1 引言
中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)專項(xiàng)TMSR（Thorium-based Molten Salt Reactor，釷基熔鹽堆）是當(dāng)今世界第四代先進(jìn)反應(yīng)堆的六個(gè)候選堆型之一，具有中子經(jīng)濟(jì)性好、放射性廢物少、功率密度高和安全性較高等優(yōu)點(diǎn)。熔鹽泵是回路和反應(yīng)堆核島的關(guān)鍵性動(dòng)力旋轉(zhuǎn)設(shè)備，在超長(zhǎng)服役的過程中，它必須能承受各種極端的事故工況。熔鹽泵的結(jié)構(gòu)力學(xué)性能不僅影響著自身的安全運(yùn)行，同時(shí)也關(guān)系著整個(gè)反應(yīng)堆回路系統(tǒng)乃至反應(yīng)堆的正常運(yùn)行。根據(jù)ASME規(guī)范和我國(guó)的核安全法規(guī)，熔鹽泵除了要滿足其自身的功能要求，還要求能承受一定的地震載荷。具體的要求是熔鹽泵在承受OBE（運(yùn)行基準(zhǔn)地震）擾動(dòng)周期內(nèi)能繼續(xù)運(yùn)轉(zhuǎn)而不發(fā)生機(jī)械損壞，同時(shí)也不會(huì)引起其他設(shè)備的損壞，在SSE（安全停堆地震）下能保持承壓邊界的完整性。
本文以現(xiàn)有TMSR回路用熔鹽泵為基礎(chǔ)，對(duì)含轉(zhuǎn)子部件在內(nèi)的熔鹽泵進(jìn)行三維有限元建模，對(duì)熔鹽泵的固有頻率和在兩類地震工況下的抗震性能進(jìn)行了相關(guān)研究。本文的研究成果為下一步核級(jí)泵的抗震設(shè)計(jì)和安全評(píng)估提供參考依據(jù)。
2 有限元模型
2.1 三維實(shí)體模型的建立
由于熔鹽泵模型較為復(fù)雜，故一般三維實(shí)體的建模工作均不考慮直接在有限元軟件中完成，而是借助建模功能更為強(qiáng)大的三維造型軟件，常用的三維造型軟件有NX、Pro/Engineer、Solid Edge、Solidworks等。TMSR工程回路用熔鹽泵初始的三維實(shí)體軟件由回路設(shè)計(jì)部門在Solid Edge中完成，其幾何模型如圖1所示。熔鹽泵模型由泵罐、連接架、轉(zhuǎn)子組件、密封組件、葉輪、泵體和液下連接架等構(gòu)成，其中泵轉(zhuǎn)子組件包含兩對(duì)滾珠軸承。電機(jī)模型包含電機(jī)和電機(jī)座。
圖1 熔鹽泵幾何模型
2.2 材料特性
TMSR當(dāng)前工程回路熔鹽泵的運(yùn)行介質(zhì)為FLiNaK熔鹽，由于熔鹽介質(zhì)的特殊性，熔鹽泵的主體結(jié)構(gòu)均采用能耐熔鹽腐蝕的C276合金，電機(jī)座和電機(jī)由于并不與熔鹽相接觸，使用材料為Incnel 316。此外，TMSR設(shè)計(jì)回路系統(tǒng)的運(yùn)行溫度在650℃以上，隨著溫度升高，材料的力學(xué)特性也會(huì)隨之變化，由于本文的計(jì)算并不涉及到熱瞬態(tài)計(jì)算，故僅選取650℃下C276和Incnel 316的力學(xué)參數(shù)進(jìn)行后續(xù)計(jì)算，兩種材料在650℃運(yùn)行溫度下的力學(xué)特性參數(shù)如表1所示。
表1 材料在650℃下的力學(xué)特性參數(shù)
2.3 有限元模型及邊界條件
在泵機(jī)組的整體計(jì)算中，一般并不對(duì)完整的軸承進(jìn)行建模和分析，而是采取等效的方式模擬滾珠對(duì)軸承座和軸套的作用。本文對(duì)轉(zhuǎn)子部件進(jìn)行簡(jiǎn)化，兩對(duì)滾珠軸承均采用SKF軸承公司等效線性化處理后的剛度進(jìn)行計(jì)算，線性化剛度參數(shù)如表2所示。
表2 熔鹽泵軸承剛度
軸套與軸承座之間的相互作用由Ansys單元庫中Combin14彈簧單元模擬，在每對(duì)軸套與軸承座圓環(huán)上各創(chuàng)建四個(gè)Combin14單元，每個(gè)單元的剛度為軸承線性化剛度的1/2。同時(shí)為了有效控制計(jì)算的規(guī)模，葉片等非關(guān)鍵部件以附加質(zhì)量單元Mass21的形式添加到整個(gè)有限元模型中。此外，為了更進(jìn)一步減小計(jì)算規(guī)模，泵與電機(jī)支座、電機(jī)支座與電機(jī)均采用剛性連接。最終熔鹽泵在Workbench平臺(tái)中采用Solid187單元進(jìn)行實(shí)體部件的網(wǎng)格劃分，整個(gè)模型的節(jié)點(diǎn)數(shù)為588374，單元數(shù)為327585，其有限元模型如圖2所示。
圖2 熔鹽泵有限元模型
由于熔鹽泵的運(yùn)行溫度在650℃左右，高溫環(huán)境下熔鹽泵內(nèi)部組件的熱匹配、熱變形釋放和管道熱位移處理等均是一個(gè)亟待解決的問題，這其中也涉及到泵支撐平臺(tái)的優(yōu)化改進(jìn)。也正是基于此原因，本文并沒有將泵支撐平臺(tái)作為有限元模型建立，僅對(duì)熔鹽泵的外部支撐做等效簡(jiǎn)化處理，即將泵罐法蘭與支撐平臺(tái)螺栓連接處的自由度全部約束。
3 模態(tài)分析
本文采用Workbench中默認(rèn)的PCG LANCZOS進(jìn)行模態(tài)提取。計(jì)算得到的一階頻率為18.7Hz，圖3和圖4分別為熔鹽泵機(jī)組和內(nèi)部轉(zhuǎn)子組件的一階振型。從圖中可以看出一階振型主要為水平X方向擺動(dòng)，擺動(dòng)幅度最大的位置在轉(zhuǎn)子組件的葉輪處，且熔鹽泵機(jī)組整體均參與了一階振型。由此可知，熔鹽泵機(jī)組整體的一階頻率和轉(zhuǎn)子組件的一階頻率相近。圖5和圖6分別為熔鹽泵機(jī)組和轉(zhuǎn)子組件的二階振型，主要體現(xiàn)為轉(zhuǎn)子組件水平Y(jié)方向的擺動(dòng)。從圖中可見熔鹽泵除轉(zhuǎn)子以外的部件未參與二階振型。泵機(jī)組第三階和第四階的頻率和振型如表3所示，從表中我們可以看出熔鹽泵第三階振型仍主要為泵轉(zhuǎn)子組件的振動(dòng)，直到第四階才同時(shí)出現(xiàn)泵機(jī)組整體和轉(zhuǎn)子組件在水平Y(jié)方向的擺動(dòng)，其中最大振動(dòng)幅度出現(xiàn)在電機(jī)位置。根據(jù)前四階模態(tài)的特征可以看出熔鹽泵機(jī)組在水平X和Y兩個(gè)方向的剛度存在一定的差異性，這與實(shí)際的結(jié)構(gòu)特征也相符合。
圖3 熔鹽泵機(jī)組一階振型
圖4 轉(zhuǎn)子組件一階振型
圖5 熔鹽泵機(jī)組二階振型
圖6 轉(zhuǎn)子組件二階振型
熔鹽泵整體前十階頻率和對(duì)應(yīng)的振型特征如表3所示�？傮w來說，在頻率較低的區(qū)間主要為泵轉(zhuǎn)子組件與電機(jī)的擺動(dòng)，在頻率較高的區(qū)間主要為泵軸與電機(jī)的彎曲振動(dòng)。值得注意的是，泵機(jī)組仍然是以泵轉(zhuǎn)子的振動(dòng)為主，水平擺動(dòng)幅度最大值多出現(xiàn)于葉輪，彎曲振動(dòng)幅度最大值多為泵主軸。
表3 熔鹽泵前十階模態(tài)及振型特征
4 抗震分析
4.1 抗震分析方法
抗震分析的主要目的是驗(yàn)證熔鹽泵在OBE和SSE兩類地震工況后，是否能保證其結(jié)構(gòu)完整性和可運(yùn)行性。常用的抗震分析方法有等效靜力法、反應(yīng)譜法和時(shí)程分析法。
等效靜力分析方法要求結(jié)構(gòu)的一階頻率大于33Hz，而反應(yīng)譜法和時(shí)程分析法則對(duì)結(jié)構(gòu)的頻率無計(jì)算方法上的要求。其中反應(yīng)譜法既考慮了結(jié)構(gòu)動(dòng)力特性和地震動(dòng)特性之間的關(guān)系，又充分運(yùn)用了靜力理論，使得復(fù)雜機(jī)構(gòu)地震作用下的響應(yīng)計(jì)算變得簡(jiǎn)單可行，但是其計(jì)算不能包含非線性因素。而時(shí)程分析法則是真正意義上的動(dòng)力分析方法，可充分考慮地震波的特征因素，如振幅、持續(xù)時(shí)間等，且可以充分考慮結(jié)構(gòu)的非線性因素，但是由于計(jì)算量龐大，一般情況下均是作為反應(yīng)譜分析方法的補(bǔ)充計(jì)算。結(jié)合前文的計(jì)算結(jié)果，熔鹽泵的一階頻率為18.7Hz，小于33Hz，再綜合考慮反應(yīng)譜法和時(shí)程分析法各自的特征，本文最后采用反應(yīng)譜法對(duì)熔鹽泵的地震響應(yīng)進(jìn)行計(jì)算。
按照規(guī)范，選取的OBE和SSE地震波譜阻尼比分別為2%和4%，水平兩個(gè)方向考慮加速度相同，兩類地震波譜分別如圖7、圖8所示。對(duì)于本文的研究對(duì)象，水平方向?yàn)閄和Y方向，豎直方向?yàn)閆。
圖7 OBE反應(yīng)譜
圖8 SSE反應(yīng)譜
4.2 抗震分析結(jié)果
4.2.1 強(qiáng)度校核
本文采用SRSS方法進(jìn)行反應(yīng)譜振型組合。OBE和SSE工況下的應(yīng)力分布云圖分別如圖9、圖10所示。從圖中可見，最大應(yīng)力均出現(xiàn)在連接架位置，兩類工況由于頻譜趨勢(shì)接近，主要差異體現(xiàn)在加速度數(shù)值上，故整體的應(yīng)力分布云圖趨勢(shì)相近。為了更進(jìn)一步考察熔鹽泵內(nèi)部的應(yīng)力情況，本文選取了SSE工況下熔鹽泵內(nèi)部的應(yīng)力云圖為研究對(duì)象，如圖11所示。從圖中可見，除了連接架位置，在轉(zhuǎn)子組件的軸承處和主軸位置整體應(yīng)力水平較其他部位較高。但熔鹽泵的整體應(yīng)力水平較低，OBE反應(yīng)譜下對(duì)應(yīng)最大應(yīng)力為12.6MPa，SSE對(duì)應(yīng)最大應(yīng)力為21.2MPa。
圖9 熔鹽泵OBE下整體應(yīng)力云圖
圖10 熔鹽泵SSE下整體應(yīng)力云圖
圖11 熔鹽泵SSE下內(nèi)部應(yīng)力云圖
泵的完整性要求主要體現(xiàn)在電機(jī)、電機(jī)座、轉(zhuǎn)子部件、連接架和泵罐等部件的強(qiáng)度問題上，本文以上述部件為研究對(duì)象，其最大應(yīng)力值如表4所示。由表可知，在當(dāng)前的OBE和SSE工況下，熔鹽泵的整體應(yīng)力水平均較低，低于材料的許用應(yīng)力。
表4 熔鹽泵在OBE和SSE兩類地震作用下關(guān)鍵部件的最大等效應(yīng)力
4.2.2 靜動(dòng)部件間隙校核
泵在地震工況下的可運(yùn)行性準(zhǔn)則要求：轉(zhuǎn)子組件在兩類地震工況下，轉(zhuǎn)動(dòng)件和靜止件變形后仍然能夠無障礙運(yùn)轉(zhuǎn)，不發(fā)生摩擦、擠壓現(xiàn)象。
針對(duì)熔鹽泵的設(shè)計(jì)，轉(zhuǎn)子組件最易與周圍部件發(fā)生接觸、碰撞的位置集中在轉(zhuǎn)子組件地震工況下?lián)隙茸畲蟮娜~輪位置和主軸與密封原件配合的部位。表5為這兩個(gè)關(guān)鍵部位在地震工況下的位移響應(yīng)值。由于與葉輪有間隙要求的泵體和與主軸配合的密封原件地震工況下的位移值均較小，與轉(zhuǎn)子組件相差一個(gè)數(shù)量級(jí)，故可考慮忽略其位移量。根據(jù)TMSR熔鹽泵設(shè)計(jì)間隙和運(yùn)行間隙數(shù)據(jù)得知，葉輪與泵體的運(yùn)行間隙為2mm，軸與密封原件的間隙為0.8mm，表5中葉輪位置和軸與密封元件位置在X、Y、Z三個(gè)方向以及總的位移響應(yīng)值均小于各自的限值，可見在兩類地震工況下轉(zhuǎn)子組件能滿足運(yùn)行需求。
表5 熔鹽泵在OBE和SSE兩類地震作用下轉(zhuǎn)子組件的位移
5 結(jié)論
本文建立了熔鹽泵三維有限元實(shí)體模型，計(jì)算得到了熔鹽泵的前十階模態(tài)和振型。在抗震分析中，采用反應(yīng)譜分析方法計(jì)算了在OBE和SSE兩類地震工況下熔鹽泵各關(guān)鍵部件的應(yīng)力響應(yīng)和轉(zhuǎn)子組件的位移響應(yīng)。計(jì)算結(jié)果表明，在當(dāng)前的地震工況下，熔鹽泵機(jī)組的整體應(yīng)力水平較低，小于材料許用應(yīng)力；熔鹽泵轉(zhuǎn)子組件和靜止件的間隙基本能滿足變形要求。由于反應(yīng)堆具體的選址和該區(qū)域的地震反應(yīng)譜尚未確定，故本文的計(jì)算僅作方法上的探討，而最終反應(yīng)堆建址地區(qū)的反應(yīng)譜確定、此確定反應(yīng)譜下的抗震研究、結(jié)合溫度場(chǎng)的應(yīng)力評(píng)定和疲勞計(jì)算將是下一步熔鹽泵計(jì)算中最重要的工作。
核心關(guān)注：拓步ERP系統(tǒng)平臺(tái)是覆蓋了眾多的業(yè)務(wù)領(lǐng)域、行業(yè)應(yīng)用，蘊(yùn)涵了豐富的ERP管理思想，集成了ERP軟件業(yè)務(wù)管理理念，功能涉及供應(yīng)鏈、成本、制造、CRM、HR等眾多業(yè)務(wù)領(lǐng)域的管理，全面涵蓋了企業(yè)關(guān)注ERP管理系統(tǒng)的核心領(lǐng)域，是眾多中小企業(yè)信息化建設(shè)首選的ERP管理軟件信賴品牌。
轉(zhuǎn)載請(qǐng)注明出處：拓步ERP資訊網(wǎng)http://www.ezxoed.cn/
本文標(biāo)題：TMSR熔鹽泵抗震性能研究
本文網(wǎng)址：http://www.ezxoed.cn/html/solutions/14019313747.html
相關(guān)文章
