目录图表
1 一体化设计

1 一体化设计

2 固态储氢容器的设计

2 固态储氢容器的设计

3 基于ANSYS的固态储氢容器应力分析

3 基于ANSYS的固态储氢容器应力分析

3.1 工作条件和结构参数

3.1 工作条件和结构参数

3.2 网格划分及无关性验证

3.2 网格划分及无关性验证

3.3 载荷施加

3.3 载荷施加

3.4 结果分析

3.4 结果分析

4 结论

4 结论

参考文献:

参考文献:

基于固态稀土储氢的储罐与重卡底盘一体化设计

郭 天磊1杨 宇宸1侯 雄坡1宋 航1杨 福胜1张 早校1吴 震1王 天亮1,2
1. 西安交通大学化学工程与技术学院,陕西 西安 710049
2. 中国一重集团有限公司,黑龙江 齐齐哈尔 161041
摘要: 氢能被认为是减少碳排放和解决化石能源枯竭的清洁能源之一。然而,储氢技术是发展氢能的瓶颈。固态储氢具有质量储氢密度高、储氢压力低和安全性高的优势,是一种理想的储氢方式。在低压和常温下固态稀土合金LaNi5具有抗毒性好、易活化和吸放氢速度快等优点,但质量储氢密度比其他常温固态储氢合金的低。将应力低且均匀的AB5型等固态储氢容器有效地适应重卡结构,实现储氢容器的异形化及与载具底盘车架的一体化设计。同时,采用轻质合金设计储氢容器,实现一体化储氢系统的轻量化。降低AB5型储氢合金的质量储氢密度低对一体化储氢系统带来的不利影响。设计Q345R、钛合金、6061铝和LA103Z等材质的储氢容器。同一工况下,容器筒体的厚度分别为6,4,7和11 mm;容器的质量分别为24.4,9.2,9.9和8.7 kg;系统的质量储氢密度分别为1.04%,1.23%,1.22%和1.24%(质量分数)。其中LA103Z材质储罐的一体化储氢系统质量最轻,储氢密度最高,低质量储氢密度的稀土储氢合金对其一体化系统储氢密度的不利影响最小。与现阶段主流的35 MPa铝内胆碳纤维全缠绕复合材料储氢瓶相比,一体化储氢系统可为重卡节省62.5 L·kg−1 H2的载货空间。针对重质Q345R储罐,进行基于稀土镧镍系,钛铁系,钛锰系,钒基固溶系等常温储氢材料一体化储氢系统的对比分析。基于钛锰系合金一体化储氢系统的储氢量(30.5 kg)和储氢密度(1.36%)最高,基于钒基固溶系合金的一体化储氢系统最有利于轻量化设计。固态储氢材料吸氢膨胀会对固态储氢容器内部产生较大的应力(膨胀应力),严重时会导致容器变形甚至开裂。建立考虑储罐内储氢材料吸氢体积膨胀和热膨胀的储氢容器应力应变模型。采用有限元法对固态储氢容器进行应力应变分析,综合考虑了储氢容器内部氢气压力,吸放氢过程中储氢合金热膨胀和体积膨胀产生的应力,储氢容器外表面的受力。分析结果表明,Q345R、钛合金、6061铝和LA103Z储氢容器的等效应力最大值分别是116.09,173.16,130.19和100.42 MPa,均未超过材料的抗拉强度和屈服强度;等效弹性应变最大值分别是0.000599,0.00156,0.00213和0.00246 mm·mm−1,应变最大值均在微米级。
关键词: 固态稀土储氢容器;一体化设计;储氢密度;应力;储氢系统
基金: 国家自然科学基金项目(52176203);陕西省重点研发计划项目(2024CY2-GJHX-19);长庆工程设计有限公司项目(2025-CQSJ-F008)
收稿日期: 2024-01-23

交通行业消耗了中国一半的石油消费量,并向大气中排放了10%的二氧化碳,其中重型卡车碳排放占汽车排放总量近一半。氢能作为21世纪最理想的清洁能源之一,发展氢燃料电池是支撑国家“双碳”战略目标的关键技术[1–5]。氢燃料电池能量密度高,电池功率大,可以解决纯电动重卡路线能量密度不足和混动重卡路线经济效益不突出的问题[6],故氢燃料电池重卡是解决全球能源问题和促进实现低碳经济的理想方案之一[7–9]。金属氢化物固态储氢技术具有较高的体积/质量储氢密度、较好的安全性以及高能效等优点,是一种极为理想的储氢方式,有望在重卡储氢系统中广泛应用,提高储氢量和安全性[10–14]。

综上所述,本文介绍了稀土固态储氢与重卡底盘车架一体化储氢系统的设计思路,分别选用Q345R、钛合金、6061铝和LA103Z金属材质设计储罐,实现固态储氢容器的轻量化。针对重质Q345R储罐,进行基于稀土镧镍系、钛铁系、钛锰系和钒基固溶系等常温储氢材料储氢系统的对比分析。建立考虑储罐内储氢材料吸氢体积膨胀和热膨胀的储氢容器应力应变计算模型,采用有限元法对固态储氢容器的应力和应变进行评估与分析。

1 一体化设计

在有限体积和质量的车载储氢系统中,提高体积储氢密度和质量储氢密度,实现氢的高密度储存。金属氢化物作为固态储氢的典型代表,体积储氢密度达到了高压气态储氢的2~5倍。但对于镧镍系、钛铁系、固溶体钒基等传统的金属氢化物储氢材料而言,其质量储氢密度仅为1%~3.5%(质量分数,下同)。AB5型LaNi5固态储氢合金材料因其具有较高的安全性和体积储氢密度而受到广大学者的关注。同时,LaNi5储氢材料具有抗毒化性能好,易活化,反应速度快等特点,可在常温和低压下完成吸放氢过程[15]。王宏和刘祖岩[16]对LaNi5晶体结构进行理论分析计算,得出LaNi5的最大储氢量为1.379%。LaNi5合金存在易粉化、易氧化等缺点,影响其储氢性能和循环稳定性。研究发现,LaNi5在O2,H2O,CO或SO环境中其吸氢能力显著降低,反应动力学变慢[17]。李珂等[18]采用溶液浸渍法对LaNi5合金粉末进行了表面包覆改性,所用包覆材料为聚甲基丙烯酸甲酯(PMMA),空气氧化后,PMMA包覆LaNi5合金的储氢量保持在1.29%,而未包覆LaNi5合金的储氢量仅为1.15%。

目前单一类型的AB5型LaNi5储氢合金质量储氢密度低,储存等量的氢气需要大量的储氢材料,将会增大储氢系统的重量,制约AB5型储氢合金在交通行业的广泛应用。由德国宝马(即巴伐利亚发动机制造厂股份有限公司)、罗伯特·博世有限公司、Hexagon Purus和Testnet Engineering等组成的“FlatHyStor——创新氢气罐系统的功能设计和测试”联盟合作成功开发了适用于燃料电池汽车底部的扁平高压氢气罐,增大了载具的储氢密度和续航能力。但目前典型的高压气态储氢容器为圆柱形或球形,很难与重卡底盘融为一体。LaNi5型固态储氢具有储氢压力低,吸放氢温度室温,安全性好等优势。应力低且均匀的AB5型等固态储氢容器能够更有效地适应重卡结构,可以实现储氢容器的异形化及与载具底盘车架的一体化。同时选用轻质罐材,实现储氢系统的轻量化。降低AB5型固态储氢材料质量储氢密度低所带来的不利影响。

重卡底盘车架广泛采用梯形车架(又称边梁式车架),由若干根横梁和相互平行的纵梁组成[19]。梯形车架弯曲刚度较大,当承载扭矩时,各部分同时产生扭转和弯曲。其优势是便于安装车厢、车身和布置其他总成,车辆的改型和改装。因此也被广泛地应用在越野汽车,特种汽车,客车等。重卡底盘梯形车架由两条开口朝内相互平行、冲压制成的槽型纵梁及若干开口槽型横梁组合而成。等长等宽车架简化了纵梁的冲压工艺并在横梁上不会产生附加扭矩。目前大多数氢能和纯电动重卡车都是在不改变燃油车的布置基础上,将储氢容器或电池通过框架的方式直接安装在驾驶室的后方,占用了重卡额外的载货空间。同时重卡的重心变高,而且,框架遮挡住了驾驶室后方的视野,存在巨大的安全隐患[20]。

固态储氢与重卡底盘车架一体化储氢装置通过把储氢容器安装在重卡底盘闲置的空间内,未占用重卡额外的载货空间。同时,降低了重卡的重心,减小重卡高速行驶时受到的离心力和对轮胎的磨损,提高重卡的防侧翻能力,延长轮胎的使用寿命。因此,一体化装置更安全可靠。固态储氢与重卡底盘车架一体化储氢装置如图1所示,该一体化装置包括集成式储氢容器、外部框架和冷却装置。外部框架由底盘车架纵梁,底盘车架横梁,受力骨架底板和受力骨架盖板组成。重卡底盘的横梁和纵梁对集成式储氢容器起固定和支撑作用。外部框架可以有效抵御行驶过程中对储氢容器的各种振动与冲击。储氢容器,底盘车架纵梁,底盘车架横梁,受力骨架底板和受力骨架盖板围成的腔体为冷却液流道。通过四重高效交叉换热单元精确调控储氢容器温度,保障吸放氢环节的正常运行。同时,一体化装置在底盘集成多个固态储氢容器,提高了重卡底盘的利用率和重卡的续航能力。

图1 固态储氢与重卡底盘车架一体化储氢装置
Fig.1 Solid-state hydrogen storage and heavy truck chassis frame integrated hydrogen storage device

图2为一体化装置管道分布,可知,在重卡横梁上设有冷却液进口,冷却液出口,加氢口和放氢口。储氢容器两端采用焊接的方式固定在重卡底盘横梁上,储氢容器直径与重卡底盘车架高度一致。装置整体的表面平齐,便于安装和集成。在相邻的固态储氢容器和换热通道之间设置接管,便于调节氢气和换热工质流动,保证高效的吸放氢速率和换热效果。

图2 一体化装置管道分布
Fig.2 Distribution of piping of integration device
(a) Front;(b) Back

2 固态储氢容器的设计

根据重卡底盘车架尺寸,各附件的布置,车架相邻纵梁横梁的间距,为了最大程度地利用重卡底盘空间,储氢容器的筒体长度设计为800 mm。重卡完全氢化时,计划每个固态储氢容器可以储存1.035 kg的氢气,在重卡底盘车架内部空间设计28个大小相同并紧凑排列的储氢容器。一体化储氢系统的总储氢达到29 kg。在Sato等[21]的研究中,当温度为303 K,压力低于1 MPa时LaNi5的储氢量为1.4%。LaNi5储氢容器的设计温度为323 K,工作压力为5 MPa,均高于Sato等的实验工况,材料的储氢密度将达到1.4%。此时每个储氢容器的内径可由式(1)求得: D i = 4 M 1 L w 1 π ( 1 - φ ) ρ a α 1 (1)式中,Di为储氢容器的内径,m;M1为单个储氢容器的储氢质量,kg;L为储氢容器筒体长度,m;α1为储氢材料的填充率;ϕ为储氢材料的孔隙率;ρa为储氢材料未吸氢时的密度,kg·m−3;w1为LaNi5储氢材料的质量储氢密度,%。计算得储存1.035 kg的氢气时储氢容器的内径为200 mm。

设计Q345R、钛合金、6061铝和LA103Z材质储氢容器。在固态储氢容器内部进行涂层处理,隔离储氢容器材料与氢气的直接接触。Q345R是屈服强度为340 MPa级的压力容器专用钢板,它具有良好的综合力学性能和制造工艺性能,在中国压力容器行业用量较大、用途较广。但Q345R具有密度大、耐腐蚀性差等缺点,钢制储氢容器容易被氢气腐蚀出现氢脆现象,导致容器在压力作用下失效,出现泄露等风险[22]。钛合金对氧化性、弱还原性、中性等介质耐腐蚀,具有重量轻、强度高、低温性能好、粘附力小等优点。在长寿命、强腐蚀性介质、轻质的储氢设备中应用较广。然后,钛合金导热系数低、化学性质活泼,在切削加工过程中,刀具附件区域不易散热和冷却易形成高温区,会对刀具造成物理和化学损坏,加工性能差。

铝合金的设计压力小于16 MPa,设计温度范围在–269~200 ℃。铝合金的密度轻,仅为钢的三分之一,具有良好的韧性、塑性、成型和焊接性能,可应用于低压储氢容器[22–23]。镁锂合金具有超轻、高比强度、高比模量等特点。温室常压下镁锂合金的刚性是钢铁的22倍,满足同等刚度所需材料重量是钢材的三分之一。而且,镁锂合金的力学性能均匀一致,可实现少甚至无切削加工,机械加工余量小,提高了生产效率。同时,镁锂合金具有优异的可焊接性能。然而,镁锂合金的价格昂贵,生产成本高。表1给出了以上4种储氢容器材料的力学性能参数[22–25]。

表1 不同储氢容器材料的力学性能参数
Table 1 Mechanical property parameters of different hydrogen storage container materials[22–25]
Materials Tensile strength (Rm)/MPa Yield strength(Rel)/MPa Allowable stress/MPa
Q345R 510 345 189
Titanium alloy 400 245 169
6061Al 262 241 97
LAl03Z 162 154 60

考虑焊接可能引起的强度削弱,承受内压单层薄壁圆筒计算厚度由式(2)求得: δ = P c D i 2 [ σ ] t φ − P c (2)式中,δ为计算厚度,mm;Pc表示计算压力,MPa;φ为焊接接头系数;[σ]t为设计温度下的许用应力,MPa。焊接接头采用双面焊对接接头和相当于双面焊的全熔透对接接头,进行局部无损检测,焊接接头系数取值为0.85。

当已求得圆筒尺寸Di,δn和δe时,按式(3)对圆筒进行应力强度校核: σ t = P c ( D i + δ e ) 2 δ e ≤ [ σ ] t φ (3)式中,δe为有效厚度,mm,δe=δn–C;δn为名义厚度,mm;C为附加厚度,mm;σt表示设计温度下圆筒的计算应力,MPa。

由式(4)求得设计条件下储氢容器的最大运行工作压力[Pw]: [ P w ] = 2 δ e [ σ ] t φ D i + δ e (4)

采用式(5~7)计算储氢容量的质量和系统质量储氢密度: M 2 = π [ ( D 0 2 ) 2 − ( D i 2 ) 2 ] L ρ (5) M 3 = π ( D i 2 ) 2 L ρ a ( 1 − φ ) α 1 (6) w = M 1 M 1 + M 2 + M 3 × 100 % (7)式中,M2为单个储氢容器的质量,kg;M3为固态储氢材料的质量,kg;D0为储氢容器筒体的外径,m;ρ为储氢容器材料的密度,kg·m−3:w为储氢装置的质量储氢密度,%。

不同固态储氢容器的设计结果以及不同固态储氢容器的重量及系统储氢密度如表2和图3所示。当重卡固态储氢容器使用传统的钢材时,钢制储氢容器耐腐蚀性差,容易被氢气腐蚀出现氢脆现象,同等工况下所需的厚度较大。4种固态储氢容器材料中,Q345R的密度最大,是其他3种材料密度的1.5~5倍。采用Q345R储罐储存29 kg氢气时,储氢系统的总质量(氢气质量、储氢合金质量和储氢罐质量)最重将达2780 kg,系统的质量储氢密度最低仅为1.04%。采用Q345R设计储氢容器增加了容器的厚度、重量和成本。钛合金具有耐腐蚀、高强度等优点,钛合金储氢容器最薄,仅为4 mm,但钛合金储氢容器加工制造困难,在工业领域应用较少。镁锂合金的密度仅为1500 kg·m−3,约为钢铁的五分之一。在设计条件下4种材料中镁锂合金的许用应力最小,承受同等氢气压力和载荷时储氢容器最厚。但采用镁锂合金时车载储氢系统装备的重量最轻,约为钢制储氢容器质量的三分之一,一体化装置的质量储氢密度最高,达到1.24%。

表2 不同固态储氢容器的设计结果
Table 2 Design results of different solid-state hydrogen storage containers
Materials Thickness/mm Density/(kg·m−3) Container mass/kg Hydrogen storage density/%
Q345R 6 7850 24.4 1.04
Titanium alloy 4 4620 9.2 1.23
6061Al 7 2730 9.9 1.22
LAl03Z 11 1500 8.7 1.24
图3 不同固态储氢容器的重量及系统储氢密度
Fig.3 Weight and system hydrogen storage density of different solid-state hydrogen storage containers

现阶段主流的35 MPa铝内胆碳纤维全缠绕复合材料储氢瓶(Ⅲ型)的体积储氢密度约为16 g·L−1[26]。如果不进行一体化设计,由式(8)求得,固态储氢装置只比35 MPa的储氢瓶节省了38.1 L·kg−1 H2的载货空间。由式(9)求得,一体化储氢系统比35 MPa的储氢瓶可为重卡节省62.5 L·kg−1 H2以上的载货空间,增加了64.4%。 V 1 = 1000 V H 1 − 1000 V H 2 (8) V 2 = 1000 V H 1 (9)式中,V1表示若不进行一体化设计,固态储氢装置比35 MPa储氢瓶节省的载货空间,L·kg−1;V2表示一体化设计比35 MPa储氢瓶节省的载货空间,L·kg−1;VH1表示35 MPa储氢瓶的体积储氢密度,g·L−1;VH2表示固态储氢装置的体积储氢密度,g·L−1。

燃料电池重卡对自重极为敏感。针对Q345R储罐,系统储氢密度仅为1.04%。储氢密度较低的LaNi5材料限制了重质Q345R储氢系统储氢密度的提升。除此之外,钛铁系、钛锰系和钒基固溶系等合金均可在室温附近吸放氢,且可逆质量储氢密度可达到1.8%~2.5%。室温下TiFe合金可以吸收大量的氢气,其储氢量为1.8%,是一种含量丰富、成本低廉的储氢材料。但TiFe基合金的动力学性能差,活化条件要求高[27]。TiMn2合金的吸放氢性能良好,储氢量约为2.0%[28]。钒基固溶系合金具有较高的储氢能力,但其氢化物过于稳定。这造成了钒基固溶系合金在常温常压下难以放氢,实际可逆储氢容量在1.9%左右[29–30]。

图4为针对重质Q345R储罐,进行基于稀土镧镍系、钛铁系、钛锰系和钒基固溶系等常温储氢材料储氢系统的对比分析。由于重卡底盘的空间受到限制,一体化储氢系统的体积是定量,质量相关的量是变量。4种常温材料中,钛锰系合金的密度较轻,质量储氢密度最高。因此,基于钛锰系合金储氢系统的储氢量(30.5 kg)和储氢密度(1.36%)最高,系统装置较轻。钒基固溶系合金的密度最轻,最有利于储氢系统的轻量化,但限制了系统储氢量(体积储氢密度)的提升。稀土镧镍系合金的质量储氢密度最低,其密度明显比其他3种材料的密度高。因此,基于稀土镧镍系合金储氢系统的装置最重,系统储氢密度最低。但具有高密度、低储氢密度的稀土镧镍系合金对系统储氢量(体积储氢密度)的不利影响较低。

图4 基于稀土镧镍系、钛铁系、钛锰系和钒基固溶系的储氢系统对比分析
Fig.4 Comparative analysis of hydrogen storage systems based on rare earth lanthanide-nickel,titanium-iron,titanium-manganese,and vanadium-based solid solution systems

3 基于ANSYS的固态储氢容器应力分析

储氢容器受到的应力不仅由气体压力引起,还受吸放氢过程中储氢材料热膨胀和体积膨胀的影响。合金吸收氢气时体积膨胀,释放氢气后体积又会收缩。吸氢过程中储氢材料的体积膨胀会对储氢容器内壁面产生挤压。膨胀的合金粉末会在储氢容器内造成局部应力集中,应力过大时会导致储氢容器变形甚至失效,这对储氢容器的强度要求更高[31]。按照膨胀应力增加的速率可以将吸放氢循环分为3个阶段:膨胀应力缓慢上升期、膨胀应力快速上升期、膨胀应力稳定期。在膨胀应力稳定期,膨胀应力达到峰值且趋于动态平衡。此时,膨胀应力占固态储氢容器总应力的比重最大,将达到60%以上[31]。根据Nasako等[32]的实验,在第120~140次吸放氢循环中,储氢容器表面的膨胀应力达到260 MPa以上,该表面应力超过了304不锈钢的屈服强度,容器发生了塑性变形。因此,考虑吸放氢过程中膨胀应力对储氢系统的影响十分有必要。由于氢化物体积随氢浓度线性增加,在氢化过程中考虑LaNi5合金热膨胀和体积膨胀的有效膨胀系数αe为[33]: α e = α t [ 1 + α t Δ T + 0.167 3 α t Δ T ( 1 − [ ρ sat − ρ s ρ sat − ρ α ] ) (10)式中,αt为LaNi5的热膨胀系数,mm−1·K−1;ΔT为温差,K;ρs表示合金材料正在吸附氢时的密度,kg·m−3;ρsat表示合金材料吸附氢饱和时的密度,kg·m−3。

考虑储氢合金材料的孔隙率,计算结果表明计算多孔固体的有效杨氏模量时,采用Phani-Niyogi关系式[34],由有效膨胀系数和有效杨氏模量得出吸放氢过程中合金颗粒膨胀产生的应力: Y e = Y 0 ( 1 − φ φ cr ) p (11) P e = a e k Y e Δ φ (12)式中,k为杨氏模量修正系数;Y0为多孔介质材料的杨氏模量,GPa;Δϕ为孔隙率差;ϕcr为杨氏模量为零时的临界孔隙率;Ye为LaNi5材料的杨氏模量,GPa;Pe为膨胀应力,GPa;p为材料特定指数。LaNi5储氢材料的热物理性能如表3所示[35]。

表3 LaNi5储氢材料的热物理性能
Table 3 Thermophysical properties of LaNi5 hydrogen storage material[35]
Parameters Values
at/(mm−1·K−1) 1.23×10−5
ΔT/K 20
Y0/GPa 140
p 2.1
ϕ 0.5
ϕcr 1

在强度分析中常采用第三强度理论(最大切应力理论)或第四强度理论(最大畸变能密度理论)校核,吸放氢过程的交变载荷疲劳分析不在本文研究范畴。第三强度理论虽然可以较好表现塑性材料塑性屈服现象,但未考虑主应力的影响,故只适用于拉伸屈服极限和压缩屈服极限相同的材料。第三强度理论主要用于设计中高压压力容器,该方法设计的压力容器安全性高。第四强度理论考虑了主应力的影响,比第三强度理论更接近实际情况,所以在强度评价中多用第四强度理来评价。ANSYS后处理中Mises等效应力遵循的就是第四强度理论,无论材料处于何种应力状态,按第四强度理论材料不发生破坏的强度条件为: σ eq = 1 2 [ ( σ 1 − σ 2 ) 2 + ( σ 2 − σ 3 ) 2 + ( σ 3 − σ 1 ) 2 ] ≤ [ σ ] (13)式中,σeq为等效应力,MPa;σ1,σ2和σ3分别为第一、二和三主应力,MPa;[σ]为许用应力,MPa。

3.1 工作条件和结构参数

图5是固态储氢容器的设计示意图,利用Space Claim建立固态储氢容器三维模型。材料的弹性模量(E)和泊松比(μ)如表4所示[22–25]。

图5 固态储氢容器设计示意图
Fig.5 Schematic diagram of design of solid-state hydrogen storage container
表4 不同储氢容器材料的弹性模量和泊松比
Table 4 Elasticity modulus and Poisson´s ratio of different hydrogen storage container materials[22–25]
Materials Elasticity modulus/GPa Poisson's ratio
Q345R 210 0.30
Titanium alloy 114 0.36
6061Al 71 0.33
LAl03Z 43 0.30

3.2 网格划分及无关性验证

固态储氢容器的几何结构及承受的载荷左右对称,故受载后的应力和应变也对称。只需建立一半部分的固态储氢容器有限元模型,简化分析过程,提高求解速度和精度。网格划分是对固态储氢容器进行有限元分析极其重要的一步。单元网格的面积越小,网格划分的就越细,计算出的分析结果就会越准确。在保证结果准确的基础上,对网格进行无关性验证,采用最少的网格数仿真计算。图6为不同储氢容器材料的网格无关性验证,可以看出,当网格数分别为26.4万、25万、79万和78万时,Q345R、钛合金、6061铝和LA103Z固态储氢容器的等效应力峰值趋于稳定,选取此网格划分方式进行后续的仿真分析。

图6 不同储氢容器材料的网格无关性验证
Fig.6 Grid-independent verification of different hydrogen storage container materials

3.3 载荷施加

ANSYS有限元分析是为了了解对模型施加载荷时,模型的各个响应。正确且有效地施加载荷并进行仿真,是得到真实仿真结果的关键一环。车厢底部与重卡底盘车架处于“面-面”接触状态,车厢与车架的刚度共同决定了接触压力的大小与分布形式。首先对车架净载荷进行定义,车厢货物简化为均布载荷,方向是垂直车架向下,直接施加在车厢与重卡底盘纵梁、固态储氢容器的接触面上。每个固态储氢容器受到的车厢货物的压力由式(14~17)求得: A 1 = 2 L 1 B 1 (14) A 2 = n L 2 B 2 (15) F = G A 2 n ( A 2 + A 2 ) (16) p = F L 2 B 2 (17)式中,G为车厢(包括货物)总重量,N;A1为车厢与车架的接触面积,mm2;L1为车厢与车架的接触长度,mm;B1为纵梁的宽度,mm;A2为车厢与储氢容器罐的接触总面积,mm2;L2为车厢与单个储氢容器罐的接触长度,mm;B2为车厢与单个储氢容器罐的接触宽度,mm;n为固态储氢容器罐的数量。根据实际尺寸、设计参数及式(14~17),可得到各项数据,结果如表5所示。

表5 车厢与货物载荷施加数据
Table 5 Carriage and cargo load applied data
Parameters Values
L1/mm 5700
L2/mm 800
B1/mm 100
B2/mm 4.4
n 20
G/N 196000
F/N 566
P/MPa 0.16

对固态储氢容器内表面施加5 MPa的氢气压力和吸放氢循环初期产生的膨胀应力,在固态储氢容器与盖板接触面上施加566 N的力。固态储氢容器底部封头处设置为固定约束,对称面截面处设置为无摩擦支撑。

3.4 结果分析

图7为不同固态储氢容器等效应力分布云图,可以看出,固态储氢容器筒体部分的应力变化比较平缓,封头处的应力值最小。接管和筒体与封头连接部分结构尺寸突变,存在应力集中现象,这些区域应力值最大。在设计时可以增大该区域厚度值或者减缓结构尺寸的突变过程以减小应力值。Q345R、钛合金、6061铝和LA103Z储氢容器的等效应力最大值分别是116.09,173.16,130.19和100.42 MPa,均未超过材料的抗拉强度和屈服强度。4种固态储氢容器中,钛合金固态储氢容器筒体处的应力值最大在110 MPa左右,LA103Z固态储氢容器筒体处的应力值最小仅在40 MPa左右,钢制和6061铝储氢容器筒体处的应力值接近,都是在70 MPa附近。筒体处的等效应力值均小于材料的许用应力,具有较大的安全系数。

图7 不同固态储氢容器等效应力分布云图
Fig.7 Cloud picture of the equivalent stress distribution of different solid-statehydrogen storage containers
(a) Q345R;(b) Titanium alloy;(c) 6061Al;(d) LA103Z

图8为不同固态储氢容器等效弹性应变分布云图。由于容器受到均匀氢气内压和膨胀应力作用,筒体向外均匀膨胀。同时筒体受到车厢的压力作用,筒体发生弯曲变形。因此,在筒体与车厢的接触区域有较大的变形。由图8可知,Q345R、钛合金、6061铝和LA103Z储氢容器的等效弹性应变最大值分别是0.000599,0.00156,0.00213和0.00246 mm·mm−1。钢制固态储氢容器的等效弹性应变值最小,钛合金、6061铝和LA103Z制固态储氢容器的等效弹性应变云图相似。反应器的等效弹性应变值都很小,应变最大值均在微米级。

图8 不同固态储氢容器等效弹性应变分布云图
Fig.8 Cloud picture of equivalent elastic strain distribution of different solid-state hydrogen storage containers
(a) Q345R;(b) Titanium alloy;(c) 6061Al;(d) LA103Z

为了获得沿固态储氢容器轴向和径向各个关键点的具体应力值,直观体现固体储氢容器应力变化情况。在筒体中径面和直径方向上创建路径,等效应力沿固态储氢容器轴向和径向变化趋势如图9和10所示。在筒体中部和封头与筒体连接处创建截面,截面处的应力分布云图如图11所示。固态储氢容器筒体部分的应力变化比较平缓,在筒体与封头连接处的地方结构尺寸突变,所以会出现应力剧增的情况,筒体最大应力值出现在固态储氢容器的中部。其中LA103Z制储氢容器轴向应力值最小,最大处应力值仅在40 MPa左右,而钛合金制储氢容器轴向应力的最大值达到了110 MPa以上,是LA103Z的两倍以上。Q345R、钛合金和6061铝固态储氢容器筒体处径向应力值呈缓慢下降趋势。筒体内壁受氢气压力作用,所以在筒体内壁处应力值最大,外壁处应力值最小。LA103Z储氢容器应力值在径向呈缓慢上升趋势,筒体外壁面应力值比内壁面的大,这与镁锂合金低密度、高比强度、高比刚度的特点相关。

图9 不同固态储氢容器等效应力沿轴向的变化趋势
Fig.9 Trend of equivalent stresses in axial direction of different solid-state hydrogen storage containers
图10 不同固态储氢容器等效应力沿径向的变化趋势
Fig.10 Trend of equivalent stresses in radial direction of different solid-state hydrogen storage containers
图11 不同固态储氢容器各个截面应力分布云图
Fig.11 Cloud figure of stress distribution in each cross-section of different solid-state hydrogen storage containers

4 结论

本文设计了固态储氢容器与重卡底盘车架一体化储氢系统,探讨采用不同金属材质的储氢容器对一体化储氢系统的质量储氢密度和重量的影响。针对重质Q345R储罐,进行基于稀土镧镍系、钛铁系、钛锰系和钒基固溶系等常温储氢材料一体化储氢系统的对比分析。并建立考虑储罐内储氢材料吸氢体积膨胀及热膨胀的储氢容器应力应变模型,得出主要结论如下:

1)通过将应力低且均匀的AB5型等固态储氢容器有效地适应重卡结构,设计储氢容器与重卡底盘车架一体化储氢系统。如果不进行一体化设计,与35 MPa气瓶相比,固态储氢装置只为重卡节省了38.1 L·kg−1 H2的载货空间。一体化储氢系统比35 MPa气瓶为重卡节省了62.5 L·kg−1 H2的载货空间,增加了64.4%。

2)设计Q345R、钛合金、6061铝和LA103Z等材质的储氢容器。同一工况下,储氢容器的厚度分别为6,4,7和11 mm;储氢容器的质量分别为24.4,9.2,9.9和8.7 kg;系统的质量储氢密度分别为1.04%,1.23%,1.22%和1.24%。LA103Z材质储罐的一体化储氢系统质量最轻,储氢密度最高。低质量储氢密度的稀土储氢合金对其一体化系统储氢密度的不利影响最小。

3)针对Q345R储罐,基进行于稀土镧镍系、钛铁系、钛锰系和钒基固溶系等常温储氢材料一体化储氢系统的对比分析。基于钛锰系合金储氢系统的储氢量(30.5 kg)和储氢密度(1.36%)最高,系统装置较轻。钒基固溶系合金的密度最轻,最有利于储氢系统的轻量化(2210 kg),但限制了系统储氢量的提升。基于稀土镧镍系合金储氢系统的装置最重(2780 kg),系统储氢密度最低(1.04%)。但具有高密度、低储氢密度的稀土镧镍系合金对系统储氢量的不利影响较小。

4)由于接管与封头之间的不连续性发生应力集中,最大应力出现在接管和封头的连接处。Q345R、钛合金、6061铝和LA103Z储氢容器的等效应力最大值分别是116.09,173.16,130.19和100.42 MPa,均未超过材料的抗拉强度和屈服强度。在设计时可以增大该区域厚度值或者减缓结构尺寸的突变过程。Q345R、钛合金、6061铝和LA103Z储氢容器的等效弹性应变最大值分别是0.000599,0.00156,0.00213和0.00246 mm·mm−1,应变最大值均在微米级。

参考文献:

[1] 

(李航,胡尊严,胡家毅,董家臣,李建秋,徐梁飞,欧阳明高,卜玉,王力军,秦志东. 新型分布式驱动液氢燃料电池重型商用车设计、分析与验证 [J]. 汽车工程,2022,44(8):1183.)

[2] 

(吴岱丰,李睿,陈鹏允,周庆,唐仁衡,肖方明. Ti-Mn基高容量固态低压储氢合金研究 [J]. 稀有金属,2023,47(3):357.)

[3] 

(刘友成,周兵,刘绍忠. NiF2掺杂对MgH2解氢性能影响机制研究 [J]. 稀有金属,2018,42(11):1193.)

[4] 

Mohammadshahi S S,Gray E M,Webb C J. A review of mathematical modeling of metal-hydride systems for hydrogen storage applications [J]. International Journal of Hydrogen Energy,2016,41:3470.

[5] 

(高亚男,张全逾,魏彦杰,郑镭. 汽车用铝合金在搅拌摩擦加工下的组织改性与力学性能分析 [J]. 塑性工程学报,2024,31(3):200.)

[6] 

(刘世闯,孙桓五,王瑞鑫,李昊,张东光. 大功率型氢燃料电池重卡动力系统匹配设计 [J]. 汽车工程,2021,43(2):196.)

[7] 

Thomas J M,Edwrds P P,Dobson G P. Decarbonising energy:the developing international activity in hydrogen technologies and fuel cells [J]. Journal of Energy Chemistry,2020,51(12):405.

[8] 

(雷超,李韬. 碳中和背景下氢能利用关键技术及发展现状 [J]. 发电技术,2021,42(2):207.)

[9] 

(马义,何特立,刘利连,李学锐,张剑. 燃料电池系统氢循环部件性能测试 [J]. 汽车科技,2023,(1):72.)

[10] 

(黄妹青,徐文强,李先明,杨康,何广利,蒋利军. 金属储氢材料的性能边界条件分析 [J]. 稀有金属,2018,42(10):1046.)

[11] 

Kaur M,Pal K. Review on hydrogen storage materials and methods from an electrochemical viewpoint [J]. Journal of Energy Storage,2019,23:234.

[12] 

Yang H,Ding Z,Li Y T,Wu P K,Hou Q H,Zhang Y,Gao B,Hou K F,Du W J,Shaw L L. Recent advances in kinetic and thermodynamic regulation of magnesium hydride for hydrogen storage [J]. Rare Metals,2023,42(9):2906.

[13] 

(徐忆明. 铜冶炼产品碳足迹核算方法与优化路径 [J]. 有色金属工程,2024,14(8):67.)

[14] 

(张志芸,张国强,刘艳秋,康启平. 车载储氢技术研究现状及发展方向 [J]. 油气储运,2018,37(11):1207.)

[15] 

Boateng E,Chen A C. Recent advances in nanomaterial based solid-state hydrogen storage [J]. Materials Today Advances,2020,6:100022.

[16] 

(王宏,刘祖岩. LaNi5最大储氢量的晶体结构分析 [J]. 稀有金属材料与工程,2004,(3):239.)

[17] 

Yang F S,Chen X Y,Wu Z,Zhang Z X,Wang S M,Wang G X,Wang Y Q. Experimental studies on the poisoning properties of a low-plateau hydrogen storage alloy LaNi4.3Al0.7 against CO impurities [J]. International Journal of Hydrogen Energy,2017,42(25):16225.

[18] 

(李珂,高云鹤,黄通文,周承商. PMMA表面包覆LaNi5合金粉末的储氢性能 [J]. 粉末冶金材料科学与工程,2022,27(6):595.)

[19] 

(翁宇. YT3621矿用车车架的有限元分析与改进设计 [D]. 西安:长安大学,2013. 1.)

[20] 

(高超,薛正成,黄露,湛时时. 一种低重心氢能重卡底盘 [P]. 中国:CN 216886261U,2022.)

[21] 

Sato T,Saitoh H,Utsumi R,Ito J,Nakahira Y,Obana K,Takagi S,Orimo S I. Hydrogen absorption reactions of hydrogen storage alloy LaNi5 under high pressure [J]. Molecules,2023, 28(3):1256.

[22] 

(郑津洋,桑芝富. 过程设备设计第五版 [M]. 北京:化学工业出版社,2021. 1.)

[23] 

(武磊,何兵,覃文东,覃铭. 异速轧制和常规轧制对7075铝板力学性能的影响 [J]. 有色金属科学与工程,2024,15(3):400.)

[24] 

(中华人民共和国国家质量监督检验检疫总局,中国国家标准化管理委员会. 压力容器:GB/T 150—2011 [S]. 北京:中国标准出版社,2011.)

[25] 

(中华人民共和国国家质量监督检验检疫总局,中国国家标准化管理委员会. 锅炉和压力容器用钢板:GB/T 713—2014 [S]. 北京:中国标准出版社,2014.)

[26] 

(徐双庆,范志超,刘孝亮. 车载高压-固态复合储氢系统储氢密度的数值分析 [J]. 压力容器,2020,9:24.)

[27] 

Zeaiter A,Nardin P,Yazdi M A P,Billard A. Outstanding shortening of the activation process stage for a TiFe-based hydrogen storage alloy [J]. Mrs Bulletin,2019,112:132.

[28] 

Zhan L J,Zhou P,Xiao X Z,Cao Z M,Piao M Y,Li Z N,Jiang L J,Li Z P,Chen L X. Numerical simulation and experimental validation of Ti0.95Zr0.05Mn0.9Cr0.9V0.2 alloy in a metal hydride tank for high-density hydrogen storage [J]. International Journal of Hydrogen Energy,2022,47(91):38655.

[29] 

Tiwari G P,Krishnamurthy N,Kumar S. Tailoring the hydrogen desorption thermodynamics of V2H by alloying additives [J]. Journal of Alloys and Compounds,2015,645(1):S252.

[30] 

(张晓飞,蒋利军,叶建华,武媛方,郭秀梅,李志念,李海文. 固态储氢技术的研究进展 [J]. 太阳能学报,2022,43(6):345.)

[31] 

Lin C K,Chen Y C. Effects of cyclic hydriding-dehydriding reactions of LaNi5 on the thin-wall deformation of metal hydride storage vessels with various configurations [J]. Renewable Energy,2012,48:404.

[32] 

Nasako K,Ito Y,Hiro N,Osumi M. Stress on a reaction vessel by the swelling of a hydrogen absorbing alloy [J]. Journal of Alloys and Compounds,1998,264:271.

[33] 

Joubert J M,Cerny R C,Latroche M. Compressibility and thermal expansion of LaNi5 and its substitutional derivatives [J]. Intermetallics,2005,13(2):227.

[34] 

Phani K K,Niyogi S K. Young modulus of porous brittle solids [J]. Journal of Material Science,1987,22(1):257.

[35] 

Dinachandran L,Mohan G. Numerical simulation of the parametric influence on the wall strain distribution of vertically placed metal hydride based hydrogen storage container [J]. International Journal of Hydrogen Energy,2015,40(16):5689.

顶部