摇板式造波机所造二维规则波的时域线性解.docVIP

摇板式造波机所造二维规则波的时域线性解 PLa 【.一印 第2期(总第121崩) 1993年4月 造眨机,镏抱喇,成似带 币国造船 SHIPBUILDINGOFCHINA No.0(Ser^JNo.121) Apt1993 摇板式造波机所造二维规则波 的时域线性解 (哈尔滨船舶工程学院) 提要 本文给出摇扳武造渡机所造二维规则玻时域线性解的样析表达武,它描述 一 个带有玻前的二维线性玻的传播过程.当摇板作简谐摇荡的时问趋于无穷大 时,该解趋干线性稳态解.时域线性解可为求解时域非线性解提供自由表面速 度势和波形的韧始值,目此具有实际应用价值. (一)引言 摇板式造渡机造=维规则渡是一个典型的二维摇荡辐射问题.在频率域中求稳态线 性解及二阶解一般可采用数值方法,但由于摇板式造波机几何上的特殊性,即摇板的平 均位置是垂直的,二维规则波的频域解可以解析表达,论文[I3做了这一工作 在少数情况下,可以写出线性时域解的解析表达式,例如wM.Lin等I984年的 论文.中给出了一个范侧.本文讨论的摇板式造波机所造二维规则波的时域线性解也有 解析表达式.它可以显示二维规则波线性解的瞬态形式考察其趋于稳态解的过程,可 为在时域中求解线性绕射和浮体运动响应提供一个二维入射波模型,也可为步进法求解 时域非线性问题提供台理的初值. (=)时域线性解的表达式. 1.线性氇态解 围1是水池沿长度方向的剖面,池深H,摇板在静水以下高度^,坐标系取法为 轴与静水面重台,y轴与摇板的平均位置重台,垂直向上 本文于l9年S月3日救到 图1水池剖面与坐标系 e中国迫船1903年 设摇板作筒谐运动.角速度=sincot?则线性稳态解为 一一4一 CkohshkoH :(T 2 cb H ko + H+ sh2 chkHo( ) H--h).3一 ch(+H)eos(k一州) +4 荤(k,,fisink.H+吨eoHsk恤.H--2e㈣osk.(tt-h)3.碱(j.+ (1) 其中.kothk~H=co/9.k.tgk.H=一∞/9. (1)式表示线性稳态解由一项沿正方向传播的行进波和无穷多项衰减波组成,在 远离造波机处,仅剩下正向行进波. 2.时域线性解的推导 如果摇板从时刻f一0起作简谐运动.角速度—s.f.则时域线性解可分解为瞬 时部分和记忆部分之和 (x,Y,f)一T(x,Y)sincot+Y(,y,f)(2) . (x,j,,lgt;-S:(x,j,,f—t)sin.dt(3) 瞬时部分的定解条件是 IV一一o(在域内) It—o(Y=0) … 一 ¨(=㈩ f三oy一0(j,=一H) 瞬为 ∑6[口-(j,+H),口_一(5)m●… 代人物面条件决定系数k,身知 6一击J:^(e+^)cos【(+H)]d ;(旦(6) Hd:L 从而 堕Oy一一 音c川(c一cⅧ 利用_F面的式子来改变(7)的形式 (7) (+肋一-告(8)l…～∞… 其中积分路径(1-是图2所示复平面上围道ct当月一..时的极IV.情况.易知无穷大 2鞋(总第121期J埂五溯等.樯板式遗渡机断墙二堆用研浊的时域线性奄阜 O 田2星平阿圈 半圆周上积分为零,故(8)式右端等于 睾: 于是 兽l…一一÷:(+Ⅷi 园为 :^(+^)ch矗(+日)d毒—争(与^shhH--chkH+chk(H--)]垒Ff膏) 所以 堕 ayl.o一一÷Jr0chk㈩ly?…… 顺便指出,毒+时,F(k)以—争为极限. 再求(,y,1).其定解条件为 +罟一.(一o) 号一0(一日) 告一0(:0) lr.6;0 詈l…一等c=. 易知 一 2g J.(+H)半(11) 满足(1o)中所有条件,其中=/磊r.这样我们就得到记忆部分 中国造船t993年 等厂:8旦chk(y+H)dkfotdr 进一步可将乱写成 一 警旦盎嚣chh(Y+H)d +警/chk(十H)等?号曲) 记 G--dF㈩?丽1丽(13) 则(12)式右端前一项积分为 —-咖ReI/Tc.~1(14) 将这个主值积分转换成复平面上正实轴和正虚轴上G的诸奇点留数来计算.实轴上奇 点为.虚轴上奇点为池和f口～它们均是一阶极点如图3所示 k O 剧a夏平皿囝22朋性点 (1)‰的贡献 Gdk--号鹄 对()式的贡献为 —一 4~kohshko … H--chko + H+ehko( , H-h)]ehk 口()r+H)siHsiD耐(15)HaL十sh2k0H,……～0 (2)的贡献 为计算厂.上的积分.把G的表达式(13)改写一下 G一()掰 fr— Gdk=~e~*[(-1)a=h+eosa..( 一 H--h, )]cosa.(Y+H)一 诸f对(14)式的贡献为 三 1 奇()抖∞一.毗(+H)sin耐(16) 由(5),(6)式可