林 磊1,陳乃斌,徐德城1,朱如東
(1.蘇州熱工研究院有限公司,江蘇 蘇州 215004;2.陽江核電有限公司,廣東 陽江 529941)
公稱直徑2 in以下小支管的振動疲勞斷裂是國內(nèi)外核電廠普遍問題之一,其失效占核電廠管道振動疲勞失效總量的90%以上[1]。一些重要小支管的斷裂將導(dǎo)致機(jī)組停機(jī),造成巨大經(jīng)濟(jì)損失。近幾年來,某CPR1000核電機(jī)組頻繁發(fā)生汽輪機(jī)高壓缸進(jìn)汽管道上GPV/GRE系統(tǒng)儀表管線斷裂事件,造成巨大經(jīng)濟(jì)損失。為緩解相關(guān)儀表管線的振動水平,核電廠開展了多項管線改進(jìn)措施[2],一定程度上降低了其振動交變應(yīng)力水平,但仍存在啟機(jī)過程低功率下振動應(yīng)力高于允許限值的問題。
小支管的振動疲勞屬于受應(yīng)力幅控制的高周疲勞,其破壞過程通常包含了裂紋萌生、擴(kuò)展和瞬斷三個階段,其中裂紋的萌生出現(xiàn)在局部最高應(yīng)力處、在最弱的及應(yīng)力最大的晶粒上形成微裂紋,然后發(fā)展為宏觀裂紋,最終導(dǎo)致疲勞斷裂[3]。疲勞曲線(又稱SN曲線)是對構(gòu)件進(jìn)行疲勞設(shè)計或評估的基礎(chǔ),基于疲勞曲線和Miner線性累積損傷準(zhǔn)則,通過計算總的疲勞累積使用因子,可以進(jìn)行有限壽命或無限壽命設(shè)計與評估[4]。
與上述正向疲勞設(shè)計的方法相反,對于小支管的振動疲勞問題,本文借助材料的疲勞曲線并基于Miner線性累積損傷準(zhǔn)則,基于實測振動應(yīng)變信號,來反推對應(yīng)于有限壽命和無限壽命的應(yīng)變范圍限值,用于全新安裝小支管在機(jī)組一次啟機(jī)過程中的實時監(jiān)測預(yù)警,當(dāng)實際應(yīng)變范圍高于限值時,啟機(jī)保駕人員即采取相應(yīng)措施,避免小支管在啟機(jī)過程中疲勞斷裂。
疲勞損傷與交變應(yīng)力水平、結(jié)構(gòu)不連續(xù)處的應(yīng)力集中情況、平均應(yīng)力水平等有關(guān)系,且是一個不斷累積、最終發(fā)生裂紋萌生的過程。當(dāng)交變應(yīng)力水平低于材料的疲勞強(qiáng)度時,理論上管道具有無限壽命。只有當(dāng)交變應(yīng)力水平高于材料的疲勞強(qiáng)度時,才會形成損傷。
結(jié)構(gòu)振動疲勞與典型低周疲勞(如壓力或溫度變化引起的疲勞)的區(qū)別在于,振動疲勞循環(huán)的快慢與結(jié)構(gòu)動力學(xué)特性相關(guān),在結(jié)構(gòu)幾何特征不變、但激振力變化后,結(jié)構(gòu)的振動響應(yīng)幅值和頻率也將發(fā)生變化,則相同時間內(nèi)形成的累積損傷大大不同;而當(dāng)激振力保持不變,但振動持續(xù)時間變化后,相同振動應(yīng)力在不同時間內(nèi)形成的累積損傷也發(fā)生變化。本文的研究中,使用某機(jī)組上一次啟機(jī)過程儀表管的實測振動應(yīng)變數(shù)據(jù),來計算下一次啟機(jī)的應(yīng)變范圍限值,需作出如下假設(shè)。
(1)下次啟機(jī)前汽輪機(jī)高壓缸等相連的主要設(shè)備未經(jīng)過全檢等大范圍檢修工作,設(shè)備的振動特性與上一循環(huán)基本保持一致;
(2)下一次啟動的功率曲線與上一次啟機(jī)功率曲線保持一致。
小支管為細(xì)長型結(jié)構(gòu),呈現(xiàn)出梁的橫向振動特征,其橫截面主要受彎矩作用,其外表面以軸向變形為主,軸向應(yīng)變幅值最大。但由于小支管為三維走向,除受振動彎矩作用外,還會受到相鄰管段運(yùn)動帶來的剪切力和扭矩,因此其變形比較復(fù)雜,并非理想的純彎曲變形。針對不同受力情況下管道的疲勞損傷評估,在ASME BPVC標(biāo)準(zhǔn)中提出了兩種方法:1)主應(yīng)力方向不變的情況,可對單軸應(yīng)力通過雨流統(tǒng)計進(jìn)行循環(huán)計數(shù)和累積使用因子計算;2)主應(yīng)力方向變化的情況,可采用如最大最小應(yīng)力法進(jìn)行循環(huán)計數(shù)和累積使用因子計算,該方法不考慮加載次序,具有極大的保守性[6]。而ASME OM Part3標(biāo)準(zhǔn)中對振動應(yīng)力的表達(dá)式,則是僅考慮橫向彎曲變形。為降低實際監(jiān)測評估時的復(fù)雜性,本文提出第(3)個假定,即假定小支管以橫向純彎曲振動為主,僅考慮其軸向應(yīng)力應(yīng)變。
對于上述第(1)條假設(shè),根據(jù)以往多臺機(jī)組相似位置小支管的測試經(jīng)驗,在啟機(jī)過程中其振動時域和頻域特征基本相似。對于上述第(2)條假設(shè),通過對比歷次啟機(jī)功率曲線,可以確定每次啟機(jī)曲線基本一致,但部分功率平臺持續(xù)時間會有差異。綜合考慮后,確定機(jī)組并網(wǎng)至50%功率平臺結(jié)束的總時間T0=53 h。為防止由于實測應(yīng)變范圍遠(yuǎn)大于53 h應(yīng)變范圍限值而導(dǎo)致的短時間疲勞斷裂風(fēng)險,同時保證啟機(jī)保駕人員有充足的時間進(jìn)行機(jī)組后撤處理,增加T1=2 h和T2=6 h作為補(bǔ)充條件。
因應(yīng)變幅等于應(yīng)變范圍的二分之一,故純彎曲振動管道無限壽命對應(yīng)的最大軸向應(yīng)變范圍限值εr-I為960 με。由于儀表管斷裂主要發(fā)生在對接焊縫處,在進(jìn)行應(yīng)變監(jiān)測數(shù)據(jù)的分析時,還應(yīng)將實測數(shù)據(jù)乘以對接焊縫的疲勞強(qiáng)度減弱系數(shù)(或應(yīng)力集中系數(shù))2.0,才是最終評估數(shù)據(jù)。因此,用于監(jiān)測限值的應(yīng)變范圍應(yīng)取為εr-I=960/2=480 με。當(dāng)實測應(yīng)變范圍低于εr-I時,管道不存在振動疲勞失效的風(fēng)險。
由于儀表管的振動是穩(wěn)態(tài)隨機(jī)振動,其振動頻率范圍較寬,不同量級的振動造成的疲勞損傷程度也不同,且大部分振動循環(huán)的交變應(yīng)力幅值是低于材料疲勞強(qiáng)度的,不會造成疲勞損傷。也即在給定時間內(nèi),管道的最大交變應(yīng)力可以超過材料疲勞強(qiáng)度,但并不會引起疲勞失效。因此,采用無限壽命的應(yīng)變范圍作為啟機(jī)過程的后撤預(yù)警值會過于保守,從而導(dǎo)致不必要的機(jī)組狀態(tài)后撤。而基于疲勞累積使用因子來倒推最大允許應(yīng)變范圍,則可給出能夠保證一定安全運(yùn)行時間的有限壽命應(yīng)變范圍限值。
1.3.1 振動應(yīng)變基本數(shù)據(jù)
根據(jù)已有監(jiān)測數(shù)據(jù)的分析,啟機(jī)過程中儀表管線振動過大的原因為所連接主管道存在高頻殼壁振動,振動頻率以400 Hz以上的高頻成分為主。圖1為所關(guān)注儀表管線的外形圖,實測振動應(yīng)力較大的位置分別為圖中所示根部焊縫和下游焊縫。啟機(jī)過程管線各位置的振動水平及頻率特征并不相同,其根部更接近高壓缸進(jìn)汽管,相應(yīng)的振動頻率與主管基本一致;下游除存在高頻成分外,還伴有低頻振動。本文分別以兩個位置的測試數(shù)據(jù)為基礎(chǔ)數(shù)據(jù),計算應(yīng)變范圍限值。
圖1 儀表管線布置示意圖Fig.1 Layout of typical instrumental pipe
啟機(jī)過程中,在并網(wǎng)至50%功率期間儀表管線的振動應(yīng)力要普遍高于后續(xù)工況,在機(jī)組功率達(dá)到50%左右后,管線振動水平快速下降,管線離開危險工況。而在并網(wǎng)至50%功率期間,儀表管線的振動應(yīng)變趨勢為先增大后緩慢降低,其中電功率120~138 MW期間的振動應(yīng)力水平相對最高,如圖2所示。因此,為了提高保守性,并降低限值計算的復(fù)雜度,取電功率120~138 MW的軸向振動應(yīng)變?yōu)榛緮?shù)據(jù),用于進(jìn)行監(jiān)測限值的計算,取樣時間長度t0為2 100 s,在并網(wǎng)至50%功率的全過程中,均以該取樣數(shù)據(jù)進(jìn)行計算。
圖2 軸向振動應(yīng)變監(jiān)測數(shù)據(jù)Fig.2 Axialvibration strain measurement curve
1.3.2 計算方法
根據(jù)Miner線性累積損傷準(zhǔn)則,對于總運(yùn)行時間為T的n個應(yīng)力循環(huán),假設(shè)第i個應(yīng)力循環(huán)的交變應(yīng)力幅為Si,對應(yīng)的循環(huán)次數(shù)Ni,SN曲線上Si所對應(yīng)的允許循環(huán)次數(shù)為Nai,那么第i個應(yīng)力循環(huán)所產(chǎn)生的疲勞使用因子ui=Ni/Nai。計算所有n個應(yīng)力循環(huán)的累積疲勞使用因子CUF,即:
當(dāng)U<1.0,則表明在壽期內(nèi)不會發(fā)生疲勞裂紋萌生,預(yù)測總的疲勞壽命(指裂紋萌生前)為LF=T/U;當(dāng)U≥1.0,表明會發(fā)生疲勞裂紋的萌生。因此,為保證在有限運(yùn)行時間內(nèi)管線最大應(yīng)力處不發(fā)生裂紋萌生,則應(yīng)使總的疲勞使用因子U<1.0,從保守的角度出發(fā),取U=0.8作為其限值。
應(yīng)變范圍的計算流程如圖3所示,分別為:
圖3 有限壽命對應(yīng)應(yīng)變范圍限值計算流程Fig.3 Procedure of limited-life strain range calculation
(1)以實測管線取樣時間t0內(nèi)的軸向應(yīng)變數(shù)據(jù)為基準(zhǔn),按照胡克定律將應(yīng)變轉(zhuǎn)化為應(yīng)力,并考慮局部疲勞強(qiáng)度減弱系數(shù)的影響;
(2)采用雨流統(tǒng)計法對振動應(yīng)力數(shù)據(jù)進(jìn)行統(tǒng)計,得到交變應(yīng)力幅向量[S]和循環(huán)次數(shù)向量[N];
(3)對(2)中得到的交變應(yīng)力幅和循環(huán)次數(shù)向量,采用Miner線性累積損傷準(zhǔn)則,計算疲勞累積使用因子u0;
(4)擴(kuò)展得到T時間內(nèi)的總累積使用因子U=T×u0/t0;
(5)若U<0.8,則對基本應(yīng)變數(shù)據(jù)進(jìn)行等比例放大,直到U=0.8,獲得此時對應(yīng)的最大應(yīng)變范圍εL0-T;
(6)考慮對接焊縫的應(yīng)力集中系數(shù)2.0,取有限壽命對應(yīng)的應(yīng)變范圍限值為εL-T=εL0-T/2,即可用作實測應(yīng)變數(shù)據(jù)的監(jiān)測限值。
(7)針對T0=53 h、T1=2 h和T2=6 h,分別計算三級有限壽命應(yīng)變范圍限值:
· 保證53 h安全運(yùn)行的應(yīng)變范圍限值εL-53 h;
· 保證6 h安全運(yùn)行的應(yīng)變范圍限值εL-6 h;
· 保證2 h安全運(yùn)行的應(yīng)變范圍限值εL-2 h。
1.3.3 計算結(jié)果
計算采用MATLAB程序?qū)崿F(xiàn),結(jié)果如表1所示。從表中可以看出,根據(jù)根部焊縫實際測試數(shù)據(jù)計算得到的應(yīng)變范圍限值要低于下游焊縫實測數(shù)據(jù)得到的限值,這與上文所分析的兩處頻率成分差異相一致。下游焊縫的振動頻率低,在相同時間內(nèi)允許的振動幅值就較大,而上游焊縫振動頻率高,允許的振動幅值自然要小。從保守的角度出發(fā),機(jī)組啟機(jī)時取根部焊縫的應(yīng)變范圍限值作為監(jiān)測預(yù)警值。對比上文計算的無限壽命應(yīng)變范圍限值,能夠保證機(jī)組一次安全啟機(jī)的允許應(yīng)變范圍接近無限壽命應(yīng)變范圍的2倍,因此采用有限壽命的應(yīng)變范圍限值能夠降低機(jī)組后撤或停機(jī)概率,在保證儀表管運(yùn)行安全性的同時提高經(jīng)濟(jì)性。
表1 有限壽命對應(yīng)的應(yīng)變范圍限值計算結(jié)果
對于全新安裝小支管的一次啟機(jī)過程,基于實測振動應(yīng)變數(shù)據(jù),采用雨流統(tǒng)計法和Miner線性累積損傷準(zhǔn)則,計算了保證無限壽命、53 h、6 h和2 h安全運(yùn)行的應(yīng)變范圍限值,分別為:480 με、856 με、1028 με、1160 με。
由于管道振動特征各不相同,本文所計算應(yīng)變范圍限值中,無限壽命應(yīng)變范圍限值可普遍適用,而有限壽命應(yīng)變范圍限值則僅能用于本文所述特定電廠、特定管道。但本文所提出的應(yīng)變范圍限值計算方法,對于穩(wěn)態(tài)隨機(jī)或穩(wěn)態(tài)周期振動的管道具有普遍適用性。通過有限壽命應(yīng)變范圍限值的計算,可用于管道實際振動狀態(tài)的監(jiān)測和預(yù)警、減少不必要的機(jī)組狀態(tài)后撤或停機(jī),并可為管道的更換或維修計劃提供參考,具有較大的實際工程應(yīng)用價值。