【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬

【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖1

Code_Aster是法國電力集團(tuán)(EDF)研發(fā)的一款開源有限元仿真軟件,適用于固體力學(xué)、熱學(xué)和聲學(xué)等物理現(xiàn)象,具體細(xì)分為靜力學(xué)、動力學(xué)、土-結(jié)構(gòu)相互作用、流固耦合、傳熱學(xué)等方面的仿真分析,具備廣闊的應(yīng)用場景。Code_Aster通過核工業(yè)認(rèn)證,滿足法國核安全局和英國核監(jiān)管部門要求,在各工業(yè)領(lǐng)域尤其是能源電力領(lǐng)域有大量的工程和研發(fā)應(yīng)用案例。


01

研究背景

汽輪機是核電站常規(guī)島中最關(guān)鍵的設(shè)備之一,負(fù)責(zé)將二回路水蒸氣的熱能轉(zhuǎn)化為機械能。而套環(huán)是汽輪機的一種鋼制旋轉(zhuǎn)件,安裝在轉(zhuǎn)子軸的兩端,用于固定感應(yīng)器繞組的線圈頭,正常運行工況下,不僅需要承擔(dān)部分汽輪機的自重載荷,還需承受汽輪機葉片轉(zhuǎn)動過程產(chǎn)生的高溫和機械振動(1500轉(zhuǎn)每分鐘)。

【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖2

基于這種情況,細(xì)小的缺陷也很可能導(dǎo)致汽輪機發(fā)生故障,最終造成停堆。由于汽輪機運行工況的周期性特點,汽輪機套環(huán)易發(fā)生疲勞損傷現(xiàn)象,從而導(dǎo)致套環(huán)產(chǎn)生細(xì)微裂紋,或是導(dǎo)致已有的裂紋擴展。

法國電力公司為了評估套環(huán)的安全性和生命周期,使用Code_Aster的疲勞計算模塊進(jìn)行了靜力加載的裂紋擴展分析和循環(huán)加載的疲勞裂紋擴展分析。


02

理論基礎(chǔ)—斷裂和疲勞擴展的原理


   
斷裂的基本類型有三種:  
  • 張開性裂紋(Ⅰ型)
  • 滑開型裂紋(Ⅱ型)
  • 撕開型裂紋(Ⅲ型)
材料和結(jié)構(gòu)的斷裂與否常由應(yīng)力強度因子(Facteurs d‘intensité de contraintes)與材料的標(biāo)準(zhǔn)臨界值K   c對比來判斷,而Ⅰ型裂紋在實際工程應(yīng)用是最容易斷裂的裂紋類型。
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖3  
斷裂的三種類型  
應(yīng)力強度因子計算公式:  
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖4  
能量釋放率公式
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖5  
 對于線彈性材料:  
如果K < KIC,材料不發(fā)生斷裂;  
如果K > KIC,材料發(fā)生斷裂,裂紋擴展。  
疲勞問題的載荷通常以周期性的形式存在,此時的應(yīng)力強度因子也隨加載規(guī)律而呈現(xiàn)周期性變化,材料的裂紋擴展則需要用另一種方法來計算:  
疲勞裂紋擴展速率  【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖6  
即用裂紋長度的增量和交變應(yīng)力的循環(huán)次數(shù)增量的比值來定量計算,表示交變應(yīng)力每循環(huán)一次裂紋長度的平均增量。  
根據(jù)Paris提出的公式   【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖7  
來推得每次循環(huán)的應(yīng)力強度因子的峰值   【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖8  
和材料的疲勞斷裂臨界值   【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖9   對比
判斷材料的裂紋是否會發(fā)生擴展(僅適用于短裂紋)。
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖10  
周期性加載下的應(yīng)力強度因子變化曲線示意圖  
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖11  
Paris公式  

03

模型建立

根據(jù)幾何和邊界條件的對稱性,將汽輪機套環(huán)簡化為1/8模型,通過Salome中的HOMARD模塊將網(wǎng)格作密化處理,最終采用了1501970單元數(shù)量,包含286338節(jié)點數(shù)的TETRA4四面體網(wǎng)格。
   
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖12  
汽輪機套環(huán)的幾何模型  
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖13  
汽輪機套環(huán)的網(wǎng)格  
汽輪機套環(huán)的裂紋萌生和擴展風(fēng)險主要來源有兩點:  
  • 啟堆和停堆期間的汽輪機軸速(離心力)變化,這將導(dǎo)致汽輪機套環(huán)的鼓脹和軸向裂紋發(fā)生;
  • 交變電流產(chǎn)生的熱量引起線圈受熱膨脹,由內(nèi)而外地對套環(huán)進(jìn)行擠壓,這將導(dǎo)致汽輪機套環(huán)容易發(fā)生周向裂紋。
對比兩種導(dǎo)致裂紋產(chǎn)生的原因,熱膨脹引起的周向裂紋概率遠(yuǎn)大于軸速變化引起的軸向裂紋。  
因此在本研究中僅考慮了周向裂紋,選取其常見的表征類型Ⅰ型(張開型)裂紋。  
計算模型主要考慮:  
 
  • 復(fù)雜的熱力學(xué)邊界條件;
  • 接觸和摩擦問題(套環(huán)和線圈等);
  • 不同尺寸的橢圓形裂紋缺陷,使用擴展有限元方法(X-FEM,該方法無需畫出裂紋形狀的網(wǎng)格)。
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖14  
汽輪機套環(huán)的邊界條件  
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖15  
橢圓形裂紋缺陷  

04

計算結(jié)果分析


   

   
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖16    
汽輪機套環(huán)的高斯應(yīng)力云圖      
code_aster的CALC_G模塊提供了計算應(yīng)力強度因子的功能,能夠直接輸出觀測點上的應(yīng)力強度因子。  
在本研究中,我們對套環(huán)最容易產(chǎn)生裂紋的部位(槽)建立了X-FEM單元的橢圓張開型裂紋。下圖為不同尺寸的裂紋和不同方向上的應(yīng)力強度因子曲線:
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖17  
各方向上裂紋尖端的應(yīng)力強度因子  
從圖中可觀察到,左右兩側(cè)的裂紋存在應(yīng)力集中現(xiàn)象,應(yīng)力強度因子最大,即裂紋擴展速率最快;而隨著初始裂紋尺寸的增大,應(yīng)力強度因子的最大值也在增大。  
針對交變電流產(chǎn)生的周期性熱載荷所產(chǎn)生的疲勞問題,裂紋擴展分析主要分為以下步驟:  
  • 應(yīng)力強度因子的計算與周期性的熱量輸出峰值相關(guān):,建立應(yīng)力強度因子與熱量的映射關(guān)系。

【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖18  
應(yīng)力強度因子峰峰值– 熱量變化映射關(guān)系圖  
  • 根據(jù)已有的應(yīng)力強度因子峰峰值計算材料壽命Nf

【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖19  
  • 迭代計算隨循環(huán)次數(shù)而變化裂紋尖端位置,得到下圖(初始橢圓幾何尺寸a=40mm,b=10mm)。
【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖20  
不同尺寸的橢圓裂紋的裂紋尖端位置  
圖中黑色實線為初始的裂紋位置,紅線和黑色圓形線為實驗觀測得到的結(jié)果,綠線和藍(lán)線為計算得到的結(jié)果。  
結(jié)果顯示裂紋最終在套環(huán)邊界處停止,和   真實情況的裂紋相當(dāng)吻合。   在應(yīng)力集中的左右兩側(cè),真實情況的裂紋位置在計算的裂紋外側(cè),表明此處的真實裂紋擴展速率比計算的更大,計算結(jié)果相對保守。  

05

結(jié)論

該案例主要對汽輪機套環(huán)進(jìn)行了斷裂-疲勞裂紋擴展分析,考慮了由交變電流引起的周期性熱載荷,利用測得的交變電流的頻率,可以推斷裂紋箍在循環(huán)載荷下的疲勞壽命,并從而研究最佳的汽輪機運行控制策略,為汽輪機套環(huán)部件的壽命評估提供了重要分析依據(jù)。   

近期Code_Aster培訓(xùn)信息

Code_Aster動力學(xué)分析專題培訓(xùn)


#

培訓(xùn)時間

2021年3月17日~19日

報名截止時間:2021年3月16日

本次培訓(xùn)名額為30人,名額有限,報名從速。

#

培訓(xùn)地點

遠(yuǎn)算培訓(xùn)教室(杭州市西湖區(qū)千島湖智谷大廈)

#

擬定培訓(xùn)內(nèi)容

1、Code_Aster軟件進(jìn)行固體力學(xué)仿真計算的基礎(chǔ)操作和功能介紹

2、Code_Aster軟件針對結(jié)構(gòu)動力學(xué)模態(tài)&頻響問題的基本解決方案和實踐指導(dǎo)

3、Code_Aster軟件針對瞬態(tài)問題的非線性和模態(tài)疊加解決方案對比和實踐指導(dǎo)

4、Code_Aster軟件在工業(yè)領(lǐng)域的應(yīng)用實例

#

報名方式

微信掃描下方二維碼線上報名,我們會在一個工作日內(nèi)聯(lián)系您!

【EDF開源CAE案例】Code_Aster對汽輪機套環(huán)的疲勞損傷與裂紋擴展模擬的圖21


更多資訊可登錄格物CAE官方網(wǎng)站

https://cae.yuansuan.cn/

遠(yuǎn)算在bilibili、技術(shù)鄰、知乎定期發(fā)布課程視頻等內(nèi)容

敬請關(guān)注

技術(shù)鄰.jpg


登錄后免費查看全文
立即登錄
App下載
技術(shù)鄰APP
工程師必備
  • 項目客服
  • 培訓(xùn)客服
  • 平臺客服

TOP

8
6
1