基于等几何分析与半解析框架的激光粉末床熔融高效热模拟新方法

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：Efficient thermal simulation in metal additive manufacturing via semi-analytical isogeometric analysis

为解决激光粉末床熔融(LPBF)热模拟中传统有限元法(FEM)计算成本高、解析法边界处理难的问题，研究人员开展了基于等几何分析(IGA)的半解析热模拟研究。该研究通过将温度场分解为解析点源解与IGA修正场，实现了对复杂几何零件的高效、高精度温度预测，显著提升了计算效率。

　　

在金属增材制造（Additive Manufacturing, AM）的竞技场上，激光粉末床熔融（Laser Powder Bed Fusion, LPBF）技术凭借其制造复杂几何结构的能力备受青睐。然而，在微观层面，激光与粉末的相互作用犹如一场“热力闪电战”——激光束以极快的速度移动，在微米级的熔池内产生瞬间的极高温度梯度。想要在计算机中精准复现这一物理过程，对工程师来说一直是个巨大的挑战。

传统的有限元法（Finite Element Method, FEM）虽然应用广泛，但为了捕捉激光带来的剧烈温度变化，往往需要在激光光斑附近布置极其细密的网格（通常要求网格尺寸与光斑半径相当，甚至细化到5–25微米）。这种“杀鸡用牛刀”的策略导致计算量激增，使得零件尺度的仿真耗时漫长，甚至难以完成。虽然动态自适应网格（如扫描级重网格）能缓解这一问题，但频繁的网格更新和矩阵重组本身又带来了新的计算开销。

另一方面，纯解析方法（如点热源解）虽然能快速算出温度场，但它假设材料是无限大的，完全忽略了零件实际边界（如散热面、尖角）的影响。当激光扫描到零件边缘时，这种“无限大”的假设会导致计算结果严重失真。半解析方法试图取二者之长，用解析解处理激光主效应，用数值解补充边界条件。但现有的半解析方法在复杂边界处往往“水土不服”，要么被迫在边界处进行精细网格加密，要么依赖繁琐的“镜像热源法”，难以处理真实零件中常见的复杂轮廓和尖角。

面对这一困境，Yang Yang、Ye Ji、Matthias M?ller和Can Ayas团队在《Computer Methods in Applied Mechanics and Engineering》上发表的研究中，提出了一种融合等几何分析（Isogeometric Analysis, IGA）的半解析新框架，旨在打破计算效率与几何精度之间的壁垒。

关键技术方法

本研究核心是将温度场分解为已知的解析点源解（满足控制方程）与待求的IGA修正场（满足边界条件）。针对修正场，建立了基于非均匀有理B样条（NURBS）的弱形式控制方程，利用IGA的高阶连续性基函数进行空间离散，并采用隐式θ-时间积分方案推进求解。该方法利用IGA精确的几何描述能力，避免了扫描过程中的重网格操作，实现了对复杂几何边界的高效处理。

研究结果与发现

1. 方法框架构建

研究人员构建了一种新颖的数学“分身术”。他们将总温度场（T）拆解为两部分：一部分是描述激光热输入的解析点源解（T_a），它自带捕捉陡峭梯度的能力；另一部分是负责“修正”边界条件的互补场（T_c），通过IGA求解。由于T_c场的变化相对平缓，IGA网格可以完全摆脱激光尺度的束缚，采用较粗的离散化即可满足精度要求，从而大幅减少了自由度数量。

2. 复杂几何适应性验证

针对传统镜像法难以处理的几何特征（如锐角、变截面），该方法展现了强大的鲁棒性。在包含尖角的几何算例中，IGA凭借其精确的几何描述能力，无需像传统FEM那样通过加密网格来“拟合”曲线或尖角，有效避免了在边界处因解析解失效而被迫进行的局部加密，确保了计算效率的稳定性。

3. 计算效率对比

数值实验表明，在达到相同预测精度的情况下，该半解析IGA方法所需的计算资源显著低于标准FEM。得益于IGA更高的单位自由度精度和避免了扫描级重网格，该方法在处理复杂轮廓零件时，实现了计算效率的数量级提升，为零件级快速热模拟提供了可行路径。

结论与意义

这项研究成功地将等几何分析（IGA）与半解析热模拟相结合，为LPBF热管理提供了一种高效、高精度的新工具。其重要意义在于：一是突破了传统数值方法在复杂几何边界处的处理瓶颈，利用IGA的精确几何表达避免了繁琐的镜像源构造；二是通过解耦激光尺度与网格尺度，大幅降低了计算成本，使快速预测零件制造过程中的温度历史、优化工艺参数成为可能，对推动增材制造的数字化设计与质量控制具有重要价值。