弹性壳结构静力与动力分析的光滑粒子法

本文通过采用移动最小二乘函数作为近似函数 和完全拉格朗日方程作为近似方程来改善光滑粒子法的稳定性和数值精度; 在此基础上, 提出了壳结构静力分析的光滑粒子法, 并完善了壳结构动力分析方法; 最后, 采用国际公认的壳结构的标准测试模型对静力和动力问题分别进行了验证, 所得结果与已有数据吻合良好, 证明了本文数值模型的有效性和可靠性, 为光滑粒子法进一步在裂纹、破碎等非线性壳结构中的应用提供参考. 关键词： 弹性壳 静力与动力分析 光滑粒子法 完备性和稳定性     

舰船水下接触爆炸数值研究

计及船体外板的材料强度,基于SPH算法建立水下接触爆炸数值模型,提出混合密度、广义光滑长度等解决数值模拟中计算崩溃的问题,对舰船的水下接触爆炸进行数值模拟,计算结果和冲击波基本理论及接触爆炸试验吻合良好.结果表明,接触爆炸时钢结构会产生圆形冲塞破口,其破口大小和炸药的厚度接近,当板与炸药的正则化厚度达1.8时,板仅会发生较大的塑性变形而不会冲塞破坏.并得到钢的反射系数约为1.2,及接触爆炸钢板的近似变形公式.     

自由上浮气泡运动特性的光滑粒子流体动力学模拟

基于虚功原理, 在Hu X Y等和Grenier N等的研究结果基础上推导了多相流光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)控制方程, 采用精度较高的黏性力和表面张力模型, 发展了一套适用于具有大密度比和大黏性比界面的多相流SPH方法. 首先, 通过施加人工位移修正, 适当背景压力和异相界面力, 使得计算全程粒子分布相对均匀, 改善了界面处的失稳现象, 防止了异相界面处粒子的非物理性穿透; 在此基础上, 利用方形流体团振荡模型对表面张力模型进行了验证, 数值结果与解析解甚为吻合; 然后采用上浮气泡经典数值算例对比研究了不同黏性力计算方法、不同核函数的适用性以及人工位移修正的效果; 最后, 对单个气泡的上浮、变形、撕裂以及垂向两个气泡的追赶、融合等现象进行了模拟, 初步揭示了气泡上浮过程中各种有趣物理现象的细节过程和动力学机理.     

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

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

Investigation on charge parameters of underwater contact explosion based on axisymmetric SPH method

This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.