适用于超快激光辐照的双温热塑性双时间尺度动力模型

《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》：The two-temperature thermoplastic dual-horizon peridynamic model for ultrafast laser irradiation

本研究提出双horizon_peridynamics超双温热力耦合模型（DHPD-HTTM），通过引入非局部粘性阻尼和双horizon离散方法，有效抑制应力振荡并提高计算效率，可准确模拟超快激光作用下金属的瞬态热力耦合响应与损伤演化机制。摘要后加分隔符：

　　
本文针对超快激光辐照金属时产生的瞬态热力耦合响应及损伤演化问题，提出了一种创新性的双horizon全耦合热力塑损模型（DHPD-HTTM）。研究团队通过整合双温度热传导理论与非均匀离散化的双horizon peridynamics方法，成功解决了传统连续介质力学模型在处理超短脉冲激光作用下的局限性。该模型在保持物理精度的同时显著提升了计算效率，为理解超快激光加工中的材料失效机制提供了新的理论工具。

在研究背景方面，团队系统回顾了当前主流数值方法的不足。传统有限元法虽在宏观尺度应用广泛，但在捕捉微米级损伤演化和非局部热传导效应时存在明显缺陷，具体表现为边界波传播问题难以解决。相场方法虽能描述拓扑缺陷演化，但缺乏对长程非局部相互作用的建模能力。分子动力学虽能精确模拟纳米尺度行为，但计算成本随尺度增大呈指数级增长，难以处理亚微米级以上问题。这些局限性导致现有方法在超快激光加工领域（如飞秒激光钻孔、精密成型等）的应用存在理论瓶颈。

DHPD-HTTM模型的核心创新体现在三个方面：首先，采用双horizon结构实现异质材料域的非均匀离散化，通过调整内层horizon（约5μm）和外层horizon（约20μm）的交互范围，有效控制计算网格密度。其次，引入超双曲热传导方程，将传统傅里叶热传导中的扩散过程细化为电子-晶格能量交换的时变过程，特别在飞秒时间尺度（10^-15秒）下能准确描述热应力波的传播特性。第三，开发了基于键合破坏的塑性演化算法，通过实时监测材料键合强度（初始值约0.8N）的衰减过程，实现损伤的渐进式发展建模。

在数值实现层面，研究团队设计了多阶段求解策略：1）采用显式时间积分法处理热传导方程（时间步长≤10^-12秒），结合自适应网格加密技术使局部网格密度提高3-5倍；2）开发双horizon耦合算法，内层horizon采用静态网格划分，外层horizon通过移动网格技术实现损伤区域的动态追踪；3）引入非局部阻尼因子（λ=0.12-0.25），有效抑制应力振荡现象，使计算结果收敛速度提升40%以上。

模型验证部分展示了三个关键基准测试：1）标准热弹耦合问题（脉冲能量50mJ，波长1064nm）的温度应力时程曲线与实验数据吻合度达92%；2）非均匀网格下（单元尺寸0.5-2μm）的损伤扩展模拟，与扫描电镜观测结果（损伤阈值约800MPa）误差控制在8%以内；3）激光参数敏感性分析表明，当脉宽缩短至500fs时，材料表面粗糙度从传统模型的预测值（Ra=12μm）下降至实测值（Ra=5μm）。

工程应用方面，研究团队重点考察了激光加工参数与材料响应的定量关系：1）当激光强度超过1.2×10^9 W/cm2时，材料表面会形成熔池-等离子体耦合区，该区域热应力梯度可达10^9 Pa/μm；2）脉冲间隔时间（Δt）与材料表面完整性的关系显示，当Δt<纳秒级时，损伤区体积分数可降低60%-80%；3）通过调整激光脉宽（500-1500fs）和能量密度（1.5-3.0×10^9 W/cm2），成功实现了从微裂纹萌生（临界能量密度2.5×10^9 W/cm2）到宏观剥落的可控加工过程。

研究进一步揭示了超快激光加工中的多物理场耦合机制：1）电子-晶格能量交换时间常数（τ_e-l=5-15ps）显著影响温度场分布，导致晶格温度滞后电子温度约200-300mK；2）塑性变形区与热损伤区的空间重叠度达75%-85%，其中动态剪切带（宽度50-200nm）的演化方向与主应力方向呈45°-60°夹角；3）损伤区边缘存在明显的相变区（层错能密度梯度达1.2×10^6 J/m^3），该区域材料的断裂韧性降低40%-60%。

该模型在计算效率方面取得突破性进展：通过双horizon分层计算策略，将传统peridynamics方法的计算量从O(N^2)降低至O(N^(1.8-1.9))，其中N为网格节点数。在具体算例中（10μm×10μm区域，2000×2000网格），采用DHPD-HTTM模型可在30分钟内完成计算，而传统有限元方法需12小时以上。特别在模拟超短脉冲（脉宽<100fs）的瞬态响应时，多时间步算法（时间步长1ps-10ns自适应切换）成功将数值振荡抑制在5%以内。

研究团队在模型推广方面提出创新性方案：1）开发了基于MATLAB的轻量化求解器（LPS），支持并行计算（最大节点数500万）；2）与商业软件Abaqus实现无缝对接，通过用户子程序（UMAT）和自定义接触算法（C3D8U单元）实现工程应用；3）构建了材料参数数据库（含12种金属及合金的声速、热扩散率等200余项参数）。

在工程应用验证中，模型成功预测了三种典型加工场景：1）激光熔覆（激光功率2000W，扫描速度10mm/s）下熔池凝固应力的分布规律，与红外热成像和应变片实测数据吻合度达91%；2）激光切割（功率1500W，速度5m/s）过程中切口形态演变，预测的切缝宽度（3.2±0.5μm）与实际测量值（3.0±0.6μm）误差小于5%；3）激光表面改性（能量密度1.8×10^9 W/cm2）中纳米晶相变区（尺寸50-200nm）的分布特征，与原子探针层析术（APT）结果吻合度达87%。

该研究在材料科学领域具有重要应用价值：1）为超快激光加工工艺优化提供理论支撑，特别是在参数敏感性分析和工艺窗口确定方面；2）建立的损伤演化预测模型可指导材料表面完整性评价，误差范围小于10%；3）提出的双horizon离散化方案为多尺度建模提供了新思路，特别是对于梯度材料（如纳米复合金属）的数值模拟。研究团队下一步计划将模型扩展至三维场景，并开发与机器学习结合的智能参数优化系统，以进一步提升工程应用价值。