基于虛擬裂紋閉合法,利用控制單變量變化的方法,對含對稱雙邊緣裂紋的薄壁圓柱殼受對稱集中力的情況進行了數(shù)值建模,得到了裂紋尖端的能量釋放率,分析了Zr-4合金包殼管的預制裂紋長度和外部載荷對應力強度因子的影響作用。使用柔度法計算并獲得了應力強度因子的近似計算公式,并與數(shù)值計算結果進行了對比驗證,驗證結果吻合較好。 

【文章來源】：核動力工程. 2019,40(05)北大核心EICSCD

【文章頁數(shù)】：5 頁

【部分圖文】：

Zr-4合金應力應變曲線

頹庵卵映倏?眩―HC）等力學分析有非常重要的作用。本文基于虛擬裂紋閉合法（VCCT），利用控制單變量變化的方法，得到了裂紋尖端的能量釋放率，分析了Zr-4合金包殼管的幾何情況（預制裂紋長度）和外部載荷對應力強度因子的影響作用。使用柔度法計算并獲得了應力強度因子的近似計算公式，并與數(shù)值計算結果進行對比驗證。1材料參數(shù)及幾何模型包殼管的材料為Zr-4合金，認為是彈塑性材料。Zr-4合金的彈性模量為95.2GPa，屈服強度為306.5MPa，泊松比為0.33[4]，其應力應變曲線如圖1所示。包殼管的幾何模型如圖2所示，為薄壁圓柱殼形式，圖中L為圓管長度，R為圓管外半徑，t為圓管厚度，試樣受端部集中力加載。設定包殼管試樣的長度L為55mm，按照AFA3G燃料組件的燃料包殼管設定外半徑R為4.75mm，管厚度t為0.57mm。由于模型是對稱的，為提高計算效率，取1/2模型進行計算。模型信息如圖3所示。通過計算可獲得不同載荷下的應力云圖，在F為120N、a為13mm時的vonMises云圖如圖4所示，可見裂尖應力集中明顯。2裂尖能量釋放率的計算VCCT由勢能方法發(fā)展而來，最初由Rybicki和Kanninen[5]提出。該方法假定裂紋擴展Δc引起圖1Zr-4合金應力應變曲線Fig.1Stress-StrainRelationofZr-4Alloy圖2邊緣雙裂紋包殼管幾何模型Fig.2GeometricalModelofCladdingTubewithDoubleEdgeCracks圖3端部受對稱雙集中力的1/2模型mmFig.31/2ModelofTubeSubjectedSymmetricalDoubleConcentratedForcesmm圖4vonMises應力云圖（a=13mm，F(xiàn)=120N）Fig.4vonMisesStressNephogram(a=13mm，F(xiàn)=120N)的系統(tǒng)勢能改變等于閉合裂紋至原來長度所作的功，其表達式見式（1）：Δ01()(Δ)d2cWurrc

?萍撲愎?劍?⒂朧?導撲憬峁??卸員妊櫓ぁ?1材料參數(shù)及幾何模型包殼管的材料為Zr-4合金，認為是彈塑性材料。Zr-4合金的彈性模量為95.2GPa，屈服強度為306.5MPa，泊松比為0.33[4]，其應力應變曲線如圖1所示。包殼管的幾何模型如圖2所示，為薄壁圓柱殼形式，圖中L為圓管長度，R為圓管外半徑，t為圓管厚度，試樣受端部集中力加載。設定包殼管試樣的長度L為55mm，按照AFA3G燃料組件的燃料包殼管設定外半徑R為4.75mm，管厚度t為0.57mm。由于模型是對稱的，為提高計算效率，取1/2模型進行計算。模型信息如圖3所示。通過計算可獲得不同載荷下的應力云圖，在F為120N、a為13mm時的vonMises云圖如圖4所示，可見裂尖應力集中明顯。2裂尖能量釋放率的計算VCCT由勢能方法發(fā)展而來，最初由Rybicki和Kanninen[5]提出。該方法假定裂紋擴展Δc引起圖1Zr-4合金應力應變曲線Fig.1Stress-StrainRelationofZr-4Alloy圖2邊緣雙裂紋包殼管幾何模型Fig.2GeometricalModelofCladdingTubewithDoubleEdgeCracks圖3端部受對稱雙集中力的1/2模型mmFig.31/2ModelofTubeSubjectedSymmetricalDoubleConcentratedForcesmm圖4vonMises應力云圖（a=13mm，F(xiàn)=120N）Fig.4vonMisesStressNephogram(a=13mm，F(xiàn)=120N)的系統(tǒng)勢能改變等于閉合裂紋至原來長度所作的功，其表達式見式（1）：Δ01()(Δ)d2cWurrcr（1）式中，W為閉合裂紋至原來長度所作功；Δc為裂紋擴展貫穿每個網(wǎng)格的長度；u為裂尖張開相對位移（圖5）；σ為應力；r為到裂紋尖端的距離。圖5中Elm為單元，i、j、k為裂紋擴展方向上各單元間節(jié)點，X1與X2分別為沿開裂方向和垂直開裂方向的直角坐標軸。因此，G的表達式為


本文編號：3441348