海上横向补给高架索系统[1]的设计受制于绞车负荷能力、甲板高度、舷侧布置和补给品载荷等因素,在输运过程中应确保补给品不与海面和船体相碰。在实际作业过程中,船舶运动与高架索张力间存在显著的动力学耦合作用[2],波浪补偿装置控制系统存在一定的滞后效应[3],这使得高架索系统的精确动态模拟非常复杂。在工程上,常采用静力学方法对高架索的线形进行分析,以保证海上补给作业的安全可靠。张良欣[4, 5]等采用分段受力分析的方法得到了高架索挠度表达式,分析了补给品的位置和重量、张力、补给距离对挠度的影响。卢永锦[6]利用弹性力学原理建立了海上横向补给高架索系统的静态分析模型,并通过样机试验验证了其方法的可靠性。余建星等[7]基于上述理论,给出了考虑波峰波谷影响的实际算例,并从高架索张力、钢索挠度和波浪补偿性能等角度提出了适用于小型船舶补给的高架索海上补给装置[8]。
上述文献或假设高架索线形为以悬挂点为拐点的2段直线,或在考虑高架索曲线挠度的同时忽略了高架索自重的变化,因此计算结果并不精确,只适用于集中载荷远大于高架索自重时的估算。本文基于悬链线方程,全面考虑了高架索曲线挠度对高架索线形的影响,得到了更精确的高架索系统静力学分析模型,探讨了高架索线形与高架索自重和张力等参数的关系,并给出了关于高架索最大挠度和最低位置的一般结论。
1 力学分析模型 1.1 基本假设1)忽略船体在波浪中的运动;
2)假设高架索为质量均匀分布的柔性构件,仅受轴向拉伸作用,忽略其轴向弹性变形;
3)假设高架索水平张力恒定;
4)假设补给品在高架索上的运动缓慢平稳。
基于上述假设,可以将高架索上的补给品输运过程视为准静态过程,开展静力学分析。
1.2 静力学分析在补给品输运过程中,悬架索系统可能存在如图 1和图 2所示的2种线形:线形I和线形II,其中点A为接收点,点B为输送点,点O为补给品悬挂点。
图 3为高架索系统为线形I时的布置情况和受力分析图。高架索单位长度质量为w,接收点A和输送点B的横向间距为L,垂向间距为H。选择点A为系统坐标原点,距原点XG处为货物悬挂点O,补给品重G,OA间的垂向距离为ZG。以点O为界,高架索被分成S1和S2两段,每段均可视为单独的悬链线。假设FA,FAH,TAV,FB,FBH和TBV分别为点A和B处的轴向张力、水平张力和垂直张力,高架索任意位置处横截面上的轴向张力为T,与水平面的夹角为θ。
在S1段内,定义原点为A的局部坐标系。取距点A水平距离为xs的一段高架索做力学分析,其长度为S。
$ \sum{{{F}_{x}}}={{F}_{A}}\cos {{\theta }_{A}}=T\cos \theta \text{,} $ | (1) |
$ \sum{{{F}_{y}}=T\sin \theta -{{F}_{A}}\sin {{\theta }_{A}}-w\cdot S=0}. $ | (2) |
由于货物载荷仅在垂直方向,所以高架索各截面处的轴向张力水平分量相同,如式(1),记作H0。则式(2)可写成:
$ {{H}_{0}}\tan \theta -{{H}_{0}}\tan {{\theta }_{A}}-w\cdot S=0, $ | (3) |
两边取微分,并引入$\text{d}S=\sqrt{1+{{\tan }^{2}}\theta }\text{d}x$,可得:
$ \frac{\text{d}\left(\tan \theta \right)}{\sqrt{1+{{\tan }^{2}}\theta }}=\frac{w}{{{H}_{0}}}\text{d}x, $ | (4) |
对上式两边进行积分,可得式(5)。
$ \left\{ \begin{matrix} \theta =\arctan \left( sh\left( \frac{w}{{{H}_{0}}}{{x}_{s}}+{{C}_{1}} \right) \right), \\ \begin{array}{*{35}{l}} z=\frac{{{H}_{0}}}{w}ch\left( \frac{w}{{{H}_{0}}}{{x}_{s}}+{{C}_{1}} \right)+{{C}_{2}}, \\ S=\frac{{{H}_{0}}}{w}sh\left( \frac{w}{{{H}_{0}}}{{x}_{s}}+{{C}_{1}} \right)+{{C}_{3}}\text{.} \\ \end{array} \\ \end{matrix} \right. $ | (5) |
其中系数C1、 C2和C3取决于悬链线的边界条件。
上述方法同样适用于S2段悬链线方程的推导,区别仅在于局部坐标系原点需取在点O处。
由此可以得到S1段和S2段上点O处的切线水平夹角。
再以点O处的补给为对象,进行力学分析,可得式(6):
$ \begin{align} & {{H}_{0}}sh\left( \frac{w}{{{H}_{0}}}{{X}_{G}}+ash\left( \frac{w{{Z}_{G}}}{2{{H}_{0}}}\cdot \frac{1}{sh\left( \frac{w{{X}_{G}}}{2{{H}_{0}}} \right)} \right)-\frac{w{{X}_{G}}}{2{{H}_{0}}} \right)+G= \\ & {{H}_{0}}sh\left( ash\left( \frac{w\left( H-{{Z}_{G}} \right)}{2{{H}_{0}}}\cdot \frac{1}{sh\left( \frac{w\left( L-{{X}_{G}} \right)}{2{{H}_{0}}} \right)} \right)-\frac{w\left( L-{{X}_{G}} \right)}{2{{H}_{0}}} \right), \\ \end{align} $ | (6) |
当高架索系统为线形II时,同理可得:
$ \begin{align} & {{H}_{0}}sh\left( ash\left( \frac{w\left( H-{{Z}_{G}} \right)}{2{{H}_{0}}}\cdot \frac{1}{sh\left( \frac{w\left( L-{{X}_{G}} \right)}{2{{H}_{0}}} \right)} \right)-\frac{w\left( L-{{X}_{G}} \right)}{2{{H}_{0}}} \right)+ \\ & {{H}_{0}}sh\left( ash\left( \frac{-w{{Z}_{G}}}{2{{H}_{0}}}\cdot \frac{1}{sh\left( \frac{w{{X}_{G}}}{2{{H}_{0}}} \right)} \right)-\frac{w{{X}_{G}}}{2{{H}_{0}}} \right)=G. \\ \end{align} $ | (7) |
根据式(6)和式(7),通过数值方法可以得到ZG的值,进而得到2段悬链线的数学表达式和高架索各点处的受力情况。
1.3 力学模型验证选取文献[5]中的算例进行对比计算,计算模型的主要参数列于表 1。
图 4为整个补给过程中悬挂点处的挠度曲线,两者几乎完全重合。这表明在当前补给条件下,无论是本文还是文献中的方法均能准确给出悬挂点处的位移。
但当高架索张力较小(H0=15 kN)时,2种方法的区别非常明显。
尤其当悬挂点靠近接收点时,高架索从线形I过渡到线形II,其悬挂点挠度迅速增大,更容易发生补给品落水事故。
文献中的方法忽略了高架索变形带来的自重增加,因此低估了悬挂点位置和高架索变形,在工程实际中可能会带来安全问题,有必要采用本文提出的悬链线方法开展静力学分析。
由式(6)和式(7)可知,高架索的线形主要取决于索端横向间距L、索端垂向间距H、补给品重G、高架索单位长度质量W和水平张力H0。其中高架索索端横向间距和垂向间距多根据水动力特性[9]和补给俯角[4]要求选取,水平张力多根据补给品和架索重量选择。本文主要探讨架索自重和水平张力对高架索线形的影响,不作特别说明的情况下,高架索参数仍取表 1中的数值。
2.1 架索单位长度质量对线形的影响图 6为架索单位长度质量变化时不同位置悬挂点处的挠度大小。架索变重相当于补给品重量增加,此时相同悬挂点处挠度增大。当架索自重较小时,挠度较大区域集中在架索中间段,越靠近架索两端,挠度越小;但当架索自重较大时,靠近接收点处的挠度也明显增加。
图 7为相应的不同位置悬挂点处的垂向位置。由于补给端与接收端间存在高度差,所以挠度大的地方并不一定垂向位置低。在文中给出的架索自重变化范围内,垂向位置较低区域始终靠近补给品接收端;随着架索自重增加,相同悬挂点处垂向位置逐渐变低。
图 8和图 9进一步给出了架索单位长度质量不同时架索挠度最大值、垂向位置最低位置及其对应的横向位置。随着架索单位长度质量增加,其最大挠度线性增加,但其所在横向位置基本不变,均为架索中间区域;其最低位置呈抛物线减小趋势,但其所在横向位置均在靠近接受端0~4 m区间内。
在船舶实际补给过程中,水平张力是最主要的控制变量。图 10和图 11分别给出了水平张力变化时不同位置悬挂点的挠度大小。显然,仅当H0/G < 3时,在靠近接收端处架索会发生很大的挠度变形,补给品垂向位置偏低,易发生事故。
图 12为水平拉力不同时架索挠度最大值及其对应的横向位置。当张力H0/G < 3时,架索发生明显变形且挠度随着张力的增加缓慢减小;随着张力的进一步提高,架索挠度迅速减小并逐渐趋于稳定值,且其对应横向位置始终在补给端;当到达H0=20 kN时候,架索挠度已经可以忽略不计。
图 13为水平拉力不同时架索垂向最低位置及其对应的横向位置。当张力H0/G < 4时,架索最低位置低于接收点,补给品可能存在落水的危险。当H0/G > 4时,架索垂向最低位置趋于稳定,悬挂点在补给品输运过程中基本均位于水面以上,因此工程常取用H0/G=4 [1]。在整个水平张力变化范围内,架索最低点对应的横向位置始终位于距接收点0~4 m范围内,且当H0/G > 4.5,接收点即为架索的最低点。
本文从悬链线方程出发,建立了计入架索长度和自重变化的海上横向补给高架索系统静力学分析模型,并通过与文献数据的对比,验证了本文方法的可靠性和普适性。基于上述静力学方程,定量研究了架索单位长度质量和水平张力等参数对高架索线形的影响。
研究结果表明,高架索在输运过程中可能存在2种线形,架索单位长度质量和水平张力的变化会显著影响架索线形,其挠度最大处和垂向位置最低处并不一致。在保证架索张力的前提上应尽量减小单位长度质量,当水平张力与补给品重量比值大于4时,架索线形已经趋于稳定。
[1] | 严梅剑. 海上航行横向补给装置选型设计[J]. 船舶, 2004(1):51-53. |
[2] | 邓凯, 李红涛, 余建星. 高架索航行补给中船舶在波浪中的运动性能研究[J]. 船舶力学, 2009, 13(2):217-225. |
[3] | 赵小兵, 赵伟, 马良荔, 等. Ian航行横向补给高架索张力及张力波动控制参数仿真[J]. 舰船电子工程, 2007, 27(5):115-117. |
[4] | 何学军, 张良欣, 任爱娣. 海上横向补给系统高架索参数影响分析[J]. 海军工程大学学报, 2009, 21(5):52-56, 66. |
[5] | 马仁刚, 张良欣, 何学军. 横向干货补给系统高架索的静力学分析[J]. 舰船科学技术, 2012, 34(2):47-49, 80. |
[6] | 卢永锦. 海上横向补给高架索道系统数学模型研究[J]. 中国造船, 1996(1):19-22. |
[7] | 阎宏生, 余建星, 郭海涛. 海上轻型高架索补给系统的力学分析及安全性研究[J]. 海洋技术, 2006, 25(2):39-43. |
[8] | 余建星, 李红涛. 高架索海上补给装置在小型船舶补给上的应用[J]. 海洋技术, 2005, 24(2):59-62. |
[9] | 施平安, 何立居, 叶家玮. 两船波浪中并靠补给相对运动模型试验研究[J]. 船海工程, 2012, 41(5):26-30. |