孫曉燕， 鄭津洋， Dryver R.Huston， 秦 權(quán)， 王海龍

（1.浙江大學(xué)建筑工程學(xué)院 杭州，310058） （2.浙江大學(xué)材料化學(xué)學(xué)院 杭州，310027）（3.佛蒙特大學(xué)工程學(xué)院 佛蒙特州，美國 VT05405） （4.清華大學(xué)土木水利學(xué)院 北京，100084）

引 言

隨著石油等不可再生能源日益緊張，作為清潔可再生能源的氫燃料在能源領(lǐng)域中越來越得到重視，其生產(chǎn)、儲運(yùn)、使用都成為現(xiàn)階段的重要研究方向。其中，氫氣儲存技術(shù)是整個(gè)系統(tǒng)中最基本、最重要的技術(shù)之一［1］，高壓儲氫目前被普遍認(rèn)為是極具發(fā)展前景的儲氫技術(shù)。氫能汽車技術(shù)的開發(fā)實(shí)施使其有望代替燃油汽車成為未來交通工具，移動式的高壓氫氣儲存氣瓶在汽車行駛過程中的安全運(yùn)營也得到越來越多的關(guān)注?，F(xiàn)階段關(guān)于纖維纏繞壓力儲氫容器的結(jié)構(gòu)分析已經(jīng)較為系統(tǒng)［2-5］，相應(yīng)的生產(chǎn)、設(shè)計(jì)和安全評定準(zhǔn)則已形成行業(yè)技術(shù)規(guī)范［6-7］，其儲運(yùn)和使用過程中的安全和風(fēng)險(xiǎn)也被研究和關(guān)注［8-9］。但是，對于纖維纏繞壓力儲氫容器的模態(tài)分析和損傷識別研究還處在探討階段。以往傳統(tǒng)的安全評估方法建立在定期檢測的基礎(chǔ)上，不能滿足現(xiàn)有的工作環(huán)境改變和使用要求的提高帶來的安全行駛需要，這就給纖維纏繞壓力儲氫容器提出了在線監(jiān)測、實(shí)時(shí)識別損傷和失效預(yù)警的技術(shù)要求。

Nathanael等［10］針對碳纖維絲纏繞的罐體結(jié)構(gòu)開展了沖擊荷載作用下的模態(tài)試驗(yàn)，基于空間傳感器測試得到的三向正交加速度響應(yīng)確定頻響函數(shù)，并在此基礎(chǔ)上研究了通過非線性波調(diào)制譜進(jìn)行裂紋損傷識別的可行性。研究認(rèn)為即使考慮工作環(huán)境的噪聲干擾，高靈敏的多向量測傳感器也可以使得纖維纏繞罐體的裂紋識別成為可能。Sun等［11］針對高壓儲氫容器基于動態(tài)參數(shù)的損傷識別和泄漏預(yù)警技術(shù)進(jìn)行了理論研究，針對環(huán)向裂紋和徑向裂紋兩種疲勞裂紋和不同的裂紋貫穿深度，考慮了汽車行駛過程中的噪聲和路面不平度引起的振動干擾，以有限元仿真分析為基礎(chǔ)，探討了采用基于頻率和振型的損傷識別指標(biāo)進(jìn)行損傷識別的可行性。研究表明，頻率對儲氫容器的局部損傷并不敏感，導(dǎo)致基于頻率的識別指紋在損傷前后變化完全被環(huán)境噪聲和路面不平度振動所湮沒，而基于振型的損傷識別指標(biāo)可以實(shí)現(xiàn)裂紋損傷的識別和定位。

筆者基于0.074m3車用纖維纏繞儲氫容器開展了模態(tài)試驗(yàn)，針對反復(fù)充裝引起的容器桶身疲勞裂紋損傷狀態(tài)，探討了基于模態(tài)參數(shù)進(jìn)行損傷識別的可行性。該研究對車用氫氣儲存氣瓶的在線損傷識別研究奠定了基礎(chǔ)，對建立實(shí)時(shí)監(jiān)測及氫氣泄漏預(yù)警系統(tǒng)具有借鑒參考意義。

1 試驗(yàn)研究

車用復(fù)合纖維纏繞高壓儲氫容器基本結(jié)構(gòu)［6］主要由金屬內(nèi)襯、纏繞在內(nèi)襯外的增強(qiáng)復(fù)合材料層、包覆在增強(qiáng)復(fù)合材料層外表面的保護(hù)層、與容器一端封頭極孔連接的閥座等幾部分構(gòu)成，其結(jié)構(gòu)如圖1所示。

圖1 車用纖維纏繞高壓儲氫容器結(jié)構(gòu)構(gòu)造Fig.1 The structure detail of automobile hydrogen onboard storage tank

本次試驗(yàn)研究針對0.074m3的車用纖維纏繞高壓儲氫容器（材料性能見表1，2）開展，其結(jié)構(gòu)構(gòu)造從內(nèi)至外依次為鋁合金內(nèi)襯層 （厚度5mm）、碳纖維層（8.5mm）和玻璃纖維層（3.5mm）。其運(yùn)行使用中的約束狀態(tài)如圖2所示。該車用纖維纏繞高壓儲氫容器工作壓力為40MPa，質(zhì)量為38kg，容器桶身直徑為400mm，桶身長度為400mm。

模態(tài)試驗(yàn)的開展針對車用纖維纏繞高壓儲氫容器的完好狀態(tài)、輕微損傷和疲勞開裂3種狀態(tài)進(jìn)行。為減免在模擬氫氣反復(fù)充裝的疲勞試驗(yàn)過程中發(fā)生爆炸的風(fēng)險(xiǎn)，本次試驗(yàn)的充壓載體為水，疲勞試驗(yàn)在專用的防爆試驗(yàn)室進(jìn)行，如圖2所示。第1輪疲勞加載幅值為2MPa～35MPa，循環(huán)次數(shù)為15 000次，觀察儲氫容器的外觀及變形，施加第2輪疲勞荷載，加載幅值為2MPa～44MPa，循環(huán)次數(shù)為10 500次。此時(shí)按照設(shè)計(jì)計(jì)算分析儲氫容器存在輕微疲勞裂紋，但結(jié)構(gòu)仍可完好承受工作壓力，表明裂紋沒有貫穿，該狀態(tài)為模態(tài)試驗(yàn)中的輕微損傷狀態(tài)。然后將結(jié)構(gòu)反復(fù)疲勞加載直至裂紋貫穿，結(jié)構(gòu)破裂，此狀態(tài)為疲勞開裂狀態(tài)。為了模擬高壓儲氫容器在汽車運(yùn)行狀態(tài)的實(shí)際約束情況，如圖3所示，將容器采用兩根與地剛結(jié)的環(huán)形支架固定。在調(diào)試試驗(yàn)階段，根據(jù)振動信號與模態(tài)參數(shù)分析，將支架剛度調(diào)整并加高，以實(shí)現(xiàn)全罐體布置測點(diǎn)，并增加懸掛約束的試驗(yàn)工況進(jìn)行對比研究，如圖4所示。

表1 鋁合金內(nèi)襯與碳纖維材料參數(shù)Tab.1 Material parameters of aluminum alloy liner and carbon fiber

表2 碳纖維／環(huán)氧樹脂復(fù)合材料性能Tab.2 Properties of composite material carbon fiber／expoxyresin

圖2 疲勞加載中的儲氫容器Fig.2 Hydrogen on－board storage tank under working condition

圖3 車用纖維纏繞儲氫容器工作示意Fig.3 Hydrogen on－board storage tank under fatigue loading

圖4 模態(tài)試驗(yàn)的約束示意Fig.4 Constrain for model test

1.1 測試方法

錘擊法具有數(shù)據(jù)采集快、分布頻域較寬、所需儀器少、操作簡單等特點(diǎn)，比較適合于纖維纏繞儲氫罐體模型這樣的復(fù)雜空間結(jié)構(gòu)的模態(tài)分析，因此在試驗(yàn)中采用力錘激振、加速度傳感器采集響應(yīng)的方法。當(dāng)用力錘激振模型時(shí)，與力錘相連的力傳感器將力信號通過電荷放大器放大后傳輸?shù)絼討B(tài)信號分析系統(tǒng)，而響應(yīng)點(diǎn)處的加速度響應(yīng)則通過壓電陶瓷傳感器來測量，同樣通過電荷放大器傳回動態(tài)信號分析系統(tǒng)。經(jīng)過FFT分析，得到傳遞函數(shù)，然后用模態(tài)分析軟件對傳遞函數(shù)進(jìn)行分析，識別出結(jié)構(gòu)的模態(tài)參數(shù)。

本次試驗(yàn)采用多點(diǎn)激振一點(diǎn)測試的方法，以加速度響應(yīng)測點(diǎn)作為參考點(diǎn)，激振點(diǎn)為移動測點(diǎn)，即移動力錘進(jìn)行激振。為了提高測試精度，減少噪聲影響，在敲擊點(diǎn)上多次敲擊，經(jīng)過多段平均得到總體平均的頻率響應(yīng)函數(shù)。此外，在進(jìn)行信號譜分析時(shí)，對于力信號采用加力窗函數(shù)處理，響應(yīng)信號采用指數(shù)窗函數(shù)處理。由于單向激振法可能遺漏某階由兩個(gè)或兩個(gè)以上方向合成的振型，因而在本項(xiàng)目模型試驗(yàn)中對試驗(yàn)方法進(jìn)行了改進(jìn)，采用單點(diǎn)單向激振，然后集總平均求得模型的空間振動模態(tài)，避免了模態(tài)漏測并確保試驗(yàn)數(shù)據(jù)的完整可靠。

1.2 模態(tài)試驗(yàn)原理

設(shè)基座模型結(jié)構(gòu)在外力作用下時(shí)域內(nèi)的運(yùn)動方程為

其中：M，C，K分別為結(jié)構(gòu)的質(zhì)量、阻尼、剛度矩陣；x為位移列陣；f為外力列陣。

將式（1）兩邊進(jìn)行Laplace變換，考慮所有的振動模態(tài)并利用振動模態(tài)的正交性，可得到位移響應(yīng)與外力在頻域內(nèi)的關(guān)系式為

其中：Φr為第r階固有頻率下的振動模態(tài)；ω為激振力頻率；j為復(fù)數(shù)記號；Mr，Cr，Kr為第r階固有頻率下的模態(tài)質(zhì)量、阻尼、剛度。

式（2）表明，基座模型上各部分的振動情況X為各階振動模態(tài)的Mr，Cr，Kr，ω及Φr的組合，即模態(tài)疊加原理。

將（2）式表示為矩陣形式，即

其中：Hij是指在j點(diǎn)激振與i點(diǎn)響應(yīng)之間的位移－力的傳遞函數(shù)，可以表示為

由式（4）出發(fā)并拓寬Hij，可將傳遞函數(shù)矩陣的第j列的列陣表示為

由式（5）可知，只要測量傳遞函數(shù)矩陣中的第j列，就可識別出系統(tǒng)的各階模態(tài)參數(shù)。同理可推，只要測量傳遞函數(shù)矩陣中的第i行，也可識別出系統(tǒng)的各階模態(tài)參數(shù)，并不需要測出傳遞函數(shù)矩陣中的每一個(gè)元素。所以，對測試模型可以采用“固定激勵(lì)點(diǎn)移動測量點(diǎn)法”或“固定測量點(diǎn)移動激勵(lì)點(diǎn)法”來實(shí)現(xiàn)模態(tài)分析。在本基座的模型試驗(yàn)過程中，考慮到加速度傳感器安裝利用膠水粘結(jié)，移動不便，而敲擊點(diǎn)力錘移動則比較方便，所以采用“固定測量點(diǎn)移動激勵(lì)點(diǎn)法”來實(shí)現(xiàn)模態(tài)分析。

1.3 傳遞函數(shù)的測量

結(jié)構(gòu)試驗(yàn)的模態(tài)分析首先要測量傳遞函數(shù)，然后用模態(tài)分析軟件對傳遞函數(shù)進(jìn)行分析，識別出結(jié)構(gòu)的模態(tài)參數(shù)。傳遞函數(shù)（頻率響應(yīng)函數(shù)）是線性動力系統(tǒng)的固有特性，與激振函數(shù)的類型無關(guān)。整個(gè)試驗(yàn)過程中，在儀器選擇和試驗(yàn)技巧上遵循國家標(biāo)準(zhǔn)《機(jī)械導(dǎo)納的試驗(yàn)確定》——基本定義與傳感器（GB11349.1－2006）和沖擊激振法（GB11349.3－2006）。

1.4 模態(tài)參數(shù)分析

模態(tài)參數(shù)識別就是用數(shù)值方法得到與實(shí)測的傳遞函數(shù)曲線相吻合的擬合曲線，從而得到模態(tài)參數(shù)。這個(gè)數(shù)值逼近過程通常用最小二乘法來計(jì)算，以確定一條最優(yōu)的擬合曲線。分析儀器采用南京安正軟件工程有限責(zé)任公司開發(fā)的振動及動態(tài)信號采集分析系統(tǒng)CRAS V6.2，模態(tài)分析軟件采用 MaCras。CRAS振動及動態(tài)信號采集分析儀為16通道FFT分析儀，在精確性和通用性方面具有很高的標(biāo)準(zhǔn)。

1.5 試驗(yàn)系統(tǒng)和儀器設(shè)備

高壓儲氫容器罐身布置了270個(gè)測點(diǎn)，如圖5所示，其中1～30節(jié)點(diǎn)和240～270節(jié)點(diǎn)為約束節(jié)點(diǎn)。本次模態(tài)試驗(yàn)采用型號為KD1002的壓電加速度傳感器，其質(zhì)量為16g，電荷靈敏度為2.265pC／ms－2，諧振頻率為40kHz，使用頻率為1kHz～18kHz，最大可測值為2×104ms－2。采用型號為 CL－YD－302的力錘錘擊激勵(lì)，其參考靈敏度為3.82pC／N。上述試驗(yàn)儀器在測試前已進(jìn)行了檢查和校正，以保證試驗(yàn)結(jié)果的正確性。

圖5 高壓儲氫容器的罐身測點(diǎn)示意Fig.5 Measuring points on body of high－pressure hydrogen storage tank

2 試驗(yàn)結(jié)果分析

經(jīng)過嘗試性試驗(yàn)，發(fā)現(xiàn)懸掛約束由于激勵(lì)后達(dá)到穩(wěn)態(tài)時(shí)間與再次激勵(lì)時(shí)間交迭而導(dǎo)致交互干擾，致使動態(tài)響應(yīng)的測試精度降低，剛架支撐與懸掛約束的頻響函數(shù)如圖6所示。此外，高壓儲氫容器作為汽車燃料在行駛中固定在汽車后備箱處，采用的約束形式與試驗(yàn)中的支架約束相近，因此基于模態(tài)參數(shù)的損傷識別研究針對雙支架約束狀態(tài)開展進(jìn)行。

圖6 剛架支撐與懸掛約束的頻響函數(shù)對比Fig.6 Comparison of frequency response functions of tanks with rigid support asuspension constraints

2.1 完好狀態(tài)、輕微損傷狀態(tài)和疲勞破裂狀態(tài)的頻響函數(shù)對比

圖7 完好、輕微損傷、疲勞破裂狀態(tài)下儲氫容器的頻響函數(shù)對比Fig.7 Comparison of frequency response functions of intact，slightly damaged and fatigue cracked tanks

通過直接對比完好狀態(tài)、輕微損傷狀態(tài)和疲勞破裂狀態(tài)的頻響函數(shù)，如圖7所示，可以發(fā)現(xiàn)局部的疲勞裂紋對損傷前后的頻響函數(shù)有一定影響但并不敏感，對于局部的微細(xì)裂紋損傷則無法判定損傷是否發(fā)生并確定損傷位置。

2.2 完好狀態(tài)、輕微損傷狀態(tài)和疲勞破裂狀態(tài)的測試振型對比

由于結(jié)構(gòu)損傷影響其剛度，從而導(dǎo)致振動參數(shù)有所變化，因此可以通過直接對比結(jié)構(gòu)在損傷狀態(tài)的測試振型與完好狀態(tài)的各階振型［12］，以振型差異大的節(jié)點(diǎn)作為潛在損傷位置［13］。為避免氫氣爆炸風(fēng)險(xiǎn)，本次儲氫容器高壓反復(fù)充壓疲勞試驗(yàn)以水為加載體，試驗(yàn)至疲勞斷裂，記錄表面滲水點(diǎn)，如圖8所示。

圖8 高壓疲勞測試斷裂滲水點(diǎn)Fig.8 Seepage point during high－pressure fatigue testing

對比損傷前后測試信息，可以發(fā)現(xiàn)振型差異較大的節(jié)點(diǎn)為67，103，107，210，判斷它們?yōu)槲⒘鸭y可能發(fā)生位置，但也有一些差異節(jié)點(diǎn)不能排除是噪聲干擾和測試誤差產(chǎn)生的，如圖9所示。

圖9 儲氫容器輕微損傷前后前三階振型差比較Fig.9 The first，second and third order vibration mode differences of intact and slightly damaged tanks

對于疲勞破裂狀態(tài)，損傷節(jié)點(diǎn)振型變化突出，可以識別損傷位置為67，103和231，如圖10所示。試驗(yàn)測試得到的高壓儲氫罐體的疲勞測試斷裂損傷點(diǎn)在104附近。振型直接比較方法是基于某單一振型信息，誤差比較大。采用該方法直接比較識別損傷還缺乏客觀的判斷依據(jù)，對于其他存在振型差異的節(jié)點(diǎn)難以判斷是否存在損傷。由于缺乏對多層高壓儲氫器罐體微損傷的測試信息，仍然難以明確其損傷部位。

圖10 儲氫容器疲勞破裂前后前三階振型差比較Fig.10 The first，second and third order vibration mode differences of intact and fatigue cracked tanks

3 結(jié)束語

筆者針對0.074m3車用纖維纏繞儲氫容器開展了模態(tài)試驗(yàn)，由于車用儲氫容器在使用過程中反復(fù)充裝，疲勞裂紋成為其常見失效模式?；诓煌瑩p傷狀態(tài)前后頻響函數(shù)、振型進(jìn)行損傷識別和定位的研究，結(jié)果表明，對于疲勞裂紋這類局部損傷，基于振型的識別方法可以實(shí)現(xiàn)損傷判定并預(yù)測損傷位置。盡管直接對比損傷前后的各階振型差可以識別損傷，但是由于該方法缺乏客觀判據(jù)，并且受測試誤差和環(huán)境噪聲干擾，容易產(chǎn)生誤判?，F(xiàn)階段多向傳感器和高精度測試技術(shù)已經(jīng)成熟，車用纖維纏繞儲氫容器的在線監(jiān)測、實(shí)時(shí)識別損傷和氫氣泄露預(yù)警的硬件基礎(chǔ)已經(jīng)具備，更多的研究和分析應(yīng)投入到對于纖維復(fù)合結(jié)構(gòu)的局部損傷更為敏感的損傷識別指標(biāo)中去。

［1］ Andrew H W，Salvador M A，F(xiàn)rancisco E L，et al.Delivery of cold hydrogen in glass fiber composite［J］.International Journal of hydrogen.2009，34：9773－9780.

［2］ Zu Lei，Zhang Dinghua，Xu Yingqiang，et al.Integral design and simulation of composite toroidal hydrogen storage tanks［J］.International Journal of Hydrogen，2012，37：1027－1036.

［3］ Liu Pengfei，Xu Ping，Zheng Jinyang.Artificial immune system for optimal design of composite hydrogen storage vessel［J］.Computational Materials Science，2009，47：261－267.

［4］ 魏喜龍，孫銀寶，陳日東，等.纖維纏繞壓力容器幾何非線性分析［J］.纖維復(fù)合材料，2007（1）：24－27.

Wei Xilong，Sun Yinbao，Chen Yuedong，et al.Geometrical nonlinear analysis of filament－winding vessel［J］.Fiber Composites，2007（1）：24－27.（in Chinese）

［5］ Xu Ping，Zheng Jinyang，Liu Pengfei.Finite element analysis of burst pressure of composite hydrogen storage vessels［J］.Materials and Design，2009，30：2295－2301.

［6］ 中國國家標(biāo)準(zhǔn)化管理委員會.GB／T 6058—2005纖維纏繞壓力容器制備和內(nèi)壓試驗(yàn)方法［S］.2005.

［7］ 中國國家標(biāo)準(zhǔn)化管理委員會.GB／T19624—2004在用含缺陷壓力容器安全評定［S］.2004.

［8］ Camara S，Bunsell A R，Thionnet A，et al.Determination of lifetime probabilities of carbon fiber composite plates and pressure vessels for hydrogen storage［J］.International Journal of Hydrogen 2011，36：6031－6038.

［9］ 鄭津洋，開方明，劉仲強(qiáng)，等.高壓氫氣儲運(yùn)設(shè)備及其風(fēng)險(xiǎn)評價(jià)［J］.太陽能學(xué)報(bào)，2006，27（11）：1168－1174.

Zheng Jinyang，Kai Fangming，Liu Zhongqiang，et al.Risk assessment and control of high pressure hydrogen equipment［J］.Acta Energiae Solaris Sinica，2006，27（11）：1168－1174.（in Chinese）

［10］Nathanael C Y，Douglas E A，Matt T.Muiti－dimensional sensing for impact load and damage evaluation in a carbon filament wound canister［J］.Material Evaluation，2008（7）：756－762.

［11］Sun Xiaoyan，Dryver R H，Qin Quan，et al.Damage detection and leakage alert of fiber composite wrapped tank for high pressure hydrogen storage［C］∥SPIE Smart Structures／Non－destructive Testing.San Diego：Published by SPIE，2008.

［12］劉洋，段忠東，周道成.基于模態(tài)綜合技術(shù)的結(jié)構(gòu)有限元模型修正［J］.振動、測試與診斷，2009，29（3）：287－291.

Liu Yang，Duan Zhongdong，Zhou Daocheng.Finite element model updating of complex structure using component mode synthesis technique［J］.Journal of Vibration，Measurement ＆ Diagnosis，2009，29（3）：287－291.（in Chinese）

［13］黃維平，王曉燕.基于改進(jìn)CHC遺傳算法的結(jié)構(gòu)損傷診斷研究［J］.振動、測試與診斷，2007，27（3）：232－235.

Huang Weiping，Wang Xiaoyan.Structural damage diagnosis based on improved CHC algorithm［J］.Journal of Vibration，Measurement ＆ Diagnosis，2007，27（3）：232－235.（in Chinese）