2、热仿真基础：热传导、对流、辐射基本理论，热阻网络模型

在PA（功率放大器）版图热仿真中，一切结温预测与裕量复核都建立在传热学三大基本方式与热阻网络模型之上。PA芯片的功耗密度极高，热流从晶体管结区出发，依次穿过硅衬底、焊接层、热沉、PCB铜皮与散热器，最终通过空气或机壳带走。这一路径上的每一段材料、每一个界面都会贡献热阻，而热阻的串并联关系决定了结温的最终数值。本章将系统梳理热传导、对流、辐射的工程适用场景，并重点讲解如何将物理传热过程抽象为热阻网络模型，为后续的仿真参数设置提供理论依据。

上图展示了PA芯片从结区到环境的热流路径与对应的热阻网络。实际工程中，热阻网络并非简单的单链串联，而是包含并联支路（如PCB背面散热、侧面辐射）与多维扩展路径的复杂网络。理解这一模型是进行仿真参数设置的前提——仿真软件中的材料导热率、界面热阻、对流换热系数等参数，本质上都是在为这个网络中的每个节点赋值。

2.1 热传导：固体内部的主要传热机制

核心原理：热传导是固体内部或相互接触的固体之间，由于温度梯度引起的能量传递。其基本规律由傅里叶定律描述：

q = -k × A × (dT/dx)

其中，q为热流速率（W），k为材料导热率（W/m·K），A为垂直于热流方向的截面积（m²），dT/dx为温度梯度（K/m）。在PA版图仿真中，硅的导热率约为130~150 W/m·K（随温度升高而下降），铜约为385~400 W/m·K，FR4基材仅约0.3 W/m·K。这一数量级差异决定了热流优先沿铜皮和过孔扩散，而非穿过FR4层。

工程判断方法：在设置仿真参数时，需要关注以下三个关键点：

• 导热率的温度依赖性：硅的导热率在125°C时比25°C时下降约30%，仿真中应使用随温度变化的曲线，而非恒定值。
• 界面热阻：芯片与热沉之间的焊接层、导热硅脂层，其厚度和接触面积直接影响界面热阻。典型焊接层厚度为50~100μm，导热硅脂的导热率约1~5 W/m·K。
• 热扩散角：热流从芯片底部向下传播时，会以约45°的角度扩散。仿真中若使用一维热阻计算，会高估结温；应使用三维模型或设置合理的扩散热阻。

小型场景化说明：假设一个PA芯片功耗为10W，结到壳热阻Rth_jc=0.5°C/W，壳到热沉界面热阻Rth_cs=0.3°C/W，热沉到环境热阻Rth_sa=4°C/W。若环境温度为25°C，则结温Tj=25+10×(0.5+0.3+4)=73°C。若将界面热阻从0.3降至0.1°C/W，结温降至71°C，仅改善2°C；但若将热沉热阻从4降至2°C/W，结温降至53°C，改善20°C。这说明热阻网络中数值最大的环节是优化重点。

2.2 热对流：散热器与空气的热交换

核心原理：热对流是流体（空气或液体）流经固体表面时发生的热量传递。牛顿冷却定律描述为：

q = h × A × (Ts - Ta)

其中，h为对流换热系数（W/m²·K），A为换热面积，Ts为固体表面温度，Ta为流体温度。自然对流时h约为5~25 W/m²·K，强制风冷时可达20~100 W/m²·K，液冷则更高。PA功放通常采用自然对流或低风速强制风冷，h值一般在10~50 W/m²·K之间。

仿真参数设置要点：

1. 换热系数的取值：自然对流时，垂直表面的h约5~10，水平表面约10~15；强制风冷需根据风速计算，2m/s风速下h约20~30。仿真软件中可设置为恒定值，也可使用CFD模块自动计算。
2. 散热器效率：翅片散热器的实际换热面积远大于基板面积，但翅片根部温度高、顶部温度低，存在效率因子η（通常0.6~0.9）。仿真中应使用有效换热面积A_eff = η × A_total。
3. 辐射的耦合：在自然对流且表面温度较高（>60°C）时，辐射散热量可占总散热量的20%~30%。仿真中应同时开启辐射模型，否则会低估散热能力。

案例拆解：某5G基站PA模块，功耗25W，采用铝挤散热器（基板100mm×100mm，翅片高20mm，共10片）。初始仿真设置h=15 W/m²·K，计算结温为98°C，超出85°C的规格要求。分析发现，自然对流下该散热器的实际h约为8 W/m²·K，且未考虑辐射。重新设置h=8并开启辐射模型后，结温仿真值为92°C，仍超标。进一步优化：增加翅片高度至30mm、数量至14片，并将表面做黑色阳极氧化（辐射率0.9），最终仿真结温为83°C，满足要求。实测结果84.5°C，与仿真偏差1.8%，验证了参数设置的准确性。

2.3 热辐射：高温表面的不可忽略项

核心原理：热辐射是物体通过电磁波传递能量的方式，不需要介质。斯蒂芬-玻尔兹曼定律描述为：

q = ε × σ × A × (Ts⁴ - Ta⁴)

其中，ε为表面发射率（黑体为1，光亮铝约0.05，黑色阳极氧化约0.85），σ为斯蒂芬-玻尔兹曼常数（5.67×10⁻⁸ W/m²·K⁴）。由于温度以四次方出现，当表面温度与环境温差较大时，辐射散热显著增加。

工程判断方法：

• 当表面温度低于60°C时，辐射散热量通常小于总散热量的10%，可忽略。
• 当表面温度超过80°C且处于自然对流环境时，辐射必须纳入仿真。
• PA功放外壳若为金属喷漆或阳极氧化，ε取0.7~0.9；若为光亮铝，ε仅0.05~0.1，辐射散热几乎可忽略。

常见误区：许多工程师在仿真中忽略辐射，认为“辐射在低温下不重要”。但在PA自然对流散热场景中，散热器表面温度常在70~90°C，辐射散热量可达总散热量的25%以上。忽略辐射会导致结温仿真值偏高5~10°C，从而造成不必要的散热器尺寸加大或误判热设计不达标。

2.4 热阻网络模型：从物理到数学的抽象

核心原理：热阻网络模型将传热路径上的每个环节抽象为热阻，类似于电路中的电阻。热阻定义为温差与热流的比值：

Rth = ΔT / P

单位°C/W。对于一维稳态热传导，Rth = L/(k×A)，其中L为传热路径长度。热阻网络可以是串联、并联或串并联混合结构，通过电路类比法求解结温。

PA版图典型热阻网络结构：

建模方法：在PA版图热仿真中，热阻网络模型有两种应用方式：

1. 快速估算：使用上述典型值进行手算或Excel计算，在方案阶段快速评估结温裕量。精度约±20%，足以判断热设计方向。
2. 仿真校准：在CFD仿真或有限元仿真中，将仿真结果反推为等效热阻，与实测数据对比，校准模型参数。例如，实测结温85°C、功耗20W、环境25°C，则总热阻Rth_ja=(85-25)/20=3°C/W，可分解为各环节热阻。

案例拆解：某PA模块在版图设计阶段，工程师使用热阻网络模型进行预评估。芯片规格书给出Rth_jc=0.8°C/W，散热器供应商提供Rth_sa=3.5°C/W（自然对流），界面采用0.2mm导热硅脂（k=3 W/m·K，面积100mm²），计算Rth_cs=0.2/(3×0.0001)=667°C/W——这个数值明显异常。检查发现，导热硅脂的面积应为芯片面积（如5mm×5mm=25mm²），而非散热器面积。修正后Rth_cs=0.2/(3×0.000025)=2667°C/W，仍然很大。进一步分析，实际接触面积因表面平整度仅约60%有效，且硅脂厚度可能仅0.05mm。最终取Rth_cs=0.15°C/W，总热阻Rth_ja=0.8+0.15+3.5=4.45°C/W。在25°C环境、20W功耗下，结温Tj=25+20×4.45=114°C，超出105°C的规格。该案例说明，热阻网络模型中的每个参数都需要仔细核实，尤其是界面热阻的计算，面积和厚度的误判会导致数量级错误。

总结复盘：热阻网络模型的价值在于将复杂的传热问题转化为直观的电路类比，便于工程师快速定位热设计瓶颈。但在使用中需注意：① 热阻值随温度变化，高温下硅导热率下降，热阻增大；② 并联路径不可忽略，PCB横向散热在小型化PA中占比可达15%~30%；③ 瞬态热仿真需引入热容（Cth），此时网络模型扩展为RC网络，用于分析脉冲功率下的结温波动。