[发明专利]水下爆炸作用下舰船整体弹塑性运动响应预测方法及系统

摘要

本发明公开了一种水下爆炸作用下舰船整体弹塑性运动响应预测方法及系统，方法包括步骤：S1、将舰船简化为等截面自由船体梁；S2、将水下爆炸作用下自由船体梁载荷压力划分为五个时间阶段；S3、求解自由船体梁在水下爆炸五个时间阶段内的运动位移；S4、求解船体梁的弯矩；S5、求解船体梁正向塑性运动参数；S6、求解船体梁反向弹性卸载运动参数；S7、求解船体梁反向塑性运动参数；S8、求解船体梁正向弹性运动参数。本发明还提供了实现上述方法的系统。本发明在综合考虑冲击波和气泡载荷对舰船总体的冲击作用的基础上分析舰船塑性运动过程，可以准确且简便地实现对水下爆炸作用下舰船整体弹塑性运动变形的工程预报。

主权项

一种水下爆炸作用下舰船整体弹塑性运动响应预测方法，其特征在于，所述方法包括步骤：S1、将舰船等效为等截面自由船体梁，由此将对舰船运动响应的预测等效为对自由船体梁运动响应的预测；S2、将水下爆炸作用下自由船体梁载荷压力过程划分为五个时间阶段，采集冲击波压力峰值Pm、气泡脉动阶段负压峰值Pb、气泡第一次脉动压力峰值Ps；所述五个时间阶段为：0≤t＜t1、t1≤t＜t2、t2≤t＜t3、t3≤t＜t4、t4≤t＜t5；其中t1＝θ，θ为冲击波衰减常数；P0为炸药处静水压力，Patm为大气压，c为水中声速，r0为装药半径，R为爆距；me为装药当量，ρw为水的密度，g为重力加速度；t5=3290r0P00.71+t4;]]>S3、求解自由船体梁在水下爆炸五个时间阶段内的运动位移w(x,t)：其中，x为自由船体梁任意点处的横坐标值，k2＝t2‑t1，为自由船体梁第一阶固有振形，其中ζ1为振形幅值，μ1为梁运动频率参数，l为自由船体梁长度；为五个时间阶段的积分常数，ω1为自由梁的一阶振形固有频率；φ1为积分常数，m为考虑附连水质量的单位梁长度质量，p(x)为自由船体梁不同时间阶段压力分布特征函数，对于0≤t＜t1、t1≤t＜t2时间阶段，对于第t2≤t＜t3、t3≤t＜t4、t4≤t＜t5时间阶段，p(x)=(1-2x-ll)·exp[-8(2x-ll)2+4(2x-ll)3]+(2l-2xl)1.5·(2x-ll);]]>S4、根据船体梁在五个时间阶段的弹性运动位移w(x,t)，求解船体梁在不同阶段的速度和弯矩值；S5、当所述船体梁中点弯矩绝对值超过船体梁塑性极限弯矩Ms绝对值时，船体梁进入正向塑性运动，求解船体梁正向塑性运动时的运动幅值H'(t)和相对转动的转角α(t)：H··′(t)=3Ms-(lξ3+3ξ4)p(t)(lξ1+3ξ2)mα··(t)=24(ξ1ξ4-ξ2ξ3)p(t)-24Msξ1(lξ1+3ξ2)ml2,]]>ξ1＝λ1(l)‑λ1(l/2)，ξ2＝λ2(l)‑λ2(l/2)‑λ1(l)·l/2，ξ3＝η1(l/2)‑η1(l)，ξ4＝η2(l/2)‑η2(l)+η1(l)·l/2；∫0xp(x)dx=η1(x),∫0x∫0xp(x)dxdx=η2(x);]]>S6、当所述船体梁正向塑性变形达到最大即α(t)达到最大值时，其正向塑性变形阶段结束，进入反向弹性卸载，根据公式α(t)(x-l/2)+φ1(x)H′(t)=φ1(x)H(t)α·(t)(x-l/2)+φ1(x)H·′(t)=φ1(x)H·(t)]]>求解船体梁反向弹性卸载过程的运动幅值H(t)；S7、在船体梁反向弹性变形过程中，当梁中点弯矩绝对值超过塑性极限弯矩Ms绝对值时，船体梁进入反向塑性运动，求解梁进入反向塑性变形后的运动幅值H1'(t)和相对转动的转角α1(t)：H··1′(t)=3(-Ms)-(lξ3+3ξ4)p(t)(lξ1+3ξ2)mα··1(t)=24(ξ1ξ4-ξ2ξ3)p(t)-24(-Ms)ξ1(lξ1+3ξ2)ml2;]]>S8、当所述船体梁反向塑性变形达到最大即α1(t)达到最大值时，其反向塑性变形阶段结束，进入正向弹性卸载，根据α1(t)(x-l/2)+φ1(x)H1′(t)=φ1(x)H1(t)α·1(t)(x-l/2)+φ1(x)H·1′(t)=φ1(x)H·1(t)]]>求解船体梁正向弹性卸载阶段的运动幅值H1(t)。