1
論 文I
UDC :627
.
24 :624.
074.
4日本建 築 学 会 構 造 系 論 文 報 告 集 第428号
・
1991年10月Joumal of Strucし
.
Consしr、
Engng,
AIJ,
No.
428,
0cし,
1991海洋
に
お け
る
浮
遊式 弾性
円
筒曲面構造物
の
振動性状
’
流
体 調和
表
面 波 を受
け る 場合
RESPONSE
CHARACTERISTICS
OF
FLOATING
TYPE
OF
ELASTIC
CYLINDRICAL
STRUCTURE
IN
OCEAN
SPACE
,
Under
harmonic
surface wave of water福 住 忠 裕
* , 日 下部 馨
* * ,野 添 久 視
* * * ,堯 天 義 久
* * **TadUh
加FUKUSUMJ
,.
ぬo厂% 丿ζこfSAKABE
,Hisashi
NOZOE
, and }「oshihisa
GYOTEN
Cylindrical
she 旦1s
are usedin
oc 尹an space as resevoirs andfloats
to support platforms andbhildings
.
Dynamic
response characteristics of large scaled floating type of elastic cylindricalstructure
has
not yetbeen
dlarified.
In
this paper,
a cylindrical shellstrucuture which consists of an elastic shell wall,
bottom
plate and mooring springs on its circumference is studied , and not only motiQn as a whole sturcture but also distribution of stresses and displacements of the elasticplate and shell are investigated
.
from
止e elastic structural poi1ゴt of view,
.
Solutions
of motion ofthe elastic structure system and the f}uid are obtained as theoretical solution
』
consistently.
Dyna
皿ic【esponses of the structure anCl thefluid
to theharmonic
water wave aredetermined
andthe effect of the thickness of the shell and plate and the st 五ffness of the mooring spring o皿 the re
−
sponses aredescribed
.
、
Kegwords :aylin‘irical shell
,
floating
st厂a砌 re,
water zvave respanse,
(〜ffshore
structure「
円 筒 シェ ル,
浮 遊 式構 造物,
.
波浪 応答,
海 洋構 造 物 L 序 海 域に は,
従 来 漁 業や 工業にか か わ る多.
くの施 設が建 造さ れて き た。
近 年,
海洋空 間は さ ら に広 範 多岐に わ た る利 用を目指して開 発さ れ る方 向に あ り,
海上 空 港や海 上都 市 建 設が実現す る時期を 迎え て い る。
こ の種の海 洋 構造物の多く は,
必然 的に大 規 模 構 造 物 となる こと と予 想さ れ,
その動 力 学 的 特 性の解 明は急 務であ る。
海 洋 構 造 物一
流 体系に関しては従来 より多くの研究があ る もの の,
大 抵は構 造 物を剛 体と して扱う もの であ る。 し か し 構 造 物が大 規 模で薄 肉とな る場 合に は弾 性 変形 が 顕 著に.
生じる可 能 性も あり,
その弾 性変形 を考慮した解析を行 う 必要が あ る。
剛構造物の場合は,
波のマ クロ的 外 力 と して の波 力 (構 造 物 全 面に わた る波圧 の積分値 )と の関 係か ら応 答が決まる一
方,
弾 性 構 造 物の場 合,
構 造 物 各 部 位に分 布 する ミ クロ的 外 力である波 圧レ ベ ル で の解析 を 行 う必 要が ある。
固 定 式 弾 性 構 造物に関しては 登坂IJ,
西 村z )が流 体 波 動と3次 元弾性
体 お よ び シェ ル との連 成 問 題の定 式 化を行い,
そ の取り扱いにつ い て述べてい る。
各 種 固 定 式 弾 性 構 造 物を取り扱う研 究と して は,
濱 本,
田中4}¶
6レおよび筆 者7}に よるもの と して,
海 底に固 定さ れ た 円筒シェ ル の動 的 応 答に関す る解 析 的研 究がある。 また,
遠 藤ほ か8) に よっ て は同シェ ルが流 体の波を受け る場 合につ いて, 弾性変形に伴う波圧分 布に関して の実 験 的 研 究 が行わ れ てい る。
さ ら に, 西村ほ か pv °] に よっ て は,
同じ く海 底固定の弾 性球殻お よ び円錐シェ ル を対 象と し た 固有値問 題,
波浪・
地 震 応 答に関 する解析 的・
実 験 的 研 究が行わ れ ている。
』
固 定 式構造 物に比べ て浮 遊 式 構 造 物の特 徴とし て は , 深 水域で の設 置が可 能なこ と
,
周 辺諒
域 環 境に及ぼ す影 響が 小さい ことが あげら れ る。
松井S〕 は弾 性 構 造 物 を扱 う もの で はないが,
形 状お よ び形 式 が 任 意の構造物を扱 うことの できるハ イブ リッ ド型 解 析 法に よる定式化を示 し た。 各 種 浮 遊 式 弾 性 構 造 物の う ち2次 元 平 板 を扱う も の と して,Wenll
)は 底面に 分布 係 留 索を有 する矩 形 平 板 につ いて,
ま た,
田中ほ か12 }は円 形 平 板につ い て固 有 振 動お よび強制振 動に関する研 究を行っ てい る。 最 もシ ン プル で基本 的な3次 元 構 造 物と しては 円柱 体が あ げ ら れ る。
円柱 体は内 部 空間 を直接利用, あ るい は上部構造物 * 神戸大 学 工 学部建築学科 牌 神 戸大 学工学 部 建 築 学 科 * * * 神戸大学大学院 自然 科 学 研 究 科 * * *一 神 戸大 学 名 誉 教 授DepL of Architecture
,
Faculty Qf Engng,
Kobe UniversityDept
.
of Architecture,
Faculty of Engng,
Kobe UniversityGraduate Schoo)of Science and Technolbgy
、
Kobe UniversityEmeri
しusProf.
KobeUniversity
の 支 持 機 構と して利 用さ れ る。
Garretti3
}は有 限 水 深に おけ る剛な円柱 状 浮 体の強 制 波 力 を求め,
井 島ほ か14)は 同 浮 体の運 動につ いて も解 析 し, さらに合田ほ か15 ]は同 浮 体の波 力と運 動に関す る実 験 も含め た研 究 を 行っ た。
円柱状 浮体を弾性体と し て扱っ た研究は筆 者ら を除い て はこれ まで他には身 受け ら れ ない。筆 者ら は,
円筒シェ ル側壁,
円 形底板,
係留索か ら構成さ れ る浮遊式弾性円 筒 シェ ル構造物の 調和 波 応 答に関 す る研 究1?)−
19)ig
行っ て きてお り, 本 論文は円筒シェ ル, 円板, 流体の基 礎 微 分 方 程 式に, 有 限フー
リェ 変 換 法, 有 限ハ ンケル変 換 法 を適 用し理 論 解 を 得るもの で あり,
さ らに, 数 値 解 析に よっ て は構 造 物 構 成 要 素である シェ ル側 壁,
底 板,
係 留 索 剛 性が,
構 造 物の全 体 的な運 動と波 力 並びに構 造 物 各 部で の変 形,
応 力および波 圧応答 分布に及ぼ す影 響につ い て調べ,
弾 性円筒 構造物と 流体との動的 相互作用の 解 明を図ろ う と す る もの であ る。2.
弾 性 円 筒 構 造 物の基 礎 式お よ び解の誘 導 Fig.
1に示 すと おり,
絶 対 座 標には構 造 物 静止時の重 心 点G
を原 点に し たx
(鉛直上向 ),y ,
z
直交座標系 を用い る。一
様 水 深 L の海 域に お い て,
構 造 物 自重下 で の吃 水 をL
,とする円 筒シェ ル構 造 物が下 端 外 周に沿 う線 形ば ねによ り係 留さ れてい る もの と す る。
構 造物は 左 方 向 か ら 円 振 動 数 ω のHarmonic
な流 体入射 波を 受 けて運 動する場 合, 構 造 物の側 壁 弾 性 円 筒シェ ル,
円形 底 板, 係 留 索および構 造 物 周 辺 流 体の理 論 解の誘 導につ い て述べ る。
2.
1 円 筒シェ ル の解 円 筒シェ ル (半 径 a,
高さH ,
板 厚hs,
面密 度ρshs,
ヤ ング係 数 Es,
ボア ソ ン比 Vs)は等 方・
等質 弾性 体で 薄 肉で ある と し,
シェ ル理論 式に は,
Kirchhoff−
Love の仮 定に基づ くDonnel
型の式を適 用す る。 こ の場 合,
シェ ル下 端を原 点とする円 筒 座 標 x (鉛直上向 ),
θ,z 各 方 向の シェ ル の運 動 方 程 式は次式で あ る。 乢.
エ十Nbl,
e/α一
Pshs
U,
tt=
0 ハle,
θ/α十 ハlxe,
エ
ー
ρs/ZsIく趾=
0 ヨ Vθ/α十Mx、
xエー
21鴫儲 θ/α 十ルfe
,
ee/α2一
ρ$んs四ε厂Ps=0
た だ し,
……・
…・
…・
(1) 君π= ∂F
/∂x,F,
e= ∂F
/∂θ,1
「尸 ∂F
/∂t
・
………・
・
…・
・
…・
…・
…・
……
(2 )Fig
.
1 Analyzed model of structure・
fluid system一
108
一
−.
He Mex ’ ,ク
”
か■
一一嚊_
= = 「帽
騨
HxMxe ロ、。・
1
『
、
丶一一
6罵 漏Fig
.
2
Stresses acting on the element oI the sheUシェ ル の x, θ, 2 方 向の絶対変位
U ,y ,
W
は,
剛 体 (一
様 直線 的 )変位および相 対 (非 直線 的 )変位の和で表し,
次 式の と お り変 数 分 離 形でお く :U
(工,
θ,
の= Σ [U(コC)十 εlnVa ]COS nθ・
e剛 rt=
oぷ
y(x,
θ,
t)=
Σlv
(x)+ε1π[一
γ+(コr−
x。)姻i
n=
1・
sin nθ・
eiWtの
罪(x,
θ,t
)=
Σ 傾エ)一
ε in卜 γ 十(x−
Xo>Ψ]} η=
0・
COS nθ・
e星ω¢…・
……・
・
・
・
……・
一 一 …・
…
(3
) こ こ に.
1 (for n=
1)ε・; 。 (
f
。 , 冗≠1
ガ
… ’
… … ’
’
”
(4) γとv
は水 平 方 向 剛 体 変 位と剛 体回転 角であり,
x。は 底板 中心線か ら底 板一
シェ ル構造物全体の重心G
点まで の 距 離で ある。 本 解 析で は鉛 直 方 向 変 位 をCosine
型 フー
リエ 変 換す る た めに,
鉛 直方 向 剛 体 変 位 を 表す一
様 変 位の項を u(x)中に含んで い る。 シェ ル面 外 方 向に作 用 する流 体 圧 力Ps
(流 体は圧 縮 力 を 正 と する)お よ び シェ ル の応 力と してNx ,
Mx
(Fig.
2参 照 )は変 位 と 同 様, 周 方 向 波 数の重ね合わ せ と して次 式の形で表す :co Ps(即
,
θ,
t)=
Σ⊃Ps(」)。
COS n θ・
et4dt n±
Oco
Nx
(認,
θ,
t)= Σ nx(C)。
cos nθ・
etωt n=
om
Mx
(X,
θ,
t
)= Σ Mx 〈x}・
Cos nθ・
ettet η=
o…
(5} (1)式に現れ る応 力の変 位 表 示,
およびせ ん断 力 と置 換せ ん断力は下式,
N 。
= 〔B
、/α)(αU
=
+v,IV,
,+レsW )Ne =
(Bs
/a)(aVsU,
エ十V,
e 十W
)Nx
θ=
(Bs
/2ω (1−
Vs)〔U,
e十a蕨=
M苫=一
(05/α 2 )(a!Wx
エ
十 Vsl概θ θ> Mθ=一
く1)s/αt)(四θθ十レ3α2四苫工) Mxθ=
(Ds/a)(1−
Vs)VV=
θQ
、
c=
脇,
。
−
M。
。,
e/αQe
=Me,
θ/α一Mxe.
xVx
=Qx
− Mxe1
θ/αVe
= N =e−
Mxe/α一 ……
(6)である。 シェ ル の材 料 減 衰は履 歴 型で減 衰 比を
ds
と し,
上式 中の ひずみ剛度 Bs お よ び曲げ剛 度D。 は次式で表 す;Bs=
Eshs
(1
十2ids
)/(1一
μ舊),
DS=B
ε峠 ノ12
・
一 …・
・
………・
……・
・
……・
・
(7 ) こ こに,Es ,
Vs は シェ ル の ヤ ング係数,
ポァソ ン比で あ る。 偏微分方 程 式であ る運 動 方 程 式 (1 )式の最 初の式に 有 限フー
リエCosine
変換 Cll
を , 同 じく第 2式, 第3
式に有限 フー
リエSine
変 換Sl
トを 施す と , m=
1,
2…
の 時は下 式の とお り シェ ル 変 位の像関 数 uc(am ),
vs (am ),
w9 (αm )に関 す る 3元 連 立 代 数 方程 式が得ら れ, また,
m=
・
Oの 時は像関数uc(α。)が単 独に得ら れる ;ー
33
3 し
ま
じ
444
2.
2 しヨ ハ 乃 ー
M
渦Y
S
ー
8100ー
= oo 轟 22BBO uc(α m) vs(αin) ωs(αm) 00
BSi
・
({λ!lト← 1} Mlx /)十 εmiqA 十lqpl
……一 ・
・
(8 ) こ こ に, {鮒 ,lx
,1
は シェ ル上 下 端 境 界の 変位と応 力 ベ ク トル で あ り本 解の 未定
定数で あ る。・
lqil
,
{昭 は 構造物剛体変位に よ る慣 性 力,
流 体 圧ベ ク トル であ り,
これ ら は下 式で表 示さ れ る。1
疋1ト1
η。
(0),
v(0),
ω (0),
凋0
)1
「・
・
・
…
(9
)lx21
=
ln
エ(H ), v(H), w (H), mx (H)} 「Claurl
lqt
}=− b
ρshs (itO
)2Sl
一
γ十(x−
x。}邸l
Slγ一
一
(x一
コじo)『}lq
鋪=
io
,
0,一
かP
壼(α皿)} 7・
・
・
・
・
・
・
・
・
………・
…
(10
> さ らに,
マ ト リックスの要素は次式で表される :]
An
= a’a轟+dnt− bp
、h
、ω ’,
A、
,=−
ena αm ,A
、S− −
ysaam , A 、、=
n2 +dα 2 轟一
bρ、h、 ω 2 ,A
、、1・
n,
’
A
、、=
1+ cげ +α ’ ・k
)2− bp
、h、b、 2,
.
B .
=− b,
B、
2』−
dαn,
B、、=
da2a。
,
B33=
Cαzα m[aZα揚十(2−
Vs)n2],
B34=b
α m………一 ・
・
一 ……・
…
(ll) な.
お,
α皿
= m π/H
,b
= α t /Bs
, c= ん§/12α 2,
d =
(1−
Vs)/2, e=
(1 十 Vs)/2……・
…・
一 ・
…
(12) 逆変換に よっ て原関 数u は次 式で書き表さ・
れ る ;・ω
一
盈芳
ガ(… )… s ・,…一 …………
(13) こ こ に,
1 (for m = 0> ε皿=
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(14)2
(for
肌 ≠0
) 同 じ く,
原 関 数 v,
w は, フー
リエ 変 換 適 用 領 域 全 域 で の評 価 を 可能と す る一
般 表 現an〕 を用い て表す :’
・(x)一朞
葱
{
vs(・・)+士
[(−
1)・ v(H
)−
v(o
}]}
・
s・・・… +晋
・(・)+ 〃静
(・)ωω
一
葺
葱{
げ (a・
)+毒
[(一
・}・w (H
)−
w(o
)]}
・
・i
・・mX +lil
’
・(・)+ H静
(・)・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(15) 2。
2 円形 底 板の解 タン クの底 板 (半 径 α,
板 厚h
,,
面 密 度 ρphp )は,
等 方・
等 質の弾 性円板と する。
円板は面 内に剛セあり面 内変形は無い もの とし,’
せ ん断 力によ る変 形は考 慮し な い。
この場 合,
円板の運 動 方 程式
は次 式である :D・
(
∂2 ∂ ∂1 ∂r2+ 。∂r + 。・ ∂θ・)
(
鼻
・者
・諾
,、)
u・
∂!Up
…
r・
…
r・
・
・
・
・
・
・
・
・
…
『
・
…
7『
…
− 7
(16)=P
,一
ρ尸h
尸 ∂t2 円板の x,
z 方 向変 位U
。,
Wp を次 式で表す ;ゆ
防(r,θ,t)=
Σ [Up(r)十εlnVr ]COS nθ・
eSWt n=
oWK θ
,
・t)一
、、。[r+ X。 V ]C・S θ・
e ‘te’
t・
曁
・
・
『
9・
・
・
・
・
…
辱
一・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(17) 底 板 面 外 垂 直方 向に作 用 する流 体 圧 力Pp
お よ び底 板 応 力の例とし て板の曲げモー
メ ン トM.p は変位と同様, Pp(r,
θ,
t)=
Σ Pp(r)・
cos nθ・
eiWt n!
o…
(18)M7P
(7,
θ,
t}= =Σ]Mrp (7)・
COS 7Lθ・
eiblt n¥
e と表す。
こ こ に応 カー
変 位 関係はM ,p(r)
一
“
・
Dpl AuK r)一
(≒
ンP}(
訴
一
穿
)
u,Cr
)}
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(19) こ こ に,
A・・(r)一
(
dt
d
n2dr
・+ 7d7一
π戸
)
u・
(r)・
一 一
(・・) で ある。 底 板の材 料 減 衰も シェ ル と同嘩
履 歴 型で減 衰 比 を dpとし て,
曲げ 剛度 D,を次式で表す。
D,;E
,hXl
十2id
,)/[12(1−、
レ多)]…・
・
………
(21
) こ こ に,E
,,
レp は底板の ヤング系数,
ボアソ ン比であ る。 有 限Hankel 変 換,
晒 )}・
∬
・・
Jn
(Slr
)・
・dr
・F
・(SD
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
9・
・
…
(22) を 運動 方 程 式 (16)式に適 用すれば,
DpSEaJn−
1(s
‘α) [− SIUP
{a)+Au尸(a) uF’ (Si)=
DpSt+(iω) ’ ρFhp +P 鍔(S、)ノD,一
伽 ) ’ P。h。・
Hlr}ε、。
VID ,]…・
……・
・
…………・
………・
・
(23) と な り変 位の像関数 u#(.
Si)の表 現が得ら れ る。
こ こ に,
一
109
一
St
はJn
(Sia
)の 正根 (for
i=
1,
2…
)で あり,
Up (a )お よ びAu
,(α)は未 定 定 数で あ る。 逆変換に よっ て変位 の原 関 数ぱ次式で表さ れ る :u・(・)
一
(
看ゲ
・・(ω+鮎 [
u”s
・)+ α
響
ω 蝋 ・)]
・(・、・)…・
・
…・
・
(24
) こ こ ・,
酪誓
隔 (s
、・)12
…・
一 一 …一 ・
・
…・
(・・) 2.
3 係留索の解 構 造 物は弾性係 留 索に よっ て係 留さ れ る もの と し,
係』
留 索の剛 性を変化さ せ ること に よっ て自 由 浮 体から定 位 値 固 定 浮 体まで を カバー
す る弾 性 係 留 状 態を取り扱うこ との で き る モデルを用い る こと と する。
こ こ で はFig.
3 に示 すと お り構 造 物 下 端 円 周に沿っ て内 外 2系 列の係 留 索が底 板 中心に対し放 射 状に配 置さ れ るもの とする。
ま た,
索は 無 質 量,
線形弾 性ば ねであり,
軸力 だけ を生 じ る もの と し, 流体に よ る作用 は 無視で き る もの とする。 索の上 端 部 (シェ ルお よび 円板 端 部)にお け る x,
z 方 向 変 位 耽,
蹴 お よ び内 外 (q=
1,
2)の係 留 索の円周 方 向 分 布 軸 方 向 力 Nteは,
次 式で表さ れ る :co U、(θ,t)
=
Σ (u(0)+ε1謬 α)cos nθ・
e‘wt n=
o W,(θ,t〕=
εln(γ十XeV )COS θ・
eiωtm
NL9(θ
,
t)=
Σ n、,’
cos nθ・
e ‘b’t n=
DtS−一
一
一
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(26
) こ こ に,
n、P
・
=
λ、・
EA 副[u(0}+ε、。ea
]cos α。一
(−
oq
ε]n〔γ十Xo ψ)sin ae}/Sg・
・
・
・
・
・
・
・
・
…
一・
…
r
(27> で あ り,
α。 は索の 鉛 直 軸との な す角, Se・
=
L
、ICOS
ae は索長で あ る。
λL を係留索の剛 性パ ラメー
タとす るも の と して, λ,・
EA
,。 で円 周方 向 単 位 長さ当た りの分布ば ね剛性を 表す もの と する。
構 造 物の受ける鉛 直・
半 径 方 向 分 布 応 力 (係 留 索の分 布 反 力 )は次 式と な る :N广 乢 c°s α 一 乢
『
°s α・.
.
……….
(28
) ハrrL=
=−
NL,
sin α1十IVL, Sln α z 構 造 物 全 体 系に対 する X,Z
方 向 合 力F
。L,
F
。L お よ びY
軸回 りの モー
メ ン トTr
。 は,
次の と おり係 留 索の分 布 反 力の積分で表され る : FXL− 一
ズ
π膕 ・,
&・イ
・αC・S ・
・
d
・T
・L− 一
ズ
π
鳳 ・+ X・N・・.
)・C・S ・・
d
・・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(29
> 積分の 結果 は集 約 ばね 剛性を 用い た次 式の形で表され る。1
;
一一
[
驢
1
:
]
一
llO
一
εonδ εm γ’
COS θ・
ettOt・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(30) εln、
V・
COS θ こ こに δは鉛 直方 向 剛 体 変 位 (u{O
)atn=0
)を表す。 本 論 数値解 析で は内外索の 剛性お よび 傾 斜 角は同一
(λL、=
λL,=
AL,
EA
,、
=EAIa
=EA
,,
α1竺
α2= a}で ある もの と し,
水面 レ ベ ルの静 的変動 η。 によっ て, 鉛 直 方 向力Fs
= rra2 η 。・
Psg を生じ る剛 性 を 索 基 準 剛 性Ks
に設 定 す る :K
、= πa2ρFg…………・
・
……・
…・
・
…・
…・
・
……
(31 ) こ こ に ρF は 流体密度,
g は重 力の加 速 度で あ る。
こ の 場 合,
基準分布ば ね 剛性EA 。
は次式で表さ れる。
EA
,=
ρFgLt α/(4 cos3 α)・
…
『
・
・
・
・
…
「
・
・
・
・
・
・
・
・
…
(
32
) 索 基 準 剛性を用い て集約ばね剛性を表せば次式と な る。
Kxx
; λ 、Ks,
K
。 、=ALKs・
tan2
a/2Kze=
λ,K
.・
Xotan2
α/2…・
……
(33 )Ke
θ=
λLK3・
(τ乞tanZ
a十a2}/2 2.
4 流 体の基 礎 式および解 構 造 物の径 は流体 粒 子運 動の径に比べ て十分大き く, 波高が 比較 的小さい もの と す れ ば,
非 線 形 流 体 力は無 視 で き,
微 小振幅波理論が適 用で き る15) 。 また,
流 体は非 粘性,
無 渦 流,
非圧縮 性と仮定 す る と,
流 体 速 度ポテ ン シャ ル理論が適用でき る。 i )流体の基 礎 式 :・・ φ
一
(撃
・磐
・・、,・一
・一 ・
一 ……
・・4 ) 速 度 ポテンシャル φ を,ゆ
φ(r,
θ,x, t)=
Σ ψ(r,x)’
cos n θ・
eiω
t・
・
…
(35
) n=
O で表す ものと す れば,
流体基礎式は次式 と な る。
耐
払
一
籌
・ + … x−
…一 ・
……・
一
(36
)ii
)流体の速度ポテンシャル 流 体の解 析は, 領 域 分 割 法 (構 造 物 下部・
外周 領 域に 分 割 )によるもの とし,
各 領 域で の解 を 次に示す。
な お,
流体解 析に お け る鉛直座標π の原点は海底 面とする。ii
)−
1 構 造 物 下 部 領 域の速 度 ポテ ンシャル q、 構 造 物 下 部領域 (r≦α,0
≦x ≦L
,に お け る流体ポテ ン シャ ル g、につ いては,
流 体 基 礎 式を半 径 方 向に有 限Hankel
変換 (inO
≦r≦α),
鉛直方 向に有 限フー
リェCosine
変換 (in O≦コc≦LI )を行う。
関 数 F (r,
x)に 対す る上 記の変換は次 式で定 義さ れ る。f
, L’
xaF
(r,・)・
Jn
(S
・r)… S ・・X・
rd … x ≡FHC
(Si,
ait)・
…一 …………一 …・
・
・
・
……
(37) こ こ Fこ,
ait=h
π/」L
置・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
一・
一・
・
・
・
・
・
・
…
(38) 像関 数ゲC(S
‘,
ak )は下 式と な る。e・HC(・i
・
・aD − 、纛
・1・Jn−
i(・i・)・・(a,
・Cb
十 ψ佐(
SE,
0
)一
(−
1)露 嬬 (S
¢,
L
,)}・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(39> こ こ で,
gc(α,
ak ),
ψ鋒(∫‘,
O),
ψ態(S
‘,
Ll)は, 境 界 値;ψ(a,x),
ψ、
x{r,
0),
g,
.(r,
L,)の像 関 数であっ て 本 解の未 定 定 数で あ る。
これ ら構 造 物 下 部 領 域 速 度 ポテ ンシ ャ ル に関す る像 関 数の逆 変 換に よっ て, 次 形 式で原 関 数 が 表 される。
・・i(r
,
X)−
i
・
藷
・げ(a… )’
・ ・S ・ ・X+
盞
・
癒
み(銅撫 [
・・C(s
… a・)+
響
‘ω・
げ・・,・・]
… S・
a・・X・
1山
…(… 。)=tio
・・gc
(・… )’
・ ・S・… 9・pe(r,
o)一毒
i
赱
qX(s・・o)・
」 (s・r} ・・。(r.
・Li
)〒敲
郷
・
L1
)・
J
・(・s
・ r)・
・
「
…
(40)ii
}−
2、
構造物外周域 流 体の速 度ポ テンシャル ePtl構造物外 周 領 域 (
for
α≦r,
0≦x≦L
)の流体ポテン シャル ep、tと し て は,
自 由 表 面・
海 底 面・
無 限遠 放 射 条 件 を満 足する解Z2)としt
次 式を用い る:姦
(醐 」{
en(一
・略
η恥
・)・ 醐副
’
cosh (ko
〕m
cos (h
丿x) 。。sh (κ。
L
) +i
,A °Kn
(k
’r) 。 。s(κ、L)・
……・
…………・
……・
…
(41) こ こ に,
上式 第 1項は 入射 波 (振 幅 η, 波 数 h。)の ポ テンシャ ル であ る。A
。,
ん は構造 物と流 体の接 触 面 に お ける境 界 条件に よ汐決 定さ れ る未 定 定 数で,
Jn,
Hn, Kn は n 次の第1種ベ ッセル 関数, 第 2種ハ ンケル 関 数,
第2種 変 形ベ ッ セ ル関 数であ る。h
。,
島は次 式 を 満 足 する波 数で ある : ω t=gh
。
・
tanh (κ。ム),
ω』−
gic,・
tan(h
’L
)
一・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
一
一・
・
・
…
(42 ) iii)・
波圧 と波力 シェ ル側 面波 圧 Ps,
円板 底 面 波圧p尸は,
線 形 化さ れ たベル ヌー
イ式 によっ て表さ れ る :P。(x)
=一
ρ市ω・
ψ“(αゆ )+9
[u(Ol+・1。v
α]}pKr )
=一
ρ。lia
・・
q,(r,
Li)+9 [駕P(r)十 EinVr ]}・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
9…
9・
一・
・
・
・
・
・
・
…
(43} 上 式 第2項は静 水 圧 的 波 圧 成 分で あ る。 構造物全 体 系に 対す る鉛 直,
水 平,
回 転 波力Px
。,
.
PZF,
TrF
は,
流 体に 接 触する シェ ル側 面As (0≦ θ≦2π,
L
,≦x ≦L
), 円板 底 面 ん (0≦θ≦2π,
0≦r ≦α〉で の波圧積 分で表さ れ る。
Fig
.
3 Stresses ofthemoQring springs,
pLate and shell actingat the edge of the structureP
・F一
蕉
恥 翩P
・F−−
fX
.
Rs
’
d
・s θ・
αd
・d
・T
・・一
{
ff
. P・’
・(x−
x・1
・
dx
+疏
耳…司
・
COS θ・
d
θ・
・
…
一
一
tt−t・
−t…
tt…
∵・
−s・
t・
tt・
tt・
(44) 積 分の結 果,
Pxs は n=
0め 時に,
またPZF
お よ びTVF
は n=
1の時に の み非 零であり,
次 式で表さ れ る。
P
・F−−
1
ル
・P… ,(・,Li
)…d
・+・吶 δ}
・
・…P
・F−
(
fX
,.i
・PF・
・”(喞 ・ ・s・・伽 ・ト
・… T・F− 一
{
∬
励 … u(・,X>・
・C・S ・ ・(X−
X・)・・C・・・
ft
”‘・PF・
e・ ・(丗 ・ ・S ・θ・
・’・・de
・・a・ PF・[・・/・+L・〈・・
−
x。
)叫
… 1・
一…
一一…
一一・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(45) こ こ に, XB は底 板か ら浮心 ま での距 離 (XE=Le
/2
)で あ る。
上 式の とお り静 水 圧 的 復元力の項は線形 理論に よっ て い る ため, 水平波 力へ の寄 与は な いが,鉛直波力,
回転モー
メ ン トに対する寄 与がある15 )。
2.
5 境界条 件 本解 析 上の境界条 件は以下に示す とお りであ る。
i
)シェ ル上下端に お け る応 力・
変位に関す る条件 : か レ 如…
==
駕 0 囲 咲 呼 エ エ ー ー η η 駕 ω v、(H
}=0…・
………・
・
…
,
M 。(H
)=
0− ・
・
…・
……
, v(o)=
o・
・
・
・
・
・
・
・
・
・
・
・
・
…
∵・
,
w,
。
(O)=−
u,.
r(α)・
……
,
ii)Fig.
3に示す構造物 隅 角 部 分の鉛 直 方 向の釣 合 条 件 (Σ 凡≡
0),
お よ び同 部 分の y軸ま わ り回 転モー
メ ン トの釣 合条 件 (Σ臨;
0);n
=
(・)−
v・
・(・)・耐膿
’ 醐d
・.
一
伽 )’ ρ,h
、hp
[u,(a)+ε1。Va
]=0…・
・
…・
・
・
・
…
Mx (O)+η 。。(α)一
伽 )悔ん。ん、(h
}+鵬 /12・
(u。.
。(α)+ε 1。の=
0…・
・
…………・
・
…・
……・
・
こ こ に,
v。
p(α)は底 板 端 部 置 換せ ん断 力である。
iii)底 板 全 体の水 平 力の釣 合 条 件 (Σ Fz=
0, at n=
1>;ズ秘
・ ・S θ一Ve
・… ]+・FZL一
蕉
一
IL1
一
Ps
・ ・s θd
・・}
・d
・一
{・・) ・ p・… α ・ (・・… )・
… t−
・・
………・
・
……・
……・
…一
iv
)シェ ルー
底 板 構 造 物 全 体 系のy
軸 回り回 転モー
メ ン トの釣 合 条 件 (ΣTv=
・
O,
at n=1
); Tvs+ Tv。+ T,L+ T,F=
0…・
…………・
…一 ・
…
ここに,Trs,
Tvp
は シェ ル と底 板の回 転モー
メ ン トで あっ て次 式で表さ れ る :乃广 ・ ・
イτ
[σ・・
α C・S 卿 ・・
・…一
四、ビcos θ)(x−
x、)]αd
θdx − ・
一
(a }T
・p−−
P・h・ズ∬
陬 ・ C・S θ+(−
V・.
tt・
sin θ十Wp,
tE・
cos θ}Xo]rdrd θ……・
(b
) v )構 造物 外周流 体に対する自 由表面の条件,
お よ び海 底 面に お け る鉛直方 向流体 粒子 速度零の条 件であり,
こ れ ら は既に (41
)式の解で満足 してい る :il
、i,
tt+gdirt
,
x=
O・
・
…・
,ip
、1,
x=
O……
vi )構 造 物 下 部 領域流 体に対す る底板面に お け る流 体 粒 子と底 板の鉛直方 向運 動 速 度一
致の 条 件, お よ び海 底 面 にお け る流 体 粒 子鉛直方向速度が零の条 件 :il
、,
x=Unt……
, φ,,
x=
O……
vii)r= a 境 界 面 上で o≦x ≦L の領 域におい て,
シェ ル壁と構 造 物 下 部 流 体の半 侵 方 向 運 動速度が構造物 外周 流 体の半 径 方 向 運 動 速 度に一
致の条件;徽
:
1
:
欝
訓
…・
……・
…・
…・
・
Viii)構 造 物 下 部 流 体 圧 力と外 周 流 体 圧 カー
致の条 件;dii
,
,=ilti
,
t (inO
≦x ≦L
、)・
…・
……一 ・
…・
・
…
以上,〜
の境界条件と ともに, それ に対応する く未 定 定 数〉 をFig.
4に示す。
2.6
弾 性円筒 構造物の解 以 上の円 筒シェ ル,
円板,
係 留索お よ び流体の解を 用 い,
境 界 条 件 を表せ ば次 形 式の 複 素 連 立 方 程 式を得る。 [ノ4
]1Xl
=
IBI
・
・
・
・
・
・
・
・
・
…
噛
曁
…
9・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
《46) こ こに,
EA
],
IB
}は,
そ れ ぞれ係 数マ ト リッ ク ス,
荷重 列ベ ク トル で あ り, {X }は未 知ベ ク トル で あ る。 本 解 未 定 定 数に は,
以 下の ものが含ま れ る(Fig.
4参 照)。 i)シェ ル両 端 弾 性変位と応 力 関 係の 8元 ;nx (O ),
v(0
), ω(0
), m =(0
), nエ(H
), v(H
), w(H
),
Mr (H
),
ii
) 円 板 関 係の 2元;Up (α),
△伽 (α),
iii) 剛 体 変 位 関係の Z元 ;γ,
ψ,
iv)構 造 物 下 部 流 体 ePI関 係の 1+K
+ 21’
元 ;ψc (α, α。
)〜
ψ c( α,’
α κ), 鵐(S
、,0
)〜
撮 (Sn,0
),
轍(Si,
L
,)〜
epM(Sl
・
,
L
,)。
こ こ に,
K ,1
’
は 鉛直方向フー
リエ 級 数,
半 径 方 向ベ ッ セ ル級 数型の解に 対す る数 値 解 析上の算入項 数である。
v ) 外 周 流 体 ψ 皿 関係の1
+J
元 ;A
。〜A
,であ り,J
は同上 項 数である。 以 上,
解析上算入 さ れ る全 未 定 定 数の数は 14+K
+2 ∬’
+J
元で ある が, こ の内, 境界条件 , , , , よ り直 接 単 独に零と決ま る もの お よび境 界 条 件 , より,
式 上 消 去で きるものは省 くことができ,
最 終 的に は境界条件 , , , , , , , , に対 応 す る ∬’ +9元の {Xl
(n=
1の時 ), ある い は境 界 条 件,
が 無 用な n≠1の時に は1
’
+7元の }淵 を連 立 方 程 式と し て解くことに よっ て未 定定数 を決定す6“
:iXl
; [A
]一
匸ヨBI ・
・
曾
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(47 ) }Xl が求ま る と,
引き続き構 造 物と流 体 応 答 値は決 定 で き る。
3.
数 値 解 析 結 果 と考 察 無 限 級 数 型の本 解を用い た数 値 解 析上,
精 度の よい結 果 を 得るに は計 算 算入項 数はで き る限り多い方が よい が,
こ こ で は計算項 数を多く採っ て得ら れ る最終 的な値 に対 して 3 % 程 度の誤 差 に収 まる有限 項 数 を採る こ と と し,
周方向波 数は n=
O〜
5 (剛体解 析で剛体の応 答の みを求め る場 合はn=0,
1),
シェ ル の鉛 直方 向フー
リェ 級 数の項 数はm=
O−
30, また, 構 造 物 下部流体に関す る項 数はK =
30,
1’
=
15,
外 周 流 体に関す る項 数はJ
= 5と した。
3,
1 剛円柱 構造物の応答i
)自 由浮体の応 答 円 柱 浮 体モ デル 〔半 径 a=
O.
lm,
水 深L
=O.
4m,
底 板一
重心 問 距 離 x。=qL
/2,
流 体 密 度 ρF=
ユ000
kgf
/m3 /g >に設 定し た場 合に対 する調 和 波 応 答 解 析 並 びに実 験が合田 ほ かに よっ て扱わ れ た (文 献15)。 そ の φ1い 驚 +9φII.
x=
0rL
llilll
:
lii
}
1
ヤ
ー
≒一
ΣTゾ 。 (Ψ> L Xo 。 z ΣFz・
0 〈7) 「三 u(0)・
Up(ω 〈nk (0)〉 vω }・0 〈v (0)〉 蝋0》・
O 〈w ω )〉 叱罵ω )=
−
Up,
.(q) 〈mx( 〉 Σ Fx=
0 〈Up (ω 〉 Σ My・
0F
〈△Up (α)〉 、、1
φ・
・
ノ U・・
t (9,
xH (S二,
L1)〉’
x1
φL .・
φ11.
. , w.
t・
φILr (A亅〉I
or/ φ1
,
.・
0(?
.
罵 H (Si,
D)〉 φい=
φ一 (乎じ (α・
哩
/ φH.
rOト
ー 一一一一一
〇一一一 一一一
爿
Fig
.
4 Boundary conditions and corresponding unknown constant $ involved in the analysisうち吃水比 (q
=
L,/L >を変化させ た もの をModel l〜
Model3
と して考 察す る。 回 転 慣 性モー
メ ン トとし ては,
構 造物 質量 が浮体 高さ中 央 部に集 中してディ ス ク 状に分布す る と し た 場合 (1三 勾,
中実シ リンダー
とし て全 体 的に分 布す る と し た場 合 (∬=
lc),
側 壁と底 板は 同一
面 密 度 (Pshs=
p,h
,}を持つ もの と して構 造 物 周 辺 に質 量が面 的に分 布する とし た場 合 (1=
lsρ)の 3ケー
ス を想 定し,
Icを基 準 値とする形で Table 1に示し て、
い る。 な お,
本 解 析では シェ ルと底 板は共 通に剛性Eh
を 大きな値と な る よ うに設 定 (h =
0.
002m,
E=
2.
1× ユOPkgf/me ) する ほ か,
u= O.
3,
ds
・tdp=
0とし,
また , 浮体高さ は同 上文献に揃 えてH =
qL に設 定 し た。Fig.
5に解析結果と して重 心 点 水平変位周波数応 答曲線む
を示す と と もに, 同上文献か らの引用 実 験 結 果 (プロ ッ ト白 丸, 黒丸は入射 波 高がほ ぼ2cm , 1cm の時の結果〉 を 示す。
』
な お,
振 動 数に は無 次 元 振 動 数の 自乗 値 ξ= ω ’L /g
を使 用する。
本 解 析 値は実 験 値に大 略の一
致を 示す。 解 析 結 果に現れている共 振ピー
クは実 験 結 果に は 見ら れない が,
実 際 上は流 体の粘 性め影 響があるた め と 考え られ る。Table
2に は本 解固有 振 動 数 (共 振 応答曲 線か ら算定)と と もに, 同 上 文 献の 自由 振 動実験
値 (自 由 減 衰 記 録か らの 読 取 値で,Fig.
5にAExp.
の マー
ク で印す),・
同 じく剛円柱と して の解 析 結果 を示す 。 文献 15)に よ る実 験で は回 転 慣 性モー
メ ン トが不 明で あ る が, 本 報の パ ラ メ トリック解析か ら,
回 転 慣 性モー
メ ン トの 大き い結 果が実 験 結 果に よ く一
致 し.
て お り,
実 験 時の 回 転 慣 性モー
メ ン トの 値は ∬=
tc(
Model 1,
2の時 )か ら1=lsp
(Model
3
の時 )付近 に あっ たこと と推 定さ れ る。
ii
)係 留 浮 体の応 答Fig.
6に は索剛性1
パ ラメー
タλ,を変化 (λL=O,
1,
5,
Tablel Parameters used in this analysis
Model
1
2
3
q
3
!4
1
!2
14
IC
(xlO『
2 ¶s2} 〔is
+IP
}ノIC
Id
/Ic
0
.
9621
.
380
.
25
0
.
3741
.
50
’
0.
43
0.
1071
.
580
.
75
1
・・Ma2
/4
,1
・・
H [
02+H213
]ノ4,
ISP
・1
、+1
,{a) こ こに
、
Ip
=ρpπhp
〔12〔a?+hp213
)14
十Mpx
。2・
Is
=ρ sπH {
a22 (a22十H213
}−
a12 (ai2十H213
}}ノ4
+
M
・{II
/2−
x・
)2,
・i・ q−h
・12
, a 2・q+h
・12
,
x。・
H
Ms
/2H
,H
=πa2L2 ρ F,Ms
=2
π aHhsPs ,Mp
=π 〔】 2hp ρP.
{
b
} 1000)さ せ て得 ら れ る鉛 直お よび水 平 方 向の変 位 δ=
δ/η,
γ=
γ/ηと波 力PXF=PXF
/πα2ηρFg,
PZF=PXF
/2 aL 、ηp,・
g の周波 数応答 解 析 結 果を示す。 自 由 浮 体 (λL=
0 )の状態か らi 係 留索剛性 を増 大さ せて い く と共 振 振 動 数が増大し変
位 振幅は低 下してい く。 鉛 直・
水 平 波 力 応答 曲線に おいて は, 係留索剛 性の増 大に連れて波 力 共 振 ピー
ク は鈍 化 し,
ほ ぼ定位置固定 (λ,=
1000)の場 合 には,
ピー
クのない滑ら か な応 答 曲 線になる。
実 験 結 果 と しての自由浮 体の変位応 答 結果
および定 位 置 固 定 時の 波力 応 答 結果は解 析 結果に「
よ く一
致 して お り,
本解 析は 妥 当な もの と判 断さ れ る。
係 留索剛性が変化す れ ば構 造 物の変 位お よ び波 力は顕 著な変 化 を示 し, 係畜
索剛性の 応 答へ の影 響は大き い。
・
3.2
弾 性 構 造 物の応 答水 深 半 径比
L
/a=
4の海 域に シェ ル高’
さ半 径 比 H/a =3,
吃 水 比 qul /2 と じて,
シェ ル と 底板 材料は同一
(鋼 )で あ り,・
ヤン グ率E ’
s=E
,= 2.
1
×loiOkgf/m2,
’
ボ アッ ソ ン比 Vs・=
・
Vp=0.3
, 減 衰比d5=
俳=0
であ る場 合Table 2 Natural freqtiencies ξ
,
,
ξr in vertical and rocking vibration for model 1〜
3・
ξ×in
vertical vib.
ξYin rocki vib.
Hode11Mode12 凹ode13Hodel
1Hodel
2Hode13
Thisel
母stic analysis1=IdI
=
ICI =ISP
1.
101
.
101.
10
1.591
.
591.
59
2.
772
.
772.
770
.
530
.
250.
231.
250
.
680
.
483.
0
α2.
401
.
65Exp 。
result、
(Ref 。
15)Rigid
analysis (Ref.15
)1
.
06LO61.
551.
462.
792.
65
0.
300.
310
.
660.
50 1.
711.
55
3
2ttS
、
1
0
ISP
Ic
Idl9
.
r,。e4 %Se。 。 %
Hodel
1
(q・314
)tE
:P.
2
3
21tS1
e ξ 0 34
ISP
!
,
1
・1
111レ
1
・ ヴ 。。。’
G.
。
5.卩
■
Hode12
〔
q
・1
/2
〕・
,
△.
tEXP
2
33
21
トt
ξD
4
Hodel
3
ISP
(q・1
/4
} ムー1
‘ ,ノ
L
ノ1
・ oo■
,
o騨
9r.
り9
ユ EK ,.
23
4
ξ’
4Fig
.
5
Variation of horizontaldisplacement
response γwith rotational iner【ia momen し1
(’
AExp.
;Experimenta且resu [t of natura 【frequency
,
Ref.
15}43210
一
δ嚇
ん..
_
2、
臥、
’ ノ改
432一
δ 1 2 3 ‘ξ 0ー
・ ・訊
畢
1
唱
ー
\
卩 % 、0
7
−−ー
ハー
OO1 4 3 21 2.
o1.
51、
0 : s一
幽
}一一
一一、
_ 、・
o o.
5 OO δ4
」 4321o一
δ 23
4ξ
0
1 a) }畳odel 1 (q=
3
/4
〕 2.
0 島ll
昌 昌ハ
/1
\ 1 2 3 4 4■
321 ξ o7rIi1111
1
tl
.Yl
.・
\
こ
轟
一
1 2 3 4 b》阨del
1,
51 』 oo 0.
5\
.
°
% 6メ
一
一
一
冒
一
一
〇
一
一
〇
■
畳
bo
ノ「
一
『
緬
・
陣
0 〔 05221
ξ 2 3 4 43一
γ、
・
、
(
×・
丶 ゆ 丶 丶 、〆 。
伽
ノ
, ’ー
2ユ 1
.
0G.
5 2 3 4\
畿
聾
芟
、
.
、
、
2.
o1.
51.
oo.
51
レ
v/ e 1 2.
0 ユ.
51.
o0.
50竺
ー 巳魂
ノ,
0 7。
σ Z II ご一
ヒ
ロ
ロ
P ー,
−・
ー・
ノ「
PZF 昌 !lll
』
拶一 1} ’
r ξ 2 3 4 2.
o1.
51.
oo.
5!
l
導
睾
ξ 1 2 3 4ξ ・ 1234 ξ ゜ 123 、 ξ ゜
123 ・ c)岡Ddel 3 (q
・
ユ!4
)Experimental r巳sults (Ref
.
15
〕 O eR巳sll1しs by this analysis ;
一
λL=O,一一一一
λL:1,
一・
一
λし=
5,一
一
一
一
一
λしニ1000
Variation of displacement responses δ
,
γ,
and wave forces Pxr,
PzF with 【igidity ofthe mooring spring λLfordifferent
draft ratio qξ 丶