東海第二発電所 工事計画審査資料
資料番号

工認-094 改2

提出年月日

平成30年7月26日

NT2 補③ Ｖ-5-9 R2

Ⅴ-5-9 計算機プログラム（解析コード）の概要・ＴＯＮＢＯＳ目次

補③ Ｖ-5-9 R2

1.

はじめに

1

1.1 使用状況一覧

22.

解析コードの概要

33.

解析手法

4

3.1 基礎方程式

4

3.2 飛来物運動方程式の時間積分法

6

3.3 評価条件

73.4 適用範囲

84.

解析フローチャート

95.

解析コードの検証

10

5.1 フジタスケールとの比較

10

5.2 米国Grand Gulf 原子力発電所への竜巻来襲事例

11

5.3 佐呂間竜巻での車両飛散事例

136.

妥当性確認

197.

参考文献

19

1.

はじめに
本資料は，添付書類において使用した計算機プログラム（解析コード）TONBOSについて説明するものである。

補③ Ｖ-5-9 R2

本解析コードを使用した添付書類を示す使用状況一覧，解析コードの概要を以降に記載する。

11.1 使用状況一覧
使用添付書類
竜巻の影響を考慮する施設及び固縛対象物の選定

NT2 補③ Ｖ-5-9 R2

Ⅴ-1-1-2-3-2

バージョン

2Ver.3

2.

解析コードの概要
コード名

項目

TONBOS

使用目的

竜巻により発生する飛来物の飛散距離等の評価

開発機関

一般財団法人 電力中央研究所

開発時期

2013年

使用したバージョン

Ver.3
TONBOS（以下，「本解析コード」という。）は，一般財団法人電力コードの概要

中央研究所にて開発・保守されているプログラムである。
空気中の物体が受ける抗力，揚力による運動を計算することで，竜巻による風速場の中での飛来物の軌跡を解析することができる解析コードであり，飛来物の飛散距離等の算出が実施できる。
仮定する風速場は，フジタモデルDBT-77（DBT:Design BasisTornado）とする。

NT2 補③ Ｖ-5-9 R2【検証（Verification）】
本解析コードの検証の内容は次のとおりである。
・自動車飛散解析において，フジタスケールの各スケールに対応する被災状況とおおむね合致した結果が得られた。
・パイプ飛散解析において，Grand Gulf 原子力発電所への竜巻襲来事検証（Verification）
及び例とおおむね合致した結果が得られた。
・自動車やトラックの飛散解析において，佐呂間竜巻での車両飛散事
例とおおむね一致した結果を得られた。

妥当性確認（Validation）【妥当性確認（Validation）】
TONBOSを竜巻により発生する飛来物の飛散距離等の評価に使用することは次のとおり妥当である。
・本解析コードは，竜巻により発生する飛来物の飛散距離等の評価を目的に開発されたコードであり，使用目的が合致している。
・評価は妥当性を確認している範囲内で行うようにしている。

3

2.1

TONBOS Ver.3の特徴「TONBOS」は，竜巻による風速場は地面付近の風速場をよく表現できているフジタモデル

DBT-77（DBT:Design Basis Tornado）が移動するものとして定め，飛来物を質点系モデルとして，飛来物が（流体）抗力，揚力及び重力を受けて運動する軌跡に対する運動方程式を解くことにより，飛来物の飛散距離，飛散高さ及び飛散速度（以下，「飛散距離等」という。）を算定する解析コードである。
3.

解析手法

3.1 基礎方程式(1) 竜巻の風速場
フジタモデルDBT-77 における接線風速等の関係式については，Fujita Work Book(1)と同NT2 補③ Ｖ-5-9 R2様に以下の式で定義する。

無次元座標

r R/Rm

接線風速

Vθ Fr(r)Fh(z)VmFr(r)

r(r 1)z Z/Hizk0

Fh(z)1/r (r 1)(z 1)

exp( k(z 1)) (z 1)0半径方向風速Vrtanα0Vθtanα0ν
1 21 ν2r(v r 1)

Vθtanα0(r 1)A(1 z1.5 )Vz(z 1)B 1 exp( k(z 1)) (z 1)
7上昇風速

(r v)
2

83 ηVmA(16z 6 7z 3 )(z 1)
28 1 ν2
ηVm B exp( k(z 1))
2 exp( k(z 1)) (z 1)k(1 ν2 )k 0 ,k,ν,η, A,Bは定数

図3-1Vθ
VrVzVmRm

接線方向風速半径方向風速（内向きが正）
上昇風速
最大接線風速外部コア半径フジタモデルの概要

また，上記以外の定数として，内部コアの半径Rνと外部コアの半径Rmの比ν＝Rν /Rm，流入層高さHi，その他定数k0，k，A，B についても，Fujita Work Book(1)に基づく以下の式及び値を用いる。
内部コア半径と外部コア半径の比
流入層高さ

(1)

0.9 0.7exp( 0.005R m )Hi ηRm4

ここでη 0.55(12)

(2)

k0

その他定数1/6 ， k0.03 ， A0.75 ， B0.0217(3)

時刻t＝0 で竜巻中心は原点に位置するものとし，竜巻の風速場は飛来物の影響を受けないと仮定すると，時刻t における竜巻の風速場Vw の各風速成分（u,v,w）は以下のようになる。
uvw

Vrcosθ V sinθVrsinθ V sinθ
Vz

Vtr
00(4)

ここで，Vtr はx軸方向の竜巻の移動速度，θは竜巻中心と各点を結ぶ線分がx軸となす角であり，竜巻中心からの距離rは以下で定義されるものである。
r(x - Vtr t) 2(5)y2

なお，静止する観測者からみた最大水平風速VD は，x軸方向の竜巻の移動速度VtrをVmに加えた以下で計算される。
VDVm

(6)

V tr

NT2 補③ Ｖ-5-9 R2(2) 飛来物の運動
飛来物の運動は，既往研究(2)〜(6)と同様に方向性がない平均的な抗力（抗力係数一定）と重力のみが外力として作用するものとした。ただし，地面付近では地面の存在により地面の損座により流れが非対称になること（地面効果）による揚力加速度L を考慮した。この場
合，物体飛行の空中での運動方程式は以下のようになる。
dV M
dt

ρ CDA
VW
2 m

V M (V W

V M ) (gL)k(7)

ここで，各変数は以下で定義する。
VM：飛来物の速度ベクトル
VW：風速ベクトル（フジタモデルでの風速場と移動速度の和）
ρ：空気密度CDA：飛来物の各方向の抗力係数と見附面積の積の平均値
m：飛来物の質量g：重力加速度k：鉛直上向きの単位ベクトル
また，揚力加速度 は，既往の風洞実験の結果（7）（8）を考慮して，物体高さdの物体にかかる揚力が，物体底面が地面から3dの高度であるときに消滅するものとした。具体的には，地面からzの距離（高度）にある物体に作用する揚力加速度L を以下の関数形でモデル化する。（Z：物体底面の高度（=z-d/2））

5

1 CDAρ
VW2mL(8)2V M x,y f(Z/d)

ここでf(Z/d)は，既往の風洞実験の結果(7)を参考に，以下で定義する。
1 (z/3d)f(Z/d)/ 10(Z/d)

0Z
3d3dZ

(9)

3.2 飛来物運動方程式の時間積分法
時刻t における竜巻の風速場Vw＝(u,v,w)は，式(4)で与えられるため，式(7)を時間積分することによって，飛来物の速度と位置の時刻歴を求めることができる。一定加速度法では，時刻τから時刻τ＋⊿τにおける加速度が一定と仮定する。つまり，dV M(t)
dt

(10)

A(τ)

ここで，加速度ベクトルA(τ)は，以下を表す。

NT2 補③ Ｖ-5-9 R2A(τ)ρ CDAVW(τ)
2 m

V M(τ)

V W(τ)V M (τ)(gL)k(11)

式(11)を時間積分すると，以下を得る。
VM(t)

VM(τ)

A(τ)(t(12)τ)

さらに式(12)を時間積分すると，以下の飛来物位置XM (t)を得る。
X M(t)X M(τ)

VM(τ)(t -τ)A(τ)(t τ) 22(13)

従って，時刻t＝τ＋⊿τにおける飛来物の速度と位置は以下で与えられる。
V M(tΔτ)V M(τ)

A(τ)ΔτX M(tΔτ)X M(τ)

VM(τ)Δτ(14)

A(τ)Δτ 22(15)式(14)，(15)の右辺は時刻t＝τにおける既知の速度・位置の関数であるため，陽解法として時間積分を実施する。

6

3.3 評価条件
以下の初期条件を設定する。
・竜巻の風速場（図3-2）
・竜巻の最大水平風速（Vm）が85 m/sとなる風速場を設定する。
・竜巻の移動速度（Vtr）は15 m/sとする。
・竜巻風速半径（外部コア半径：Rm）は30 mとする。
初期高さ（物体の設置高さ。地面設置の場合は初期高さ0 m）に51×51 個の物体を配置す

ｙる。

2Rm
本評価（多点配置）

Vｍ

ＮＲＣガイド（1点配置）

第3象限のみ例示

Rmｘ
NT2 補③ Ｖ-5-9 R2VＴ

図3-2 竜巻の風速場及び初期物体配置

7

・飛散距離等
図3-3に示すとおり，飛散距離は初期位置からの距離とし，飛来物が地表面に落ちた場合は，それ以上移動しないものとする。飛散高さは，初期高さからの相対高さとする。
2,601（51×51） 個の算出結果のうち，飛散距離，飛散高さ及び飛散速度の最大値を，解NT2 補③ Ｖ-5-9 R2

析結果として採用する。

図3-3 飛散距離及び飛散高さ
3.4 適用範囲
解析は以下の条件で実施することとする。
・竜巻の最大風速：100 m/sとする。
・竜巻の最大接線風速：最大風速との比が0.85となる85 m/sとする。
・竜巻中心の移動速度：最大風速との比が0.15となる15 m/sとする。
・竜巻の最大接線風速半径：30mとする。
・空力パラメータ：制限なし。
・解析時間刻み：0.01 秒以下。
・重力加速度，空気密度：それぞれ9.80665 m/s2，1.22 kg/m3とする。

8

4.

解析フローチャート
図4-1に解析フローチャートを示す。

ＳＴＡＲＴ

入力データの読み込み（最大風速，最大水平風速，竜巻中心の移動速度，最大風速半径，空力パラメータ，NT2 補③ Ｖ-5-9 R2

重力加速度，空気密度等）

飛来物の飛散距離等の算定

解析結果の出力

ＥＮＤ

図4-1 解析フローチャート

95.

解析コードの検証

5.1 フジタスケールとの比較
フジタスケールは，竜巻等の突風により発生した建築物や車両等の被害状況から，当時の竜巻風速を推定するために考案された指標である。このフジタスケールで示されている自動車の被災状況を表5-1に示す。
ここで，TONBOSによって，各スケールに対応する最大風速(69 m/s，92 m/s，116 m/s)におけるフジタモデルによる自動車飛散解析を行う。その結果を表5-2に示す。
TONBOSによる自動車飛散解析の結果は，各スケールに対応する自動車の被災状況とおおむね合致していると考えられる。
表5-1 フジタスケールで示されている自動車の飛散状況
フジタ風速

スケール

[m/s]

F250-69

F370-92

自動車の被災状況
cars blown off highway（自動車が道路からそれる）
cars lifted off the ground（自動車が地面から浮上する）

NT2 補③ Ｖ-5-9 R2

cars thrown some distances or rolled considerable
F493-116distances（自動車がある距離を飛ばされる，又は，かなりの距離を転がる）

表5-2 TONBOSによる自動車の飛散解析結果（CDA/m=0.0052 m2/kg）(9)
計算結果

フジタ

最大竜巻

竜巻

スケール

水平風速

接線速度

移動速度

速度

距離

高さとの対応

[m/s]

[m/s]

[m/s]

[m/s]

[m]

[m]

F26959101.01.40F392791323341.1

F4116991742593.1

105.2 米国Grand Gulf原子力発電所への竜巻来襲事例
1978年4月17日に米国のミシシッピー州にて建設中のGrand Gulf原子力発電所にF3の竜巻が来襲した。主な被害として，建設中の冷却塔内部に設置されていたコンクリート流し込み用のクレーンが倒壊し，冷却塔の一部が破損したことが挙げられる。また，竜巻によりトレーラーが台から剥がれ移動したことや，直径8〜10インチの木が折れた事例等も確認されており，図5-1は，竜巻による飛来物の飛散状況が定量的に分かる事例として，資材置き場のパイプの飛散状況を示したものである。なお，通過時の竜巻規模はF2 であったと考えられている。このパイプはコンクリート・石綿製で，長さは8フィート，直径（内径）は8インチであった。このパイプの飛散状況に対して，TONBOS を用いた飛来解析を行った。その計算条件は過去の記録

NT2 補③ Ｖ-5-9 R2に基づき表5-3のとおりとする。

図5-1 Grand Gulf原子力発電所資材置き場におけるパイプの散乱状況(10)
被害状況

・パイプを収納した木箱（一部は二段重ね）は浮上せずに転倒し，パイプが周辺7 m〜9 mに散乱。
（Pieces of pipe were scattered over the area, but none traveled morethan 25-30 ft. The pipe joints are 8 in. dia x 8 ft long.(24)）

11表5-3 Grand Gulf原子力発電所の竜巻によるパイプ飛散の再現をする上での計算条件(9)
竜巻条件

飛来物条件

初期配置

竜巻の最大風速67 m/s

最大接線風速53.6 m/s

移動速度13.4 m/s

コア半径45.7 m/s

直径（外径）

9 inch（0.2286 m）

物体高さ

0.229 m密度1700 kg/m3

飛行定数CDA/m

0.0080 m2/kg

・物体個数51×51 個，竜巻半径の4 倍を一辺とする正方形内（x, y=[-2Rm,+2Rm ]）に等間隔配置
・設置高さ1 m（パイプが収納されていた木箱が2 段重ねで配置されていた状況を想定。）

計算結果を表5-4に示す。TONBOSによるパイプ飛散解析の結果は,パイプがほとんど飛散せず，木NT2 補③ Ｖ-5-9 R2箱が倒れた影響で散乱したと思われる状況とおおむね合致している。
表5-4 Grand Gulf原子力発電所のパイプの飛散計算結果(9)風速場モデル
フジタモデル

計算結果

初期物体
高さ

飛散距離

飛散高さ＊１

最大水平速度1 m

1.2 m

0.0 m4.9 m/s

注記 ＊1：初期物体高さからの飛散高さ。

12

5.3 佐呂間竜巻での車両飛散事例
2006年11月7日に北海道網走支庁佐呂間町に発生した竜巻（以下「佐呂間竜巻」という。）により，4tトラックが約40 m移動したことが報告されている(11)。被災状況を図5-2に示す。この事例では被災時に4tトラックに乗員2名が乗車しており，4tトラックの初期位置と移動位置が分かっている（図5-2左上画像の②）。また，4tトラックの他に2 台の自動車（図5-2左上画像の③と⑥）について，初期位置と被災後の移動位置が分かっている。このように竜巻被災前後で車両等の位置が明確になっている事例は極めてまれである。なお，竜巻飛来物の再現計算は，竜巻が頻発する米国でもほとんど実施されていない。この理由としては，来襲した実際の竜巻特性を精度よく計測・推測することが困難であることや自動車等の移動前後の位置が不明

NT2 補③ Ｖ-5-9 R2

確な場合が多いことが挙げられる。

図5-2 佐呂間竜巻(2006.11.7)による被災状況（工事事務所敷地内の車両被災）（文献(11)の写真に竜巻被害の方向をで加筆）

ここでは，TONBOSによるフジタモデルを風速場として用いた車両（4tトラック，乗用車）の飛散評価を行い，実際の被害状況と比べて妥当な結果となるかどうかの確認を行う。
方法としては，下記の2通りとする。
(a) 竜巻特性や飛来物（4tトラック，乗用車）の状況を現実的に設定した場合の再現解析(b) 東海第二発電所に適用する飛散解析手法による検証

13

(a) 竜巻特性や飛来物の状況を現実的に設定した場合の再現解析(ⅰ) 4tトラックの飛散解析
再現解析の条件として，入手可能なデータ(11)(12)に基づき，合理的と考えられる竜巻特性条件と飛来物（4tトラック）の条件を表5-5のように設定する。初期配置の条件として，配置個数は１個とし，竜巻が遠方から近づく状況設定としている。また，風速60 m/s
以下では浮上しない設定となっている。その上で，竜巻との距離を合理的な範囲で変化させ，佐呂間竜巻の再現性を確認する。
車両と竜巻中心との距離を18 m，20 m，22 m とした場合の解析結果を表5-6及び図5-3に示す。車両の軌跡は竜巻中心との相対位置関係に敏感であるが，各ケースとも飛散方向が実際の移動方向とおおむね合致しており，特に車両と竜巻中心との距離を20 m とした
ケース2 では飛散距離もほぼ正確に再現されている。このように，TONBOSによる4tトラックの飛散解析の結果で，飛来物が地上に設置された状況からの飛散挙動が再現できることが確認できる。
表5-5 佐呂間竜巻の4t トラックの計算条件(9)NT2 補③ Ｖ-5-9 R2

竜巻条件

飛来物条件

竜巻の最大風速

92 m/s

最大接線風速70 m/s

移動速度22 m/s

コア半径20 m

車種不明のため，三菱

車両長さ

8.1 m

ふそうPA-FK71Dの仕様

車両幅2.24 mを採用

車両高さ2.5 m車両質量4000 kg
0.0056 m2/kg

飛行定数CDA/m
初期配置

・物体個数：1台
・竜巻は遠方から物体に近づくが，風速60 m/s 以下では浮上しない
・設置高さ0 m

表5-6 佐呂間竜巻での4tトラックの飛散計算結果(9)計算結果（フジタモデル）

解析

車両と竜巻中心とケースの距離

飛散距離

飛散高さ

最大水平速度

１22 m

45.4 m2.8 m

25.8 m/s

２20 m

35.5 m2.3 m

22.2 m/s

３18 m

25.9 m1.7 m

18.8 m/s

14図5-3 TONBOSによるトラック飛散の再現解析結果(9)（ⅱ）乗用車の飛散解析
白い乗用車（図5-2の⑥）の被災事例を対象として，物体を１点初期配置した条件で最大水平速度等を計算する。
NT2 補③ Ｖ-5-9 R2

白い乗用車の計算条件について，表5-7に示す。
表5-7 佐呂間竜巻の白い乗用車の計算条件
竜巻条件

飛来物条件

表5-5と同様
白い自動車

車両長さ4.4 m

トヨタカローラ

車両幅1.7 mを仮定

車両高さ

1.5 m
0.0097 m2/kg

飛行定数CDA/m
・物体個数：1台
初期配置

・竜巻は遠方から物体に近づくが，風速60 m/s 以下では浮上しない
・設置高さ0 m

白い乗用車と竜巻中心との距離を，18 m，20 m，22 m とした場合の解析結果を表5-8及び図5-4に示す。飛散距離についてはケース1 でおおむね合致している。
飛散方向については，飛び出し方向はおおむね合致しているものの，最終的な着地点には多少のずれが生じている。これは乗用車（白）が建物に近接して駐車していたため，この建物の倒壊の影響を受けて飛散方向のずれが生じたものと推定される。
なお，赤い乗用車（図5-2の③）について評価した場合は，竜巻中心との距離が大きいため飛散しない解析結果となる。ただし，実際には，赤い乗用車は全壊・飛散したプレハブ建物（軽量鉄骨造２階建て，図5-2のA）のすぐ下流側に駐車しており，そのがれきの影響を受けて一緒に移動したものと考えられる。

15表5-8 佐呂間竜巻での白い乗用車の飛散計算結果
計算結果（フジタモデル）

白い乗用車と竜巻

ケース

中心との距離飛散距離

飛散高さ

最大水平速度

１22 m

51.9 m3.6 m

28.9 m/s

２20 m

43.5 m3.4 m

24.7 m/s

３18 m

34.7 m2.9 m

21.1 m/s

NT2 補③ Ｖ-5-9 R2解析

図5-4 TONBOSによる白い乗用車飛散の再現解析結果

16

(b) 東海第二発電所に適用する飛散解析手法による検証(13)
ここでは，東海第二発電所に適用する飛来物評価法の竜巻条件・物体初期配置条件で前述の佐呂間竜巻における4tトラック及び白い乗用車の被災事例を評価し，佐呂間竜巻での実際の被災状況（移動距離等）との結果を比較する。
（ⅰ） 4tトラックの飛散解析
計算条件について表5-9に示す。竜巻条件としては，設計竜巻の最大風速を92 m/sとし，その他の特性量については，ガイドに例示されている方法に従い，移動速度Vtを14 m/s（最大風速の15 %），竜巻コア半径Rm を30 m とする。
表5-9 東海第二発電所に適用する飛来物評価法の計算条件

竜巻条件

飛来物条件

設計竜巻風速92 m/s

最大接線風速78 m/s

移動速度14 m/s

コア半径30 m

表5-5と同様

NT2 補③ Ｖ-5-9 R2・物体個数51×51 個，竜巻半径の4 倍を一辺とする正方形内（x, y=[-2Rm,初期配置+2Rm]）に等間隔配置
・設置高さ0m

表5-10に実際の被災状況と，東海第二発電所に適用する飛来物評価法の結果の比較を示す。また，図5-5に被災後の4tトラックの状況を示す。
TONBOSによる飛散評価結果として，4tトラックの最大飛来物速度は36 m/s，最大飛散高さは3.6 m，最大飛散距離は63.4 m となる。
実際の4tトラック飛散距離は約40 mであり，TONBOSによる飛散距離の評価結果はこれを上回る。また，飛散高さや最大水平速度については，直接の比較はできないものの，4tトラックの乗員2 名が存命であったこと，被災後の4t トラックがほぼ元の外形をとどめていること等から，東海第二発電所に適用する飛来物評価法で飛散解析をした場合でも，実際の被災状況と比較して妥当な結果となるものと考えられる。

17

表5-10 実際の被災状況と「東海第二発電所の飛来物評価法」との結果の比較（4tトラックの場合）
風速場モデル
フジタモデル（地上）

飛散距離

飛散高さ＊ 1

63.4 m

3.6 m

最大水平速度
36.0 m/s(毎時130 km)
被災後もほぼ元の外形をと実際の被害状況

約40 m

乗員2名が幸いにも存命で救出どめていることが示されてされ，搬送先の病院で聞き取おりり調査に応じており(11) ，トラ(11)

ックが地面から3.6 m以上の高解析で得られた最大速度

所から落下したとは考えにく（約130 km/h）を超える飛い。

散速度であったとは考えに，実際の飛来物速度は本くい。

NT2 補③ Ｖ-5-9 R2

注記 ＊1：初期物体高さからの飛散（浮上）高さ

図5-5 竜巻による被災後の4t トラックの様子(12)(14)（ⅱ） 乗用車（白）の飛散解析
4tトラックの場合と同様に，東海第二発電所に適用する飛来物評価条件で白い乗用車の飛散解析を行った場合の結果を表5-11に示す。
白い乗用車の場合も，TONBOSによる評価が実際の被災状況を包含する結果となっている。

18

表5-11 実際の被災状況と「東海第二発電所の飛来物評価法」との結果の比較（白い乗用車の場合）
風速場モデル
フジタモデル（地上）
実際の被災状況

飛散距離

飛散高さ＊ 1

最大水平速度

82.3 m4.2 m

44.1 m/s

約50 m−−

注記 ＊1：初期物体高さからの飛散高さ。
6.

妥当性確認
TONBOSを竜巻により発生する飛来物の飛散距離等の評価に使用することは次のとおり妥当である。
・本解析コードは，竜巻により発生する飛来物の飛散距離等の評価を目的に開発されたコードであり，使用目的が合致している。
・評価は妥当性を確認している範囲内で行うようにしている。

7.

参考文献

NT2 補③ Ｖ-5-9 R2(1)

Fujita, T. T., Workbook of tornadoes and high winds for engineering applications,U. Chicago, 1978

(2)

Simiu, E. and Cordes, M., Tornado-Borne Missile Speeds, NBSIR 76-1050, 1976.

(3)

Maruyama, T., Simulation of flying debris using a numerically generated tornadolike vortex. Journal of Wind Engineering and Industrial Aerodynamics, vol.99(4)，pp.249-256，2011.

(4) 東京工芸大学（2011）：平成21〜22 年度原子力安全基盤調査研究（平成22 年度）竜巻による原子力施設への影響に関する調査研究，独立行政法人原子力安全基盤機構(5)

U.S. Nuclear Regulatory Commission, Regulatory Guide 1.76: Design-Basis Tornado
and Tornado Missiles for Nuclear Power Plants, Revision 1, March 2007.

(6) 丸山敬，数値的に生成された竜巻状の渦の性質，平成２１年度京都大学防災研究所研究発
表要旨(7)

EPRI, Wind field and trajectory models for tornado propelled objects, report NP2898, 1978.

(8)

Lei, C., Cheng, L. and Kavanagh, K., Re-examination of the effect of a planeboundary on force and vortex shedding of a circular cylinder, J. of Wind
Engineering and Industrial Aerodynamics, Vol.80，pp.263-286, 1999.

(9) 日本保全学会 原子力規制関連事項検討会， 2015：軽水型原子力発電所の竜巻影響評価
における設計竜巻風速および飛来物速度の設定に関するガイドライン(JSM-NRE-009)(10) Fujita, T. T., and J. R. McDonald, Tornado damage at the Grand Gulf, Mississippi
nuclear power plant site: Aerial and ground surveys, U.S. Nuclear Regulatory
Commission NUREG/CR-0383, 1978.

19

(11) 札幌管区気象台:平成18年11月7日から9日に北海道(佐呂間町他)で発生した竜巻等の突風.
災害時気象調査報告, 災害時自然現象報告書, 2006年第1号,2006.
(12) 奥田泰雄,喜々津仁密,村上知徳,2006年佐呂間町竜巻被害調査報告.建築研究所災害調査,
46, 2006.
(13) 江口譲，杉本聡一郎，服部康男，平口博丸, 原子力発電所での竜巻飛来物速度の合理的評
価法 （Fujita の竜巻モデルを用いた数値解析コードの妥当性確認），日本機械学会論文集， Vol.81, No.823, 2015.

補③ Ｖ-5-9 R1E(14) 土木学会 平成18年11月北海道佐呂間町竜巻緊急災害調査報告書

0