收稿日期: 2020-04-24.
作者简介: 董仁义(1979-),男,博士,工程师,研究方向为船舶与海洋结构物设计与制造
Research on displacement of ship underwater pipeline induced by multiple loads
China Ship Development and Design Center,Wuhan 430064, China
0 引言 舰船管系承担着输送各种流体介质的任务,特别是舰船水下管系,在设计时受到船体内部空间的限制,具有较为复杂的空间特性且距离长跨度较大,不同介质的管路往往布置较为紧凑。另外,不同介质的管系几何尺寸、安装固定以及工作压力、使用工况各不相同,介质管系往往受到多个负载的同时作用。
多负载一方面是来自管系的固有物理特性,比如管路质量产生的重力荷载、输送介质的温度所导致温度荷载、水下管系连同舷外具备初始压力的管系内压力荷载;另一方面是来自管路的使用工况。比如文献[1]中,介绍了舰船液压、疏水、海水冷却等管系,由于流体流速较高、流量较大,在进行管路流体截至和开启过程中,改变了流体的方向或速度,在惯性作用下,产生流体瞬变现象。流体瞬变产生的压力增值荷载和时程压力荷载,会使管道发生的剧烈强迫振动,使管系的位移产生变化。这不仅会与其他管系产生碰撞干涉,损伤管系的防腐、保温、隔振层,更严重的是会造成管路与支撑结构表面产生微动磨损,使支撑刚度逐渐下降从而降低管系的固有频率,当其与管系泵组的压力脉动频率接近或者重合,就产生流固耦合振动,进而影响整个管系的性能和降低设备的使用寿命,甚至可能导致舰船管系功能失效。
在对单一要素比如弯头[2-3]、单一管道位移分析[4-5]、单一舰船系统仿真分析[6-7]和流体瞬变对管路激振分析[1]的基础上,结合文献[8]组合载荷对舰船管路的静动态分析,本文建立一个典型的具有复杂空间走向和安装固定特性的舰船水下管系,对不同载荷作用下的管系位移进行了仿真研究,分析了多荷载对舰船管系的影响。
1 舰船管系荷载类型及模型 舰船管系可能承受的荷载类型有,重力荷载:包括管道自重、防腐保温阻尼层重、介质重等;位移荷载:包括管道热胀冷缩位移、端点附加位移、支撑沉降等;压力载荷:包括内压力和外压力;瞬变流载荷:如安全阀起跳或阀门快速启闭时的压力冲击;两相流脉动载荷;压力脉动荷载;机械振动荷载等。文献[1]进行了的舰船水系统水锤特性仿真与试验研究,都是针对流体瞬变引起的动荷载。
本文选取的舰船管系是舰船水下连通舷外的疏水管路,在注水过程中具有重力荷载、热胀冷缩位移荷载、压力(内压力)荷载,在疏水管系的应变分析中,这几种属于静荷载。另外还有瞬变流荷载,是偶然性动荷载,其瞬时压力值(Pressure transient)和压力增值(Pressure Rise)采用下式计算:
$ \Delta p = \rho \cdot a \cdot \Delta \nu{\text{。}} $
|
式中:
$ \Delta p$
为管路瞬时压力值和压力增值;
$ \rho $
为流体密度;
$ a$
为流体内声速;
$ \Delta \nu$
为阀关闭后流体流速的变化值。另外,流体内声速由下式来计算:
$ a = \sqrt {\dfrac{{K/\rho }}{{1 + (K/E)\left( {D/e} \right)}}} {\text{。}} $
|
其中:
$\ \rho $
为流体密度;
$ K$
为体积模量;
$ E$
为管材的弹性模量;
$ D$
为管路直径;
$ e$
为管路壁厚。当
$ \Delta p$
作为压力增值的时候,本文将其用于管系增加的内压力荷载。
2 管系仿真模型 2.1 管系描述 某段典型的舰船空间管系的走向如图1所示,管系模型各个节点的编号和坐标如表1所示。管路公称通径300 mm,实际外径为323.85 mm,壁厚为9.525 mm,最小屈服强度为206.84 N/mm2,管材密度7833 kg/m3。A00点和A18点固定,流体假定为海水,密度为1028 kg/m3,由A18点流入,沿管路至A00点,同时在图2中标注出各管点的标号和管路的长度。在管系点A00处有一管路截至阀,A00点截止阀突然关闭将产生沿管线传递的压力冲击波。
表 1(Tab. 1)
表 1 典型舰船空间管系算例节点坐标
Tab. 1 Point name and coordinate of space pipeline
节点代号 |
坐标/mm |
|
节点代号 |
坐标/mm |
X
|
Y
|
Z
|
X
|
Y
|
Z
|
A00 |
0.00 |
0.00 |
0.00 |
A09 |
6000.00 |
2000.00 |
−24000.00 |
A01 |
0.00 |
0.00 |
−6000.00 |
A10 N |
6000.00 |
2000.00 |
−25542.80 |
A02 N |
0.00 |
0.00 |
−6542.80 |
A10 |
6000.00 |
2000.00 |
−26000.00 |
A02 |
0.00 |
0.00 |
−7000.00 |
A10 F |
6000.00 |
1542.80 |
−26000.00 |
A02 F |
457.20 |
0.00 |
−7000.00 |
A11 N |
6000.00 |
−542.80 |
−26000.00 |
A03 N |
5542.80 |
0.00 |
−7000.00 |
A11 |
6000.00 |
−1000.00 |
−26000.00 |
A03 |
6000.00 |
0.00 |
−7000.00 |
A11 F |
6000.00 |
−1000.00 |
−26457.20 |
A03 F |
6000.00 |
0.00 |
−7457.20 |
A12 |
6000.00 |
−1000.00 |
−28000.00 |
A04 |
6000.00 |
0.00 |
−8000.00 |
A13 |
6000.00 |
−1000.00 |
−34000.00 |
A05 |
6000.00 |
0.00 |
−14000.00 |
A14 |
6000.00 |
−1000.00 |
−40000.00 |
A06 N |
6000.00 |
0.00 |
−15542.80 |
A15 |
6000.00 |
−1000.00 |
−46000.00 |
A06 |
6000.00 |
0.00 |
−16000.00 |
A16 N |
6000.00 |
−1000.00 |
−51542.80 |
A06 F |
6000.00 |
457.20 |
−16000.00 |
A16 |
6000.00 |
−1000.00 |
−52000.00 |
A07 N |
6000.00 |
1542.80 |
−16000.00 |
A16 F |
6457.20 |
−1000.00 |
−52000.00 |
A07 |
6000.00 |
2000.00 |
−16000.00 |
A17 |
7000.00 |
−1000.00 |
−52000.00 |
A07 F |
6000.00 |
2000.00 |
−16457.20 |
A18 |
13000.00 |
−1000.00 |
−52000.00 |
A08 |
6000.00 |
2000.00 |
−18000.00 |
|
|
|
|
|
表 1 典型舰船空间管系算例节点坐标
Tab.1 Point name and coordinate of space pipeline
|
2.2 管系固定方式 图2中结合舰船管系的特点,给出了管系的固定方式,包含固定支座和移动支座,刚性支座和弹性支座。节点A00和A18为固定座,A01,A04,A05,A08,A09,A13,A15,A17为垂向限位支架(vertical stop),限制管路上下移动与管路上下的间隙为0;A08和A09是弹簧吊架(spring hanger),A12和A14是导向支架(guide),限制管路的径向运动与管路上下左右的间隙为0,允许管路的轴向移动。
3 不同荷载对管系位移作用仿真 空间管系固定时,已将支座与管系间隙设置为0,并且不考虑摩擦力的作用,进行静力荷载分析时忽略了系统中所有非线性约束,进行静荷载的线性分析。本文的静荷载包括重力荷载(GR)、温度荷载和内压力荷载。在本算例中,由于在舰船疏水管路工作时,介质的温度变化不大,并且初始设定为20 ℃,由温度变化引起的热胀冷缩导致的管系位移变化非常小,也可以忽略。初始的内压力载荷设定为水下10 m,管系内海水静压为0.1 MPa。当关闭A00处的管系阀门,瞬变流产生压力增值
$ \Delta p$
,在确定内压力荷载时,考虑压力值P2为
$ \Delta p$
加上初始静压0.1 MPa后的值,流体瞬变产生的动荷载的时程分析过程中的位移用代号M1表示。
管路的位移特性分析如图3~图7所示,前8阶振动频率见表2所示。
表 2(Tab. 2)
表 2 管系模态频率计算结果
Tab. 2 Modal frequency result of pipeline
阶数 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
振动频率/Hz |
1.0066 |
1.3830 |
2.2020 |
3.2012 |
3.9077 |
4.8575 |
7.2185 |
8.0470 |
|
表 2 管系模态频率计算结果
Tab.2 Modal frequency result of pipeline
|
表3给出了管系节点在不同荷载下的沿坐标轴位移值,GR+P2和GR+M1是为了方便对比分析所做的代表性的位移叠加,分别表示最大的静力荷载所产生的位移,重力载荷与动荷载所产生的位移。
表 3(Tab. 3)
表 3 管系节点在多荷载下的位移
Tab. 3 Pipeline displacement under multiple loads
节点代号 |
荷载 组合 |
位移/mm |
|
节点代号 |
荷载 组合 |
位移/mm |
X
|
Y
|
Z
|
X
|
Y
|
Z
|
A00 |
GR
|
0.00 |
0.00 |
0.00 |
A09 |
GR
|
−1.23 |
0.00 |
−0.02 |
|
P2 |
0.00 |
0.00 |
0.00 |
|
P2 |
0.49 |
0.00 |
0.81 |
|
M1 |
0.00 |
0.00 |
0.00 |
|
M1 |
32.32 |
0.00 |
46.64 |
|
GR+P2 |
0.00 |
0.00 |
0.00 |
|
GR+P2 |
−0.74 |
0.00 |
0.78 |
|
GR+M1 |
0.00 |
0.00 |
0.00 |
|
GR+M1 |
31.09 |
0.00 |
46.62 |
A01 |
GR
|
−0.01 |
0.00 |
0.00 |
A10N |
GR
|
−1.23 |
−0.33 |
−0.02 |
|
P2 |
0.37 |
0.00 |
−0.33 |
|
P2 |
0.30 |
0.08 |
0.72 |
|
M1 |
10.47 |
0.00 |
0.02 |
|
M1 |
27.85 |
0.16 |
46.64 |
|
GR+P2 |
0.36 |
0.00 |
−0.33 |
|
GR+P2 |
−0.93 |
−0.25 |
0.70 |
|
GR+M1 |
10.46 |
0.00 |
0.02 |
|
GR+M1 |
26.62 |
−0.17 |
46.62 |
A02N |
GR
|
−0.01 |
−0.02 |
0.00 |
A10F |
GR
|
−1.03 |
−0.43 |
0.04 |
|
P2 |
0.45 |
0.00 |
−0.36 |
|
P2 |
0.23 |
0.07 |
0.70 |
|
M1 |
12.39 |
0.07 |
0.02 |
|
M1 |
23.19 |
0.14 |
46.39 |
|
GR+P2 |
0.44 |
−0.02 |
−0.36 |
|
GR+P2 |
−0.80 |
−0.36 |
0.74 |
|
GR+M1 |
12.38 |
0.05 |
0.02 |
|
GR+M1 |
22.17 |
−0.29 |
46.43 |
A02 F |
GR
|
−0.01 |
−0.31 |
0.00 |
A11N |
GR
|
−0.17 |
−0.43 |
0.05 |
|
P2 |
0.57 |
0.00 |
−0.26 |
|
P2 |
0.10 |
−0.05 |
0.76 |
|
M1 |
14.65 |
0.27 |
2.80 |
|
M1 |
7.27 |
0.14 |
44.54 |
|
GR+P2 |
0.56 |
−0.31 |
−0.26 |
|
GR+P2 |
−0.07 |
−0.48 |
0.80 |
|
GR+M1 |
14.65 |
−0.04 |
2.80 |
|
GR+M1 |
7.10 |
−0.28 |
44.59 |
A03N |
GR
|
−0.01 |
−0.32 |
−0.01 |
A11F |
GR
|
0.01 |
−0.33 |
−0.02 |
|
P2 |
0.85 |
−0.01 |
1.60 |
|
P2 |
0.05 |
−0.06 |
0.73 |
|
M1 |
14.66 |
0.29 |
43.41 |
|
M1 |
2.92 |
0.26 |
44.31 |
|
GR+P2 |
0.84 |
−0.33 |
1.59 |
|
GR+P2 |
0.06 |
−0.40 |
0.71 |
|
GR+M1 |
14.65 |
−0.03 |
43.40 |
|
GR+M1 |
2.93 |
−0.08 |
44.29 |
A03 F |
GR
|
−0.01 |
−0.01 |
−0.01 |
A12 |
GR
|
0.00 |
0.00 |
−0.02 |
|
P2 |
0.98 |
0.00 |
1.70 |
|
P2 |
0.00 |
0.00 |
0.65 |
|
M1 |
17.22 |
0.03 |
46.29 |
|
M1 |
0.00 |
0.00 |
44.30 |
|
GR+P2 |
0.97 |
−0.02 |
1.69 |
|
GR+P2 |
0.00 |
0.00 |
0.63 |
|
GR+M1 |
17.21 |
0.02 |
46.28 |
|
GR+M1 |
0.00 |
0.00 |
44.28 |
A04 |
GR
|
−0.02 |
0.00 |
−0.01 |
A13 |
GR
|
−0.01 |
0.00 |
−0.02 |
|
P2 |
1.07 |
0.00 |
1.67 |
|
P2 |
0.00 |
0.00 |
0.32 |
|
M1 |
19.64 |
0.00 |
46.29 |
|
M1 |
3.99 |
0.00 |
44.27 |
|
GR+P2 |
1.04 |
0.00 |
1.66 |
|
GR+P2 |
−0.01 |
0.00 |
0.30 |
|
GR+M1 |
19.62 |
0.00 |
46.28 |
|
GR+M1 |
3.98 |
0.00 |
44.25 |
A05 |
GR
|
−0.12 |
0.00 |
−0.01 |
A14 |
GR
|
0.00 |
0.00 |
−0.02 |
|
P2 |
1.50 |
0.00 |
1.34 |
|
P2 |
0.00 |
0.00 |
−0.01 |
|
M1 |
38.27 |
0.00 |
46.29 |
|
M1 |
0.00 |
0.00 |
44.24 |
|
GR+P2 |
1.38 |
0.00 |
1.33 |
|
GR+P2 |
0.00 |
0.00 |
−0.03 |
|
GR+M1 |
38.15 |
0.00 |
46.29 |
|
GR+M1 |
0.00 |
0.00 |
44.22 |
A06N |
GR
|
−0.15 |
−0.19 |
−0.01 |
A15 |
GR
|
0.00 |
0.00 |
−0.02 |
|
P2 |
1.50 |
−0.04 |
1.26 |
|
P2 |
−0.26 |
0.00 |
−0.34 |
|
M1 |
40.78 |
0.14 |
46.30 |
|
M1 |
7.35 |
0.00 |
44.20 |
|
GR+P2 |
1.35 |
−0.23 |
1.25 |
|
GR+P2 |
−0.26 |
0.00 |
−0.36 |
|
GR+M1 |
40.63 |
−0.05 |
46.29 |
|
GR+M1 |
7.34 |
0.00 |
44.18 |
A06 F |
GR
|
−0.41 |
−0.25 |
−0.05 |
A16N |
GR
|
0.00 |
−0.22 |
−0.02 |
|
P2 |
1.49 |
−0.02 |
1.23 |
|
P2 |
−0.41 |
0.00 |
−0.65 |
|
M1 |
42.02 |
0.12 |
46.35 |
|
M1 |
2.28 |
0.43 |
44.16 |
|
GR+P2 |
1.08 |
−0.27 |
1.18 |
|
GR+P2 |
−0.41 |
−0.22 |
−0.67 |
|
GR+M1 |
41.61 |
−0.12 |
46.30 |
|
GR+M1 |
2.28 |
0.21 |
44.14 |
A07N |
GR
|
−0.99 |
−0.25 |
−0.06 |
A16F |
GR
|
0.00 |
−0.02 |
−0.02 |
|
P2 |
1.48 |
0.04 |
1.25 |
|
P2 |
−0.36 |
0.00 |
−0.64 |
|
M1 |
45.40 |
0.12 |
46.58 |
|
M1 |
0.01 |
0.17 |
41.18 |
|
GR+P2 |
0.50 |
−0.21 |
1.19 |
|
GR+P2 |
−0.36 |
−0.02 |
−0.66 |
|
GR+M1 |
44.41 |
−0.13 |
46.52 |
|
GR+M1 |
0.01 |
0.15 |
41.17 |
A07 F |
GR
|
−1.23 |
−0.19 |
−0.02 |
A17 |
GR
|
0.00 |
0.00 |
−0.02 |
|
P2 |
1.45 |
0.06 |
1.22 |
|
P2 |
−0.33 |
0.00 |
−0.57 |
|
M1 |
47.12 |
0.11 |
46.65 |
|
M1 |
0.01 |
0.00 |
36.58 |
|
GR+P2 |
0.22 |
−0.14 |
1.20 |
|
GR+P2 |
−0.33 |
0.00 |
−0.59 |
|
GR+M1 |
45.90 |
−0.08 |
46.63 |
|
GR+M1 |
0.01 |
0.00 |
36.56 |
A08 |
GR
|
−1.23 |
0.00 |
−0.02 |
A18 |
GR
|
0.00 |
0.00 |
0.00 |
|
P2 |
1.27 |
0.00 |
1.14 |
|
P2 |
0.00 |
0.00 |
0.00 |
|
M1 |
45.49 |
0.00 |
46.65 |
|
M1 |
0.00 |
0.00 |
0.00 |
|
GR+P2 |
0.05 |
0.00 |
1.11 |
|
GR+P2 |
0.00 |
0.00 |
0.00 |
|
GR+M1 |
44.26 |
0.00 |
46.63 |
|
GR+M1 |
0.00 |
0.00 |
0.00 |
|
表 3 管系节点在多荷载下的位移
Tab.3 Pipeline displacement under multiple loads
|
表3中所示的管系各个节点的位移数据与图3~图7对应。节点A00和A18为固定座,图表中所示3个方向的位移均为0;A01,A04,A05,A08,A09,A13,A15,A17为垂向限位支架,图表中所示的Y方向的位移为0;A12和A14是导向支架,图表中所示X和Y方向的位移为0。限制管路的径向运动与管路上下左右的间隙为0,允许管路的轴向移动。在本算例中,可以看到流体瞬变所产生的动荷载是引起管路位移变化的最主要的因素,这与实际情况相符,证明对于管路振动特性的仿真可以作为管路设计参考。
4 结 语 舰船管系在工作过程中受到多荷载的作用,通过建立典型的空间管系,分析不同荷载对管路位移的影响。针对工作过程中介质温度不高,温差变化和初始工作压力较小的管系,温度变化导致的热胀冷缩和初始内压力荷载引起的位移变化相比动载荷引起的位移变化,可忽略其位移影响,管系位移主要是动载荷造成的。
舰船水下管系流体瞬变引起的管系压力增量远大于管系静水压力,有必要作为内压力静载荷进行位移影响分析,特别是其扩展应用至材料刚性较小、结构呈薄壁几何特性并且工作介质温度极高或温差范围大的工业管系。多荷载对舰船管系位移影响研究方法和结论,需要全面考虑对管系设计因素,结果对于设计舰船管系空间管路设计有一定的指导意义,后续可以进一步考虑其他因素,比如高温、高压、蒸汽介质、管系恒定或周期的受力,支座和管系的非线性分析等来丰富荷载类型。