[设计与研究]悬链线方程的求解及其应用胡灵斌(上海三航工程船舶公司上海200032唐军(上海船舶研究设计院上海200032[关键词]悬链线;方程[摘要]本文结合实例介绍了悬链线方程的求解方法及其在船舶设计中的应用。[]U664.4+1[文献标识码]A[]1001-9855(200401-0017-04Solemnestandapplicationofcatenary'scurveequationHuLingbinTangJunKeywords:catenary'scurve;equationAbstract:Solutionandapplicationofcatenary'scurveequationinshipdesignispresentedinthepaperwithexamples.1前言船舶设计者都知道按舾装数法来确定锚、系泊设备,而对于用抓持力法来确定锚、系泊设备则可能由于其计算的局限性而令广大船舶设计人员感到陌生。然而,如灯船、电站船、起重船等工程船的锚泊设备用抓持力法来确定为好,则悬链线方程的应用是必不可少的;而对于较为复杂的多根锚链泊碇的浮船坞,其外力通过专用程序计算得到后确定其链力也要用到悬链线方程。2悬链线方程的求解众所周知,对应图1有以下著名的悬链线方程。y=asecTchxa-1+tgTshxa=Rw1cosTexwR+e-xwR2-1+tgTexwR-e-xwR2(1图1锚泊悬链线式中y——水深,m;R——作用在船上的外力,t;x——锚泊点至挂链点的水平距离,m;T——起锚角,°;w——每米链重,t/m。从以上公式中可以看出y是R,x,T,w的函数,即y=f(R,x,T,w。已知R,x,T,w就可以求得水深y。而实际上船舶抛锚时,通常设置有足够长度的锚链,以便有一段躺底锚链,使锚不产生上拔力,这样,T=0;而如果又选定了锚链,则每米链重w也为17[收稿日期]2003-4-29[作者简介]胡灵斌(1962-,男,工程师,上海交通大学在读工业工程研究生,从事船舶研究设计工作。唐军(1946.8-,男,汉族,江苏无锡人,高级工程师,从事船舶总体设计与研究工作。常量;又如果作用在船上的外力R也已求得,这样,公式y=f(R,x,T,w就可写成y=f1(x。但是船在某水域抛锚,则水深y反而是已知的。所以往往要从y来求x,显然,从公式y=f1(x要推导出它的反函数x=f2(y来求x是很困难的。只能用不同的x值代入公式(1去凑y,当计算得到的y恰好等于水深h时,x就是我们要求的值。这的确很繁琐,好在我们现在已处在计算机时代,我们可以利用Ex-cel软件的强大功能或编制某计算程序,让计算机来完成。求得x后,接下来再接以下公式求悬链长。l=atgTchxa-1+secTshxa=Rw1tgTexwR+e-xwR2-1+1cosTexwR-e-xwR2(2对于用链索泊碇的浮船坞,由于受水域限制,锚不能抛得很远;而采用沉石,这就允许有适量的上拔力,但不能太大,要求出该上拔力的大小,这也可以编制另一个程序求解悬链线方程得到。在浮船坞的链索泊碇工程中碰到的另一个问题是:有时要求出在测定的挂链力T下,且已知水深y,想要求出使锚不产生上拔力,锚至少应该抛多远(x=?,即要求得到x=f3(y,R这一关系式。而R又是T1的水平分力,显然,从公式(1要推导出x=f4(y,T1也很困难。我们也可以编制某计算程序,程序框图见图2(a,利用计算机每秒钟几千万次的计算速度,让其按一定规律从较大的水平阻力悬链线方程曲线开始进行迭代计算搜索,搜索过程见图2(b,当搜索到某一点,这一点的y恰好等于设定水深h,T1等于测定的挂链力T。这时,x就是我们要求的最小抛锚距离。同样,接下来就可以按公式(2求得悬链长。3悬链线方程的应用实例一用抓持方法,即根据锚及锚索在水底抓住泥土时产生的抓持力与作用在船舶上的外力平衡条件来确定锚及锚素的尺寸。以某船为例:船长97m,船宽15.8m,型深7.7m,吃水3.8m。风阻力:R1=0.006dV2(F′+0.31FkN=0.006×0.124×30(107+0.31×61.6=83.93kN图2搜索计算程序框图式中d——空气密度,为0.124kg.s2/m4V——风速,海船V=30m/sF′——上层建筑在舯剖面上的投影面积,F′=107m2F——船体干舷部分在舯剖面上的投影面积,F=61.6m2水流阻力:R2=0.01fKv1.83kN=0.01×0.18×1500×1.51.83=5.70kN式中f——阻力系数,取0.18K——船体浸水面积,1500m2v——水流速度,取v=1.5m/s波浪漂流阻力:R3=34kdgB(14H1/32cos0°×10-2kN=0.75×0.7×104.61×9.8×15.8×(0.25×1.52×1.0×10-2=11.96kN式中K——系数,0.7d——水的密度,104.61kg.s2/m4g——重力加速度9.8m/s2B——船宽15.8mH1/3——有义波高1.5m18作用在船上的外力:R=R1+R2+R3=83.93+5.70+11.96=101.59kN...