【論 文
1
UDC :551
.
46 :624.
042 :627.
223.
6日本建築 学 会 構造 系 論 文報 告 薬 塘 428 号
・
1991 年10月Jolirnal ot Stluct
.
ConstT.
Engng,
AIJ,
Ne.
428,
0ct,
1991曳航 中
の
コー ン
形
状構
造
物
に
作
用
す
る
流体 力
と
波 浪 強 制 力
お
よ び
運
動特性
に
関
す
る
研
究
第
1
報
流
体
力
と波 浪強 制 力
STUDY
ON
THE
CHARACTERISTICS
OF
THE
且YDRODYNAMIC
FORCES
,
WAVE
EXCITING
FORCES
AND
MOTIONS
ON
CONE
TYPESTRUCTURE
WITH
A
CONSTANT
FORWARD
SPEED
.
IN
WAV
耳
S
Part
l The
hydrodynamic
forces
and the wave excitingforces
藤
澤
康
雄
* ,増 田 光
一
* * ,前 田 久 明
* * *Yasuo
FUJISAWA
,
Koichi
MASUDA
andHisaaki
A4AEDA
The
authorsdeveloped
a c6mputer programby
using thr6qdimensional
sourcedistribution
method considering the
fo
[wa [d
speedbased
onlinear
potential theory.・
The
authors conducted model experiments on a.
cone typed structureduring
towing and non・
toWing conditions and com.
pared calculated results with experimental resu 且ts regarding hydrodynamic
forces
,
wave excitingforces
and the effects offorward
speed and clarified the characteristics of these results.
The
calculated results and experimental results with a constantforward
speedin
waves are shown a good dorrelation within thelinear
potential theory.
Kegwerds :
forward
sPeed,
COiteC
コPed
strttcture,
threedimenSiopal
scurcedistri
’
butiony
methed,
hydrody
namic 〆
force
,
wave exCitingforce
前 進 速 度
,
コー
ン形 状 構 造 物.
,
3次 元 特 異 点 分 布 法,
流 体 力,
波 浪 強 制 力 1.
緒 言 大 型 海 洋 構 造 物,
特に,
海 洋 石 油・
ガス掘 削 /生産 用 重 力 式プラッ トフ ォー
ムは,
ドライドック並び にフロー
ティ ングサ イ トで の建 造 後, 設 置 海 域まで 曳 航さ れ沈 設・
設置さ れ ることが 特 徴の一
っ である。 現 在,
稼 働 中 のプラッ トフ ォー
ムの事 例で は,
曳 航 時に作 用 する流 体 力と波浪
強 制 力お よ び運 動 特 性の評 価は,
水 槽 実 験に よ り行われ て い る のが実 状であり,
系 統だっ た理 論,
解 析 手 法 等の研 究はほ と ん ど見 受け ら れ ない。し か し な が ら,
建 設 段 階の一
環で ある曳 航 段 階におい て,
曳 航 時の プ ラッ トフ t一
ム に 作 用 する流 体力と波浪強制力お よ び運 動特性を 正確に把 握す ること は,
プラッ トフォー
ム の設 計に も大 き な 影 響 を与え, その安 全 性お よ び稼 働 時の コ ス トダウン等につ な が る た め非 常に重 要な検討 項 目であ る。一
般に,
波 浪中を一
定の前 進 速 度 を 持っ て移 動する3
次 元任意 形 状浮体に作 用 す る流 体 力と波 浪 強 制 力お よ びこれ らに支 配さ れ る運動の解析は,
線形理論の範 囲に おいて も著 し く複 雑である。
通常,
波 浪 中を前 進 速 度を 持っ て航行する の は船 舶に限 られ , 船 舶工学の分野にお いて多くの研究がな さ れ てい る。
船 舶 工 学で は一
般 的に,
船体動 揺 理論に おい て前 進速度を考 慮 し た簡 便な方 法と して,
ス ト リップ法M,
あるい は細 長 体 理 論2,等で取 り 扱わ れて い る が,
これ らの理 論で は,
大 型 海 洋 構 造 物の よ う な3
次元物 体の 3次 元 影 響および前 進 速 度 影 響を十 分に考 慮す るこ とは困 難であ る。 そこで最 近の高 速 大容 量の大 型 電 子計算機の出現に よ’
り,
前 進 速 度のある場 合 の 3次元定 常 動 揺問 題 に対 して,
線形 理論の範 囲で数 値 計 算 が 実 施さ れ て い る。
小林3i,Inglis
&Price
“・
5 ),
Chang6
)らは,
積 分 方 程 式の核 関 数に前 進 速度を考 慮し 713 次元Green
関数
を用い て,
任 意の 3次元 物体 (半 没 楕円船 型 ),
船 体 形 状 (Ship m 。del
Series
60)にっい て物 体 表 面に特 異 点 を分 布する特 異 点 分 布 法を用いて 前 進速度 影 響 等の研 究を行っ て い る
。
しか し な が ら,
重 力式 プラ ッ トフ ォー
ム のよ う な大型海洋構造物を曳 航 す * (株 )大 林 組 技 術 研 究 所・
工 修 * * 日本 大 学 理工学 部 海 洋 建 築工学 科 助 教 授・
工博 1 # 東 京大 学生産 技 術 研 究 所 教授・
工 博Technical Resarch Institure
,
Obayashi Corporation,
M.
Eng.
Assoc
,
Pro{.
,
Dept.
ofOceanic
Architecヒure & Engineering,
CollegedScience & Technology
,
Nihon Univ,
,
D匸.
Eng.
幽
Prof
.
,
1nstitute of Industrial Science,
Univ,
o 正Tokyo,
Dr,
Eng.
る場 合の前 進速度影 響等を系統的に検 討し た研 究 事 例は 海 外を含め非 常に少ない の が 現状であ り, い まだ解 明さ れて いない部 分が多く残さ れ て お り極め て興 味深く, 重 要な研 究 課 題であ る
.
。 そこ で本研究で は,
前進 速度影 響 を考 慮して, 曳 航 時の重 力 式 海洋構 造 物に働く流 体 力 と 波 浪 強 制 力お よ び その運 動 応 答の評 価 法を開 発す ること を目的とす る。
著者ら は,
既にケー
ソ ンタワー
タイ プ・
プラ ッ トフ ォー
ム に関して線形ポテンシャル理論に基づ く流 体 力と波 浪 強 制 力の評価法お よび 前 進 速 度 影 響 等に つ い て理論 的かつ 実 験 的 検 討 を行い一
連の成 果7)β 吃 得 た。
著 者 ら は,
これ らの研 究 成 果 を踏まえ て,
重 力 式 プ ラッ トフォー
ム の一
つ である氷 海 域 用コー
ンタ イフ・
フ ラッ トフt一
ム (以 下, コー
ン形 状 構 造 物 )を取り上 げ,
曳 航 時の流 体 力と波 浪 強 制 力の評 価 法とし て の線 形ポテ ンシャル理 論に基づ く3
次 元 特 異 点 分 布 法の適 用 性,
前 進 速 度 影 響等に関し て検 討 を行っ た。 また,
コー
ンタ イ ブモ デルに よ る曳 航 実 験を実 施して,
前 進 速 度 を 有する 場合の流体 力と波浪強制力お よび運 動 応 答 を計 測し, 実 験結果と 理論計算 結果を比 較・
検 討し た。
さ ら に,
コー
ン タ イプの よ う な特 殊な構 造 形 状のた め強 制 動 揺の振 幅 影響 お よび入射 波 高 変 化に よ る流 体 力, 波 浪 強 制 力の実 験 を行い, コー
ン形 状 構 造 物に作 用 する流 体 力および波 浪 強 制 力の線 形 範 囲お よ び非 線 形 影 響を明 確に し た。 そ こ で, 本 論 文では その第 1報と して,
曳 航 中の コー
ン形 状 構 造 物に作 用す る流 体 力と波 浪 強制力,
前進 速度影 響 お よび線形範囲 と非線形 影 響につ い て述べ る。
2.
理 論 解 析 法 2.
1 基 礎 方 程 式お よ び境界条件流 体は, 非 粘 性, 非 圧 縮 性
,
非 回 転の完 全 流 体を仮定 するので,
速 度 ポ テン シャ ル の存 在が保 証 されて いる。 ま た。
自 由表 面の変 位お よび運 動の振 幅は微 小で あり,
これ らの 高 次の微 小量 は無 視で きる もの と 仮 定する。 コー
ン形 状 構 造 物は,
自由 表 面 を横 切る物 体 面の なす 角 が一
般に鉛 直と異な る ため水 線に非 線 形 性の強い特 異 性 を有 する可 能 性がある。 し か し,
本 論 文の理 論 解 析で は 自由 表 面を横 切る物 体 面と水線との非線形性は少ない も の と仮定す るこ とにす る。
以 上の仮 定によ り,
理 論 解 析 法は線 形 理 論の範囲内で議論す る。 座標系は空間固 定 座 標系 (X
,y
,Z
)と物体座標系 (x, y,2)を 理 論 解 析に お い て適 宜 用いる もの とする。
また, 原 点 を静 止 水 面 上に取 りZ
軸は上向きを 正と す る。
構 造 物は, 無 限 水 深 中 をX
軸正方向に一
定の前進 速度U
。で進むもの と する。
境 界 面上に立てた法線は,
流体 領域か ら外向きを正と す る。
Fig.
1に座 標 系を示 す。 こ こ に,
SH
は物体 表面の境界,SF
は自 由 表 面の境 界,SR
は無 限遠方の境界を示す。 ま た,
入射 波は X 軸の正 方 向か ら負 方 向へ 進 行し て くる もの と する。 さ らに, 構 造 物は x,
z 平 面に対して対 称Fig
.
1 Coordinate systemn & n
囗
.
転であり, 構 造 物の運 動は surge, heave
,
pitchにっ い て 考え ること と する。
し たがっ て, 構 造 物の運 動に伴 う流 体 運 動もx,
z 平 面に対 して対 称と なる。
前進 速度を有す る構造物ま わ りの全 速度ポテンシャ ル は,
(1
)式の よ うに表せ る。
φ(x,
y,
z ;t
) = rUeX
+il
。
(x, y, z〕 +ip
(x,
y
, z)e’
‘tUet…・
・
『
……
(1 ) こ こに,
右 辺 第1項は構 造 物の前 進に伴 う一
様 流U
,の 速 度ポテ ン シャル,
第 2項は,
時 間に依 存せ ず一
定 速 度 で静 止 水 面上 を前 進す る と きの構 造 物 回りの流れ場 を表 す定 常 ポテンシャ ル , 第3
項は 入射 波お よび構造物に よ る非 定 常 速 度 ポテン シャ ル である。
また,
非 定 常の 速度 ポテンシャル は, 入射波の速度ポテン シャ ル φ‘と構造 物に よ り 入射 波が か く乱さ れ るdiffraction
ポテン シャ ル φ‘ お よび 構 造 物の運 動に伴 うradiation ポテン シャ ル φ。とに分けて表せ る。 すな わち,
非定 常の速 度ポテンシャ ル を φ=
φ‘+ilel
+ φ。…・
……・
・
…・
・
………・
一
(2) と する と,
(2) 式は,
Laplace 方 程 式お よ び物 体 表 面 条 件,
無 限 遠 方 条 件,
さ らに下 式の線 形 化し た自 由表 面 条 件を満 足し な け れ ば な ら ない。(
叫,
∂一
酬 否ア)
2 φ・ ・{
謬
一
q
・・S
・・
7…一
(・) 2,
2 前進速度を考慮 し たGreen
関数 無 限水深中を単一
振 幅の強さで調和 振動す る特異点が一
定速 度U
。で前 進す る場 合, 自 由 表 面 条件お よび無 限 遠 方 条 件 を 満 足するGreen
関 数の表 示は, 種々 の もの が示さ れ て い るが, 本 研 究で は,Wehausen
&Lait・
one9 } の表示に従うと 下 式の ように与え ら れ る。 な お, 前 進 速度を考慮し たGreen
関数の 導 出お よ び物理的な 検 討につ い て の詳 細は流 体 研 究グルー
プL°),
高 木・
大 楠11)の論 文に詳しい 。 G(x.,
Yp.
Zp;XQ,
YQ,
ZQ)一
÷÷ 努
鷹
゜
°
∫(翩 ・d
θ・
鷹
・・e・・1
・・dkd
・+囂
・・…剛
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(4) こ こ で,
X,,
Y
尸,
Z
.は領 域 内の任 意の点を,
XQ,
YQ,
2Q は構 造 物上 の点を, θ,
’
h
は 積分変 数,L
,,
L2
は関 数f
(θ,
k
)の特 異点を避け る た めの積分 経 路を そ れ ぞ れ表 す。 ま た,
嬬=
(Xp−
XQ)2十(Y
,−
YQ) 2 十(Z
尸十(− 1
> MZQ )z(m = =
1,
2 )・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
一・
・
…
(5
)’
ノ…
h
・一
[,、.
,。讎
。。S 、,r
…一
:
・
・(・) h(θ,
k)=
k exp [h
{Zp
十 ZQ十i
(Xp−
XQ〕cos θ1
]・
cos [h
(y
,−
YQ)sin θ]・
tt……一
(7) である。 2,
3 積 分 方 程式 上記の Green 関 数 を 用 い て構 造 物 表 面上の 点Q
(XQ,
YQ,
ZQ)に σ(Q
)の一
定 強 さの 吹 出 しが 分 布 して いると き,
構造 物まわ りの点尸
(X。,
y。,
Z。〉に おける非 定 常ポ テン シヤル φ(P)’
は,
R,
Brardi2
)の表示に従っ て次の よ うに表せ る。φ(・)
−
41 。蕉
・(ω・(・.
Q
>・・・
岩
野
・・Q
) ∂G
黷
Q
) 吻・
…・
・
(・) こ こで,
(8 )式の第 2項は線 積 分の 項で あ る。
こ の 線 積分項は, 物 体 と 自 由表 面との交 線に表れ る が こ の交 線におい て本 来 流 速は特 異 性を示す。UrselP3
トに よ れ ば,
Kelvin−Neumann
問 題にお け る線 積分項を考 慮し』
た計 算 で は,
交 線の 特 異 性が弱 く なっ た解 (least singular solution )と して意味を持つ と し てい る。一
方,
交 線で の弱い特 異 性という 条 件は物理的に あまり意 味がない]4 ) とい う意 見も あ る。
ま た, 柏 木 IS)に よ れ ば 非 定 常 問題に おいて も線 積分 項 を 含 んだ解は,Ursell
と同様に特異 性の 弱いleast
singuLar solution と考え ら れ ることを 示 唆し た が,
エ ネル ギー
保存 則は満 足され ないこ と を示 し た。
この原 因と し て線 形 理 論の定 式 化に問 題が あ る の か, 他の物 理 的 制 約 条件が あ るの か今 後の検 討 課 題と してい る。 以 上の ような理由に よ り, こ の線積 分 項に は物理 的 条 件および定 式 化 等に不 明な点が あ ると思わ れる の で,
本 研 究で は (8)式の線 積分 項 を省 略し て計 算を行っ た。 (8) 式の線積 分 項を省略 し た式 を構 造 物 表 面上で法 線 方 向に微 分す る と,
∂
grc
’)一一
岩[
・π σ(・礁
・(Q
) ∂eg
(;
S
・
Q
)d ・]
・
・
・
・
・
…
t−・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
∵・
・
…
(9) の よ うに表す こと がで きる。
ま た,
(9)式の左 辺は,
物 体 表 面の境 界 条 件で あ る か ら, σ に関す る第 2種Fredholm
型の積 分 方 程 式と な るの で これに より吹 出し 分 布 関 数 σ を 決 定す れ ば 速度 ポテン シャ ル は求まる こ とになる。 2.
4 流 体 力お よ び波浪 強制力(9)式 より速 度 ポ テン シャ ルが求ま れば流 体 力お よ び波 浪 強 制 力は
,
変動圧 力 を構造 物 表 面 上で積 分 するこ とに より求め ること がで き る。Bernoulli
の圧力 方 程 式 におい てX
方 向の速 度に,
前 進 速 度 UDを考慮 して微 小 量の2
次 以 上の項を省 略し, 静 水 圧の項 を 除く と 圧力方 程 式は,
P
− 一
職 φ・・仏器
・
………
∵鹽
・
一 一 …
(・・〉 で表せ る。
こ こで, 単 位 速 度の速 度ポ テン シャル φ‘を 導入し, 【adiatiQn ポテ ン シャ ル を用い て,
(10
)式か ら 得られ た圧 力 を構 造 物表面上で積分す ると 流 体 力Ft
∫ は,Fv
− 一
・蕉
{
聯 ・孀略劉
黔
・・ (i,
ノ=
1,
3,
5 )・
・
一
・
曁
…
9・
…
一・
・
・
・
・
…
(11) と表せ る。
Fり は,
j
モー
ドの動 揺に よ る iモー
ド方 向 の流 体 力 を表 す。 (11
)式を時 間項も含めて実 数 部 と虚 数部に分 け る と,
」モー
ドの動 揺に よるi
モー
ド方 向の 付加質量AM
‘, , 造 波 減 衰 力Ntjは,
AM
、,一一
礁
{
畆一
鵬
窪階
d
・ (i,
ノ=
1,
3,
5 )・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(12)Nv
−
・唖蕉
{
甑 ・農籌
}
咎
d
・ (i,
j
;
1,
3,
5)・
・
・
・
・
・
・
・
・
・
・
・
…
t・
・
・
・
…
(13) と表せ る。
こ こ で,iPi
。
,
φtCは,
それぞ れ φ‘の実 数 部と 虚 数 部である。
また,i
方 向に作 用 す る波浪 強制力は,
E,
一一
・蕉
隔
晦磯 隴
・・一
・蕉
{
甜 砺磯
}
髪
・・ (‘=1,
3,
5
)・
・
・
・
・
・
・
・
…
1・
・
・
・
・
・
・
・
・
…
:・
(14
> と表すこと がで き る。 (14 >式第1
項は,
入射波に よ る 圧 力 を 物 体 表 面で積 分 して 求め ら れ る Froude−
Krylovforce
で あり,
第2
項はdiffraction
ボ テ ン シヤ ル に基 づ くdiffraction
forceで ある。
Z,
5Green
関 数の数 値 計 算 精 度の確 認 (4)式に よ り求めたGreen
関 数の計 算 精 度の確 認の ため, 小 林に よっ て与え られ たGreen 関 数の 計 算 結 果 3) と 本 計 算 結 果の比 較を行っ た。 その 結 果をFig.
2
に示 す。
Fig
.
2は,
PointSource
を (xQ,
Yo,
2e ); (0.
0,
0.
O,
−
1,
0
)に置い た と きの一
10<X
,<10, Yp=
0, Zp=
0にお け る
Green
関 数の実 数 部 と虚 数 部 を 比 較 した もの で2
.
1.
Ccs一
1、
一
2.
一
3,
3、
2、
。 c.
1
’
ll
Grttn「
s funCtien G〔β再
O.
1216,
FL=
O.
5,
ξ=
O.
2)Ls
「
Grten’
51凵
目
⊂
[lonG 〔β
=
02739,
’剛
=
05,
ξ篇
02〕Fig
.
2 The real and the lmag章nary part of green・
s functionあ り, β〈1/4の 場 合とβ≧1/4の 場 合 (β
=
ω。 U。/g
) にっ い て検 討 した。 ま た,
こ の 計 算に は 1/ri, 1/r、の 成 分は含 まれて いない。
Fig
.
2
は,
フ ルー
ド数F 。
=o.
5
(F
。=
・
u。IV
]
giiJ
, L :物体の代 表 長さ)β=
0.
2236,
ξ=0.
2 (ξ一
ω:/g)の場 合と,
F。=
O.
5, β=
0.
2739, ξ=
O.
3の計 算 結 果で あ る。
両 計 算は,
周 波 数を固 定し て,
前 進 速 度 U。の変 化につ いて β=
1/4付 近で の Green 関 数の挙 動を確 認し た もの で あ る。 本 計 算 結 果と小 林の計 算 結果と は良く一
致 して お り, β=
1/4で特 異性を持っGreen
関 数の計 算 精 度の確 認が で き た。
同様に, 本 計 算 結 果 とInglis& PTiceの計 算 結 果4) を比 較し両 者が良く一
致して いること を確 認し た。
以 上の検 討よ り,
本計算 に用い たGreen
関 数の計 算 精 度 を確認する こと が で き た。3.
実験概 要 模型実験は,東京大学工学部 船舶航海性 能 試 験 水 槽 (長 さ50m,
幅30 m,
水 深 2.
5m )で実 施し た。 供 試 模 型は,
Fig.
3に示す よ うに想 定 実 機の 1/200 縮尺 模 型を用い た。
模 型の曳 航 時 喫 水は O.
1m , 構 造 物の前進方 向と 入 射 波の相 対 角 度は G度,
実 機 曳 航 速 度は,O,3,6
ノ ッ トであり,
相 似 則に合わせ て フ ルー
ド数を 凡;
0.
O,0.
055,0.
11に換 算して曳 航 中 強 制 動 揺 実 験,
曳 航 中 波 浪強制 力実験を行っ た。 ま た,
コー
ンタイプの よ うな構 造 形 式 を考え る 場 合,
その独特な形 状の た め強 制 動 揺 振一
Fig
.
3
Coロe typed structure mode1tllli 七: 幅変化お よび入射 波 高 変 化による流 体 力および 波 浪 強 制 力へ の影 響 並びに非 線 形 影 響が どの程 度ある の か把 握す る た め, 非 曳 航 時 (前 進 速 度 Oknot )およ び 曳 航 時の実 験 を行い,
流 体 力お よ び波 浪 強 制 力の線 形 範 囲 と非 線 形 影 響を実験によ り明確に し た。
非曳航 時の強制動揺振幅 変化お よび入射 波 高 変 化の実 験は, 日本大 学 理 工学 部海 洋建築工学 科 平面 造 波 水 槽 (長 さ19m ,
幅7m ,
水 深 1.
Om )で実 施 し た。 以 下に そ れ ぞ れの実験の概 要 を示 す。 (1 )曳航 中 強 制 動 揺 実 験お よ び波 浪 強 制 力 実 験 実験は模型の曳航 時 喫 水 を 合わ せ た後,surge ,heave.
,
pitchの3
モー
ドに対し,
ポテンショメー
タの 片 振 幅 を そ れ ぞ れ 1.
5cm,
1.
5cm,
回 転 中 心 角 を O.
03491 radian とし,
規 則 的な強 制 動 揺 周 期を0.
7〜
3.
O sec , また,
過 渡 的な強 制 動 揺 試 験で は,
出 会い周 期 を0.
7−
2.
ユsec と し,
曳 航 速9e
Fn…
O.
0,
0.
055.
O.
11 の 3ケー
ス につ い て実 験を行い そ れ ぞ れロー
ドセ ル で反 力 を, ポ テンショ メー
タで 変 位を計 測 し た。 波 浪 強 制 力 実 験は ロー
ドセル を模型に固定し,
そ の真 横にサー
ボ式 波高計 を設 置して入射 波と波浪強 制 力の位 相 差 を直 接 計 測で き る よ う に し た。 そ して,
曳 航 速 度Fn=
0.
0,
0.
055,
0.
11,
構造物の 前進方向と 入射 波の相 対 角 度が0度とな る よ う 曳航し,
波 周 期0.
7− 3.Osec
の規 則 波お よ び出 会い周me
O.
7〜2.1
sec の過渡水 波を用い て実 験を行い, surge ,heave
, pitchの 3モー
ドの波 浪 強 制 力お よび入射 波 高を計 測し た
。
(2) 非 曳 航 時および曳 航 時の強 制 動 揺 振 幅 変 化お よび 入射 波 高 変 化に よ る流 体 力と波 浪 強 制 力実 験 非 曳 航 時の流体 力 実 験で は
,
模 型の喫 水を合わ せ た後surge
,
heave
モー
ドで はポテ ン シ ョ メー
タの片振 幅を そ れ ぞ れ0.25cm
,0.
5cm
,0.
75 cm,
ユ.
Ocm ,1.
5cm ,
2.
Ocm,3.
O
cm と し た。
また, pitchモー
ドで は回転中 心 角 をそ れ ぞ れ 1.
0 度, 2.
0 度, 3.
0度とし,
規 則 的な 強 制 動 揺 周 期をO.
8〜1.
5
sec と し , ロー
ドセ ル で反 力 を,
ポテンショ メー
タで変位 を 計 測 した。
ま た,
入射 波 高 変 化に よ る波浪強制力実験で は, ロー
ドセ ルを模型に 固定し そ の真 横に サー
ボ 式 波 高 計 を設 置 し て入射濠と波 浪 強制力の 位 相 差 を 直 接 計 測で きるよ うに し た。 そ し て入射 波 高を それ ぞれ 2
、
Ocm , 3.
Ocm,
4.
Ocm,6.
Ocm
と し, 波周期D
.
8〜
1.
5sec の規 則 波を用いて波 浪 強 制 力 を計 測した。
ま た,
曳 航 時の実 験で は,
上 記の同じ実 験条件で前 進 速 度3.
0,
6,
0 knots (F
。;O,
055,
0.
ll) で実 施し た。
4,
結果お よび 考察層
本論 文で は, 曳 航 中の コ
ー
ン形 状 構 造 物に作 用す る流 体力 と 波 浪 強 制 力, 前 進 速 度 影 響お よ び線形 ポテンシャ ル理論の適用性 等につ い て計 算 結 果と実験結果 を 比較し 蓬て考 察 し た。
数 値 計 算で は前 進 速 度を考 慮し た3次元特 異 点 分 布 法に よる解 析 プロ グラム を開 発し た。
解析は構 造 物 表 面を円周 方 向,
深 さ 方 向 お よ び 底 面 部 を 分割する が,
要 素 分 割 数の違いに よ る計算精度の影 響 を確 認する た め要素数を変え て各前 進 速 度に お け るsurge モー
ドの 波 浪 強 制力の 解の収束 性を確 認した。Fig.
4に そ の結 果 を示す。 コー
ン形状構造 物は, 本 研 究で考 慮し ている 前 進 速度では全要素数 を500以 上 分 割 すれ ば安 定し た解が 得ら れ ること が 分 かっ た。
ただし, 本 論 文の計算では,
線 形 理 論の適用性も検討する ため総 計 1440 め要 素に分 割して,
付 加質量,
造波減衰係 数,
波 浪 強 制 力につい て 数 値 計 算を行っ た。 数 値計算に用い た要 素 分 割 をFig.
5 に示す。 な お,
数 値 計 算は (株 〉大林組技 術 研 究 所の スー
パー
コ ンピュー
タSX −
1EA を用いて実 施し た。 また,
本 論 文は線形 範 囲 内にお ける問 題 を対 象と し て い るが、
コー
ン形 状 構造物の よ うに水 線 面と構 造 物 側 壁の なす角 ;が鉛 直と異な り,
こ の ことに起 因 する流 体 力お よび波浪 強 制 力の非 線 形 影 響が考え られる ため強制動揺 振 幅変化 お よび入射 波高変化に よる実 験を行い, コー
ン形 状 構 造 物に作用 す る流 体 力お よ び波 浪 強 制 力の線 形 範 囲.
を明ら かに し線形 理論 解 析 結 果 と比 較・
検 討し た。
(1 > 曳航 時に作 用する流 体 力Fig.
6〜Fig.
8
に前 進 速 度 有り および前準
速 度 無しの場 合の surge
,
heave
, pitch モー
ドに おける付 加 質 量の計 算 結 果と実験結果を示す
。
図中 実 線, ○ 印は前 進 速 度 無し (0
ノッ ト,F
。;O.0
)の場合,
点線,
4
印は前 進 速 度3 ノ ッ ト (Fn=
0.
055)の場 合 お よ び一
点 鎖 線,
口 印は前 進 速 度6 ノッ、
「
ト (Fn =
0.
11 )の 場 合の計 算 結果 と実験結果 を そ れ ぞ れ示す。
o o.
8 64 駄 似 コ 矗 ミ
一
凵 0 2 O Q O 2 4 6 B lO 匹2 14冥
LO軍
MeshFig
.
4 Convergence of wave exciting force results of surge mode 〉 隻 \ = Σ く 】.
5 1.
0 0.
5 0、
0Fig
.
5 Meshes for cone typed mode 且Fig
.
61
『
D 2.
D 3 D 4、
D 5』
D KeLAdded mass coefficients of surge
付 加質量係 数
Fig
.
6よ りsurge モー
ドにお け るFn
=O.
O,0.
055,
0
.
11の実 験結果と計算結果は良い一
致を示す。 前 進 速 度の違い による surge の付加 質量は若 干の前進速度 影 響 が見 受け られる。 特に,
付 加 質量係 数は低 周 波 数 側か ら 付 加 質量の ピー
クとなる K。L=
1.
8近 傍で は,
Fn=
0.
O,
0.
11,
0.
055の順で大き く’
な り,
KeL=
1.
8以 上で は Fn=
0.
055,
0.
ll,
0,
0の順で大き くな る傾 向を示 し た。
ま た,
本 構 造 物が独 特の形 状の た めK
。L =3.
5〜
4.
2近一 123
一
4
.
0 03 D2 〉 ミ 靄 Σ く 1、
O e.
D Fig.
74.
o 3.
o 02 」 〉 隻 \ 踞 Σ く 1、
0 1,
0 2』 3.
0 4.
0 5.
O KeLAdded mass coefficients of heave
0
.
0 9 0 2.
0 3 0 4 0KeL
Fig
.
8 Added mass ceefficients of pitch傍で計 算 結 果が ゼロに近い値 を示 してお り実 験 結 果に も こ の傾 向がうか が わ れ る。 こ の結 果 よりsurge モ
ー
ドの 付 加 質 量の推 定は,
線 形ポテ ン シ ャ ル理 論で か な り良く 評 価で きそ うで あ る。 Fig.
7よ り,
heave
モー
ドに お け る計 算 結 果 と実 験 結 果は良く一
致 して お り,heave
モー
ドに よる付 加 質 量は 前 進速度の違いに よる影 響は あ ま り な い と思わ れ る。
heave
の付 加 質 量の推 定は,
線 形ポ テンシャ ル理 論でか な り評 価でき る。
Fig.
8より,pitch
モー
ドにお け るそれぞ れの前 進 速 度に よる計 算 結 果と実 験 結 果を比 較 する と,
両 者は ほ ぼ 良い一
致 を 示して い る。
特に,Fn=
O.
055の場 合の計 算 結 果は, 良く実 験 結 果を説 明し て い る。
また,
F。=
O,
0 お よびFn=O,11
の場 合,K
,L =1,
8
以 下でFn=O.
11
の 場 合は実 験 結 果よ り計 算 結 果の方 が 大き く,F 。
; O.
0
の 場 合ではこ れ が若干下ま わ る傾 向を示し てい る。
今回の 実 験 範 囲 全 般で は, 前 進 速 度 影 響 が 見ら れ,Fn=0,
055
, Fn=
0,
11,
Fn=
O.
Oの順で pitch の付 加 質 量は小さ く な る こと がわ かっ た。 e > 咲 \ = Z s > ミ 羇 Z 1.
5 1 0.
5 0.
D I.
D 2.
D 3.
0 4.
0 5、
O KeLFig
.
g Wave damping coefficients of surge 1.
5.
o0
.
5O
、
O LO 20 30 4.
0 5』 KeLFig
.
10 Wave damping coefflcients of heave造 波減衰係数
Figs.
9−
11に前進 速 度有りお よ び前進速度無 し の場合の surge
,
heave,
pitchモー
ドにお け る造 波 減 衰 係 数の数値 計算 結果と実験結果を示す。 図中 各 線お よ び記 号 は
,
上 記 付 加 質 量 係 数の 各線お よび 記 号 説 明と同 様であ る。
Fig.
9に surge モー
ドにお ける造 波 減 衰 係 数の計 算結 果と実 験 結 果を示 す。 計 算 結 果は実 験 結 果 を よく説 明し て いる といえる。 surge の造 波 減衰係 数に は前 進 速 度 影 響が顕著に現れており,K
。L
=
3.
0以 下で はFn
=
0.
11,
0.
0,
0,
055の 順に造 波 減 衰 係 数が大き くな り,
K。L・
=
3.
0 以上で は これ が逆 転する現象を示し て い る。 こ の傾 向は実 験 結 果に も 同様に現れて い る。
Fig.
10にheave
モー
ドにお ける計 算 結 果 と 実 験 結 果 を示す。
計 算 結 果と実 験 結 果は ほぼ同様な傾 向 を示し て お り, 顕 著な前 進 速 度 影 響が ない こ と が わ か る。
特に,
計 算 結 果の Fn=
O.
0で はK 。
L
・
=
2.
5で波 無し点 が 現れ て い る。
実 験 結 果で は同 じ くF。
=
0,
11の結果に K。
L=
2.
5で波 無し点が現れて いる。智
〉 隻 \ 旨 Z 3.
0 2、
o 1.
O 。窄
…ぎ
塋1
}
il
含
≡i
ヨロじ T △ △
合
ダ
ミ
・
/
’ ロ /d
今/ / /’
1’
!_一
〆 ! A ノ!
!ノ
o O OG 1』 2.
0 3.
0 40 KeLFig
.
11 Wave damping coefficients Qf pitch o.
8 50 40 』 矗 ミ.
凵 o.
2 副‘
龠
曹
曹
叩 O ム ロ E【
5
1 FOO1
吾
o.
t
σ 二二b
!こミノ
,
〆ク! ! 0 △ o 00 Fig.
12
1.
D 20 30 ヰ.
0 5.
O K,LWave exciting forcE coefficiemts of surge
Fig
.
11に pitchモー
ドに お け る計 算結果 と実 験 結果 を示す。 計 算 結 果と実 験 結 果を 比較する と, 両者は定 性 的に 同様な傾 向を示し てい る が, 計 算 結 果は実 験 結 果 よ り も下ま わっ た値 を示し た。
こ の理 由につ い て は粘 性 減 衰およ び自由表 面の か く乱による非 線 形 影 響 等の た め と 思わ れ る。
計 算 結 果 を比 較して み ると,
pitchモー
ドの 造 波 減 衰 係 数に は顕 著な前 進 速 度 影 響が現れてお り,
K
。L
= 2.
5
以下では前 進 速 度の速い順に値が大き く な る。 (2 ) 曳航時に作用 す る波 浪 強 制 力 Figs.
12〜
14 に前 進 速 度有り お よ び前進 速度無 しの場合の surge
,
heave,
pitchモー
ドに お け る波 浪 強 制 力の計算結果と実験結果を 示す
。
図 中の各 線お よ び記 号は,
上記の説 明と 同 じである。
Fig.
12に surge モー
ドに おける計 算 結 果と実 験 結 果 を示 す。
計 算 結 果と実 験 結 果は,
ほ ぼ良い一
致を示し て い る。
K。L=
2.
5の ピー
ク近 傍では,
Fn羃
O.
0,
Fn
=O,
055,
F
。1O.
11の前 進 速 度の遅い 順に大き く なっ て い る。
ま た,
計算結果に も 同様な傾向が見ら れ る。
特に,
顕 著な 0.
8 0.
5 40 畠 ミ ぼ D.
2。
守
咢
1塋li
謡
含
≡ D △ Oレ
必
叶 丶丶
丶 O 丶 」 丶 丶LN
.
、 、 \ 、 、 0.
0 卩 1 コじMT Fig.
13o.
B 0.
6 脳 」 島 隻 \ n 凵 r勹
o.
2 1、
0 2.
0 3.
0 4.
D 5.
O KeLWave exciting force coefficients of heave
。
呈
nEぎ
隻llii
昌
≡!f2 ?・
寸
胃 な蝉
〆、
o一
彡
1’
t甼
M 0、
0 1.
0 2.
D 3.
0 4.
0.
50 Kε
LFig
.
14 Wave exciting、
foTce coefflcients of pitch前 進 速 度 影 響は見られない
。
Fig,
13 に heave モー
ド に お け る計 算 結果と実 験結果 を示 す。 計 算 結 果と実 験 結 果はほぼ良い一
致 を示 して い る。K 。
L ;
2.
5以 下で は,
Fn =0.
0,
Fn
・
O.
11,
Fn =
O.
.
055 の順で 大き く な る が,K
。L =2.
5以上で はFn=
0.
11,
F。!
・
O.
055,
F
。=
0.
0の順に大き く な る。 また,
LK 。
L #2.
5
近 傍で波浪強制力は ゼロ に近い値を示す。
こ れ はheave
の造 波 減 衰 係 数の波 無し周 波 数に対 応 し て いる。 また,
顕 著 な前 進 速 度 影 響は見 られない。
Fig.
14に pitchモー
ドに おける計 算 結 果と実 験 結 果 を 示す。
計 算 結 果と実 験 結 果 を 比 較する と, 両 者は ほぼ 同 様な傾 向 を示 す。K
。L =
2.
1以 下ではF
。’
=
O.
11,Fn
=
0.
0,
F。
=
O.
055の順で 大き く な り,
K。
L=.
2.
1以 上 で は Fn=
O.
0,
Fn=
O.
11,
Fn=
O;55の順に大き く な る。 実 験 値も同様な傾向を 示 す。
ま た,
顕 著な前進 速度 影 響 は見られ ない。
(3
) 曳 航 時 お よ び 非 曳 航 時の強 制 動 揺 振 幅 変 化によ る 流体力Figs.15−
20に 曳 航 時お よび 非 曳 航 時の 場 合の一
125
一
・ ×
8.
m ロ 0.
▼
OO、
m。
OmL 丶 ≡ 」 OO,
N OO一
、
』門
三 : L.
[
Oロ
.
面 Fig,
15 OlX00.
い
OO.
Ψ
O ロ 』ミ
,
n ≡ 」 囗 O.
N OO一
、
巨
。
一
閃
−
二 二 」ロ
O.
O x10齟
■
5.
00 10.
00 し/BTotal
and linear hydrodynamic forces of sttrge(Kβ
=
=
O,
9019 and Fn=
0.
0,
0.
055 and O.
11)Fig
.
16×
10’
■
5.
00 10.
OO L/BTotal and linear hydrodynamic forces of suTge
(KA= ユ
.
4585 and Fn=
O.
O,
0.
055 and O.
11) surge,
heave,
pitchモー
ドにお け る 流体力の一
次成 分と二次 成分の 和の最 大 値 (以 下
Tota1
の流 体 力 )お よび一
次 成分の流体力 (以 下Linear の流 体 力 )の結 果 を示 す。 図中○印, △ 印および×印, ◇印は,
曳 航 時 (F
。 =0.
055 , 0.
11)の 無 次 元 周 波 数 K8=
0.
9019,
Ks; 1.
4585 (KB=
ωヲg ・
Bl 波 数h
:h =
ω2/g ,
B :B;
L/2 ) お よび‘ ● 印, ▲ 印は,
非 曳 航 時のKe=
O.
9019,
KB=
1.
4585 のTotal
と Linearの流 体 力の実 験 値 をそれ ぞれ 示す。 なお, 図 中 実 線, 破 線お よ び一
線鎖 線は , Fn;O.
055,0.11
お よ び 0.
0の同 無 次 元 周 波 数のLinear
の 流体力の計算結果 をそれぞれ示す。
surge モー
ドの流体 力 Fig.
15,
16に,Ks=O.
9019
,K
.=
1.
4585 お よ び強 制動 揺 片 振 幅ユ
.
Ocm , 1.
5cm , 2.
Ocm ,
3.
Ocm
の場合のsurge モ
ー
ドの流 体 力 を示す。
両図 よ りFn=
0.
0の場 合,
両 周 波 数 と も強 制 動 揺 振 幅が大き く な るに従っ てTQtal
の流 体 力お よ び Linear の流 体 力は大き く な り, 両 者は ほ ぼ同 じ値で ある。
強 制 動 揺 片 振 幅が 3cm の場 合, 両 周 波 数 と もTotalとLinear の流 体 力に若 干の差 を 生じ 高 次の流 体 力 成 分の影 響が見られ る。
両 周 波 数の計 算 結 果とLinear
の実 験 結 果 を比 較し ても, 両者は よく一
致門
ロ
α
氏
丶=
一
」磊
一
←
t
x 亅o’
1 ×10ー
コ
・
00 10.
00 し/BFig
.
17 Tota正and [inear hydrodynamic forces of heave(κ8
=
0.
9019 a既d Fπ
=
0.
0,
0.
055 and O.
11) xlO’
o『
00円
b 研 」 \一
≡ 」 O囗
閏
O冖
一
、
閂
嚠
計 7 ≡」
0 Fig.
18 x 凪0’
3
5.
00 ]0,
D口 L〆BTDtai and linear hydTodynamic forces of heave
(Ks=1
.
4585 and Fn= O.
0,
0.
055 and O,
11>して おり強 制 動 揺の振幅 影 響は ほと んどない こと が分か る。 また
,
Fn=
0.
055,
0.
11の場 合,
両 図よ り,
強 制 動 揺 片 振幅の 変化にか か わ ら ず,
両 周 波 数ともTotal とLinear
の流体 力に差を生じず,
Total
お よ びLinear
の流体力の強 制 動 揺 振 幅変化に よ る影 響は見ら れ ない。 ま た. 二次 成分の流体力の影 響も ほ と ん ど ない こと が分か る
。Linear
の計算 結果と実験 結果の一
致も よ く,
高次 の流体力 成分 によ る非線形影 響は あ ま り生じ ない こと が 分かっ た。 以上 よ り,
surge モー
ドの流 体 力は,
曳 航 時 お よび非 曳 航時で も強 制 動 揺 片振幅 lcm,
1.
5cm,
2 cm まで は高 次の流体力に よ る非線形 影 響は な く,
強 制 動 揺の振 幅影 響 も少 ない こと を 確 認 した。
heave
モー
ドの流 体 力 Fig.
ユ7,
18に K.・
=
O.
9019 と 1.
4585,
強 制 動 揺 片『
振 幅0.
5cm,
L O cm,
1.
5cm,
2.
O cm の場 合の heave モー
ドの流体力を示す。 両 図よ りF
。=
O.
0
の場 合,
両 周 波 数の結果は強制動 揺 振幅が大き く な る とTotal
お よびLinear
の流体力は増加し, 強 制動揺 片 振幅3cm
で高次 の流 体 力 成 分に よ る非 線 形 影 響 が若 干 見ら れ るだけ で, 後は 同 じ値と な る。
ま た,
計算結 果と実 験 結果 を比較す る と,
両周波 数と も強制動揺 片振幅 0.
5cm の場 合に,
・ 艮 OO
.
α
X OO.
凹
b−
丶
[
二
巴一
.
腕
.
『 ≡一
」▼
O氏
\
一
=
「
」一
.
−
u
一
警
}
「
;
」 OO、
O Fig.
19 ゜ 且 OO.
四 X ×10.
t
L,
00 2.
00 ヨ.
00 L/BTotal and [in巳ar hydrodynamic fQrces Qf pitc卜
(Ks
=
0.
9019 and F町
=
0.
O,
0.
055「
and O.
11>]
.
o囗し/B
2
.
OO翌品
置
Fig
.
20 Tota1 and li皿ear hydrodynamic正orces of pitch(Ke=1
.
4585 and Fn=
=O.
0,
0.
055 and O.
11) 計算 値よりも小さ な値を示して い る がheave
モー
ドの 流 体力に強 制動 揺の振 幅 影響は ほと ん どな い こと が分か る。 また,F
。=
o.
055,
0、
11 の場 合,
Fig.
18 よりt Fn=
O, 055 で 強 制 動 揺 片 振 幅 O.
5cm,
2.
Ocm の場 合の実 験 結 果が小さ く なっ て いる が,
強 制動揺 片 振 幅1.
Ocm,
1,
5cm
の場 合は計 算 結 果とも良く一
致して お り,
Total
とLinear
の流体力に も差が なく強 制 動 揺の振 幅 変化に よ るheave
モー
ドの流 体 力へ の影 響は少ない と 思 われ る。
Fn=
O.
llの場 合, Totalお よ び Linearの流 体 力は ほ とん ど 変わ らずheave
モー
ドの流 体 力の非 線 形 影 響 は見られ ない。
ま た,
両図と・
も計 算 結 果と実 験 結 果も良 く一
致 して い る の で, 前 進 速 度 を 有する場 合の heave の流 体 力には非 線 形影響 が ない ことが分か っ た。 pitch モー
ドの流体 力Fig.
19,
20 に・
K
,=o.
9019 と1.
4585, 回転 中心角 度1,2,3
度の場 合の pitchモー
ドの流 体 力 を示す。
両 図 よ りFn・
O.
O
の 場合,
他の モー
ド と 同様に,
両周波 数 と も強 制動揺振 幅が大き く な る と流 体 力は増 加し,Tot.
ai・
とLinear
の流 体 力は同じ値 を示す。 また,
両 周 波 数 と も計 算 結 果と実 験 結 果を比 較すると,
pitchモー
ドの 流体力の 強 制 動 揺の 振 幅 影 響 は ほ と ん ど ない こ とが分か OXOOO
.
旧
門
ロ ロ.
9一
ロ O.
の
n
」 ω 匹\
一
二
」一
.
匚
凾
‘
〔
縄 L + = 」〕
x 】o’
1
000 100 200 3、
DO H/xFig
.
21Total
alld Iinear wave exciting 正orces of s皿rge(K
.
=
0,
9019 and Fn=
0.
0,
0.
055 and O.
11) OXロ
0O 門ロ
O.
n[
n
」 m ミ0
口
.
0【
ロ
ロ
.
頃一
=」
.
「
。
回
〔
記 止舮
=
」一
旧・
聖、
1565 Fn 歐 卩.
馳 [ ooTotal ● Lime証
▲一.
一
0055 丁帆aIoLユ
nE旺
△一
01崖 b臨1xUn 閑 「 ◇一
一
一
!!
/ /
麺
/’
/q /卩
ノ
〆
ノ
//
・ ツ「
、
彡
ル
ア
一
男
/ノ
ノ
ノ
’
!
〃
.
〃
。、
。。L 。。
2
.
。。3
.
。。,
、
。。5
.
茜イ ゜.
【
HAFig
.
22 Total and tinear wave exciting forces of surge(K,
=
1.
4585 and Fn=
=
O.
0,
0,
055 and O.
11) る。
また,
高次の流 体 力 成 分に よ る非線形影 響もほ と ん どな い こと が分か る。F
。=
0.
055,
0.
llの場 合,
両図 よ りTotal
と Linearの流 体 力に は差が見ら1 れず,
計 算 結 果と実 験 結 果 も良く一
致してお り,
前 進 速度を有す る場 合の pitchモー
ドの流体 力には非 線 形影 響 が 見 ら れ ない層
こと が分かっ た。 (4 )曳 航 時お よび非 曳 航 時の入射 波 高 変 化による波 浪 強 制 力Fig.
2ト 26 に曳航 時お よ び非 曳航時の場 合の surge ,heave,
pitchモー
ドに お・
ける波 浪 強 制 力の一
次 成 分と 二 次 成分の和の最 大 値 (以下Total
の波 浪 強 制 力 )お よ び一
次 成 分の 波 浪 強 制 力.
(以 下Linearの波 浪 強 制 力 ) の結 果 を示 す。
図 中 各 種 印および各 種 線は流 体 力で説 明 し た もの と 同 じ で ある。
surge モー
ドの波 浪 強制力Fig
.
21,22
に,
K
β=o
,.
9019
と1.
4585,
入 射波 高LO
cm,
2.
Ocm,
3.
O cm,
4.
O cm の 場合の surge モー
ドの 波 浪 強 制 力を示す。 両図 よ り, F。; O.
Oの場合, 両 周 波 数 と も波 形 勾配 が 大き くな る に従っ てTota1
お よ びLinear
の波 浪 強 制 力は増 加し両 者は ほ ぼ同じ値 を 示す。
両 周 波 数と もTotal とLinear の波 浪 強 制 力に は,
高次一 127 一
×
聰0’
58]
.
国門
」2
丶 ロ O.
の
,
90.
V
ロ O =r
し・
匹
霊冖
昌」
+一
F
」〕 .
N 0・
OO「
Fig.
23 OXロ
O Oロ
.
囗
0ロ
ロn
」 n 」 丶一
= 」 〇 〇、
V gO冨
oσ
[
胃 」+
=
」
一
N x10.
’
1.
DO 2.
00 ヨ.
00 HハTotal and lnear wave exciting forces of heave
(Ke
=
O.
9019 and Fn=O.
0,
0.
055 a皿d O,
11) 0.
001 GO2 0DH /x3 DO4,
ao xlO』
:
5.
GOFig
.
24 Total and 駈ロear wave exciting forces of heave(K,
≡
1.
4585 and Fn=
0.
0,
0.
055 and O.
11) もL °『
t:
00,
」 m ミ.
守
一
= 」 0 ロ,
N〔
窩 L + 二」
一
xlO”
O、
00 1.
00 2、
OD ヨ.
00 H/λFig
.
25
Total a皿d
linear
wave exciting {orces o 正piIch(κo
耳
0.
9019 a琢d F昂三
〇.
0,
0.
055 a皿d O.
ll) xlo幽
呂
旨口
OO.
田
ロ
O.
笛
OO.
噂 00.
Jo 乢 丶 =.
」鬩
。
ロ
[
窓 軋 +一
一
巴.
四
0.
001
.
002
.
。0300
亀
.
0。5
も
1
『:
H〆λFig
.
26 Total and 旦inear wave exciting forces Df pitch(κ