简要介绍了物质点法(material point method)的离散原理,通过引入Johnson-Cook材料模型和Mie-Gruneisen状态方程,将其用于超高速碰撞问题的分析中,并编制了相应三维物质点法程序MPM3D。该方法避免了拉格朗日格式因网格畸变产生的数值困...简要介绍了物质点法(material point method)的离散原理,通过引入Johnson-Cook材料模型和Mie-Gruneisen状态方程,将其用于超高速碰撞问题的分析中,并编制了相应三维物质点法程序MPM3D。该方法避免了拉格朗日格式因网格畸变产生的数值困难,也克服了欧拉格式材料界面跟踪问题以及因非线性对流扩散项而引起的数值困难。利用该程序对Taylor杆高速碰撞问题和空间碎片防护超高速碰撞问题进行了数值模拟,所得数值结果与实验结果基本吻合,验证了程序的正确性,说明了物质点法在分析超高速碰撞问题时相对于有限元法的优势。展开更多
基于物质点方法(material point method,MPM)理论框架,提出了处理饱和多孔介质与固体间动力接触问题的新方法。其中饱和多孔介质的动力学响应通过文献[1]中发展的耦合物质点方法进行分析,单相固体的力学行为由传统单相物质点方法进行预...基于物质点方法(material point method,MPM)理论框架,提出了处理饱和多孔介质与固体间动力接触问题的新方法。其中饱和多孔介质的动力学响应通过文献[1]中发展的耦合物质点方法进行分析,单相固体的力学行为由传统单相物质点方法进行预测。通过本文提出的接触算法使二者相结合,在保证饱和多孔介质与固体间不存在相互穿透的前提下,允许饱和多孔介质与固体间的相互滑动,以预测整个接触/碰撞系统的动力学响应。同时进行了数值算例计算,通过算例验证了此方法的正确性,展示了此方法有效性。展开更多
提出用于饱和多孔介质动力学响应分析的耦合物质点方法(Coupling material point method)。采用u–p形式控制方程对饱和多孔介质进行数值模拟,建立了耦合物质点方法的弱形式离散求解方程,阐述了耦合物质点方法压强场边界条件的处理方式...提出用于饱和多孔介质动力学响应分析的耦合物质点方法(Coupling material point method)。采用u–p形式控制方程对饱和多孔介质进行数值模拟,建立了耦合物质点方法的弱形式离散求解方程,阐述了耦合物质点方法压强场边界条件的处理方式,通过引入边界压强层近似描述指定压强边界,并给出了算法的实施过程。通过数值算例,验证了所提出的耦合物质点方法用于饱和介质动力学分析的正确与有效性。展开更多
物质点法MPM(Material Point Method)是无网格方法之一。它是在质点网格法(PIC)基础上发展而来的一种新数值方法,它利用了欧拉法和拉格朗日法两者的优点,计算物质点在冲击载荷下的应力和应变;通过物质点来跟踪材料体的变形和破损,而在...物质点法MPM(Material Point Method)是无网格方法之一。它是在质点网格法(PIC)基础上发展而来的一种新数值方法,它利用了欧拉法和拉格朗日法两者的优点,计算物质点在冲击载荷下的应力和应变;通过物质点来跟踪材料体的变形和破损,而在整个计算过程中背景网格始终固定不变,避免了重新划分网格。本文应用MPM法计算三维爆炸焊接问题,在爆轰载荷作用下的飞板和基板的金属动态变形过程进行了三维数值模拟,并且对飞板的碰撞点速度和爆轰压力变化进行了计算分析。展开更多
利用MPM(material point method)计算方法,以爆炸焊接作为实例,对炸药滑移爆轰过程和金属飞板与基板之间的碰撞变形进行了数值模拟.并将无网格MPM法的数值计算结果与近似解析公式的计算结果进行了比较,二者基本吻合,有力地证明了无网格...利用MPM(material point method)计算方法,以爆炸焊接作为实例,对炸药滑移爆轰过程和金属飞板与基板之间的碰撞变形进行了数值模拟.并将无网格MPM法的数值计算结果与近似解析公式的计算结果进行了比较,二者基本吻合,有力地证明了无网格MPM法在求解爆炸冲击问题中的有效性和健壮性.展开更多
A new multi-mesh contact algorithm for three-dimensional material point method is presented. The contact algorithm faithfully recovers the opposite acting forces between colliding bodies. Collision procedures between ...A new multi-mesh contact algorithm for three-dimensional material point method is presented. The contact algorithm faithfully recovers the opposite acting forces between colliding bodies. Collision procedures between regular bodies and/or rigid bodies are treated within the same framework. Multi-value of momentum and mass are defined on every node to describe the contact/sliding/separation procedure. Both normal and tangential velocities of each particle at the contact surface are calculated in respective individual mesh. A Coulomb friction is applied to describe the sliding or slipping between the contacting bodies. The efficiency of the contact algorithm is linearly related to the number of the contacting bodies because the overlapped nodes are labeled by sweeping the material particles of all bodies when the nodal momentum and mass are formed at every time step. Numerical simulation shows that our contact algorithm possesses high accuracy and low numerical energy dissipation, which is very important for solving collision problems.展开更多
The physical fields in porous materials under strong shock wave reaction are very complicated. We simulate such systems using the grain contact material point method. The complex temperature fields in the material are...The physical fields in porous materials under strong shock wave reaction are very complicated. We simulate such systems using the grain contact material point method. The complex temperature fields in the material are treated with the morphological characterization. To compare the structures and evolution of characteristic regimes under various temperature thresholds, we introduce two concepts, structure similarity and process similarity. It is found that the temperature pattern dynamics may show high similarity under various conditions. Within the same material, the structures and evolution of high-temperature regimes may show high similarity if the shock strength and temperature threshold are chosen appropriately. For process similarity in materials with high porosity, the required temperature threshold increases parabolically with the impact velocity. When the porosity becomes lower, the increasing rate becomes higher. For process similarity in different materials, the required temperature threshold and the porosity follow a power-law relationship in some range.展开更多
为了解决当炸弹在近场爆炸时爆轰波驱动破碎的弹片共同作用于混凝土墙壁的过程中所涉及的多物理场计算和多相介质耦合分析等问题,利用物质点法(material point method,MPM)不需要考虑物质间的分界面、耦合条件自动满足等特点,应用无网格...为了解决当炸弹在近场爆炸时爆轰波驱动破碎的弹片共同作用于混凝土墙壁的过程中所涉及的多物理场计算和多相介质耦合分析等问题,利用物质点法(material point method,MPM)不需要考虑物质间的分界面、耦合条件自动满足等特点,应用无网格MPM法对两种类型的炸弹(带金属外壳和不带金属壳)产生的爆炸场、弹片破碎和混凝土墙壁的破坏进行数值模拟。数值结果表明,无网格MPM法是计算多相介质爆炸效应的一种有效的算法。展开更多
文摘简要介绍了物质点法(material point method)的离散原理,通过引入Johnson-Cook材料模型和Mie-Gruneisen状态方程,将其用于超高速碰撞问题的分析中,并编制了相应三维物质点法程序MPM3D。该方法避免了拉格朗日格式因网格畸变产生的数值困难,也克服了欧拉格式材料界面跟踪问题以及因非线性对流扩散项而引起的数值困难。利用该程序对Taylor杆高速碰撞问题和空间碎片防护超高速碰撞问题进行了数值模拟,所得数值结果与实验结果基本吻合,验证了程序的正确性,说明了物质点法在分析超高速碰撞问题时相对于有限元法的优势。
文摘基于物质点方法(material point method,MPM)理论框架,提出了处理饱和多孔介质与固体间动力接触问题的新方法。其中饱和多孔介质的动力学响应通过文献[1]中发展的耦合物质点方法进行分析,单相固体的力学行为由传统单相物质点方法进行预测。通过本文提出的接触算法使二者相结合,在保证饱和多孔介质与固体间不存在相互穿透的前提下,允许饱和多孔介质与固体间的相互滑动,以预测整个接触/碰撞系统的动力学响应。同时进行了数值算例计算,通过算例验证了此方法的正确性,展示了此方法有效性。
文摘提出用于饱和多孔介质动力学响应分析的耦合物质点方法(Coupling material point method)。采用u–p形式控制方程对饱和多孔介质进行数值模拟,建立了耦合物质点方法的弱形式离散求解方程,阐述了耦合物质点方法压强场边界条件的处理方式,通过引入边界压强层近似描述指定压强边界,并给出了算法的实施过程。通过数值算例,验证了所提出的耦合物质点方法用于饱和介质动力学分析的正确与有效性。
文摘物质点法MPM(Material Point Method)是无网格方法之一。它是在质点网格法(PIC)基础上发展而来的一种新数值方法,它利用了欧拉法和拉格朗日法两者的优点,计算物质点在冲击载荷下的应力和应变;通过物质点来跟踪材料体的变形和破损,而在整个计算过程中背景网格始终固定不变,避免了重新划分网格。本文应用MPM法计算三维爆炸焊接问题,在爆轰载荷作用下的飞板和基板的金属动态变形过程进行了三维数值模拟,并且对飞板的碰撞点速度和爆轰压力变化进行了计算分析。
基金The project supported by the Science Foundation of Laboratory of Computational Physics,Science Foundation of China Academy of Engineering Physics,and National Natural Science Foundation of China under Grant Nos.10702010,10775018,10472052,and 10604010
文摘A new multi-mesh contact algorithm for three-dimensional material point method is presented. The contact algorithm faithfully recovers the opposite acting forces between colliding bodies. Collision procedures between regular bodies and/or rigid bodies are treated within the same framework. Multi-value of momentum and mass are defined on every node to describe the contact/sliding/separation procedure. Both normal and tangential velocities of each particle at the contact surface are calculated in respective individual mesh. A Coulomb friction is applied to describe the sliding or slipping between the contacting bodies. The efficiency of the contact algorithm is linearly related to the number of the contacting bodies because the overlapped nodes are labeled by sweeping the material particles of all bodies when the nodal momentum and mass are formed at every time step. Numerical simulation shows that our contact algorithm possesses high accuracy and low numerical energy dissipation, which is very important for solving collision problems.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10702010, 10775018, and 10771019)Science Foundation of Laboratory of Computational Physics and Science Foundation of China Academy of Engineering Physics (Grant Nos. 2009A0102005 and 2009B0101012)
文摘The physical fields in porous materials under strong shock wave reaction are very complicated. We simulate such systems using the grain contact material point method. The complex temperature fields in the material are treated with the morphological characterization. To compare the structures and evolution of characteristic regimes under various temperature thresholds, we introduce two concepts, structure similarity and process similarity. It is found that the temperature pattern dynamics may show high similarity under various conditions. Within the same material, the structures and evolution of high-temperature regimes may show high similarity if the shock strength and temperature threshold are chosen appropriately. For process similarity in materials with high porosity, the required temperature threshold increases parabolically with the impact velocity. When the porosity becomes lower, the increasing rate becomes higher. For process similarity in different materials, the required temperature threshold and the porosity follow a power-law relationship in some range.
文摘为了解决当炸弹在近场爆炸时爆轰波驱动破碎的弹片共同作用于混凝土墙壁的过程中所涉及的多物理场计算和多相介质耦合分析等问题,利用物质点法(material point method,MPM)不需要考虑物质间的分界面、耦合条件自动满足等特点,应用无网格MPM法对两种类型的炸弹(带金属外壳和不带金属壳)产生的爆炸场、弹片破碎和混凝土墙壁的破坏进行数值模拟。数值结果表明,无网格MPM法是计算多相介质爆炸效应的一种有效的算法。