水下爆炸近场非均熵流的特征线差分解法

提出了一种改进的特征线差分方法,根据力学基本方程推导出2族带有熵变修正项的特征线方程,并将质点迹线方程补充为第3族特征线方程,利用这3个特征线方程可以对水下爆炸冲击波及波后的非均熵流场进行求解。由于特征线差分法无需引入人工粘性项处理冲击波间断问题,所以减少了人为干预的影响因素。最后,应用这种改进的特征线法对质量分数比w(TNT)∶w(RDX)=40∶60的TNT/RDX炸药球的水下爆炸问题进行了数值模拟。在比例距离1R/R04的近场范围内,计算得到的冲击波峰值压力与文献实验数据的最大相对误差在13.1%以内,进一步证明了这种改进的特征线差分方法的正确性。     

柱状含铝炸药水下爆炸近场的特征线法研究

基于之前提出的一种含熵变项的特征线法,通过控制非等熵流中的能量释放来刻画铝粉燃烧的影响,结合简单Chapman-Jouguet模型和JWL-Miller状态方程,计算了柱形含铝炸药水下爆炸的近场参数。对比模拟结果与实验数据,发现这种特征线法可以较好地预测近场冲击波的传播迹线、爆轰产物的膨胀轨迹以及内部压缩波的反射过程。结果表明,这种特征线法可用于含铝炸药水下爆炸的近场计算,进一步可用于评估含铝炸药性能或计算水下能量输出。     

爆轰产物状态方程的水下爆炸反演理论研究

由于水下爆炸过程中爆轰产物的信息以水中压缩波形式向外传递,本文旨在探讨如何利用水下爆炸试验数据确定爆轰产物的状态方程。相较于常规圆筒试验,水下爆炸试验具有装置简单成本低、装药尺寸限制少、测定压力范围更广的特点,更适用于大药量或非理想炸药的现场测试。本文从水中冲击波轨迹和波后压力时程曲线出发,发展了由冲击波及其波后流场还原水气界面的逆特征线算法,以及根据水气界面确定爆轰产物状态方程的遗传算法。与水下爆炸正演结果对比,发现逆特征线法可以较准确地还原水气界面的位置和压力参数,有效压力下限可达2 MPa,远低于圆筒试验的测试下限0.1GPa。根据水气界面的反演结果,遗传算法也能稳定地优化JWL方程参数,8种常用炸药的等熵衰减压力误差在爆压至0.01 GPa的区间内都小于3%。结果表明,利用本文的逆特征线算法和遗传算法,理论上可以反演出压力范围较宽、准确性较高的爆轰产物状态方程。     

水下目标在水下爆炸作用下冲击响应的时频特征

基于某水下目标的抗水下爆炸实验数据,利用小波变换良好的时频局部化性质对监测的水下目标内部装置冲击加速度信号进行时频特征分析,得到了不同频带上冲击信号振动分量的时间历程曲线和能量分布。结果表明,基于冲击信号小波时频特征分析的频带能量可以同时反映冲击振动的强度、频率和持续时间对目标损伤的影响,同时结合不同频带的时间历程曲线可以获得冲击波和二次压力波在不同频带上的分布和衰减的细节信息,由此可确定冲击波和二次压力波各自对目标毁伤的影响。     

球对称形式物质点法及其在水下爆炸中的应用

利用物质点法求解三维远场水下爆炸冲击响应问题时,爆炸冲击波到达结构物之前传播过程的计算会消耗大量的计算资源,影响求解效率。针对这一问题,并考虑到水下爆炸冲击波在自由场中的传播特点,提出了一维球对称形式的物质点法,对球形炸药水下爆炸过程进行了数值模拟,并与Cole公式以及DYNA计算结果进行了比较,结果吻合较好,验证了一维球对称形式物质点法的正确性。在此基础上,提出了基于物质点法的重映射算法,将爆炸冲击波在自由场中的传播过程在一维球对称模型中进行计算,并将计算结果映射到三维模型中,以便继续完成计算过程,有效避免了计算资源在冲击波传播过程计算上的消耗。     

水下爆炸非均熵二维定常流的三族特征线解法

针对二维定常可压缩超声速非等熵柱状流,提出一种特征线差分解法,通过在沿马赫线的相容方程中添加沿流线的熵变项以描述非等熵效应,得到等熵流和非等熵流均适用的三族特征线方程组。根据水下爆炸近场特点,建立无限长柱状装药的定常模型,将三族特征线方程组用有限差分法离散求解,通过构造合适的网格保证计算格式可以数值上收敛,由此编制程序并计算几种柱状炸药的水下爆炸近场冲击波。对比有限元模拟结果和实验结果发现,特征线差分法可以比较准确地捕捉冲击波形状并计算冲击波后流场,从而验证了所提出的三族特征线差分法的准确性。     

基于小波分析的实船水下爆炸船体响应特征

为了研究水下爆炸条件下船体冲击振动响应时频特征,针对某实船非接触水下爆炸实验冲击响应测试实验数据,基于小波分析及能量统计方法对响应信号进行时频特性分析,得到了实船非接触水下爆炸冲击振动响应的时频分布和能量分布。分析结果表明,采用基于小波变换的时频分析方法,可以成功获得船体冲击响应信号不同频率段下的强度、能量和作用时间等时频细节信息,包括响应信号各频段冲击峰值、衰减过程、振动能量及其在全频率段上所占的分数。通过对小波频段能量统计以及冲击强度分析发现,冲击响应能量频段分布较广,主甲板及以下甲板全频段振动能量的80%以上在312.5 Hz以上,上层建筑甲板平台各频段冲击振动能量分数向低频段转移。     

基于区域分解的并行有限元程序及其在水下爆炸中的应用

基于MPI标准定义的消息传递接口实现了显式动力学有限元程序的并行计算.通过基于Hilbert空间填充曲线的区域分解算法,实现各区域的独立运算和区域之间共享数据的相互通信.利用程序对水下爆炸自由场中的冲击波传播规律进行了数值模拟,并通过对不同进程数下的并行加速比进行测试,验证了程序的并行效率.     

特征线法在三维弹／粘塑性波数值解析中的应用

根据广义特征理论建立了求解三维非线性应力波传播问题的特征关系式,并采用特征线法对纵向冲击载荷下弹／粘塑性矩形截面棒的三维应力波传播过程进行了数值模拟。     

基于二阶双渐近法的双层圆柱壳在水下爆炸作用下的鞭状运动

针对潜艇在水下爆炸载荷下的鞭状运动,从波动方程出发,推导了二阶双渐近法后期近似,并结合声固耦合法初步解决了双层圆柱壳的内域问题,然后将其与显式有限元耦合形成了圆柱壳结构水下爆炸流固耦合分析方法。通过简单算例验证了本文分析方法的有效性和精度。最后,基于此方法分析了双层圆柱壳结构在水下爆炸载荷下的总体响应特性以及周期比和爆距比对其影响规律。     

近自由面水下爆炸冲击载荷特性三维数值模拟

基于三维SPH 方法,对传统链表搜索算法进行了变光滑长度改进,并提出了具有较好稳定性的多相物质交界面的处理方法,模拟了三维无限域水下爆炸问题,验证了改进的三维SPH 方法模拟水下爆炸问题的可行性和有效性。在此基础上,建立了水下爆炸三维数值模型,模拟了近自由面水下爆炸过程,研究了冲击波传播特征、自由面下压力场和能量场特性以及水柱的产生过程。结果表明:自由面可将冲击波压力峰值和压力冲量最大衰减到1/3和1/7;爆深的增加会导致压力比和冲量比的等值柱面曲率变小,产生的水柱也逐渐由破碎的喷柱向高而窄的水冢过渡。     

夹芯圆柱壳在水下爆炸载荷作用下的抗冲击特性

采用声固耦合方法对夹芯圆柱壳和等质量的普通圆柱壳在爆炸载荷作用下的应变、速度和加速度进行有限元计算。结果表明:夹芯防护层对爆炸冲击波可起到较好的衰减作用,即通过芯层的塑性变形,耗散了冲击过程中产生的大部分能量,对里面的圆柱壳体起到较好的保护作用,由于夹芯防护层的存在,与等质量的普通圆柱壳相比,夹芯圆柱壳能够承受更强的爆炸冲击波,降低结构的整体变形。     

修正的虚拟介质方法用于水下爆炸冲击薄板的模拟

近期发展的修正的虚拟介质方法(MGFM)已经成功应用于求解多介质流动问题.本文研究了近体水下爆炸中强激波对薄板冲击的演变及薄板的动力响应过程.为了真实地考虑不同介质在界面处的非线性相互作用,本文将MGFM推广应用于处理流体与弹性薄板的相互作用.研究发现尽管流体和薄板结构的模拟基于各自不同的求解方法,但是用修正的虚拟介质...     

水下爆炸气泡射流载荷计算的新方法研究

由于存在强冲击、多相界面复杂运动等强非线性效应,目前对于水下爆炸气泡运动的计算方法无法给出较为可信的水射流运动特征及其载荷形式.本文基于多相可压缩流体的five-equation计算模型,引入界面函数和界面密度压缩技术,采用5阶WENO重构与HLLC近似Riemann求解器进行空间离散,时间离散采用3阶TVD Rung...     

基于小波变换的水下连续爆炸声信号特征分析

为了研究水下连续爆炸声信号的特征,利用Mallat算法,采用离散小波变换对水下连续爆炸声信号进行了分层提取分析,讨论了水下连续爆炸声信号在各频带的能量分布状况,采用Welch方法实现了对水下连续爆炸声信号的功率谱特征提取,并采用离散小波变换对声信号进行时频谱特性分析。结果表明,水下连续爆炸声信号具有很强的声功率,声压级可以达到190 dB以上,声持续时间较长,频率范围宽、声信号的能量主要集中在频率48 kHz以下,其中在低频段能量更大,这些特点使其有望成为水声干扰源。

     

水下接触爆炸下防雷舱舷侧空舱的内压载荷特性

采用模型实验方法,研究了近自由面水下接触爆炸下防雷舱舷侧空舱的内压载荷特性。根据实验模型的破坏结果和压力测试结果,分析了水下爆炸产物与防雷舱舷侧空舱的相互作用过程以及水下爆炸产物的压力变化规律。研究表明:防雷舱舷侧空舱的载荷可分为冲击波载荷、准静态压力载荷和负压载荷3种,防雷舱舷侧空舱的破坏主要由冲击波载荷和准静态压力载荷造成,并且准静态压力载荷的比冲量是冲击波载荷的数倍,而负压载荷对防雷舱舷侧空舱破坏的影响可忽略不计。