HZSM-5 synthesized by montmorillonite and its application for production of hydrogen via steam reforming of dimethyl ether
-
摘要:
以亚熔岩活化蒙脱土为原料，采用传统水热法、晶种导向法、蒸汽辅助晶化法和气溶胶晶化法制备了具有不同织构和酸性的HZSM-5，并以其为固体酸与商品化的Cu/ZnO/Al2O3物理混合组成双功能催化剂，并考察其二甲醚水蒸气重整(SRD)反应性能。通过XRD、FTIR、SEM、N2低温吸附-脱附和NH3-TPD对样品进行表征。结果表明：合成方法影响了HZSM-5的晶体尺寸、织构及酸性等，进而影响了以其为固体酸的双功能催化剂的SRD性能。传统水热法合成的TH-HZSM-5因其具有适宜的酸量及相对大的表面积、孔容和孔径，以其为固体酸与Cu/ZnO/Al2O3以质量比为1∶1组成的双功能催化剂表现出较好的SRD反应性能，在压力为0.1 MPa、温度为300℃、空速为
5000 mL/(g·h)的反应条件下，二甲醚初始转化率和氢气收率分别达到了100%和91%。Abstract:HZSM-5 with different texture and acidity was prepared using sub-molten salt activated montmorillonite as raw material by traditional hydrothermal method, seed-directed method, steam-assisted crystallization method and aerosol crystallization method. And the HZSM-5 used as solid acid was physically mixed with commercial Cu/ZnO/Al2O3 to obtain bifunctional catalysts for steam reforming of dimethyl ether (SRD) reaction. The samples were systematically characterized by XRD, FTIR, SEM, N2 adsorption-desorption at low temperature, and NH3-TPD techniques. The results showed that the crystal size, textural and acidity properties of HZSM-5 could be effectively regulated by the synthesis method, and then affecting the SRD performance of the corresponding bifunctional catalyst. TH-HZSM-5 synthesized by the traditional hydrothermal method exhibited suitable acid content and relatively large surface area, pore volume and pore size, and thus bifunctional catalyst composed of TH-HZSM-5 as solid acid and Cu/ZnO/Al2O3 with a mass ratio of 1∶1 showed better SRD reaction performance, and the dimethyl ether conversion and H2 yield reached 100% and 91% under the conditions of reaction temperature 300℃, pressure 0.1 MPa, space velocity
5000 mL/(g·h), respectively.-
Keywords:
- dimethyl ether /
- steam reforming /
- hydrogen production /
- HZSM-5 /
- montmorillonite
-
随着全球能源结构的转型和对清洁能源需求的不断增长，氢能作为一种清洁、高效的能源载体，正逐渐成为新能源领域的研究热点[1-3]。然而，氢能的高效储存和运输是实现其广泛应用的关键技术挑战之一。在众多储氢技术中，纤维缠绕复合材料高压储氢技术因技术成熟度和成本优势，被广泛应用于航空航天和交通运输领域[4-7]。复合材料缠绕层作为载荷承担的主体，是决定气瓶结构承载能力和安全性的重要因素。爆破压力是高压储氢容器设计的关键指标之一。为提高复合压力容器的结构效率，降低储氢成本，通常在设计阶段使用有限元方法对缠绕结构的爆破行为进行分析并对缠绕层进行优化。目前储氢气瓶的有限元分析中，直接将气瓶缠绕层看作不同角度纤维层的堆叠，基于层合结构数值模型对纤维缠绕储氢气瓶的结构性能进行评估[8-10]。但纤维束在螺旋缠绕过程中会产生交叉起伏，直接采用层合数值模型难以考虑交叉起伏形态特征对缠绕层力学性能的影响[11-15]，影响爆破压力预测结果的准确性。
已有研究发现纤维束的交叉起伏会降低复合材料的力学性能[16-18]。Shen等[19]建立了一种介观模型来研究纤维交叉起伏对缠绕复合材料刚度的影响，证明了纤维交叉起伏会产生显著影响纤维缠绕复合材料(Filament-wound composites，FWC)的弹性参数。Shen提出，缠绕纤维束厚度也是影响力学性能的一个关键因素。Henry等[20]研究了纤维缠绕圆筒在压缩载荷下损伤行为，发现纤维起伏会降低纤维缠绕圆筒的压缩强度。Chang等[21]研究了缠绕圆筒在拉伸、扭转和多轴载荷下的损伤响应，结果表明纤维损伤首先出现在缠绕交叉起伏区域，并沿此区域扩展。Chang等[22]进一步建立缠绕圆筒全尺寸模型，基于缠绕特征区域赋予材料属性，结果证明这种模型可以模拟缠绕圆筒的损伤过程。Hameed等[23]提出了一种模拟缠绕图案特征的压力容器分区建模方法。与传统的有限元模型相比，该模型可以捕捉缠绕图案对应变分布的影响。Arellano等[24]则采用数字图像相关技术(DIC)对缠绕平板应变分布进行了监测，验证了分区建模方法的合理性。在此基础上，肖磊等[25]使用数值仿真和实验手段相结合的方法进一步对比了FWC平板和标准层合复合材料(Standard laminate composites，SLC)平板的应变分布，结果证明缠绕结构菱形特征图案中部纤维交叉起伏区域存在明显的应变集中现象，是导致该区域损伤起始的主要原因。
在复合材料气瓶结构设计中，有限元分析是常用的手段之一[26]。当前报道的建模类型包括轴对称模型[11]和三维模型[15]两大类。根据Lekhnitskii假设[27]，材料围绕一个轴是对称的，使用轴对称模型不会影响爆破压力结果，并且可以有效提高计算效率[28]。相比轴对称模型，三维模型可以考虑纤维缠绕变角度等特征，开展基于三维尺度的水压爆破和低速冲击分析[29-30]，可提供更高的分析精度。三维壳单元模型通过在壳结构上定义复合材料厚度和方向模拟气瓶结构[31-35]，从而降低气瓶建模的难度。壳单元在建模难度和计算效率上具有优势[36-38]，但三维实体单元[39-42]更能反映缠绕层的几何特征，更适合于精确计算和捕捉复合材料层间的损伤行为。为了研究结构失效机制，许多损伤失效模型被应用在气瓶有限元分析中，包括Hashin失效准则[43]、Puck失效准则[31]、Hashin–Rottem失效准则[35]等。为了提高三维实体模型的计算效率，减少结构模型的尺寸，借助周期性边界条件，建立气瓶的1/4[44]、1/12[8]和1/72[39]模型，是当前常用的有效手段。
在已有的文献报道中，气瓶多基于简化的轴对称模型和三维模型，将缠绕层假设为层合结构进行分析，忽略了纤维缠绕结构中的纤维交叉起伏对纤维方向强度的影响。但对缠绕结构的研究已经表明，纤维交叉起伏会影响缠绕结构的力学性能。因此在复合材料气瓶结构分析中，忽略纤维交叉起伏影响会减低预测精度，难以准确分析缠绕结构的失效模式和强度。虽然通过建立包含纤维束形态的气瓶介观模型，可以直接考虑纤维束交叉起伏特征，但与平板模型和圆筒模型相比，气瓶模型更为复杂，在复合材料气瓶结构分析中使用该方法会显著增加建模难度，并降低计算效率。
为了评估纤维束介观交叉起伏形态的影响，实现气瓶爆破行为的准确预测，并兼顾计算效率，本文提出一种考虑纤维强度折减效应对气瓶层合模型进行修正的爆破失效分析方法(折减分析方法)，并与不考虑纤维强度折减效应的爆破失效分析方法(传统分析方法)进行了对比。首先，采用实验和数值模拟相结合的方法研究纤维交叉起伏特征对缠绕平板拉伸行为的影响机制，探索缠绕参数对纤维拉伸强度的影响规律。然后，通过对传统层合模型分析方法中各角度缠绕层纤维方向拉伸强度进行修正，发展了一种考虑纤维交叉起伏影响的Ⅳ型气瓶爆破行为预测数值分析方法。最后，分别使用折减分析方法和传统分析方法，开展3种不同铺层Ⅳ型气瓶的爆破失效分析，并与对应铺层的气瓶水压爆破试验结果进行对比，验证本文提出的Ⅳ型气瓶爆破失效分析方法的准确性。
1. 纤维强度折减效应研究
在纤维螺旋缠绕过程中，会形成一个以纤维交叉起伏为特征的区域，如图1所示。由纤维交叉引起的纤维束起伏可能会影响复合材料的拉伸强度。纤维强度折减效应研究选择了一个含纤维交叉起伏的特征区域作为研究对象，制备了纤维缠绕的FWC平板试样，用来模拟典型的纤维束交叉起伏特征。同时制备了相同角度的标准层合结构的SLC平板试样作为对照组。沿纤维的+α°方向，对两类试样进行了单向拉伸的实验和仿真的对比分析，获得纤维交叉起伏结构的强度折减规律。
1.1 实 验
首先进行FWC平板和SLC平板试样制备。如图2所示，将预浸料裁成长400 mm、宽6 mm的预浸带，模拟缠绕过程交叉铺贴获得缠绕平板，铺贴过程不考虑张力作用，采用真空袋压法制备了4组缠绕角度为±10°、±20°、±30°和±40°的FWC平板试样。作为对比，制备了4组相同角度的SLC平板试样。试样整体厚度为0.43 mm，共两层，其中试样长度方向沿着+α°纤维束方向，实际平板试样的纤维束夹角为2α°。
![]()
图 2 固化成型和拉伸试样FWC—Filament-wound compositesFigure 2. Curing process and tensile specimens在电子万能试验机(LD23.104，力试科学仪器有限公司)上对试样进行拉伸试验，参考标准ASTM D3039[45]进行加载，并使用DIC设备(Vic-2D，Correlated Solutions Europe)监测加载过程中试样表面的位移场和应变场，如图3所示。为确保DIC方法监测应变场的准确性，在试验前先对试样表面进行散斑喷涂。然后进行试样和设备安装，在加载过程中设置拍摄频率为1 Hz，拉伸位移速度为2 mm/min。拉伸载荷通过试验机传感器获得，实际位移通过DIC的非接触式引伸计获得。
![]()
图 3 实验装置和加载方式DIC—Digital image correlation technologyFigure 3. Experimental equipment and loading apparatus1.2 有限元分析
基于ABAQUS/Explicit求解器建立了三维有限元模型，采用Hashin失效起始准则和渐进损伤演化方式模拟FWC平板和SLC平板试样的拉伸失效行为。为了在模型中准确建立纤维交叉起伏的细节特征，使用光学显微镜拍摄并测量了不同缠绕角度下FWC平板试样截面参数，如图4(a)所示。将测量的起伏角度作为评估缠绕结构起伏特征的关键参数。同时为了建立双层厚度的纤维缠绕结构平板模型，进行缠绕层厚度对纤维折减系数影响的研究，进一步制备了双层厚度的缠绕结构试样，并测量了这些试样交叉区域的起伏角度，如图4(b)所示。图4(c)、图4(d)为单层和双层FWC试样的纤维起伏角度随缠绕角度的变化规律，单层厚度试样平均起伏角度分别为7.6°、7.7°、8.1°和8.6°，双层厚度试样的平均起伏角度分别为9.7°、10.6°、11.2°和12.5°。
在Abaqus有限元软件中建立了FWC平板和SLC平板的有限元模型。SLC平板的有限元模型和边界条件设置如图5所示，在左侧加载区域的上下表面施加固定约束，在右侧加载区域施加位移加载，加载速率与试验一致。样品的拉伸区域长度为100 mm，为了避免边界处的应力集中导致样品过早失效，在拉伸区域的左侧和右侧各7 mm长度区域内未设置材料失效，而在拉伸区域中间的86 mm长度区域内考虑了材料失效。
![]()
图 4 起伏角度(α)随缠绕角度变化：(a)单层纤维的起伏形态；(b) 2层纤维的起伏形态；(c)单层纤维变化规律；(d) 2层纤维变化规律Figure 4. Variations of undulation angle (α) with the winding angle: (a) Undulation feature of single layer; (b) Undulation feature of 2 layers; (c) Variation curve of single layer; (d) Variation curve of 2 layers![]()
图 5 标准层合复合材料(SLC)平板模型和加载条件Figure 5. Model of standard laminate composites (SLC) plate and its loading conditions根据显微镜观察到的FWC平板试样截面交叉起伏几何特征，建立了FWC平板试样的细节模型，如图6所示。模型中的纤维起伏角度与测量的纤维起伏平均角度相同，在交叉起伏边缘和纤维束之间的间隙存在，并使用树脂对这些间隙进行了填充。图7为建立的FWC平板的有限元模型，边界条件和加载方式与试验保持一致。在FWC平板左侧宽50 mm区域上下表面施加固定约束，在FWC平板右侧宽50 mm区域上下表面设置向右的拉伸位移载荷，加载速率与试验一致。为避免边界应力集中导致样品过早失效，试样中间宽86 mm区域的复合材料设置了材料失效，两侧7 mm区域未设置材料失效。
![]()
图 6 FWC平板中纤维交叉起伏形态：(a)实验观测截面图；(b)建模示意图Figure 6. Undulation of fiber bundle in FWC plate: (a) Experimental image; (b) Modeling image单向纤维束材料参数的弹性模量、泊松比和极限强度等参数通过试验测试获得，试验测试参考ASTM测试标准[45-49]，详见表1。
![]()
图 7 不同缠绕角度FWC平板有限元模型：(a) ±10°；(b) ±20°；(c) ±30°；(d) ±40°Figure 7. Finite element (FE) models of FWC plates with different winding angles: (a) ±10°; (b) ±20°; (c) ±30°; (d) ±40°表 1 单向纤维束材料参数Table 1. Material properties of unidirectional fiber bundleItems Value Longitudinal modulus, E11/GPa 125.4 Transverse modulus, E22 =E33/GPa 7.7 In-plane shear modulus, G12=G13/GPa 3.8 Out-of-plane shear modulus, G23/GPa 4.8 Major Poisson's ratio, μ12 = μ13 0.33 Through-thickness Poisson's ratio, υ23 0.35 Longitudinal tensile strength, XT/GPa 2.18 Longitudinal compressive strength, XC/GPa 1.2 Transverse tensile strength, YT/MPa 60 Transverse compressive strength, YC/MPa 140 Density of laminate, ρ/(kg·m−3) 1600 Tensile fracture energy of fiber, Gft/(N·mm−1) 133 Compressive fracture energy of fiber, Gfc/(N·mm−1) 40 Tensile fracture energy of matrix, Gmt/(N·mm−1) 0.6 Compressive fracture energy of matrix,Gmc/(N·mm−1) 2.1 Elastic modulus of resin, E/GPa 3.0 Density of resin, ρr/(kg·m−3) 1200 Poisson's ratio of resin, μ 0.3 采用Hashin失效准则对纤维和基体的起始损伤进行判定，采用基于能量的线性本构描述材料损伤起始后的演化行为。平板拉伸失效分析和气瓶的爆破分析采用相同的失效准则，如下式：
纤维拉伸失效(σ11 > 0)：
(σ11XT)2⩾ (1) 纤维压缩失效(σ11 < 0)：
{\left(\frac{{\sigma }_{11}}{{X}_{\mathrm{C}}}\right)}^{2}\geqslant 1 (2) 基体拉伸失效(σ22 > 0)：
{\left(\frac{{\sigma }_{22}}{{X}_{\mathrm{T}}}\right)}^{2}+{\left(\frac{{\sigma }_{12}}{{S}_{12}}\right)}^{2}+{\left(\frac{({\sigma }_{23}{)}^{2}+{\sigma }_{22}{\sigma }_{23}}{{S}_{23}}\right)}^{2}\geqslant 1 (3) 基体压缩失效(σ22 < 0)：
\left(\frac{\sigma_{22}}{S_{22}}\right)^2+\left(\frac{\sigma_{22}+\sigma_{33}}{Y_{\mathrm{C}}}\right)\left(\left(\frac{Y_{\mathrm{C}}}{2S_{23}}\right)^2-1\right)+ \quad \frac{({\sigma }_{23}^{2}-{\sigma }_{22}{\sigma }_{23})}{{S}_{23}^{2}}+{\left(\frac{{\sigma }_{12}}{{S}_{12}}\right)}^{2}+{\left(\frac{{\sigma }_{13}}{{S}_{13}}\right)}^{2}\geqslant 1 (4) 面内损伤演化：
{d}_{i}=\mathrm{m}\mathrm{a}\mathrm{x}\left\{0,\mathrm{m}\mathrm{i}\mathrm{n}\left\{1,{\varepsilon }_{\mathrm{e}\mathrm{q}}^{\mathrm{f}}\frac{{\varepsilon }_{\mathrm{e}\mathrm{q}}-{\varepsilon }_{\mathrm{e}\mathrm{q}}^{0}}{{\varepsilon }_{\mathrm{e}\mathrm{q}}({\mathrm{\varepsilon }}_{\mathrm{e}\mathrm{q}}^{\mathrm{f}}-{\varepsilon }_{\mathrm{e}\mathrm{q}}^{0})}\right\}\right\} (5) 其中：XT、XC分别为纤维方向拉伸和压缩强度；YC为垂直于纤维方向压缩强度；S22为垂直于纤维方向拉伸强度；S12、S13、S23为剪切强度；σ11、σ22分别为纤维方向和垂直于纤维方向的正应力；σ12、σ13和σ23为剪切应力； {\varepsilon }_{\mathrm{e}\mathrm{q}}^{0} 为复合材料结构失效起始状态下的等效应变； {\varepsilon }_{\mathrm{e}\mathrm{q}}^{\mathrm{f}} 为复合材料结构完全失效状态下的等效应变； {\varepsilon }_{\mathrm{e}\mathrm{q}} 为复合材料结构当前状态下的等效应变[50]。
1.3 结果对比与分析
图8为0.43 mm厚度4组不同缠绕角的FWC平板和SLC平板拉伸实验的位移-载荷曲线。在加载过程中加强片和试件之间会发生滑动，为了使位移值更接近实际值，位移值使用DIC系统的电子引伸计功能测量，载荷值采用试验机记录的加载力。从图中发现，FWC平板的极限强度低于相同角度的SLC平板试件，FWC平板的承载能力显著降低。与相同角度的SLC平板对比，FWC平板刚度分别降低了3.70%、2.93%、1.50%和0.55%，均小于5.00%。FWC平板刚度降低并不明显，因此本文重点讨论纤维强度折减效应。
图9和图10分别为实验和有限元分析中得到的FWC平板和SLC平板的长度方向应变场（εyy）。两者结果相比，FWC平板的实验和模拟应变场都表现出应变集中。相反，SLC平板的应变场没有显示出应变集中，尽管位移增加，SLC平板的表面应变保持着相对均匀的分布。有限元与实验的表面应变结果具有良好的一致性。
![]()
图 8 FWC平板和SLC平板的位移-载荷试验曲线：(a) ±10°；(b) ±20°；(c) ±30°；(d) ±40°Figure 8. Experimental displacement-load curves of FWC and SLC plate specimens: (a) ±10°; (b) ±20°; (c) ±30°; (d) ±40°![]()
图 9 数字图像相关技术(DIC)监测FWC平板和SLC平板长度方向应变(εyy)结果：(a) ±10°；(b) ±20°；(c) ±30°；(d) ±40°Figure 9. Distribution of vertical strain (εyy) of FWC and SLC plates monitored by DIC: (a) ±10°; (b) ±20°; (c) ±30°; (d) ±40°![]()
图 10 FWC平板和SLC平板长度方向仿真应变结果：(a) ±10°；(b) ±20°；(c) ±30°；(d) ±40°Figure 10. Numerical predicted distribution of vertical strain fields of FWC and SLC plates: (a) ±10°; (b) ±20°; (c) ±30°; (d) ±40°图11为两种平板各角度仿真(FEA)和实验(EXP)的位移-载荷曲线及FWC平板的纤维拉伸失效云图。结果表明，仿真与实验所得的平板拉伸刚度高度一致，且在拉伸极限强度方面的误差较小，这验证了数值分析与实验结果之间具有很好的一致性。右侧云图为FWC平板纤维拉伸失效过程仿真结果。图11中载荷达到 a点时，FWC平板少数单元应力水平达到纤维拉伸起始失效判据，单元开始失效，单元刚度出现退化；随着继续载荷达到b点，失效面积逐渐扩展；最终加载至c点，缠绕交叉区域多数单元完全失效，FWC平板失去承载能力。
![]()
图 11 仿真和实验位移-载荷曲线和FWC平板纤维拉伸失效演化：(a) ±10°；(b) ±20°；(c) ±30°；(d) ±40°dft—Variable of fiber tensile failure state; FEA—Simulation; EXP—ExperimentFigure 11. Simulation and experimental displacement-load curves and the evolution of fiber tensile failure in FWC plate: (a) ±10°; (b) ±20°; (c) ±30°; (d) ±40°图12为
4000 N载荷下FWC平板+α°层拉伸方向应变云图和纤维方向应力云图。可知，当载荷均为4000 N时，FWC平板在纤维交叉起伏区域均出现了明显的应变集中和应力集中现象，且随着缠绕角度的增大，应变集中和应力集中效果越明显，拉伸方向最大应变从0.0127 增加到0.0143 ，纤维方向最大应力从1713 MPa增加到1986 MPa。同时，应变集中和应力集中的区域与纤维起伏的区域重合。据此可知，+α°层的起伏是导致应变集中和应力集中产生的原因，使纤维过早达到其极限拉伸强度，进而导致FWC平板的承载能力下降。缠绕角度增加，纤维起伏角度越大，承载能力下降效果越明显，纤维强度折减效应越显著。由此分析，当相同角度缠绕层厚度增加时，纤维起伏角度也会增加，并使承载能力进一步下降，纤维强度折减效应更加显著。![]()
图 124000 N载荷下FWC平板+α°层云图：(a)拉伸方向应变；(b)纤维方向应力LE—Logarithmic strain; LE11—LE in the tensile direction; S11—Stress in the fiber direction; S—StressFigure 12. FWC plate +α° layer nephogram under4000 N: (a) Strain in the tensile direction; (b) Stress in the fiber direction1.4 纤维强度折减系数
图13为有限元分析两种厚度试样的拉伸失效载荷结果。可以发现，所有FWC平板的拉伸失效载荷都低于对应的SLC平板。定义FWC平板与SLC平板的最大失效载荷之比为纤维强度折减系数，则纤维强度折减系数的变化规律如图14所示。单层板的不同缠绕角度的强度折减系数分别为0.93、0.91、0.90和0.89。双层平板不同缠绕角度的纤维强度折减系数分别为0.80、0.78、0.70和0.68。随着缠绕角度变大，强度折减系数逐渐降低。同时缠绕层厚度的增加，相同缠绕角度的纤维强度折减系数也显著降低。
![]()
图 13 有限元分析预测的FWC平板和SLC平板拉伸失效载荷对比：(a)单层厚度；(b)双层厚度Figure 13. Comparison of tensile failure loads of FWC and SLC plates predicted by finite element analysis (FEA): (a) 1-layer thickness; (b) 2-layers thickness![]()
图 14 纤维强度折减系数随厚度和缠绕角度变化Figure 14. Variation of strength reduction factors with winding angle and thickness根据图14中的纤维强度折减系数随缠绕角度和厚度的变化规律，建立纤维强度折减系数随厚度和角度的线性拟合公式，如下式所示：
\kappa =\frac{({t}_{\mathrm{f}}-{t}_{2})}{({t}_{1}-{t}_{2})}({\kappa }_{1}-{\kappa }_{2})+{\kappa }_{2} (6) 式中：κ表示纤维强度折减系数；t1和t2表示本文中使用的单层和双层预浸料的厚度；tf表示缠绕单层的厚度； {\kappa }_{1} 和 {\kappa }_{2} 的表示式分别为： {\kappa }_{1} =0.94−
0.0013 α和 {\kappa }_{2} =0.85−0.0044 α；α表示当前的缠绕角度。当缠绕制品的纤维束带宽为6 mm，缠绕螺旋层厚度在0.4~0.8 mm，缠绕角度范围在±10°~±40°之间，可以使用式(6)对纤维拉伸强度进行折减分析。本文的气瓶数值分析中，使用该经验公式对气瓶不同缠绕层的拉伸强度进行折减，开展进一步的气瓶爆破失效分析，并与未考虑强度折减的分析方法进行了对比。
2. 基于纤维强度折减效应的气瓶数值分析
基于本文提出的纤维强度折减效应，对3种不同铺层的IV型气瓶进行爆破行为预测，分别使用不考虑纤维强度折减效应的传统分析方法和考虑纤维强度折减效应的折减分析方法对气瓶爆破压力和爆破失效位置进行数值分析预测。
图15为气瓶数值分析框架。在气瓶数值分析中，首先建立包含变厚度分层信息的气瓶三维模型，并对气瓶划分网格。然后计算缠绕层不同位置的变角度信息后，将角度信息赋予在模型网格上。随后将考虑纤维强度折减效应的VUMAT材料参数赋予在气瓶缠绕层几何模型上，并赋予内胆对应的材料属性，继续设置模型的接触属性和边界条件，最后提交计算进行气瓶的渐进损伤有限元分析。
![]()
图 15 基于纤维强度折减效应的气瓶数值分析框架Figure 15. Numerical analysis framework for cylinders based on the effect of strength reduction effect2.1 铺层信息和材料参数
使用图16所示的9 L-IV型气瓶内胆进行铺层设计。筒身为高密度聚乙烯(HDPE)材料，两端金属封头(BOSS)为6061-T6铝合金材料。筒身长度400 mm，外径164 mm，厚度5 mm，采用椭圆形封头，长轴为82 mm，短轴为53 mm，极孔半径为23 mm。为了验证本文提出的折减分析方法在不同爆破失效模式下的预测精度和适用性，设计了3种不同铺层，铺层信息如下：
A：[±14°/±90°/±90°/±25°/±90°/±35°/±90°/±14°]；
B：[±14°/±90°/±90°/±25°/±90°/±35°/±90°/±90°]；
C：[±14°/±90°/±25°/±90°/±35°/±35°/±90°/±14°]。
![]()
图 16 9 L-IV型内胆尺寸HDPE—High density polyethyleneFigure 16. Dimensions of the 9 L-IV type liner螺旋层纤维方向拉伸强度XT取折减后对应强度值，一个循环的螺旋层厚度为0.5 mm，根据对纤维强度折减系数的经验公式计算，最终±14°、±25°及±35°螺旋缠绕层的纤维强度折减系数分别为0.90、0.88和0.86，对应折减后的螺旋层纤维拉伸强度分别为
2250 MPa、2200 MPa和2150 MPa。本文所用IV型气瓶材料中的复合材料缠绕层[51]、两端金属瓶口(BOSS)[52]和筒身高密度聚乙烯(HDPE)内胆[53]的材料属性如表2所示。表 2 IV型气瓶材料力学性能参数Table 2. Type IV cylinder material mechanical parametersItem Value Longitudinal modulus, E11/GPa 154 Transverse modulus, E22 =E33/GPa 11.4 In-plane shear modulus, G12=G13/GPa 4.8 Out-of-plane shear modulus, G23/GPa 3.8 Major Poisson's ratio, μ12 = μ13 0.3 Through-thickness Poisson's ratio, υ23 0.33 Longitudinal tensile strength, XT/GPa 2.5 Longitudinal compressive strength, XC/GPa 1.2 Transverse tensile strength, YT/MPa 70 Transverse compressive strength, YC/MPa 180 Density of laminate, ρ/(kg·m−3) 1600 Tensile fracture energy of fiber, Gft/(N·mm−1) 133 Compressive fracture energy of fiber, Gfc/(N·mm−1) 40 Tensile fracture energy of matrix, Gmt/(N·mm−1) 0.6 Compressive fracture energy of matrix,Gmc/(N·mm−1) 2.1 Elastic modulus of HDPE, E/GPa 1.1 Poisson's ratio of HDPE, μ 0.38 Yield strength of HDPE, σs/MPa 22.9 Ultimate strength of HDPE, σb/MPa 25 Fracture elongation of HDPE, δ/% >600 Elastic modulus of BOSS, E/GPa 69 Poisson's ratio of BOSS, μ 0.324 Yield strength of BOSS, σs/MPa 298 Ultimate strength of BOSS, σb/MPa 330 Fracture elongation of BOSS, δ/% 12 Note: BOSS—Bolted opening support structure. 2.2 模型与边界条件
建立气瓶的1/360模型，单元类型为C3D8R，模型共计
26318 个单元，如图17所示。在气瓶剖面施加周期性边界条件，根据气瓶爆破试验工况建立边界条件：在气瓶左侧接头处设置为沿轴向的固定约束，避免发生刚体位移；内胆的塑料部分和金属BOSS部分的接触属性设置为Tie，模拟两者之间的固定效果；内胆外表面与缠绕层内表面的接触属性设置为Tie，模拟两者之间的粘接效果。气瓶内表面施加100 MPa压力(P)，基于Hashin失效准则进行分析，基于能量的线性本构描述材料损伤起始后的演化行为，实际爆破压力以缠绕层纤维贯穿损伤为判断依据。![]()
图 17 IV型气瓶模型：(a)整体模型和载荷；(b)封头；(c)过渡区域P—PressureFigure 17. Type IV cylinder model: (a) Global model and loading; (b) Dome; (c) Transition region3. 气瓶制备与水压测试
为了验证两种预测方法的准确性，根据3种缠绕线型制备了对应的IV型气瓶，并进行了水压爆破测试。
制备过程如图18所示，采用湿法缠绕工艺，对IV型气瓶进行缠绕制备。内胆主体材料为高密度聚乙烯(High-density polyethylene，HDPE)，BOSS材料为6061-T6铝合金。碳纤维缠绕层材料为光威TZ700-24 K级碳纤维增强树脂基复合材料，树脂为博汇EpoTech®425型环氧树脂。气瓶内胆使用滚塑工艺制成，通过滚塑模具将铝合金BOSS和聚乙烯内胆滚塑一体成型。根据湿法缠绕工艺要求，缠绕前先将树脂和固化剂按比例混合，倒入到胶槽中，纤维束经过胶槽滚轮浸润树脂后，连接到丝嘴位置，然后缠绕到瓶身位置。缠绕路径由一台六轴五联动数控缠绕机(SLW01.6-500-4/1-5000CNC，湖南江南四棱数控有限公司)控制，通过 CADWIND软件编写的环向缠绕和螺旋缠绕程序，将纤维束按照缠绕层顺序缠绕在气瓶表面。缠绕完成后，将气瓶安装到旋转固化炉中，按照树脂固化曲线设置温度，经固化成型后，得到碳纤维增强树脂基体(Carbon fiber reinforced polymer，CFRP)复合材料储氢气瓶试样。
如图19(a)所示，根据标准GB/T 15385—2022《气瓶水压爆破试验方法》[54]，使用水压爆破试验机(EHM-8102，深圳市恩普达工业系统有限公司)对IV型气瓶进行爆破测试。爆破前首先将在气瓶内部灌满水并排除内部空气，防止在爆破瞬间内部空气造成的碎片飞溅，损坏试验设备或者危害人身安全。然后使用高压水管连接爆破试验机和气瓶，水管和气瓶之间使用转接头连接，用生胶带和橡胶圈保证密封性。将气瓶放置在钢筒中，防止爆破碎屑飞出。最后启动试验机，加压开始爆破测试。测试装置如图19(a)所示。
![]()
图 18 IV型气瓶制造CFRP—Carbon fiber reinforcedpolymer; PE—PolyethyleneFigure 18. Type IV cylinder manufacturing![]()
图 19 气瓶水压实验及爆破失效形式：(a)试验装置；(b) A铺层；(c) B铺层；(d) C铺层Figure 19. Pressure test and burst modes of type IV cylinders: (a) Burst test setup; (b) Layup-A; (c) Layup-B; (d) Layup-C根据图19(b)~图19(d)爆破试验结果，A铺层气瓶的爆破压力为62.86 MPa，爆破位置在封头与筒身过渡的肩部，B铺层气瓶的爆破压力为49.21 MPa，爆破位置在BOSS区域，C铺层气瓶的爆破压力为52.95 MPa，爆破位置在筒身区域。
4. 结果与讨论
图20为气瓶爆破失效前纤维应力云图。图20(a)~图20(c)为传统分析方法的结果，图20(d)~图20(f)为折减分析方法的结果。对传统分析方法的结果进行分析，根据图20(a)可知，A气瓶在爆破失效前纤维方向最大应力出现在筒身段环向层位置，当内部压力加载至66.4 MPa，环向层纤维方向应力值达到
2493.36 MPa，接近环向层复合材料纤维拉伸强度。图20(b)中，B气瓶在爆破失效前纤维方向最大应力出现在靠近两端金属BOSS附近的最内层14°螺旋层，56.6 MPa压力下该位置纤维方向应力为2486.25 MPa，接近未折减14°螺旋层复合材料纤维拉伸强度。图20(c)中，C气瓶在爆破失效前纤维方向最大应力在筒身段环向层位置，压力达到53.0 MPa时，该位置纤维方向应力为2492.30 MPa，接近环向层复合材料纤维拉伸强度。![]()
图 20 气瓶爆破失效前纤维应力云图：(a) A铺层-传统分析方法；(b) B铺层-传统分析方法；(c) C铺层-传统分析方法；(d) A铺层-折减分析方法；(e) B铺层-折减分析方法；(f) C铺层-折减分析方法Figure 20. Fiber stress nephogram of cylinder before failure: (a) Layup A-traditional method; (b) Layup B-traditional method; (c) Layup C-traditional method; (d) Layup A-reduction modified method; (e) Layup B-reduction modified method; (f) Layup C-reduction modified method根据图20(d)~图20(f)，对折减分析方法的结果进行分析。A气瓶在爆破失效前最大应力同样出现在筒身中部最内层的环向缠绕层，当内压达到57.6 MPa时，最大应力为
2156.69 MPa，未达到环向层复合材料纤维拉伸强度。对于B气瓶，在爆破失效前，压力达到52.0 MPa时，缠绕层纤维方向最大应力出现在金属BOSS附近最内层14°螺旋层，最大应力为2236.68 MPa，接近折减后14°螺旋层复合材料纤维拉伸强度。对于C气瓶，在爆破失效前，最大应力出现在筒身中部最内层的环向缠绕层，内压达到53.0 MPa，纤维方向最大应力为2492.30 MPa，接近环向层复合材料纤维拉伸强度。图21为气瓶爆破失效时纤维损伤云图，将纤维损伤贯穿缠绕层作为判断气瓶爆破失效的依据。其中图21(a)~图21(c)为传统分析方法的结果，图21(d)~图21(f)为折减分析方法的结果。对A气瓶进行分析，根据图21(a)、图21(d)结果发现，使用传统分析方法预测的爆破压力为68 MPa，爆破发生在筒身位置。使用折减分析方法预测A气瓶爆破压力为59.6 MPa，气瓶爆破发生在筒肩位置。对B气瓶进行分析，根据图21(b)、图21(e)结果发现，传统分析方法爆破压力为56.8 MPa，爆破发生在BOSS位置。使用折减分析方法，B气瓶爆破压力为52.2 MPa，爆破发生在BOSS位置。对C气瓶进行分析，根据图21(c)和图21(f)结果发现，两种方法的爆破压力都是56.4 MPa，爆破位置都位于筒身位置。和图19中水压爆破试验的结果对比，试验中A气瓶破坏位置在筒肩位置，和折减分析方法破坏位置相同，传统分析方法则破坏于筒身位置。试验中B气瓶破坏发生在BOSS位置，与两种分析方法预测的爆破位置相同。试验中C气瓶破坏发生在筒身位置，与两种分析方法预测的爆破位置相同。
表3给出了3种铺层气瓶的实验与仿真爆破压力和爆破失效位置的结果对比。从表中可以发现，对A气瓶，折减分析方法得到的爆破压力误差为5.19%，爆破位置预测正确；传统分析方法得到的爆破压力误差为8.18%，爆破位置预测错误。对B气瓶，折减分析方法得到的爆破压力误差为6.07%，爆破位置预测正确；传统分析方法得到的爆破压力误差为15.4%，爆破位置预测正确。对C气瓶，两种方法预测的结果相同，爆破压力误差均为6.52%，爆破位置预测正确。
![]()
图 21 气瓶爆破失效纤维拉伸损伤云图对比：(a) A铺层-传统分析方法；(b) B铺层-传统分析方法；(c) C铺层-传统分析方法；(d) A铺层-折减分析方法；(e) B铺层-折减分析方法；(f) C铺层-折减分析方法Figure 21. Comparison of fiber tensile failure nephogram of cylinder: (a) Layup A-traditional method; (b) Layup B-traditional method; (c) Layup C-traditional method; (d) Layup A-reduction modified method; (e) Layup B-reduction modified method; (f) Layup C-reduction modified method表 3 气瓶爆破失效结果对比Table 3. Comparison of burst failure results of cylindersNumber Method Pressure/MPa Burst location Error/% A Test 62.86 Transition region − Traditional method 68.00 Cylinder body +8.18 Reduction modified method 59.60 Transition region −5.19 B Test 49.21 BOSS − Traditional method 56.80 BOSS +15.42 Reduction modified method 52.20 BOSS +6.07 C Test 52.95 Cylinder body − Traditional method 56.40 Cylinder body +6.52 Reduction modified method 56.40 Cylinder body +6.52 综合以上结果分析，对于A气瓶，使用传统分析方法预测时，其爆破预测位置在筒身段，各角度缠绕层强度相同，最大应力出现在筒身段。而考虑纤维强度折减效应后，在环向层应力未达到强度极限之前，筒肩部位的螺旋层就已经达到了经折减后的纤维拉伸强度，导致预测的爆破位置发生改变，爆破压力数值降低，更接近实验测试值。对于B气瓶，由于相对于A气瓶减少了一个±14°螺旋层，导致气瓶BOSS部位附近的螺旋层非常薄弱，成为爆破失效起始位置，增加的±90°环向层导致筒身冗余设计，因此气瓶爆破压力较小。对于两种分析方法，其预测的失效部位一致，均为BOSS区域，由于折减分析方法的螺旋层纤维强度进行了折减，因此其预测的爆破压力更低，误差更小。对于C气瓶，由于相对于A气瓶减少了一个±90°环向层，因此环向层的承载能力降低，而封头部位承载能力较强，气瓶在筒身段发生爆破。在两种分析方法中，其筒身段环向层的纤维拉伸强度相同，因此两种方法的预测结果一致。因此，本文提出的考虑纤维强度折减效应的IV型气瓶爆破失效分析，能更精准预测螺旋层主导的气瓶爆破压力和爆破形式，对于IV型气瓶的轻量化设计具有重要的指导意义。
5. 结 论
以提高气瓶分析精度为目的，本文进行了纤维强度折减效应研究，并发展了考虑纤维强度折减效应的IV型气瓶折减分析方法。得到以下结论：
(1)基于纤维缠绕结构平板和标准层合结构平板拉伸的实验和仿真对比，发现了缠绕结构纤维方向强度折减规律。随着缠绕角度和缠绕层厚度的增加，折减系数减小，纤维方向承载能力降低，使用线性拟合获得本文气瓶缠绕层纤维强度随缠绕角度和厚度变化的折减经验公式；
(2)开展了3种铺层气瓶的水压爆破试验，并与两种预测方法的预测结果进行对比。结果发现，折减分析方法较传统分析方法具有更高的爆破压力预测精度，对A铺层和B铺层气瓶的爆破压力预测误差分别从+8.18%和+15.42%降低至5.19%和+6.07%；
(3) 3种铺层的爆破失效位置预测结果中，折减分析方法较传统分析方法具有更高的预测准确度。采用纤维强度折减系数对Ⅳ型气瓶螺旋层纤维拉伸强度进行修正，可有效提升气瓶爆破失效行为分析结果的合理性。
-
![]()
图 2 不同方法合成HZSM-5的XRD图谱
Figure 2. XRD patterns of HZSM-5 synthesized by different methods
TH-HZSM-5, SD-HZSM-5, SAC-HZSM-5, AC-HZSM-5—HZSM-5 synthesized by traditional hydrothermal, seed-directed, steam-assisted crystallization and aerosol crystallization method, respectively
![]()
图 3 不同方法合成HZSM-5的FTIR图谱
Figure 3. FTIR spectra of HZSM-5 synthesized by different methods
![]()
图 4 不同方法合成HZSM-5的SEM图像
Figure 4. SEM images of HZSM-5 synthesized by different methods
![]()
图 5 不同方法合成HZSM-5的N2吸附等温线
Figure 5. N2 adsorption isotherm of HZSM-5 synthesized by different methods
![]()
图 6 不同方法合成HZSM-5分子筛的NH3-TPD图
Figure 6. NH3-TPD spectra of HZSM-5 synthesized by different methods
![]()
图 7 不同方法合成的HZSM-5+Cu/ZnO/Al2O3双功能催化剂的二甲醚(DME)转化率(a)和H2收率(b)
Figure 7. Dimethyl ether (DME) conversion (a) and H2 yield (b) of bifunctional catalysts composed of HZSM-5 synthesized by different methods and Cu/ZnO/Al2O3
![]()
图 8 不同方法合成的HZSM-5+Cu/ZnO/Al2O3双功能催化剂的含碳产物选择性
Figure 8. Selectivity of carbon containing products of bifunctional catalyst composed of HZSM-5 synthesized by different methods and Cu/ZnO/Al2O3
![]()
图 9 不同质量比的TH-HZSM-5+Cu/ZnO/Al2O3双功能催化剂的DME转化率(a)和H2收率(b)
Figure 9. DME conversion (a) and H2 yield (b) of bifunctional catalysts composed of TH-HZSM-5 and Cu/ZnO/Al2O3 with different mass ratios
![]()
图 10 不同质量比的HZSM-5+Cu/ZnO/Al2O3双功能催化剂的含碳产物选择性
Figure 10. Selectivity of carbon containing products of bifunctional catalyst composed of HZSM-5 and Cu/ZnO/Al2O3 with different mass ratios
表 1 不同方法合成的HZSM-5分子筛的织构特征
Table 1 Summary of the textural properties of HZSM-5 synthesized by different methods
Sample SBET/(m2·g−1) Vtotal/(cm3·g−1) dave/nm Smicro/(m2·g−1) Vmicro/(cm3·g−1) TH-HZSM-5 399 0.30 3.03 315 0.15 SD-HZSM-5 380 0.26 2.77 288 0.13 SAC-HZSM-5 416 0.30 2.86 335 0.16 AC-HZSM-5 422 0.33 3.11 371 0.17 Notes: SBET—Brunner-Emmet-Teller (BET) specific surface area; Vtotal—Total pore volume; dave—Average pore size; Smicro—Microporous surface area; Vmicro—Microporous volume. 
下载: 导出CSV
表 2 不同方法合成的HZSM-5的组成及硅铝比
Table 2 Composition and molar ratio of SiO2/Al2O3 of HZSM-5 synthesized by different methods
Sample SiO2/wt% Al2O3/wt% Na2O/wt% Other impurity/wt% SiO2/Al2O3 molar ratio TH-HZSM-5 89.86 3.63 0.66 5.85 42 SD-HZSM-5 93.20 3.43 0.92 2.45 46 SAC-HZSM-5 92.30 3.42 0.74 3.54 46 AC-HZSM-5 92.49 3.40 0.82 3.29 46
下载: 导出CSV
表 3 不同方法合成HZSM-5的酸量
Table 3 Amount of acid site of HZSM-5 synthesized by different methods
Sample Amount of desorbed NH3/(mmol·g−1) Weak acid Strong acid Total acid TH-HZSM-5 0.13 0.09 0.22 SD-HZSM-5 0.08 0.04 0.12 SAC-HZSM-5 0.22 0.18 0.40 AC-HZSM-5 0.24 0.20 0.44
下载: 导出CSV
-
[1] LI J, WU T, CHENG C, et al. A review of the research progress and application of key components in the hydrogen fuel cell system[J]. Processes, 2024, 12(2): 249-277. DOI: 10.3390/pr12020249
[2] YANG M, MEN Y, LI S L, et al. Enhancement of catalytic activity over TiO2-modifed Al2O3 and ZnO-Cr2O3 composite catalyst for hydrogen production via dimethyl ether steam reforming[J]. Applied Catalysis A: General, 2012, 433-434: 26-34. DOI: 10.1016/j.apcata.2012.04.032
[3] 高天宇, 赵永华, 郑择, 等. 酸活化蒙脱土在二甲醚水蒸气重整制氢中的应用[J]. 燃料化学学报, 2021, 49(10): 1498-1503. GAO Tianyu, ZHAO Yonghua, ZHENG Ze, et al. Acid activation of montmorillonite and its application for production of hydrogen via steam reforming of dimethyl ether[J]. Journal of Fuel Chemistry and Technology, 2021, 49(10): 1498-1503 (in Chinese).
[4] WEI S, LI C, REN H, et al. Design and optimization of hydrogen production system model for dimethyl ether self-heating steam reforming[J]. International Journal of Hydrogen Energy, 2024, 49: 1450-1467. DOI: 10.1016/j.ijhydene.2023.08.369
[5] YANG W W, MA X, TANG X Y, et al. Review on developments of catalytic system for methanol steam reforming from the perspective of energy-mass conversion[J]. Fuel, 2023, 345: 128234. DOI: 10.1016/j.fuel.2023.128234
[6] LI J, ZHANG Q J, LONG X, et al. Hydrogen production for fuel cells via steam reforming of dimethyl ether over commercial Cu/ZnO/Al2O3 and zeolite[J]. Chemical Engineering Journal, 2012, 187: 299-305. DOI: 10.1016/j.cej.2012.01.126
[7] GAO T Y, ZHAO Y H, ZHANG Q J, et al. Zinc oxide modified HZSM-5 as an efficient acidic catalyst for hydrogen production by steam reforming of dimethyl ether[J]. Reaction Kinetics, Mechanisms and Catalysis, 2019, 128: 235-249. DOI: 10.1007/s11144-019-01642-5
[8] LONG X, SONG Y H, LIU Z T, et al. Insights into the long-term stability of the magnesia modified H-ZSM-5 as an efficient solid acid for steam reforming of dimethyl ether[J]. International Journal of Hydrogen Energy, 2019, 44: 21481-21494. DOI: 10.1016/j.ijhydene.2019.06.177
[9] NISHIGUCHI T, OKA K, MATSUMOTO T, et al. Durability of WO3/ZrO2-CuO/CeO2 catalysts for steam reforming of dimethyl ether[J]. Applied Catalysis A: General, 2006, 301: 66-74. DOI: 10.1016/j.apcata.2005.11.011
[10] LUO H X, ZHAO Y H, ZHANG Q J, et al. The role of promoters in Cu/acid-MMT catalysts for production of hydrogen via steam reforming of dimethyl ether[J]. Journal of Chemical Technology and Biotechnology, 2023, 98: 718-725. DOI: 10.1002/jctb.7275
[11] FAUNGNAWAKIJ K, KIKUCHI R, SHIMODA N, et al. Effect of thermal treatment on activity and durability of CuFe2O4-Al2O3 composite catalysts for steam re-forming of dimethyl ether[J]. Angewandte Chemie International Edition in English, 2008, 47: 9314-9317. DOI: 10.1002/anie.200802809
[12] LIU Y, HAN S Y, GUAN D D, et al. Rapid green synthesis of ZSM-5 zeolite from leached illite clay[J]. Microporous and Mesoporous Materials, 2019, 280: 324-330. DOI: 10.1016/j.micromeso.2019.02.027
[13] HAN S Y, LIU Y, YIN C R, et al. Fast synthesis of submicron ZSM-5 zeolite from leached illite clay using a seed-assisted method[J]. Microporous and Mesoporous Materials, 2019, 275: 223-228. DOI: 10.1016/j.micromeso.2018.08.028
[14] LIU H Y, YUE Y Y, SHEN T, et al. Transformation and crystallization behaviors of titanium species in synthesizing Ti-ZSM-5 zeolites from natural rectorite mineral[J]. Industrial and Engineering Chemistry Research, 2019, 58(27): 11861-11870. DOI: 10.1021/acs.iecr.9b01826
[15] PAN F, LU X C, YAN Y, et al. Synthesis of nano/micro scale ZSM-5 from kaolin and its catalytic performance[J]. Kinetics and Catalysis, 2017, 58(5): 541-548. DOI: 10.1134/S0023158417050184
[16] 蔡玉福, 周艳军, 路君凤, 等. 碱活化蒙脱土负 载铁类芬顿体系去除亚甲基蓝[J]. 复合材料学报, 2023, 40(8): 4601-4612. CAI Yufu, ZHOU Yanjun, LU Junfeng, et al. Removal of methylene blue by Fenton-like system with alkali-activated montmorillonite supported iron catalyst[J]. Acta Materiae Compositae Sinica, 2023, 40(8): 4601-4612(in Chinese).
[17] ZHAO Y H, GAO T Y, WANG Y J, et al. Zinc supported on alkaline activated HZSM-5 for aromatization reaction[J]. Reaction Kinetics, Mechanisms and Catalysis, 2018, 125: 1085-1098. DOI: 10.1007/s11144-018-1426-9
[18] 王荧光, 桂建舟, 张晓彤, 等. 纳米ZSM-5分子筛的合成与表征[J]. 光谱实验室, 2005, 22(2): 225-229. DOI: 10.3969/j.issn.1004-8138.2005.02.001 WANG Yingguang, GUI Jianzhou, ZHANG Xiaotong, et al. Synthesis and characterization of nanosized ZSM-5 zeolite[J]. Chinese Journal of Spectroscopy Laboratory, 2005, 22(2): 225-229 (in Chinese). DOI: 10.3969/j.issn.1004-8138.2005.02.001
[19] XIAO W Y, WANG F, XIAO G M. Performance of hierarchical HZSM-5 zeolites prepared by NaOH treatments in the aromatization of glycerol[J]. RSC Advances, 2015, 5: 63697-63704. DOI: 10.1039/C5RA07593A
[20] LUO C W, HUANG C, LI A, et al. Influence of reaction parameters on the catalytic performance of alkaline-treated zeolites in the novel synthesis of pyridine bases from glycerol and ammonia[J]. Industrial and Engineering Chemistry Research, 2016, 55: 893-911. DOI: 10.1021/ie504934n
[21] SHIRAZI L, JAMSHIDI E, GHASEMI M R. The effect of Si/Al ratio of ZSM-5 zeolite on its morphology, acidity and crystal size[J]. Crystal Research and Technology, 2008, 43: 1300-1306. DOI: 10.1002/crat.200800149
[22] ZHANG S G, HIGASHIMOTO S, YAMASHITA H, et al. Characterization of vanadium oxide/ZSM-5 zeolite catalysts prepared by the solid-state reaction and their photocatalytic reactivity. In situ photoluminescence, XAFS, ESR, FT-IR, and UV-Vis Investigations[J]. The Journal of Physical Chemistry B, 1998, 102: 5590-5594. DOI: 10.1021/jp981230r
[23] DUMITRIU D, BÂRJEGA R, FRUNZA L, et al. BiOx clusters occluded in a ZSM-5 matrix: Preparation, characterization, and catalytic behavior in liquid-phase oxidation of hydrocarbons[J]. Journal of Catalysis, 2003, 219: 337-351. DOI: 10.1016/S0021-9517(03)00216-1
[24] THOMMES M, KANEKO K, NEIMARK A V, et al. Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution (IUPAC Technical Report)[J]. Pure and Applied Chemistry, 2015, 87: 1051-1069. DOI: 10.1515/pac-2014-1117
[25] FAUNGNAWAKIJ K, TANAKA Y, SHIMODA N, et al. Influence of solid-acid catalysts on steam reforming and hydrolysis of dimethyl ether for hydrogen production[J]. Applied Catalysis A: General, 2006, 304: 40-48. DOI: 10.1016/j.apcata.2006.02.021
[26] SEMELSBERGER T A, OTT K C, BORUP R L, et al. Generating hydrogen-rich fuel-cell feeds from dimethyl ether (DME) using physical mixtures of a commercial Cu/Zn/Al2O3 catalyst and several solid-acid catalysts[J]. Applied Catalysis B: Environmental, 2006, 65(3-4): 291-300. DOI: 10.1016/j.apcatb.2006.02.015
-
其他相关附件
-
目的
二甲醚(DME)水蒸气重整(SRD)制氢是解决质子交换膜燃料电池(PEMFC)燃料氢气供应问题的一条合理且实用的途径。SRD是由DME水解生成甲醇和甲醇水蒸气重整生成氢气两步连续反应构成，DME水解可能是SRD的速控步骤。因此，开发成本低且高效的具有DME水解功能的酸性催化剂尤为重要。
方法以亚熔岩活化处理的廉价蒙脱土为原料，分别采用传统水热法、晶种导向法、蒸汽辅助晶化法和气溶胶晶化法制备一系列HZSM-5分子筛(分别记为TH-HZSM-5、SD-HZSM-5、SAC-HZSM-5、AC-HZSM-5)，以其为固体酸与商品化的Cu/ZnO/AlO进行物理混合制备双功能催化剂，并用于SRD反应。对所制备的材料进行了采用XRD、FT-IR、SEM、N吸附-脱附、NH-TPD和XRF等表征分析，详细探究了由不同方法合成HZSM-5构成的双功能催化剂的SRD反应性能。
结果四种方法都成功合成了以微孔为主HZSM-5分子筛，尽管不同方法合成的HZSM-5有着相似的酸强度，但酸量略有不同，酸量从高到低依次为AC-HZSM-5>SAC-HZSM-5>TH-HZSM-5>SD-HZSM-5。所有双功能催化剂的DME初始转化率都比较高，具有最多酸量的AC-HZSM-5基双功能催化剂诱发的DTH和MTH等副反应较快，导致催化剂失活迅速，反应11 h内DME转化率下降了23%，H收率最低(60%)。TH-HZSM-5具有适中的酸量以及大的比表面积、孔容和孔径，反应产物能够快速的扩散出去，进而抑制积碳的发生，因此表现出了较好的稳定性，在反应11 h时H产率最高(85%)。所有双功能催化剂都表现出较低的CH选择性，C产物在SD-HZSM-5基双功能催化剂上的选择性约为4%，随着不同方法制备HZSM-5分子筛酸量的增加，促进了DTH和MTH副反应的发生，导致其相应双功能催化剂的C产物选择性逐渐增高，其增高的程度与酸量增加的程度正相关。由于AC-HZSM-5具有较高的酸量，基于SRD的两步反应机理，DME在固体酸上水解为甲醇的速率大于甲醇在Cu/ZnO/AlO上的重整速率，从而表现出较高的甲醇选择性。当采用其他方法制备的HZSM-5的酸性减弱时，DME的水解速度减慢，所以甲醇的选择性减小。甲醇的选择性与相应固体酸的酸量是成正比的。在当前的反应条件下，Cu/ZnO/AlO催化剂对甲醇重整制氢反应具有良好的催化活性，同时极大的抑制了逆水气变换(RWGS)反应的发生，使得所有催化剂均表现出较高的CO选择性(63-89%)和较低的CO选择性(<4%)。
结论以亚熔岩活化蒙脱土为原料，采用传统水热法、晶种导向法、蒸汽辅助晶化法和气溶胶晶化法都能成功制备HZSM-5。但合成方法影响了HZSM-5的晶体尺寸、织构及酸性等，进而影响了以其为固体酸的双功能催化剂的SRD性能。传统水热法合成的TH-HZSM-5因其具有适宜的酸量以及相对大的表面积、孔容和孔径，以其为固体酸与Cu/ZnO/AlO以质量比为1:1组成的双功能催化剂表现出较好的SRD反应性能，在空速为5000 mL/(g·h)、压力为0.1 Mpa、温度为300 ℃的条件下反应11 h，DME的转化率从100%降低到95%，H收率从91%降低到85%，表现出良好的稳定性。
-
质子交换膜燃料电池(PEMFC)因其环境友好性、工作温度低、高效率和高可靠性而被认为是新一代清洁高效的动力能源之一。但由于PEMFC的燃料氢气的密度较低，难以被压缩和液化，而且易燃易爆，致使其在输送和储存上都存在着极大的不便，这限制了PEMFC的进一步发展。分析现有技术，发现通过车载重整富氢化合物“氢载体”生产氢气将是一条合理且实用的解决该问题的途径。在众多的氢载体中，由于二甲醚(DME)无毒无害，氢含量高，具有类似气体的特性和液体储存密度，以及较低的重整反应温度，是一种理想的重整富氢原料。通常，二甲醚水蒸气重整制氢反应(SRD)包括两步连续反应，第一步是DME在固体酸催化剂(氧化铝、分子筛等)上进行水解反应生成甲醇，第二步是甲醇在金属催化剂(贵金属、Cu基、Zn基催化剂等)上进行水蒸气重整反应生成CO2和H2，由于两个反应需要不同的催化活性中心，所以SRD催化剂必须同时含有这两个活性中心，是一种双功能催化剂。另外研究表明DME水解可能是SRD的速控步骤。因此，提高DME水解反应速率是加快整个反应的关键，这意味着固体酸催化剂的选择尤为重要。HZSM-5具有良好的水热稳定性，使用HZSM-5作为固体酸催化剂，DME水解在较低温度(<300 ℃)就能表现出较好的效果，但是商品化的HZSM-5通常是以含硅、铝的化学试剂合成，成本较高。以富含硅铝的天然黏土如伊利石、累托石和高岭土等代替含硅、铝的化学试剂均可成功合成HZSM-5分子筛，极大降低了HZSM-5的成本。
而蒙脱土(MMT)作为天然黏土，其含有大量的Si和Al元素，并且MMT在我国储量丰富、价格低廉。鉴于此，以MMT为原料，分别采用传统水热法、晶种导向法、蒸汽辅助晶化法和气溶胶晶化法制备了TH-HZSM-5、SD-HZSM-5、SAC-HZSM-5、AC-HZSM-5，并以其为固体酸与商品化的Cu/ZnO/Al2O3进行物理混合制备双功能催化剂，详细探究了双功能催化剂的SRD反应性能。采用XRD、FT-IR、SEM、N2吸附-脱附、NH3-TPD和XRF对HZSM-5进行表征。结果表明不同的合成方式可以显著影响HZSM-5的结构、织构和酸性等性质，进而影响HZSM-5基双功能催化剂的SRD性能。其中采用晶种导向法仅用8 h就成功合成了HZSM-5，以传统水热法合成的TH-HZSM-5为固体酸的双功能催化剂表现出最佳的SRD性能，在压力为0.1 MPa、温度300 ℃、空速 5000 mL/(g·h)的反件下，反应10 h时的DME转化率和H2收率分别为95%和84%，并且反应稳定性良好。
不同合成方法 HZSM-5 的 XRD 图 (a) TH-HZSM-5, (b) SD-HZSM-5, (c) SAC-HZSM-5, (d) AC-HZSM-5
Cu/ZnO/Al2O3和不同方法合成 HZSM-5 组 成双功能催化剂的 DME 转化率和 H2收率
计量
- 文章访问数: 120
- HTML全文浏览量: 64
- PDF下载量: 5
