冲击作用下空洞塌缩问题的物质点模拟.pdfVIP

第八扁伞困爆炸力学学术会议论集 冲击作用下空洞塌缩问题的物质点模拟 许爱国 潘小飞 张广4财 朱建士 北京应用物理’j计算数学研究所计算物理实验空高于 100088 摘要：空涮的冲击塌缩是材料老化、炸药起爆和声敛发光等研究领域所共蚓关注的问题。奉文使用l捌体力学领域新 近发展起来的物质点方法进行模拟研究。重点考察空涮塌缩的对称性．。j冲击强度、其‘I三界面的关系和纯流体力学模型所 忽视的一些动力擎、热力学特征。在冲击速度较高的情况下，研究J，空洞内射流产生以及材料表面附近空洞内的射流冲 出材料表面的过程’j规律；在冲击速度较低的情况下，发现空洞在一定条件下不能完伞塌缩，月．叫。能会发生近似均匀、 对称塌缩的情形，研究__r对称性转变的规律。发现，随着塌缩方式的改变，“热点”在材料内的分布规律发生相应的改 变。最后使用物质点方法研究了HMX材料中窄涮1竽租i对冲击起爆的影响。窄洞的塌缩使“热点”的温度火幅度提升。由 1i变。 关键词：窄洞塌缩物质点法材料老化炸药起爆 1．引言 空洞塌缩问题的研究可以追溯到十九世纪后期。当时人们发现汽船推进器中的空洞具有很强的侵 蚀作刚。随后，不同应刚背景下的空洞塌缩问题相继被发现和提出。从而，空洞塌缩问题的研究便引 起人们的普遍关注。随着研究的进行、认识的深入，空洞塌缩带来的新现象也不断涌现。现在，空洞 的冲击塌缩问题已经成为水力学、材料老化、炸药起爆、惯性约束聚变、卢致发光和化jJ：等领域所共 同关注的课题【1．3】。 空洞塌缩问题的多数理论研究为流体力学描述。在这种描述中，边界条件的描述是一个核心问题。 一种简单描述为将注入气泡的“液体舌头”视为边界的儿何流效应：另一种描述为冲击作川下来流方向 泡鼙的破碎形成射流。需要指出的是，在研究(金属)材料中空洞的冲击塌缩时，流体力学描述忽视了 塑性应变、强化效应以及与变形历史相关的一些动力学特征。由丁．塌缩过程涉及到很强的非线性，其理 论描述非常复杂，所以除实验外，这方面的研究在很人程度上依赖丁．数值模拟。 冲击波的主要效应为使得受载荷作用的自由表面进行加速，朝对面静，I}_=的腔蹙运动。由于几何效应， 粒子速度发生汇聚，更高的粒子速度会山现。另外，冲击波在传播过程中，可能会压塌相邻的空洞；如 果作用在I矧体表面，‘损伤和进一步的成腔行为可能会产生。对于这类涉及到材料大变形的问题，纯拉氏 和纯欧氏算法均遇到了一些不易克服的困难：拉氏方法能够很好的跟踪物质的界面，但在处理人变形时 会遇到网格畸变的问题；在欧氏方法中，网格不随物体的变形而变化，但是它不易跟踪物质的界面。为 充分发挥欧氏和拉氏算法的优点、避免其缺点，本文将使用最近发展起来的一种混合方法一物质点法来 模拟材料中空洞问题的冲击塌缩过程。 2．物质点方法 ^， 在该方法由流体力学的PIC发展而来。其中，连续物体被离散成”P个离散的拉格朗日质点，每个 质点携带有质繁聊，、鹿力op、速度V、密度砟等信息：背景网格采用欧拉法描述，网格结点f携带有 速度Vi利加速度a，的信息。在每一个时间步，将物质点P携带的物质信息映射到背景网格结点f，在背 景网格’1，点z上求解动量方程，获得网格’仃点的速度、加速度；然后再映射同物质点处，更新下一时刻 冲击作用下夺洞塌缩问题的物质点模拟 物质点所携带的物质信息。在每一个时间步均采川原始背景网格，冈此避免了冈网格畸变而产生的数值 困难。_卜．面我们简要地叙述其基本原理。 在物质点法中，每个物质点的质量不随时间变化，所以质量连续性方程。 竺+厕．v：o (1) 。 dt 自动满足，其中∥代表密度，V代表速度。拉氏力学平衡方程为： ．(2)