《Journal of Materials Research and Technology》：Shock-induced spallation behavior and deformation mechanisms of nanotwinned Cu: MD simulations and partial experimental validation
编辑推荐：
本研究结合分子动力学（Molecular Dynamics, MD）模拟与部分实验表征，系统研究了[111]织构纳米孪晶铜（nanotwinned Cu）在不同加载速度和加载方向下的冲击诱导层裂行为及潜在损伤机制。模拟结果揭示，在低冲击速度下存在显著的层裂各向
本研究结合分子动力学（Molecular Dynamics, MD）模拟与部分实验表征，系统研究了[111]织构纳米孪晶铜（nanotwinned Cu）在不同加载速度和加载方向下的冲击诱导层裂行为及潜在损伤机制。模拟结果揭示，在低冲击速度下存在显著的层裂各向异性行为，而随着冲击强度的增加，加载方向的影响逐渐减弱。在0°和45°加载条件下，孔洞主要沿晶界形核与长大，呈现层状合并特征；相反，在90°加载条件下，孔洞主要沿孪晶界形核，呈独立球形分布。X射线衍射（XRD）模拟与实验衍射结果一致表明，位错滑移主导动态变形过程，而非孪生和晶粒转动。透射电子显微镜（TEM）观察进一步证实了冲击加载后强烈的位错累积、层错形成以及位错-孪晶界相互作用。由于纳米孪晶的强化效应，织构纳米孪晶铜的动态性能与单晶铜相当。在45°加载条件下，裂纹扩展偏离最敏感路径，导致其具有显著的层裂强度；而在90°加载条件下，高密度位错的相互作用阻碍了位错的进一步运动，导致其层裂强度最低。
冲击诱导纳米孪晶铜层裂行为与变形机制研究解读
一、研究背景与科学问题
在高速交通、航空航天、国防技术、核能应用及天体物理等领域，材料经常承受高速冲击所引起的高温、高压和高应变率极端环境。在此条件下，理解材料的失效阈值和损伤机制对于评估和优化其性能至关重要。已有研究表明，材料经历典型的高速平面冲击时，主要经历压缩、卸载和拉伸三个阶段。对回收样品的微观分析证实，金属材料的失效过程涉及孔洞的萌生、长大和合并，最终导致层裂。然而，初始微观结构特征（如晶粒尺寸、晶体学织构和孪晶界特性）及其在冲击加载过程中的演化对位错活动、孔洞演化和层裂强度的耦合影响仍认识不足。
一方面，即使借助高时空分辨率的X射线原位诊断技术，对单晶、多晶和纳米晶等典型微观结构对层裂影响的实时观察仍面临重大挑战。另一方面，在缺乏原位微观诊断工具的情况下，分子动力学（Molecular Dynamics, MD）模拟已被广泛用于研究微观结构对层裂机制的影响，同时可通过X射线衍射模拟实现模拟与实验之间的桥梁连接。
近年来，合成的纳米孪晶铜展现出独特的织构和性能，具有超高强度和高导电性。已有研究发现其屈服强度与加载方向密切相关，可实现三种位错模式的有效切换。然而，关于纳米孪晶铜层裂机制的研究仍十分有限，特别是压缩和卸载引起的微观结构演化对动态拉伸过程中层裂行为的影响，与准静态拉伸失效研究相比鲜有涉及。
二、研究内容与结论概述
基于前期研究，本研究结合MD模拟和部分冲击回收实验，系统研究了织构纳米孪晶铜在动态加载下的层裂行为和变形机制，重点关注孔洞形核与长大，以及加载方向和孪晶界对层裂的影响。研究结果表明，在低冲击速度下存在明显的层裂各向异性，而随着冲击强度增加，加载方向的影响逐渐减弱。45°加载条件下裂纹扩展偏离最敏感路径，导致最高的层裂强度；90°加载条件下高密度位错的相互作用阻碍位错进一步运动，导致最低的层裂强度。由于纳米孪晶的强化效应，织构纳米孪晶铜的动态性能与单晶铜相当。该研究为具有特殊结构材料的层裂机制提供了有价值的见解，有助于深入理解微观结构特征对极端条件下损伤抗性的影响。
三、关键技术方法
研究人员构建了具有[111]织构的纳米孪晶铜模型，采用Voronoi法则构造20×20×40 nm的模拟盒，包含8个晶粒，孪晶层厚度设为2 nm，经旋转后模型包含约135万个原子。随后沿X、Y、Z方向复制并按照0°、45°和90°三种加载角度裁剪，最终模型尺寸为20×20×360 nm，包含约1224万个原子。采用共轭梯度法和NPT系综退火进行结构优化，使用Mishin等人提出的嵌入原子法（Embedded Atom Method, EAM）势函数描述原子间相互作用，通过活塞法施加冲击加载，选取0.1、0.3、0.5、1.0、1.5和2.5 km/s六种加载速度。所有MD模拟由LAMMPS模拟器完成。微观结构分析采用Ovito软件进行共近邻分析（Common Neighbor Analysis, CNA）、位错分析（Dislocation Analysis, DXA）、构造表面网格（Constructed Surface Mesh, CSM）和模板匹配（Polyhedral Template Matching, PTM）。X射线衍射（XRD）图谱使用GPU加速多色衍射（GPU-Accelerated Polychromatic Diffraction, GAPD）代码模拟。实验方面，通过直流电沉积制备[111]织构纳米孪晶铜样品，使用14 mm口径单级轻气炮进行冲击加载-卸载实验，飞片为高导电无氧铜（OFHC Cu），样品通过软回收舱回收。
四、研究结果
3.1 纳米孪晶铜的冲击响应
通过自由表面速度（Free Surface Velocity, FSV）波剖面分析发现，在0.1 km/s冲击速度下，冲击波约在70-80 ps到达自由表面，自由表面速度迅速增加至0.2 km/s，但由于冲击速度较低，未显示卸载和再加载波形。在0.3 km/s冲击下，0°和45°加载方向在压缩阶段呈现明显的弹塑性转变点，而90°加载未观察到该转变。45°加载下卸载持续时间最长（约60 ps），波形演化最缓慢。随着冲击速度增加至1.0 km/s以上，各加载方向间的波形差异开始减小，在1.5 km/s和2.5 km/s时波形曲线趋于收敛。时间-位置-法向应力云图分析表明，样品对冲击波的响应可分为四个阶段：冲击压缩阶段、卸载阶段、拉伸阶段和复杂状态阶段。在0.5 km/s冲击下，峰值拉应力分别为0°加载11.11 GPa、45°加载12.05 GPa和90°加载11.63 GPa；在1.0 km/s冲击下，峰值拉应力分别为0°加载12.12 GPa、45°加载12.83 GPa和90°加载11.86 GPa。45°加载在三种取向中表现出最高的峰值拉应力。
3.2 纳米孪晶铜的层裂强度
研究人员根据自由表面速度回拉，利用公式计算了有效层裂强度。体声速对冲击速度高度敏感，呈正相关关系。回拉速度不随冲击速度呈简单单调依赖。三种加载方向呈现一致趋势：层裂强度在研究的冲击速度范围内随冲击速度增加而增大。最低层裂强度出现在0.3 km/s冲击速度下，最高出现在2.5 km/s。三种加载方向中，45°加载层裂强度最高，0°加载次之，90°加载最低。在2.5 km/s和45°加载条件下，织构纳米孪晶铜的层裂强度达到峰值14.59 GPa。在0.1 km/s冲击速度下未观察到层裂强度，因为未形成孔洞。与单晶铜在1.0 km/s冲击下层裂强度12.2 GPa相比，织构纳米孪晶铜的层裂强度与之相当且高于纳米晶铜。纳米孪晶的强化效应可归因于三个因素：位错屏障作用、集中位错束的存在和共格孪晶界的热稳定性。
3.3 孔洞与微观结构演化
3.3.1 冲击速度的影响
在0.5 km/s冲击下，样品在135 ps时呈现显著的结构损伤，形成明显裂纹。0°加载下孔洞沿晶界呈层状合并，层状裂纹平行于晶界；90°加载下孔洞主要呈独立球形分布，且位错密度显著高于0°和45°加载。随着冲击速度增加至1.0 km/s，所有加载方向下位错密度和孔洞数量均大幅增加。应变率随冲击速度增加而增大，高应变率导致材料内位错密度显著上升，激活更多滑移系，形成更多形核位点，促进孔洞形核与长大。0°加载0.5 km/s冲击下的切片分析揭示了孔洞演化的两个阶段：快速孔洞形核阶段和孔洞长大与合并阶段。剪切应变集中在晶界处，孔洞周围呈现明显的应变局部化。
3.3.2 冲击方向的影响
在0°加载下，加载方向平行于孪晶界，孔洞在晶界处形核并沿晶界长大，断裂主要发生在晶间区域。在45°加载下，孔洞主要沿晶界分布并沿该方向长大，孔洞倾向于沿冲击加载方向向晶粒内部扩展，裂纹扩展路径被迫偏离最敏感轨迹，增加了裂纹扩展所需能量，这解释了45°加载下层裂强度最高的原因。在90°加载下，孔洞主要沿孪晶界形核并长大，晶界对孔洞长大的影响显著降低。90°加载下损伤程度明显小于0°和45°加载。X射线衍射模拟显示，初始样品具有明显的(111)衍射环，确认了(111)织构。层裂后(111)衍射环发生轻微变形，但未发现新衍射环，表明晶粒转动不是冲击断裂过程中的主要变形机制。衍射斑点向外移动并变得弥散和拉长，表明微观结构内形成了缺陷。未观察到破碎衍射斑点，表明断裂过程中未发生孪生。0°和45°加载下出现衍射斑点消失和新斑点出现，表明发生了塑性变形引起的再取向；90°加载下衍射斑点基本未出现明显消失或出现，表明晶格再取向有限。
3.4 微观结构定量分析
孔洞体积分数演化曲线表明，孔洞体积分数随卸载时间延长而增加。0.5 km/s冲击下，0°加载的孔洞体积始终高于45°和90°加载，45°加载的孔洞体积分数最低。1.0 km/s冲击下，三种加载方向的孔洞体积分数差异减小，但0°加载仍略高。位错密度分析显示，0.5 km/s冲击下位错密度先增加后降低，90°加载下位错密度始终最高。1.0 km/s冲击下位错密度在初始记录时间已处于高位，随后逐渐降低，90°加载下仍保持最高。六方密排（Hexagonal Close-Packed, HCP）原子含量分析显示，0.5 km/s冲击下0°和90°加载的HCP原子含量总体下降，而45°加载的HCP原子含量在层裂过程中总体增加并逐渐超过其他方向。1.0 km/s冲击下各加载方向的HCP原子含量均随层裂时间延长而逐渐下降，45°加载下HCP原子含量最高，表明更多层错参与变形。90°加载下位错密度远高于0°和45°加载，高密度位错间的相互作用阻碍位错进一步运动，降低了材料的协调变形能力，从而导致层裂强度最低。
五、部分实验验证与讨论
实验制备的纳米孪晶铜样品厚度约59.3 μm，厚度均匀性良好。截面扫描电子显微镜（SEM）形貌显示沿厚度方向生长的柱状晶结构，与MD模拟中构建的织构纳米孪晶铜模型高度一致。电子背散射衍射（EBSD）分析显示Σ3孪晶界占全部晶界的83.5%，表明材料中存在高密度共格孪晶界。反极图（IPF）图显示孪晶主要分布在晶粒内部，属于晶内纳米孪晶。极图分析确认了强烈的(111)择优取向。平均晶粒尺寸约0.8 μm，范围为0.2～5.1 μm，表明"纳米孪晶铜"主要指晶粒内部的纳米尺度孪晶片层间距，而非晶粒尺寸本身。
冲击回收样品的XRD图谱显示所有衍射峰均可标定为面心立方（Face-Centered Cubic, FCC）铜相，未出现新衍射峰，表明冲击加载过程中未发生相变。回收样品仍表现出强烈的Cu(111)衍射特征，表明原始[111]织构在冲击后基本保留。不同加载方向下可观察到峰强度和峰展宽的明显变化，表明动态变形过程中产生了残余应变和高密度缺陷。0°和45°加载下衍射峰强度变化更明显，90°加载下衍射特征相对稳定，表明垂直于加载方向的孪晶界可部分抑制严重的晶格再取向和缺陷演化。这些实验观察与XRD模拟结果一致。
六、研究结论
本研究通过MD模拟与部分实验表征相结合，系统研究了[111]织构纳米孪晶铜在不同加载速度和加载方向下的冲击诱导层裂行为和变形机制。模拟结果揭示了低冲击速度下明显的各向异性层裂行为，而加载方向的影响随冲击强度增加而逐渐减弱。0°和45°加载下孔洞主要沿晶界形核与长大，呈层状合并；90°加载下孔洞主要沿孪晶界形核，呈独立球形分布。XRD模拟与实验衍射结果一致表明位错滑移主导动态变形过程，而非孪生和晶粒转动。TEM观察进一步证实了冲击加载后强烈的位错累积、层错形成和位错-孪晶界相互作用。由于纳米孪晶的强化效应，织构纳米孪晶铜的动态性能与单晶铜相当。45°加载下裂纹扩展偏离最敏感轨迹，赋予其显著的层裂强度；90°加载下高密度位错的相互作用阻碍位错进一步运动，导致其层裂强度最低。
