原子—有限元方法的研究铁中裂纹低温脆性解理扩展.pdfVIP

原子一有限元方法研究铁中 裂纹低温脆性解理扩展 郭雅芳 王崇愚 (钢铁研究总院) 摘要运用原子一有限元方法(FEAt)研究了裂纹在低温加载条件下的脆性解理扩 0 1l边界I型裂纹扩展 展。原子模拟结果表明，低温和一定外载下bcc．Fe中f1O)(o 为裂尖原子键断裂与层错或李晶扩展相伴随的过程，裂纹低温脆性解理扩展的最大速 度可达到1125m／s，即约为06Va。相应的有限元分析表明，裂尖存在较大的能量和 应力集中。这是导致裂尖原子键断裂及弹性孪晶形成的原因。 关键词 原子一有限元裂纹扩展孪晶形成 一定应力和载荷条件下金属材料的断 力学模拟中也观测到层错及挛晶的生成 裂过程贯穿宏观、介观、微观尺度及原子、 [8-101。但是，人们对层错和孪晶的形成和运 电子层次，从裂尖原子键的断裂到宏观材 动与脆性解理过程的关系还缺乏足够认 料的破损，空间尺度跨越约107量级，是典 识。 型的多尺度现象。所以，为了深入探讨材 本文采用原子模拟技术与有限元方法 料断裂、破损的本质，进一步提高材料的 相结合的原子．有限元方法，研究了bcc-Fe 0 1 l边界I型裂纹的动态扩展， 强度，人们一直在探索运用宏微观相结合 中f1O}O 的方法研究材料的断裂问题，这其中包括 同时给出相应的能量和应力分析。 采用原予模拟技术与传统的有限元方法相 结合的原子．有限元方法”月。 1 研究方法及模型 对于体心立方金属如铁，其断裂强度 强烈依赖于环境的温度条件，这就是我们 所常说的脆．韧转变(BDT)。由实验早已得 到bcc．Fe在低碳多晶情况下脆．韧转变温度 为250Kf”，即低于脆．韧转变温度时材料为 脆性断裂，高于脆．韧转变温度时材料为韧 性断裂。运用原子模拟方法，人们已模拟 得到低温时bcc-Fe中裂纹脆性解理扩展过 程，表明bec．Fe是本征脆性金属【4jI。另外， 图1原子．有限元模型 ofFEAt 人们在材料的低温脆性断裂实验中已普遍 Fi91．Model 观测到挛晶的存在【6t7】。随后的裂尖分子动 612 本文采用的原子．有限元模型如图1所 在实际模拟过程中，有限元区的大小 示。模型共分为：乏个区，即原子区，又可 为2u×2u，相当子纳米晶的尺寸。垂直 称为有限单元区l，及有限单元区II和有限 F纸面厚度取无穷大，为平面应变条件。 原子区含160×160×6层原子，尺寸为 单元区III。原子区即分子动力学(MD) 区，由离散的粒子，即原子组成，处于模 O0321．t×0．023H，共含76632个原子，厚 型内部的裂尖区域。同时，原子区也为有 度方向采用周期性边界条件。初始裂纹通 限元的有机组成部分，即所有原子均与有 过挖去原子区2层原子而产生，一共挖去 限元中的网格节点重合，并形成有限元区 168个原子，裂纹在原子区长度为28层原 I。有限单元区II为精细的有限单元区，其 子。同时，有限元区也产生相应的初始裂 网格的大小由原子尺度向外逐渐加大，使 缝。本文研究bcc-Fe中{100}01l边界I 得有限单元在靠近原子区的过渡较为平 型裂纹，其空间构型如图2所示，即裂纹 缓。而有限单元区III为粗网格区，在此区 表面为(010)面，裂纹前缘为[10H