技术分享丨PCB电热仿真方法及实例分析

allanwei1984 · 发表于 2022-5-17 08:41

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

本帖最后由 Heaven_1 于 2022-5-17 19:30 编辑

如今电子产品及其相关产业正在追求小型化、集成化、高频高速率，电路板上的元器件功率越来越大，同时高密度封装和组装使得元器件的散热空间越来越小，导致热流密度急剧增高。因此作为元器件散热载体的电路板，需要详细分析其散热特性，保证产品的性能和可靠性。

由于电热效应相互影响，电路板上的走线、平面、过孔等产生的焦耳热会影响温度分布，而且如今的电子产品设计中电流越来越大，会进一步增加电热耦合的影响。反过来，温度升高会导致电阻的增大，带来更高的直流压降从而影响元器件供电。因此这两个问题必须同时得到模拟和解决。

板级的传热问题包含了丰富的对流与传导现象，集成电路、封装、电路板、散热器之间的散热主要是热传导问题，而上述因素和环境（气体或液体）之间的散热则是热对流问题。因此，对于板级热分析来讲，不仅要同时分析电气和热物理领域，更要兼顾热传导分析和热对流分析，需要的是一个多物理场的仿真解决方案。

FEA（有限元分析）求解器是用于热传导电热耦合分析的，该方法以传热系数为边界条件，以简化的方式考虑对流和辐射效应，详细模拟固体内部的传导问题，可以在短时间内获得高精度的热传导分析结果。另一种CFD（计算流体动力学）求解器用于热对流和热辐射模拟，通过流体流动的实际模拟（如风扇吹过PCB上的空气）进行对流和辐射的详细建模，但该方法在处理传导问题时，要求尽量简化设计，所以达不到FEA的求解精度和效率。因此我们可以同时使用上述两种仿真方法进行热分析工作，达到优势互补的目的。

FEA求解器主要解决的区域

CFD求解器主要解决的区域

本文章以cadence的Celsius Thermal Solver作为仿真工具，利用它的FEA-CFD电热仿真流程实施分析工作，流程图如下所示：

导入PCB文件，进行电与热相关的参数设置，运行FEA仿真，得到包含PCB各区域功率耗散的简化仿真模型，再导入到CFD求解器中，添加风扇、机箱等结构实施CFD仿真，结果保存为一个CFD模型，代表设置环境（自然环境、风冷或水冷等）下的真实传热系数，再回到FEA求解器中导入CFD模型作为边界条件，重新执行电热仿真，最后得到精确的电热仿真结果。

下面以一个PCB热仿真实例，详细说明仿真步骤。

① FEA求解器提取简化模型

启动Cadence Celsius Thermal Solver，选择并打开Solid Objects Simulation for Layered Structures模块，该模块就是基于FEA求解仿真层状结构（PCB板）。6 ^ u5 R+ R% q1 R

选择菜单栏中Tools > Options > Edit Options，在Translator一栏中使能下图中的两项，保证板框按照design outline、器件尺寸按照place bound定义。& Y J# k) R% b! p9 v: D

在软件界面左侧向导中选择Load Layout File，导入.brd格式的PCB文件，将PCB文件转换为.spd格式的仿真文件。$ \) k- e4 s" U+ C1 Z& {, D8 ?

在向导中使能Enable E/T Co-Simulation Mode选项，选择电热耦合仿真模式。
点击Check Stackup，打开层叠管理界面，检查并设置层叠厚度和材料属性，其中每一层的fill-in填充材料务必设置。此外，点击Pad Stack打开焊盘编辑界面，设置所有的焊盘过孔材料为copper。8 f- v/ ]5 Q0 h2 W

设置电气参数，包括VRM和Sink的定义，通过使能的网络（电源VCC和地网络GND）和器件自动生成VRM和Sink，其中VRM电压为5V，板中所有用电负载电流为2A（参数均为举例）。
" A2 Z' P( K7 x' ^3 Q% S- v

设置环境温度和边界条件，环境温度为25℃，环境风速为0，PCB板水平放置。这里没有设置自然环境的换热系数，因为该步骤主要是为了提取用于CFD放置的功率耗散简化模型，后面用CFD计算出真实的换热系数后，会再回到FEA环境中重新计算。$ V" v# ?+ X, x2 G* y

勾选作为发热件的所有器件，设置每个器件的尺寸参数，包括长宽、厚度，点击Calc Outline and Position，软件会自动根据place bound得到器件的尺寸。8 z: Y5 x$ |# [0 w# E' X

设置每个发热器件的属性，可以是热阻模型参数、也可以是定义元器件的等效材料，这里设置器件的材料为MoldingCompound。
6 x3 f) k% t/ _3 C( _9 e2 ?

设置每个发热器件的功率值，所有器件功率如下图。
E) _7 W0 j$ p. \

点击向导中的Set up Thermal Simulation，打开Thermal Simulation Options，选择Generate Simplified CFD Model，勾选Generate CFD Model选项代表仿真结束后会输出一个供CFD仿真的模型，设置NX和NY两个值，表示将PCB划分为多少个计算功率耗散的区域，区域越多，模型越精细，这里设置为10X10。, c/ [4 \; x0 c5 V7 ^9 K+ \& [

设置计算资源，可选分布式计算或单机计算，设置计算机核心数量和内存，检查是否存在错误，保存仿真文件，然后点击Start Simulation进行仿真。仿真结果的3D热分布图如下所示，可见器件U99的最高温度已经达到170℃，这个温度是不可接受的，因此下面我们在CFD求解器中使用风扇降低这个温度。- z% Q- O$ X8 |% T4 [* ?

刚才的仿真结束后，会在工程目录下生成一个.xml格式的简化模型，用于后面的CFD仿真。
8 `4 f( C/ D2 D9 E& e

② CFD求解器提取CFD模型

Fluid Flow Simulation模块，该模块就是基于CFD求解仿真流体。
; h4 U: ?- g' ?6 f9 o* Z

Modeller> Import Simplifed Model,导入前面生成的.xml格式的简化模型，导入的模型包含PCB和元器件的简化块结构、材料属性、功耗、终端定义等。) n% v$ t! l# O. o0 D, ^

Tools> Edit Library > Material，从库中选择aluminum铝材料并点击Export to Design输出到工程模型里。
s% J( v" h: `8 H2 z

From Design栏中的所有102种材料，右键并选择Bulk Edit，设置Emissivity参数为0.9。该步骤在本案例强迫对流（风扇）应用中并不重要，但如果是仿真自然对流的情况则十分重要。/ E; M+ J) [9 z* X) U

CFD Component下的上下前后左右所有Chassis部件，勾选Installed，说明加入机箱结构，并设置材料、厚度、尺寸和位置参数如下。7 w$ D3 `+ N0 a1 d; m

Chassis下面，点击Chassis_Back右键，选择Add Axial Fan，在机箱Back一侧添加2个风扇，风扇参数设置如下。& O" A3 ~; S/ J- ?. u9 w- K

Chassis_Front右键，选择Add Opening，在机箱Front一侧添加4个通风孔，通风孔参数设置如下。
, I9 }# h( `9 m$ E' R

在CFD Component下面选择Test_Chamber,设置边界盒子，尺寸如下。
% R5 R) F6 T5 N# a8 j; l

Simulation> Options > Simulation Options，设置静态分析中最大迭代数为200，环境温度设置为25℃。
4 ^. f8 N' C. r1 X1 N$ s( w3 g* n$ B

Solver Options，设置CPU最大使用率为67%，Turbulence Model选择Standard K-e model（自然对流使用层流，强制对流使用标准K-e模型），设置Convergence Control为Residual Only。
# I! F8 P4 f9 Y- r

保存工程，点击红色三角形符合开始仿真，仿真结束可以右键Result Summary查看结果。在工程目录下***_EX_CFD文件夹中会自动生成一个包含强制对流的CFD模型。
7 h5 X. ? J3 Y3 t3 S5 B

③ 导入CFD模型到FEA求解器中实施最终电热仿真

FEA求解器的PCB板仿真中，打开之前仿真过设置好的demo.spd文件，在workflow中点击Set up Thermal Simulation，选择Generate Simplified CFD Model界面，取消Generate CFD Model的勾选。
6 w. S/ O* F+ _ x5 }1 W& F7 @

选择Setup Heat Transfer Coefficients界面，使能Use Defined CFD File选项，点击Browse，指向刚才CFD仿真生成的.cfd文件，点击Auto-match by Terminal Name，这样通过CFD仿真得到的、真实准确的换热系数就应用到PCB和元器件表面作为边界条件了。
" o" z$ Q; P4 U% E5 x( Z0 H