發布日期：2025-12-20 15:34:23
引言
鈦合金以其卓越的高比強度、優異的耐腐蝕性、出色的低溫性能及良好的焊接性能，在深海環境結構中占據了舉足輕重的地位 [1]，特別是Ti6Al4V ELI鈦合金展現了其獨特的材料優勢，廣泛應用于載人潛水器的球艙制造中 [2]。然而，深海服役的鈦合金結構需同時承受巨大的靜水壓力和復雜的交變載荷，這些因素極易誘發微小裂紋的萌生并推動其擴展，對深海結構的長期穩定性和安全性構成了嚴峻挑戰。在鈦合金結構的疲勞壽命中，裂紋萌生階段占據了絕大部分 [3]。因此，深入探究Ti6Al4V ELI鈦合金構件的疲勞裂紋萌生壽命，對于精確評估服役環境下鈦合金結構的使用壽命，以及科學制定其安全檢查周期具有重要意義。
目前，針對Ti6Al4V ELI鈦合金疲勞壽命分析的方法主要分為兩類：①基于疲勞累積損傷理論的疲勞安全壽命計算方法；②基于線彈性斷裂力學的損傷容限計算方法。疲勞安全壽命計算方法利用簡單試件的應力 - 壽命曲線進行構件壽命估算。高焱等 [4] 基于機械工程研究協會疲勞強度設計規范 FKM (Forschungs Kuratorium Maschinenbau) 得到了適用于增材制造 TA15 鈦合金結構的疲勞壽命計算方法，且計算精度較好。LIN 等 [5] 的研究表明，與應力 - 壽命評估方法相比，應變 - 壽命方法能夠更準確地描述結構在周期性過載下的局部塑性變形。損傷容限計算方法是基于線彈性斷裂力學建立的，該方法對于穩定擴展階段的長裂紋擴展壽命預測比較準確。吳連生等 [6]、ZHAO 等 [7]、BIAN 等 [8]、YU 等 [9-10]、陶蒙等 [11] 在三維斷裂理論框架下，實現了不同厚度Ti6Al4V ELI鈦合金試件在不同載荷下的疲勞裂紋擴展壽命和裂紋長度的準確預測，提高了深海耐壓結構的損傷容限設計精度。張海英等 [12] 通過試驗獲取了不同狀態激光選區熔化 Ti6Al4V 鈦合金的裂紋擴展性能數據，證明了超過一定范圍的裂紋尺寸和應力強度因子會導致表面機械加工處理對裂紋擴展性能幾乎無影響。WANG 等 [13] 對深海載人潛水器中使用的Ti6Al4V ELI鈦合金進行了結構級駐留疲勞試驗，驗證了載人潛水器壓力艙體在駐留疲勞載荷下的結構完整性。然而，當前針對Ti6Al4V ELI鈦合金疲勞壽命預測的研究仍存在一定的局限性。疲勞安全壽命計算易受結構幾何參數、連接方式、工藝方法及加載順序等多重因素的影響，導致預測結果出現偏差。而損傷容限計算方法則無法涵蓋裂紋萌生階段的壽命預測。CHAN [14] 通過大量試驗和數值模擬發現，材料微結構對于疲勞裂紋的萌生和擴展起著重要的作用，如晶粒尺寸的均勻性很大程度上會影響材料的疲勞壽命。Ti6Al4V ELI 鈦合金結構的裂紋萌生壽命深受材料微觀結構特征和局部應力狀態的共同影響，它們在結構微缺陷和微裂紋的損傷及演化過程中起重要作用，使得在宏觀特性和載荷水平相同的情況下，不同微結構試件的疲勞壽命呈現出顯著的差異性。因此，需要采用晶體塑性 (Crystal Plasticity, CP) 有限元模型來模擬微觀結構對材料疲勞壽命的影響。在晶體塑性理論中，晶體材料的塑性變形主要源自滑移面上的位錯滑移。通過預設晶粒尺寸的成核與生長關系，在不均勻微觀結構中對有限小裂紋的擴展進行模擬，從而可以較為準確地估計初始裂紋的分布情況。這種方法不僅能夠精確預測裂紋的萌生過程，還能夠深入揭示材料疲勞破壞的微觀機制 [15]。
為了準確預測裂紋萌生壽命，需要引入疲勞指示因子 (Fatigue Indicator Parameter, FIP)，并將其與疲勞裂紋萌生相關的微觀力學行為和驅動力緊密聯系起來。自 MANONUKUL 等 [16] 在 2004 年提出晶體塑性有限元框架下的 FIP 概念以來，該領域取得了快速發展，形成了多種基于累積塑性滑移 (Plastic Strain Accumulation, PSA)、應變能耗散 [17]、應變統計學理論 [18] 及駐留滑移帶能量 [19] 的疲勞指示因子。其中，累積塑性滑移 [20-23] 和累積能量耗散 [24] 是 2 種應用較多的疲勞指示因子。MCDOWELL 等 [25-27] 引入 FIP，將細觀力學與基于微觀結構的驅動力結合，描述晶粒疲勞裂紋萌生行為。LIU 等 [28] 將對 TC4 合金的試驗與晶體塑性理論相結合，得出柱面滑移處高度累積的塑性應變和應力集中導致裂紋萌生，裂紋通常在界面處成核的結論。然而，現有的預測模型大多聚焦于均勻晶粒結構，忽視了真實微觀結構的復雜性，這導致預測結果與實際情況之間存在顯著差異。因此，有必要通過電子背散射衍射 (Electron Back Scatter Diffraction, EBSD) 技術重構Ti6Al4V ELI鈦合金的真實微觀結構，并在此基礎上建立晶體塑性有限元模型。
針對具有雙相組織的Ti6Al4V ELI鈦合金，運用 EBSD 技術構建了一個基于其真實微觀結構的晶體塑性循環本構模型。為了更精確地評估微觀短裂紋的成核壽命，模型采用了累積塑性滑移和累積能量耗散 2 種疲勞指示因子。首先，本文對光滑試樣開展了單軸拉伸試驗和疲勞試驗。其次，將試驗結果與有限元模擬結果進行比較，確定了晶體塑性模型中的關鍵參數，包括初始硬化模量和臨界分解切應力等，這些參數的精確獲取為后續的疲勞壽命預測奠定了堅實的基礎。隨后，利用疲勞指示因子與循環次數之間存在的近似線性關系，對不同微結構Ti6Al4V ELI鈦合金試樣在不同應變幅下的疲勞裂紋萌生壽命進行了預測。最后，將預測結果與經典的 Coffin-Manson (C-M) 定律進行了對比分析。強調了考慮材料微觀結構在疲勞壽命預測中的重要性，建立的晶體塑性循環本構模型在揭示材料疲勞行為微觀機制方面展現了獨特優勢。這一研究旨在更深入地探索Ti6Al4V ELI鈦合金的疲勞裂紋萌生壽命，從而為其耐久性和結構安全性提供有力支持。
1、試驗設計
本試驗選用具有高損傷容限性能的Ti6Al4V ELI鈦合金。為了獲取Ti6Al4V ELI鈦合金的真實微觀結構和力學參數，采用 EBSD 技術對在不同應變加載條件下的試件進行微觀結構表征，并建立與其對應的晶體塑性有限元模型。根據樣品表面質量，選擇的掃描波束電壓為 20 kV，工作距離為 15 mm，傾斜角度為 70°，總掃描時間為 45 min，放大倍率為 1000。經過 EBSD 測試，獲得了樣品的顯微組織信息。圖 1 (a)~ 圖 1 (c) 分別為晶粒的取向分布圖、尺寸分布圖及極圖。
使用圖 2 (a) 所示的 MTS Landmark 370.10 型試驗機對Ti6Al4V ELI鈦合金進行準靜態單軸拉伸試驗和疲勞裂紋萌生試驗。Ti6Al4V ELI 鈦合金的元素組成如表 1 所示。在單軸拉伸試驗中，加載速度設定為 3 mm/min，并使用常溫引伸計來精確測量試樣在拉伸過程中所經歷的變形。為確保結果的可靠性，試驗選取 3 塊試樣進行重復拉伸試驗。同時，針對含有預設缺口的試樣，進行了疲勞裂紋萌生試驗。該試驗中，加載頻率設定為 8 Hz，載荷可用三角波形加載幅值曲線描述，并確保應變比為 0。為了深入探究不同應變幅值對Ti6Al4V ELI鈦合金疲勞裂紋萌生壽命的具體影響，本文設計了多組疲勞裂紋萌生試驗，各組的應變 (ε) 分別為 0.5%、0.55%、0.6%、0.65%、0.7% 及 0.75%。
依照標準 GB/T 228.1-2021 [29]，設計了圖 2 (b) 和圖 2 (c) 所示的板狀結構試樣。這些試樣通過精密的線切割技術加工而成，具體分為單邊缺口拉伸 (Single Edge Notch Tension, SENT) 試樣和無缺口試樣 2 種類型。圖 2 (b) 為含缺口試樣，用于評估Ti6Al4V ELI鈦合金在不同應變條件下的裂紋萌生壽命。圖 2 (c) 為無缺口試樣，其主要用于校準和驗證仿真模型的參數。在試驗準備階段，首先在試件表面均勻噴涂一層白色油漆，以增強裂紋的可視性。隨后，使用菲林尺配合數字顯微鏡對缺口部位進行細致觀察，這一步驟對于準確監測裂紋的起裂位置和測量裂紋長度至關重要。為了精確獲取含缺口試樣的裂紋萌生壽命，試驗采取定期中斷的策略，即每當達到一定的循環次數后，試驗會暫停，并施加一個相當于 80% 最大水平應變的靜態拉伸載荷。這一操作有助于清晰地暴露出潛在的裂紋。接著，利用數碼顯微鏡仔細觀察缺口處的裂紋形態，并通過菲林尺精確測量裂紋的長度，其測量范圍設定為 0~0.3 mm。與此同時，詳細記錄當前應變每個裂紋長度所對應的循環次數，從而為后續的裂紋萌生壽命分析提供準確的數據。
表 1Ti6Al4V ELI鈦合金的元素組成
Tab. 1 Elemental composition of theTi6Al4V ELItitanium alloy（%）
| Ti | Al | V | Fe | O | N | C | H |
| 余量 Balance | 6.2 | 4.2 | 0.046 | 0.12 | 0.008 | 0.0013 | 0.0045 |
2、晶體塑性理論
本文采用率相關晶體塑性本構模型，其本構方程主要描述了塑性變形過程中剪切速率與剪切應力之間的關系。滑移系統的變形是由分解切應力驅動的。晶體材料的總變形梯度 F 的表達式 [30] 為：
式中， F e 為由于晶格點陣的拉伸和扭曲而產生的彈性變形梯度； F p 為由于晶體滑移而引起的塑性變形梯度。塑性速度梯度 L p 與滑移系 α 上的滑移率 γ ˙ α 的關系式為：
式中， s α 為滑移方向； m α 為滑移面法向。塑性變形是由各個滑移系下的剪切變形貢獻得到的。根據 Schmid 定律，當分解切應力的大小超出臨界值 (臨界分解切應力) 時，滑移系統將被激活。晶體滑移系上分解切應力 τ α 的表達式為：
式中，σ 為 Cauchy 應力張量。在晶體發生塑性變形時，滑移系開動以及滑移時產生的滑移剪切速率 γ ˙ α 可根據黏塑性冪法則計算 [31]，即：
式中， γ ˙ 0 為參考塑性滑移率； χ α 為背應力； g α 為某一時刻滑移系的應變硬化強度；n 為速率敏感指數，當 n 趨于無窮時，上述方程是率無關的。參數 g α 可表征材料的應變硬化過程，其變化率 g ˙ α 的表達式為：
式中， h αβ 為自硬化模量 (α=β) 或潛硬化模量 (α≠β)； γ ˙ β 為滑移系 β 上的滑移率。潛硬化模量可表示為：
式中，q 為常數； h 0 、 τ 0 和 g ∞ 分別為初始硬化模量、臨界分解切應力和飽和應力，且均為材料常數； γ α 為滑移剪切應變。
γ 為所有滑移系上的累積剪切應變，即滑移率 γ ˙ 隨時間 t 的積分，可表示為：
式 (4) 中，運動硬化的背應力χα遵循非線性演化規律，即：
式中，C、D 分別為描述滑移系線性、非線性硬化特性的材料常數。累積塑性滑移已被證明與循環變形過程中的疲勞損傷密切相關，因此常作為疲勞指示因子，用于預測低周疲勞裂紋萌生壽命。累積塑性滑移 [32] 為：
其中，
當累積塑性滑移達到臨界值pcrit時，疲勞裂紋萌生。
pcrit滿足的關系式為：
式中，Ni為從試驗中獲得的低周疲勞 (Low Cycle Fatigue, LCF) 裂紋長度 0~0.3 mm 對應的循環次數，近似等于 LCF 裂紋萌生壽命；pcycj為每循環下的預測塑性應變增量。
此外，累積能量耗散 W 是評估疲勞裂紋萌生的另一種疲勞指示因子，表示為：
與累積塑性滑移類似，累積能量耗散臨界值Wcrit滿足的關系式為：
3、有限元仿真
基于真實微觀結構的模型能夠精確地再現Ti6Al4V ELI鈦合金的實際微觀結構特征、形貌以及晶粒間的不均勻性，這對于準確描述其低周疲勞損傷行為至關重要。為了在保證模擬精度的同時降低計算復雜度，模型采用 EBSD 數據后處理軟件 Aztec Crystal，從圖 1 (a) 所示的真實雙相微觀結構圖像中裁剪出一個尺寸為 9.6 μm×8.0 μm 的代表性區域，如圖 3 (a) 所示。隨后，利用 Matlab 軟件中的 MTEX 工具箱，基于該裁剪后的區域構建一個二維的代表性體積元 (Representative Volume Element, RVE) 模型，用于深入研究循環加載下材料的變形行為。圖 3 (b) 為采用 C3D8 單元對 RVE 幾何模型進行細致的離散化處理。根據 EBSD 的結果，將 RVE 中的晶粒G1和G2指定為 β 相，而其余的晶粒則統一歸類為 α 相，并使用 Python 腳本為所有晶粒隨機分配了取向。由于 RVE 模型在宏觀尺度上可以視為一個物質點，因此為其施加了周期性邊界條件，并在 D 點施加循環位移載荷，如圖 3 (b) 所示。這樣設計是為了模擬材料在真實環境中的受力情況。圖 3 (c) 為所施加的三角波循環載荷，其中應變比均設為 0，以確保試驗條件的一致性。
為了精確獲取Ti6Al4V ELI鈦合金在晶體塑性模擬中所需的參數，對試件進行了一系列單軸拉伸試驗以及不同應變 (Δεt) 的低周疲勞試驗。圖 4 為試驗結果與仿真結果的對比。由圖 4 可知，無論是單軸拉伸還是循環加載，有限元仿真得到的應力 - 應變曲線都與試驗結果保持了高度的一致性。這里所提及的應力與應變為真應力、真應變。利用試驗中獲得的名義應力和名義應變，通過方程σtrue=σe(1+εe)可得到相應的真應力。試驗和仿真結果的一致性不僅彰顯了仿真的精確性，也充分驗證了所選用材料參數的有效性和可靠性。
從基于試驗所得的應力 - 應變曲線中，可得出一系列關鍵參數，包括初始硬化模量h0、飽和應力g∞、臨界分解切應力τ0以及與背應力相關的 C、D 參數，部分參數值如表 2 所示。對于 α 相的晶體塑性參數，參考了 BALASUBRAMANIAN 等 [33] 的研究成果；而對于 β 相的參數，則參考了 OGI 等 [34] 和 HE 等 [35] 的研究數據。在仿真過程中，設定試樣的參考塑性滑移率γ˙0為0.001s−1，滑移自硬化與潛硬化效應的比值 q 為 1.0，與背應力相關的 C、D 參數分別為 500、50。
表 2 α 相和 β 相的材料常數
Tab. 2 Material constants of α-phase and β-phase
| 滑移系類型 | 晶面 | h 0 /MPa | τ 0 /MPa | g ∞ /MPa |
| 基面 Basal | {0001} | 605.0 | 390.0 | 400.2 |
| 柱面 Prismatica | {10 1 ˉ0} | 640.0 | 391.0 | 400.2 |
| 錐面 Pyramidal | {10 1 ˉ1} | 645.0 | 391.0 | 400.2 |
| 錐面 |
{10 1 ˉ1} | 650.0 | 391.0 | 400.2 |
| β 相 β-phase | {110} | 150.0 | 80.3 | 99.5 |
表 3Ti6Al4V ELI鈦合金的晶格滑移系
Tab. 3 Lattice slip systems of theTi6Al4V ELItitanium alloy
| 滑移系統 Slip system | 晶面 - 晶向 | 數量 Number |
| 基面 Basal | {0001}<11 2 ˉ0> | 3 |
| 柱面 Prismatica | {10 1 ˉ0}<11 2 ˉ0> | 3 |
| 錐面 Pyramidal | {10 1 ˉ1}<11 2 ˉ0> | 6 |
| 錐面 |
{10 1 ˉ1}<11 2 ˉ3> | 12 |
| β 相 β-phase | {110}<111> | 12 |
鑒于疲勞試驗采用的是帶有缺口的試件，并且 EBSD 技術測量的是缺口部位的微觀結構，因此必須考慮缺口導致的應變集中對 RVE 模型可能產生的影響。為了精確評估這一影響，結合試驗數據構建了一個宏觀有限元模型。該模型首先模擬了不同應變條件下，缺口拉伸方向上的最大應變值；隨后，將這些模擬得到的最大應變值作為輸入條件，應用于 RVE 模型中。為了確保模擬過程能夠準確描述材料的力學性能，針對Ti6Al4V ELI鈦合金的特性，本研究選擇了最為合適的本構模型。通過分析Ti6Al4V ELI鈦合金的應力 - 應變曲線，確定了其彈性模量為 121.2 GPa，并采用了指數形式的數學本構模型對其進行描述。該模型的具體表達式為 σ = A ′ + B ′ε z，其中，參數 A ′ = 242.34MPa、 B ′ = 1433.69MPa和 z = 0.128分別為特定溫度條件下材料的屈服強度、應變硬化系數以及應變硬化指數。在宏觀有限元模型的設置中，對模型的下部在 x、y、z 這 3 個方向均施加了固定約束，對上部則施加了與試驗條件相匹配的單軸拉伸應變。模擬的結果如圖 5 所示，對于應變值為 0.5%、0.55%、0.6%、0.65%、0.7% 和 0.75% 的試件，其缺口拉伸方向上的最大應變分別達到了 1.26%、1.48%、1.74%、2.02%、2.32% 和 2.64%。
4、模擬結果與討論
4.1 疲勞損傷
通過 EBSD 試驗，構建了反映真實微觀結構的 RVE 模型，并利用 Abaqus 有限元分析軟件對其進行了循環變形模擬，旨在從微觀層面深入探究疲勞損傷的相關因素。圖 6 (a) 和圖 6 (c) 分別為在應變幅為 0.63% 的條件下，經過 10 周和 100 周循環后 RVE 中的累積塑性滑移分布情況。由圖 6 (a)、圖 6 (c) 可知，與宏觀上的均勻應變狀態相比，微觀層面上 RVE 內各晶粒的應變呈現出顯著的不均勻性，特別是 β 相晶粒更傾向于發生塑性變形。沿路徑 PP ′的累積塑性滑移分布進一步揭示了應變主要集中在 α 相與 β 相的界面處，如圖 6 (e) 所示。圖 6 (e) 中黑色虛線標記的位置為累積塑性滑移的峰值，這暗示了相鄰晶粒間滑移阻力的差異是上述現象出現的關鍵因素。隨著循環次數的累積，局部應變呈現出明顯的增強趨勢，如圖 6 (c) 所示。此外，圖 6 (b) 和圖 6 (d) 分別呈現了相同應變幅下，第 10 周和第 100 周循環后的累積能量耗散分布情況。值得注意的是，累積能量耗散的分布與累積塑性滑移的分布高度相似。在給定應變幅下，累積能量耗散隨著循環次數的增加而持續攀升。沿路徑 PP ′的累積能量耗散分布趨勢與累積塑性滑移分布趨勢基本一致，如圖 6 (f) 所示，這進一步驗證了兩者之間的緊密聯系。鑒于疲勞裂紋通常傾向于先在應變較高的區域萌生，累積塑性滑移值因此成為評估疲勞損傷的有效指標，能夠準確預示裂紋萌生的趨勢。圖 6 中顯示的 α 相與 β 相界面以及 β 相內部較高的累積塑性滑移值預示著這些區域將承受更大的疲勞損傷風險。
已有研究表明，局部應變往往與微觀結構中特定區域的晶格旋轉角密切相關 [36]。這種晶格旋轉角可通過計算旋轉矩陣來確定，具體表達式為θ=π180∘arccos21[trRe−1]（Re為晶格旋轉矩陣），它反映了晶格結構的局部變化。圖 7 (a) 和圖 7 (b) 分別為在應變幅為 0.63% 的條件下，經過 10 周和 100 周循環后 RVE 模型中的晶格旋轉角分布情況。通過對比分析，發現局部疲勞指示因子與局部晶格旋轉之間存在緊密關聯：晶格旋轉角較大的區域局部疲勞指示因子值也相對較高。晶格的局部旋轉不僅導致了微觀紋理的形成，還顯著提升了裂紋萌生位置的 FIP 值，且在 α 相與 β 相的界面處表現得尤為明顯。這一現象清晰地表明，晶格旋轉是促使疲勞裂紋萌生以及小裂紋進一步擴展的關鍵因素之一，這為從微觀角度深入理解和預測材料的疲勞行為提供了新的視角和依據。
為了更深入地探究疲勞裂紋萌生與局部晶格旋轉之間的內在聯系，對圖 7 (a) 標示的P1、P2、P3和P4這 4 個關鍵位置的晶格旋轉角隨時間的變化情況進行進一步分析。圖 8 比較了上述 4 個點在連續 30 個拉伸與壓縮循環周期中的晶格旋轉角變化。由圖 8 可知，從第 1 個周期至第 30 個周期，晶格在拉伸與壓縮過程中的旋轉表現出顯著的差異，呈現出明顯的拉壓不對稱性。這種不對稱性不僅加速了微觀紋理的形成，還進一步加劇了晶格旋轉的不均勻性，從而導致了局部區域累積塑性滑移和累積能量耗散的顯著增加。換言之，這種拉壓不對稱的晶格旋轉行為是促使局部疲勞損傷加劇的關鍵因素之一。值得注意的是，每個循環周期內晶格旋轉的不對稱性還可能引發材料流動行為的拉壓不對稱性，這將進一步影響材料的整體力學性能和疲勞壽命。
針對不同應變加載條件下的試件，進行了系統的微觀結構表征，并結合有限元分析技術開展了深入的模擬研究。圖 9 為在不同應變幅下，經過 80 周循環后的累積塑性滑移分布情況。由圖 9 可知，隨著應變幅的不斷增加，累積塑性滑移量大致呈上升趨勢。這表明疲勞裂紋萌生的風險也在逐步增加，進而導致材料的疲勞壽命相應縮短。此外，圖 10 為在不同應變幅下，經過 80 周循環后的累積能量耗散分布情況。由圖 10 可知，累積能量耗散的分布與累積塑性滑移的分布高度一致，從而可將累積能量耗散視為評估疲勞損傷的又一重要指標參數，它能夠有效地預測與累積塑性滑移相對應的疲勞裂紋萌生位置。模擬研究的結果進一步證實，無論是基于應變的疲勞指標還是基于能量的疲勞指標，在預測疲勞裂紋萌生位置時均展現出較高的準確性和可靠性。
4.2 疲勞壽命的預測
研究局部疲勞指示因子與循環次數之間的動態關系有助于預測Ti6Al4V ELI鈦合金的裂紋萌生壽命。圖 11 為累積塑性滑移和累積能量耗散隨循環次數遞增的變化趨勢。由圖 11 可知，循環次數與這兩者的最大值之間近似成正相關。這一發現為確定如參考塑性滑移率等關鍵材料參數提供了重要依據。圖 11 中的數據均代表整個代表性體積元疲勞指示因子的最大值，這是因為在累積塑性滑移和累積能量耗散達到峰值的區域，材料損傷的風險顯著增加。為了更精確地量化這一關系，在應變幅為 0.87% 的條件下，利用試驗壽命數據 (Ni) 和模擬至第 40 周末的累積塑性滑移增量 (pcyc)，通過式 (11) 得出累積塑性滑移的臨界值 (pcrit=311)；同樣的，利用試驗壽命數據 (Ni) 和模擬至第 40 周末的累積能量耗散增量 (Wcyc)，依據式 (13) 確定累積能量耗散的臨界值 (Wcrit=64374MJ/m3)。
根據試驗結果，將含缺口試件在不同應變 (0.5%、0.55%、0.6%、0.65%、0.7% 和 0.75%) 下的裂紋萌生壽命進行了詳細記錄，其對應的循環次數依次為 6029、5415、3920、3607、3263、2093，如圖 12 (a)~ 圖 12 (f) 所示。此外，為了更深入地理解疲勞損傷機制，并驗證有限元模擬的準確性，本文基于 2 個疲勞指示因子及其臨界值對模擬壽命進行了預測，并對試驗與模擬 2 種預測方法進行了全面評估。圖 13 對比了基于累積塑性滑移模型、累積能量耗散模型和疲勞安全壽命計算方法的疲勞壽命預測結果與低周循環疲勞試驗結果。其中，疲勞安全壽命計算依據文獻 [37] 中的Ti6Al4V ELI鈦合金應變 - 壽命關系式Δεt/2=0.0126(2Nf)−0.080+7.565(2Nf)−1.21，通過線性插值法得出不同應變幅Δε/2下的循環壽命Nrc。對比結果顯示，盡管存在一定的誤差，但基于兩種 FIPs 預測的疲勞壽命呈現出相似的趨勢，并且在準確性上高于疲勞安全壽命計算方法。這一發現驗證了圖 9 和圖 10 中累積塑性滑移和累積能量耗散分布的準確性，且符合前述提出的兩者在相同位置發生了局部化增強，并且都隨著循環次數的增加呈近似線性上升趨勢的結論。進一步分析發現，隨著應變幅的增加，基于累積塑性滑移的 FIP 相較于基于累積能量耗散的 FIP 表現出更高的預測準確性。
圖 14 為塑性應變幅值與基于累積塑性滑移和累積能量耗散這 2 種疲勞指示因子預測的疲勞壽命的對比。在對數坐標系下，這一關系呈鮮明的線性特征。在描述低周循環疲勞壽命Ni時，通常采用 Coffin-Manson (C-M) 定律，該定律將塑性應變幅值Δεp/2通過式Δεp/2=εf′(Ni)c與Ni聯系起來，其中εf′為延性系數，c 為疲勞指數。利用 FIP 進行預測的結果揭示了一個重要現象：與宏觀 C-M 模型相比，微觀結構的不同對材料的裂紋萌生壽命產生了顯著影響。這一發現強調了微觀結構在疲勞損傷過程中的重要作用。因此在進行疲勞壽命預測時，必須充分考慮材料的微觀特征。
4.3 在三維模型中的適用性
探究疲勞指示因子在三維模型中的適用性時，本研究深入分析了晶粒尺寸與網格數量對拉伸應力 - 應變曲線的影響。這對于選定一個既能使模擬曲線與試驗數據高度吻合又能確保計算效率的代表性體積元模型至關重要。借助 Neper 建模軟件，構建了體積相同但晶粒數量各異的模型，具體涵蓋 30、50、100、120 及 150 個晶粒，如圖 15 (a) 所示。對比這些多晶體模型的拉伸應力 - 應變曲線，注意到含 30 個與 50 個晶粒的模型的屈服強度與含 100 個以上晶粒的模型存在顯著差異。在彈性階段，所有模型的應力 - 應變曲線近乎重合；然而，步入塑性階段后，100、120 及 150 晶粒模型的曲線展現出高度的相似性，相比之下，30 與 50 晶粒模型的曲線則明顯偏離。基于上述事實，選定 100 晶粒模型，并對其進行了細致的網格敏感性分析。在 Neper 軟件中，分別設置模型的網格數量大約為 6 萬、10 萬、18 萬和 46 萬，得出不同網格數量的多晶體模型拉伸應力 - 應變曲線，如圖 15 (b) 所示。由圖 15 (b) 可知，網格數量約為 10 萬的晶粒模型的拉伸應力 - 應變曲線呈現出良好的一致性。這一結果意味著，在該網格密度下，模型的計算結果具有較高的穩定性，非常適用于后續的疲勞壽命預測。
采用 Neper 建模軟件，依據 EBSD 獲取的晶粒尺寸平均值，構建了三維晶粒隨機生長的代表性體積元模型。此模型遵循正態分布原則，標準差設定為 1.23，其中 10% 的晶粒設定為 β 相，其余 90% 則為 α 相。RVE 的邊長為 1 mm，內含 100 個晶粒，網格數量調整為約 10 萬個，以確保計算效率與精度。為確保計算的連貫性，對 RVE 施加周期性邊界條件，以模擬材料的真實變形過程。在此基礎上，施加了應變比為 0 的循環位移載荷，并使用第 3 節中確定的材料參數，以深入研究三維結構下的疲勞裂紋萌生壽命。
為分析不同應變幅下疲勞因子在三維結構中的適用性，本節沿用第 3 小節的宏觀模擬方法，重新模擬了 0.4% 應變條件下缺口拉伸方向對應的最大應變幅，結果為 0.45%；試驗測得的裂紋萌生壽命則為 11908 周。圖 16 (a)~ 圖 16 (d) 為應變幅值分別為 0.45%、0.63%、0.74% 和 0.87% 時，第 30 周末的累積塑性滑移分布，以表明這些應變幅值下的疲勞損傷情況。同時，圖 17 (a)~ 圖 17 (d) 為相應應變幅值下第 30 周末的累積能量耗散分布。對比分析可知，隨著應變幅值的增加，累積塑性滑移和累積能量耗散均呈顯著的上升趨勢。
圖 18 (a) 和圖 18 (b) 為不同應變幅下 2 種疲勞指示因子隨循環次數變化的規律。由圖 18 可知，在三維結構中，循環次數與疲勞指示因子同樣呈近似線性的關系。在應變幅為 0.63% 的條件下，借助試驗所得的壽命數據和模擬至第 30 周末時的累積塑性滑移增量及累積能量耗散增量，利用式 (11) 和式 (13) 分別計算出，累積塑性滑移的臨界值為 143，累積能量耗散的臨界值為 30136 MJ/m³。為了進一步驗證這些疲勞指示因子的有效性，圖 19 對比了基于這 2 個因子預測的裂紋萌生壽命、采用疲勞安全壽命計算方法得出的預測壽命，以及實際的試驗結果。由圖 19 可知，基于疲勞指示因子的預測結果與試驗結果高度吻合，且預測值相對保守；在預測準確性方面，基于疲勞指示因子的方法優于疲勞安全壽命計算方法。這一發現驗證了累積塑性滑移和累積能量耗散作為疲勞指示因子在三維結構分析中的適用性。
4.4 討論
盡管本研究在預測Ti6Al4V ELI鈦合金的疲勞裂紋萌生壽命方面取得了初步成果，但仍存在諸多值得深入探討和改進的問題，特別是在模型的維度、疲勞指示因子的線性假設以及試驗觀察方法上。本研究主要基于 EBSD 構建了Ti6Al4V ELI鈦合金的疲勞裂紋萌生壽命預測模型，但這些模型主要局限于二維和簡化的三維問題。盡管二維模型在計算效率和模型復雜性方面具有明顯優勢，但它無法全面考慮三維微觀結構對疲勞行為的影響。實際上，材料的疲勞行為是一個復雜的三維問題，涉及晶粒的空間分布、相界面等多種因素。因此，未來的研究應致力于開發真實微觀結構的三維模型，以更精確地模擬實際條件下的疲勞過程。在構建這些三維模型時，如何在計算復雜性和物理真實性之間找到平衡將是一個重大挑戰。對于高度復雜的微觀結構，可能需要借助先進的數值算法和高性能的服務器進行計算，以更準確地捕捉各個晶粒和相界面之間的相互作用。此外，本研究假設累積塑性滑移和累積能量耗散這 2 種疲勞指示因子與循環次數存在近似線性關系。然而，這 2 種疲勞指示因子在高疲勞加載條件下可能表現出非線性特性，因而上述線性假設雖然簡化了模型，但可能導致對疲勞行為的預測不夠精確。因此，未來的研究應考慮引入非線性疲勞指示因子，以更準確地描述材料在疲勞載荷下的真實響應。盡管本研究得出的結果與試驗結果之間表現出良好的吻合性，為模型的有效性提供了有力支持，但將裂紋長度在 0~0.3 mm 假設為試驗壽命對應的裂紋萌生壽命可能存在一定的局限性。為了更真實地反映裂紋萌生壽命，后續研究需要采用更先進的技術，如高分辨率成像技術或實時監測技術，來觀察裂紋的萌生過程并測量所需數據。這些技術將有助于驗證模型的預測能力并優化模型，從而提高疲勞壽命預測的準確性和可靠性。
5、結論
本研究致力于開發一個針對雙相微觀結構Ti6Al4V ELI鈦合金的晶體塑性模型，以深入探究該材料的低周疲勞行為，并精確預測其裂紋萌生壽命。得出如下主要結論：
為雙相Ti6Al4V ELI鈦合金構建了一個基于真實微觀結構的循環 CP 本構模型。此模型不僅考慮了背應力、分解切應力和滑移阻力的動態演變，還結合試驗數據確定了相關的 CP 參數。此外，采用 2 種疲勞指示因子 (累積塑性滑移和累積能量耗散)，以準確表征疲勞損傷程度。
研究顯示，疲勞損傷更易在 α 相與 β 相的界面處萌生，這一現象與相鄰晶粒間滑移阻力的差異密切相關。局部晶格旋轉不僅加速了微觀紋理的形成，還導致了累積塑性滑移和累積能量耗散的局部集中，從而加速了疲勞裂紋的產生。值得注意的是，每個循環中晶格旋轉的不對稱性引發了流動行為中拉伸與壓縮的不對稱性。
隨著應變幅的增大，累積能量耗散和累積塑性滑移均呈顯著上升趨勢，而疲勞裂紋萌生壽命則相應縮短。通過與試驗結果的對比發現，基于累積塑性滑移的疲勞指示因子相較于基于累積能量耗散的疲勞指示因子具有更高的預測準確性。與經典的 C-M 模型相比，基于疲勞指示因子的預測結果更清晰地展示了Ti6Al4V ELI鈦合金因微觀結構差異而導致的裂紋萌生壽命的顯著差異。此外，這 2 種疲勞指示因子在三維結構中的適用性也得到了驗證，其預測結果不僅準確，而且偏于保守，為工程應用提供了可靠的安全裕量。
參考文獻 (References)
[1] CARRION P E, SHAMSAEI N, DANIEWICZ S R, et al. Fatigue behavior of Ti-6Al-4V ELI including mean stress effects [J]. International Journal of Fatigue, 2017, 99: 87-100.
[2] WANG F, CUI W C, PAN B B, et al. Normalised fatigue and fracture properties of candidate titanium alloys used in the pressure hull of deep manned submersibles [J]. Ships and Offshore Structures, 2014, 9 (3): 297-310.
[3] LIU C H, XU X, SUN T Z, et al. Microstructural effects on fatigue crack initiation mechanisms in a near-alpha titanium alloy [J]. Acta Materialia, 2023, 253: 118957.
[4] 高焱，金海波，李遠杰，等。基于 FKM 的增材制造鈦合金結構疲勞壽命分析方法研究 [J]. 機械強度，2024, 46 (4): 984-991. GAO Yan, JIN Haibo, LI Yuanjie, et al. Research on fatigue life analysis method of additive manufacturing titanium alloy structure based on FKM standard [J]. Journal of Mechanical Strength, 2024, 46 (4): 984-991. (In Chinese)
[5] LIN C W, JU C P, LIN J C. A comparison of the fatigue behavior of cast Ti-7.5Mo with c. p. titanium, Ti-6Al-4V and Ti-13Nb-13Zr alloys [J]. Biomaterials, 2005, 26 (16): 2899-2907.
[6] 吳連生，于培師，韋朋余，等。基于三維理論的 TC4ELI 鈦合金疲勞裂紋擴展研究 [J]. 船舶力學，2022, 26 (9): 1354-1362. WU Liansheng, YU Peishi, WEI Pengyu, et al. Fatigue crack growth of TC4ELI titanium alloy based on three-dimensional theory [J]. Journal of Ship Mechanics, 2022, 26 (9): 1354-1362. (In Chinese)
[7] ZHAO Y X, YU P S, WU L S, et al. Corrosion-induced deceleration-to-acceleration of fatigue crack growth for deep-seaTi6Al4V ELItitanium alloy [J]. Engineering Fracture Mechanics, 2023, 281: 109160.
[8] BIAN J K, YU P S, ZHAO Y X, et al. Fatigue crack growth of marine 980 high-strength steel: from standard specimen testing to three-dimensional curved crack life prediction [J]. Ocean Engineering, 2024, 296: 116950.
[9] YU P S, GUO W L. An equivalent thickness conception for prediction of surface fatigue crack growth life and shape evolution [J]. Engineering Fracture Mechanics, 2012, 93: 65-74.
[10] YU P S, GUO W L. An equivalent thickness conception for evaluation of corner and surface fatigue crack closure [J]. Engineering Fracture Mechanics, 2013, 99: 202-213.
[11] 陶蒙，于培師，趙軍華。基于晶體塑性與內聚力模型的 Cu-Ni-Si 合金晶界起裂研究 [J]. 機械強度，2024, 46 (5): 1184-1191. TAO Meng, YU Peishi, ZHAO Junhua. Study on crack initiation at grain boundary of Cu-Ni-Si alloy based on crystal plasticity and cohesive zone model [J]. Journal of Mechanical Strength, 2024, 46 (5): 1184-1191. (In Chinese)
[12] 張海英，李剛，臧偉鋒，等。激光選區熔化 Ti-6Al-4V 鈦合金疲勞裂紋擴展性能研究 [J]. 機械強度，2023, 45 (6): 1355-1360. ZHANG Haiying, LI Gang, ZANG Weifeng, et al. Fatigue crack growth performance of Ti-6Al-4V titanium alloy by selective laser melting [J]. Journal of Mechanical Strength, 2023, 45 (6): 1355-1360. (In Chinese)
[13] WANG L, YE C, SUN C Q, et al. Experimental investigation on compressive dwell fatigue behavior of titanium alloy pressure hull for deep-sea manned submersibles [J]. Ocean Engineering, 2024, 303: 117646.
[14] CHAN K S. Roles of microstructure in fatigue crack initiation [J]. International Journal of Fatigue, 2010, 32 (9): 1428-1447.
[15] LU X C, ZHAO J F, YU C, et al. Cyclic plasticity of an interstitial high-entropy alloy: experiments, crystal plasticity modeling, and simulations [J]. Journal of the Mechanics and Physics of Solids, 2020, 142: 103971.
[16] MANONUKUL A, DUNNE F P E. High- and low-cycle fatigue crack initiation using polycrystal plasticity [J]. Proceedings of the Royal Society of London Series A: Mathematical, Physical and Engineering Sciences, 2004, 460 (2047): 1881-1903.
[17] KORSUNSKY A M, DINI D, DUNNE F P E, et al. Comparative assessment of dissipated energy and other fatigue criteria [J]. International Journal of Fatigue, 2007, 29 (9/10/11): 1990-1995.
[18] ZHANG K S, JU J W, LI Z H, et al. Micromechanics based fatigue life prediction of a polycrystalline metal applying crystal plasticity [J]. Mechanics of Materials, 2015, 85: 16-37.
[19] SANGID M D, MAIER H J, SEHITOGLU H. A physically based fatigue model for prediction of crack initiation from persistent slip bands in polycrystals [J]. Acta Materialia, 2011, 59 (1): 328-341.
[20] YUAN G J, ZHANG X C, CHEN B, et al. Low-cycle fatigue life prediction of a polycrystalline nickel-base superalloy using crystal plasticity modelling approach [J]. Journal of Materials Science & Technology, 2020, 38: 28-38.
[21] SWEENEY C A, O’BRIEN B, DUNNE F P E, et al. Strain-gradient modelling of grain size effects on fatigue of CoCr alloy [J]. Acta Materialia, 2014, 78: 341-353.
[22] YUAN G J, WANG R Z, GONG C Y, et al. Investigations of micronotch effect on small fatigue crack initiation behaviour in nickel-based alloy GH4169: experiments and simulations [J]. International Journal of Fatigue, 2020, 136: 105578.
[23] YUAN G J, WANG R Z, ZHU W B, et al. Experimental and simulated investigations of low cycle fatigue behavior in a nickel-based superalloy with different volume fractions of δ phase [J]. International Journal of Fatigue, 2021, 153: 106411.
[24] CRUZADO A, LUCARINI S, LLORCA J, et al. Microstructure-based fatigue life model of metallic alloys with bilinear Coffin-Manson behavior [J]. International Journal of Fatigue, 2018, 107: 40-48.
[25] MCDOWELL D L, DUNNE F P E. Microstructure-sensitive computational modeling of fatigue crack formation [J]. International Journal of Fatigue, 2010, 32 (9): 1521-1542.
[26] PRZYBYLA C P, MUSINSKI W D, CASTELLUCCIO G M, et al. Microstructure-sensitive HCF and VHCF simulations [J]. International Journal of Fatigue, 2013, 57: 9-27.
[27] MCDOWELL D L. Simulation-based strategies for microstructure-sensitive fatigue modeling [J]. Materials Science and Engineering: A, 2007, 468: 4-14.
[28] LIU W C, HUANG J, LIU J W, et al. Experimental and crystal plasticity modelling study on the crack initiation in micro-texture regions of Ti-6Al-4V during high cycle fatigue tests [J]. International Journal of Fatigue, 2021, 148: 106203.
[29] 全國鋼標準化技術委員會。金屬材料 拉伸試驗 第 1 部分：室溫試驗方法：GB/T 228.1-2021 [S]. 北京：中國標準出版社，2021: 12-15. National Technical Committee of Steel Standardization. Metallic materials: tensile testing: part 1: method of test at room temperature: GB/T 228.1-2021 [S]. Beijing: Standards Press of China, 2021: 12-15. (In Chinese)
[30] LIU L Y, YANG Q S, LIU X, et al. Crystal cracking of grain-gradient aluminum by a combined CPFEM-CZM method [J]. Engineering Fracture Mechanics, 2021, 242: 107507.
[31] MATSUMOTO H, YOSHIDA S, SASE R, et al. Crystallographic rotation of Ti-6Al-4V alloy with a fiber texture component under isothermal forging through experimental and CPFEM analyses [J]. Materialia, 2023, 30: 101839.
[32] ZHOU D W, WANG X W, YANG X Y, et al. A generalized physical-based failure indicator parameter used in crystal plasticity model to predict fatigue life under low cycle fatigue and creep-fatigue loadings [J]. International Journal of Fatigue, 2023, 166: 107290.
[33] BALASUBRAMANIAN S, ANAND L. Plasticity of initially textured hexagonal polycrystals at high homologous temperatures: application to titanium [J]. Acta Materialia, 2002, 50 (1): 133-148.
[34] OGI H, KAI S, LEDBETTER H, et al. Titanium’s high-temperature elastic constants through the hcp-bcc phase transformation [J]. Acta Materialia, 2004, 52 (7): 2075-2080.
[35] HE T, TANG Y X, WU H B, et al. Research and analysis of TC4 titanium alloy cutting based on crystal plasticity theory [J]. Journal of Physics: Conference Series, 2024, 2785 (1): 012131.
[36] ZHANG K, WU X, DAVIES C H J. Effect of microtexture on short crack propagation in two-phase titanium alloys [J]. International Journal of Fatigue, 2017, 104: 206-220.
[37] 王雷，王琨，李艷青，等. TC4ELI 鈦合金低周疲勞性能研究 [J]. 鈦工業進展，2018, 35 (2): 17-21. WANG Lei, WANG Kun, LI Yanqing, et al. Low-cycle fatigue properties of TC4ELI titanium alloy [J]. Titanium Industry Progress, 2018, 35 (2): 17-21. (In Chinese)
（注，原文標題：基于微結構特征的Ti6Al4V_ELI鈦合金疲勞裂紋萌生壽命晶體塑性預測方法）
tag標簽:Ti6Al4V,ELI鈦合金,深海服役,疲勞安全性評估