楊智榮,楊再明,高海波,陳祖志,林治國(guó)

(1. 中國(guó)特種設(shè)備檢測(cè)研究院 國(guó)家市場(chǎng)監(jiān)管總局特種設(shè)備安全與節(jié)能重點(diǎn)實(shí)驗(yàn)室，北京 100029;2. 武漢理工大學(xué)能源與動(dòng)力工程學(xué)院，湖北 武漢 430063)

0 引 言

在礦物能源枯竭的背景下，全球天然氣貿(mào)易活躍，市場(chǎng)越來(lái)越廣闊[1]。CNG罐車(chē)作為天然氣主流的運(yùn)輸渠道具有成本低、操作簡(jiǎn)單的特點(diǎn)[2]。但其具有較高的安全隱患，且大多由于氣瓶因腐蝕或老化所引發(fā)的事故。傳統(tǒng)方法是對(duì)氣瓶進(jìn)行定期檢驗(yàn)，但其檢驗(yàn)流程復(fù)雜，且不能動(dòng)態(tài)監(jiān)測(cè)氣瓶的狀態(tài)。對(duì)此較多專家在無(wú)損檢測(cè)的基礎(chǔ)上，提出了在線監(jiān)測(cè)，如文獻(xiàn)[3]介紹了超聲波檢測(cè)在高溫壓力容器在線監(jiān)測(cè)上的應(yīng)用，并提出應(yīng)注意的問(wèn)題和測(cè)試條件，文獻(xiàn)[4]也指出了壓力容器聲發(fā)射檢測(cè)的發(fā)展方向，即在線監(jiān)測(cè)、聲發(fā)射信號(hào)的模式識(shí)別和人工神經(jīng)網(wǎng)絡(luò)模式識(shí)別分析。然而大多方法都基于傳統(tǒng)的檢驗(yàn)方法，因此其監(jiān)測(cè)設(shè)備繁雜、監(jiān)測(cè)參數(shù)復(fù)雜，并不適合監(jiān)測(cè)移動(dòng)式設(shè)備，因此高效的檢測(cè)及監(jiān)測(cè)技術(shù)成為當(dāng)前研究的重點(diǎn)。

中國(guó)特種設(shè)備檢測(cè)研究院聯(lián)合浙江大學(xué)、合肥通用機(jī)械等29家單位共同申報(bào)的“十三五”國(guó)家重點(diǎn)研發(fā)計(jì)劃“公共安全風(fēng)險(xiǎn)防控與應(yīng)急技術(shù)裝備”專項(xiàng)項(xiàng)目“移動(dòng)式承壓類特種設(shè)備風(fēng)險(xiǎn)防控與治理關(guān)鍵技術(shù)研究”，并提出移動(dòng)式承壓設(shè)備風(fēng)險(xiǎn)監(jiān)管平臺(tái)監(jiān)測(cè)，可實(shí)時(shí)監(jiān)測(cè)氣瓶運(yùn)輸中的狀態(tài)，保證車(chē)輛安全性，本文主要針對(duì)其中氣瓶損傷監(jiān)測(cè)內(nèi)容進(jìn)行研究，通過(guò)有限元方法，模擬氣瓶在運(yùn)輸中的狀態(tài)，定量制造氣瓶損傷，監(jiān)測(cè)氣瓶運(yùn)輸中由于損傷導(dǎo)致的各參數(shù)變化，得出合適的監(jiān)測(cè)參數(shù)，并結(jié)合了仿真結(jié)果、案例統(tǒng)計(jì)以及現(xiàn)場(chǎng)氣瓶水壓實(shí)驗(yàn)結(jié)果得出合適的監(jiān)測(cè)位置，為之后現(xiàn)場(chǎng)實(shí)驗(yàn)以及平臺(tái)搭建做準(zhǔn)備。

1 在線監(jiān)測(cè)平臺(tái)

目前，由于智能傳感、物聯(lián)網(wǎng)等技術(shù)的發(fā)展，為壓力容器在線監(jiān)測(cè)技術(shù)提供了新的監(jiān)測(cè)方法。本文提出多傳感CGN氣瓶損傷在線監(jiān)測(cè)平臺(tái)，其構(gòu)成如圖1所示，主要包括數(shù)據(jù)采集模塊、診斷模塊、無(wú)線通信模塊以及電源模塊。其監(jiān)測(cè)流程為：由傳感器采集氣瓶的特征參數(shù)，之后將數(shù)據(jù)傳輸?shù)皆\斷模塊，通過(guò)對(duì)數(shù)據(jù)進(jìn)行處理、融合以及診斷，最后通過(guò)無(wú)線通信模塊將診斷結(jié)果傳輸?shù)揭苿?dòng)設(shè)備。其主要采用多傳感網(wǎng)絡(luò)技術(shù)、多傳感信息融合技術(shù)以及數(shù)據(jù)處理等。

圖1 在線監(jiān)測(cè)平臺(tái)架構(gòu)示意圖

本文主要對(duì)數(shù)據(jù)采集模塊進(jìn)行研究，選定合適的監(jiān)測(cè)參數(shù)與位置，通過(guò)有限元數(shù)值仿真方法對(duì)氣瓶運(yùn)輸過(guò)程中的狀況進(jìn)行模擬分析，并分別對(duì)無(wú)損以及有損氣瓶進(jìn)行特征參數(shù)監(jiān)測(cè)，從多方面考慮，得到監(jiān)測(cè)參數(shù)以及監(jiān)測(cè)位置。

2 無(wú)損氣瓶數(shù)值仿真計(jì)算

對(duì)大容積無(wú)縫鋼瓶進(jìn)行建模，分別對(duì)有損氣瓶以及無(wú)損氣瓶進(jìn)行結(jié)構(gòu)分析，得出氣瓶在給定靜力載荷作用下的響應(yīng)[5]，主要包括結(jié)構(gòu)分析、模態(tài)分析，為隨機(jī)振動(dòng)仿真做準(zhǔn)備。

2.1 無(wú)損氣瓶模型建立

本文以大容積無(wú)縫鋼瓶作為研究對(duì)象，其材料為4130X，材料性能如表1所示[6]。由于在運(yùn)輸過(guò)程中，每個(gè)氣瓶相互獨(dú)立，所以只需對(duì)單獨(dú)氣瓶建模分析，其尺寸如圖2所示。通過(guò)Soildworks進(jìn)行建模得出模型如圖3所示。

表1 4130X材料性能

圖2 長(zhǎng)管拖車(chē)氣瓶（單位：mm）

圖3 氣瓶模型圖

2.2 模型結(jié)構(gòu)分析

結(jié)構(gòu)分析主要作用是確定氣瓶結(jié)構(gòu)在給定靜力載荷下是否安全，其原理如下[7]：

式中：K——?jiǎng)偠认禂?shù)矩陣；

x——位移矢量；

F——力矢量。

分析前需對(duì)模型進(jìn)行網(wǎng)格劃分，本次劃分采用掃掠法生成單元為20節(jié)點(diǎn)的六面體單元（SOLID186），網(wǎng)格劃分結(jié)果如圖4所示。

圖4 部分網(wǎng)格劃分結(jié)果

設(shè)置邊界條件為：重力加速度（9.8 m /s2）、氣瓶?jī)?nèi)壓（20 MPa）、氣瓶?jī)啥斯潭s束[8]。

由于氣瓶應(yīng)存在氣閥，會(huì)在工作壓力下對(duì)瓶口產(chǎn)生軸向應(yīng)力，為簡(jiǎn)化模型將其直接等效為兩端瓶口表面軸向拉應(yīng)力（-9.79 MPa）[9]。其中軸向拉應(yīng)力計(jì)算公式如下：

式中：Ps——?dú)馄績(jī)?nèi)壓（20 MPa）；

Di——瓶口內(nèi)直徑；

Do——瓶口外直徑。

最終得出通過(guò)網(wǎng)格無(wú)關(guān)性驗(yàn)證[10]的結(jié)構(gòu)分析結(jié)果如圖5所示，由于極限載荷判斷氣瓶結(jié)構(gòu)失穩(wěn)的主要評(píng)定標(biāo)準(zhǔn)，本次采用極限載荷法對(duì)氣瓶安全性進(jìn)行評(píng)定[11]，由于氣瓶最大應(yīng)力（537.18 MPa）未超過(guò)材料屈服極限657 MPa，說(shuō)明在此載荷下氣瓶處于安全狀態(tài)。

圖5 結(jié)構(gòu)分析應(yīng)力云圖

3 有損氣瓶數(shù)值仿真計(jì)算

氣瓶的損傷類型大致分為10類，分別是裂紋、腐蝕、鼓包、凹陷、瓶口螺紋損傷、熱損傷、壁厚超標(biāo)、熱損傷、硬度超標(biāo)、其他缺陷（表面不平整、分層、結(jié)疤等）。根據(jù)氣瓶產(chǎn)生損傷的原因可將其歸為3類：在制造中產(chǎn)生的損傷、氣瓶在使用中的損傷以及其他損傷。對(duì)于制造中缺陷可以通過(guò)出廠時(shí)檢驗(yàn)篩選出來(lái)，而使用中產(chǎn)生缺陷，只能通過(guò)定期檢驗(yàn)得出，并不能保證氣瓶的安全性，因此本次在線監(jiān)測(cè)平臺(tái)建立主要針對(duì)氣瓶運(yùn)行中的損傷。

3.1 有損氣瓶模型建立

由GB 13004—1999《鋼質(zhì)無(wú)縫氣瓶定期檢驗(yàn)與評(píng)定》[12]可知?dú)馄亢穸仁呛饬繗馄康闹匾蛩?，而且氣瓶損傷部位厚度低于最小允許壁厚（設(shè)計(jì)壁厚的95%）時(shí)則氣瓶報(bào)廢，所以部分有損氣瓶損傷參數(shù)如表2所示，其損傷等級(jí)范圍為0.1～1.87 mm。

表2 損傷等級(jí)

3.2 隨機(jī)振動(dòng)仿真

對(duì)于長(zhǎng)管拖車(chē)各部件來(lái)說(shuō)，一般不會(huì)由于靜強(qiáng)度的問(wèn)題而破壞的，其破壞往往是在行駛過(guò)程中由于動(dòng)強(qiáng)度不足而破壞，或由于長(zhǎng)期受到周期性載荷的作用而引起疲勞失效。結(jié)合長(zhǎng)管拖車(chē)行駛的實(shí)際情況，更多的是動(dòng)態(tài)載荷是由于路面凹凸不平所產(chǎn)生的隨機(jī)振動(dòng)，因此本文模擬振動(dòng)的方法是隨機(jī)振動(dòng)仿真實(shí)驗(yàn)，進(jìn)而得出各參數(shù)與損傷程度的關(guān)系，并選出變化較為明顯的參數(shù)作為在線監(jiān)測(cè)平臺(tái)的監(jiān)測(cè)參數(shù)。

首先需確定隨機(jī)振動(dòng)實(shí)驗(yàn)激勵(lì)，對(duì)于車(chē)輛而言，行駛過(guò)程中主要產(chǎn)生兩類激勵(lì)，分別為離散事件激勵(lì)與路面隨機(jī)激勵(lì)，其中第二類表示各種等級(jí)的高速公路路面。結(jié)合長(zhǎng)管拖車(chē)行駛的實(shí)際情況，本文采取隨機(jī)激勵(lì)。由文獻(xiàn)[13]可知，車(chē)輛隨機(jī)激勵(lì)一般由路面不平度的空間頻率位移功率譜密度表示。其計(jì)算式如下所示：

式中：n——空間頻率,m?1；是波長(zhǎng) λ的倒數(shù)，且0.011 m?1

Gq(n0)——參考空間頻率n0下的路面功率譜密度值，稱為路面不平度系數(shù)，其取值如表3所示；

表3 路面8級(jí)各路面不平度系數(shù)取值

w——頻率指數(shù)，取w=2。

由于長(zhǎng)管拖車(chē)具有一定車(chē)速，所以需要將空間頻率功率譜密度Gq(n0)換算為時(shí)間頻率功率譜密度。時(shí)間頻率和空間頻率存在如下關(guān)系：

其中u為長(zhǎng)管拖車(chē)的速度，根據(jù)式（3）進(jìn)一步推導(dǎo)得

本次設(shè)定長(zhǎng)管拖車(chē)以80 km/h(22.3 m/s) 的速度在B級(jí)路面上行駛，利用Matlab通過(guò)諧波疊加法[14]得出其功率譜密度作為隨機(jī)振動(dòng)的激勵(lì)，如圖6所示。

圖6 B級(jí)路面功率譜密度

4 監(jiān)測(cè)參數(shù)與位置選取

對(duì)于監(jiān)測(cè)參數(shù)的選擇主要是通過(guò)隨機(jī)振動(dòng)實(shí)驗(yàn)，并監(jiān)測(cè)氣瓶由于損傷導(dǎo)致的各參數(shù)變化，選擇較為明顯的參數(shù)作為監(jiān)測(cè)參數(shù)。而監(jiān)測(cè)位置的選擇主要結(jié)合仿真結(jié)果、案例分析以及氣瓶爆破實(shí)驗(yàn)的結(jié)果進(jìn)行分析，得出較為合適監(jiān)測(cè)位置。

4.1 監(jiān)測(cè)參數(shù)選取

得出隨機(jī)振動(dòng)實(shí)驗(yàn)的激勵(lì)后，對(duì)兩端固定約束輸入上述振動(dòng)激勵(lì)，且X、Y、Z3個(gè)方向都加載B級(jí)路面的功率譜密度，對(duì)有損氣瓶損傷位置各參數(shù)進(jìn)行監(jiān)測(cè)，最終得出損傷位置的垂直加速度以及應(yīng)力變化較為明顯，且與損傷等級(jí)成正比，其變化如圖7和圖8所示。

圖7 應(yīng)力與損傷等級(jí)關(guān)系

圖8 垂直加速度與損傷等級(jí)關(guān)系

4.2 監(jiān)測(cè)位置選擇

4.2.1 仿真結(jié)果

隨機(jī)振動(dòng)實(shí)驗(yàn)的應(yīng)力云圖如圖9所示，從圖中可看出其最大應(yīng)力點(diǎn)出現(xiàn)在瓶肩處，而且由圖10可看出，X、Y、Z方向最大變形也為瓶肩位置，所以可得出瓶肩處為缺陷多發(fā)位置。

圖9 等效應(yīng)力云圖

圖10 等效變形云圖

4.2.2 案例統(tǒng)計(jì)

根據(jù)氣瓶的損傷類型，本文取較為常見(jiàn)的損傷進(jìn)行分析，即裂紋與腐蝕。

表4是中國(guó)特種設(shè)備檢測(cè)研究院對(duì)氣瓶檢驗(yàn)開(kāi)裂部位的統(tǒng)計(jì)，可以看出裂紋在氣瓶的各個(gè)部位均有可能發(fā)生，但以筒體外表面和瓶頸內(nèi)表面裂紋最為常見(jiàn)[15]。

表4 裂紋發(fā)生部位統(tǒng)計(jì)

除此之外，中國(guó)特檢院也收集了501個(gè)長(zhǎng)管拖車(chē)氣瓶腐蝕案例，腐蝕嚴(yán)重的區(qū)域絕大多數(shù)都位于氣瓶下部靠近端部的地方，其中以靠近后端部的居多，也有靠近前端部的，也有同時(shí)在靠近兩端發(fā)生腐蝕的。通過(guò)檢驗(yàn)發(fā)現(xiàn)氣瓶實(shí)測(cè)最小壁厚大多數(shù)位于端部瓶肩上或靠近瓶肩的部位，也有少數(shù)氣瓶的腐蝕區(qū)域位于氣瓶靠近中部某一位置處，少數(shù)氣瓶筒體整個(gè)內(nèi)表面都發(fā)生了腐蝕[15]。

但一般情況下，氣瓶裂紋的產(chǎn)生主要在制造中生成，而氣瓶運(yùn)輸中發(fā)生的減薄大多由于氣體腐蝕，所以將瓶肩作為損傷多發(fā)位置是合理的。

4.2.3 氣瓶爆破實(shí)驗(yàn)

按照ISO 11439—2013[16]要求，對(duì)該規(guī)格氣瓶進(jìn)行水壓實(shí)驗(yàn)，以尋找氣瓶較為合適的監(jiān)測(cè)位置，其升壓曲線如圖11所示，圖中有兩段15 min左右的保壓階段，氣瓶爆破前后如圖12所示。由圖12（a）可知此階段目的是為測(cè)試筒體與氣瓶瓶肩處的應(yīng)變大小，以選擇合適的監(jiān)測(cè)位置。從圖中可得出氣瓶屈服壓力為50 MPa左右，且其最終爆破壓力為59.45 MPa，其符合ISO11439—2013要求。

圖11 爆破實(shí)驗(yàn)壓力曲線

圖12 氣瓶爆破前后

氣瓶最終爆破位置如圖12（b）所示，其位置在氣瓶筒體，雖氣瓶最終爆破位置不處于瓶肩位置，這是由于氣瓶瓶肩厚度（平均厚度18 mm左右）相比于筒體（17 mm）較厚，對(duì)瓶肩有抗壓能力有一定的加強(qiáng)，由表5可看出，保壓階段瓶肩的應(yīng)變相比于筒體較大，可看出瓶肩應(yīng)為易損傷區(qū)域，因此將瓶肩作為氣瓶損傷監(jiān)測(cè)部分的監(jiān)測(cè)位置是合理的。

表5 保壓階段應(yīng)變

5 結(jié)束語(yǔ)

本文對(duì)長(zhǎng)管拖車(chē)氣瓶損傷在線監(jiān)測(cè)，提出較為合理的監(jiān)測(cè)參數(shù)以及監(jiān)測(cè)位置，分別對(duì)無(wú)損氣瓶以及損傷氣瓶進(jìn)行數(shù)值仿真分析，以隨機(jī)振動(dòng)實(shí)驗(yàn)?zāi)M氣瓶在運(yùn)輸中的狀況，監(jiān)測(cè)各參數(shù)的變化，得出當(dāng)出現(xiàn)損傷時(shí)其損傷位置的等效應(yīng)力與垂直加速度都有異常變化，且與損傷等級(jí)成正比的結(jié)論。而且其瓶肩與瓶口處的變形與應(yīng)力較大，由于瓶口存在氣閥，監(jiān)測(cè)難度較大，因此將瓶肩作為在線監(jiān)測(cè)平臺(tái)的監(jiān)測(cè)位置，為之后平臺(tái)搭建以及現(xiàn)場(chǎng)監(jiān)測(cè)實(shí)驗(yàn)做準(zhǔn)備。