楊剛 歐晨希 郭嘉煒 黃思

摘要：為研究充裝完成后的碳纖維全纏繞儲(chǔ)氫氣瓶內(nèi)熱力學(xué)參數(shù)隨時(shí)間的變化規(guī)律，以氫氣為工質(zhì)，選取了3種常用規(guī)格的車用儲(chǔ)氫氣瓶（27L、52L和135L）進(jìn)行傳熱計(jì)算。對(duì)于氣瓶內(nèi)的氫氣采用飽和均質(zhì)模型;對(duì)于氣瓶中的各隔層采用多層穩(wěn)態(tài)熱傳導(dǎo)的結(jié)構(gòu)模型。根據(jù)飽和狀態(tài)下的熱平衡條件以及熱力學(xué)關(guān)系，計(jì)算得到了氣瓶內(nèi)氫氣溫度、壓力以及氣瓶碳纖維纏繞層溫度隨時(shí)間的變化規(guī)律。計(jì)算結(jié)果表明：氣瓶容積越大，充裝完成時(shí)刻的氫氣溫度以及碳纖維纏繞層溫度越高;靜置過程中氫氣溫度以及碳纖維纏繞層溫度均呈同趨勢(shì)緩慢下降，且氣瓶容積越大，下降的幅度越小，碳纖維纏繞層溫度變化小于氫氣溫度變化。氫氣壓力在充裝完成時(shí)刻均為70 MPa，氣瓶容積越大，靜置過程中氫氣壓力下降速度越小。

關(guān)鍵詞：碳纖維全纏繞儲(chǔ)氫氣瓶;熱力學(xué)參數(shù)變化;飽和均值模型;熱傳導(dǎo)

中圖分類號(hào)：TH49

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1009-9492f 2022） 02-0007-04

0 引言

氫能是21世紀(jì)最有潛力的新型能源，正在被廣泛應(yīng)用。其中氫燃料電池汽車是突出的代表，氫能安全的話題自然成為大眾關(guān)注的重點(diǎn)[1]。目前，車用高壓氣態(tài)儲(chǔ)氫容器多采用碳纖維增強(qiáng)復(fù)合材料，作為儲(chǔ)氫氣瓶的承壓材料，具有重量輕、使用壽命長、強(qiáng)度高等優(yōu)點(diǎn)，從而保證車用儲(chǔ)氫系統(tǒng)具有足夠的質(zhì)量儲(chǔ)氫密度[2]。但是由于碳纖維增強(qiáng)復(fù)合材料對(duì)溫度的敏感性較高，若氣瓶的使用溫度過高，時(shí)間久了會(huì)導(dǎo)致金屬內(nèi)膽與纖維纏繞層形成脫粘[3]，甚至引起瓶體宏觀變形、碳纖維失效，給氣瓶帶來很大的安全隱患[4]。

近年來，國內(nèi)外學(xué)者相繼開展了碳纖維全纏繞儲(chǔ)氫氣瓶內(nèi)部熱力學(xué)參數(shù)變化的研究。伊藤裕一等[5]開展了氫氣瓶快充溫升的數(shù)值模擬研究，采用3D幾何模型，應(yīng)用理想氣體狀態(tài)方程對(duì)高壓氫氣進(jìn)行求解;EtienneWer-len等[6]對(duì)氫氣快充過程進(jìn)行了數(shù)值計(jì)算并推導(dǎo)出充裝完成時(shí)氣瓶內(nèi)部溫度與壓力等的解析解;Dicken等[7]對(duì)快充溫升進(jìn)行了數(shù)值分析，應(yīng)用R-K真實(shí)氣體狀態(tài)方程得出了用于絕熱狀態(tài)下快充過程氣體溫升的計(jì)算公式;劉延雷和劉格思[8]考慮氣體湍流、真實(shí)氣體效應(yīng)以及壁面?zhèn)鳠幔ㄟ^Fluent軟件建立模型對(duì)氣瓶不同充裝條件下的溫升狀況進(jìn)行預(yù)測(cè);王光緒等[9]通過Fluent軟件對(duì)不同長徑比與進(jìn)口直徑的氣瓶快速充氣溫升過程進(jìn)行研究。

上述研究主要集中在碳纖維儲(chǔ)氫氣瓶快充過程中的溫升研究，缺少氣瓶充裝后熱力學(xué)參數(shù)隨時(shí)間規(guī)律變化的研究。因此，在本文的研究中，選取Ⅲ型儲(chǔ)氫氣瓶為研究對(duì)象，對(duì)充裝完成后的碳纖維全纏繞儲(chǔ)氫氣瓶進(jìn)行傳熱學(xué)計(jì)算分析，得到氣瓶內(nèi)溫度、壓力等熱力學(xué)參數(shù)隨時(shí)間的變化規(guī)律，為車載LNC氣瓶的安全問題提供技術(shù)支持。

1 計(jì)算模型和計(jì)算方法

1.1 物理模型

圖1所示為所研究的車載儲(chǔ)氫氣瓶的計(jì)算模型。該儲(chǔ)氫氣瓶為多層結(jié)構(gòu)，內(nèi)膽為鋁合金，起密封防漏作用;包裹鋁合金內(nèi)襯層的為碳纖維纏繞層，起承壓結(jié)構(gòu)作用;外層為玻璃纖維層，起隔熱作用。氣瓶充裝結(jié)束，處于放置過程中，充氣口已封口。

1.2 計(jì)算方法

本研究采用文獻(xiàn)[10]提供的計(jì)算方法，該方法已得到相關(guān)實(shí)驗(yàn)數(shù)據(jù)的驗(yàn)證[11]。

由于所研究的碳纖維儲(chǔ)氫氣瓶容積較小，介質(zhì)較快達(dá)到熱平衡，不容易形成溫度分層，因此采用飽和均質(zhì)模型對(duì)氣瓶進(jìn)行傳熱計(jì)算。根據(jù)熱力學(xué)定律在At時(shí)間內(nèi)有如下關(guān)系[12]：

Q=C.m.△T=φ.△t

（1）式中：Q為氣瓶釋放的熱量;C為氫氣飽和狀態(tài)比熱容;m為氫氣質(zhì)量;△T為△t時(shí)間間隔前后氣瓶內(nèi)溫差;φ為漏熱率。

理想狀態(tài)下氫氣的飽和狀態(tài)比熱容計(jì)算公式[13]：

氣體質(zhì)量根據(jù)理想氣體狀態(tài)方程pV=m/M RT以及初始?jí)毫Α囟扔?jì)算可得，以鋁合金作內(nèi)襯的儲(chǔ)氫氣瓶對(duì)氣體的封閉性較好，忽略氣瓶泄漏量，認(rèn)為氣瓶在重裝完成后瓶內(nèi)氣體的質(zhì)量保持不變。且根據(jù)理想氣體狀態(tài)方程，可得氣瓶內(nèi)部在i時(shí)刻的壓力為：

式中：R為氣體常數(shù)，R=8.314 kJ/ （kg.K）;M為氫氣物質(zhì)的摩爾質(zhì)量。

為計(jì)算氣瓶的漏熱率φ作出如下假設(shè)：（1）漏熱量只考慮以熱傳導(dǎo)的方式;（2）整個(gè)氣瓶的傳熱方向?yàn)楦鞲魧拥姆ň€方向，氣瓶總熱阻是多個(gè)薄層熱阻的疊加;（3）鋁合金襯里、碳纖維/環(huán)氧復(fù)合層壓板（CFEC）和玻璃纖維/環(huán)氧復(fù)合層壓板（GFEC）的傳熱特性被認(rèn)為是各向同性的;（4）充裝完成時(shí)刻氣瓶內(nèi)壓力達(dá)到了標(biāo)稱壓力;（5）忽略瓶口對(duì)傳熱的影響;（6）認(rèn)為氣瓶充裝完成后瓶內(nèi)較快達(dá)到熱平衡狀態(tài)。

該氣瓶漏熱主要考慮網(wǎng)簡體φ1和封頭φ2[14]：

給定氣瓶內(nèi)部溫度為T，碳纖維纏繞層內(nèi)徑處溫度為Te，即以Te表征碳纖維纏繞層的溫度。根據(jù)多層網(wǎng)筒壁面的穩(wěn)態(tài)熱傳導(dǎo)計(jì)算公式可以計(jì)算得碳纖維纏繞層溫度，由式（6）得到i時(shí)刻氣瓶得漏熱率φi，且滿足關(guān)系式：

1.3 具體算例

本文選取了3個(gè)具有代表性的儲(chǔ)氫氣瓶容積27 L、52 L和135 L來研究碳纖維全纏繞儲(chǔ)氫氣瓶內(nèi)熱力學(xué)參數(shù)和時(shí)間的變化規(guī)律。各容積氣瓶初始條件如表1所示。

已知環(huán)境溫度Ts=20℃，標(biāo)稱工作壓力po=70 MPa，初始溫度P參考有關(guān)文獻(xiàn)可得[15]，給定to=0，te=600 s，△T=0.1℃，即計(jì)算從氣瓶充裝完成時(shí)刻開始至靜置10 min后結(jié)束。氣瓶物性參數(shù)如表2所示。

計(jì)算流程如圖2所示，根據(jù)圖2在Matlab上編制程序進(jìn)行運(yùn)算，由所得數(shù)據(jù)可繪制出氫氣溫度T與時(shí)間t的關(guān)系曲線圖、氫氣壓力P與時(shí)間t的關(guān)系曲線圖以及氣瓶碳纖維纏繞層溫度Te與時(shí)間的關(guān)系曲線圖。

2 計(jì)算結(jié)果分析

由于碳纖維增強(qiáng)復(fù)合材料對(duì)溫度變化的敏感性較高，根據(jù)GB/T35544-2017規(guī)定，氣瓶的使用溫度不宜過高，過高的溫度會(huì)誘發(fā)環(huán)氧樹脂剝離、碳纖維失效，故研究儲(chǔ)氫氣瓶充裝完成后的熱力學(xué)參數(shù)隨時(shí)間變化規(guī)律是十分有必要的。

2.1 瓶內(nèi)氫氣溫度T和時(shí)間t的關(guān)系

氫氣充裝完成后，瓶內(nèi)氫氣溫度高于環(huán)境溫度，使得氣瓶與外界發(fā)生熱傳遞，瓶內(nèi)氫氣溫度不斷下降。圖3所示為3種容積的碳纖維全纏繞儲(chǔ)氫氣瓶內(nèi)氫氣溫度T隨時(shí)間t的變化曲線圖。由圖可知，氣瓶容積越大，充裝完成時(shí)刻的氫氣溫度越高，這是由于容積增大后，氣瓶換熱表面積與容積的比值減小，氣瓶蓄熱能力增強(qiáng)。在隨后的靜置過程中，氫氣溫度呈緩慢下降的趨勢(shì)，且氣瓶容積越大，氫氣溫度下降的幅度越小。從氣瓶充裝完成時(shí)刻到靜置10 min的過程中，27 L氣瓶內(nèi)部溫度從74.2℃下降到了35.7℃;52 L氣瓶內(nèi)部溫度從79.9℃下降到了42.2℃;135 L氣瓶內(nèi)部溫度從83.7℃下降到了59.4℃，瓶內(nèi)氣體的平均溫度滿足安全要求[16]。

2.2 瓶內(nèi)氫氣壓力p和時(shí)間t的關(guān)系

在充裝氫氣結(jié)束時(shí)瓶內(nèi)氫氣壓力達(dá)到70 MPa，根據(jù)理想氣體狀態(tài)方程，當(dāng)氫氣溫度降低時(shí)，其壓力也會(huì)隨之降低。圖4所示為3種容積的碳纖維全纏繞儲(chǔ)氫氣瓶內(nèi)氫氣壓力P隨時(shí)間t的變化曲線圖。由圖可知，氣瓶容積越小壓力下降越快。從氣瓶充裝完成時(shí)刻到靜置10 min的過程中，27 L氣瓶內(nèi)壓力下降到了62.4 MPa，52 L氣瓶內(nèi)壓力下降到了62.7 MPa，135 L氣瓶內(nèi)壓力下降到了65.3 MPa。由此可看出儲(chǔ)氫氣瓶內(nèi)部壓力變化與氣瓶容積有較大關(guān)系。

2.3 碳纖維纏繞層溫度Te和時(shí)間t的關(guān)系

碳纖維纏繞層是以環(huán)氧樹脂作為基體，碳纖維以全纏繞的方式附在鋁合金層外表，碳纖維復(fù)合層具有很好的拉伸性能，斷裂韌性高，可有效傳遞載荷[17]，是高壓氣瓶的主要承壓結(jié)構(gòu)。但是碳纖維對(duì)溫度變化有較高敏感性，高溫條件下環(huán)氧樹脂會(huì)從樹脂粘合劑上剝離，環(huán)氧樹脂的機(jī)械性能下降[18]，造成碳纖維復(fù)合材料的失效，進(jìn)而影響到車用儲(chǔ)氫系統(tǒng)的安全[19]。

由于所研究儲(chǔ)氫氣瓶內(nèi)膽材料為鋁合金，導(dǎo)熱系數(shù)較大，隔熱性能差，碳纖維纏繞層的溫度也會(huì)隨著氣瓶內(nèi)部溫度變化。圖5所示為3種容積的碳纖維全纏繞儲(chǔ)氫氣瓶碳纖維纏繞層溫度Te隨時(shí)間t的變化曲線圖。由圖可知，氣瓶容積越大，充裝完成時(shí)刻的氣瓶碳纖維纏繞層溫度越高。在隨后的靜置過程中，氣瓶碳纖維纏繞層溫度都呈緩慢下降的趨勢(shì)，氣瓶容積越大，溫度下降的幅度越小，且碳纖維纏繞層溫度變化小于氣瓶內(nèi)部溫度變化。從氣瓶充裝完成時(shí)刻到靜止10 min的過程中，27 L氣瓶碳纖維纏繞層最高溫度從74.1℃下降到了36.3℃;52 L氣瓶碳纖維纏繞層最高溫度從79.7℃下降到了42.9℃：135 L氣瓶碳纖維纏繞層最高溫度從83.4℃下降到了59.9℃。

3 結(jié)束語

本文以氫氣為工質(zhì)，選取27 L、52 L、135 L三種規(guī)格的鋁內(nèi)膽碳纖維全纏繞儲(chǔ)氫氣瓶進(jìn)行傳熱計(jì)算，采用飽和均質(zhì)模型得到了如下結(jié)論。

（1）氣瓶容積越大，充裝完成時(shí)刻的氫氣溫度越高。氣瓶充裝完成后，氣瓶內(nèi)部溫度高于環(huán)境溫度，使得氣瓶與外界發(fā)生熱傳遞，瓶內(nèi)溫度下降，在結(jié)束充裝的10 min內(nèi)呈緩慢下降的趨勢(shì)，且氣瓶容積越大，氫氣溫度下降的幅度越小，溫度在氣瓶使用安全范圍之內(nèi)。

（2）瓶內(nèi)氫氣溫度下降，其壓力也隨之下降，氣瓶容積越大，氫氣壓力下降速度越小，且氣瓶內(nèi)壓力變化與氣瓶容積大小有較大關(guān)系。

（3）由于鋁合金隔熱性能較差，碳纖維纏繞層溫度隨氫氣溫度變化，在結(jié)束充裝的10 min內(nèi)也呈緩慢下降趨勢(shì)，氣瓶容積越大，碳纖維層溫度下降的幅度越小，且碳纖維纏繞層溫度變化小于瓶內(nèi)溫度變化。

參考文獻(xiàn)：

[1]黃正夏.深入認(rèn)識(shí)開發(fā)氫能源的戰(zhàn)略意義[J].時(shí)代汽車，2009（9）：44-46.

[2]秦玉琪，袁奕雯，楊振國.纖維纏繞儲(chǔ)氫氣瓶及燃料汽車應(yīng)用現(xiàn)狀綜述[J].中國特種設(shè)備安全，2019，35（2）：70-75.

[3]劉延雷，鄭津洋，韋新華，等，復(fù)合材料氫氣瓶快充過程溫升控制方法研究[J].太陽能學(xué)報(bào)，2012，33（9）：1621-1627.

[4] Mengxiao Li， Yunfeng Bai， Caizhi Zhang， et al.Review on the re-search of hydrogen storage system fast refueling in fuel cell vehi-cle[J]. International Journal of Hydrogen Energy， 2019. 44（21）：10677-10693，

[5]伊藤裕一.自動(dòng)車用壓縮水素容器一④急速充填④數(shù)值解析[J].自動(dòng)車研究，2006（28）： 29-32.

[6] Werlen，E.Thermal effects related to H2 fast filling in high pres-surevessels depending on vessels type and filling procedures：modeling， trials and studies[C]//European Energy Conference.September 2003.

[7]C J B Dicken，W Merida. Measured effects of filling time and ini-tial mass on the temperature distribution within a hydrogen cylin-der during refuelling[J]. Journal of Power Sources， 2006， 165（1）：324-336.

[8]劉延雷，鄭津洋，韓樹新，等，天然氣氣瓶快充過程溫升數(shù)值模擬研究[J].高校化學(xué)工程學(xué)報(bào)，2011，25（6）：966-970.

[9]王光緒，周劍秋，胡淑娟，等，高壓儲(chǔ)氫氣瓶的長徑比和進(jìn)氣口直徑對(duì)其快充溫升及溫度分布的影響[J].工程力學(xué)，2014.31（5）：225-232.

[10]楊剛，黃思，李蔚，等.車載LNG氣瓶在高溫下熱力學(xué)參數(shù)的變化研究[J].機(jī)電工程技術(shù)，2020，49（10）：21-23.

[11]夏莉，勞英杰，李蔚，等.車用LNG氣瓶升壓規(guī)律的研究[J].制冷與空調(diào)（四川），2015，29（6）：741-743.

[12] Technical Committee CEN/TC. Cryogenic vessels-Methods for-performance evaluation of thermal insulation： BS EN12213-1999[S]. BSI， London. 1999.

[13]道客巴巴，氫氣物性參數(shù)[EB/O L]. http：//www. doc88.com/p-6846670985306.html.

[14] Shih T M. Numerical heat transfer [M]. Boca Raton： CRC Press.1984.

[15]劉峻，周連勇，馬華慶，等.不同水容積的車載儲(chǔ)氫氣瓶快充溫升研究[J].化工機(jī)械，2021，48（2）：160-165.

[16]鄭津洋，別海燕，徐平，等.車用纖維全纏繞高壓儲(chǔ)氫氣瓶標(biāo)準(zhǔn)研究[J].壓力容器，2007（11）：48-56.

[17]白文怡，廖英強(qiáng)，劉勇瓊.等.碳纖維復(fù)合材料纏繞氣瓶優(yōu)化研究進(jìn)展[J].炭素，2020（1）：36-41.

[18]益小蘇，杜善義，張立同，等，復(fù)合材料工程/中國材料工程大典第10卷[M].北京：化學(xué)工程出版社，2006.

[19]劉格思.70MPa車用儲(chǔ)氫氣瓶快充溫升研究[D].杭州：浙江大學(xué)。2012.