竜巻被害 EF スケール・デブリ飛距離シミュレーター 戻る
気象災害・耐風建築

竜巻被害 EF スケール・デブリ飛距離シミュレーター

EFスケール(EF0〜EF5)の風速・デブリ(飛翔物)の質量・投影面積・形状から、抗力・飛翔臨界風速・終端速度・衝撃エネルギー・住宅被害区分をリアルタイムで計算します。竜巻被害の主因である「飛翔物の衝突」を定量的に評価できます。

パラメータ設定
EF スケール
Enhanced Fujita スケール(3秒平均風速)
デブリ質量 m
kg
投影面積 A
高さ h
m
デブリの初期高さ(地表からの高度)
形状
抗力係数 Cd を決定
旋回半径 R
m
中心圧力低下 Δp
mb
竜巻中心と外周の気圧差(mb = hPa)
計算結果
風速 v (m/s)
抗力 F_drag (N)
飛翔臨界風速 v_crit (m/s)
終端速度 v_term (m/s)
衝撃エネルギー KE (kJ)
被害区分
竜巻・デブリ・建物の可視化

中央の渦が竜巻、飛翔するデブリが軌跡を描き、右側の建物に衝突します。色は衝撃エネルギーの大きさ(緑→橙→赤)を表します。

飛翔距離 vs 風速
EFスケール別の衝撃エネルギー比較
理論・主要公式

$$F_{\text{drag}} = \tfrac{1}{2}\,\rho\,C_d\,A\,v^{2}, \qquad v_{\text{crit}} = \sqrt{\dfrac{m\,g}{0.5\,\rho\,C_d\,A}}$$

ρ=空気密度(1.225 kg/m³)、Cd=抗力係数(棒1.0、板1.5、ブロック2.0)、A=投影面積、v=風速、m=質量、g=9.81 m/s²。

$$v_{\text{term}} = \sqrt{\dfrac{2\,m\,g}{\rho\,C_d\,A}}, \qquad KE = \tfrac{1}{2}\,m\,v_{\text{impact}}^{2}$$

終端速度 v_term は自由落下時の上限速度、衝撃エネルギー KE は飛翔デブリの運動エネルギー。住宅被害は KE が 1/10/100 kJ を超えるごとに段階的に重症化する。

竜巻被害評価 EF スケール・デブリ飛距離

🙋
竜巻の強さって「EFスケール」って言うんですよね?昔の藤田スケールと何が違うんですか?
🎓
いい質問だね。藤田哲也博士が1971年に作った元のFスケールは、被害の写真から風速を読み取る方式だった。ただ「住宅全壊=F4で約100 m/s」みたいに高めに見積もる傾向があってね。米国NOAAが2007年に「Enhanced Fujita(EF)」へ改良して、建物種別ごとの被害指標(28種類)と、3秒平均風速の対応表を作り直した。例えば旧F5は117 m/s以上だったのが、EF5は90 m/s以上に見直されたんだ。
🙋
なるほど。じゃあ、竜巻の被害ってやっぱり「風そのもの」で家が壊れるんですか?
🎓
それが意外と違うんだよ。実は被害の主因は「飛翔物(デブリ)の衝突」なんだ。EF2クラス(55 m/s)でも、屋根材の2x4 lumber(約2.4 kg)が30-50 m/sで飛んでくると、窓ガラスや薄い外壁を簡単に貫通する。シミュレーターの左側で「形状=板状、質量10 kg」のデフォルトを見ると、抗力 F_drag=139 N で飛翔臨界 v_crit=46 m/s、EF2の55 m/sで余裕で浮揚するのが分かるよ。
🙋
それは怖い…!「Storm Shelter」っていうのは、そのデブリ衝撃に耐える部屋ですか?
🎓
その通り。米国ではFEMA P-361とICC 500っていう基準があって、住宅向け耐力室は250 km/h(約70 m/s)の風と、2x4 lumberを約 45 m/sで投射する衝撃試験に耐える必要がある。典型的な構造は3層CMUブロック+1/8インチ鋼板で、ドアも同等の貫通抵抗を持つ。トルネードアレー(米中部)では新築住宅に標準装備されつつあるよ。
🙋
日本でも竜巻ってあるんですよね?
🎓
あるよ、しかも結構深刻。2006年北海道佐呂間町でF3が発生して9名亡くなった事例や、2012年つくば市のF3で1名死亡・住宅大規模被害があった。米国の年1200件と比べれば少ないけど、日本では2016年から「日本版改良藤田スケール(JEF)」が運用されていて、JEF3で建物2階部分の大破などの指標が決められている。シミュレーターの「中心圧力低下 50 mb」は典型的なEF2-3クラスの値だから、日本の竜巻でも十分あり得る範囲だ。
🙋
気候変動で竜巻って増えてるんですか?
🎓
頻度自体は横ばい〜微減という研究もあるけど、「強い竜巻(EF3以上)の比率増加」と「Tornado Outbreak(1日に複数発生する集中発生)」の増加が報告されている。温暖化で対流不安定(CAPE)が強まる一方、上空風シアーは減るという相反する効果のせいで地域差が大きいんだ。日本でも台風随伴の竜巻が増える可能性が指摘されていて、耐風設計の見直しが進んでいる。Doppler weather radar(NEXRAD WSR-88D)やMobile Doppler(DOW)でリアルタイム検知の精度も上がってきたよ。

よくある質問

EF(Enhanced Fujita)スケールは、米国NOAAが2007年に旧Fujita-Pearsonスケールを改良して採用した竜巻強度の指標で、被害状況から推定3秒平均風速を逆算します。EF0は29-38 m/s、EF1は39-49、EF2は50-60、EF3は61-74、EF4は75-89、EF5は90 m/s以上で、住宅全壊・自動車飛翔・コンクリート造大破などの被害指標と関連付けられます。考案者の藤田哲也博士に由来し、日本では2016年から日本版改良藤田(JEF)スケールが運用されています。
デブリが浮揚を開始する臨界風速は、抗力と重力が釣り合う条件 F_drag = m·g から導かれ、v_crit = √(2mg / (ρ·Cd·A)) です。例えば 10 kg・投影面積 0.05 m²・板状(Cd=1.5)のデブリでは v_crit ≈ 46 m/s で、EF2クラス(50-60 m/s)から浮揚が始まります。住宅で典型的な2x4 lumber(約2.4 kg)は30-50 m/sで飛翔し、窓や壁を貫通する主要な被害メカニズムになります。
米国ではFEMA P-361とICC 500がコミュニティおよび住宅向けStorm Shelterの基準で、250 km/h(約70 m/s)の風と、2x4 lumberを約 45 m/sで投射する飛翔物試験に耐える耐力室を要求します。典型構造は3層CMU(コンクリートブロック)に1/8インチ鋼板を組み合わせたもので、ドアも同等の貫通抵抗を持たせます。日本でも一部の自治体やメーカーが鉄筋コンクリート造の「耐風シェルター」を提案しています。
米国では年間約1200件、世界では年間2000件程度の竜巻が観測され、米国だけで年80-100名の死者が出ています。近年は強い竜巻(EF3以上)の比率増加と「Tornado Outbreak(集中発生)」の増加が報告されており、温暖化に伴う対流不安定の強化との関連が議論されています。日本でも2006年北海道佐呂間F3(9名死亡)、2012年つくばF3など、深刻な被害が継続的に発生しています。

実世界での応用

耐風建築の設計:米国トルネードアレーでは、新築住宅にFEMA P-361準拠のStorm Shelterが標準装備されつつあります。3層CMU(コンクリートブロック)と1/8インチ鋼板を組み合わせた構造で、2x4 lumberを約45 m/sで投射する飛翔物試験を満たします。住宅本体の屋根構造もハリケーンタイ(金物)でトラスを基礎に連結し、屋根剥がれ起点の倒壊を防止します。

Doppler weather radarによる検知:米国NOAAは全国に約160基のNEXRAD WSR-88D(Sバンドドップラーレーダー)を配備し、メソサイクロン(低気圧性の渦)とTVS(Tornado Vortex Signature)をリアルタイム検出します。空港にはTDWR(Terminal Doppler Weather Radar)が追加配備され、Mobile Doppler(DOW: Doppler on Wheels)が研究目的で実測風速を取得します。これらにより竜巻警報の平均リードタイムは13分にまで延伸しています。

保険・損害査定への適用:米国の住宅保険会社は、被害住宅の屋根・外壁・窓の損傷写真から逆算でEFスケールを推定し、雹害・風害・洪水被害との切り分けを行います。本ツールのような飛翔物エネルギー計算は、損害請求の妥当性検証(claim verification)に使われます。日本でも風水害保険の査定で類似手法が研究されています。

気象工学・CAEとの連携:竜巻シミュレーションは、LES(Large Eddy Simulation)やTLV(Translating Large-eddy Vortex)モデルで実施されます。本ツールのような梁理論・抗力理論の概算は、詳細CFD解析の前段で「どの程度のデブリが浮揚し得るか」を当たり付ける用途に有用で、メッシュ作成・境界条件設定の参考になります。

よくある誤解と注意点

まず最大の落とし穴が、「EFスケールの風速=3秒平均風速」であることを忘れることです。EFスケールの風速値(EF5なら90 m/s以上)は、地表から10 m高さでの3秒平均風速の推定値であり、瞬間最大風速ではありません。同じ竜巻でも瞬間値は1.3〜1.5倍に達するため、構造設計では「3秒風速×ガストファクター」で換算する必要があります。本ツールでは EF_winds として中央値(EF2なら55 m/s)を採用していますが、実設計では各EFランクの上限値や瞬間最大相当を使うのが安全側です。

次に、「抗力係数 Cd を単純な定数とみなす」こと。本ツールの Cd は形状別の代表値(棒1.0、板1.5、ブロック2.0)ですが、実際の Cd は Reynolds数、デブリの姿勢(迎角)、表面粗さ、回転状態で大きく変動します。特に飛翔中のデブリはタンブリング(不規則回転)しており、瞬間的に Cd が2倍以上になることもあります。FEMA P-361の飛翔物試験で2x4 lumberを使うのは、最も貫通力の高い「先端衝突」を想定した保守的設定です。

最後に、「衝撃エネルギーだけで被害を評価しない」こと。本ツールは KE = (1/2)·m·v² で運動エネルギーを計算しますが、実際の貫通被害は「単位面積あたりのエネルギー」つまり KE/A や、貫通深さに支配されます。同じ KE でも、釘の先端のような小面積衝突は鋼板を貫通しますが、平板の側面衝突は曲げ変形で吸収されます。詳細な被害評価には、衝撃力(F=dp/dt)と接触面積、対象物の動的吸収特性を考慮した非線形FEM解析が必要です。