發(fā)布時(shí)間：2022-01-27 02:22

　　為研究某航空發(fā)動(dòng)機(jī)渦輪盤(pán)典型材料粉末高溫合金FGH96高溫下的變形特性,基于Bodner-Partom（B-P）統(tǒng)一粘塑性本構(gòu)理論對(duì)其高溫下的力學(xué)行為進(jìn)行本構(gòu)建模。開(kāi)展了550℃下的單軸拉伸及低周疲勞試驗(yàn),利用Levenberg-Marquardt算法對(duì)模型參數(shù)進(jìn)行了優(yōu)化識(shí)別。采用隱式積分算法將本構(gòu)方程離散為差分方程組,推導(dǎo)了一致切線剛度矩陣,為提高積分過(guò)程的精度和可靠性,引入以非彈性應(yīng)變?cè)隽繛槎攘康姆e分步長(zhǎng)控制策略。通過(guò)用戶子程序接口UMAT將B-P模型引入到ABAQUS有限元軟件中進(jìn)行數(shù)值模擬。結(jié)果表明:高溫下FGH96材料表現(xiàn)出一定的率相關(guān)性及循環(huán)軟化特性,模擬曲線與試驗(yàn)結(jié)果平均相對(duì)誤差在10%以內(nèi),具有較好的一致性,說(shuō)明B-P模型能夠較好地模擬FGH96合金高溫下的變形特性,驗(yàn)證了本實(shí)驗(yàn)?zāi)Ｐ团cUMAT子程序的準(zhǔn)確性。 

【文章來(lái)源】：材料導(dǎo)報(bào). 2020,34(16)北大核心EICSCD

【文章頁(yè)數(shù)】：6 頁(yè)

【部分圖文】：

有限元仿真應(yīng)力云圖:(a)單軸拉伸，(b)低周疲勞(電子版為彩圖)

通過(guò)用戶子程序UMAT接口將B-P模型植入到ABAQUS軟件中，UMAT流程圖如圖1所示。利用表1中的材料參數(shù)，對(duì)FGH96合金在550℃下的單軸拉伸及低周疲勞行為進(jìn)行數(shù)值模擬。幾何模型的約束條件為一端固定，一端受載，采用C3D8R單元對(duì)模型進(jìn)行網(wǎng)格劃分，節(jié)點(diǎn)總數(shù)為384 650，單元總數(shù)為364 536，單軸拉伸模型施加5%的位移載荷，低周疲勞模型施加應(yīng)力控制的三角波。圖2給出了FGH96合金單軸拉伸(1×10-3s-1)和低周疲勞(850 MPa)應(yīng)力云圖。圖2 有限元仿真應(yīng)力云圖:(a)單軸拉伸，(b)低周疲勞(電子版為彩圖)

圖3給出了FGH96合金在550℃時(shí)不同應(yīng)變速率(1×10-3s-1、5×10-4s-1)下的單軸拉伸及有限元模擬曲線。從圖3中可以看出，材料在彈性段受應(yīng)變速率影響不明顯，兩組不同應(yīng)變速率下的彈性變形曲線基本重合，而當(dāng)材料進(jìn)入屈服階段后不同的應(yīng)變速率對(duì)材料的形變影響比較明顯，當(dāng)應(yīng)變速率由5×10-4s-1增加到1×10-3s-1時(shí)，材料的屈服強(qiáng)度增加了約54 MPa，由此可見(jiàn)，高溫下FGH96合金表現(xiàn)出一定的拉伸率相關(guān)性。采用B-P模型對(duì)FGH96合金高溫下不同應(yīng)變速率拉伸曲線進(jìn)行了數(shù)值模擬，平均相對(duì)誤差分別為2.07%(1×10-3s-1)、3.05%(5×10-4s-1)，說(shuō)明B-P模型能夠準(zhǔn)確地表征該材料在高溫下的拉伸變形特性，同時(shí)也能對(duì)材料的率相關(guān)行為進(jìn)行較好模擬。圖4給出了FGH96合金550℃時(shí)不同應(yīng)力水平(850MPa、950 MPa)下低周疲勞循環(huán)響應(yīng)規(guī)律及有限元模擬曲線。從圖4中可以看出，隨著應(yīng)力水平的增加，材料的循環(huán)應(yīng)變范圍也增加，穩(wěn)態(tài)滯回曲線面積(即循環(huán)應(yīng)變遲滯能)也相應(yīng)增大，疲勞壽命減小;同一應(yīng)力水平下應(yīng)變范圍隨著循環(huán)數(shù)的增加而增加，表現(xiàn)出明顯的循環(huán)軟化現(xiàn)象。采用B-P模型對(duì)FGH96合金高溫下的循環(huán)應(yīng)變范圍響應(yīng)及穩(wěn)態(tài)滯回曲線進(jìn)行了數(shù)值模擬，循環(huán)應(yīng)變范圍響應(yīng)模擬平均相對(duì)誤差分別為1.46%(850 MPa)、1.52%(950 MPa)，說(shuō)明B-P模型能夠?qū)Σ牧系难h(huán)軟化現(xiàn)象進(jìn)行準(zhǔn)確模擬。穩(wěn)態(tài)滯回曲線模擬平均相對(duì)誤差分別為8.43%(850 MPa)、9.15%(950MPa)，與試驗(yàn)結(jié)果相比，模擬結(jié)果在加載和卸載終段比較平坦，出現(xiàn)“過(guò)方”現(xiàn)象，這是由于B-P本構(gòu)模型運(yùn)動(dòng)硬化變量的演化方程中缺少動(dòng)力恢復(fù)項(xiàng)，使得模型的彈性階段與穩(wěn)態(tài)流動(dòng)塑性階段迅速銜接[33]。此外，本構(gòu)模型的變形方程中缺少粘彈性變形也會(huì)造成遲滯環(huán)加/卸載過(guò)渡段出現(xiàn)“過(guò)方”現(xiàn)象[36]。因此，需要對(duì)B-P模型進(jìn)行改進(jìn)，以提高對(duì)材料變形過(guò)渡段的模擬精度。

【參考文獻(xiàn)】：
期刊論文
[1]焊接參數(shù)對(duì)FGH96慣性摩擦焊接頭組織和高溫拉伸性能的影響[J]. 楊俊,李京龍,董登科,廖江海.  航空材料學(xué)報(bào). 2019(02)
[2]FGH96激光熔化沉積成形過(guò)程對(duì)組織與性能的影響[J]. 劉帥,張雪,于海鑫,劉常升.  航空制造技術(shù). 2019(Z1)
[3]表面完整性對(duì)FGH96粉末合金高溫低循環(huán)疲勞性能的影響研究[J]. 王欣,陳星,王曉峰,宋穎剛,湯智慧,鄒金文.  稀有金屬材料與工程. 2019(01)
[4]GH4169合金的循環(huán)本構(gòu)模型研究[J]. 盧孔漢,張宏建,賈鵬超.  推進(jìn)技術(shù). 2019(02)
[5]梯度組織FGH96合金殘留枝晶區(qū)的組織特征研究[J]. 田高峰,陳陽(yáng),汪煜.  粉末冶金技術(shù). 2018(06)
[6]熱變形對(duì)FGH96高溫合金原始顆粒邊界的影響[J]. 傅豪,王夢(mèng)雅,紀(jì)箴,田高峰,秦恬,賈成廠.  粉末冶金技術(shù). 2018(03)
[7]航空發(fā)動(dòng)機(jī)高溫材料的研究現(xiàn)狀及展望[J]. 付青峰,楊細(xì)蓮,劉克明.  熱處理技術(shù)與裝備. 2018(03)
[8]熱變形工藝對(duì)FGH96高溫合金PPB析出相及組織的影響[J]. 劉昌奎,魏振偉,張佳慶,鄭真.  航空材料學(xué)報(bào). 2018(03)
[9]粉末高溫合金的研究現(xiàn)狀[J]. 郭茂文,劉春榮,鄭雪萍,劉超,劉勝林.  熱加工工藝. 2017(20)
[10]粉末冶金鎳基高溫合金FGH96高溫疲勞壽命分散性特征[J]. 苗國(guó)磊,楊曉光,石多奇.  航空動(dòng)力學(xué)報(bào). 2017(02)

博士論文
[1]粉末高溫合金材料的力學(xué)特性及其在渦輪盤(pán)上的應(yīng)用研究[D]. 楊治國(guó).南京航空航天大學(xué) 2007
[2]粘彈塑性統(tǒng)一本構(gòu)理論[D]. 馮明琿.大連理工大學(xué) 2000



本文編號(hào)：3611572