24小时服务热线13370829909
青岛六和迅捷物流有限公司
青岛专业的大型设备运输公司
电 话:86-0532-55578775
电 话:86-0532-55578076
传 真:86-0532-55578773
手 机:013370829909(王先生)
邮 箱:andowang@qdliuhe.com
地 址:青岛市崂山区海尔路61号天宝国际银座1302/1408室
------集装箱桥吊海上运输主动力计算分析
廖建清
(中华人民共和国福州海事局,福州 350015)
摘 要 集装箱桥吊海上运输属特殊重大件货物海上运输。本文应用IMO《货物系固手册编制指南》的编制原理,结合集装箱桥吊海上运输实践,推导出集装箱桥吊海上运输过程中整体及重要部位受力分析公式,为这一特殊领域的货物系固和有关强度核算打下理论基础。
关键词 桥吊 海上运输 受力分析
1 引言
集装箱桥吊海上运输具有专业性强、技术要求高、风险大的特点,它需要专业的运输船舶或专业的运输驳船,专门的积载和系固,专业的技术人员,涉及船舶稳性和强度以及专业运输管理等方面技术,保证其足以抵御海上的一般风险,从而实现这种特殊的重大件货物的海上运输。
目前国内外在这方面的技术处于业务垄断和技术保密状态,理论研究分析不多。平常航海人员通常采用设计好的固定堆垛方式或用粗略估计、非语言式思维和模糊控制理论等方式处理一般货物组件的船上积载。但对于重大件非标准货物海上运输,一定要做到心中有数,一定要用数学思维、量化分析处理船舶稳性和强度等问题。在整个运输过程中,要考虑种种因素,预先制定出一整套运输方案,且在每一个运输环节还要反复考察论证方案的正确性。所以本文先从桥吊在船上积载和主动力计算两方面入手,归纳两点理论的主动力计算方法,为这一特殊运输领域的船舶稳性、强度核算、海上拖带、风险控制提供参考。
2 集装箱桥吊海上运输的积载
2.1 桥吊海上运输要求专门的船舶
首先船舶要有适当的尺度和足够的强度,以抵御海上的一般风险,其次船舶要有能够调节稳性和吃水系统,实现门对门运输。本文只讨论第一部分的主动力计算问题。此外船上要装备起重、顶压、气割、焊接等安装设备。
2.2 运输船舶尺度要求
一是宽度要基本与桥吊的横向门跨度相一致,或偏宽、偏窄误差在2~3m之内,长度则要大大于桥吊的纵向门跨度,甲板是连续甲板,船体的横梁和肋板要对应在桥吊侧腿支点的正下方,如果不刚好的话,在船体上要加焊强力肋骨和肋板,使之处于桥吊侧腿支点的正下方,以满足船体局部强度的要求。
2.3 桥吊积载考虑的因素
桥吊的横向跨度与甲板宽度大致相同,由于桥吊的自身重量很重,单件重量在千吨以上,且重量分布不一定均匀,分布在四个门腿的重量有几百吨,积载时以下因素要认真考虑:
(1)装卸是要靠外力使桥吊从岸上移动到船甲板上或从船甲板上移到岸上的预定位置,整个过程要平稳有序进行,首选的是在门腿下安装轨道,实行轨道滑动装卸。
(2)轨道下面要有足够的衬垫板,以保证船甲板的局部强度要求。
(3)保证桥吊的重心在船舶的中纵剖面上。所以在装载之前一定要有桥吊的有关资料,如主要尺度、重量、重心位置、受风面积、主要设备保护措施等资料。
(4)要事先选定装卸的舷侧,是右舷还是左舷,以便轨道的安装和预定的装卸操作走向。这就要求在制定运输方案之前要反复考察论证。
(5)如果甲板宽度不够,或者为调整重心位置而需门腿向一侧偏移时,导致一侧门腿超出甲板宽度,这时需在超出一侧横向加接受载底座,以满足装载要求。
(6)船体各舱要专门设计成可在一定范围内调节船舶重心高度、船舶吃水和防自由液面影响的隔舱,以满足受载后船舶稳性以及装卸时根据潮位调整甲板与码头平台的随时水平。
(7)位置较低的电机、接线箱和开关箱等电器设备,以及裸露的较精密的机械零件和设备均需用油布等包裹严实,避免海水侵蚀,造成损坏。
2.4 桥吊船上积载方法
2.4.1 轨道安装方法
轨道安装方法如图1,钢轨的铺设有横向和纵向,横向的轨距与门腿横跨度一致,误差应保持在10mm内,轨道下面铺垫强力钢板,钢板直接焊接在甲板上,轨道再焊接在强力钢板上,焊接要牢固,轨道之间用连接板连接定位,轨面不平时,用垫板垫平。
2.4.2 桥吊腿侧积载要求
由于船舶在海上航行会遇到风浪的危险,船舶摇摆是造成桥吊移动的最大危险,在积载上采取两腿侧轮相互垂直积载,以减少摇摆造成滑动的危险(图1为桥吊的腿轮可以90°调整和安装,安装工艺可分别用大型千斤顶顶离轨道,脱开构件连接螺栓,旋转90°后再对准连接固定,强度上不变)。
3 装箱桥吊在甲板上的主动力计算公式
受力分析是本文的中心内容,通过受力分析,了解桥吊相对于船甲板的相对滑动力和受压力变形趋势,为下一步的桥吊系固和强度计算提供最基础的理论依据。
3.1 风浪中桥吊在甲板上的主动力分析
海上运输静态情况很少,在有风浪小倾角的情况下,即船舶倾角小于10o~15o(甲板边缘开始入水情况下另外核算)情况下,甲板上积载的货物单元主要受到物体运动的惯性力、风的作用力和浪打的波溅力等三种主动力的作用,这是普遍现象。根据力的分解原理,三种主动力可以分解为纵向、横向和垂向的分力:
Fx = M·ax + Fwx + Fsx
Fy = M·ay + Fwy + Fsy
Fz = M·g ±M·az
式中,Fx-纵向主动力;
Fy-横向主动力;
Fz-垂向主动力;
M -货物单元的质量(重量,t);
ax、 ay 、az-货物单元据其在不同货位时纵向、横向和垂向加速度(m/s2);
Fw(xy)-货物单元所受到纵向或横向的风力(kN);
Fs(xy)-货物单元所受到纵向或横向波溅力。
对集装箱桥吊在海上运输而言,纵向主动力Fx受到纵摇加速度力M·ax、纵向的风力Fwx和纵向波溅力Fsx。由于桥吊主体构架较高,浪打不到那么高,波溅力Fsx可以忽略不计,所以纵向受力Fx只剩下加速度力M·ax和风力Fwx。即 Fx = M·ax + Fwx
横向主动力Fy受到横摇加速度力M·ay、横向的风力Fwy和横向波溅力Fsy,同理横向波溅力Fsy可以忽略不计,所以横向受力Fy只剩下加速度力M·ay和风力Fwy。即Fy = M·ay + Fwy
垂向受力Fz有自身重力M·g、波浪引起的上下垂荡运动的加速度力±M·az。但要说明的是,自身重力M·g永恒向下(设向下为+),波浪引起的上下垂荡运动的加速度力±M·az有时向上,有时则向下,由低位置向上做加速度运动时,桥吊对船甲板的反作用力为正,这是我们要考虑的主要问题,所以对甲板受力的程度最大的应考虑在船体横倾到一侧且同时由低位置转向上运动的瞬间受力最大状态。即Fz= M·g + M·az
综上所述,桥吊整体主动力计算公式为:
Fx = M·ax + Fwx
Fy = M·ay + Fwy ………(1)
Fz= M·g + M·az·
风的压力可估算:
Fwx = Pw·A wx
Fwy = Pw·A wy
式中,Pw -为估算风压强,取Pw =1(kN/m2 );
A wx-露天甲板上桥吊纵向受风面积(m2);
A wy-露天甲板上桥吊横向受风面积(m2);
A wx和A wy-厂家应提供的数据。
公式(1)是本文把集装箱桥吊海上运输当作整体的货物组件进行分析的重要结论。
3.2 几个加速度值的确定
桥吊在运输过程中必须保持刚性不变,所以根据刚性物体的运动学原理和转动定律,加速度:
式中,l-是桥吊重心距船舶横纵漂心的距离;
M-是船舶所受的合外力矩;
J-是船舶的转动惯量。
合外力矩M和转动惯量J与船舶的质量、形状、运动状态有关,所以桥吊的加速度a(分纵ax、横ay、垂az)与船舶的质量、形状、运动状态有关。
为知道纵ax、横ay、垂az的值,本文重点介绍应用IMO《货物系固手册编制指南》确定加速度值的原理,确定ax、 ay、 az值,所以简捷的确定加速度值方法为:
式中,a0x-桥吊在不同装载位置时纵向基本加速度值;
aoy-桥吊在不同装载位置时横向基本加速度值;
aoz-桥吊在不同装载位置时垂向基本加速度值;
k1-船速修正系数;
k2-稳性高度修正系数。
以上基本值和修正系数是IMO在标准条件下(无限航区、全年航行、连续航行25d、船舶长度为100m和船舶航速15kn)计算出aox 、aoy 、aoz,并以表格的形式确定具体数值提供给使用者,充分体现了不同位置摇摆产生的加速度值是不一样的,位置越高加速度值越大,离船中(漂心)越远加速度值越大,与实际完全吻合。修正系数k1的意义是当实际船舶与标准船型、船速不一样时的修正系数,也体现了船速越快,急转向产生的加速度值越大,船型越短急转向时回旋半径越小,离心加速度值越大。修正系数k2的意义是稳性高度修正值,体现了初稳性高度越大,船舶摇摆的剧烈程度越厉害,加速度值也越大。以上5个数值均可在各船旗国政府认可的每艘船舶《货物系固手册》中随时查到,是可查数据,而且考虑到最不利因素的影响,也就是说实际的数值比表列的数值要小,所以根据表列数值计算的结果是安全的。
4 主动力计算公式(1)和(2)结论分析
4.1 主动力对桥吊自身构件的影响
以上桥吊的主动力计算公式(1)和(2)是建立在把桥吊整体当作质点来分析的,这种分析方法可以使复杂问题简单化,只要桥吊整体刚性保持不变,这种按质点的分析方法是科学的,满足客观了解事物的要求。但在风浪中要保证桥吊整体刚性不变,就得要求对桥吊自身各个构件进行绑扎以满足海上运输的一般风险。同时本文分析的目的在于认识整体桥吊在海上运输中产生主动力情况,而对桥吊自身各个构件的受力不作详细分析。
4.2 桥吊主动力对船舶稳性、强度影响的定性定量分析
4.2.1 桥吊主动力对船舶稳性的影响
桥吊整体重量一般在1200t以上,且重心位置高离甲板30m以上,如果一船同时运送几台桥吊,那么对船舶载重、船舶初稳心高度GM0和船舶倾覆力矩Mh将产生重大影响。根据船舶原理,载一个用下面公式(3)计算初稳性高度变动值;当装载多个集装箱桥吊(∑Pi>10%Δ)时,稳心高度新值KM1及船舶重心高度新值KG1的计算式为:
GM1= KM1-KG1(m)
KG1=(Δ·KG0+∑Pi·Zi)/Δ1 (m)
ΔGM= GM1- GM0(m)
式中,KM1 -根据排水量新值Δ1(Δ1=Δ+∑Pi)由静水力曲线图重新查取;
Pi-每件桥吊的重量(t),装载时取(+),卸载时取(-);
Zi- 每件桥吊的相应重心距基线高度(m)。
下面举例分析:
例1:某船排水量D0为38806.18t(并非满载),船长200m,甲板宽度32m,型深14m,KM0=15.113m,KG0=12.73m,单台桥吊重量为1200t,桥吊重心高度距甲板36.9m,设桥吊装在0.6L处,船速12kn,横向受风Y方向受风面积1156.3m2,Y方向受风作用高度Zm=36.84m。求(1)装船后对初稳性高度影响;(2)在风浪中对动稳性的影响有多大。
解:(1)近似的计算方法是:在只装卸一个桥吊时,一般Pi<10%Δ,忽略了KM值的少量变化,把KM值看作不变,仅仅看作是船舶重心高度变化的结果。也就是可以假设,只装一个桥吊时,先将桥吊重量P装在船舶原重心G0位置上,再将桥吊重量P自G0位置垂直上移至实际装载位置Zi,桥吊的垂移距离为(Z G0-Z i ),船舶重心也由G0移至G1。根据平行力移动原理,船舶重心移动距离为G0 G1,即初稳性高度变动值ΔGM可按下式计算:
因为Pi)
或 …………(3)
计算得:ΔGM=1200×(-12.73+36.9+14)/(38806.18+1200)=1.1449(m)
由公式(3)计算得稳性高度降低约1.1449m,初稳性高度降为:
GM1= KM0-KG0 -ΔGM,
=15.113-12.73-1.1449
=1.2381m
也就是说装完一个桥吊后GM值降为1.2381m。
(2)对动稳性方面,横向受力由公式(1)和(.2)计算
横向主动力Fy得:
Fy= M·ay + Fwy
ay = aoy·k1·k2
经查 aoy =6.8 ,k1=0.56 ,k2=1
ay=6.8×0.56×1=3.8(m/s2)
那么Fy=1200×3.8+
1156.3×1=5716.3(kn)。
增加倾覆力矩M ah:
Mah= M·ay·(36.9 +型深)+Fwy·(Zm+型深)
=1200×3.8×(36.9+14)+1156.3 ×(36.84+14)
=290959.67(knom)
定量算出桥吊的增加倾覆力矩M ah,根据船舶的稳性报告书,计算风对船体上部产生的倾覆力矩Mw ,总的倾覆力矩Mh就可以计算出来。经过核算就可分析是否满足船舶动稳性要求,见动稳性曲线图(图2)。根据动稳性平衡原理,面积a=面积b,θ2要小于进水角。以上分析只说明桥吊运输对船舶稳性影响相当大,所有的量化计算结果均应满足《IMO稳性规则》要求,这里不作重点说明。
4.2.2 腿侧受力对船体局部强度的影响
从桥吊主动力计算公式(1)中,还说明力来自纵、横、垂三个方向,并随着船舶的纵摇、横摇和垂荡桥吊各腿侧的受力各不相同,但作为海上运输风险预测,必须考虑最大危险存在的情况。下面通过船舶横风横浪向右横倾15o时(小倾角中的最大角)桥吊受力在腿侧处分解的分析,知道腿侧支点处受压力有多大,从而进一步知道桥吊腿侧压力对船体局部强度的影响有多大。
考虑风力作用中心与重心不很接近,则分开分解后落到甲板,再叠加,可推导出以下公式(4),公式如下:
Zt右压=1/2{(F1+F2)·cosα+ F3ˊ·cosδ }
Zt左压=1/2{(F2ˊ-F1ˊ)·cosβ + F3·cosγ}
Yt各个滑=1/4{(F1+F2)·sinα+( F2ˊ-F1ˊ) · sinβ+ ………(4)
(F3·sinγ+ F3ˊ· sinδ)}
式中,Zt右压-右侧腿压力;
Zt左压-左侧腿压力;
Yt各个滑-每个侧腿的平均横向滑动力;
α-桥吊重心所在纵剖面与桥吊重心和右侧腿支点连线的夹角;
β-桥吊重心所在纵剖面与桥吊重心和左侧腿支点连线的夹角;
γ-横向风力作用中心所在纵剖面与风力作用中心和左侧腿支点连线的夹角;
δ-横向风力作用中心所在纵剖面与风力作用中心和右侧腿支点连线的夹角;
F1-Fz在桥吊重心和右侧腿支点连线上的分力;
F2-横向加速度力Fay在桥吊重心和右侧腿支点连线上的分力;
F1ˊ- Fz在桥吊重心和左侧腿支点连线上的分力;
F2ˊ-横向加速度力Fay在桥吊重心和左侧腿支点连线上的分力;
F3-横向风力Fwy在风力作用中心和左侧腿支点连线上的分力;
F3ˊ-横向风力Fwy在风力作用中心和右侧腿支点连线上的分力。
以上公式(4)是本文集装箱桥吊海上运输对腿侧支点处局部受力分析的结论。
本例题的参数在表1中列出,按公式(4)计算的结果见表2。
表1 求解参数