4. 3D IC封装几何建模与简化准则
3D IC 的热结构耦合仿真，第一步不是画网格，而是把版图、封装和装配信息“翻译”成可计算的几何模型。这一步决定了后续网格规模、求解时间以及结果可信度。几何建模的核心矛盾在于：真实封装包含数百万个凸点、微凸点、TSV、走线和介质层，若全部保留，网格量会爆炸；若删减过度，热流路径和应力集中位置又会失真。因此，本节围绕“细节删除、等效建模、对称性与周期性利用、网格过渡处理”四条准则展开。
一、从版图与封装文件到几何模型
典型输入包括：芯片 GDSII/OASIS 版图、封装基板 Gerber 或 DXF、BGA 球阵列坐标、TSV 阵列文件、材料层叠表（layer stack-up）以及键合/装配公差。仿真前需要先建立统一的坐标系约定：以封装基板底面中心为原点，X/Y 为平面方向，Z 向上为正，芯片堆叠方向沿 +Z。所有文件导入后必须对齐到同一原点，否则后续对称面、周期边界都会错位。
导入后的几何通常处于“过度详细”状态。例如一个 20 mm × 20 mm 的封装基板，若保留全部 2000 个 BGA 焊球、数万条铜走线和阻焊开窗，直接划分四面体网格可能产生上千万单元。此时需要按热-结构贡献度分级处理，而不是简单按尺寸删除。
二、细节删除准则
细节删除不是“越小越删”，而是判断该特征是否影响温度梯度、热流路径或应力集中。可参考以下判据：
- 尺寸判据：特征尺寸小于目标网格尺寸的 1/3，且不位于关键热流路径上，可删除。例如阻焊层开窗、表面文字、测试焊盘。
- 热阻判据：若某结构的热阻远小于相邻主路径热阻（如小于 1%），可忽略。例如芯片表面钝化层相对于硅衬底。
- 刚度判据：若某细节对整体刚度贡献极小，且不产生局部应力集中，可删除。例如非承载的塑封料溢边。
- 保留例外：TSV、微凸点、关键互连、散热通孔、芯片边缘和角部必须保留，因为这些位置往往是热流瓶颈和应力集中区。
三、等效建模方法
对于数量庞大但几何重复的结构，如 BGA 球阵列、微凸点阵列、TSV 阵列，直接逐个建模不现实，应采用等效建模。常用方法有三种：
1. 等效热导率法：将阵列区域视为均匀材料，计算其等效导热系数。以 TSV 阵列为例，铜填充 TSV 直径 10 μm、间距 40 μm、硅基体导热系数 148 W/(m·K)、铜 400 W/(m·K)，则面内等效导热系数可近似为面积加权：
k_eq = f_Cu * k_Cu + (1 - f_Cu) * k_Si
f_Cu = π*(d/2)^2 / p^2 = π*25 / 1600 ≈ 0.049
k_eq ≈ 0.049*400 + 0.951*148 ≈ 160.3 W/(m·K)该值明显高于纯硅，说明 TSV 阵列在面内增强了导热，但沿 Z 向的等效热阻还需考虑 TSV 长度和界面层。
2. 等效热阻网络法：将焊球、微凸点、TIM 层表示为热阻节点，用 R = t/(k·A) 计算。例如一个直径 0.3 mm、高度 0.25 mm 的 SnAgCu 焊球，k ≈ 50 W/(m·K)，则单球热阻约为 0.25e-3/(50×π×0.15e-3²) ≈ 70.7 K/W。若阵列有 1000 个球并联，总热阻约 0.07 K/W。
3. 等效力学模型：对结构仿真，焊球阵列可用等效弹性模量和泊松比替代，或采用“子模型”技术：全局模型用等效层，局部关键焊球用详细模型二次分析。
四、对称性与周期性利用
3D IC 封装通常具有几何和载荷对称性。若芯片布局、TSV 阵列和热载荷均关于 XZ 平面对称，可只建 1/2 模型；若关于 XZ 和 YZ 双对称，可建 1/4 模型。对称面需施加对称边界条件：法向位移为零，面内旋转自由。热分析中对称面为绝热边界。
对于周期性阵列，如 HBM 堆叠中的重复存储层，可取一个周期单元（unit cell），在侧面施加周期性边界条件：对应面温度差相等、位移协调。这样可将 8 层堆叠简化为 1 层单元，计算量降低近一个数量级。但需注意：周期边界要求载荷也周期分布，若存在全局温度梯度或边缘效应，则不能简单使用。
五、网格过渡处理
简化后的模型仍存在尺度跨度：TSV 直径约 10 μm，封装尺寸约 20 mm，相差三个数量级。若全局用 10 μm 网格，单元数不可接受；若全局用 200 μm，TSV 区域又无法捕捉热流集中。因此必须做网格过渡。
推荐策略是“分区过渡 + 渐变比控制”：在 TSV、微凸点、芯片边缘等关键区采用 0.5–2 μm 的细网格；在基板和塑封体远场采用 50–200 μm 的粗网格；两者之间设置过渡层，相邻单元尺寸比不超过 1.5，避免刚度矩阵病态和热流失真。过渡层厚度一般取细网格尺寸的 3–5 倍。
对于热-结构耦合，还需注意网格类型匹配：热分析可用四面体，结构分析在焊球处宜用六面体或棱柱层，以减少剪切锁定。若使用同一网格，应在焊球界面处加密，并检查雅可比行列式是否为正。
原始模型：2000 个微凸点 + 500 个 BGA 球 + 完整走线，四面体单元约 1200 万，求解约 6 小时。
简化后：微凸点阵列用等效层，BGA 用热阻网络，走线删除，保留 TSV 和边缘焊球，单元约 180 万，求解约 25 分钟。
验证：关键芯片结温差异 < 1.5%，最大 von Mises 应力位置一致，数值差异 < 8%。
六、建模检查清单
完成几何简化后，建议按以下清单逐项检查，避免把错误带入求解阶段：
| 检查项 | 合格判据 | 常见问题 |
|---|---|---|
| 坐标与单位 | 所有部件原点一致，单位统一为 mm | 版图与封装错位、μm 与 mm 混用 |
| 细节删除 | 主热路径热阻变化 < 5% | 误删散热通孔或边缘焊球 |
| 等效参数 | k_eq、E_eq 有计算依据并可追溯 | 直接拍脑袋取值 |
| 对称/周期 | 对称面边界条件正确，载荷也对称 | 只几何对称、载荷不对称仍用 1/4 模型 |
| 网格过渡 | 相邻单元尺寸比 ≤ 1.5，关键区足够细 | 过渡太陡导致热流震荡 |
| 网格无关性 | 加密后关键温度变化 < 2% | 只做一次网格就下结论 |
几何建模与简化不是一次性动作，而是与网格、求解器、后处理反复迭代的过程。实际项目中，建议先建立“粗模型”快速验证热路径和边界条件，再逐步加密关键区域，最后用子模型或局部详细模型校核。这样既能控制计算成本，又能保证关键位置的精度。