摘要
质子交换膜燃料电池的边框一般通过黏结剂黏接在质子交换膜上,在保护质子交换膜的同时,还起到密封作用。燃料电池长期工作在复杂的环境下,边框和黏结剂的性能均会发生疲劳衰减,导致被剥离或破坏,进而影响燃料电池的气密性。因此,以边框和黏结剂的搭接试件为对象,结合内聚力模型和有限元软件设计了多种边框结构,研究了新结构对黏结性能的影响。结果表明:在低弹性模量的边框表面做三角形或梯形的凸形结构,将有效提升黏结结构的分离位移,但对于高弹性模量的边框,新结构反而会降低黏结性能;此外,优化结构的尺寸、形状对搭接模型黏结性能也有影响,当凸三角形底边长为1.5 mm时,搭接模型有最好的黏结性能。
在中国提出的“2030年前碳达峰,2060年前实现碳中和”的“双碳”目标指引下,质子交换膜燃料电池(PEMFC,proton exchange membrane fuel cell)凭借高效、清洁等优点,获得了极大的关注。燃料电池若要在车用方向实现大规模商业化,除了需要有较高的输出功率、较低的催化剂制作成本和氢气储运成本外,还需要有良好的耐久性和气密性,在极端工作条件下具有较长的使用寿命。目前,燃料电池密封结构的耐久性问题是限制燃料电池使用寿命的一块“短板”,为解决该问题,各密封部件的性能需要进一步提升。
燃料电池中的密封结构由密封圈和边框组成,其中,边框一般通过直接接触和粘贴绑定的方式组装在质子交换膜(PEM,proton exchange membrane)活性区域的周围。由于质子交换膜自身机械强度低,边框的存在可保护其不被过压,且与密封圈一起发挥密封作用,即阻止阴阳极之间的气体窜漏及防止气体泄漏至外界环境中。PEMFC中边框密封件若长期在高温、高湿、强酸性和氟离子存在的工作环境下,容易出现破裂,开胶等现象,这将严重降低电池的气密性,影响其正常工作,因此,若要提升燃料电池的气密性,对PEMFC中边框密封结构的研究则尤为关键。
目前,国内外对于PEMFC边框气密性的研究较少,在提升燃料电池密封性方面,大多数成果是结合有限元模型,优化密封
在对边框的研究之中,Ye
基于以上研究可知,合理选取边框材料及组装方式能够减小燃料电池组装、工作过程中在膜电极上产生的应力;此外燃料电池工作的温度、湿度、酸性条件会严重降低边框力学性能,影响使用寿命,而且湿度改变引发的质子交换膜膨胀收缩现象,将在边框密封件附近产生一定的集中应力,长期工作后会使边框和膜疲劳受损,引发边框密封件与膜开胶、质子交换膜破损等现象,降低燃料电池密封性甚至工作性能。因此,笔者以边框和黏结剂组成的试件为对象,通过建立力学模型,并设计不同的边框结构以提升边框密封件的黏结性能,来改善模型的力学表现,同时探究新的边框结构在不同边框材料上的表现,为边框材料选择提供理论指导意义。
典型的质子交换膜燃料电池由阴阳极端板、集流板、流场板和MEA通过螺栓紧固的方式组成,其中MEA包括质子交换膜、催化层和气体扩散层,它与密封结构的组成示意图如

图1 MEA密封结构与边框黏结剂试件示意图
Fig. 1 Schematic diagrams of MEA sealant structure and frame-adhesive test specimen
为了对试件结构进行准确的力学分析,需要对其中的黏结剂选择合理的力学模型。大多数文献中黏结剂的力学模型是基于Camanho等提出的双线性内聚力模
, |
, | (1) |
, | (2) |
式中:τ1和τ2是分离平面的面内方向剪切应力,τ3是分离平面法向应力,且法向位移δ3的方向决定了H(-δ3)的值,继而影响了在该方向上的受力大小;K是内聚力单元的初始刚度;d是损伤因子,当内聚力单元位移超过起始损伤位移δ后,内聚力单元的初始刚度将产生折损。混合模式下的位移由式(
, | (3) |
, | (4) |
, | (5) |
式中:δshear和δnormal分别表示面内剪切方向位移大小和平面法向位移大小;δm表示混合模式下合位移值。Gc是临界应变能量释放率,当形变过程中单元应变能累积超过该值后(或是单元位移超过失效位移δ后),单元将发生破坏,失去黏结作用,混合模式下内聚力单元的失效标准使用B-K律判断,为
, | (6) |
, | (7) |
式中:GIc和GIIc分别为内聚力单元平面上和法向上的临界应变能量释放率;r为模式比;η为B-K失效参数,在本模型中该值设置为

图2 混合模式下应力—位移曲线
Fig. 2 The stress-displacement curve in mixed-mode
基于以上黏结剂模型的建模原理,在ABAQUS中对试件进行建模仿真:模型示意图由3部分组成(如

图3 试件模型示意图
Fig. 3 Schematic diagrams of the specimen model
材料属性 | 值 | |
---|---|---|
边框 | 弹性模量E/MPa | 8 000 |
泊松比ν | 0.3 | |
内聚力单元 |
K/(N·m | 10 |
τ/MPa | 0.5 | |
τ/MPa | 0.5 | |
GIc/(N·m | 0.05 | |
GIIc/(N·m | 0.05 |
通过以上力学理论的应用和模型的搭建,对

图4 试件初始模型和优化模型
Fig. 4 Original model and optimized model of the specimen
图

图5 初始模型和优化模型的仿真计算结果
Fig. 5 The simulation results of the original model and the optimized model

图6 初始模型和优化模型牵引力—位移曲线对比图
Fig. 6 The traction-separation curves of the original model and optimized model
由于燃料电池中使用的边框及黏结剂材料种类较为多样,因此有必要对不同材料在加入优化结构后的黏结性能表现进行考察。

图7 不同边框材料、黏结剂属性的牵引力位移曲线
Fig. 7 The traction-separation curves of the models with different frames and adhesive
值得注意的是,

图8 不同临界应变能量释放率的模型牵引力位移曲线
Fig. 8 The traction-separation curves of the models with different energy release rates
对于2.1中凸三角形结构模型,三角形底边长1.5 mm,且两两之间间隔1 mm,若改变三角形底边长度(L),三角形之间的间距(D)也会发生相应变化,因此建立一系列不同三角形底边长度的模型,进一步探究优化结构的尺寸对黏结性能的影响。模型和主要的计算结果如
编号 | L/mm | D/mm | Fmax/N | Smax/mm |
---|---|---|---|---|
1-1 | 1.0 | 1.5 | 23.0 | 0.639 |
1-2 | 1.1 | 1.4 | 25.3 | 0.666 |
1-3 | 1.2 | 1.3 | 23.7 | 0.684 |
1-4 | 1.3 | 1.2 | 20.8 | 0.579 |
1-5 | 1.4 | 1.1 | 20.9 | 0.567 |
1-6 | 1.5 | 1.0 | 26.1 | 0.684 |
1-7 | 1.6 | 0.9 | 23.5 | 0.650 |
1-8 | 1.7 | 0.8 | 19.6 | 0.538 |
1-9 | 1.8 | 0.7 | 21.6 | 0.616 |

图9 不同凸三角形尺寸模型仿真结果
Fig. 9 The simulation results of models with different convex triangle sizes
为进一步提升试件的黏结性能,对黏结区域处边框的截面形状进行了改变,如

图10 边框截面形状优化结果
Fig. 10 The optimized results of the specimens with different section shapes
此外,凹三角形、凸梯形和凹梯形模型结构的尺寸对力学性能的影响亦进行了分析,其中对于相同底长(1.5 mm)的三角形和梯形结构,可将三角形视作梯形顶边长度为0的特殊形式,因此对于梯形结构,主要研究顶边(B)对模型力学性能的影响。凹三角形模型的尺寸参数和计算结果如
编号 | L/mm | D/mm | Fmax/N | Smax/mm |
---|---|---|---|---|
2-1 | 1.0 | 1.5 | 25.0 | 0.462 |
2-2 | 1.1 | 1.4 | 22.4 | 0.452 |
2-3 | 1.2 | 1.3 | 21.3 | 0.392 |
2-4 | 1.3 | 1.2 | 25.2 | 0.449 |
2-5 | 1.4 | 1.1 | 21.1 | 0.437 |
2-6 | 1.5 | 1.0 | 23.0 | 0.472 |
2-7 | 1.6 | 0.9 | 20.8 | 0.395 |
2-8 | 1.7 | 0.8 | 20.7 | 0.403 |
2-9 | 1.8 | 0.7 | 18.0 | 0.358 |
梯形形状 | 编号 | B/mm | Fmax/N | Smax/mm |
---|---|---|---|---|
凸梯形 | 3-1 | 0.2 | 23.6 | 0.620 |
3-2 | 0.4 | 20.7 | 0.614 | |
3-3 | 0.6 | 15.8 | 0.449 | |
凹梯形 | 4-1 | 0.2 | 17.6 | 0.369 |
4-2 | 0.4 | 18.2 | 0.355 | |
4-3 | 0.6 | 21.1 | 0.439 |
在ABAQUS中对边框—黏结剂试件建立了三维有限元模型,模拟了在相同牵引载荷作用下,不同边框结构的试件因黏结剂破坏而发生分离的情况,通过牵引力—位移曲线比较了它们黏结性能的差异,并对优化结构的尺寸和截面形状进行了进一步探索。结果表明引入凸形结构,能提升试件的分离位移,但会降低最大承载力;优化结构在低弹性模量的边框试件上能延长分离位移,但在高弹性模量的边框上有负面效果,而对不同属性的黏结剂均能提升分离位移;不同截面形状的边框结构对黏结性能也会产生影响,与其他截面形状相比,凸三角形的结构黏结性能最好,且当凸三角形底边长为1.5 mm时,优化结构的力学表现最佳。这对于提升边框结构在PEMFC中的密封性能有重要意义。
参考文献
张智明, 胡淞, 李昆朋, 等. PEMFC密封胶接触压力的均匀性研究及改善方法[J]. 湖南大学学报(自然科学版), 2018, 45(10): 30-37. [百度学术]
Zhang Z M, Hu S, Li K P, et al. Study on uniformity of contact pressure of PEMFCs sealant and its improvement method[J]. Journal of Hunan University(Natural Sciences), 2018, 45(10): 30-37.(in Chinese) [百度学术]
张智明, 张娟楠, 潘涛, 等. PEMFC金属双极板密封设计的有限元仿真[J]. 中国科技论文, 2017, 12(16): 1895-1901. [百度学术]
Zhang Z M, Zhang J N, Pan T, et al. Finite element method modeling on sealant design of PEMFC metal bipolar plate[J]. China Sciencepaper, 2017, 12(16): 1895-1901.(in Chinese) [百度学术]
肖文灵. 燃料电池密封结构优化设计与性能分析[D]. 大连: 大连理工大学, 2020. [百度学术]
Xiao W L. Performance analysis and optimization design of fuel cell sealing structure[D]. Dalian: Dalian University of Technology, 2020. (in Chinese) [百度学术]
Liang P, Qiu D K, Peng L F, et al. Structure failure of the sealing in the assembly process for proton exchange membrane fuel cells[J]. International Journal of Hydrogen Energy, 2017, 42(15): 10217-10227. [百度学术]
Li B, Kim Y S, Mukundan R, et al. Mixed hydrocarbon/fluoropolymer membrane/ionomer MEAs for durablity studies[J]. ECS Transactions, 2010, 33(1): 913-924. [百度学术]
Crum M, Liu W. Effective testing matrix for studying membrane durability in PEM fuel cells: part 2. mechanical durability and combined mechanical and chemical durability[J]. ECS Transactions, 2006, 3(1): 541-550. [百度学术]
Huang X Y, Solasi R, Zou Y, et al. Mechanical endurance of polymer electrolyte membrane and PEM fuel cell durability[J]. Journal of Polymer Science Part B: Polymer Physics, 2006, 44(16): 2346-2357. [百度学术]
Alizadeh E, Barzegari M M, Momenifar M, et al. Investigation of contact pressure distribution over the active area of PEM fuel cell stack[J]. International Journal of Hydrogen Energy, 2016, 41(4): 3062-3071. [百度学术]
Yang D J, Tan Y L, Li B, et al. A review of the transition region of membrane electrode assembly of proton exchange membrane fuel cells: design, degradation, and mitigation[J]. Membranes, 2022, 12(3): 306. [百度学术]
Ye D H, Zhan Z G. A review on the sealing structures of membrane electrode assembly of proton exchange membrane fuel cells[J]. Journal of Power Sources, 2013, 231: 285-292. [百度学术]
Bograchev D, Gueguen M, Grandidier J C. Stress and plastic deformation of MEA in fuel cells: Stresses generated during cell assembly[J]. Journal of Power Sources, 2008, 180: 393-401. [百度学术]
叶东浩, 李永静, 汪广进, 等. PEM燃料电池膜电极边框材料对膜干湿应力影响研究[J].武汉理工大学学报,2012,34(10):1-4. [百度学术]
Ye D H, Li Y J, Wang G J, et al. Stress effect of membrane electrode frame on the membrane dry-wet cycling in PEM fuel cell[J]. Journal of Wuhan University of Technology, 2012, 34(10): 1-4. (in Chinese) [百度学术]
叶东浩. 膜电极压缩引起的燃料电池内部应力、内阻及水传输行为研究[D]. 武汉: 武汉理工大学, 2014. [百度学术]
Ye D H. Effects of MEA compression on fuel cell stress, resistances and water transmission[D]. Wuhan: Wuhan University of Technology, 2014. (in Chinese) [百度学术]
Yue W, Qiu D, Yi P, et al. Study on the degradation mechanism of the frame for membrane electrode assembly in proton exchange membrane fuel cell[J]. International Journal of Hydrogen Energy, 2021, 46: 36954-36968. [百度学术]
Camanho P P, Davila C G, de Moura M F. Numerical simulation of mixed-mode progressive delamination in composite materials[J]. Journal of Composite Materials, 2003, 37(16): 1415-1438. [百度学术]
Banan R, Bazylak A, Zu J. Effect of mechanical vibrations on damage propagation in polymer electrolyte membrane fuel cells[J]. International Journal of Hydrogen Energy, 2013, 38(34): 14764-14772. [百度学术]