案例丨热浪中的变形之谜：用PERA SIM Mechanical揭秘PCB热翘曲

3.1 模型建立及简化

本文研究对象为电路板，根据分析目的，对电路板模型进行了简化，包括去除不必要的圆角、倒角，删除与分析相关性较低的零部件。为了保证计算精度，整个模型采用实体单元进行模拟。 

分析的整体几何模型如图1所示：

图1  电路板稳态热－固耦合分析的几何模型

3.2 网格划分

网格划分以基本保持模型原有的几何形状作为前提，为了保持网格划分后与几何的贴体性，同时结合几何的基本形状，本文采用高阶四面体网格。使用平均2mm网格、最大4mm网格进行整体网格划分，共生成85921个高阶四面体网格，共计165799个节点，生成的网格模型如下：

图2  电路板热－固耦合仿真网格

3.3 材料定义

本文利用热－固耦合分析来研究PCB板的受热翘曲问题，需要先计算PCB板的温度分布，然后再基于温度分布来求解PCB板的热翘曲。默认所有焊接在电路板上的器件与电路板之间的接触关系均为绑定接触，材料本构均为线弹性，对应的材料参数如下表1所示。

表1 电路板稳态热-固仿真分析涉及的材料参数

3.4 截面定义

采用实体单元。

3.5 边界条件

模型的边界条件为侧面与电路板安装孔的固定，具体位置如图3所示。

图3 电路板固定支撑的位置

3.6 热源施加

根据电路板及其器件的工况，共涉及5个热源，具体的热源大小与位置如图4所示：

图4  电路板热源施加图

3.7 对流换热边界施加

根据电路板及其器件的工况，一共涉及3个对流换热边界，具体的对流换热条件与位置如下图5所示：

对流换热1

对流换热2

对流换热3

图5 电路板对流换热条件施加图

4.1 计算分析设置

1）在完成热源、固定约束条件以及对流条件的指定后，创建对应的稳态热仿真，配置稳态热仿真分析工况；

2）创建一个稳态热-固耦合仿真分析工况，以稳态热作为载荷条件，以固定约束条件作为约束工况。

4.2 稳态热-固耦合计算结果

完成稳态热－固耦合配置后，提交作业进行计算，在后处理模块加载并查看计算结果。

计算得到的电路板温度分布如图6所示，从图中可以看出，电路板在热源发热以及对应的对流换热条件下，稳定的最高温度为86℃。 

图6  电路板温度分布云图

电路板受热后的结构翘曲变形（位移分布）如图7所示，升温引起的最大结构翘曲为1.72e-4 m。

图7  电路板受热后的结构翘曲