張貴濱,張振禹
?
屏蔽主泵貧鈾飛輪完整性分析
張貴濱1,2,張振禹1,2
(1.黑龍江省核主泵工程技術(shù)研究中心,哈爾濱 150066;2. 哈爾濱電氣動力裝備有限公司,哈爾濱 150066)
為避免反應(yīng)堆緊急停堆或斷電事故中堆芯損壞,需要在屏蔽主泵電機(jī)上安裝儲能飛輪,為使其在有限的空間內(nèi)產(chǎn)生較大的轉(zhuǎn)動慣量,屏蔽主泵飛輪材料采用高密度貧鈾合金金屬材料。本文主要對貧鈾飛輪的完整性進(jìn)行了研究,從應(yīng)力強(qiáng)度和斷裂力學(xué)兩方面對貧鈾飛輪進(jìn)行了分析,結(jié)果表明:在設(shè)計和額定轉(zhuǎn)速下飛輪的應(yīng)力強(qiáng)度和斷裂強(qiáng)度因子滿足設(shè)計規(guī)范的要求,并通過解析分析結(jié)果與有限元分析結(jié)果的對比表明應(yīng)力分析結(jié)果合理可信。
主泵屏蔽電機(jī);貧鈾飛輪;應(yīng)力強(qiáng)度;斷裂力學(xué)
反應(yīng)堆冷卻劑泵是壓水堆核電站一回路系統(tǒng)的關(guān)鍵設(shè)備之一,其主要功能是驅(qū)動反應(yīng)堆冷卻劑在一回路系統(tǒng)中循環(huán),將堆芯產(chǎn)生的熱量傳遞給蒸發(fā)器二次側(cè),是有效調(diào)整反應(yīng)堆溫度、保證反應(yīng)堆安全的重要設(shè)備[1,2]。為了確保反應(yīng)堆緊急停堆或全廠斷電事故等工況下堆芯的安全,要求主泵具有足夠的惰轉(zhuǎn)時間,繼續(xù)驅(qū)動反應(yīng)堆冷卻劑循環(huán),冷卻堆芯,需要在屏蔽主泵電機(jī)軸上安裝大慣量的飛輪,以保證反應(yīng)堆停堆后主泵具有足夠的惰轉(zhuǎn)能力,保證反應(yīng)堆的安全[3]。
飛輪是主泵重要的功能部件之一,屬于核安全相關(guān)部件,為了確保飛輪能夠在高轉(zhuǎn)速下長期運(yùn)行,需要對飛輪結(jié)構(gòu)的完整性進(jìn)行分析[4]。本文分別采用解析計算方法與有限元計算方法對飛輪的應(yīng)力強(qiáng)度進(jìn)行對比分析,證明了方法的正確性和可行性,并依據(jù)法規(guī)對分析結(jié)果進(jìn)行判定[5~8],以確保飛輪結(jié)構(gòu)的完整性。
1.1 飛輪結(jié)構(gòu)
屏蔽主泵飛輪位于承壓邊界內(nèi),整體浸漬在高溫高壓的水中。主泵運(yùn)行時大直徑的飛輪會產(chǎn)生大量的摩擦損耗使機(jī)組效率降低并使主泵機(jī)組的整體尺寸變大,而小直徑的飛輪又無法滿足機(jī)組所需的轉(zhuǎn)動慣量的要求[9]。為降低主泵機(jī)組的整體尺寸,以及提高機(jī)組使用效率,屏蔽主泵飛輪材料采用高密度貧鈾合金金屬材料能使其在有限的空間內(nèi)的轉(zhuǎn)動慣量達(dá)到最大,該種材料密度約為18.5g/cm3。屏蔽主泵貧鈾飛輪結(jié)構(gòu)示意圖如圖1所示,由貧鈾盤、輪轂、鍵、側(cè)蓋板和外殼組成。貧鈾盤采用圓環(huán)結(jié)構(gòu),為方便進(jìn)行無損檢測,貧鈾飛輪由9個等厚的貧鈾盤疊壓而成,采用平鍵進(jìn)行定位。貧鈾盤與輪轂之間采用過盈連接:在飛輪兩個端面及外表面分別設(shè)置不銹鋼蓋板及屏蔽外殼,與飛輪輪轂構(gòu)成封閉腔體,避免貧鈾合金受到反應(yīng)堆冷卻劑的腐蝕。
圖1 屏蔽主泵貧鈾飛輪結(jié)構(gòu)示意圖
1.2 飛輪材料特性
表1飛輪材料特性
輪轂貧鈾盤 Sy/MPa245400 Su/MPa465670 ρ/(g/cm3)7.7518.5 E/MPa2.05×105 1.6×105 v0.30.3
注:材料特性中屈服強(qiáng)度與抗拉強(qiáng)度均為設(shè)計溫度下的材料特性參數(shù)。
1.3 飛輪設(shè)計準(zhǔn)則[10,11]
飛輪完整性分析包括應(yīng)力強(qiáng)度分析和斷裂力學(xué)分析。應(yīng)力強(qiáng)度分析包括在額定轉(zhuǎn)速(1500r/min)下,飛輪一次應(yīng)力不應(yīng)超過最小屈服強(qiáng)度的1/3;設(shè)計轉(zhuǎn)速(1875r/min)下,飛輪一次應(yīng)力不應(yīng)超過最小屈服強(qiáng)度的2/3。斷裂力學(xué)分析是飛輪在有最大假想裂紋的工況下,判斷高應(yīng)力區(qū)裂紋處的應(yīng)力強(qiáng)度因子是否小于臨界應(yīng)力強(qiáng)度因子。
2.1 飛輪應(yīng)力強(qiáng)度解析法分析
2.1.1 飛輪分析模型
飛輪應(yīng)力強(qiáng)度解析法采用沿軸橫截面建立的軸對稱模型進(jìn)行分析,如圖2所示,r1,r2分別為輪轂的內(nèi)外徑;r3為貧鈾盤外徑。
2.1.2 飛輪貧鈾盤應(yīng)力強(qiáng)度解析分析
在設(shè)計轉(zhuǎn)速和額定轉(zhuǎn)速工況下,貧鈾盤自身旋轉(zhuǎn)產(chǎn)生的離心力和外部壓力綜合作用下產(chǎn)生的平均環(huán)向應(yīng)力:
設(shè)計轉(zhuǎn)速計算公式[12]:
額定轉(zhuǎn)速計算公式[12]:
圖2 貧鈾飛輪軸對稱分析模型
2.1.3 飛輪輪轂應(yīng)力強(qiáng)度解析分析
在設(shè)計轉(zhuǎn)速和額定轉(zhuǎn)速工況下,輪轂自身旋轉(zhuǎn)產(chǎn)生的離心力產(chǎn)生的平均環(huán)向應(yīng)力:
設(shè)計轉(zhuǎn)速計算公式[12]:
額定轉(zhuǎn)速計算公式[12]:
2.1.4 應(yīng)力強(qiáng)度解析法分析結(jié)果
按照上述分析計算方法對飛輪貧鈾盤及輪轂在設(shè)計轉(zhuǎn)速及額定轉(zhuǎn)速下的一次應(yīng)力強(qiáng)度進(jìn)行計算,結(jié)果見表2。
表2 飛輪一次應(yīng)力解析分析結(jié)果
2.2 飛輪應(yīng)力強(qiáng)度有限元法分析
采用有限元分析程序ANSYS 14.5對貧鈾飛輪結(jié)構(gòu)進(jìn)行強(qiáng)度計算分析,主要計算了設(shè)計轉(zhuǎn)速和額定轉(zhuǎn)速工況下貧鈾盤和輪轂的應(yīng)力。
2.2.1 飛輪有限元分析模型
在建立模型的過程中各部件均采用實體模型,網(wǎng)格由六面體單元(SOLID186)組成,根據(jù)結(jié)構(gòu)及其載荷的對稱性,模型取整體軸向結(jié)構(gòu)的1/9進(jìn)行分析,計算模型如圖4、圖5所示。
圖4 貧鈾盤有限元分析模型
圖5 輪轂有限元分析模型
2.2.2 應(yīng)力強(qiáng)度有限元分析結(jié)果
利用有限元分析軟件對飛輪貧鈾盤及輪轂在設(shè)計轉(zhuǎn)速及額定轉(zhuǎn)速下的應(yīng)力強(qiáng)度進(jìn)行計算,設(shè)計轉(zhuǎn)速下貧鈾盤和輪轂的應(yīng)力強(qiáng)度云圖如圖6和圖7所示,結(jié)果見表3。
圖6 貧鈾盤有限元分析結(jié)果
圖7 輪轂有限元分析結(jié)果
表3 飛輪一次應(yīng)力有限元分析結(jié)果
分別采用解析分析方法與有限元分析方法計算了飛輪貧鈾盤與輪轂在設(shè)計轉(zhuǎn)速及額定轉(zhuǎn)速工況下的一次應(yīng)力,并依據(jù)飛輪應(yīng)力設(shè)計準(zhǔn)則對分析結(jié)果進(jìn)行了評定,結(jié)果證明飛輪一次應(yīng)力均滿足強(qiáng)度要求。同時對解析分析方法與有限元分析方法得到的結(jié)果進(jìn)行了對比,計算結(jié)果基本一致,證明了兩種方法的可行性。飛輪應(yīng)力分析結(jié)果對比與評定見表4。
表4 飛輪應(yīng)力分析結(jié)果對比與評定 單位:MPa
方法設(shè)計轉(zhuǎn)速工況額定轉(zhuǎn)速工況 貧鈾盤一次應(yīng)力輪轂一次應(yīng)力貧鈾盤一次應(yīng)力輪轂一次應(yīng)力 計算值限值評定計算值限值評定計算值限值評定計算值限值評定 有限元法130266.6合格1581.6合格78133.3合格9.8163.2合格 解析法132.5266.6合格15.581.6合格78.8133.3合格9.9163.2合格 誤差%1.93.311
根據(jù)線彈性斷裂力學(xué)分析原理,首先假定在飛輪最危險部位存在一個最大的假想裂紋,并計算該部位的應(yīng)力強(qiáng)度因子,然后將值與臨界應(yīng)力強(qiáng)度因子相比較,判斷是否能防止延性斷裂。
4.1 裂紋的位置和大小
最大假想缺陷為垂直于最大應(yīng)力方向(縱向)的尖銳表面缺陷。對于貧鈾盤鍵槽處的環(huán)向拉應(yīng)力是最大的,為305MPa,因此假設(shè)裂紋的位置在鍵槽的應(yīng)力最大點位置?,F(xiàn)按飛輪內(nèi)圓存在一個橢圓形表面裂紋進(jìn)行計算,初始裂紋尺寸取為預(yù)期漏檢缺陷的最大尺寸i=1.6mm。
4.2 應(yīng)力強(qiáng)度因子計算[14~15]
將計算對象簡化成半無限大體具有表面半橢圓形片狀裂紋的計算模型,其應(yīng)力強(qiáng)度因子表達(dá)式為:
式中E()是以為參數(shù)的第二類完整橢圓積分
式中:a——半橢圓裂紋的短半軸;
C——半橢圓裂紋的長半軸。
當(dāng)<<時,/0,有=1(相當(dāng)于單邊貫穿裂紋)
應(yīng)力強(qiáng)度因子
本文對貧鈾飛輪結(jié)構(gòu)的完整性進(jìn)行了分析,包括應(yīng)力強(qiáng)度分析和斷裂力學(xué)分析,得到以下結(jié)論:
(1)在應(yīng)力強(qiáng)度分析中采用解析分析方法和有限元分析方法計算了貧鈾飛輪貧鈾盤及輪轂在設(shè)計轉(zhuǎn)速和額定轉(zhuǎn)速下的一次應(yīng)力,一次應(yīng)力值均小于設(shè)計準(zhǔn)則中的應(yīng)力限值,滿足應(yīng)力強(qiáng)度要求。
(2)將解析法分析結(jié)果與有限元法分析結(jié)果對比,得出兩種方法計算結(jié)果基本一致的結(jié)論,證明飛輪一次應(yīng)力分析結(jié)果是合理可信的。
(3)斷裂力學(xué)分析是在貧鈾盤鍵槽位置假想有一個橢圓形表面裂紋,計算該部位的應(yīng)力強(qiáng)度因子,小于臨界應(yīng)力強(qiáng)度因子,并有一定的余量。
[1] 李夢啟.三代核電反應(yīng)堆冷卻劑泵屏蔽電機(jī)的技術(shù)特點[J]. 大電機(jī)技術(shù), 2009, 8(10): 1-4.
[2] 陳興江, 羅志遠(yuǎn)等. 核電站軸封型主泵電機(jī)飛輪結(jié)構(gòu)分析與研究[J]. 大電機(jī)技術(shù), 2015, (241): 30-33.
[3] 陳剛. 核電屏蔽電機(jī)主泵大慣量飛輪技術(shù)[J]. 防爆電機(jī), 2016, 51(192): 23-26.
[4] RG 1.14-1975, 反應(yīng)堆冷卻劑泵飛輪完整性[S].
[5] 葛磊, 王德忠, 張繼革. 屏蔽電機(jī)主泵慣性飛輪完整性分析[J]. 上海理工大學(xué)學(xué)報, 2012, 34(6): 584-588.
[6] 李笑天, 鞠紅惠. 主循環(huán)泵慣性飛輪完整性分析[J]. 原子能科學(xué)技術(shù), 2007, 41(4): 1-6.
[7] 張貴濱. 冷卻劑泵電機(jī)飛輪完整性分析[J]. 民營科技, 2010(6):16-16.
[8] 張強(qiáng), 李昊, 張立梅. 柴油機(jī)飛輪的應(yīng)力分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報, 2002, 33(4):131-132.
[9] 崔海燕, 李天斌. 核主泵飛輪結(jié)構(gòu)及設(shè)計要求探討[J]. 通用機(jī)械, 2015(9):79-82.
[10] 上海核工程研究設(shè)計院.反應(yīng)堆冷卻劑泵設(shè)備規(guī)范書(屏蔽泵) Rev.1[R]. 上海: 2014.
[11] 宋艷文, 李再赫, 于令新,等. 核電站反應(yīng)堆冷卻劑泵電動機(jī)轉(zhuǎn)子飛輪的設(shè)計與評定[J]. 大電機(jī)技術(shù), 1999(5):1-3.
[12] 電指(DZ)-1965, 大型電機(jī)機(jī)械計算公式[S].
[13] 宋俊杰, 劉茂江譯. 材料力學(xué)手冊[M]. 河北: 人民出版社, 1982.
[14] 李慶芬等. 斷裂力學(xué)及其在工程上的應(yīng)用[M]. 哈爾濱:工程大學(xué)出版社, 2007.
[15] 程靳, 趙樹山. 斷裂力學(xué)[M]. 北京: 科學(xué)出版社, 2006.
Integrity Analysis on Depleted Uranium Flywheel of the Shield-Pump
ZHANG Guibin, ZHANG Zhenyu
(1. Engineering Research Center of RCP Heilongjiang Province;2. Harbin Electric Power Equipment Co., Ltd.,Harbin150066, China)
To avoid damage of reactor core in scram of the reactor or the event of loss of motive power,the flywheels have to be mounted on canned-motor pump to store potential energy mechanically. Simultaneously, shied pump flywheel is made out of a heavy depleted uranium alloy, which has higher inertia in confined space. The integrity analysis of depleted uranium flywheel is concerned in the aspects of structure strength and fracture mechanics. The analysis results show that the stress and fracture ductility are both satisfactory to the requirements of regulations under the rated speed as well as the design speed. It can be concluded that the result is credible due to the consistency between the results of analytical method and finite element method.
pump canned-motor; depleted uranium flywheel; stress intensity; fracture mechanics
TM343
A
1000-3983(2017)02-0047-04
2016-10-13
張貴濱(1981-),男,2005年畢業(yè)于吉林大學(xué)機(jī)械專業(yè),現(xiàn)從事核主泵、大中型交直流電動機(jī)結(jié)構(gòu)強(qiáng)度計算方面的工作,工程師。
審稿人:李志和