铁的冲击相变及其对层裂影响的数值模拟研究.pdfVIP

第九届全国冲击动力学学术会议论文集 铁的冲击相变及其对层裂影响的数值模拟研究 裴晓阳 贺红亮 李平 于继东 (中国工程物理研究院流体物理研究所，四川绵阳621900) 摘要：本文对铁的冲击相变及其对层裂的影响进行了数值模拟。模拟中分别采用一个相变模型和损伤度函数模型来 描述铁的冲击相变和层裂。自由面速度剖面和层裂片厚度的计算结呆和实验吻合，数值模拟再现了由于相变导致反 常的多层层裂的物理过程。结果分析表明，对于发生相变的层裂区域，层裂强度明显提高，相变动力学对层裂破坏 特性有着显著的影响。 关键词：相变，数值模拟，多层层裂 1 引言 众所周知，当冲击压力超过13卯口时，铁将发生口§占多形性结构相变(从体心立方结构到 密排六方结构)。Bancroft等人…最早从波剖面测量发现了该现象，很快被静高压实验所证实埋1。随 后，Barker等人及众多学者睇刊对铁的冲击相变作了一系列波剖面测量，得剑了精确的铁的雨贡纽数 据。Andrews【78l计算了热力学平衡态情况下铁的雨贡纽状态，Hayes旧。101则提出了一个非平衡态模型 力学、热力学及动力学过程进行了精彩评述。最近，文献u}131报道了通过x光衍射对单晶铁冲击相 5‘。 变的首次直接测量，实验结果证实了之前的分子动力学模拟¨41 在高于相变压力时，铁有可能发生光滑层裂u…，普遍认为这个现象和冲击相变后逆相变造成的 稀疏冲击波相关，稀疏冲击波使密度、声速等物理量发生突变。许多学者…屯刮通过爆轰驱动加载、 飞片撞击及激光加载等加载技术对相变对层裂的影响进行了广泛的实验研究。 尽管如此，仍然还存在很多疑问。本文主要关注铁在发生相变情况下的层裂特性，以及相变动 力学对层裂的影响。针对铁的相变，建立了一个相变模型，包括混合相模型、热力学势模型以及相 变动力学模犁。对于层裂特性则采用一个损伤度函数模型进行描述。利用所建立的模型对铁的相变 压力之上的层裂实验进行了数值模拟，并分析了由于相变对层裂过程的影响。 2物理模型 本文采用一个二维动力有限元程宇——-DEFEL对相关实验进行了数值模拟，物理模型包括一个 相变模型和一个损伤度函数模型来描述铁的动态响应特性。 A．相变模型 首先把每一个材料微元都看成是口相和s相的混合相，采用两相状态方程性卜驯的思想来来描述 铁的相变，比容和比内能满足如下关系： y=(1一石)K+x％ (1) E=(1一x)巨+x岛 下表l和2 式中：x代表s相的质量分数(O≤x≤1)，y和E分别代表混合相的比容和比内能， 分别代表口相和占相。两相之间满足压力平衡和热平衡： 尸(y，E)=墨(K，巨)=￡(％，岛) (3) 丁(y，E)=五(K，巨)=五(砭，易) (4) 425 对于每个单相采用施P—G，|锄P施，z状态方程进行描述，以雨贡纽状态为参考线，状态方程形 式如下： 尸一名=芳(E—EH) (5) 2风爵南 式中：尸为压力，／为格林奈森系数，风未初始密度，c0和兄表述线性的冲击波速度和粒子速度 关系参数，77=l—y／K代表压缩比。图1给出了本文计算的铁的雨贡纽线和文献的对比结果，计 算参数列于表l。 2∞