1. 引言
码头是供船舶停靠的水工建筑物,是港口的重要的组成部分。随着我国经济快速稳定发展,河道航运条件不断改善,船舶逐渐趋于大型化,同时一些较早建成的码头存在一定的老化问题,这对于码头的长久安全服役有着巨大挑战。在船舶停靠码头的时候不可避免的会产生撞击力作用在码头结构上,从而产生动力学响应,当产生的响应超过了码头的当前安全承载力,可能会发生不可预估的事故损失。对服役多年的码头进行安全评估是十分必要的,但仅仅将船舶撞击码头的过程看作是一个静力作用过程 [1],是不能完全保证结构的安全。目前高桩码头结构分析方法主要分为平面计算方法、空间计算方法和物理模型试验方法,考虑到码头的实际情况,本文根据相关规范中的参数建立详细合理的船舶和码头的三维模型 [2],通过模拟船撞事故来对码头进行动力学分析,对现有码头结构的安全评估以及后续高桩码头结构的设计都有一定的指导价值。本文采用软件LS-DYNA对船舶碰撞码头的动态过程进行有限元模拟,LS-DYNA主要采用拉格朗日算法,在结构撞击模拟中被广泛应用,十分适合用来分析船舶碰撞码头的响应状况 [3]。
2. 有限元计算模型的建立
2.1. 高桩码头结构模型
本文以镇江某高桩码头为研究对象建立有限元模型进行数值分析,码头的上部结构主要包括平台面板、横梁、纵梁、前边梁、系船梁及靠船构件,桩则包括直桩和斜桩。
由于码头受到船舶撞击后的动力响应受接岸结构等的影响相对较小,并且码头纵向长度长,各个结构分段间的相互约束较弱,选取一个结构分段作为独立空间受力体系进行分析 [4]。为了简化计算,选取码头的四跨五榀结构作为计算分析对象。同时考虑到码头结构构件较多,模型只保留主要构件,包括桩、横梁、纵梁、系船梁、平台面板以及撞击位置处的靠船构件橡胶护舷,略去了撞击点以外的靠船构件与橡胶护舷。同时为了减小计算量,采取刚度等效原则对码头结构构件的截面简化,同时以等效弹性模量模拟钢筋作用。
码头结构和橡胶护舷选用SOLID164单元模拟,均设置为弹性材料“MAT_ELASTIC”。码头结构材料参数:密度为2350 kg/m3,弹性模量为30 GPa,泊松比为0.167。橡胶护舷材料参数:密度为1150 kg/m3,弹性模量为8 MPa,泊松比为0.36。图1中红色部分为钢筋–混凝土结构,蓝色部分为撞击处的橡胶护舷。模型整体采用六面体单元进行扫略网格划分 [5],动力计算共包括149,917个节点,87,162个单元。建立的码头有限元模型见图1。
Figure 1. Finite element model of wharf structure
图1. 码头结构有限元模型
船舶对码头上部结构的撞击是一个短暂的受力过程,受岸坡土体的影响相对较小。因此参照地基土的实际情况,根据JTS 167-4-2012《港口工程基桩规范》 [6] 中的规定,计算得出码头桩基的受弯嵌固点深度为4.2 m。将固结位移边界条件施加在嵌固点以下的桩基,以此来代替土体对码头桩基的约束作用。
2.2. 船舶模型
研究选取5千吨级散货船为设计船型,由于不考虑船体变形,对船舶进行了简化,未细致模拟全船结构。船舶选用SHELL163单元 [7],设置为弹性材料“MAT_ELASTIC”,通过Mass质量点单元施加船体质量 [8]。船舶材料参数:密度为7800 kg/m3,弹性模量为210 GPa,泊松比为0.3。船舶的有限元模型见图2,共包括4599个节点和3170个单元。
在船舶与码头的碰撞过程中,周围的水体对船舶的运动起到了一定的阻碍作用,水体吸收了撞击船舶的部分动能。因此为了表征船舶因水的阻力而减少的部分动能,本文采用附加质量法,以0.03倍的排水量作为附加质量 [9]。
Figure 2. Finite element model of the ship
图2. 船舶有限元模型
3. 撞击仿真计算结果及其分析
3.1. 撞击力时程分析
失控船舶垂直撞击码头时,撞击区域在极短时间内会产生巨大的撞击力,造成被撞部位的损伤破坏,继而整个码头也会发生一定的变形和损伤。通过有限元计算得到的撞击力时程曲线见图3。
Figure 3. Time history curve of impact force (Unit: N)
图3. 撞击力时程曲线(单位:N)
图3中,红色曲线A为船舶与橡胶护舷的接触力,绿色曲线B为船舶与码头结构的接触力,蓝色曲线C为橡胶护舷与码头结构的接触力。在0.15 s左右时船舶首先开始接触护舷,撞击力不断增大,受到橡胶护舷的吸能缓冲作用,此时传递给码头结构的撞击力约为护舷受力的三分之一。随着船舶继续撞击,在1.32 s左右时橡胶护舷失效,船舶开始接触码头结构,此时撞击力迅速增大,达到峰值3914.1 kN后又迅速减小,在3.42 s左右时撞击力降低为0。碰撞过程持续时间约为3.27 s,撞击力曲线图形为抛物线形。
3.2. 码头结构位移时程分析
为分析码头在受到船舶撞击后,从左侧至右侧的位移变化,在码头前端设置5个追踪单元如图4所示,从左至右依次为A至E点。A点为追踪单元6049,B点为追踪单元10,838,C点为追踪单元26,947,D点为追踪单元16,785,E点为追踪单元21,580。
Figure 4. Tracking elements of front end
图4. 码头前端追踪单元示意图
通过图5(a)的横向位移时程曲线可以看出,在船舶开始接触橡胶护舷时,码头在横向上产生微小的位移。在橡胶护舷失效、码头开始受到撞击后,位移大小明显增大,同时伴随着显著的扭转,从左至右追踪单元的横向位移逐渐增大。随后可以观察到码头在沿着船舶入射方向运动的同时,发生了往复的扭转振动。由图5(b)纵向位移时程曲线可以看出码头前端整体性较好,在整个撞击过程中没有较大的位移差异。
为分析码头从前端至后端的位移变化,在码头的右侧上部结构设置9个追踪单元如图6所示。从左至右依次为A至H点。A点为追踪单元54,080,B点为追踪单元52,388,C点为追踪单元50,696,D点为追踪单元49,004,E点为追踪单元47,312,F点为追踪单元45,620,G点为追踪单元43,928,H点为追踪单元42,236。
(a)(b)
Figure 5. Displacement time history curve of front end (Unit: m)
图5. 码头前端位移时程曲线(单位:m)
Figure 6. Tracking elements of right end
图6. 码头右侧追踪单元示意图
码头右侧横向位移时程曲线见图7(a),可以看出码头离岸侧和近岸侧位移规律基本一致,曲线波形与图5(a)一致,其上波动为码头结构的扭转振动所致。纵向位移时程曲线见图7(b),可以明显发现码头受撞击后发生的扭转振动反映在曲线上的波动情况,离岸侧和近岸侧在扭转部分的位移规律差异明显,呈现出几乎完全相反的形式。其原因为梁板式高桩码头的梁系结构相互牵制作用大,码头整体性较好。
(a)(b)
Figure 7. Displacement time history curve of right end (Unit: m)
图7. 码头右侧位移时程曲线(单位:m)
3.3. 码头位移方向变化分析
位移矢量图能够完整反映碰撞过程中码头结构位移的大小及其方向的变化。码头结构在撞击力作用下,在8个不同时刻的位移矢量图如图8所示。
Figure 8. Displacement vector diagram of at different time (Unit: m)
图8. 不同时刻码头结构位移矢量图(单位:m)
船舶与码头初步接触后,在1.8 s时可以发现撞击点处靠船构件的位移矢量最大,位移方向与船舶的入射方向一致。同时可以看到码头近岸侧开始出现位移变化,靠近撞击点的位移明显大于远离撞击点处的位移,即码头结构在右侧开始受到撞击时,出现了绕左侧离岸端逆时针旋转的趋势。2.1 s至2.6 s的位移矢量图可以明显观察到码头结构在平动同时,伴随着显著的扭转现象,这期间撞击力达到峰值,可以看到码头结构江侧远离撞击点位置的位移变化不明显,码头绕左侧逆时针扭转的程度显著增大。随后码头结构的位移继续增大,在4.4 s左右时码头结构的位移达到峰值。5.4 s时船舶已与码头分离,码头位移值逐渐减小,同时码头开始绕右侧顺时针扭转,且扭转中心逐渐向右侧过渡。随后码头将在扭转一定程度后回正,持续重复这一旋转过程直至停止。这是由于船舶脱离码头后,码头还具有动能,因此码头将通过不断振动直至能量完全耗散后停止。
3.4. 码头结构等效应力分析
通过等效应力云图可以看出船舶撞击码头结构的过程中应力的传递过程。码头受到撞击后不同时刻的等效应力云图见图9。
Figure 9. Equivalent stress nephogram of the wharf (Unit: Pa)
图9. 码头结构等效应力云图(单位:Pa)
0.87 s时刻为船舶与护舷接触阶段,这期间应力主要集中在撞击位置附近的构件与近岸侧的两根斜桩中,此时等效应力值较小,护舷吸收了大部分的能量。1.68 s时橡胶护舷已经失效,码头受到的撞击力开始迅速增大,应力此时分布于离岸侧直桩,在1.71 s时已迅速由被撞排架的横梁传递至中部的桩。2.07 s时撞击力达到峰值,应力传递至离岸侧斜桩。在1.68 s至2.07 s撞击力迅速增大的过程中,应力峰值在0.4 s内迅速由8 MPa增大到22 MPa。由2.31 s和2.7 s时码头上部的梁结构的应力云图可以看出,在撞击力到达峰值至船舶脱离码头期间,应力仍不断从撞击处通过前端梁和纵梁传递至其他排架。船舶脱离码头后,码头结构仍将进行持续的伴随扭转的振动,这种运动形式导致桩与上部结构的连接处会出现应力集中的现象,等效应力的分布将在3.42 s与4.17 s的两种分布状态间不断切换,此时离岸侧斜桩桩顶的等效应力峰值已接近40 MPa,远超过混凝土强度设计值C30。
通过对码头结构受撞击过程的等效应力云图进行分析,可以发现码头结构在受到撞击时应力主要由离岸侧桩斜度较大的斜桩承担,此处为码头的薄弱位置。码头在受到撞击后的振动过程中,应力主要分布于桩与横梁连接处。因此在发生船撞事故后,除更换护舷、维修靠船构件外,应重点对撞击点附近的排架中的斜桩和各桩与上部结构连接处进行质量检测。
3.5. 码头结构加速度时程分析
图10为码头结构在受到撞击过程中横向和纵向的加速度–时间关系图,其中红色曲线A为码头横向加速度时程曲线,绿色曲线B为码头纵向加速度时程曲线。
Figure 10. Acceleration time history curve of wharf structure (Unit: m/s2)
图10. 码头结构加速度时程曲线(单位:m/s2)
通过图10可知,码头结构在受到撞击后,加速度时程曲线的振幅不断扩大,且在1.5 s左右护舷失效后陡增,说明振动随着碰撞的发生愈演愈烈。在船舶与码头分离后横向振动减弱,加速度幅值开始逐渐降低。横向和纵向两个方向的加速度时程曲线呈现出相同的规律性,这说明结构在这两个方向上的振动形式一致,曲线形状类似于正弦函数。
4. 结论
本文通过对失控船舶撞击高桩码头的动态过程进行数值模拟分析,得到以下结论:1) 码头上部结构受到撞击后,码头最高点处的位移最大,位移自上至下逐渐减小。不同排架之间,受撞击排架内的应力远大于其他排架;2) 码头受到撞击后会出现较为明显的扭转,撞击侧在撞击力作用下会绕着另一侧转动,在船舶脱离后逐渐恢复至最初位置。由于码头仍具有动能,码头将不断振动直到能量完全衰减后停止;3) 码头应力最开始产生于撞击位置,应力通过系船梁、前端梁等构件传递至其他横梁,再通过横梁传递至桩和纵梁。在撞击过程中,应力主要由码头离岸侧的斜桩承担。在船舶脱离后码头不断振动的过程中,桩与上部结构的连接位置极为薄弱。当船撞事故发生后,应重点对这些位置进行质量检测。
NOTES
*通讯作者。