2. 封装结构建模基础

封装电磁热结构联合仿真的第一步,是把真实封装体转化为仿真器能够求解的模型。这一步决定了后续所有计算的物理可信度:几何过度简化会丢失关键效应,材料属性取值不当会让结果偏离实际,网格过粗则无法捕捉梯度,过细则让计算量失控。本节围绕几何简化、材料属性、网格划分与精度权衡四条主线,说明建模阶段应当做什么、怎么做,以及如何判断模型是否够用。

一、几何简化:保留主导路径,剔除次要细节

封装体通常包含芯片、焊料层、铜柱/凸点、基板、塑封料、散热盖、焊球等结构,还伴随大量倒角、通孔、走线、标记和圆角。若全部按原始CAD导入,往往产生数十万甚至上百万个面,网格划分困难且求解规模不可控。几何简化的原则是:保留对电磁场、热流和应力传递起主导作用的连续路径,剔除对目标物理量影响低于设定阈值的局部特征。

具体操作上,可按以下优先级处理:

  • 必须保留:芯片有源区、焊料/铜柱连接层、基板铜箔层、散热路径上的金属块、影响电感回路的键合线或走线。
  • 可合并:多层薄铜箔若厚度远小于热扩散长度,可等效为一层各向异性材料;多个小焊球可等效为均质焊料层。
  • 可删除:倒角、丝印、定位孔、非功能圆角、外部标记。
  • 可对称化:若结构、载荷和边界条件均对称,可取1/2或1/4模型,但需确认电磁激励是否破坏对称性。

例如,某倒装芯片封装中,芯片尺寸为5 mm×5 mm×0.2 mm,下方有100个直径80 μm、高度60 μm的铜柱。若逐个建铜柱,网格量巨大;若等效为厚度60 μm、铜体积分数约35%的均质层,则热导率和电导率可按体积分数加权估算。这种简化在热分析中误差通常小于5%,但在高频电磁分析中可能丢失铜柱间的互感细节,因此是否简化必须结合分析类型判断。

二、材料属性:区分常数、温度相关与各向异性

封装涉及金属、半导体、聚合物和焊料等多类材料,其属性不能一律按常温常数处理。建模时应至少明确三类属性:

  • 电磁属性:电导率σ、相对介电常数εr、损耗角正切tanδ、磁导率μr。铜的电导率常取5.8×10⁷ S/m,但高频下需考虑趋肤效应,必要时用表面阻抗边界替代实体铜。
  • 热属性:热导率k、比热容Cp、密度ρ。硅的k约148 W/(m·K),铜约400 W/(m·K),塑封料仅0.5~1.5 W/(m·K),焊料约50 W/(m·K)。
  • 结构属性:杨氏模量E、泊松比ν、热膨胀系数CTE。硅的CTE约2.6 ppm/K,铜约17 ppm/K,塑封料约8~15 ppm/K,焊料约21 ppm/K。

温度相关属性在热-结构联合仿真中尤其重要。例如焊料在接近熔点附近E和屈服强度显著下降,若按常温常数计算,热应力会被高估。各向异性也不可忽略:塑封料在注塑流动方向与垂直方向的热导率和CTE可相差2倍以上,层压基板的X/Y方向与Z方向热导率差异更大。建议在材料库中为每种材料标注“常数/温度相关/各向异性”三种状态,并在仿真日志中记录取值来源。

教学示例:某塑封料在25℃时k=0.8 W/(m·K),在125℃时k=1.1 W/(m·K)。若芯片功耗2 W、结温预计从25℃升到125℃,用常数k=0.8会高估热阻约15%,导致结温预测偏高。此时应使用温度相关曲线或至少取工作区间中值。

三、网格划分策略:按梯度分布而非均匀加密

网格是几何与物理场的离散载体。封装仿真的网格难点在于尺度跨度大:芯片厚度0.2 mm,焊料层0.06 mm,而基板可达10 mm以上,跨尺度比超过100:1。均匀加密会让单元数爆炸,合理做法是按物理量梯度分布网格。

具体策略包括:

  • 边界层网格:在铜-塑封料、焊料-芯片界面附近,热流和电流梯度最大,应布置3~5层薄单元,首层厚度可取特征长度的1/10。
  • 扫掠网格:对芯片、基板等规则六面体区域,用扫掠网格生成结构化单元,质量高且数量少。
  • 局部细化:仅在焊点、铜柱、键合线根部等应力/电流集中区细化,其余区域保持较粗网格。
  • 过渡控制:相邻区域网格尺寸比不宜超过2:1,避免出现畸形单元导致收敛困难。

网格数量并非越多越好。以热分析为例,当网格从5万增加到20万时,最高结温变化若小于0.5%,说明已进入网格无关区间;若变化仍大于2%,则需继续细化或检查边界层设置。电磁分析中还需关注趋肤深度δ=√(2/(ωμσ)),铜在1 GHz时δ约2 μm,若网格尺寸远大于δ,则表面电阻计算失真。

封装建模四要素与精度权衡关系 几何简化 保留主导路径 合并薄层/小焊球 删除倒角/标记 对称面取1/2或1/4 影响:面数、网格量 材料属性 σ、εr、tanδ、μr k、Cp、ρ E、ν、CTE 温度相关/各向异性 影响:物理场幅值 网格划分 边界层3~5层 扫掠六面体 局部细化焊点/铜柱 过渡比≤2:1 影响:单元数、收敛性 精度权衡 网格无关性验证 结温变化<0.5% 趋肤深度δ校核 计算资源预算 目标:够用即可 网格无关性验证示意(教学示例) 网格数量(万) 最高结温(℃) 0 5 10 15 20 25 120 125 130 135 网格无关区间:结温变化 < 0.5% 继续加密收益有限,可停止细化 粗网格:结温偏高 梯度未捕捉 注:曲线为教学示意,实际数值随封装结构、功耗与边界条件变化

四、精度权衡:以目标误差反推建模粒度

建模精度不是越高越好,而是要与分析目标匹配。工程上常用“目标误差反推法”:先明确关注量(如结温、插入损耗、最大等效应力)允许的误差范围,再决定几何、材料和网格的粒度。

例如,若热分析目标是结温预测误差不超过3℃,而材料属性不确定性已贡献约2℃,则网格和几何简化带来的额外误差应控制在1℃以内。此时可通过一次粗网格和一次细网格的对比,确认网格误差是否满足要求。若电磁分析关注回波损耗S11,则需保证端口网格和传输线截面网格能分辨趋肤深度,否则S11相位误差会明显增大。

实际操作中,建议按以下顺序迭代:

  1. 用简化几何和常数材料建立基线模型,快速求解,观察物理量量级是否合理。
  2. 对关注区域局部加密,检查结果变化是否收敛。
  3. 引入温度相关材料或各向异性,比较与常数模型的差异。
  4. 若差异超过目标误差,则保留高精度设置;若差异很小,则回退到简化模型以节省计算资源。
常见误区:为了“更准”而全局加密网格,导致单元数从10万涨到200万,求解时间增加数十倍,但结温仅变化0.2℃。正确做法是只在梯度大的界面和集中区加密,其余区域保持粗网格。

五、建模检查清单

完成建模后,可用以下清单快速自检:

检查项 判断标准 不满足时的处理
几何是否保留主导热/电/力路径 删除次要特征后,关注量变化小于目标误差的1/3 恢复被删特征或改用等效层
材料属性是否覆盖工作温度区间 关键材料k、E、CTE随温度变化超过10%时需用曲线 替换为温度相关数据
界面网格是否分辨梯度 边界层首层厚度≤特征长度的1/10 增加边界层数或减小首层厚度
网格无关性是否验证 加密后关注量变化<目标误差 继续细化或检查边界条件
电磁趋肤深度是否分辨 导体网格尺寸≤δ/2 使用表面阻抗边界或加密导体表层

几何简化、材料属性、网格划分和精度权衡四者相互制约。简化越激进,网格越容易划分,但可能丢失关键物理;材料越精细,结果越可信,但数据获取成本越高;网格越密,梯度捕捉越好,但计算资源消耗越大。建模的本质是在目标误差约束下寻找最经济的组合,而不是追求单一环节的极致精度。