滑移爆轰问题无网格MPM法数值模拟

利用MPM(material point method)计算方法,以爆炸焊接作为实例,对炸药滑移爆轰过程和金属飞板与基板之间的碰撞变形进行了数值模拟.并将无网格MPM法的数值计算结果与近似解析公式的计算结果进行了比较,二者基本吻合,有力地证明了无网格MPM法在求解爆炸冲击问题中的有效性和健壮性.     

层裂问题的交错网格物质点法研究

交错网格物质点法(Staggered Grid Material Point Method, SGMP)有效地消除了物质点法(Material Point Method, MPM)的跨网格数值误差,是模拟冲击爆炸等涉及大变形与材料破坏问题的一种有效数值分析方法.为了研究层裂问题,论文将质点失效模型引入交错网格物质点法中,对Hopkinson杆碰撞、飞片撞击靶板碰撞和混凝土受爆炸载荷冲击等层裂问题进行了数值模拟研究,结果表明SGMP可以很好地模拟层裂问题.     

无网格MPM法三维爆炸焊接数值模拟

物质点法MPM(Material Point Method)是无网格方法之一.它是在质点网格法(PIC)基础上发展而来的一种新数值方法,它利用了欧拉法和拉格朗日法两者的优点,计算物质点在冲击载荷下的应力和应交;通过物质点来跟踪材料体的变形和破损,而在整个计算过程中背景网格始终固定不变,避免了重新划分网格.本文应用MPM法计算三维爆炸焊接问题,在爆轰载荷作用下的飞板和基板的金属动态变形过程进行了三维数值模拟,并且对飞板的碰撞点速度和爆轰压力变化进行了计算分析.     

土中爆炸冲击混凝土道面的数值模拟

炸药在多相介质中的爆炸是气态爆轰产物和固态介质之间的耦合作用问题，对其进行数值模拟一般不考虑炸药而直接将爆轰压力加在介质上，而耦合方法对炸药和其周围的介质进行离散，炸药和介质采用统一的控制方程，通过界面条件实现相互作用。任意拉格朗日欧拉（Arbitary Lagrangian and Eularian）法，是最近兴起的适用于计算固体大变形问题的数值计算方法。将其引入爆炸力学计算，可以避免由于物质的大变形引起的单元畸变和计算困难。同时混凝土的损伤模型，更加符合混凝土在高应变率加载条件下的力学特性和破坏机理。将这三种方法引入数值模拟，计算结果与相关试验结果有较好的吻合。     

结构爆炸毁伤的浸没多介质有限体积物质点法

结构在爆炸载荷作用下的毁伤现象涉及强非线性激波、固体结构极端变形和破坏破碎、强流固耦合, 给数值计算方法带来了极大的困难与挑战. 针对结构爆炸毁伤问题, 建立了浸没多介质有限体积物质点法(iMMFV-MPM), 采用基于黎曼求解器的多介质有限体积法(MMFVM)模拟爆炸产物和空气的多介质流体, 采用物质点法(MPM)模拟固体结构, 并将提出的基于拉格朗日乘子的连续力浸没边界法(lg-CFIBM)扩展到多介质流体中以处理流固耦合边界条件. 该算法在每个时间步严格满足流固耦合界面处的速度边界条件及动量守恒方程, 不需要重构流固耦合界面, 能够有效地模拟近场爆炸下爆炸产物与结构的相互作用、激波与结构的相互作用和演化以及结构的动态断裂和拓扑变化. 利用iMMFV-MPM对近场爆炸下方形钢筋混凝土靶板的失效模式、外爆载荷下建筑物的毁伤现象以及多腔室内爆炸试验进行了模拟, 模拟结果与相关实验数据吻合良好, 验证了所建立的流固耦合算法的有效性及精度.      

极端变形问题的物质点法研究进展

物质点法采用了拉格朗日和欧拉双重描述,将物体离散为一组在固定于空间的网格(欧拉描述)中运动的质点(拉格朗日描述),有效地综合了拉格朗日法和欧拉法的优点,是分析冲击爆炸等极端变形问题的一种有效方法. 本文系统总结了本课题组近年来针对冲击、爆炸和流固耦合等极端变形问题在物质点法算法方面的研究进展,并简要介绍了开发的三维显式并行物质点法数值仿真软件MPM3D及其在超高速碰撞、侵彻、爆炸、边坡失效、金属切削和流固耦合等问题中的应用.     

基于场分离的水下爆炸流固耦合计算方法

基于谱单元方法对水下爆炸流固耦合机理进行了研究,针对总场模型中入射波在流体网格中传播会产生失真现象的问题,采用场分离技术,通过理论推导建立了水下爆炸流固耦合问题的散射场数值模型.经过编程计算比较得出,散射场数值模型有效解决了入射波在流体网格中传播的失真问题,比总场模型计算精度更高,可较好地应用于水下爆炸流固耦合问题的求...     

物质点法的理论和应用

物质点法采用质点离散材料区域, 用背景网格计算空间导数和求解动量方程,避免了网格畸变和对流项处理, 兼具拉格朗日和欧拉算法的优势, 非常适合模拟涉及材料特大变形和断裂破碎的问题. 本文详细论述了物质点法在基本理论、算法和软件开发方面的进展, 包括广义插值物质点法、接触算法、自适应算法、并行算法、与其他算法的杂交和耦合等. 系统地总结了物质点法在超高速碰撞、冲击侵彻、爆炸、动态断裂、流固耦合、多尺度分析、颗粒材料流动和岩土失效等一系列涉及材料特大变形问题中的应用,展示了其相对于传统数值计算方法的优势.     

混凝土中爆炸模拟的数值方法比较

应用AUTODYN商业软件对装药在混凝土靶介质中的爆炸过程进行数值模拟,对比了混凝土中爆炸模拟的不同方法。将各数值方法模拟得到的混凝土毁伤特征尺寸与相关实验数据进行比较,研究了各数值方法的准确性及描述混凝土毁伤特征的可靠性参数。结果发现：Euler算法计算的压力最大,SPH和Lagrange耦合算法计算的压力最小;Euler算法和SPH算法计算的压力衰减过程会发生比较显著的扰动。以不同的平均网格尺寸对数值模拟结果的影响分析了各种数值方法的稳定性。比较靶体上同一测点的压力时间历程表明网格变化对测点处压力影响不大。通过应用各种数值方法模拟再现混凝土中的爆炸过程,揭示了各种方法的优缺点,明确了现有商业软件模拟混凝土中爆炸的适用性。     

大体积混凝土水下接触爆炸破坏分区特征分析

为探究大体积混凝土水下接触爆炸破坏分区特征，基于水中爆炸冲击波与混凝土的相互作用过程分析，建立了综合考虑爆炸冲击波冲击破碎和爆轰产物准静态压拉致裂机制的混凝土爆炸破坏分区计算方法，并与有限元数值模拟和试验实测数据开展对比验证。结果表明：与空中接触爆炸相比，水对爆轰产物膨胀起抑制作用，使得爆炸荷载持时增加、作用于周围介质的冲量增大；采用建议的环向压碎判据计算破碎区，并将开裂区分为动态压裂、准静态压裂和准静态拉裂区的计算方法能够很好地预测混凝土水下接触爆炸破坏分区范围；炸药类型和起爆水深一定时，混凝土的抗拉强度和抗压强度比对开裂区范围起重要影响。     

地下坑道对其临界震塌爆距处钻地武器爆炸载荷的动力响应

利用非线性显式动力有限元程序,采用多物质流固耦合计算方法,就GBU-28钻地弹在地下坑道临界震塌爆距处爆炸时,对地下直墙拱坑道的动力响应进行数值模拟。根据围岩动力稳定性和混凝土动态强度判据,结合模拟结果,分析衬砌结构与围岩的相互作用。钻地弹在直墙圆拱断面的坑道临界震塌爆距处爆炸时:围岩处于临界破坏状态,但混凝土衬砌结构处于稳定状态;拱顶的应力峰值明显,且柱状装药情况下,爆炸近区的应力较集团装药情况下的大;拱肩位置出现应力集中;围岩与衬砌结构特征位置处的相互作用载荷与对应质点的振动速度相互耦合,基本成对应的关系。     

回收战斗部破片的新型爆炸容器及应用

介绍一种回收战斗部破片的新型爆炸容器装置。该装置采用水介质和可拆卸内衬防护板结构抵御破片侵彻，能够反复使用且有效回收爆炸试验后的破片，较好地反映战斗部的破碎状况。利用该装置进行了典型弹丸爆炸试验，研究了战斗部破片特性，证实了爆炸容器在战斗部破碎性试验中的应用是可行的。     

无网格MPM法三维前处理系统设计

无网格法在解决爆炸冲击、裂纹扩展、高速碰撞以及大变形等问题时,由于与网格无关,能够避免有限元法因网格畸变计算精度降低或计算终止的矛盾,是目前求解固体力学问题的有效方法之一,也是近年来研究的热点。前处理对于无网格法数值计算是非常重要的环节,基于面向对象程序设计思想,应用C 和OpenGL开发了三维图形交互方式的MPM法前处理系统,提出了应用递归四叉树法、扫描表示法和CSG法划分任意形状非等间距的MPM法质点单元,为建立通用的无网格前处理平台提出了有效的解决方案。     

Coupling between finite element method and material point method for problems with extreme deformation

As a Lagrangian meshless method, the material point method (MPM) is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method (FEM) for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with experimental data and the efficiency of this method is higher than that of both FEM and MPM.     

混凝土中爆炸应力波衰减规律的数值模拟研究

基于Kong-Fang混凝土材料模型和LS-DYNA的多物质ALE算法,开展混凝土中爆炸波衰减规律的数值模拟研究。首先,基于已有实验数据对材料模型参数和数值算法的可靠性进行了验证,在此基础上分析球形装药在混凝土自由场中爆炸波衰减规律,利用量纲分析和数值模拟拟合了球形装药在混凝土自由场中近区爆炸波峰值应力计算公式并明确其适用范围;然后,分析装药埋深对混凝土中装药正下方不同距离处爆炸波峰值应力分布的影响,建立了耦合系数与装药埋深和测点距离之间的定量关系。结果表明:Kong-Fang混凝土材料模型可实现对混凝土中爆炸波传播衰减规律的高精度数值模拟;定义混凝土中装药质量系数和耦合常数,可定量描述装药埋深和测点距离对峰值应力耦合系数的影响;建立的混凝土中近区爆炸波峰值应力计算公式可较准确地快速预测不同装药埋深、不同测点距离和不同混凝土强度时爆炸波峰值应力。研究结果可为混凝土结构抗爆设计和爆炸毁伤评估提供参考。