4. 仿真前处理实战

封装翘曲与残余应力仿真的精度,往往在几何前处理阶段就已经被决定。CAD/EDA 工具输出的模型面向制造与布线,而有限元求解器需要的是“干净、封闭、可网格化、参数可控”的实体几何。两者之间的差距,就是前处理要填补的部分。本节以典型的倒装芯片 BGA 封装为例,走一遍从原始模型到可求解几何的完整链路,重点讲清楚每一步“为什么做、做到什么程度、怎么判断做完了”。

一、先明确仿真目标,再决定几何粒度

前处理不是盲目地“把模型修干净”,而是围绕分析目标做取舍。翘曲仿真关注的是整体刚度与热失配,残余应力仿真关注的是材料界面附近的应力集中。因此:

  • 芯片、基板、塑封体、焊球这些承担主要刚度和热膨胀失配的部件必须保留实体;
  • 基板内部的铜布线层,若只关心整体翘曲,可等效为均质正交各向异性层,不必逐层建铜;
  • 若关心焊点附近的界面应力,则焊球、IMC 层、underfill 的几何必须保留,且不能过度简化;
  • 键合线、标识字符、倒角、小孔等对刚度贡献极小的特征,一律删除。

这一步的产出是一份“几何清单”:哪些部件保留、哪些等效、哪些删除,以及每个保留部件的尺寸来源(CAD 实测、数据手册还是假设值)。清单确定后,后续清理才有判断依据。

二、导入与单位、坐标系统一

CAD/EDA 模型常见的第一个坑是单位。STEP 文件通常以毫米为单位,但部分 EDA 导出以微米或密尔为单位,导入后模型可能放大或缩小 1000 倍。导入后第一件事是测量一个已知尺寸,例如基板边长,与数据手册核对。若不符,按比例缩放。

坐标系统一同样关键。建议约定:Z 轴垂直于基板平面,向上为正;原点放在基板几何中心。这样后续施加对称约束、读取翘曲量(Z 向位移)时符号含义一致。若原始模型原点在角点,应平移几何而非依赖求解器后处理补偿,避免约束与载荷施加位置出错。

三、几何清理:缝合、去重、消除微小特征

导入的实体常出现面与面之间缝隙、重复面、自相交、微小边。处理顺序建议为:

  1. 缝合(Sew):把公差内的相邻面缝合成封闭实体。缝合公差一般取模型最小特征尺寸的 1/100~1/1000,例如最小特征 0.05 mm 时取 5×10⁻⁵~5×10⁻⁴ mm。公差过大会把本应分离的部件粘在一起,过小则缝不上。
  2. 去重与删除小面:删除面积小于阈值(如 1×10⁻⁴ mm²)的面和长度小于 1 μm 的边,这些特征无法划分网格,只会拖慢求解。
  3. 修复自相交:用布尔运算或面修复工具消除自相交,否则网格划分会在此处失败或产生畸变单元。
  4. 检查实体有效性:确认每个部件都是封闭实体(solid),而非壳体(shell)或面片集合。
判断清理是否完成,可用三个指标:实体数量与几何清单一致;无自由边、无重复面;最小边长远大于网格目标尺寸的 1/5。

四、简化与等效:把“制造细节”变成“力学等效”

简化不是随意删减,而是保持力学响应等价。常见做法:

  • 圆角与倒角:半径小于 0.05 mm 的圆角直接删除,对整体刚度影响可忽略;但焊球与焊盘连接处的圆角若影响界面应力,应保留。
  • 过孔与布线:基板内过孔阵列可等效为各向异性材料属性,用铜占比计算等效弹性模量和热膨胀系数,而不是逐个建孔。
  • 薄层合并:多层铜箔与介质层可合并为一层等效材料,厚度取各层之和,属性按体积分数加权。
  • 对称利用:若几何与载荷均关于 XZ、YZ 平面对称,可只建 1/4 模型,在对称面施加对称约束,显著降低规模。

简化后必须做一次“质量守恒”核对:等效层的等效密度乘以体积,应接近原始各层质量之和;若偏差超过 5%,说明等效参数或厚度取值有误。

五、参数化建模:让几何可扫描、可复用

封装仿真的价值在于比较不同设计或工艺参数。因此几何不应是“死”的,而应由参数驱动。建议把以下量设为参数:芯片尺寸、基板尺寸与厚度、塑封体厚度、焊球直径与间距、underfill 厚度、各层材料厚度。参数命名要能自解释,例如 sub_thk、die_size_x、solder_d。

# 参数化建模示例(伪代码,单位 mm)
die_size_x   = 8.0      # 芯片 X 向尺寸
die_size_y   = 8.0
die_thk      = 0.10     # 芯片厚度
sub_size_x   = 14.0     # 基板 X 向尺寸
sub_thk      = 0.40     # 基板厚度
mold_thk     = 0.60     # 塑封体厚度
solder_d     = 0.25     # 焊球直径
solder_pitch = 0.50     # 焊球间距
# 由参数生成几何,修改参数即可重建模型

参数化后,几何重建应能一键完成。每次改参数后,仍需重新检查实体有效性与最小边,因为参数变化可能引入新的微小特征。

六、前处理流程与判断节点

把上述步骤串起来,可以用下面的信号链式流程图表示。图中每个节点都对应一个可检查的判据,只有通过判据才进入下一步,否则回到清理或简化环节。

封装CAD/EDA模型 → 仿真几何前处理链路 CAD/EDA 原始模型 STEP / Parasolid / ODB++ 单位与坐标系统一 Z 轴垂直基板,原点居中 几何清理 缝合 / 去重 / 去小特征 简化与等效 各向异性 / 对称 / 合并 参数化建模 尺寸/厚度/间距驱动 几何有效性检查 封闭实体 / 无自由边 网格划分与求解 输出可求解几何 检查不通过:返回清理/简化 参数变更后重建:重新检查实体有效性与最小边 关键判据: ① 实体数量与几何清单一致;② 无自由边、无重复面;③ 最小边 ≥ 网格目标尺寸的 1/5; ④ 等效层质量偏差 ≤ 5%;⑤ 参数重建后几何仍为封闭实体。

七、常见问题与处理建议

实际前处理中,以下问题出现频率最高:

现象可能原因处理建议
导入后模型尺寸差 1000 倍单位不一致(mm 与 μm)测量已知尺寸核对,按比例缩放
缝合后部件粘连缝合公差过大减小公差至最小特征的 1/1000
网格划分失败存在微小边或自相交删除小特征,修复自相交
翘曲量符号与预期相反Z 轴方向或原点定义不一致统一约定并检查约束施加面
参数修改后结果跳变重建引入新微小特征每次重建后重新执行有效性检查

前处理的本质,是把“制造视角的模型”翻译成“力学视角的模型”。翻译过程中每一次简化都要有依据,每一次修复都要有判据。做到这一点,后续的网格划分和求解才会稳定,翘曲与残余应力的结果才具备可比性和可复现性。