發(fā)布日期：2025-12-20 15:34:23
引言
鈦合金以其卓越的高比強(qiáng)度、優(yōu)異的耐腐蝕性、出色的低溫性能及良好的焊接性能，在深海環(huán)境結(jié)構(gòu)中占據(jù)了舉足輕重的地位 [1]，特別是Ti6Al4V ELI鈦合金展現(xiàn)了其獨特的材料優(yōu)勢，廣泛應(yīng)用于載人潛水器的球艙制造中 [2]。然而，深海服役的鈦合金結(jié)構(gòu)需同時承受巨大的靜水壓力和復(fù)雜的交變載荷，這些因素極易誘發(fā)微小裂紋的萌生并推動其擴(kuò)展，對深海結(jié)構(gòu)的長期穩(wěn)定性和安全性構(gòu)成了嚴(yán)峻挑戰(zhàn)。在鈦合金結(jié)構(gòu)的疲勞壽命中，裂紋萌生階段占據(jù)了絕大部分 [3]。因此，深入探究Ti6Al4V ELI鈦合金構(gòu)件的疲勞裂紋萌生壽命，對于精確評估服役環(huán)境下鈦合金結(jié)構(gòu)的使用壽命，以及科學(xué)制定其安全檢查周期具有重要意義。
目前，針對Ti6Al4V ELI鈦合金疲勞壽命分析的方法主要分為兩類：①基于疲勞累積損傷理論的疲勞安全壽命計算方法；②基于線彈性斷裂力學(xué)的損傷容限計算方法。疲勞安全壽命計算方法利用簡單試件的應(yīng)力 - 壽命曲線進(jìn)行構(gòu)件壽命估算。高焱等 [4] 基于機(jī)械工程研究協(xié)會疲勞強(qiáng)度設(shè)計規(guī)范 FKM (Forschungs Kuratorium Maschinenbau) 得到了適用于增材制造 TA15 鈦合金結(jié)構(gòu)的疲勞壽命計算方法，且計算精度較好。LIN 等 [5] 的研究表明，與應(yīng)力 - 壽命評估方法相比，應(yīng)變 - 壽命方法能夠更準(zhǔn)確地描述結(jié)構(gòu)在周期性過載下的局部塑性變形。損傷容限計算方法是基于線彈性斷裂力學(xué)建立的，該方法對于穩(wěn)定擴(kuò)展階段的長裂紋擴(kuò)展壽命預(yù)測比較準(zhǔn)確。吳連生等 [6]、ZHAO 等 [7]、BIAN 等 [8]、YU 等 [9-10]、陶蒙等 [11] 在三維斷裂理論框架下，實現(xiàn)了不同厚度Ti6Al4V ELI鈦合金試件在不同載荷下的疲勞裂紋擴(kuò)展壽命和裂紋長度的準(zhǔn)確預(yù)測，提高了深海耐壓結(jié)構(gòu)的損傷容限設(shè)計精度。張海英等 [12] 通過試驗獲取了不同狀態(tài)激光選區(qū)熔化 Ti6Al4V 鈦合金的裂紋擴(kuò)展性能數(shù)據(jù)，證明了超過一定范圍的裂紋尺寸和應(yīng)力強(qiáng)度因子會導(dǎo)致表面機(jī)械加工處理對裂紋擴(kuò)展性能幾乎無影響。WANG 等 [13] 對深海載人潛水器中使用的Ti6Al4V ELI鈦合金進(jìn)行了結(jié)構(gòu)級駐留疲勞試驗，驗證了載人潛水器壓力艙體在駐留疲勞載荷下的結(jié)構(gòu)完整性。然而，當(dāng)前針對Ti6Al4V ELI鈦合金疲勞壽命預(yù)測的研究仍存在一定的局限性。疲勞安全壽命計算易受結(jié)構(gòu)幾何參數(shù)、連接方式、工藝方法及加載順序等多重因素的影響，導(dǎo)致預(yù)測結(jié)果出現(xiàn)偏差。而損傷容限計算方法則無法涵蓋裂紋萌生階段的壽命預(yù)測。CHAN [14] 通過大量試驗和數(shù)值模擬發(fā)現(xiàn)，材料微結(jié)構(gòu)對于疲勞裂紋的萌生和擴(kuò)展起著重要的作用，如晶粒尺寸的均勻性很大程度上會影響材料的疲勞壽命。Ti6Al4V ELI 鈦合金結(jié)構(gòu)的裂紋萌生壽命深受材料微觀結(jié)構(gòu)特征和局部應(yīng)力狀態(tài)的共同影響，它們在結(jié)構(gòu)微缺陷和微裂紋的損傷及演化過程中起重要作用，使得在宏觀特性和載荷水平相同的情況下，不同微結(jié)構(gòu)試件的疲勞壽命呈現(xiàn)出顯著的差異性。因此，需要采用晶體塑性 (Crystal Plasticity, CP) 有限元模型來模擬微觀結(jié)構(gòu)對材料疲勞壽命的影響。在晶體塑性理論中，晶體材料的塑性變形主要源自滑移面上的位錯滑移。通過預(yù)設(shè)晶粒尺寸的成核與生長關(guān)系，在不均勻微觀結(jié)構(gòu)中對有限小裂紋的擴(kuò)展進(jìn)行模擬，從而可以較為準(zhǔn)確地估計初始裂紋的分布情況。這種方法不僅能夠精確預(yù)測裂紋的萌生過程，還能夠深入揭示材料疲勞破壞的微觀機(jī)制 [15]。
為了準(zhǔn)確預(yù)測裂紋萌生壽命，需要引入疲勞指示因子 (Fatigue Indicator Parameter, FIP)，并將其與疲勞裂紋萌生相關(guān)的微觀力學(xué)行為和驅(qū)動力緊密聯(lián)系起來。自 MANONUKUL 等 [16] 在 2004 年提出晶體塑性有限元框架下的 FIP 概念以來，該領(lǐng)域取得了快速發(fā)展，形成了多種基于累積塑性滑移 (Plastic Strain Accumulation, PSA)、應(yīng)變能耗散 [17]、應(yīng)變統(tǒng)計學(xué)理論 [18] 及駐留滑移帶能量 [19] 的疲勞指示因子。其中，累積塑性滑移 [20-23] 和累積能量耗散 [24] 是 2 種應(yīng)用較多的疲勞指示因子。MCDOWELL 等 [25-27] 引入 FIP，將細(xì)觀力學(xué)與基于微觀結(jié)構(gòu)的驅(qū)動力結(jié)合，描述晶粒疲勞裂紋萌生行為。LIU 等 [28] 將對 TC4 合金的試驗與晶體塑性理論相結(jié)合，得出柱面滑移處高度累積的塑性應(yīng)變和應(yīng)力集中導(dǎo)致裂紋萌生，裂紋通常在界面處成核的結(jié)論。然而，現(xiàn)有的預(yù)測模型大多聚焦于均勻晶粒結(jié)構(gòu)，忽視了真實微觀結(jié)構(gòu)的復(fù)雜性，這導(dǎo)致預(yù)測結(jié)果與實際情況之間存在顯著差異。因此，有必要通過電子背散射衍射 (Electron Back Scatter Diffraction, EBSD) 技術(shù)重構(gòu)Ti6Al4V ELI鈦合金的真實微觀結(jié)構(gòu)，并在此基礎(chǔ)上建立晶體塑性有限元模型。
針對具有雙相組織的Ti6Al4V ELI鈦合金，運用 EBSD 技術(shù)構(gòu)建了一個基于其真實微觀結(jié)構(gòu)的晶體塑性循環(huán)本構(gòu)模型。為了更精確地評估微觀短裂紋的成核壽命，模型采用了累積塑性滑移和累積能量耗散 2 種疲勞指示因子。首先，本文對光滑試樣開展了單軸拉伸試驗和疲勞試驗。其次，將試驗結(jié)果與有限元模擬結(jié)果進(jìn)行比較，確定了晶體塑性模型中的關(guān)鍵參數(shù)，包括初始硬化模量和臨界分解切應(yīng)力等，這些參數(shù)的精確獲取為后續(xù)的疲勞壽命預(yù)測奠定了堅實的基礎(chǔ)。隨后，利用疲勞指示因子與循環(huán)次數(shù)之間存在的近似線性關(guān)系，對不同微結(jié)構(gòu)Ti6Al4V ELI鈦合金試樣在不同應(yīng)變幅下的疲勞裂紋萌生壽命進(jìn)行了預(yù)測。最后，將預(yù)測結(jié)果與經(jīng)典的 Coffin-Manson (C-M) 定律進(jìn)行了對比分析。強(qiáng)調(diào)了考慮材料微觀結(jié)構(gòu)在疲勞壽命預(yù)測中的重要性，建立的晶體塑性循環(huán)本構(gòu)模型在揭示材料疲勞行為微觀機(jī)制方面展現(xiàn)了獨特優(yōu)勢。這一研究旨在更深入地探索Ti6Al4V ELI鈦合金的疲勞裂紋萌生壽命，從而為其耐久性和結(jié)構(gòu)安全性提供有力支持。
1、試驗設(shè)計
本試驗選用具有高損傷容限性能的Ti6Al4V ELI鈦合金。為了獲取Ti6Al4V ELI鈦合金的真實微觀結(jié)構(gòu)和力學(xué)參數(shù)，采用 EBSD 技術(shù)對在不同應(yīng)變加載條件下的試件進(jìn)行微觀結(jié)構(gòu)表征，并建立與其對應(yīng)的晶體塑性有限元模型。根據(jù)樣品表面質(zhì)量，選擇的掃描波束電壓為 20 kV，工作距離為 15 mm，傾斜角度為 70°，總掃描時間為 45 min，放大倍率為 1000。經(jīng)過 EBSD 測試，獲得了樣品的顯微組織信息。圖 1 (a)~ 圖 1 (c) 分別為晶粒的取向分布圖、尺寸分布圖及極圖。
使用圖 2 (a) 所示的 MTS Landmark 370.10 型試驗機(jī)對Ti6Al4V ELI鈦合金進(jìn)行準(zhǔn)靜態(tài)單軸拉伸試驗和疲勞裂紋萌生試驗。Ti6Al4V ELI 鈦合金的元素組成如表 1 所示。在單軸拉伸試驗中，加載速度設(shè)定為 3 mm/min，并使用常溫引伸計來精確測量試樣在拉伸過程中所經(jīng)歷的變形。為確保結(jié)果的可靠性，試驗選取 3 塊試樣進(jìn)行重復(fù)拉伸試驗。同時，針對含有預(yù)設(shè)缺口的試樣，進(jìn)行了疲勞裂紋萌生試驗。該試驗中，加載頻率設(shè)定為 8 Hz，載荷可用三角波形加載幅值曲線描述，并確保應(yīng)變比為 0。為了深入探究不同應(yīng)變幅值對Ti6Al4V ELI鈦合金疲勞裂紋萌生壽命的具體影響，本文設(shè)計了多組疲勞裂紋萌生試驗，各組的應(yīng)變 (ε) 分別為 0.5%、0.55%、0.6%、0.65%、0.7% 及 0.75%。
依照標(biāo)準(zhǔn) GB/T 228.1-2021 [29]，設(shè)計了圖 2 (b) 和圖 2 (c) 所示的板狀結(jié)構(gòu)試樣。這些試樣通過精密的線切割技術(shù)加工而成，具體分為單邊缺口拉伸 (Single Edge Notch Tension, SENT) 試樣和無缺口試樣 2 種類型。圖 2 (b) 為含缺口試樣，用于評估Ti6Al4V ELI鈦合金在不同應(yīng)變條件下的裂紋萌生壽命。圖 2 (c) 為無缺口試樣，其主要用于校準(zhǔn)和驗證仿真模型的參數(shù)。在試驗準(zhǔn)備階段，首先在試件表面均勻噴涂一層白色油漆，以增強(qiáng)裂紋的可視性。隨后，使用菲林尺配合數(shù)字顯微鏡對缺口部位進(jìn)行細(xì)致觀察，這一步驟對于準(zhǔn)確監(jiān)測裂紋的起裂位置和測量裂紋長度至關(guān)重要。為了精確獲取含缺口試樣的裂紋萌生壽命，試驗采取定期中斷的策略，即每當(dāng)達(dá)到一定的循環(huán)次數(shù)后，試驗會暫停，并施加一個相當(dāng)于 80% 最大水平應(yīng)變的靜態(tài)拉伸載荷。這一操作有助于清晰地暴露出潛在的裂紋。接著，利用數(shù)碼顯微鏡仔細(xì)觀察缺口處的裂紋形態(tài)，并通過菲林尺精確測量裂紋的長度，其測量范圍設(shè)定為 0~0.3 mm。與此同時，詳細(xì)記錄當(dāng)前應(yīng)變每個裂紋長度所對應(yīng)的循環(huán)次數(shù)，從而為后續(xù)的裂紋萌生壽命分析提供準(zhǔn)確的數(shù)據(jù)。
表 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、晶體塑性理論
本文采用率相關(guān)晶體塑性本構(gòu)模型，其本構(gòu)方程主要描述了塑性變形過程中剪切速率與剪切應(yīng)力之間的關(guān)系。滑移系統(tǒng)的變形是由分解切應(yīng)力驅(qū)動的。晶體材料的總變形梯度 F 的表達(dá)式 [30] 為：
式中， F e 為由于晶格點陣的拉伸和扭曲而產(chǎn)生的彈性變形梯度； F p 為由于晶體滑移而引起的塑性變形梯度。塑性速度梯度 L p 與滑移系 α 上的滑移率 γ ˙ α 的關(guān)系式為：
式中， s α 為滑移方向； m α 為滑移面法向。塑性變形是由各個滑移系下的剪切變形貢獻(xiàn)得到的。根據(jù) Schmid 定律，當(dāng)分解切應(yīng)力的大小超出臨界值 (臨界分解切應(yīng)力) 時，滑移系統(tǒng)將被激活。晶體滑移系上分解切應(yīng)力 τ α 的表達(dá)式為：
式中，σ 為 Cauchy 應(yīng)力張量。在晶體發(fā)生塑性變形時，滑移系開動以及滑移時產(chǎn)生的滑移剪切速率 γ ˙ α 可根據(jù)黏塑性冪法則計算 [31]，即：
式中， γ ˙ 0 為參考塑性滑移率； χ α 為背應(yīng)力； g α 為某一時刻滑移系的應(yīng)變硬化強(qiáng)度；n 為速率敏感指數(shù)，當(dāng) n 趨于無窮時，上述方程是率無關(guān)的。參數(shù) g α 可表征材料的應(yīng)變硬化過程，其變化率 g ˙ α 的表達(dá)式為：
式中， h αβ 為自硬化模量 (α=β) 或潛硬化模量 (α≠β)； γ ˙ β 為滑移系 β 上的滑移率。潛硬化模量可表示為：
式中，q 為常數(shù)； h 0 、 τ 0 和 g ∞ 分別為初始硬化模量、臨界分解切應(yīng)力和飽和應(yīng)力，且均為材料常數(shù)； γ α 為滑移剪切應(yīng)變。
γ 為所有滑移系上的累積剪切應(yīng)變，即滑移率 γ ˙ 隨時間 t 的積分，可表示為：
式 (4) 中，運動硬化的背應(yīng)力χα遵循非線性演化規(guī)律，即：
式中，C、D 分別為描述滑移系線性、非線性硬化特性的材料常數(shù)。累積塑性滑移已被證明與循環(huán)變形過程中的疲勞損傷密切相關(guān)，因此常作為疲勞指示因子，用于預(yù)測低周疲勞裂紋萌生壽命。累積塑性滑移 [32] 為：
其中，
當(dāng)累積塑性滑移達(dá)到臨界值pcrit時，疲勞裂紋萌生。
pcrit滿足的關(guān)系式為：
式中，Ni為從試驗中獲得的低周疲勞 (Low Cycle Fatigue, LCF) 裂紋長度 0~0.3 mm 對應(yīng)的循環(huán)次數(shù)，近似等于 LCF 裂紋萌生壽命；pcycj為每循環(huán)下的預(yù)測塑性應(yīng)變增量。
此外，累積能量耗散 W 是評估疲勞裂紋萌生的另一種疲勞指示因子，表示為：
與累積塑性滑移類似，累積能量耗散臨界值Wcrit滿足的關(guān)系式為：
3、有限元仿真
基于真實微觀結(jié)構(gòu)的模型能夠精確地再現(xiàn)Ti6Al4V ELI鈦合金的實際微觀結(jié)構(gòu)特征、形貌以及晶粒間的不均勻性，這對于準(zhǔn)確描述其低周疲勞損傷行為至關(guān)重要。為了在保證模擬精度的同時降低計算復(fù)雜度，模型采用 EBSD 數(shù)據(jù)后處理軟件 Aztec Crystal，從圖 1 (a) 所示的真實雙相微觀結(jié)構(gòu)圖像中裁剪出一個尺寸為 9.6 μm×8.0 μm 的代表性區(qū)域，如圖 3 (a) 所示。隨后，利用 Matlab 軟件中的 MTEX 工具箱，基于該裁剪后的區(qū)域構(gòu)建一個二維的代表性體積元 (Representative Volume Element, RVE) 模型，用于深入研究循環(huán)加載下材料的變形行為。圖 3 (b) 為采用 C3D8 單元對 RVE 幾何模型進(jìn)行細(xì)致的離散化處理。根據(jù) EBSD 的結(jié)果，將 RVE 中的晶粒G1和G2指定為 β 相，而其余的晶粒則統(tǒng)一歸類為 α 相，并使用 Python 腳本為所有晶粒隨機(jī)分配了取向。由于 RVE 模型在宏觀尺度上可以視為一個物質(zhì)點，因此為其施加了周期性邊界條件，并在 D 點施加循環(huán)位移載荷，如圖 3 (b) 所示。這樣設(shè)計是為了模擬材料在真實環(huán)境中的受力情況。圖 3 (c) 為所施加的三角波循環(huán)載荷，其中應(yīng)變比均設(shè)為 0，以確保試驗條件的一致性。
為了精確獲取Ti6Al4V ELI鈦合金在晶體塑性模擬中所需的參數(shù)，對試件進(jìn)行了一系列單軸拉伸試驗以及不同應(yīng)變 (Δεt) 的低周疲勞試驗。圖 4 為試驗結(jié)果與仿真結(jié)果的對比。由圖 4 可知，無論是單軸拉伸還是循環(huán)加載，有限元仿真得到的應(yīng)力 - 應(yīng)變曲線都與試驗結(jié)果保持了高度的一致性。這里所提及的應(yīng)力與應(yīng)變?yōu)檎鎽?yīng)力、真應(yīng)變。利用試驗中獲得的名義應(yīng)力和名義應(yīng)變，通過方程σtrue=σe(1+εe)可得到相應(yīng)的真應(yīng)力。試驗和仿真結(jié)果的一致性不僅彰顯了仿真的精確性，也充分驗證了所選用材料參數(shù)的有效性和可靠性。
從基于試驗所得的應(yīng)力 - 應(yīng)變曲線中，可得出一系列關(guān)鍵參數(shù)，包括初始硬化模量h0、飽和應(yīng)力g∞、臨界分解切應(yīng)力τ0以及與背應(yīng)力相關(guān)的 C、D 參數(shù)，部分參數(shù)值如表 2 所示。對于 α 相的晶體塑性參數(shù)，參考了 BALASUBRAMANIAN 等 [33] 的研究成果；而對于 β 相的參數(shù)，則參考了 OGI 等 [34] 和 HE 等 [35] 的研究數(shù)據(jù)。在仿真過程中，設(shè)定試樣的參考塑性滑移率γ˙0為0.001s−1，滑移自硬化與潛硬化效應(yīng)的比值 q 為 1.0，與背應(yīng)力相關(guān)的 C、D 參數(shù)分別為 500、50。
表 2 α 相和 β 相的材料常數(shù)
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
| 滑移系統(tǒng) Slip system | 晶面 - 晶向 | 數(shù)量 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 技術(shù)測量的是缺口部位的微觀結(jié)構(gòu)，因此必須考慮缺口導(dǎo)致的應(yīng)變集中對 RVE 模型可能產(chǎn)生的影響。為了精確評估這一影響，結(jié)合試驗數(shù)據(jù)構(gòu)建了一個宏觀有限元模型。該模型首先模擬了不同應(yīng)變條件下，缺口拉伸方向上的最大應(yīng)變值；隨后，將這些模擬得到的最大應(yīng)變值作為輸入條件，應(yīng)用于 RVE 模型中。為了確保模擬過程能夠準(zhǔn)確描述材料的力學(xué)性能，針對Ti6Al4V ELI鈦合金的特性，本研究選擇了最為合適的本構(gòu)模型。通過分析Ti6Al4V ELI鈦合金的應(yīng)力 - 應(yīng)變曲線，確定了其彈性模量為 121.2 GPa，并采用了指數(shù)形式的數(shù)學(xué)本構(gòu)模型對其進(jìn)行描述。該模型的具體表達(dá)式為 σ = A ′ + B ′ε z，其中，參數(shù) A ′ = 242.34MPa、 B ′ = 1433.69MPa和 z = 0.128分別為特定溫度條件下材料的屈服強(qiáng)度、應(yīng)變硬化系數(shù)以及應(yīng)變硬化指數(shù)。在宏觀有限元模型的設(shè)置中，對模型的下部在 x、y、z 這 3 個方向均施加了固定約束，對上部則施加了與試驗條件相匹配的單軸拉伸應(yīng)變。模擬的結(jié)果如圖 5 所示，對于應(yīng)變值為 0.5%、0.55%、0.6%、0.65%、0.7% 和 0.75% 的試件，其缺口拉伸方向上的最大應(yīng)變分別達(dá)到了 1.26%、1.48%、1.74%、2.02%、2.32% 和 2.64%。
4、模擬結(jié)果與討論
4.1 疲勞損傷
通過 EBSD 試驗，構(gòu)建了反映真實微觀結(jié)構(gòu)的 RVE 模型，并利用 Abaqus 有限元分析軟件對其進(jìn)行了循環(huán)變形模擬，旨在從微觀層面深入探究疲勞損傷的相關(guān)因素。圖 6 (a) 和圖 6 (c) 分別為在應(yīng)變幅為 0.63% 的條件下，經(jīng)過 10 周和 100 周循環(huán)后 RVE 中的累積塑性滑移分布情況。由圖 6 (a)、圖 6 (c) 可知，與宏觀上的均勻應(yīng)變狀態(tài)相比，微觀層面上 RVE 內(nèi)各晶粒的應(yīng)變呈現(xiàn)出顯著的不均勻性，特別是 β 相晶粒更傾向于發(fā)生塑性變形。沿路徑 PP ′的累積塑性滑移分布進(jìn)一步揭示了應(yīng)變主要集中在 α 相與 β 相的界面處，如圖 6 (e) 所示。圖 6 (e) 中黑色虛線標(biāo)記的位置為累積塑性滑移的峰值，這暗示了相鄰晶粒間滑移阻力的差異是上述現(xiàn)象出現(xiàn)的關(guān)鍵因素。隨著循環(huán)次數(shù)的累積，局部應(yīng)變呈現(xiàn)出明顯的增強(qiáng)趨勢，如圖 6 (c) 所示。此外，圖 6 (b) 和圖 6 (d) 分別呈現(xiàn)了相同應(yīng)變幅下，第 10 周和第 100 周循環(huán)后的累積能量耗散分布情況。值得注意的是，累積能量耗散的分布與累積塑性滑移的分布高度相似。在給定應(yīng)變幅下，累積能量耗散隨著循環(huán)次數(shù)的增加而持續(xù)攀升。沿路徑 PP ′的累積能量耗散分布趨勢與累積塑性滑移分布趨勢基本一致，如圖 6 (f) 所示，這進(jìn)一步驗證了兩者之間的緊密聯(lián)系。鑒于疲勞裂紋通常傾向于先在應(yīng)變較高的區(qū)域萌生，累積塑性滑移值因此成為評估疲勞損傷的有效指標(biāo)，能夠準(zhǔn)確預(yù)示裂紋萌生的趨勢。圖 6 中顯示的 α 相與 β 相界面以及 β 相內(nèi)部較高的累積塑性滑移值預(yù)示著這些區(qū)域?qū)⒊惺芨蟮钠趽p傷風(fēng)險。
已有研究表明，局部應(yīng)變往往與微觀結(jié)構(gòu)中特定區(qū)域的晶格旋轉(zhuǎn)角密切相關(guān) [36]。這種晶格旋轉(zhuǎn)角可通過計算旋轉(zhuǎn)矩陣來確定，具體表達(dá)式為θ=π180∘arccos21[trRe−1]（Re為晶格旋轉(zhuǎn)矩陣），它反映了晶格結(jié)構(gòu)的局部變化。圖 7 (a) 和圖 7 (b) 分別為在應(yīng)變幅為 0.63% 的條件下，經(jīng)過 10 周和 100 周循環(huán)后 RVE 模型中的晶格旋轉(zhuǎn)角分布情況。通過對比分析，發(fā)現(xiàn)局部疲勞指示因子與局部晶格旋轉(zhuǎn)之間存在緊密關(guān)聯(lián)：晶格旋轉(zhuǎn)角較大的區(qū)域局部疲勞指示因子值也相對較高。晶格的局部旋轉(zhuǎn)不僅導(dǎo)致了微觀紋理的形成，還顯著提升了裂紋萌生位置的 FIP 值，且在 α 相與 β 相的界面處表現(xiàn)得尤為明顯。這一現(xiàn)象清晰地表明，晶格旋轉(zhuǎn)是促使疲勞裂紋萌生以及小裂紋進(jìn)一步擴(kuò)展的關(guān)鍵因素之一，這為從微觀角度深入理解和預(yù)測材料的疲勞行為提供了新的視角和依據(jù)。
為了更深入地探究疲勞裂紋萌生與局部晶格旋轉(zhuǎn)之間的內(nèi)在聯(lián)系，對圖 7 (a) 標(biāo)示的P1、P2、P3和P4這 4 個關(guān)鍵位置的晶格旋轉(zhuǎn)角隨時間的變化情況進(jìn)行進(jìn)一步分析。圖 8 比較了上述 4 個點在連續(xù) 30 個拉伸與壓縮循環(huán)周期中的晶格旋轉(zhuǎn)角變化。由圖 8 可知，從第 1 個周期至第 30 個周期，晶格在拉伸與壓縮過程中的旋轉(zhuǎn)表現(xiàn)出顯著的差異，呈現(xiàn)出明顯的拉壓不對稱性。這種不對稱性不僅加速了微觀紋理的形成，還進(jìn)一步加劇了晶格旋轉(zhuǎn)的不均勻性，從而導(dǎo)致了局部區(qū)域累積塑性滑移和累積能量耗散的顯著增加。換言之，這種拉壓不對稱的晶格旋轉(zhuǎn)行為是促使局部疲勞損傷加劇的關(guān)鍵因素之一。值得注意的是，每個循環(huán)周期內(nèi)晶格旋轉(zhuǎn)的不對稱性還可能引發(fā)材料流動行為的拉壓不對稱性，這將進(jìn)一步影響材料的整體力學(xué)性能和疲勞壽命。
針對不同應(yīng)變加載條件下的試件，進(jìn)行了系統(tǒng)的微觀結(jié)構(gòu)表征，并結(jié)合有限元分析技術(shù)開展了深入的模擬研究。圖 9 為在不同應(yīng)變幅下，經(jīng)過 80 周循環(huán)后的累積塑性滑移分布情況。由圖 9 可知，隨著應(yīng)變幅的不斷增加，累積塑性滑移量大致呈上升趨勢。這表明疲勞裂紋萌生的風(fēng)險也在逐步增加，進(jìn)而導(dǎo)致材料的疲勞壽命相應(yīng)縮短。此外，圖 10 為在不同應(yīng)變幅下，經(jīng)過 80 周循環(huán)后的累積能量耗散分布情況。由圖 10 可知，累積能量耗散的分布與累積塑性滑移的分布高度一致，從而可將累積能量耗散視為評估疲勞損傷的又一重要指標(biāo)參數(shù)，它能夠有效地預(yù)測與累積塑性滑移相對應(yīng)的疲勞裂紋萌生位置。模擬研究的結(jié)果進(jìn)一步證實，無論是基于應(yīng)變的疲勞指標(biāo)還是基于能量的疲勞指標(biāo)，在預(yù)測疲勞裂紋萌生位置時均展現(xiàn)出較高的準(zhǔn)確性和可靠性。
4.2 疲勞壽命的預(yù)測
研究局部疲勞指示因子與循環(huán)次數(shù)之間的動態(tài)關(guān)系有助于預(yù)測Ti6Al4V ELI鈦合金的裂紋萌生壽命。圖 11 為累積塑性滑移和累積能量耗散隨循環(huán)次數(shù)遞增的變化趨勢。由圖 11 可知，循環(huán)次數(shù)與這兩者的最大值之間近似成正相關(guān)。這一發(fā)現(xiàn)為確定如參考塑性滑移率等關(guān)鍵材料參數(shù)提供了重要依據(jù)。圖 11 中的數(shù)據(jù)均代表整個代表性體積元疲勞指示因子的最大值，這是因為在累積塑性滑移和累積能量耗散達(dá)到峰值的區(qū)域，材料損傷的風(fēng)險顯著增加。為了更精確地量化這一關(guān)系，在應(yīng)變幅為 0.87% 的條件下，利用試驗壽命數(shù)據(jù) (Ni) 和模擬至第 40 周末的累積塑性滑移增量 (pcyc)，通過式 (11) 得出累積塑性滑移的臨界值 (pcrit=311)；同樣的，利用試驗壽命數(shù)據(jù) (Ni) 和模擬至第 40 周末的累積能量耗散增量 (Wcyc)，依據(jù)式 (13) 確定累積能量耗散的臨界值 (Wcrit=64374MJ/m3)。
根據(jù)試驗結(jié)果，將含缺口試件在不同應(yīng)變 (0.5%、0.55%、0.6%、0.65%、0.7% 和 0.75%) 下的裂紋萌生壽命進(jìn)行了詳細(xì)記錄，其對應(yīng)的循環(huán)次數(shù)依次為 6029、5415、3920、3607、3263、2093，如圖 12 (a)~ 圖 12 (f) 所示。此外，為了更深入地理解疲勞損傷機(jī)制，并驗證有限元模擬的準(zhǔn)確性，本文基于 2 個疲勞指示因子及其臨界值對模擬壽命進(jìn)行了預(yù)測，并對試驗與模擬 2 種預(yù)測方法進(jìn)行了全面評估。圖 13 對比了基于累積塑性滑移模型、累積能量耗散模型和疲勞安全壽命計算方法的疲勞壽命預(yù)測結(jié)果與低周循環(huán)疲勞試驗結(jié)果。其中，疲勞安全壽命計算依據(jù)文獻(xiàn) [37] 中的Ti6Al4V ELI鈦合金應(yīng)變 - 壽命關(guān)系式Δεt/2=0.0126(2Nf)−0.080+7.565(2Nf)−1.21，通過線性插值法得出不同應(yīng)變幅Δε/2下的循環(huán)壽命Nrc。對比結(jié)果顯示，盡管存在一定的誤差，但基于兩種 FIPs 預(yù)測的疲勞壽命呈現(xiàn)出相似的趨勢，并且在準(zhǔn)確性上高于疲勞安全壽命計算方法。這一發(fā)現(xiàn)驗證了圖 9 和圖 10 中累積塑性滑移和累積能量耗散分布的準(zhǔn)確性，且符合前述提出的兩者在相同位置發(fā)生了局部化增強(qiáng)，并且都隨著循環(huán)次數(shù)的增加呈近似線性上升趨勢的結(jié)論。進(jìn)一步分析發(fā)現(xiàn)，隨著應(yīng)變幅的增加，基于累積塑性滑移的 FIP 相較于基于累積能量耗散的 FIP 表現(xiàn)出更高的預(yù)測準(zhǔn)確性。
圖 14 為塑性應(yīng)變幅值與基于累積塑性滑移和累積能量耗散這 2 種疲勞指示因子預(yù)測的疲勞壽命的對比。在對數(shù)坐標(biāo)系下，這一關(guān)系呈鮮明的線性特征。在描述低周循環(huán)疲勞壽命Ni時，通常采用 Coffin-Manson (C-M) 定律，該定律將塑性應(yīng)變幅值Δεp/2通過式Δεp/2=εf′(Ni)c與Ni聯(lián)系起來，其中εf′為延性系數(shù)，c 為疲勞指數(shù)。利用 FIP 進(jìn)行預(yù)測的結(jié)果揭示了一個重要現(xiàn)象：與宏觀 C-M 模型相比，微觀結(jié)構(gòu)的不同對材料的裂紋萌生壽命產(chǎn)生了顯著影響。這一發(fā)現(xiàn)強(qiáng)調(diào)了微觀結(jié)構(gòu)在疲勞損傷過程中的重要作用。因此在進(jìn)行疲勞壽命預(yù)測時，必須充分考慮材料的微觀特征。
4.3 在三維模型中的適用性
探究疲勞指示因子在三維模型中的適用性時，本研究深入分析了晶粒尺寸與網(wǎng)格數(shù)量對拉伸應(yīng)力 - 應(yīng)變曲線的影響。這對于選定一個既能使模擬曲線與試驗數(shù)據(jù)高度吻合又能確保計算效率的代表性體積元模型至關(guān)重要。借助 Neper 建模軟件，構(gòu)建了體積相同但晶粒數(shù)量各異的模型，具體涵蓋 30、50、100、120 及 150 個晶粒，如圖 15 (a) 所示。對比這些多晶體模型的拉伸應(yīng)力 - 應(yīng)變曲線，注意到含 30 個與 50 個晶粒的模型的屈服強(qiáng)度與含 100 個以上晶粒的模型存在顯著差異。在彈性階段，所有模型的應(yīng)力 - 應(yīng)變曲線近乎重合；然而，步入塑性階段后，100、120 及 150 晶粒模型的曲線展現(xiàn)出高度的相似性，相比之下，30 與 50 晶粒模型的曲線則明顯偏離。基于上述事實，選定 100 晶粒模型，并對其進(jìn)行了細(xì)致的網(wǎng)格敏感性分析。在 Neper 軟件中，分別設(shè)置模型的網(wǎng)格數(shù)量大約為 6 萬、10 萬、18 萬和 46 萬，得出不同網(wǎng)格數(shù)量的多晶體模型拉伸應(yīng)力 - 應(yīng)變曲線，如圖 15 (b) 所示。由圖 15 (b) 可知，網(wǎng)格數(shù)量約為 10 萬的晶粒模型的拉伸應(yīng)力 - 應(yīng)變曲線呈現(xiàn)出良好的一致性。這一結(jié)果意味著，在該網(wǎng)格密度下，模型的計算結(jié)果具有較高的穩(wěn)定性，非常適用于后續(xù)的疲勞壽命預(yù)測。
采用 Neper 建模軟件，依據(jù) EBSD 獲取的晶粒尺寸平均值，構(gòu)建了三維晶粒隨機(jī)生長的代表性體積元模型。此模型遵循正態(tài)分布原則，標(biāo)準(zhǔn)差設(shè)定為 1.23，其中 10% 的晶粒設(shè)定為 β 相，其余 90% 則為 α 相。RVE 的邊長為 1 mm，內(nèi)含 100 個晶粒，網(wǎng)格數(shù)量調(diào)整為約 10 萬個，以確保計算效率與精度。為確保計算的連貫性，對 RVE 施加周期性邊界條件，以模擬材料的真實變形過程。在此基礎(chǔ)上，施加了應(yīng)變比為 0 的循環(huán)位移載荷，并使用第 3 節(jié)中確定的材料參數(shù)，以深入研究三維結(jié)構(gòu)下的疲勞裂紋萌生壽命。
為分析不同應(yīng)變幅下疲勞因子在三維結(jié)構(gòu)中的適用性，本節(jié)沿用第 3 小節(jié)的宏觀模擬方法，重新模擬了 0.4% 應(yīng)變條件下缺口拉伸方向?qū)?yīng)的最大應(yīng)變幅，結(jié)果為 0.45%；試驗測得的裂紋萌生壽命則為 11908 周。圖 16 (a)~ 圖 16 (d) 為應(yīng)變幅值分別為 0.45%、0.63%、0.74% 和 0.87% 時，第 30 周末的累積塑性滑移分布，以表明這些應(yīng)變幅值下的疲勞損傷情況。同時，圖 17 (a)~ 圖 17 (d) 為相應(yīng)應(yīng)變幅值下第 30 周末的累積能量耗散分布。對比分析可知，隨著應(yīng)變幅值的增加，累積塑性滑移和累積能量耗散均呈顯著的上升趨勢。
圖 18 (a) 和圖 18 (b) 為不同應(yīng)變幅下 2 種疲勞指示因子隨循環(huán)次數(shù)變化的規(guī)律。由圖 18 可知，在三維結(jié)構(gòu)中，循環(huán)次數(shù)與疲勞指示因子同樣呈近似線性的關(guān)系。在應(yīng)變幅為 0.63% 的條件下，借助試驗所得的壽命數(shù)據(jù)和模擬至第 30 周末時的累積塑性滑移增量及累積能量耗散增量，利用式 (11) 和式 (13) 分別計算出，累積塑性滑移的臨界值為 143，累積能量耗散的臨界值為 30136 MJ/m³。為了進(jìn)一步驗證這些疲勞指示因子的有效性，圖 19 對比了基于這 2 個因子預(yù)測的裂紋萌生壽命、采用疲勞安全壽命計算方法得出的預(yù)測壽命，以及實際的試驗結(jié)果。由圖 19 可知，基于疲勞指示因子的預(yù)測結(jié)果與試驗結(jié)果高度吻合，且預(yù)測值相對保守；在預(yù)測準(zhǔn)確性方面，基于疲勞指示因子的方法優(yōu)于疲勞安全壽命計算方法。這一發(fā)現(xiàn)驗證了累積塑性滑移和累積能量耗散作為疲勞指示因子在三維結(jié)構(gòu)分析中的適用性。
4.4 討論
盡管本研究在預(yù)測Ti6Al4V ELI鈦合金的疲勞裂紋萌生壽命方面取得了初步成果，但仍存在諸多值得深入探討和改進(jìn)的問題，特別是在模型的維度、疲勞指示因子的線性假設(shè)以及試驗觀察方法上。本研究主要基于 EBSD 構(gòu)建了Ti6Al4V ELI鈦合金的疲勞裂紋萌生壽命預(yù)測模型，但這些模型主要局限于二維和簡化的三維問題。盡管二維模型在計算效率和模型復(fù)雜性方面具有明顯優(yōu)勢，但它無法全面考慮三維微觀結(jié)構(gòu)對疲勞行為的影響。實際上，材料的疲勞行為是一個復(fù)雜的三維問題，涉及晶粒的空間分布、相界面等多種因素。因此，未來的研究應(yīng)致力于開發(fā)真實微觀結(jié)構(gòu)的三維模型，以更精確地模擬實際條件下的疲勞過程。在構(gòu)建這些三維模型時，如何在計算復(fù)雜性和物理真實性之間找到平衡將是一個重大挑戰(zhàn)。對于高度復(fù)雜的微觀結(jié)構(gòu)，可能需要借助先進(jìn)的數(shù)值算法和高性能的服務(wù)器進(jìn)行計算，以更準(zhǔn)確地捕捉各個晶粒和相界面之間的相互作用。此外，本研究假設(shè)累積塑性滑移和累積能量耗散這 2 種疲勞指示因子與循環(huán)次數(shù)存在近似線性關(guān)系。然而，這 2 種疲勞指示因子在高疲勞加載條件下可能表現(xiàn)出非線性特性，因而上述線性假設(shè)雖然簡化了模型，但可能導(dǎo)致對疲勞行為的預(yù)測不夠精確。因此，未來的研究應(yīng)考慮引入非線性疲勞指示因子，以更準(zhǔn)確地描述材料在疲勞載荷下的真實響應(yīng)。盡管本研究得出的結(jié)果與試驗結(jié)果之間表現(xiàn)出良好的吻合性，為模型的有效性提供了有力支持，但將裂紋長度在 0~0.3 mm 假設(shè)為試驗壽命對應(yīng)的裂紋萌生壽命可能存在一定的局限性。為了更真實地反映裂紋萌生壽命，后續(xù)研究需要采用更先進(jìn)的技術(shù)，如高分辨率成像技術(shù)或?qū)崟r監(jiān)測技術(shù)，來觀察裂紋的萌生過程并測量所需數(shù)據(jù)。這些技術(shù)將有助于驗證模型的預(yù)測能力并優(yōu)化模型，從而提高疲勞壽命預(yù)測的準(zhǔn)確性和可靠性。
5、結(jié)論
本研究致力于開發(fā)一個針對雙相微觀結(jié)構(gòu)Ti6Al4V ELI鈦合金的晶體塑性模型，以深入探究該材料的低周疲勞行為，并精確預(yù)測其裂紋萌生壽命。得出如下主要結(jié)論：
為雙相Ti6Al4V ELI鈦合金構(gòu)建了一個基于真實微觀結(jié)構(gòu)的循環(huán) CP 本構(gòu)模型。此模型不僅考慮了背應(yīng)力、分解切應(yīng)力和滑移阻力的動態(tài)演變，還結(jié)合試驗數(shù)據(jù)確定了相關(guān)的 CP 參數(shù)。此外，采用 2 種疲勞指示因子 (累積塑性滑移和累積能量耗散)，以準(zhǔn)確表征疲勞損傷程度。
研究顯示，疲勞損傷更易在 α 相與 β 相的界面處萌生，這一現(xiàn)象與相鄰晶粒間滑移阻力的差異密切相關(guān)。局部晶格旋轉(zhuǎn)不僅加速了微觀紋理的形成，還導(dǎo)致了累積塑性滑移和累積能量耗散的局部集中，從而加速了疲勞裂紋的產(chǎn)生。值得注意的是，每個循環(huán)中晶格旋轉(zhuǎn)的不對稱性引發(fā)了流動行為中拉伸與壓縮的不對稱性。
隨著應(yīng)變幅的增大，累積能量耗散和累積塑性滑移均呈顯著上升趨勢，而疲勞裂紋萌生壽命則相應(yīng)縮短。通過與試驗結(jié)果的對比發(fā)現(xiàn)，基于累積塑性滑移的疲勞指示因子相較于基于累積能量耗散的疲勞指示因子具有更高的預(yù)測準(zhǔn)確性。與經(jīng)典的 C-M 模型相比，基于疲勞指示因子的預(yù)測結(jié)果更清晰地展示了Ti6Al4V ELI鈦合金因微觀結(jié)構(gòu)差異而導(dǎo)致的裂紋萌生壽命的顯著差異。此外，這 2 種疲勞指示因子在三維結(jié)構(gòu)中的適用性也得到了驗證，其預(yù)測結(jié)果不僅準(zhǔn)確，而且偏于保守，為工程應(yīng)用提供了可靠的安全裕量。
參考文獻(xiàn) (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] 高焱，金海波，李遠(yuǎn)杰，等。基于 FKM 的增材制造鈦合金結(jié)構(gòu)疲勞壽命分析方法研究 [J]. 機(jī)械強(qiáng)度，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 鈦合金疲勞裂紋擴(kuò)展研究 [J]. 船舶力學(xué)，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] 陶蒙，于培師，趙軍華。基于晶體塑性與內(nèi)聚力模型的 Cu-Ni-Si 合金晶界起裂研究 [J]. 機(jī)械強(qiáng)度，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] 張海英，李剛，臧偉鋒，等。激光選區(qū)熔化 Ti-6Al-4V 鈦合金疲勞裂紋擴(kuò)展性能研究 [J]. 機(jī)械強(qiáng)度，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] 全國鋼標(biāo)準(zhǔn)化技術(shù)委員會。金屬材料 拉伸試驗 第 1 部分：室溫試驗方法：GB/T 228.1-2021 [S]. 北京：中國標(biāo)準(zhǔn)出版社，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]. 鈦工業(yè)進(jìn)展，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)
（注，原文標(biāo)題：基于微結(jié)構(gòu)特征的Ti6Al4V_ELI鈦合金疲勞裂紋萌生壽命晶體塑性預(yù)測方法）
tag標(biāo)簽:Ti6Al4V,ELI鈦合金,深海服役,疲勞安全性評估