燃料電池城市客車儲氫瓶疲勞分析

唐學(xué)志 趙雨東 樊 彬 王文軍

（清華大學(xué) 汽車安全與節(jié)能國家重點實驗室）

1 前言

燃料電池城市客車作為新能源汽車的代表之一，其安全性備受關(guān)注。燃料電池客車的安全包括氫安全、電安全及碰撞安全等。在3項安全問題中，氫安全最被人們所關(guān)注。氫安全是指氫的特性引起的安全問題，主要涉及高壓儲氫瓶、供氫管路和燃料電池。國外關(guān)于燃料電池客車的氫安全問題的研究主要集中在氫泄漏、著火安全性等方面，而關(guān)于儲氫瓶疲勞方面的研究未見相關(guān)報道；我國在這方面的研究也處于起步階段。

本文所研究的儲氫瓶是我國某企業(yè)生產(chǎn)的鋁合金內(nèi)膽碳纖維全纏繞復(fù)合瓶，瓶口通過螺紋連接安裝了某公司開發(fā)的質(zhì)量為6.71 kg的集成閥，集成閥和儲氫瓶隨車振動時可能會引起儲氫瓶瓶口附近疲勞破壞，為此，通過試驗測量并結(jié)合相關(guān)理論對儲氫瓶瓶口附近的疲勞強(qiáng)度進(jìn)行分析研究。

2 分析流程

根據(jù)文獻(xiàn)[1]～文獻(xiàn)[7]確定儲氫瓶疲勞分析流程為:

a.通過試驗?zāi)B(tài)分析測量裝有集成閥儲氫瓶的前幾階固有頻率及其振型，進(jìn)而找出儲氫瓶的關(guān)鍵頻率和振型；

b.測量燃料電池城市客車實際運(yùn)行時車頂?shù)恼駝忧闆r，得到車頂振動的功率譜密度（PSD）和時域數(shù)據(jù)；

c.根據(jù)實測PSD譜和時域數(shù)據(jù)，簡化得到振動試驗臺所需的評價譜 （PSD譜和正弦定頻譜），將儲氫瓶安裝在振動臺上，按照各評價譜分別激振，同時測量瓶口的動態(tài)應(yīng)變；

d.根據(jù)所測動態(tài)應(yīng)變計算出等效應(yīng)力幅及等效平均應(yīng)力，得到多軸疲勞理論下的Goodman圖，分析儲氫瓶的疲勞特性。

3 儲氫瓶試驗?zāi)B(tài)分析

3.1 主要設(shè)置

理論上試驗?zāi)B(tài)分析是用力錘敲擊物體得到輸入信號 f（t）和物體上的加速度響應(yīng)信號x（t），然后進(jìn)行頻譜變換得到頻響函數(shù) H（f）=X（f）/F（f），再根據(jù)模態(tài)參數(shù)估計算法得到物體的固有頻率、模態(tài)阻尼及模態(tài)振型。

儲氫瓶的試驗?zāi)B(tài)分析方法為，利用細(xì)鋼絲繩懸吊儲氫瓶（圖1），采用PCB公司加速度傳感器和具有硬塑料錘頭的力錘。傳感器安裝位置見圖2。儲氫瓶模態(tài)試驗的測點為64個，瓶口模態(tài)試驗的測點為39個，每個測點敲擊3次，模態(tài)參數(shù)估計算法選用系統(tǒng)實現(xiàn)算法（ERA）。

3.2 試驗?zāi)B(tài)分析結(jié)果

在儲氫瓶裝有集成閥情況下，試驗?zāi)B(tài)分析得出的瓶體前3階固有頻率和第1階振型分別如表1和圖3所示。由表1和圖3可看出，儲氫瓶的第1階固有頻率遠(yuǎn)高于正常車頂?shù)姆治鲱l率（0～75 Hz），且振動主要發(fā)生在儲氫瓶中部。

由于儲氫瓶瓶口附近的疲勞特性更關(guān)鍵，所以對包括集成閥在內(nèi)的儲氫瓶前端瓶口進(jìn)行了瓶口試驗?zāi)B(tài)分析，結(jié)果如表2和圖4所示。由表2和圖4可看出，第1階振動的最大振幅發(fā)生在集成閥的質(zhì)心附近，第1階固有頻率也遠(yuǎn)高于正常的車頂分析頻率。

表1 儲氫瓶前3階固有頻率 Hz

表2 瓶口前3階固有頻率 Hz

由上述試驗結(jié)果可知，在正常的運(yùn)行過程中，燃料電池城市客車的儲氫瓶瓶體和集成閥的第1階共振是不會發(fā)生的，儲氫瓶的疲勞應(yīng)主要受車頂振動的影響。

4 車頂振動測量

為能在振動試驗臺上重現(xiàn)燃料電池城市客車儲氫瓶口附近振動情況，以便進(jìn)行疲勞評價，對實際運(yùn)行的燃料電池城市客車儲氫瓶口附近的車頂隨機(jī)振動進(jìn)行了測量。對所測數(shù)據(jù)進(jìn)行頻域和時域處理后，得出了振動試驗臺重現(xiàn)實際振動所需要的PSD，以及能在振動試驗臺上產(chǎn)生比實際振動強(qiáng)度大的加速PSD和最大峰值正弦譜。

測試所用燃料電池城市客車作為公共汽車每天運(yùn)行于北京頤和園北宮門與人民大學(xué)東門之間，行程為15 km，且每日需從頤和園北宮門返回加氫站，行程為10 km。

如圖5所示，在1個儲氫瓶靠近瓶口的前支架附近的車頂上布置了1個加速度傳感器。共測量3組數(shù)據(jù)，第1組和第2組數(shù)據(jù)為客車在往返全程載客路線行駛時測得，第3組數(shù)據(jù)為客車返回加氫站過程中測得。儲氫瓶前支架附近車頂振動測量結(jié)果（第1組數(shù)據(jù)）見圖6。

由圖6及其它2組數(shù)據(jù)（本文略）可知，車頂?shù)恼駝又饕性?0 Hz內(nèi)。考慮到振動試驗臺難以實現(xiàn)頻率低于2 Hz的振動，以及功率譜密度中應(yīng)盡可能反映實際振動能量所在的頻率范圍，所以對2～200 Hz頻率范圍內(nèi)的功率譜密度進(jìn)行分析，得到的可用于振動試驗臺激振并重現(xiàn)車頂振動的功率譜密度見圖7。在2～200 Hz頻率范圍內(nèi)，選取功率譜密度恒為6.299×10-4g2/Hz，可得到加速功率譜密度譜。

根據(jù)測量數(shù)據(jù)及最大峰值頻率分析理論，可得出3組測量時域上的最大加速度值及其對應(yīng)的頻率，即最大峰值正弦評價譜，見表3。

表3 最大峰值正弦評價譜

5 振動重現(xiàn)及動態(tài)應(yīng)變測量

用上述測得的PSD、加速PSD和最大峰值正弦譜，在振動試驗臺上對儲氫瓶激振，使儲氫瓶瓶口附近重現(xiàn)實車運(yùn)行條件下的振動情況，或產(chǎn)生比實際運(yùn)行情況更劇烈的振動，采用應(yīng)變法測量儲氫瓶瓶口附近的動態(tài)應(yīng)力。

因儲氫瓶較長，所以在振動臺上安置儲氫瓶時，將激振基座支撐在儲氫瓶前端靠近瓶口的位置上，而瓶的后端采用2個帶軸承的小輪（簡稱車輪）實現(xiàn)鉸鏈支撐（圖8）。應(yīng)變片A、B分別貼在儲氫瓶瓶口的上、下端。

6 疲勞分析

基于多軸疲勞理論和Goodman疲勞準(zhǔn)則對儲氫瓶的疲勞情況進(jìn)行分析。

首先根據(jù)應(yīng)變片在 3 個軸（0°、45°、90°）上的應(yīng)變得到主應(yīng)力，然后計算出最大交變應(yīng)力幅和最大平均應(yīng)力，再根據(jù)多軸疲勞狀態(tài)下的等效應(yīng)力幅和等效平均應(yīng)力計算方法[7]，計算出不同評價譜下激振瓶口振動狀況的等效應(yīng)力幅Sqa和等效平均應(yīng)力Sqm，見表 4。

表4 各疲勞載荷譜振動下的Sqa和Sqm值 MPa

根據(jù)儲氫瓶瓶口材料、結(jié)構(gòu)尺寸和載荷情況以及文獻(xiàn) [8]，得出儲氫瓶構(gòu)件的等效屈服極限Sy=21.88 MPa，等效疲勞極限Sf=9.95 MPa，等效拉伸極限Su=24.88 MPa。

根據(jù)不同評價譜下激振引起儲氫瓶瓶口振動的等效應(yīng)力幅和等效平均應(yīng)力，以及儲氫瓶構(gòu)件的各疲勞極限，可作出儲氫瓶在各評價譜振動情況下的Goodman圖。圖9為功率譜密度評價譜（圖7）振動下的Goodman圖。

從圖9及表4中其它各評價譜振動的Good?man圖（本文略）可知，在各疲勞載荷譜下，瓶口上端測點A和下端測點B都在Goodman疲勞分析圖中強(qiáng)度可接受的區(qū)域內(nèi)，說明A、B測點及瓶口內(nèi)螺紋危險點都是安全的。但A點的等效應(yīng)力幅和等效平均應(yīng)力都大于B點，表明A點的疲勞強(qiáng)度比B點略低。

加速評價PSD和最大峰值正弦譜振動下（本文略），等效應(yīng)力幅和等效平均應(yīng)力變化幅度都在1～4 MPa的范圍內(nèi)，說明儲氫瓶安裝集成閥后，在瓶口疲勞中起決定因素的不是實際行駛過程中由集成閥振動引起的應(yīng)力變化，而可能是由儲氫瓶充放氣過程和充放氣次數(shù)導(dǎo)致的。

7 結(jié)束語

為分析某燃料電池城市客車儲氫瓶因懸置安裝集成閥引起的瓶口附近疲勞情況，通過測量儲氫瓶模態(tài)和實車運(yùn)行時的車頂加速度，得出在振動試驗臺重現(xiàn)振動所需的疲勞評價譜。利用評價譜在振動試驗臺上重現(xiàn)了儲氫瓶振動，測量了儲氫瓶瓶口在各種評價譜下的動態(tài)應(yīng)力，根據(jù)多軸疲勞理論分析了儲氫瓶瓶口附近的疲勞特性。分析結(jié)果表明，儲氫瓶瓶口的疲勞強(qiáng)度在可接受的區(qū)域內(nèi)，集成閥的振動不是儲氫瓶瓶口疲勞的決定因素。

1 Uday Senapati， Yogesh Dhage， Vinaya Sawant， et al.Ac?celerated Test Method for Validation of Vehicle Components Subjected to Fatigue Failure.SAE 2007-26-072.

2 Yi Zhang，Mohammad Usman.Life Prediction for Lighting Bulb Shield Designs Subjected to Random Vibration.SAE 1999-01-0705.

3 Steve Bogert.Integration of Fatigue Technology with Modal Analysis.SAE 933068.

4 Hong Su.Vibration Test Specification for Automotive Prod?ucts Based on Measured Vehicle Load Data.SAE 2006-01-0729.

5 Jack A Collins.Mechanical Design of Machine Elements and Machines.New York :John Wiley&Sons， Inc.2003.

6 H.O.福克斯，R.I.斯蒂芬斯，著，漆平生，徐桂琴，譯.工程中的金屬疲勞.北京:中國農(nóng)業(yè)機(jī)械出版社，1983.

7 Ralph I Stephens， Ali Fatemi， Robert R Stephens et al.Metal Fatigue in Engineering （second edition）.New York:John Wiley&Sons，Inc.2000.

8 John H Bickford.An Introduction to the Design and Behav?ior of Bolted Joints （Second Edition）.New York:Marcel Dekker，Inc.1990.