李 偉,周錄坤,袁 禹,方奇術(shù),喬彥龍
三代核電廠卸壓系統(tǒng)管線冰塞模擬和材料損傷特性研究
李偉1,周錄坤1,袁禹1,方奇術(shù)2,*,喬彥龍2
(1. 山東核電有限公司,山東 海陽,216600;2. 國家電投集團(tuán)電站運(yùn)營技術(shù)(北京)有限公司,北京 昌平 112209)
為驗(yàn)證三代核電廠大口徑卸壓系統(tǒng)管線實(shí)施冰塞作業(yè)的安全可靠性,模擬管線的實(shí)際工況對(duì)冰塞形成過程進(jìn)行試驗(yàn)驗(yàn)證分析,并在評(píng)估反復(fù)冰塞作業(yè)后,對(duì)管材損傷進(jìn)行試驗(yàn)研究??煽啃则?yàn)證過程包括有限元模擬計(jì)算和理化試驗(yàn)兩部分內(nèi)容?;贏NSYS軟件的運(yùn)用,對(duì)核電站卸壓系統(tǒng)管線冰塞形成過程中管道溫度場進(jìn)行模擬計(jì)算,進(jìn)而得出冰塞形成的確切時(shí)刻,以及某時(shí)域內(nèi)冰塞厚度的變化情況,以保證冰塞形成后具有一定的強(qiáng)度。理化試驗(yàn)主要分為金相組織檢驗(yàn)、顯微硬度測(cè)試、夏比沖擊試驗(yàn)、常溫拉伸試驗(yàn)和掃描電鏡斷口分析。并根據(jù)上述結(jié)果分析,對(duì)現(xiàn)有的冰塞作業(yè)進(jìn)行優(yōu)化建議,對(duì)核電廠冰塞作業(yè)程序進(jìn)行評(píng)估改進(jìn)。
核電廠;有限元冰塞模擬;材料損傷;焊接作業(yè);
在核電廠運(yùn)行維護(hù)中,設(shè)備如故障或預(yù)防性維修需解體或更換,必須對(duì)系統(tǒng)進(jìn)行隔離,經(jīng)常存在系統(tǒng)的無隔離設(shè)備,或隔離設(shè)備不可用,造成系統(tǒng)不能隔離或者由于安全和經(jīng)濟(jì)上的原因,如系統(tǒng)介質(zhì)量大、疏排和充裝時(shí)間長;放射性外逸等,介質(zhì)疏排難度大或難以滿足運(yùn)行要求,這些情況下,采用冰塞來進(jìn)行系統(tǒng)隔離是一種可行的方法。冰塞技術(shù)雖然在國外已在使用,但在國內(nèi)也只有少數(shù)核電廠使用以及部分火電廠中,如秦山、田灣、福清等。伴隨著國內(nèi)三代核電機(jī)組建成與投入運(yùn)營,冰塞技術(shù)應(yīng)用在新機(jī)組逐步從小管徑向大管徑延伸[1-5]。
AP1000核電廠引入了安全系統(tǒng)的非能動(dòng)概念,在設(shè)計(jì)中采用非能動(dòng)的理念預(yù)防嚴(yán)重事故,大大降低了人因失誤的可能性。三代核電廠自動(dòng)卸壓系統(tǒng)(簡稱ADS)是AP1000核電廠非能動(dòng)系統(tǒng)中的典型系統(tǒng),在嚴(yán)重事故工況下,能提供快速卸壓以防止堆芯融化。每個(gè)機(jī)組共設(shè)置4列,ADS管線上串聯(lián)有爆破閥和電動(dòng)隔離閥各一臺(tái),爆破閥的設(shè)計(jì)可保證事故條件下能夠可靠地開啟,且不會(huì)出現(xiàn)誤關(guān)閉。根據(jù)維修大綱要求,爆破閥每10年需要進(jìn)行解體檢修,需要將爆破閥從ADS系統(tǒng)管線拆除,由于機(jī)組停堆后該管段上方水位較低,管線及爆破閥內(nèi)部的冷卻劑無法自行排出[6]。
采用液氮法將爆破閥與電動(dòng)隔離閥之間的冷卻劑制備成冰塞,避免在拆卸時(shí)爆破閥液態(tài)冷卻劑外泄造成設(shè)備與人員的沾污,是實(shí)現(xiàn)檢修必要條件。國內(nèi)冰塞實(shí)施經(jīng)驗(yàn)表明冰塞作業(yè)沒有詳細(xì)的指導(dǎo)書對(duì)冰塞制作的整個(gè)過程給予系統(tǒng)準(zhǔn)確的指導(dǎo),作業(yè)人員對(duì)冰塞形成規(guī)律認(rèn)識(shí)不足,冷凍介質(zhì)選用、冷凍介質(zhì)用量、冷凍時(shí)間都具有隨機(jī)性,都是冰塞作業(yè)風(fēng)險(xiǎn)來源[7]。為研究大口徑實(shí)施冰塞作業(yè)的安全可靠性,有必要模擬爆破閥管線的實(shí)際工況對(duì)冰塞形成過程進(jìn)行試驗(yàn)驗(yàn)證分析,并在評(píng)估反復(fù)冰塞作業(yè)后,對(duì)管材損傷進(jìn)行試驗(yàn)研究。并根據(jù)試驗(yàn)結(jié)果,對(duì)現(xiàn)有的冰塞夾套工具進(jìn)行優(yōu)化建議,對(duì)冰塞作業(yè)程序進(jìn)行評(píng)估與改進(jìn)[8,9]。
本次冰塞綜合評(píng)價(jià)試驗(yàn)主要針對(duì)特定尺寸304L不銹鋼材料??煽啃则?yàn)證過程包括模擬計(jì)算和理化試驗(yàn)兩部分內(nèi)容。模擬計(jì)算主要是基于ANSYS軟件的運(yùn)用對(duì)核電站冰塞形成過程中管道溫度場進(jìn)行理論分析,進(jìn)而得出冰塞形成的確切時(shí)刻,以及某時(shí)域內(nèi)冰塞厚度的變化情況,以保證冰塞形成后具有一定的強(qiáng)度,為冰塞技術(shù)進(jìn)一步的應(yīng)用研究提供一定的理論依據(jù)。
理化試驗(yàn)主要分為材料冷凍前后的金相組織檢驗(yàn)、顯微硬度測(cè)試、夏比沖擊試驗(yàn)、常溫拉伸試驗(yàn)、和掃描電鏡斷口分析。通過金相組織分析觀察材料,用于分析材料在常溫、低溫下,材料組織結(jié)構(gòu)是否發(fā)生變化;顯微硬度測(cè)試用于分析材料在常、低溫條件下的硬度分布及變化情況;夏比沖擊試驗(yàn)確定304L不銹鋼焊縫的韌脆轉(zhuǎn)變溫度;拉伸試驗(yàn)?zāi)康氖潜容^材料常溫、低溫力學(xué)性能的變化;掃描電鏡斷口分析目的是觀察材料的斷口類型[4]。
管道及管道內(nèi)部的水設(shè)為恒溫49 ℃,管道中間400 mm范圍處的管外壁,加載-30 ℃冷源代替冰塞夾套作用,模擬冰塞形成過程,并保證冰塞厚度至少到達(dá)400 mm,隨后分別在距冰塞40 mm和100 mm處施加2 500 ℃的溫度熱源,進(jìn)行焊接作業(yè)模擬,觀察其對(duì)冰塞溫度場的影響。
在材料熱變化過程中,各物理性能參數(shù)不斷隨溫度變化,因此,在模擬計(jì)算溫度場時(shí)要對(duì)材料的熱物理性能參數(shù)進(jìn)行定義。采用線性插值法和外推法求得材料在某一溫度區(qū)間內(nèi)某個(gè)時(shí)刻的性能參數(shù)。此次模擬所涉及的材料的熱物理性能參數(shù)如表1所示。
依據(jù)實(shí)際管道尺寸及焊縫分布建立焊件的幾何模型。幾何尺為219.1 mm/23.01 mm/1 000 mm(外徑/厚度/長度)。為了得到合理準(zhǔn)確地得到計(jì)算結(jié)果,在建立有限元計(jì)算模型時(shí),需要在ANSYS中選擇與材料結(jié)構(gòu)相對(duì)應(yīng)的單元類型進(jìn)行網(wǎng)格劃分。本文采用間接耦合法進(jìn)行有限元溫度場計(jì)算,在進(jìn)行溫度場計(jì)算時(shí)選擇具有4個(gè)節(jié)點(diǎn)的PLANE55面單元以及具有8個(gè)節(jié)點(diǎn)的SOLID70六面體單元分別對(duì)面和體進(jìn)行網(wǎng)格劃分,這兩種單元在每個(gè)節(jié)點(diǎn)上都只有一個(gè)溫度自由度。
圖1為大口徑管道在冰塞夾套作用下,形成冰塞并且冰塞厚度逐漸增加過程中縱截面溫度場分布云圖。圖1(a)為1 000 s時(shí),可以看出管道中心400 mm范圍,少部分管道內(nèi)水的溫度達(dá)0 ℃以下。圖1(b)和圖1(c)分別為2 000 s和3 000 s時(shí),管道縱截面溫度場分布云圖,可以看出,管道內(nèi)管壁附近冰層厚度逐漸增加,直至3 822 s時(shí)形成冰塞,溫度分布由中心向兩端逐漸增加。圖1(d)為4 980 s時(shí),冰塞厚度達(dá)到400 mm,呈扁沙漏狀,溫度大部分為-10 ℃以下,冰塞兩端部分約為0 ℃,其整體溫度場分布云圖如圖2所示。
圖3為管道冰塞橫截面溫度場分布云圖,在4 980 s冰塞形成時(shí),取冰塞中心位置橫截面溫度場,可以看出溫度由中心向四周逐漸降低,最高溫度約為-1 ℃,最低溫度為-30 ℃。冰塞厚度達(dá)400 mm時(shí),在冰塞一端(距離中心位置200 mm)處取橫截面。
通過圖4可以看出,計(jì)算溫度及環(huán)境溫度設(shè)為30 ℃時(shí),在5 000 s時(shí),冰塞形成厚度達(dá)400 mm,7 000 s時(shí),冰塞厚度區(qū)域穩(wěn)定,冰塞區(qū)域溫度直到10 000 s時(shí)變化不大,因此,7 000 s時(shí),冰塞趨于穩(wěn)定。
計(jì)算溫度及環(huán)境溫度設(shè)為57 ℃,進(jìn)行加載計(jì)算,得出以下云圖及中心軸線溫度分布曲線與上述情況類似,通過溫度場分布云圖及溫度分布曲線圖可以得出,8 000 s時(shí),冰塞形成厚度達(dá)400 mm,9 000 s時(shí),冰塞厚度區(qū)域穩(wěn)定,冰塞區(qū)域溫度較10 000 s時(shí)變化不大,因此,計(jì)算溫度及環(huán)境溫度設(shè)為57 ℃時(shí),在9 000 s時(shí),冰塞趨于穩(wěn)定。
冰塞形成后,在距冰塞40 mm處進(jìn)行焊接工況模擬,整體溫度場分布云圖如圖5所示。冰塞形成后,從第6 000 s開始進(jìn)行焊接熱源加載,焊接速度為2.3 mm/s,焊接總時(shí)長為300 s,從6 000 s時(shí)開始進(jìn)行焊接。6 150 s時(shí),焊接過程進(jìn)行一半,靠近焊接處的冰塞溫度場發(fā)生變化,受焊接過程影響,呈非對(duì)稱分布,如圖5(a)和圖5(c)所示。6 300 s時(shí),焊接過程完成,管壁溫度升高,使管壁附近冰塞少部分融化,冰塞右端靠近焊接處呈外凸?fàn)?,如圖5(d)所示。
圖1 冰塞形成過程縱截面溫度場云圖
圖2 冰塞厚度達(dá)400 mm時(shí)管道溫度場云圖
圖3 t=4 980 s冰塞橫截面溫度場分布云圖
圖4 中心軸線溫度分布曲線圖
在距冰塞100 mm處進(jìn)行焊接工況模擬時(shí),管道及冰塞溫度場分布特征相似。由于焊接位置相對(duì)較遠(yuǎn),因此對(duì)冰塞的加熱作用減小,冰塞融化量降低。
取兩種焊縫距離模式下的冰塞中心點(diǎn)溫度及時(shí)間數(shù)據(jù),作溫度隨時(shí)間變化曲線圖,如圖6所示。通過比較,在6 000 s施加焊接后,中心點(diǎn)受冰塞夾套作用溫度繼續(xù)降低,6 300 s時(shí),焊接完成,在焊接完成后的200 s內(nèi),溫度依然降低,直至6 564 s時(shí),冰塞中心點(diǎn)溫度開始上升,隨焊接熱輸入的傳遞,溫度不斷上升,計(jì)算至9 300 s,冰塞中心點(diǎn)溫度約為-2.2 ℃,冰塞仍未完全融化??梢钥闯鰞煞N焊縫分布情況進(jìn)行焊接工況模擬,溫度變化區(qū)別甚微。
圖5 距冰塞40 mm處焊接模擬溫度場分布云圖
圖6 冰塞中心溫度變化曲線圖
金相分析依據(jù)GB/T 13298—1991《金屬顯微組織檢驗(yàn)方法》進(jìn)行,將樣品置于顯微鏡下觀察,分析材料宏觀和微觀組織結(jié)構(gòu),分析冰凍對(duì)組織造成的變化。
如圖7、圖8所示,通過冷凍前后304L不銹鋼焊接試樣母材、焊縫中心區(qū)兩個(gè)典型部位的進(jìn)行組織分析,對(duì)比其金相組織可以發(fā)現(xiàn),晶粒形態(tài)、大小、組織基本類似,可以確認(rèn),低溫冷凍未對(duì)304L不銹鋼焊接試樣的金相組織產(chǎn)生明顯的影響。
采用顯微硬度儀,載荷100 g,加載時(shí)間15 s,對(duì)冷凍前和冷凍后的材料進(jìn)行硬度測(cè)試,結(jié)果如表2所示。
圖7 未冷凍的304L不銹鋼焊接試樣不同部位的金相組織
圖8 冷凍的304L不銹鋼焊接試樣不同部位的金相組織
圖8 冷凍的304L不銹鋼焊接試樣不同部位的金相組織(續(xù))
表2 冷凍前后試樣不同部位的硬度
通過未冷凍和冷凍后304L不銹鋼焊接試樣母材、焊縫中心兩個(gè)典型部位的顯微硬度的分析,冷凍前后硬度變化不大,冷凍未對(duì)304L不銹鋼焊接試樣的顯微硬度產(chǎn)生明顯的影響。
每種材料分別在焊縫區(qū)和母材部位取樣。試驗(yàn)時(shí)從0 ℃開始,對(duì)于每類樣品每隔5 ℃進(jìn)行一次沖擊試驗(yàn),測(cè)定材料的沖擊吸收能量,最后擬合出材料的沖擊能量和溫度的關(guān)系曲線,從而確定出材料的韌脆轉(zhuǎn)變溫度。試驗(yàn)最低溫度到-60 ℃。
如圖9所示,從沖擊功的數(shù)據(jù)來看,在試驗(yàn)范圍內(nèi),304L不銹鋼母材的沖擊功無明顯變化,表明304L不銹鋼母材在試驗(yàn)溫度的范圍內(nèi)沒有發(fā)生韌脆轉(zhuǎn)變;從宏觀形貌上,0 ℃到-60 ℃的試樣沖擊后斷口端都有明顯的塑性變形,呈梯形斷面,斷口末端有剪切唇,印證了304L不銹鋼母材在試驗(yàn)范圍內(nèi)未發(fā)生韌脆轉(zhuǎn)變。
圖9 304L不銹鋼母材的沖擊功隨溫度變化曲線
從微觀形貌上,從0 ℃一直到-60 ℃,304L不銹鋼母材試樣沖擊后斷口上都存在明顯的韌窩,這是塑性斷裂的典型特征。進(jìn)一步確認(rèn)了在試驗(yàn)范圍內(nèi),304L不銹鋼母材未發(fā)生韌脆轉(zhuǎn)變(見圖10)。
如圖11所示從沖擊功的數(shù)據(jù)來看,在試驗(yàn)范圍內(nèi),304L不銹鋼焊縫中心部分的沖擊功無明顯變化,表明304L不銹鋼焊縫中心在試驗(yàn)溫度的范圍內(nèi)沒有發(fā)生韌脆轉(zhuǎn)變;但比較而言,304L不銹鋼焊縫中心的沖擊功低于母材的沖擊功,表明焊縫中心的韌性比母材差。從宏觀形貌上,0 ℃到-60 ℃的試樣沖擊后斷口端都有明顯的塑性變形,呈梯形斷面,斷口末端有剪切唇,印證了304L不銹鋼焊縫中心在試驗(yàn)范圍內(nèi)未發(fā)生韌脆轉(zhuǎn)變。
圖10 304L不銹鋼母材的斷口微觀形貌隨溫度的變化
圖11 304L不銹鋼焊縫中心部分的沖擊功隨溫度變化曲線
如圖12所示,從微觀形貌上,從0 ℃一直到-60 ℃,304L不銹鋼焊縫中心試樣沖擊后斷口上都存在明顯的韌窩,這是塑性斷裂的典型特征。進(jìn)一步確認(rèn)了在試驗(yàn)范圍內(nèi),304L不銹鋼焊縫中心未發(fā)生韌脆轉(zhuǎn)變。
按照GB/T 2651《焊接接頭拉伸試驗(yàn)方法》加工試樣,進(jìn)行室溫和低溫拉伸性能試驗(yàn)。從304L不銹鋼母材試樣的常溫和-40 ℃低溫拉伸結(jié)果來看,試樣的部位斷裂有一定隨機(jī)性,從斷口微觀形貌看均表現(xiàn)出塑性斷裂的特征。
如圖13所示,低溫拉伸性能試驗(yàn)依據(jù)GB/T 2651《焊接接頭拉伸試驗(yàn)方法》,測(cè)量材料的拉伸強(qiáng)度、屈服強(qiáng)度、延伸率及斷面收縮率。常溫拉伸和低溫拉伸(-40 ℃)試驗(yàn)所得結(jié)果與304L不銹鋼的規(guī)定值進(jìn)行比較,評(píng)價(jià)其拉伸性能優(yōu)劣。由于焊接件在組織和性能上都存在不均勻性,屈服強(qiáng)度s、斷后伸長率(%)、斷面收縮率(%)的數(shù)值僅供參考。
圖12 304L不銹鋼焊縫中心的斷口微觀形貌隨溫度的變化
圖13 304L不銹鋼母材試樣常溫和-40 ℃低溫拉伸后的宏觀形貌和斷口微觀形貌
從表3的測(cè)試結(jié)果可以看出,不管是常溫拉伸還是低溫拉伸,兩種材料母材和焊接試樣的拉伸強(qiáng)度保持在拉伸性能要求的范圍內(nèi),母材的屈服強(qiáng)度、斷后伸長率、斷面收縮率也接近或達(dá)到拉伸性能要求,焊接試樣的屈服強(qiáng)度、斷后伸長率、斷面收縮率也接近或達(dá)到拉伸性能要求。但是從試驗(yàn)結(jié)果來看,-40 ℃時(shí)材料的拉伸強(qiáng)度和屈服強(qiáng)度均大于對(duì)應(yīng)材料的常溫拉伸強(qiáng)度的屈服強(qiáng)度,說明隨溫度降低,材料的抗拉性能增強(qiáng)。焊縫區(qū)域的斷后伸長率較母材有所降低,說明焊縫的存在對(duì)材料的塑性變形能力產(chǎn)生一定影響,但仍基本滿足要求。
表3 304L拉伸性能要求及-40 ℃下的測(cè)試結(jié)果
通過有限元數(shù)值模擬計(jì)算以及對(duì)304 L材料冷凍前后理化試驗(yàn),得出以下結(jié)果:
(1)通過有限元數(shù)值模擬,能夠模擬冰塞形成時(shí)間及溫度分布,對(duì)實(shí)際作業(yè)中冷凍介質(zhì)用量確定、冷凍時(shí)間控制具有重要的指導(dǎo)意義。
(2)冰塞作業(yè)區(qū)附近焊接或切割作業(yè)所引起熱輸入,對(duì)形成的冰塞影響較小,不會(huì)對(duì)形成的冰塞造成較大影響。
(3) 304L不銹鋼鋼材料母材、焊縫中心的沖擊功在冷凍前后均沒有發(fā)生明顯突變,即沒有發(fā)生韌脆轉(zhuǎn)變,母材的沖擊韌性要好于焊縫區(qū)域的沖擊韌性,焊縫區(qū)域更容易發(fā)生塑性斷裂。
(4)多次冷凍試驗(yàn)發(fā)現(xiàn),304L材料的拉伸強(qiáng)度和屈服強(qiáng)度完全滿足材料性能要求,說明實(shí)施冰塞時(shí)可在同一部位進(jìn)行多次冰塞,而不會(huì)對(duì)材料產(chǎn)生很大影響[10, 11]。
[1] 徐昶,王沁宇,肖朱強(qiáng),等. 冰塞技術(shù)在福清核電項(xiàng)目的研究與應(yīng)用[J]. 能源與環(huán)境,2015(04):27-29.
[2] 彭峰. 淺談冰塞作業(yè)管理模式[J]. 科技創(chuàng)新與應(yīng)用,2013(3):157.
[3] 朱友學(xué). 冰塞技術(shù)在核電廠的應(yīng)用[J]. 黑龍江科學(xué),2015(06):7.
[4] 杜軍林. 謝海濤. 冰塞法在600 MW火電機(jī)組上的應(yīng)用[J]. 科學(xué)技術(shù)創(chuàng)新,2018(32):66-67.
[5] 王長征,田亞,王保田. 冰塞法處理600MW機(jī)組密封油系統(tǒng)異常問題的探討[A]. 全國火電大機(jī)組(300 MW級(jí))競賽第36屆年會(huì)論文集(上冊(cè))[C],2007.
[6] 花龍. PCS系統(tǒng)安全殼噴淋隔離閥維修可執(zhí)行性分析[J]. 山東工業(yè)技術(shù),2019(09):194.
[7] 王賀宇. 冰塞技術(shù)在核電廠應(yīng)用的理論研究[J]. 中國核工業(yè),2008(05):136-141.
[8] 劉媛. 冰塞行成過程中管道溫度場理論研究[J]. 科技資訊,2012(06):242.
[9] 謝林君,沈昕,林韓波,等. 工業(yè)管道液氮冰塞試驗(yàn)及其應(yīng)力分布研究[J]. 低溫工程,2020(6):38-43.
[10] G J Howard,蔡伯淳,黃志剛. 管道試壓時(shí)冷凍封堵技術(shù)取得的進(jìn)展[J]. 低溫與特氣,1983(3):49-51.
[11]楊斐. 冰塞技術(shù)應(yīng)用風(fēng)險(xiǎn)與對(duì)策[J]. 中國核電,2012(2):183-189.
Simulation of Ice Plugs in Pipelines of Pressure Relief Systems in Third-generation Nuclear Power Plants and Study on Material Damage Characteristics
LI Wei1,ZHOU Lukun1,YUAN Yu,F(xiàn)ANG Qishu2,*,QIAO Yanlong2
(1. Shandong Nuclear Power Co.,Ltd.,Yantai of Shandong Prov. 265116,China; 2. SPIC Power Station Operation Technology(Beijing)Co.,Ltd.,Beijing 112209,China)
In order to verify the safety and reliability of the ice plug operation in the large-diameter pressure relief system pipeline of the third-generation nuclear power plant, the actual working conditions of the pipeline are simulated to carry out experimental verification and analysis of the ice plug formation process. After repeated ice plug operations are evaluated, the experimental study on the damage of the pipe is carried out. The reliability verification process includes two parts: finite element simulation calculation, and physical and chemical tests. Based on the application of the ANSYS software, the temperature field of the pipeline during the formation of the ice plug in the pressure relief system of the nuclear power plant is simulated and calculated, and then the exact moment of the ice plug formation and the change of the thickness of the ice plug in a certain time domain are obtained to ensure the formation of the ice plug. Afterwards, it has a certain strength. The physical and chemical tests include microstructure examination, microhardness test, Charpy impact test, normal temperature tensile test, metallographic structure inspection, and scanning electron microscope fracture analysis. Based on the analysis of the above results, optimization suggestions are made for the existing ice jam operation, and the ice jam operation procedure of the nuclear power plant is evaluated and improved.
Nuclear power plant; Finite element ice plug simulation; Material damage; Welding operations
TL48
A
0258-0918(2023)05-1049-10
2022-08-17
國家科技重大專項(xiàng)資助“核電站Living-PSA 和在線風(fēng)險(xiǎn)監(jiān)測(cè)與管理技術(shù)研究”(2014ZX06004-003)
李 偉(1983—),男,山東濱州人,中級(jí)工程師,碩士研究生,現(xiàn)主要從事核電廠管閥類設(shè)備檢修方面研究
方奇術(shù),E-mail:fangqishu@spic.com.cn