2次元定常热传导
2次元定常热传导的理论基础
2次元热传导的基础
2次元的定常热传导是如何从1次元扩展的？
温度在2个方向变化的问题。电子基板面内温度分布、金型断面温度场、建筑墙体热桥评估等是代表例。
控制方程
等向性材料的2次元定常热传导方程如下。
若 $k$ 为常数且无发热，则化为Laplace方程 $\nabla^2 T = 0$。
有解析解吗？
对于矩形区域，各边指定温度的情况下，用变量分离法可得到级数解。例如3边为0℃、1边为 $T_0$ 的正方形区域
这个解对于FEM代码的验证基准非常重要。
形状系数法
2次元的几何学效果用形状系数 $S$ 以1次元方式处理。
例如，地中埋设管的热损失可用 $S = 2\pi L / \cosh^{-1}(D/r)$ 计算。教科书（Incropera等）中总结了主要形状的 $S$ 值。
用形状系数就能手计算2D问题了吧。
是的。但对于复杂形状，形状系数未知，所以用FEM求解。从FEM结果逆推形状系数的值，在类似设计中重新使用也是实务中经常做的。
Laplace方程与解析解的优美性
2次元定常热传导由∇²T=0（Laplace方程）描述。1822年Fourier给出的矩形板解析解用sin级数表示，至今仍被用作验证基准。Laplace方程与电势、流速势具有相同形式，19世纪物理学家称赞这种统一性为"自然深层的优美"。
2次元定常热传导的数值计算方法
FEM的2次元离散化
2D问题的FEM用什么单元？
2D热传导有三角形单元和四边形单元。
| 单元 | 节点数 | 温度分布 | 精度 |
|---|---|---|---|
| 3节点三角形 | 3 | 线性 | 低（温度梯度为常数） |
| 6节点三角形 | 6 | 二次 | 高 |
| 4节点四边形 | 4 | 双线性 | 中 |
| 8节点四边形 | 8 | 二次 | 高 |
3节点三角形内温度梯度为常数，需要细分才能得到精度。实务推荐6节点三角形或8节点四边形。
四边形和三角形哪个更好？
节点数相同时四边形精度更高。但复杂形状的自动网格生成中，三角形（四面体）更容易生成。Ansys的网格划分器选择六面体优先（hex-dominant）可以取得平衡。
有限差分法（2次元）
等间距网格的中心差分形式为
当 $\Delta x = \Delta y$ 时变成著名的5点模板。可以用Excel或自编Python代码验证。
用自编代码验证商用求解器的结果，这个方法很踏实呢。
特别是2D验证只需几十行代码，用Excel VBA或Python就能写，没有不做的理由。
FDM中心差分的二阶精度
2D有限差分法的中心差分格式具有空间二阶精度，网格宽度减半时误差变为1/4。1960年代IBM 704上NASA进行的火箭喷管2D热分析在128×128网格上也需要数小时，而现代PC可在0.01秒内完成相同计算。50年间运算速度提高了10⁹倍以上。
2次元定常热传导的实务应用
2D模型的活用场景
3D分析已成主流，为什么还要用2D模型？
计算成本差几个数量级。在进行数十个案例的参数化研究或断面方向温度梯度占主导的问题中，2D最优。
实务例：PCB基板的热桥分析
多层PCB的断面建模时，各层（铜箔、预浸料、芯板）用不同的 $k$ 值作为平行带状排列。铜箔率80%的L1层与20%的L3层等效热导率差异很大。
| 层 | 厚度 [um] | 面内k [W/(mK)] | 面直k [W/(mK)] |
|---|---|---|---|
| L1 (铜80%) | 35 | 318 | 1.1 |
| 预浸料 | 100 | 0.3 | 0.3 |
| L2 (铜50%) | 35 | 199 | 0.55 |
| 芯板 | 800 | 0.3 | 0.3 |
面内和面直相差2位数以上呢。
PCB基板是极端的异向性材料。通过2D断面分析评估面直方向的温度梯度和过孔效应，面内方向用扩散电阻补正。这也是FloTHERM和Icepak内部建模使用的手法。
对称性的活用
利用对称性可大幅减少计算规模。
- 1轴对称: 1/2模型，对称面为隔热条件
- 2轴对称: 1/4模型
- 周期对称: 仅建模1个重复单位
对称面为隔热条件的原因？
对称面两侧温度场呈镜像对称时，横跨对称面的温度梯度为零，即热流量为零。从物理角度看是显而易见的推论。
CPU芯片的2D温度图
Intel在2003年的Pentium 4"Prescott"核心（90 nm工艺，TDP 84 W）上大量使用2D定常热分析进行热点管理。芯片上的最大温差超过20℃，通过2D图谱识别电流密度集中区域并重新设计布线层。这种手法后来成为多核设计的标准。
2次元定常热传导的软件比较
各工具的2D分析
各个软件怎么进行2D热传导分析？
主要工具的2D热解析单元一览。
| 工具 | 单元 | 备注 |
|---|---|---|
| Ansys Mechanical | PLANE55(4节点), PLANE77(8节点) | 选择平面/轴对称 |
| Abaqus | DC2D4(四边形), DC2D3(三角形) | DCAX4为轴对称 |
| COMSOL | 2D component → Heat Transfer | 自动网格化任意形状 |
| Ansys Fluent | 2D/Axi-symmetric mesh | 基于有限体积法 |
轴对称模型用2D表现3D，很方便呢。
完全同意。圆柱形散热器或旋转体的温度场用轴对称2D足够精度。多数情况下计算成本仅为3D的1/100以下，却能得到相当结果。
COMSOL设置示例
COMSOL进行2D定常热传导的步骤:
1. Model Wizard → 2D → Heat Transfer in Solids
2. Geometry: 用Rectangle等创建区域
3. Materials: 内置材料或用户自定义
4. Heat Transfer: 设置Thermal insulation(默认)、Temperature、Heat Flux等
5. Mesh → Physics-controlled (Normal)
6. Study → Stationary
COMSOL的优势是在GUI上可直接编辑弱形式，自定义方程添加容易。
COMSOL适合研究用途呢。
伴随实验验证的研究中非常易用。另一方面，要自动运行大量设计案例时，Ansys Workbench的Design Explorer工作流更成体系。
COMSOL Multiphysics的2D热模块
COMSOL在2005年从FEMLAB更名，大幅改进了2D热传导分析的GUI。特别是参数化扫描功能使得材料λ可在一行设置中改变100通，复合材料最优设计研究者数量激增。目前全球180多个国家在使用。
2次元定常热传导的先端研究
等价异向性模型
2D的先端话题有哪些？
将多层PCB和层合复合材料作为等价异向性连续体进行2D建模的手法在实务中很重要。面内和面直的有效热导率为
面内用算术平均（并联），面直用调和平均（串联）计算。
XFEM/GFEM的应用
含有材料界面或裂纹的2D热传导问题中应用扩展有限元法（XFEM）。即使网格不与界面对齐，也能表现界面处的温度跳跃。
不用与界面对齐很轻松呢。
特别在裂纹或剥离扩展的问题上有效。Abaqus中用*XFEM关键字可使用，但热分析的应用目前还处于研究阶段。
AI网格优化
最近有用强化学习自动化2D热传导问题网格自适应细分的研究。自动检测温度梯度大的区域，用最少单元数达到目标精度。
AI给我们划网格的时代来了呀。
Ansys 2025R1搭载了Mesh Morpher AI功能，实现了边界层网格自动优化。2D问题计算成本低，是AI网格优化的理想测试环境。
等温线与流线的直交性
与2D势流相同，2D定常热传导的等温线与热流束线始终直交。利用这一性质的"流束图法"由Schlichting在1940年代推广，在FEM普及前的1970年代，工学部学生们用手绘此图来推算复杂形状的热流束，是当时的实用技术。
2次元定常热传导的故障排除
常见问题
请教2D热传导容易踩的坑。
1. 平面模型与轴对称模型的混淆
问题: Ansys PLANE55应指定轴对称却设为平面时，圆柱形状的热阻会完全不同。平面问题按单位厚度计算，轴对称考虑2πr。
对策: 必须检查KEY OPT设置。Abaqus因要素名不同（DC2D4 vs DCAX4），容易避免出错。
2. 异向性材料的坐标系
异向性材料坐标系错误会怎样？
材料坐标系与全局坐标系不一致时，面内k和面直k会互换。PCB基板中 $k_{in}=30$、$k_{through}=0.5$ W/(m K) 互换会导致温度上升数十倍变化。Ansys Workbench中必须逐要素检查Coordinate System的分配。
3. 直角拐角特异点
直角拐角的内角部分理论上温度梯度（热流束）趋于无穷大。网格越细，数值越大，不收敛。
对策: 不用局部最大热流值，用一定区域平均值进行设计。或加微小圆角消除特异性。
与结构分析的应力特异点一样的处理方法呢。
完全一样。实物中温度特异点会被圆角或倒角缓解，所以分析模型中加微小R（约0.1 mm）是现实的做法。
角部特异性与网格依赖性
2D形状的锐重入角处理论上热流束无穷大，出现特异解。对FEM中该部分细分，温度梯度越来越大，不收敛。Abaqus手册（2018版）推荐对重入角附近应用特殊要素（四分点要素）。
价值
更详细
错误