三维定常热传导
三维定常热传导的理论基础
三维定常热传导的基础
三维热传导分析和二维相比有什么不同？
温度场在三个方向上变化，所以形状简化在这类问题上无效。引擎缸体、模具、电子机箱等现实中的许多问题本质上是三维的。
控制方程
三维定常热传导方程为
展开后得
对于各向异性材料，有 $k_x \neq k_y \neq k_z$。
三维的解析解基本上不存在吧？
对于给定每个面都具有简单边界条件的长方体，存在三重级数解，但在实际应用中数值求解是必需的。使用FEM固体单元（四面体、六面体）进行离散化是标准方法。
边界条件的种类
在三维问题中，需要对每个面设置不同的边界条件。
| 面 | 条件 | 例子 |
|---|---|---|
| 发热面 | 热流密度 q [W/m2] | IC发热面 |
| 散热面 | 对流 h, T∞ | 散热片外表面 |
| 接触面 | 接触导热系数 | 螺栓连接部分 |
| 对称面 | 绝热 (q=0) | 对称利用 |
| 固定温度 | T = const | 冷却水面 |
对不同的面进行不同条件设置是三维的优势。
没错。在一维中只能整体处理，而在三维中可以对各个面和各个区域设置单独条件，得到物理上更准确的模型。
三维热传导方程的一般形式
三维定常热传导∇·(λ∇T)+q̇=0，当λ是各向同性均质且无内部发热时，可化为拉普拉斯方程。Green（1828）引入的"格林函数"成为三维泊松方程求解的基础，后来也被应用于亥姆霍兹方程和电磁学。Green本人是一位自学的面包师之子。
三维定常热传导的数值计算方法
FEM单元的选择
三维热分析应该用哪种单元？
三维热传导单元的对比。
| 单元类型 | 节点数 | 温度分布 | 精度 | 网格生成 |
|---|---|---|---|---|
| 4节点四面体 (TET4) | 4 | 线性 | 低 | 自动网格容易 |
| 10节点四面体 (TET10) | 10 | 二次 | 高 | 自动网格容易 |
| 8节点六面体 (HEX8) | 8 | 三线性 | 中～高 | 需要结构化网格 |
| 20节点六面体 (HEX20) | 20 | 二次 | 非常高 | 生成困难 |
在Ansys中，SOLID70（HEX8）、SOLID87（TET10）、SOLID90（HEX20）是代表性的热单元。在Abaqus中对应DC3D4、DC3D10、DC3D8、DC3D20。
实务中TET10比较稳妥吧？
考虑到自动网格生成的兼容性，TET10的通用性较强。但如果能生成六面体网格的话，HEX8可以减少单元数。Ansys Meshing的Multizone法或Sweep法可以优先生成六面体网格。
大规模问题的处理
三维问题单元数会剧增。当超过100万单元时的对策：
子模型技术具体怎么做？
先求解粗网格全局模型的定常解，然后在关注区域的边界面上转移温度，再用细网格对局部模型求解。Ansys Workbench中有Submodel命令可以简化设置。即使只对全体的10%的区域细分10倍，计算成本也只有全局细化的约1/100。
FEM四面体单元的诞生
热分析中不可或缺的四面体单元由Turner等在1960年代受美国航空航天局委托进行研究时提出。最初用于结构分析，热分析的应用整备则是在1970年代由Wilson和Bathe完成的。如今在Abaqus/Standard和Ansys Mechanical中使用的C3D10（10节点二次四面体）的起源也可以追溯到这个时代。
三维定常热传导的实务应用
形状简化的指导原则
如果直接用CAD数据的话单元数会爆炸吧？
CAD形状清理是三维热分析成败的关键。
| 简化项目 | 效果 | 注意 |
|---|---|---|
| 移除微小圆角 | 单元数减少50% | 0.5mm以下为标准 |
| 薄壁部分壳化 | 不需要厚度方向单元 | 厚度方向梯度很重要的情况不适用 |
| 省略螺栓孔 | 避免局部细化 | 仅限于不影响热路径的孔 |
| 利用对称性 | 1/2～1/8模型 | 边界条件也必须对称 |
用SpaceClaim（Ansys）或Design Modeler简化吗？
就是。Ansys Discovery Live可以边改形状边实时看到温度场，可以快速判断哪些形状特征对温度场有影响。
网格划分策略
三维定常热传导的网格划分方针：
- 发热源附近：最小单元尺寸 = 发热源尺寸的1/5以下
- 远离区域：可以粗化（温度梯度小）
- 材料界面：沿界面网格划分（节点共享或Tied Contact）
- 薄壁部分：厚度方向至少3层
收敛验证用什么标准？
如果关注的是最高温度，将网格加密一倍后最高温度变化在1%以内就足够了。但局部热流密度收敛较慢，需要进行3阶段以上的网格验证。
汽车发动机缸盖的热分析
丰田在2020年GR雅力士开发中利用三维定常热分析，优化了铝铸造缸盖燃烧室周围的温度分布。相比之前的设计，燃油喷嘴附近的最高温度降低了约25℃，同时实现了爆震耐性和热效率的平衡。分析模型的节点数超过200万，8核PC的计算时间约4小时。
三维定常热传导的软件比较
工具对比
进行三维定常热传导分析时，用哪个工具比较好？
按用途选择很重要。
| 用途 | 推荐工具 | 理由 |
|---|---|---|
| 结构零件温度场 | Ansys Mechanical、Abaqus | 热-结构耦合容易 |
| 电子设备 | Ansys Icepak、FloTHERM、6SigmaET | 零部件数据库完整 |
| 多物理场 | COMSOL | 电磁热耦合等优势强 |
| 研究与教育 | COMSOL、OpenFOAM | 可定制性高 |
FloTHERM和Icepak有什么区别？
FloTHERM由Mentor Graphics（现Siemens）开发，采用结构化网格。SmartPart（自动参数化零部件）是其强项。Icepak是Ansys公司出品，采用非结构化网格，基于Fluent求解器。近来Icepak被集成到Ansys Electronics Desktop中，与电磁场解析（HFSS、Maxwell）的耦合变得容易了。
计算规模的参考值
| 模型规模 | 单元数 | 计算时间(16核) | 内存 |
|---|---|---|---|
| 小型零部件 | ～10万 | ～1分钟 | ～2GB |
| 电路板组件 | 10万～100万 | 1～10分钟 | 2～16GB |
| 机箱整体 | 100万～1000万 | 10分钟～1小时 | 16～64GB |
定常热传导是线性的所以计算量比较小吧。
没错。同样的模型做非定常分析的话计算量会按时间步数倍增。能用定常求解就用定常，这是铁则。
Siemens Star-CCM+的热分析功能
Siemens Digital Industries的Star-CCM+前身是1999年CDadapco开发的STAR-CD。2015年强化了在同一求解器中进行三维热传导和流体耦合的"FV-FEM共享网格"功能，随后在电动车电池包三维热管理分析中的应用激增。2017年被Siemens收购。
三维定常热传导的先进研究
3D拓扑优化
三维热传导的最新趋势是什么？
三维热传导拓扑优化受到关注。以"温度场均匀化"或"最大温度最小化"为目标函数，优化材料分布。
用SIMP法（Solid Isotropic Material with Penalization）设 $k(\rho)=\rho^p k_0$。Ansys Mechanical和COMSOL都有实现。
优化后的形状能制造出来吗？
金属3D打印（SLM/EBM）的进步使得复杂的内部流道和格子结构成为可制造的。GE航空发动机喷嘴就是典型例子。
均质化方法的多尺度分析
要将微观结构（纤维排列、多孔性）的影响反映到三维宏观模型中，可以用均质化方法。用RVE（Representative Volume Element）求有效导热系数张量，应用于宏观模型。
Abaqus或COMSOL能做吗？
Abaqus可以通过Python脚本结合*ORIENTATION实现各向异性k在单元级的分配。COMSOL有Homogenization插件可以直接执行RVE分析。Digimat等专用工具也很有效。
机器学习代理模型
三维热传导的设计空间探索中，用PINN（物理信息神经网络）的研究在增加。生成数百个FEM学习数据，用神经网络瞬间预测温度场。
只用数百个学习数据就能用？
定常热传导是线性且解光滑，所以用相对较少的学习数据就能建立高精度代理。DeepXDE（PINN框架）和NVIDIA Modulus都可用。
AMD 3D V-Cache的热问题
AMD Ryzen 7 5800X3D（2022年）采用3D堆叠缓存技术"3D V-Cache"，但堆叠的SRAM芯片阻碍了CPU核心的散热，最高温度达到106℃。AMD通过三维定常热分析优化了堆叠芯片间的TIM（热界面材料）厚度和导热系数，并实施了固件对策以最小化热节流。
三维定常热传导的故障排除
常见问题
三维热分析有什么常见问题？
1. 接触面的节点不匹配
在组件模型中，两个零件接触面的网格不对齐时，如果Ansys Mechanical中Contact/Target对未正确设定，接触面就会出现温度不连续。
对策：启用Bonded Contact的Thermal Contact。通常将接触导热系数设定为10000 W/(m2K)以上，以实现实质的完全接触；或输入实测值。
2. 材料属性的单位转换
用mm单位系计算时有什么注意？
Ansys Mechanical中mm单位系的情况：
| 物理量 | SI单位 | mm单位系 |
|---|---|---|
| 长度 | m | mm |
| 温度 | K (或 degC) | K (或 degC) |
| 导热系数 | W/(m K) | mW/(mm K) → 数值相同 |
| 传热系数 | W/(m2 K) | mW/(mm2 K) → ×10^-6 |
| 发热密度 | W/m3 | mW/mm3 → ×10^-9 |
特别是传热系数和发热密度的换算容易出错。h=10 W/(m2K)在mm单位系中是10e-5 mW/(mm2K) = 0.00001。弄错这里温度会偏3个数量级。
3. 网格不足导致的温度场不自然
温度云图中出现等温线锯齿状时，说明单元太粗。TET4单元特别容易在粗网格下出现阶梯状温度分布。改用TET10或加细网格。
TET4真的不行啊。
利用模型对称性减少计算成本
当三维热分析模型收敛困难时，首先尝试利用对称面缩小到1/2或1/4模型。某涡轮机制造商使用周期对称边界条件（Abaqus的*CYCLIC SYMMETRY命令）将从360度模型的72小时计算时间缩减到单个扇形段的2小时，实现了36倍的加速。
更多
错误