摘要: |
水中结构在遭受水下爆炸作用时,在各种程度上都会受到水下爆炸气泡载荷的破坏作用,气泡运动引起的滞后流以及脉动压力引起结构的总体破坏,气泡坍塌形成的射流,引起结构的局部破坏。近年来的研究表明,气泡载荷对水中结构物造成严重的损伤,气泡与自由表面、结构等相互作用成为目前国际上研究的热点,但仍有许多现象和本质没有被揭示。为此,本论文从基本现象入手,重点研究气泡在流场中的膨胀、坍塌、回弹、迁移以及射流形成等重要动力学行为,结合实验数据研究气泡与自由表面、复杂弹塑性结构的相互作用,旨在揭示气泡与复杂界面的耦合特性及其对结构的毁伤机理。
本文从理论分析、实验技术、数值方法等方面回顾了国内外关于气泡动力学特性研究进展。且重点综述边界积分法(BIM)及其相关技术在气泡动力学特性模拟中的应用与发展,回顾讨论气泡模拟的发展过程及其关键技术;最后综述了近边界气泡的研究方法及进展,并在上述基础上提出了一些尚需进一步解决的问题,为指导本文的后续工作奠定了基础。
本文基于势流理论分别建立了轴对称及三维气泡动力学模型,包括平面单元和曲面单元,用边界积分法求解,得到气泡及其它边界的变形及位置。在计算奇异积分时,给出了几种不同的方法,如利用Green公式的积分性质等通过重构双层位势的主值积分消除双层奇异积分的奇异性。再如,利用平面单元的几何特性消除奇异性,并对此进行了讨论,给出了计算引起的误差,为边界积分法中奇异积分的计算提供了参考。
在环状气泡阶段,本文在前人研究成果的基础上分别建立了轴对称涡面模型、轴对称涡环模型以及三维涡环气泡模型。对比不同模型的优缺点及其适用范围,总结了环状气泡模型建立的详细步骤及其必须注意的关键技术,且在气泡运动模拟过程中引入数值光顺技术及弹性网格技术,避免因网格扭曲而导致的数值发散。最终开发全套可视化轴对称及三维气泡动态特性模拟程序,该程序用于解答在各种设想情况下,有关气泡动态特性的问题。
基于本文开发的计算程序,分析了不同模型、不同单元类型的计算精度,以及网格、时间步等因素对计算结果的影响,验证了本文数值方法的稳定性。分别将轴对称模型及三维模型与Rayleigh-Plesset方程的精确解及实验数据、实验照片进行了对比分析,分析表明,计算结果与Rayleigh-Plesset模型及实验数据吻合很好。并分析了气泡在重力场中的运动特性,计算出气泡运动引起滞后流速度及压力随时间的变化。在此基础上,专门研究了浮力参数δ、强度参数ε、气体参数γ等特征参数的变化,对气泡的周期、射流速度的影响,此外,本文还研究了药量和水深对气泡动态行为的影响,并总结了规律。本文系统地研究了近自由面气泡的动态特性,包括近自由面气泡的环状回弹及自由面的水冢现象;距离参数γ,、浮力参数δ、强度参数ε对气泡动态特性的影响;多气泡与自由面耦合作用研究,包括同相气泡与异相气泡的相互作用研究等。在研究过程中,讨论了基于开尔文冲量理论的Blake准则的适用范围,并解释了Blake准则失效的原因。此外,本文还发现多气泡特别是异相气泡之间存在强烈的抑制效应。
从气泡与壁面相互作用的基本现象入手,基于本文开发的程序系统地研究了,当弱浮力与强Bjerknes力,强浮力与弱Bjerknes力以及浮力与Bjerknes力相等时,刚性壁面附近气泡的动力学特性,其中包括水平壁面及倾斜壁面,并将各种工况的计算结果与Blake准则进行对比分析及讨论,得出偏射流方向及壁面压力与气泡的特征参数有密切的关系。此外,本文还模拟和讨论了气泡与结构、自由面三者之间的强非线性耦合作用。
计及结构的弹塑性,将边界元法(BEM)与有限元法(FEM)耦合提出了气泡、弹塑性结构以及自由面三者之间,复杂耦合动力学问题的计算方法,并开发了全套三维计算程序,计算值与实验值之间误差在10%以内,且通过计算发现,气泡载荷具有不对称性,同时总结了滞后流的衰减及变化规律。应用本文开发的程序,对平板,圆柱筒等简单结构进行分析,计算了滞后流、脉动压力以及射流等载荷对结构的毁伤,包括有、无自由面等多种工况,计算结果表明,气泡载荷对结构造成巨大的毁伤,特殊工况威力更大。特别地,以水面舰船为例,将三维计算程序工程化,应用于复杂弹塑性结构,最终形成一套完整的适合于工程应用的,三维气泡与弹塑性结构耦合动力学分析程序,旨在为相关舰船气泡动力学研究提供参考。
|