金属增材制造的热分析
金属增材制造热的理论基础
概述
老师！ 今天我们要讲金属增材制造的热分析对吧？ 它是什么样的？
激光PBF/DED的温度场预测。熔池、残留应力。
控制方程
听到这里，我终于明白金属增材制造热分析为什么重要了！
离散化方法
这个方程在计算机上具体怎样求解呢？
使用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构建总体刚度方程。
矩阵求解算法
矩阵求解算法具体是怎样的？
通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解线性方程组。对于大规模问题，预处理迭代法很有效。
| 求解法 | 分类 | 内存使用 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中等规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小～中等规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说有限元法阶段如果草率的话，后面会很吃亏！ 我会铭记在心的！
商用工具中的实现
那么做金属增材制造热分析的话，有什么软件可以用呢？
| 工具名称 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
供应商系统及产品整合历史
各个软件的发展历史是不是很有戏剧性？
Ansys Mechanical (旧ANSYS Structural)
请告诉我关于"Ansys Mechanical"的信息！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys Parametric Design Language）。
现在的所属：Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEA具体是怎样的？
1978年由HKS (Hibbitt, Karlsson & Sorensen)开发。2005年被Dassault Systèmes收购，并整合到SIMULIA品牌。
现在的所属：Dassault Systèmes SIMULIA
听到这里，我终于明白开发为什么重要了，终于理清了！
COMSOL Multiphysics
请告诉我关于"COMSOL Multiphysics"的信息！
1986年在瑞典成立。最初作为MATLAB联动的FEMLAB开始，后来改名为COMSOL。擅长多物理领域。
现在的所属：COMSOL AB
哇~，开发的故事超级有意思！ 还想听更多。
文件格式和互操作性
不同软件之间交换数据时有什么注意事项吗？
| 格式 | 扩展名 | 类别 | 概述 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | 遵守ISO 10303的3D CAD数据交换格式。形状+PMI支持。 |
| IGES | .igs/.iges | 中立CAD | 早期CAD数据交换规范。曲面数据兼容性有问题。正在向STEP迁移。 |
在不同求解器间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载与边界条件的表现差异。尤其是高阶单元或特殊单元（内聚单元、用户定义单元等）往往无法在求解器间直接转换。
格式看起来很简单，但实际上深藏玄机呢。
实务注意事项
有什么教科书里没有的"现场智慧"吗？
网格收敛性验证、边界条件的合理性检验、材料参数的敏感性分析非常重要。
哇，金属增材制造热分析真是非常深奥呢… 但在老师的讲解下，我整理得差不多了！
好的，很有进度！实际动手操作是最好的学习方式。有不明白的地方随时来问我。
金属3D打印的热历史问题
金属AM的热历史与1984年查尔斯·胡尔获得SLA专利的光固化截然不同。在粉床融合（PBF）中，每一层都要经历急热急冷（冷却速率10³~10⁷℃/s）的反复循环。这种独特的热历史产生残留应力和特殊的微细组织（柱状晶）。2016年Lawrence Livermore国家实验室公开的"DigiMat-AM"是世界首个大规模AM热分析框架。
金属增材制造热的数值计算方法
数值方法详解
具体用什么算法来求解金属增材制造热分析呢？
离散化的表述
用形状函数 $N_i$ 近似未知量：
用数式表示是这样的。
基本方程的离散形式
用数式表示是这样的。
嗯，只看公式有点摸不着头脑… 它表示什么呢？
连续体的控制方程离散化后，得到下述代数方程组：
这里 $[K]$ 是全体刚度矩阵（或等价的系统矩阵），$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。
啊，这样就明白了！连续体控制方程就这样变成了方程组啊。
单元技术
听说过"单元技术"，但可能理解得不太准确…
| 单元类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中 |
积分方案
积分方案具体是怎样的？
听到这里，我终于明白单元类型为什么重要了，终于理清了！
收敛性和稳定性
如果不收敛了，首先要检查什么？
收敛速度：对于二阶单元，误差以 $O(h^2)$ 阶递减（光滑解的情况）
看似简单的网格细分，实际上深藏玄机呢。
求解器设置建议
具体用什么算法来求解金属增材制造热分析呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理方法 | ILU(0) or AMG | 根据问题规模而定 |
| 最大迭代次数 | 1000 | 不收敛时需重新检查设置 |
| 内存模式 | In-core | 尽可能采用 |
线性单元 vs 2阶单元
热传导分析中线性单元通常能给出足够的精度。在温度梯度急剧的区域（热冲击等）推荐使用2阶单元。
热流通量的评估
从单元内的温度梯度计算。与节点应力一样，通常需要平滑处理。
对流-扩散问题
当Peclet数高（对流主导）时需要迎风稳定化（SUPG等）。纯热传导问题无需考虑。
非定常分析的时间步长
热扩散的特征时间 $\tau = L^2 / \alpha$（$\alpha$：热扩散率）应该远大于时间步长。对于急剧温度变化，自适应时间步长控制很有效。
非线性收敛
由温度相关材料性质引起的非线性通常相对温和，Picard迭代（直接替代法）通常足够。对于强非线性的辐射，推荐使用Newton法。
定常分析的判定
全节点温度变化小于阈值（$|\Delta T| / T_{max} < 10^{-5}$等）时判定为收敛。
金属增材制造热的实务应用
实践指南
老师，请教我"实践指南"！
金属增材制造热分析的实务分析流程和注意点详细讲解。
分析流程
从最初阶段教我！ 首先应该做什么？
1. 预处理 (前处理)
- CAD数据导入与形状简化
- 材料特性定义
- 网格生成（单元类型·大小的确定）
- 边界条件和荷载条件的设置
2. 求解 (求解)
- 求解器设置（求解方法、收敛基准、输出控制）
- 作业投入与计算执行
- 收敛监控
3. 后处理 (后处理)
- 结果的可视化（位移、应力、其他物理量）
- 结果验证与妥当性确认
- 报告编制
网格生成最佳实践
怎样判断网格的好坏？
单元品质指标
请告诉我"单元品质指标"！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度降低 |
| 雅可比比 | 1.0 | > 0.3 | 单元退化 |
| 翘曲 | 0° | < 15° | 精度降低 |
| 歪斜度 | 0° | < 45° | 收敛性恶化 |
| 锥度比 | 0 | < 0.5 | 精度降低 |
网格密度的确定
网格密度的确定具体是怎样的？
边界条件设置指南
听说边界条件设置错的话全部都白做…
啊，这样啊！ 过约束注意原来就是这个机制啊。
商用工具分别实现步骤
有很多软件对吧？ 各自的特点请告诉我！
| 工具名称 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
Ansys Mechanical (旧ANSYS Structural)
请告诉我关于"Ansys Mechanical"的信息！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys Parametric Design Language）。
现在的所属：Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEA具体是怎样的？
1978年由HKS (Hibbitt, Karlsson & Sorensen)开发。2005年被Dassault Systèmes收购，并整合到SIMULIA品牌。
现在的所属：Dassault Systèmes SIMULIA
听到这里，我终于明白开发为什么重要了，终于理清了！
常见失败与对策
新手经常犯什么错误？ 想提前知道！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、不恰当的边界条件 | 改善网格、检查约束条件 |
| 应力异常大 | 应力特异点、网格依赖 | 避开特异点、局部网格细分 |
| 位移不现实 | 材料常数错误、单位系统混乱 | 确认输入数据 |
| 计算时间过长 | 不必要的细分、低效的求解 | 网格优化、并列计算 |
质量保证检查表
有什么教科书里没有的"现场智慧"吗？
哇，金属增材制造热分析真是非常深奥呢… 但在老师的讲解下，我整理得差不多了！
好的，很有进度！实际动手操作是最好的学习方式。有不明白的地方随时来问我。
波音AM钛合金零件解析
2018年，波音公司在使用GE Additive制造的EBM机器进行钛合金（Ti-6Al-4V）支架的变形分析时，采用3DExperience Simulia AM进行评估。预测无需支撑结构即可制造，实际制造成功。零件重量比传统锻造品减少55%，同时尺寸精度达到±0.2mm，成为AM行业教科书式的成功案例。
金属增材制造热的软件比较
商用工具比较
有很多软件对吧？ 各自的特点请告诉我！
金属增材制造热分析对应的主要商用CAE工具的功能比较和各产品的历史背景详细说明。
支持工具列表
那么做金属增材制造热分析的话，有什么软件可以用呢？
| 工具名称 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
Ansys Mechanical (旧ANSYS Structural)
请告诉我关于"Ansys Mechanical"的信息！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys Parametric Design Language）。
现在的所属：Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEA具体是怎样的？
1978年由HKS (Hibbitt, Karlsson & Sorensen)开发。2005年被Dassault Systèmes收购，并整合到SIMULIA品牌。
现在的所属：Dassault Systèmes SIMULIA
听到这里，我终于明白开发为什么重要了，终于理清了！
COMSOL Multiphysics
请告诉我关于"COMSOL Multiphysics"的信息！
1986年在瑞典成立。最初作为MATLAB联动的FEMLAB开始，后来改名为COMSOL。擅长多物理领域。
现在的所属：COMSOL AB
Ansys Fluent
接下来是Ansys Fluent的话题吧。什么内容？
Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。
现在的所属：Ansys Inc.
啊，这样啊！ 开发就是这样的机制啊。
功能比较矩阵
预算和时间都有限，最划算的是哪一个？
| 功能 | Ansys Mechanical | Abaqus | COMSOL | Fluent |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并列计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换的风险
转换的风险具体是怎样的？
啊，这样啊！不同工具间的转换就是这样的机制啊。
许可证形式
听说过"许可证形式"，但可能理解得不太准确…
| 工具 | 许可 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 价格高但有官方支持 |
| OpenFOAM | GPL | 免费但支持为收费 |
| COMSOL | 节点锁定/浮动 | 按模块单位购买 |
| Code_Aster | GPL | EDF开发的OSS求解器 |
选择指南
最后，选哪一个的判断标准请告诉我？
金属增材制造热分析工具选择时需考虑以下几点：
哇，金属增材制造热分析真是非常深奥呢… 但在老师的讲解下，我整理得差不多了！
好的，很有进度！实际动手操作是最好的学习方式。有不明白的地方随时来问我。
Amphyon（Additive Works）的战略
德国初创公司Additive Works（2015年成立）开发的AM专用仿真软件Amphyon最早在业界实现了固有应变法，2017年在Formnext展览会上获得大奖。2020年被Ansys收购并整合为"Ansys Additive Print"。当前与Ansys Simulation Platform的工作流集成，在航空航天、医疗领域快速推广。
金属增材制造热的前沿研究
前沿话题和研究动向
金属增材制造热分析领域今后怎样发展？
金属增材制造热分析的最新研究动向和先进手法详细分析。
最新的数值方法
接下来是最新的数值方法的话题吧。什么内容？
嗯，只看公式有点摸不着头脑… 它表示什么呢？
高性能计算 (HPC) 的支持
| 并列化方法 | 概述 | 适用求解器 |
|---|---|---|
| MPI (区域分割) | 分布式内存型。大规模问题标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并列 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。特别对显式法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并列 | 大规模HPC环境 |
金属增材制造热的故障排除
故障排除
常见错误和对策
老师在金属增材制造热分析中熬过夜间调试吗？（笑）
1. 收敛失败
收敛失败具体是怎样的？
症状：求解器在指定迭代次数内无法收敛而异常终止
可能的原因：
- 网格品质不足（过度变形的单元）
- 材料参数设置不恰当
- 初始条件不合适
- 非线性过强（荷载步不足）
对策：
- 实施网格品质检查（纵横比、雅可比）
- 确认材料参数的单位系
- 将荷载分为多个步（增加子步数）
- 放宽收敛判定基准（注意精度）
也就是说收敛失败地方如果草率的话，后面会很吃亏！ 我会铭记在心的！
2. 非物理的结果
接下来是非物理的结果的话题吧。什么内容？
症状：应力/位移/温度等物理上不现实
可能的原因：
- 边界条件设置错误
- 单位系混乱（SI单位与工程单位混合）
- 不恰当的单元类型选择
- 应力特异点的存在
对策：
- 确认反力合计（力的平衡）
- 检查单位系的一致性
- 重新考虑单元类型的适切性
- 特异点消除或子模型化
前辈说"收敛失败一定要好好做"的含义现在明白了。
3. 计算时间超过
计算时间超过具体是怎样的？
症状：计算耗时远超预期
对策：
- 优化网格的粗细分布
- 活用对称性（1/2、1/4模型）
- 优化求解器设置（迭代法、预处理选择）
- 活用并列计算
4. 内存不足
请告诉我"内存不足"！
症状：内存溢出错误
前辈说"收敛失败一定要好好做"的含义现在明白了。
对策：
- 使用核外求解
- 减少网格规模
- 确认使用64位版求解器
- 增加内存分配
哇~，收敛失败的话题，超级有意思！ 还想听更多。
Nastran代表错误
代表错误具体是怎样的？