SAC305焊料合金全寿命单轴棘轮行为的本构与损伤模型研究

研究概述

SAC305（Sn-3.0Ag-0.5Cu）无铅焊料因其优异的润湿性、可靠的力学性能和抗疲劳性能，已成为电子封装领域替代含铅焊料的材料。随着电子设备向微型化和轻量化方向发展，焊点的尺寸不断减小，在运输和服役过程中承受的热循环载荷、振动、冲击等复杂载荷导致其易于发生疲劳失效。在非对称循环载荷下，焊点会产生塑性应变持续累积的"棘轮效应"，加速材料失效，严重影响电子封装的可靠性。传统的统一黏塑性本构模型基于结构均匀性假设，仅能描述棘轮安定行为，无法模拟材料从稳态到加速失效的全寿命过程。因此，考虑材料内部微裂纹和微缺陷演化的损伤耦合本构模型对于准确预测焊料寿命至关重要。

针对上述问题，天津大学化工学院陈刚、王宇燕、李鹏、崔云、施守文、杨静、许维玲、林强等人在Mechanics of Materials期刊（2022年第171卷）发表了题为"Constitutive and damage model for the whole-life uniaxial ratcheting behavior of SAC305"的研究论文。该研究通过系统的单轴棘轮实验，建立了基于连续损伤力学的非线性损伤耦合本构模型，以描述SAC305焊料在不同加载条件下的全寿命棘轮行为。损伤变量从循环加载过程中卸载段的弹性模量退化中提取，损伤演化基于损伤耗散势推导，并通过归一化寿命和塑性应变能密度进行修正，使模型能够准确描述不同加载条件下的损伤演化。该研究为电子封装焊料在复杂服役环境下的寿命预测和完整性评估提供了重要的理论基础和方法支撑。

实验方法与核心发现

研究采用商业SAC305合金（Sn-3.0Ag-0.5Cu），试样在0.85倍熔点温度下热处理3小时后随炉冷却，并在室温下时效超过20天以确保力学性能均匀性。实验在微型电磁力疲劳测试系统上进行，配合数字图像相关技术（DIC）测量应变值，系统研究了应力加载速率、平均应力和应力幅值对SAC305焊料棘轮行为的影响。共设计了9组不同加载条件的实验（应力速率2~200 MPa/s、平均应力8~22 MPa、应力幅值6~22 MPa），涵盖了拉-拉和拉-压两种棘轮加载模式。基于统一蠕变-塑性本构框架，建立了考虑屈服面半径为零的率相关本构模型，并引入连续损伤力学中的非线性损伤演化律，通过弹性模量退化定义损伤变量D，最终利用MATLAB拟合实验数据确定了模型参数。

图1 试样几何形状（单位：mm）

本构模型框架基于von Mises屈服准则，由于SAC305焊料无明显弹性阶段，其屈服面半径定义为0。总应变率分解为弹性应变率和黏塑性应变率两部分，其中弹性部分遵循Hooke定律，黏塑性部分由流动法则定义。背应力分解为短程和长程两部分，分别采用不同的演化方程描述。引入损伤变量D后，有效应力概念被应用于修正屈服函数、弹性本构关系、流动法则和背应力演化方程，构成完整的损伤耦合本构模型。损伤变量的定义基于卸载弹性模量的退化：D = (E₀ - E_eff)/E₀，其中E₀为初始弹性模量，E_eff为各循环的有效弹性模量。损伤演化律由连续损伤力学的损伤耗散势推导，并结合归一化寿命进行修正以简化损伤演化方程。

图2 SAC305焊料在不同应力加载速率下的棘轮响应（σₐ = 8 MPa，σₘ = 20 MPa，σ̇ = 200/20/2 MPa/s）

图3 不同平均应力下的棘轮响应（σ̇ = 20 MPa/s，σₐ = 8 MPa，σₘ = 22/20/18 MPa）

图4 不同应力幅值下的棘轮响应：(a) σₘ = 20 MPa, σₐ = 10/8/6 MPa；(b) σₘ = 8 MPa, σₐ = 22/20 MPa

主要发现如下：

1. 应力加载速率对棘轮行为的影响

SAC305焊料的棘轮响应对应力加载速率具有显著的率相关性。随着应力加载速率增大（2→20→200 MPa/s），棘轮应变累积速率逐渐减小，疲劳寿命显著增加（60→1336→29647周次）。棘轮应变累积曲线可分为三个阶段：初始瞬态累积阶段、棘轮应变率稳态阶段和棘轮应变率加速阶段。在较高加载速率下可观察到明显的初始瞬态阶段，而低速率下该阶段不明显。SAC305焊料的大部分寿命处于棘轮应变率稳态阶段。在加速阶段，微裂纹的产生促进了棘轮应变的持续累积，最终导致断裂失效。率相关性与位错微结构、晶粒尺寸和界面连接强度密切相关：较高加载速率意味着更快的位错积累速率，阻碍进一步塑性变形；较低加载速率下位错积累密度小，材料有一定程度的回复和再结晶，增大了棘轮应变。

2. 平均应力和应力幅值对棘轮行为的影响

平均应力的增大使棘轮累积的第一阶段更加明显，棘轮累积速率增大，经过较少循环周次后即发生断裂失效。应力幅值对拉-拉和拉-压两种棘轮模式的响应均很敏感。随着应力幅值增大，棘轮累积速率增大，棘轮寿命减小。在拉-压模式下（σₘ = 8 MPa，σₐ = 22/20 MPa），尽管平均应力较小，但由于较大的应力幅值导致每个循环中塑性变形较大，棘轮应变仍然显著累积。

3. 基于弹性模量退化的损伤变量定义

本研究创新性地从循环加载过程中卸载段的弹性模量退化中提取损伤变量D。初始损伤值D₀设为0，临界损伤值Dcr取所有加载条件下的平均值0.49，与先前研究中微裂纹开始聚合形成宏观裂纹的Dcr = 0.5相吻合。损伤演化律由连续损伤力学的损伤耗散势推导得到：D = D₀ + (Dcr - D₀){1 - [1 - ln(Nnorm/Nnorm_th)/ln(Nnorm_cr/Nnorm_th)]^α}，其中Nnorm为归一化寿命，Nnorm_th和Nnorm_cr分别为损伤起始和材料失效时的归一化寿命值，α为损伤指数。归一化寿命的引入简化了不同加载条件下损伤演化律的比较，提高了计算效率。

图5 偏应力空间中应力组成示意图：(a) 屈服应力σᵧ ≠ 0；(b) 屈服应力σᵧ = 0

图6 不同加载条件下SAC305焊料损伤变量的实验值与模型拟合曲线（实验编号a~i）

4. 塑性应变能密度修正损伤指数

不同加载条件下损伤指数α存在差异，但与首周循环的塑性应变能密度w之间存在幂律关系：α = A·w^B，其中A = 0.23604，B = 0.21966，回归决定系数R² = 0.93。较大的α意味着每个循环中更大的损伤值和更高的材料失效倾向。该方法通过能量法修正损伤演化参数，提供了一种新的损伤评估方法，能够区分不同加载条件下损伤演化的速度。能量法具有温度无关性和加载路径积分的优势，适合解决热-力疲劳问题。

图7 损伤指数α与塑性应变能密度w的幂律拟合关系

5. 损伤耦合本构模型的全寿命棘轮模拟

损伤耦合本构模型能够准确描述SAC305焊料在不同应力加载速率、平均应力和应力幅值下的全寿命棘轮行为。与不考虑损伤的模型相比，损伤耦合模型可以模拟棘轮应变的稳态上升趋势和加速阶段，而无损伤模型仅表现为"棘轮安定"和小幅棘轮累积，无法描述加速失效阶段。尽管在棘轮应变率加速阶段的预测值略低于实验值，但模型能够准确捕捉从稳态阶段到加速阶段的拐点，这对于焊料寿命预测具有重要意义。由于SAC305焊料的大部分寿命处于稳态阶段，稳态阶段的持续时间被采用作为失效判据。模型同时能够模拟拉-拉和拉-压两种形式的棘轮响应。

图8 损伤耦合模型与无损伤模型的模拟结果对比（σₘ = 20 MPa，σₐ = 8 MPa，σ̇ = 2 MPa/s）

图9 不同应力下SAC305焊料棘轮行为的应力加载速率实验结果与模拟结果对比

图10 不同应力下SAC305焊料棘轮行为的平均应力实验结果与模拟结果对比

图11 不同应力下SAC305焊料棘轮行为的应力幅值实验结果与模拟结果对比

结论与工程启示

该研究建立了损伤耦合本构模型来预测SAC305焊料在室温下的全寿命棘轮行为，通过不同加载速率、平均应力和应力幅值的系统实验验证了模型的有效性。主要结论包括：

（1）损伤变量通过各循环卸载弹性模量的演化来表征。损伤演化律经归一化寿命修正，并引入塑性应变能密度项，该改进损伤模型能够模拟不同加载条件下损伤变量的演化规律。

（2）不同加载条件下损伤演化的速度可通过首周循环的塑性应变能密度来区分，为损伤评估提供了一种新的能量方法。该方法具有温度无关性和加载路径积分的优势。

（3）SAC305焊料的棘轮累积与应力加载速率、平均应力和应力幅值密切相关。损伤耦合本构模型能够成功捕捉焊料棘轮行为对应力加载速率、平均应力和应力幅值变化的敏感性，模拟结果与实验结果吻合良好。

工程启示方面，该研究提出的基于弹性模量退化提取损伤变量并结合塑性应变能密度修正的方法，为电子封装焊料的可靠性评估提供了一条实用路径。稳态阶段到加速阶段的拐点作为寿命预测判据，具有较高的工程应用价值。然而，该模型目前仅适用于单轴加载条件，对多轴加载的适用性有待进一步验证。此外，能量法的温度无关性特征使其有潜力拓展至热-力耦合疲劳问题的分析，这对于电子封装焊料在复杂温度场下的寿命评估具有重要意义。该方法框架也可推广至其他材料的损伤演化研究，为蠕变-棘轮-疲劳交互作用的探索提供了有效手段。