基板级热解析
基板级热的理论基础
概述
老师！今天是关于基板级热解析的，对吧？那是什么？
多层PCB的各向异性热导模型。铜布线图案引起的有效热导率评估。部件布置优化的应用。
老师的解释很清楚！各向异性导热模型的困惑都消散了。
管制方程
离散化方法
计算机如何实际求解这个方程？
使用有限元法（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 |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
供应商系譜与产品整合历程
各个软件的来历，是不是有些曲折的故事？
Ansys Mechanical (旧ANSYS Structural)
请告诉我关于"Ansys Mechanical"的事！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys参数化设计语言）。
现在所属：Ansys Inc.
COMSOL Multiphysics
请告诉我关于"COMSOL Multiphysics"的事！
1986年在瑞典成立。以MATLAB联动的FEMLAB起步，后改名为COMSOL。多物理领域强项。
现在所属：COMSOL AB
老师讲得很清楚！工具名称的困惑消散了！
Ansys Fluent
接下来是Ansys Fluent的话题，对吧。内容是什么？
由Fluent Inc.开发。2006年被Ansys收购。非结构网格式通用CFD求解器。
现在所属：Ansys Inc.
哦~，开发历史的话，真是太有趣了！想听更多。
文件格式与互操作性
在不同软件间传递数据时有什么注意点？
| 格式 | 扩展名 | 种类 | 概述 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | ISO 10303兼容的3D CAD数据交换格式。支持形状+PMI。 |
| IGES | .igs/.iges | 中立CAD | 初期CAD数据交换规范。曲面数据兼容性有问题。正向STEP迁移。 |
| VTK | .vtk/.vtu | 可视化 | 可视化工具包格式。ParaView等使用。 |
不同求解器间转换模型时，需注意单元类型的对应关系、材料模型的兼容性、荷载边界条件的表示差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）在求解器间通常无法直接转换。
原来格式看似简单，实际上非常深邃啊。
实际应用中的注意事项
教科书里没有的"现场智慧"之类的有吗？
网格收敛性确认、边界条件妥当性验证、材料参数敏感性分析非常重要。
基板级热解析的全貌已理解！明天起在实际工作中加以重视。
不错！实际动手才是最好的学习。有不懂的随时问我。
PCB的热导率各向异性很大
印刷电路板（PCB）的热导率在面内方向（XY）为20~30 W/m·K，厚度方向（Z）为0.3~0.5 W/m·K，各向异性达约100倍。这源于铜布线层（热导率385 W/m·K）和玻璃环氧树脂（FR-4、热导率0.25 W/m·K）的层积结构。热解析中广泛采用Delphi法（1995年）从层比率计算等效热导率。
基板级热的数值计算方法
数值方法详情
具体用什么算法求解基板级热解析？
等等，基板级热解析对来说，是说这样的情况也能用？
离散化的制定
用形状函数 $N_i$ 近似未知量：
用如下方式表示。
基础方程的离散形式
用如下方式表示。
嗯~，光看公式有点不清楚… 代表什么呢？
连续体的管制方程离散化后，得到下列代数方程组：
这里$[K]$是整体刚性矩阵（或等效系统矩阵），$\{u\}$是未知节点变量向量，$\{F\}$是外力向量。
哦，是那回事！连续体的管制方程就是这样离散的啊。
单元技术
听过"单元技术"这个词，但可能没彻底理解…
| 单元类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二阶 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二阶 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二阶 | 6/15 | 中~高 | 中 |
积分方案
积分方案具体是什么意思？
听完了，总算明白单元类型为什么这么重要！
收敛性和稳定性
不收敛了首先检查什么？
收敛速度：二阶单元误差以$O(h^2)$的阶数减小（平滑解的情况）
原来网格细化看似简单，实际上非常深邃啊。
求解器设置建议
具体用什么算法求解基板级热解析？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理方法 | ILU(0) 或 AMG | 根据问题规模 |
| 最大迭代次数 | 1000 | 未收敛则重新设置 |
| 内存模式 | In-core | 尽可能 |
线性单元 vs 2阶单元
热传导解析中线性单元通常也能获得足够精度。温度梯度急变区域（热冲击等）建议使用2阶单元。
热流评估
从单元内温度梯度计算。与节点应力相同，可能需要平滑处理。
对流-扩散问题
Pe数高（对流主导）时需要迎风稳定化（SUPG等）。纯热传导问题无需。
非定常解析的时间步长
热扩散特征时间$\tau = L^2 / \alpha$（$\alpha$：热扩散率）相比设置足够小的步长。急剧温度变化可用自动时间步长控制。
非线性收敛
温度相关物性值的非线性通常较弱，Picard迭代（直接代入法）多数情况足够。放射的强非线性推荐牛顿法。
定常解析的判定
全节点温度变化小于阈值（$|\Delta T| / T_{max} < 10^{-5}$等）时判定收敛。
基板级热的实际应用
实践指南
老师，请教我"实践指南"！
讲解基板级热解析的实际应用和注意点。
也就是说，基板级热解析的实际处理偷工减料会后患无穷。记住了！
分析流程
从一开始教我！该从哪里开始？
1. 预处理 (Pre-processing)
- CAD数据导入和形状简化
- 材料特性定义
- 网格生成（单元类型·尺寸决定）
- 边界条件和荷载条件的设置
2. 求解 (Solving)
- 求解器设置（解法、收敛基准、输出控制）
- 作业投入和计算执行
- 收敛监控
3. 后处理 (Post-processing)
- 结果可视化（位移、应力、其他物理量）
- 结果验证和妥当性确认
- 报告制作
网格生成最佳实践
网格的好坏怎样判断？
单元品质指标
请教我"单元品质指标"！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 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 |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
Ansys Mechanical (旧ANSYS Structural)
请告诉我关于"Ansys Mechanical"的事！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys参数化设计语言）。
现在所属：Ansys Inc.
COMSOL Multiphysics
请告诉我关于"COMSOL Multiphysics"的事！
1986年在瑞典成立。以MATLAB联动的FEMLAB起步，后改名为COMSOL。多物理领域强项。
现在所属：COMSOL AB
老师讲得很清楚！工具名称的困惑消散了。
常见故障及对策
初学者容易犯的错误有吗？事先知道很重要！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、不恰当的边界条件 | 改进网格、重新检查约束条件 |
| 应力异常大 | 应力奇点、网格依赖 | 避免奇点、局部网格细化 |
| 位移非现实 | 材料常数误差、单位系不一致 | 确认输入数据 |
| 计算时间过长 | 不必要的细化、低效求解法 | 网格优化、并行计算 |
质量保证检查清单
教科书里没有的"现场智慧"之类的有吗？
基板级热解析的全貌已理解！明天起在实际工作中加以重视。
不错！实际动手才是最好的学习。有不懂的随时问我。
GPU板的铜铺面积是关键
搭载NVIDIA A100的主板上，通过将GPU芯片正下方的PCB铜铺面积最大化，可将主板内的热扩散抵抗降低40%，这是2021年ITHERM会议报告的设计事例。在实际应用中，制作整个主板的"铜充填率地图"，同时评估热学和电学（接地平面）的平衡的做法已成定式。
基板级热的软件比较
商用工具比较
各种软件都有呢？各自的特点教我！
讲解基板级热解析对应的主要商用CAE工具的功能比较，及各产品的历史背景。
等等，基板级热解析对来说，是说这样的情况也能用？
支持工具列表
那么，做基板级热解析需要什么样的软件？
| 工具名称 | 开发者/现所属 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
Ansys Mechanical (旧ANSYS Structural)
请告诉我关于"Ansys Mechanical"的事！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys参数化设计语言）。
现在所属：Ansys Inc.
COMSOL Multiphysics
请告诉我关于"COMSOL Multiphysics"的事！
1986年在瑞典成立。以MATLAB联动的FEMLAB起步，后改名为COMSOL。多物理领域强项。
现在所属：COMSOL AB
从这里听下去，总算明白开发历史为什么重要！
Ansys Fluent
接下来是Ansys Fluent的话题，对吧。内容是什么？
由Fluent Inc.开发。2006年被Ansys收购。非结构网格式通用CFD求解器。
现在所属：Ansys Inc.
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题，对吧。内容是什么？
由CD-adapco开发。2016年被西门子收购，并入Simcenter品牌。多面体网格是特色。
现在所属：Siemens Digital Industries Software
哦，是那回事！开发历史就是这么个意思。
功能比较矩阵
预算和时间都有限，性价比最强的是哪个？
| 功能 | Ansys Mechanical | COMSOL | Fluent | Star-CCM+ |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换时的风险
转换时的风险具体是什么意思？
哦，是那回事！不同工具间的模传递就是这么个意思。
许可证形式
听过"许可证形式"这个词，但可能没彻底理解…
| 工具 | 许可 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高额但有官方支持 |
| OpenFOAM | GPL | 免费但支持收费 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的OSS求解器 |
选择指南
最后到底选哪个，判断标准教我？
基板级热解析的工具选择应考虑以下因素：
基板级热解析的全貌已理解！明天起在实际工作中加以重视。
不错！实际动手才是最好的学习。有不懂的随时问我。
Mentor FloTHERM XT占据业界首位
在基板级热解析领域，Mentor Graphics（现Siemens EDA）的FloTHERM XT多年来保持市场首位，ECAD数据的直接导入和庞大的元器件库（超50,000点）是其优势。竞争对手ANSYS Icepak在2019年加强了流体-结构耦合，在大型数据中心用户中普及。国内方面，Software Cradle的Thermofluid在汽车电装品制造商中有较多导入实绩。
基板级热的前沿研究
前沿话题与研究趋势
基板级热解析领域今后怎样演进？
看基板级热解析最新研究动向和先进手法。
最新的数值手法
接下来是最新数值手法的话题，对吧。内容是什么？
嗯~，光看公式有点不清楚… 代表什么呢？
高性能计算 (HPC) 的应对
| 并行化方法 | 概述 | 应用求解器 |
|---|---|---|
| MPI (领域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU利用。特别对显式法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
基板级热的故障排除
故障排除
也就是说，基板级热解析能做好的话，差不多就可以了？
常见错误及对策
老师也在基板级热解析里通宵调试过吗？（笑）
1. 收敛失败
收敛失败具体是什么意思？
症状：求解器在指定迭代次数内无法收敛，异常终止
可能原因：
- 网格品质不足（过度扭曲的单元）
- 材料参数设置不当
- 初始条件不适当
- 非线性太强（荷载步不足）
对策：
- 进行网格品质检查（纵横比、雅可比）
- 确认材料参数的单位系
- 将荷载分为多个步骤（增加子步数）
- 放松收敛判定基准（但注意精度）
也就是说，收敛失败这里偷工减料会后患无穷。记住了！
2. 非物理的结果
接下来是非物理的结果的话题，对吧。内容是什么？
症状：应力/位移/温度等出现物理上非现实的值
可能原因：
- 边界条件设置错误
- 单位系混淆（SI单位与工程单位）
- 不恰当的单元类型选择
- 应力奇点存在
对策：
- 确认反力的总合（力的平衡）
- 确认单位系的一致性
- 重新检讨单元类型的适切性
- 奇点消除或子模型化
前辈说"收敛失败一定要好好处理"的意思明白了。
3. 计算时间超过
计算时间超过具体是什么意思？
症状：计算耗时远超预期
对策：