热过孔
热过孔的理论基础
概述
老师！ 今天是讲热过孔吧？这是什么东西呢？
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 |
| ANSYS Fluent | ANSYS Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | 西门子数字工业软件 | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
供应商系谱与产品整合历史
这些软件的发展历史，是不是也有戏剧性的故事呢？
ANSYS Mechanical (旧ANSYS Structural)
请给我讲讲"ANSYS Mechanical"！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys参数化设计语言）。
现属于：ANSYS Inc.
ANSYS Fluent
下面讲ANSYS Fluent是吧。是什么内容呢？
由Fluent Inc.开发。2006年被ANSYS收购。非结构网格的通用CFD求解器。
现属于：ANSYS Inc.
经过这样的讲解，我终于理解为什么说"一定要学好开发历史"了！
Simcenter STAR-CCM+
下面讲Simcenter STAR是吧。是什么内容呢？
由CD-adapco开发。2016年被西门子收购，纳入Simcenter品牌。多面体网格是其特点。
现属于：西门子数字工业软件
哇，开发历史的讲解超有趣！ 还想再听听呢。
文件格式与互操作性
在不同软件之间交换数据时有什么需要注意的吗？
| 格式 | 扩展名 | 类型 | 概述 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | 符合ISO 10303的3D CAD数据交换格式。形状+PMI支持。 |
| IGES | .igs/.iges | 中立CAD | 早期的CAD数据交换规范。曲面数据兼容性有问题。逐渐转向STEP。 |
在不同求解器间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）往往无法在求解器间直接转换。
原来如此啊…文件格式看似简单，其实内涵非常深。
实务注意事项
有没有教科书里没讲的"现场的智慧"呢？
网格收敛性的验证、边界条件的合理性检证、材料参数的敏感性分析都特别重要。
哇，热过孔竟然这么深啊…但经过老师的讲解，我渐渐理清了思路！
嗯，进展不错啊！ 亲身动手实践才是最好的学习。有不懂的随时问我。
热过孔的导热率因充填材料而变化10倍
PCB的热过孔（用于导热的贯通孔）中空（空气充填）时实效导热率约为0.1 W/m·K，但铜镀层充填（实心过孔）可达385 W/m·K。内径250μm、铜镀层厚25μm的中空过孔与同规格实心过孔的热阻差约为15倍，考虑到成本因素，通常采用环氧银膏充填（导热率3～8 W/m·K）作为折中方案。IPC-7093规范（2011年）已成为热过孔设计的业界标准。
热过孔的数值计算方法
数值方法详述
具体用什么算法求解热过孔呢？
也就是说在有限元法这里下功夫不够的话，后面会吃大亏。我会铭记在心！
离散化的定式化
用形状函数 $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迭代（直接代入法）已足够。放射的强非线性则推荐牛顿法。
定常分析的判定
以全节点温度变化小于阈值（$|\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 |
| ANSYS Fluent | ANSYS Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | 西门子数字工业软件 | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
ANSYS Mechanical (旧ANSYS Structural)
请给我讲讲"ANSYS Mechanical"！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys参数化设计语言）。
现属于：ANSYS Inc.
ANSYS Fluent
下面讲ANSYS Fluent是吧。是什么内容呢？
由Fluent Inc.开发。2006年被ANSYS收购。非结构网格的通用CFD求解器。
现属于：ANSYS Inc.
老师的讲解非常好懂！ 软件名称的困惑终于解开了。
常见失误及对策
初学者容易犯什么错误呢？ 提前知道比较好！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、边界条件不适当 | 网格改善、约束条件审视 |
| 应力异常偏大 | 应力奇点、网格相关 | 奇点回避、局部网格细化 |
| 位移非现实 | 材料常数错误、单位不一致 | 确认输入数据 |
| 计算时间超长 | 不必要的细化、求解效率低 | 网格最优化、并行计算 |
品质保证检查清单
有没有教科书里没讲的"现场的智慧"呢？
哇，热过孔竟然这么深啊…但经过老师的讲解，我渐渐理清了思路！
嗯，进展不错啊！ 亲身动手实践才是最好的学习。有不懂的随时问我。
LED基板标准设计中通常有100个过孔
搭载高功率LED（5～10W/芯片）的MCPCB上，LED焊盘正下方通常配置径0.3mm的热过孔，密度为100～200个/cm²，这是标准设计。日亚化学照明用LED应用笔记（2019年版）的测试数据显示，将过孔密度从64个/cm²增加到144个/cm²时Rth可降低40%。由于孔数成本按孔数比例计费，如何以最少孔数实现目标Rth的最优化是设计的核心。
热过孔的软件比较
商用工具比较
有那么多不同的软件呢？ 各自的特点都给我讲讲！
讲解支持热过孔分析的主要商用CAE工具的功能比较以及各产品的历史背景。
原来如此啊…支持热过孔分析看似简单，其实内涵非常深。
支持工具列表
那么做热过孔分析的话都有什么软件可以用呢？
| 工具名称 | 开发商/现属 | 主要文件格式 |
|---|---|---|
| ANSYS Mechanical (旧ANSYS Structural) | ANSYS Inc. | .cdb, .rst, .db, .ans, .mac |
| ANSYS Fluent | ANSYS Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | 西门子数字工业软件 | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
ANSYS Mechanical (旧ANSYS Structural)
请给我讲讲"ANSYS Mechanical"！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys参数化设计语言）。
现属于：ANSYS Inc.
ANSYS Fluent
下面讲ANSYS Fluent是吧。是什么内容呢？
由Fluent Inc.开发。2006年被ANSYS收购。非结构网格的通用CFD求解器。
现属于：ANSYS Inc.
经过这样的讲解，我终于理解为什么说"一定要学好开发历史"了！
Simcenter STAR-CCM+
下面讲Simcenter STAR是吧。是什么内容呢？
由CD-adapco开发。2016年被西门子收购，纳入Simcenter品牌。多面体网格是其特点。
现属于：西门子数字工业软件
COMSOL Multiphysics
请给我讲讲"COMSOL Multiphysics"！
1986年在瑞典成立。始于与MATLAB连携的FEMLAB，后改名为COMSOL。多物理场解析见长。
现属于：COMSOL AB
啊，原来如此！ 开发历史就是这样的机制。
功能比较矩阵
预算和时间都有限，性价比最高的是哪个呢？
| 功能 | ANSYS Mechanical | Fluent | STAR-CCM+ | COMSOL |
|---|---|---|---|---|
| 基础功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换时的风险
转换时的风险，具体是指什么呢？
啊，原来如此！ 不同工具间的模型转换就是这样的机制。
许可形式
听说过"许可形式"，但还没真正理解…
| 工具 | 许可 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持有偿 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的OSS求解器 |
选型指南
最后到底选哪个，有判断标准吗？
热过孔工具选型时应考虑以下因素：
哇，热过孔竟然这么深啊…但经过老师的讲解，我渐渐理清了思路！
嗯，进展不错啊！ 亲身动手实践才是最好的学习。有不懂的随时问我。
伊比德和新光电气在热过孔实装中处于国内首位
包含热过孔在内的高密度多层基板国内最大制造商是伊比德（岐阜县）和新光电气工业（长野县），两公司在Intel CPU封装基板全球市场中占有一定份额。伊比德的"FCBGA"封装基板截至2023年广泛应用于AI GPU封装。在TSV制造方面，台湾的ASE（日月光集团）作为最大的外包代工企业，而Samsung铸造厂、TSMC和UMC则采用自主一体化制造。
热过孔的先进研究
先进话题与研究动向
热过孔这一领域今后会如何发展呢？
让我们来看看热过孔领域最新的研究动向和先进手法。
前辈说"热过孔的先进手法一定要做好"，现在我理解它的意思了。
最新的数值方法
下面讲最新数值方法是吧。是什么内容呢？
只看公式我不太明白…这是表示什么呢？
高性能计算 (HPC) 的适配
| 并行化手法 | 概述 | 支持求解器 |
|---|---|---|
| MPI (领域分割) | 分布式内存型。大规模问题的标准方案 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU利用。特别在显式格式中有效 | LS-DYNA、Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
热过孔的故障排除
故障排除
也就是说在热过孔分析相关的地方不够用心的话，后面会吃大亏。我会铭记在心！
常见错误及对策
老师也曾经为热过孔彻夜调试过吗？（笑）
1. 收敛失败
收敛失败，具体是指什么呢？
症状：求解器在指定迭代回数内未收敛而异常结束
可能原因：
- 网格品质不足（过度扭曲的单元）
- 材料参数设置不当
- 初始条件不适当
- 非线性过强（缺乏足够的荷载步）
对策：
- 执行网格质量检查（纵横比、雅可比）
- 确认材料参数的单位系统
- 将荷载分成多个步骤（增加子步数）
- 放松收敛判定基准（但注意精度）
也就是说在收敛失败这里不够用心的话，后面会吃大亏。我会铭记在心！
2. 非物理的结果
下面讲非物理的结果是吧。是什么内容呢？
症状：应力/位移/温度等出现物理上不现实的数值
可能原因：
- 边界条件设置错误
- 单位系混杂（SI单位与工学单位混用）
- 单元类型选择不当
- 应力奇点存在
对策：
- 确认反力合计（力的平衡）
- 确认单位系的一致性
- 重新审视单元类型的适切性
- 奇点消除或局部子模型
前辈说"收敛失败一定要做好"，现在我理解它的意思了。