全封装氧化锆基氧传感器的力学应力致失效行为的计算与实验研究

《Results in Engineering》：Computational and experimental study on the mechanical-stress induced failure behavior of fully packaged zirconia-based oxygen sensors

　　

研究人员针对汽车与工业领域广泛应用的氧化锆基氧传感器在长期服役中受循环热载荷与机械应力作用导致可靠性不足的问题，建立了包含全部封装组件及其与陶瓷芯片相互作用的电-热-力耦合模型。该模型揭示铂加热元件与周围陶瓷间的热失配应力是决定结构完整性的主导因素：在6.2 V加热电压下，陶瓷芯片峰值温度达703.1°C，Z形加热电阻的几何过渡区最大第一主应力为601.7 MPa。与既往仅开展温度层级验证的研究不同，本工作从三个独立维度完成模型验证：热成像验证（峰值温度偏差0.4%）、三点弯曲试验结合威布尔（Weibull）可靠性分析（样本量n=9）、加速失效试验结合裂纹位置统计验证（样本量n=29）。参数敏感性分析证实铂（Pt）与钇稳定氧化锆（YSZ）的热膨胀系数（CTE）失配对应力产生起主导作用。通过全局对称电极重新设计的结构优化方案，将峰值应力降至444.5 MPa（降幅26.1%），安全系数从1.20提升至1.62。这种融合全封装传感器多物理场建模、多层次实验验证、敏感性分析与机理导向优化的集成方法，将当前传感器可靠性评估从定性的芯片级分析推进至定量、具备统计支撑的封装级评价。

本研究由南京大学的Wei Ding、Ran Li、Xihua An、Peng Li、Can Lu与Baoming Wang合作完成，发表于《Results in Engineering》。氧化锆基氧传感器凭借高灵敏度、宽工作温区与耐久性，被广泛应用于汽车尾气监测、工业流程控制及医疗诊断等高价值场景。然而其服役环境伴随剧烈热循环与机械振动，长期可靠性常因热机械疲劳、材料性能退化与封装残余应力而受限，现有研究多采用简化或半封装模型，对封装引入的机械约束考虑不足，导致实际失效预测精度有限。为此研究人员构建了全封装结构的电-热-力耦合数值体系，并结合多层级实验完成了机理解析与结构优化。

关键技术方法包括：采用COMSOL Multiphysics建立全封装氧传感器的三维多物理场耦合模型，完整纳入金属基座、密封粉末环、聚四氟乙烯（PTFE）封装件及陶瓷芯片的相互作用；通过网格独立性验证确保数值稳定性，采用稳态分离式求解器依次计算焦耳热、温度场与热应力分布；实验层面构建三组独立验证体系，包括6.2 V工况下的红外热成像测温、n=9的三点弯曲强度测试及威布尔统计分析、7.8 V过压加速失效试验结合n=29的裂纹位置统计；基于参数敏感性分析量化材料属性与边界条件对应力的影响权重，并以全局对称电极设计为优化路径开展结构改进。

研究结果如下：

建模与仿真分析部分，研究人员首先明确了传感器多层陶瓷结构（含氧化铝绝缘层、YSZ固态电解质层、铂电极与加热元件）与全封装组件的几何尺寸及材料参数（密度、比热容、热导率、弹性模量、泊松比、热膨胀系数）。建立的电-热-力耦合数学模型包含电流场（焦耳热效应，考虑铂电阻率随温度变化）、固体传热场（传导、对流与表面间辐射）及固体力学场（线性弹性本构与热膨胀耦合）。网格独立性验证表明，网格数量从541437增至2031474时，峰值温度、最大位移与von Mises应力变化率分别仅为0.02%、0.08%与1.02%，证明数值解收敛可靠。仿真结果显示，6.2 V工况下陶瓷芯片峰值温度为703.1°C，位于加热元件附近；铂加热层最大第一主应力达398.3 MPa，远高于上下氧化铝绝缘层的226.9 MPa与224.6 MPa；Z形加热电阻宽窄段过渡区的几何不连续引发显著应力集中，全芯片最大第一主应力达601.7 MPa。参数敏感性分析表明，YSZ与Pt的热膨胀系数（CTE）对应力影响最显著，二者失配是应力产生的主导因素；YSZ热导率的变化主要影响温度分布，而Al₂O₃与Pt的热导率及CTE变化对应力影响微弱；对流换热系数在5.0~8.0 W/(m²·K)范围内变化时，峰值温度与应力仅有微小波动。

多物理场模型的实验验证部分，温度场验证显示拆解传感器芯片的实测峰值温度为719°C，仿真值为722°C，偏差仅0.4%，且不同加热电压下仿真与实验温度趋势一致。应力相关验证中，n=9的三点弯曲试验测得陶瓷芯片平均抗弯强度为762.4 MPa（标准差22.0 MPa，变异系数2.89%），服从Weibull分布（形状参数m=44.2，尺度参数λ=771.9 MPa），仿真所得601.7 MPa工作应力对应的断裂概率低至1.65×10^-5，具备充足安全裕度。裂纹位置统计中，n=29的加速失效试验显示裂纹集中于距芯片尖端3.46~4.45 mm区域，正面与背面裂纹平均距离分别为4.109 mm与3.839 mm，与仿真预测的Z形电阻过渡区应力集中位置高度吻合，验证了失效定位的准确性。

工程应用与优化方法部分，研究人员发现原电极引线非对称布局导致左右铂分布失衡，引发附加弯曲力矩。通过将4条非对称引线改为3条引线加1个表面焊盘的中心线对称布局，实现了铂分布的镜像平衡。优化后芯片最大第一主应力降至444.5 MPa（降幅26.1%），安全系数从1.20提升至1.62，峰值温度仅上升2°C，电化学界面功能未受影响，仅制造复杂度略有增加。

讨论与结论部分，研究人员将本工作与既有研究对比指出，既往模型多采用1/4或1/2芯片简化几何，缺乏封装约束与实验验证，本工作的全封装模型因纳入机械约束，所得应力水平更符合实际工况。最终结论可归纳为五点：全封装电-热-力耦合模型成功预测了Z形加热电阻过渡区的应力集中现象；三层级的实验验证体系确保了模型可靠性，裂纹位置一致性证实了失效机理的准确性；Pt-YSZ热膨胀系数失配是应力主导因素，且应力集中位置对参数扰动具有鲁棒性；对称电极优化有效降低了26.1%的峰值应力，显著提升安全裕度；本模型为静态断裂风险评估框架，未来可通过纳入温度依赖属性与循环载荷进一步拓展至疲劳寿命预测。该研究提出的集成方法为多层陶瓷器件的可靠性评估提供了可迁移的工程范式。