邵 寧, 金 丹
(沈陽化工大學 機械工程學院, 遼寧 沈陽 110142)
022Cr17Ni12Mo2不銹鋼與碳素鋼相比因具備高抗拉強度、低屈服點、極好的塑性和韌性,而且在焊接性能和冷彎成型工藝性能方面表現也尤為突出,所以,在各種儲槽、塔器、反應釜等壓力容器制造領域被廣泛應用.
當前已有大量學者在有限元模擬試驗方面進行了相關研究.金丹[1]等通過ANSYS軟件,采用Von Mises屈服準則、多線性隨動硬化律對Mod.9Cr-1Mo鐵素體鋼V形缺口件進行有限元模擬分析,其模擬試驗結果顯示,缺口根部在各個路徑下的模擬分析中均表現出明顯的應力集中現象;蔣春松[2]等通過Abaqus軟件針對022Cr17Ni12Mo2不銹鋼的彎曲過程進行了微觀有限元模擬,其模擬試驗結果顯示,數值模型表現為接觸面的壓應力分布為中間小,向邊緣地帶程遞增趨勢,且最大值區域為相對靠近邊緣地帶;張星[3]等通過ANSYS軟件建立022Cr17Ni12Mo2不銹鋼納米壓痕試驗的有限元模型,并分析了TiN薄膜與基體的結合性能,其模擬試驗結果顯示,結合性能明顯高于純鈦薄膜,并通過這一模擬試驗驗證了此有限元模擬方法的可行性;周強等[4]通過有限元模擬研究了不同缺口形式對材料的低周疲勞壽命的影響,討論了圓角半徑及缺口深度對應力集中與疲勞壽命的影響程度;李江華[5]通過Abaqus軟件,采用非線性隨動硬化及各向同性強化混合模型,針對316L不銹鋼進行應變控制有限元分析,其模擬試驗結果可以較好地體現試驗中的初始循環硬化現象,并同時驗證了該混合模型的可行性.
本文針對022Cr17Ni12Mo2不銹鋼半圓形缺口件和U形缺口件進行拉-拉疲勞試驗,采用Abaqus軟件進行單軸彈塑性有限元模擬,并將模擬結果與試驗結果進行對比分析.
試驗材料為022Cr17Ni12Mo2不銹鋼,材料的化學成分如表1所示.經過1 080 ℃的保溫再經水冷處理后,加工成標準板材試樣,試樣的尺寸(單位:mm)及形狀如圖1和圖2所示(厚度為1 mm),試件的制備參考文獻[6],理論應力集中系數參照應力集中手冊[7]得到,均為2.4,處理后屈服強度σy=247 MPa,拉伸強度σb=564 MPa,楊氏模量E=200 GPa,泊松比μ=0.3.試驗在EUF-1020電子萬能疲勞試驗機上進行,三角波形控制,標距段3.5 mm,正應力水平下降半壽命載荷的25 %即為失效.試驗條件及試驗結果如表2所示.

表1 316L不銹鋼的化學成分

圖1 半圓形缺口試件

圖2 U形缺口試件

表2 疲勞試驗條件及試驗結果
在定義材料屬性的過程當中,需要將該模型解析為非線性隨動硬化模型和各向同性模型,進而求解其各自參數.在針對非線性隨動硬化模型中,其背應力為[5]:

(1)
εp為塑性應變;γ、c為隨動硬化參數,γ為塑性變形增加時隨動強化模量減小的比率,c為隨動硬化模量初始值.
在針對各向同性模型分析中,屈服面半徑為:
r(p)=Q(1-e-bp)
(2)
p為累積塑性應變;b、Q為等向硬化參數,b為隨塑性應變發展屈服面變化的比率,Q為屈服面的最大變化值.
表3為對半圓形缺口件和U形缺口件在同一應力水平加載條件下試驗數據擬合結果.

表3 材料的彈性及硬化常數
取半圓形缺口件和U形缺口件的標距段進行模擬試驗,將中軸型算法定義為網格屬性,此算法可在兩種缺口形式的試件缺口根部利用最小網格進行過渡.其中,半圓形缺口件的模型單元數為2 796,模型節點為4 987;U形缺口件的模型單元數為3 587,模型節點數為6 831.具體模型實例如圖3、圖4所示.

圖3 半圓形缺口件有限元模型及網格

圖4 U形缺口件有限元模型及網格
將疲勞試驗結果和有限元模擬結果繪制于同一張圖中,其面載荷為442 MPa,循環周次為200圈.對比結果如圖5、圖6所示.
對比分析試驗結果與有限元模擬結果,對于半圓形缺口件,模擬結果與試驗結果的最大應力值誤差為2.1 %,最大應變值誤差為1.8 %;對于U形缺口件,模擬結果與試驗結果的最大應力值誤差為3.6 %,最大應變值誤差為9.7 %.

圖5 半圓形缺口件模擬結果與試驗結果對比

圖6 U形缺口件模擬結果與試驗結果對比
針對022Cr17Ni12Mo2不銹鋼半圓形缺口件和U形缺口試件進行了應力控制低周疲勞試驗和有限元模擬試驗,并對比其結果,得到以下結論:
(1) 對022Cr17Ni12Mo2不銹鋼進行有限元分析時,采用非線性隨動硬化模型以及各向同性模型;將應力轉換為面載荷后運用Abaqus中的Function功能進行載荷加載來實現試驗要求,結果表明材料屬性定義及邊界條件的應用均合理;
(2) 對比模擬路徑結果與試驗路徑可以看出:有限元模擬能很好地實現試驗路徑,模擬得到的應力應變滯回線與試驗結果存在著不同程度的誤差,但就整體而言模擬結果比較滿意.
[1] 金丹,田大將,王巍,等.非比例載荷下缺口件疲勞壽命預測[J].工程力學,2014,31(10):212-215.
[2] 蔣春松,彭金方,沈明學,等.316L不銹鋼構件彎曲微動的有限元模擬及其疲勞壽命預測[J].機械工程材料,2013,37(8):81-84.
[3] 張星,王鶴峰,袁國政,等.基于納米壓痕試驗的316L不銹鋼表面鈦、TiN薄膜結合性能的有限元模擬[J].機械工程材料,2013,37(9):90-95.
[4] 周強.低周疲勞下切口類缺陷的疲勞特性研究[D].蘭州:蘭州理工大學,2013:43-46.
[5] 李江華.DSA效應對316L不銹鋼多軸低周疲勞特性的影響[D].沈陽:沈陽化工大學,2016:33.
[6] 吳啟舟,李江華,緱之飛,等.022Cr17Ni12Mo2不銹鋼缺口試樣的疲勞壽命預測[J].機械工程材料,2015,39(12):55-58
[7] 西田正孝.應力集中[M].李安定,譯.北京:機械工業出版社,1983:57.