研究概述

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），涵蓋了拉-拉和拉-壓兩種棘輪加載模式?；诮y一蠕變-塑性本構框架，建立了考慮屈服面半徑為零的率相關本構模型，并引入連續損傷力學中的非線性損傷演化律，通過彈性模量退化定義損傷變量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，回歸決定系數R2 = 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焊料的棘輪累積與應力加載速率、平均應力和應力幅值密切相關。損傷耦合本構模型能夠成功捕捉焊料棘輪行為對應力加載速率、平均應力和應力幅值變化的敏感性，模擬結果與實驗結果吻合良好。

工程啟示方面，該研究提出的基于彈性模量退化提取損傷變量并結合塑性應變能密度修正的方法，為電子封裝焊料的可靠性評估提供了一條實用路徑。穩態階段到加速階段的拐點作為壽命預測判據，具有較高的工程應用價值。然而，該模型目前僅適用于單軸加載條件，對多軸加載的適用性有待進一步驗證。此外，能量法的溫度無關性特征使其有潛力拓展至熱-力耦合疲勞問題的分析，這對于電子封裝焊料在復雜溫度場下的壽命評估具有重要意義。該方法框架也可推廣至其他材料的損傷演化研究，為蠕變-棘輪-疲勞交互作用的探索提供了有效手段。