4. 仿真前处理实战
封装翘曲与残余应力仿真的精度,往往在几何前处理阶段就已经被决定。CAD/EDA 工具输出的模型面向制造与布线,而有限元求解器需要的是“干净、封闭、可网格化、参数可控”的实体几何。两者之间的差距,就是前处理要填补的部分。本节以典型的倒装芯片 BGA 封装为例,走一遍从原始模型到可求解几何的完整链路,重点讲清楚每一步“为什么做、做到什么程度、怎么判断做完了”。
一、先明确仿真目标,再决定几何粒度
前处理不是盲目地“把模型修干净”,而是围绕分析目标做取舍。翘曲仿真关注的是整体刚度与热失配,残余应力仿真关注的是材料界面附近的应力集中。因此:
- 芯片、基板、塑封体、焊球这些承担主要刚度和热膨胀失配的部件必须保留实体;
- 基板内部的铜布线层,若只关心整体翘曲,可等效为均质正交各向异性层,不必逐层建铜;
- 若关心焊点附近的界面应力,则焊球、IMC 层、underfill 的几何必须保留,且不能过度简化;
- 键合线、标识字符、倒角、小孔等对刚度贡献极小的特征,一律删除。
这一步的产出是一份“几何清单”:哪些部件保留、哪些等效、哪些删除,以及每个保留部件的尺寸来源(CAD 实测、数据手册还是假设值)。清单确定后,后续清理才有判断依据。
二、导入与单位、坐标系统一
CAD/EDA 模型常见的第一个坑是单位。STEP 文件通常以毫米为单位,但部分 EDA 导出以微米或密尔为单位,导入后模型可能放大或缩小 1000 倍。导入后第一件事是测量一个已知尺寸,例如基板边长,与数据手册核对。若不符,按比例缩放。
坐标系统一同样关键。建议约定:Z 轴垂直于基板平面,向上为正;原点放在基板几何中心。这样后续施加对称约束、读取翘曲量(Z 向位移)时符号含义一致。若原始模型原点在角点,应平移几何而非依赖求解器后处理补偿,避免约束与载荷施加位置出错。
三、几何清理:缝合、去重、消除微小特征
导入的实体常出现面与面之间缝隙、重复面、自相交、微小边。处理顺序建议为:
- 缝合(Sew):把公差内的相邻面缝合成封闭实体。缝合公差一般取模型最小特征尺寸的 1/100~1/1000,例如最小特征 0.05 mm 时取 5×10⁻⁵~5×10⁻⁴ mm。公差过大会把本应分离的部件粘在一起,过小则缝不上。
- 去重与删除小面:删除面积小于阈值(如 1×10⁻⁴ mm²)的面和长度小于 1 μm 的边,这些特征无法划分网格,只会拖慢求解。
- 修复自相交:用布尔运算或面修复工具消除自相交,否则网格划分会在此处失败或产生畸变单元。
- 检查实体有效性:确认每个部件都是封闭实体(solid),而非壳体(shell)或面片集合。
四、简化与等效:把“制造细节”变成“力学等效”
简化不是随意删减,而是保持力学响应等价。常见做法:
- 圆角与倒角:半径小于 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 # 焊球间距
# 由参数生成几何,修改参数即可重建模型
参数化后,几何重建应能一键完成。每次改参数后,仍需重新检查实体有效性与最小边,因为参数变化可能引入新的微小特征。
六、前处理流程与判断节点
把上述步骤串起来,可以用下面的信号链式流程图表示。图中每个节点都对应一个可检查的判据,只有通过判据才进入下一步,否则回到清理或简化环节。
七、常见问题与处理建议
实际前处理中,以下问题出现频率最高:
| 现象 | 可能原因 | 处理建议 |
|---|---|---|
| 导入后模型尺寸差 1000 倍 | 单位不一致(mm 与 μm) | 测量已知尺寸核对,按比例缩放 |
| 缝合后部件粘连 | 缝合公差过大 | 减小公差至最小特征的 1/1000 |
| 网格划分失败 | 存在微小边或自相交 | 删除小特征,修复自相交 |
| 翘曲量符号与预期相反 | Z 轴方向或原点定义不一致 | 统一约定并检查约束施加面 |
| 参数修改后结果跳变 | 重建引入新微小特征 | 每次重建后重新执行有效性检查 |
前处理的本质,是把“制造视角的模型”翻译成“力学视角的模型”。翻译过程中每一次简化都要有依据,每一次修复都要有判据。做到这一点,后续的网格划分和求解才会稳定,翘曲与残余应力的结果才具备可比性和可复现性。