多相介质爆炸冲击响应物质点法数值模拟

为了解决当炸弹在近场爆炸时爆轰波驱动破碎的弹片共同作用于混凝土墙壁的过程中所涉及的多物理场计算和多相介质耦合分析等问题,利用物质点法(material point method, MPM)不需要考虑物质间的分界面、耦合条件自动满足等特点,应用无网格MPM法对两种类型的炸弹（带金属外壳和不带金属壳）产生的爆炸场、弹片破碎和混凝土墙壁的破坏进行数值模拟。数值结果表明,无网格MPM法是计算多相介质爆炸效应的一种有效的算法。     

超高速碰撞问题的三维物质点法

简要介绍了物质点法（material point method）的离散原理,通过引入Johnson-Cook材料模型和Mie-Gruneisen状态方程,将其用于超高速碰撞问题的分析中,并编制了相应三维物质点法程序MPM3D。该方法避免了拉格朗日格式因网格畸变产生的数值困难,也克服了欧拉格式材料界面跟踪问题以及因非线性对流扩散项而引起的数值困难。利用该程序对Taylor杆高速碰撞问题和空间碎片防护超高速碰撞问题进行了数值模拟,所得数值结果与实验结果基本吻合,验证了程序的正确性,说明了物质点法在分析超高速碰撞问题时相对于有限元法的优势。     

爆炸焊接复合界面波与温度场物质点法分析

爆炸焊接是一门双金属复合工程技术,在炸药爆轰载荷驱动下,飞板高速冲击基板时,两金属板复合界面处在高温高压作用下材料发生塑性流动并形成周期性波状界面,波状界面的形成与复合界面处的材料熔化和变形直接相关。本文应用物质点法对爆炸焊接界面波的形成和界面温度场进行数值模拟,同时开展双金属爆炸焊接实验,并结合物质点法的三维数值模拟对爆炸焊接界面波的形态和界面材料高温软化进行分析。     

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

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

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

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

非均质固体炸药冲击起爆的物质点法

将点火增长方程引入到物质点法中,编制了相应的计算程序MPM-SDT,并求解了冲击起爆问题.通过与传统算法和实验结果的比较,验证了物质点法应用于冲击起爆问题的可行性.物质点法避免了Lagrange法的网格畸变缠结和Euler法的对流误差.最后,对钨、铜、钢等材料破片撞击屏蔽炸药问题进行了多组数值模拟.3种材料对屏蔽炸药的...     

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

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

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

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

爆炸焊接界面波物质点法三维数值模拟

基于冲击动力学和爆炸焊接理论,采用物质点法对爆炸焊接界面波的形成进行三维数值模拟。通过数值模拟结果与爆炸焊接实验结果的对比,对复合界面材料的塑性流动变形以及界面波形成的机理进行探讨。结果表明:界面波是因为在碰撞点处的金属材料发生熔化并产生涡旋流动形成的;同时也说明采用物质点法模拟爆炸焊接界面波的形成是可行的。     

无网格法的理论及应用

详细论述了近年来迅速发展的无网格法的理论基础及其在各个领域内的应 用. 无网格法网格依赖性弱, 避免了传统的有限元、边界元等基于网格的数值方法 中可能出现的网格畸变和扭曲, 在一些有限元、边界元等方法难以较好处理的领域体现 出独特的优势. 以加权余量法为主线归纳了已有的30多种无网格法, 各类 无网格法的主要区别在于使用了不同的加权余量法和近似函数. 详尽介绍 了各种无网格近似方案(包括移动最小二乘近似、核近似和重构核近似、单位分 解近似、径向基函数近似、点插值近似、自然邻接点插值近似等)和无网格法 中常用的各类加权余量法(伽辽金格式、配点格式、局部弱形式、加权最小二乘 格式和边界积分格式等), 并讨论了数值积分方法和边界条件的处理等问题. 在 此基础上较系统地总结了无网格法在冲击爆炸、裂纹传播、超大变形、结 构优化、流固耦合、生物力学和微纳米力学等领域的应用, 展示了无网格法相 对于传统数值方法的优势.      

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

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

流固耦合不可压物质点法及其在晃动问题中的应用

作为一种混合拉格朗日欧拉法,物质点法在流固耦合问题中具有重要的应用前景。对于自由液面的流动问题,基于物质点法框架已建立了弱可压物质点法和完全不可压物质点法,但在处理流固耦合问题时遇到了困难。弱可压物质点法由于采用可压缩状态方程,导致求解时间步长过小,压力振荡严重,产生了非物理的飞溅现象;完全不可压物质点法基于投影算法和不可压条件,消除了弱可压物质点法的压力振荡,提高了时间步长,但难以处理移动边界问题。基于变分形式的投影算法提出了一种新型流固耦合不可压物质点法,得到了体积加权的压力泊松方程PPE(Pressure Poisson Equation),解决了完全不可压物质点法无法处理不规则边界和移动边界的问题。采用流固耦合不可压物质点法研究了运动刚体容器中的液体晃动问题,并与已有实验和数值结果进行对比,验证了算法的正确性和精度。     

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

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

A multiscale material point method for impact simulation

To better simulate multi-phase interactions involving failure evolution, the material point method (MPM) has evolved for almost twenty years. Recently, a particle-based multiscale simulation procedure is being developed, within the framework of the MPM, to describe the detonation process of energetic nano-composites from molecular to continuum level so that a multiscale equation of state could be formulated. In this letter, a multiscale MPM is proposed via both hierarchical and concurrent schemes to simulate the impact response between two microrods with different nanostructures. Preliminary results are presented to illustrate that a transition region is not required between different spatial scales with the proposed approach.     

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.     

无网格法研究进展及其应用

从加权残量法的角度出发，系统地总结了现有各种无网格法的基本格式，阐明了无网格法的特点，论述了无网格法的研究进展，给出了无网格法在碰撞、动态裂纹扩展、金属加工成型、流体力学以及其它领域中的应用。