电子封装热仿真：热应力分析_（5）.热应力的基本概念与计算方法.docxVIP

PAGE1

PAGE1

热应力的基本概念与计算方法

热应力的定义

热应力是指由于温度变化引起的材料内部应力。在电子封装中，由于不同材料之间的热膨胀系数（CTE，CoefficientofThermalExpansion）不同，当温度变化时，不同材料之间的热膨胀量不一致，从而产生热应力。这种应力可能会导致封装结构的变形、裂纹甚至失效，因此在设计和优化电子封装时，热应力分析是至关重要的一步。

热应力的产生机制

温度变化引起的材料热膨胀

当温度变化时，材料会发生活动，这种活动可以用热膨胀系数来描述。热膨胀系数定义为材料在单位温度变化时长度的变化率。不同材料的热膨胀系数差异越大，温度变化时产生的热应力也就越大。

材料之间的热膨胀不匹配

在电子封装中，通常涉及多种材料，如硅芯片、焊料、基板、封装外壳等。这些材料的热膨胀系数不同，当温度变化时，材料之间的热膨胀量不一致，从而产生热应力。例如，硅芯片的热膨胀系数约为2.6ppm/°C，而常用的基板材料如FR-4的热膨胀系数约为17ppm/°C。因此，当温度上升时，FR-4的膨胀量远大于硅芯片，导致硅芯片受到压缩应力，而当温度下降时，硅芯片的收缩量小于FR-4，导致硅芯片受到拉伸应力。

约束条件

在实际封装结构中，材料之间的约束条件也会影响热应力的大小。例如，焊点将芯片固定在基板上，当温度变化时，焊点会限制芯片的自由热膨胀或收缩，从而产生额外的热应力。这些约束条件可以通过边界条件在热仿真中进行建模。

热应力的计算方法

线性热应力计算

对于简单的几何形状和材料，可以使用线性热应力计算公式来估算热应力。线性热应力计算公式如下：

σ

其中：-σthermal是热应力-E是材料的弹性模量-α是材料的热膨胀系数-Δ

例子：简单线性热应力计算

假设有一个硅芯片，其尺寸为10mm×10mm×1mm，热膨胀系数为2.6ppm/°C，弹性模量为130GPa。当温度从25°C升高到150°C时，计算产生的热应力。

#定义材料属性

E_silicon=130e9#硅的弹性模量，单位为Pa

alpha_silicon=2.6e-6#硅的热膨胀系数，单位为1/°C

delta_T=150-25#温度变化，单位为°C

#计算热应力

sigma_thermal=E_silicon*alpha_silicon*delta_T

#输出结果

print(f产生的热应力为:{sigma_thermal/1e6}MPa)

非线性热应力计算

对于复杂的几何形状和材料，尤其是涉及多材料和多层结构的电子封装，线性热应力计算方法可能不够准确。此时需要使用非线性热应力计算方法，考虑材料的非线性特性、几何非线性和接触非线性等因素。

有限元方法（FEM）

有限元方法是目前最常用的非线性热应力计算方法。通过将复杂的结构分解为多个小的单元，每个单元的应力和应变可以通过线性或非线性方程组求解，最终得到整个结构的热应力分布。

例子：使用有限元方法进行热应力分析

假设有一个三层封装结构，包括硅芯片、焊料层和基板。使用有限元软件（如ANSYS、ABAQUS等）进行热应力分析。

模型建立：在有限元软件中建立三层封装结构的几何模型。

材料属性：定义每个材料的热膨胀系数和弹性模量。

边界条件：设置温度变化的边界条件。

求解：运行仿真，得到热应力分布结果。

#示例代码：使用Python和FEniCS库进行简单的热应力分析

#导入必要的库

fromfenicsimport*

#定义几何模型

mesh=UnitSquareMesh(10,10)

V=VectorFunctionSpace(mesh,P,1)

#定义材料属性

E_silicon=130e9#硅的弹性模量，单位为Pa

nu_silicon=0.22#硅的泊松比

alpha_silicon=2.6e-6#硅的热膨胀系数，单位为1/°C

E_solder=30e9#焊料的弹性模量，单位为Pa

nu_solder=0.33#焊料的泊松比

alpha_solder=23e-6#焊料的热膨胀系数，单位为1/°C

E_substrate=18e9#基板的弹性模量，单位为Pa

nu_substrate=0.3#基板的泊松比

alpha_substrate=17e-6#基板的热膨胀系数，单位为1/°C

#定义温度变化

delta_T=150-25#温度变化，单位为°C

#定义弹性张量

defepsilon(v):

returnsym(grad