水下爆炸气泡作用下船模的动态响应.docxVIP

水下爆炸气泡作用下船模的动态响应 这是一个非常复杂的问题，包括模拟水下爆炸负荷、爆炸气泡的动态性能、固结界面分析、水下结构的非线性动态响应等。目前,对于这类三维结构动态问题的主要研究方法是边界元/有限元结合的方法。主要采用流固耦合方程来表征耦合界面上的载荷,而不需要对水域进行建模,因此减少了大量的计算量。这种方法主要经历了平面波理论,曲面波理论,虚质量理论到双渐近(DAA)理论的发展过程。特别是由Geers所提出的DAA理论,综合并改进了前面几种方法,在高、中、低频段都有较好的精度。DAA方法经过众多学者们的改进与发展,目前是一种解决流固耦合问题的非常有效的计算方法。 水下爆炸中对舰船主要产生两种载荷:一是瞬态的冲击波,其压力峰值很高但持续时间很短,通常造成水中结构的局部破坏;二是气泡脉动载荷,气泡中大约含有47%的爆炸能量,具有非常强的破坏力。这一阶段的明显特征是达到的压力峰值较冲击波的压力峰值低,而且持续时间很长,约为几百毫秒,会造成水中结构的总体响应和局部响应。 对于气泡的动力学,众多学者已经进行了大量的研究工作[13,14,15,16,17,18,19]。对于气泡载荷作用下的舰船的动态响应,目前的主要研究方法都是将舰船简化为船体梁,研究船体梁的总体鞭状响应[20,21,22,23,24,25,26]。然而,对于气泡作用下三维全船结构的动态响应的细节分析,开展的研究还很少。 本文致力于建立一套有限元方法与DAA方法相结合的计算程序,用以研究水下爆炸气泡作用下水面舰船的动态响应。首先,建立一个考虑迁移效应,自由面效应和气泡阻力的气泡模型和一个水面舰船结构模型,并阐述了流固耦合与结构响应相关的理论分析;然后以该船模作为算例,研究了气泡作用下船体模型的总体响应和局部响应。比较了不同位置的加速度,速度和位移时程曲线。最后详细讨论了气泡作用下船体模型动态响应的特征与机理。 1 理论背景 1.1 流固耦合方程的建立 一个弹性结构受到外部激励的运动方程可以表示为: Μs??x+Cs˙x+Κsx=F(1)Msx??+Csx˙+Ksx=F(1) 其中:Ms,Cs,Ks分别为N×N阶的质量矩阵,阻尼矩阵和刚度矩阵。N×1阶向量x,˙x,??xx,x˙,x??分别为结构的位移,速度和加速度。列向量F代表外力。N为结构的自由度总数量。 对于浸没在流体的结构所受的外部激励可以表示为: F=-GAf(pi+ps) (2) 式中:pi和ps分别为湿表面上入射流和辐射流作用下产生的节点压力向量,Af为湿表面单元的面积对角矩阵,G为结构节点力和湿表面节点力的坐标转换矩阵。 当一个水中结构如水面舰船或潜艇浸没在无限声学流体介质中,结构湿表面的壳单元的控制方程可以由DAA方法来表示。对于流体中的弹性结构,DAA方法表示的结构表面流体的运动即为各正交流体边界模态的线性组合。一阶双渐近DAA方程的矩阵形式表达为: Μf˙ps+ρcAfps=ρc˙us(3)Mfp˙s+ρcAfps=ρcu˙s(3) 式中:Mf为N×N阶的流体质量矩阵;ρ和c分别为流体密度和流体中声速;us是辐射流的流体法向速度向量。 根据在流固耦合湿表面上,结构和流体的法向速度相等的条件,得到: GΤ˙x=ui+us(4)GTx˙=ui+us(4) 式中:ui为入射流的流体法向速度向量。 将式(2)代入式(1),式(4)代入式(3),得到了下面的流固耦合方程组: Μs??x+Cs˙x+Κsx=-GAf(pi+ps)(5)Msx??+Csx˙+Ksx=?GAf(pi+ps)(5) Μf˙ps+ρcAfps=ρc(GΤ??x-˙ui)(6)Mfp˙s+ρcAfps=ρc(GTx???u˙i)(6) 为了求解的方便,将式(5)和式(6)改写成对称形式。令ps=GΤ˙y?ps=GTy˙?这样式(5)和式(6)可以改写成: Μs??x+Cs˙x+Κsx=-GAf(pi+GΤ˙y)(7) ΜfGΤ??y+ρcAfGΤ˙y=ρc(GΤ??x-˙ui)(8) 将式(8)左右两边同乘以GAfM-1f,得到: GAfGΤ??y+ρcGAfΜ-1fAfGΤ˙y=ρcGAf(GΤ??x-˙ui)(9) 这样式(7)和式(9)就可以联立进行数值耦合求解。 1.2 阶龙格库塔法 冲击波之后,大约47%的爆炸能量仍然存留在气泡中。假设流体为无旋,不可压缩的,气泡中心在自由水面下深度d处。由于本文考虑的是中远场爆炸,假设船体对气泡没有影响,并且气泡在运动过程中保持球形。为考虑气泡的迁移效应,自由面效应和气泡阻力的影响,采用Vernon的无量纲方程组: σ=˙x?(10a) λ=˙ζ?(10b) ˙σ=-3δ(2δ-βx)[σ2x(1-2βx3δ)-λ26x+ζxζ0-(γ-1)kx3γ+1+βx4δ2(Cdλ2