CN114896719A - 一种深冷高压储氢气瓶设计方法 - Google Patents

一种深冷高压储氢气瓶设计方法 Download PDF

Info

Publication number
CN114896719A
CN114896719A CN202210439056.1A CN202210439056A CN114896719A CN 114896719 A CN114896719 A CN 114896719A CN 202210439056 A CN202210439056 A CN 202210439056A CN 114896719 A CN114896719 A CN 114896719A
Authority
CN
China
Prior art keywords
hydrogen storage
storage cylinder
pressure hydrogen
cryogenic high
fiber
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210439056.1A
Other languages
English (en)
Inventor
倪中华
严岩
张加俏
赵晓航
张帆远航
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202210439056.1A priority Critical patent/CN114896719A/zh
Publication of CN114896719A publication Critical patent/CN114896719A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/26—Composites
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/45—Hydrogen technologies in production processes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computing Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Filling Or Discharging Of Gas Storage Vessels (AREA)

Abstract

本发明公开了一种深冷高压储氢气瓶设计方法，包含：根据加注终态、供氢控制、瓶体结构及传热模型方程，得到深冷高压储氢气瓶服役工况参数；计算层合板低温下所受的目标参数域，校核现有层合板力学性能，构建深冷高压储氢气瓶设计用复合材料性能数据库；对缠绕层进行力学分析；根据缠绕层力学分析的结果，进行缠绕铺层优化设计，得到最优缠绕工艺参数。本发明实现了深冷、高压双重极端工况作用下的储氢气瓶可靠性设计，将为车载深冷高压储氢装备的研制提供全新的设计方法，推进深冷高压储氢技术在重载交通运输装备上的示范应用，对助力“双碳”发展目标、升级绿色交通装备具有重要的理论价值与应用前景。

Description

一种深冷高压储氢气瓶设计方法
技术领域
本发明属于液氢存储领域，具体涉及到一种深冷高压储氢气 瓶设计方法。
背景技术
高效安全的车载储氢技术是交通运输装备关注的重点。深冷 高压氢以突出的储氢效率、安全性，成为运输装备使用清洁能源 的最优解决方案。深冷高压储氢是指利用绝热、耐压气瓶将氢以 超临界态储存在低温(20-50K)、高压(35MPa)复合工况下， 与车载液氢相比，具有无损维持时间长、加注速度快、耐压性能 高等显著优势，能够解决目前车载液氢研究面临的漏热蒸发率高、 加注损耗大、供氢压力低等难点。然而，深冷、高压两种极端工 况的共同作用，对储氢气瓶材料性能、结构设计和试验条件提出 了极高要求，是目前车载储氢装备领域的前沿热点。
发明内容
鉴于现有技术中的缺陷或不足，期望提供一种深冷高压储氢 气瓶设计方法，通过建立热力学集总参数模型预测深冷高压储氢、 供氢全过程工况，获得气瓶设计的边界载荷，通过复合材料层力 学分析和缠绕层优化设计得到最优的深冷高压储氢气瓶设计方案。
本发明为了实现上述目的采取的技术方案是：
一种深冷高压储氢气瓶设计方法，其特征在于，包含以下步 骤：
步骤一、根据加注终态、供氢控制、瓶体结构及传热模型方 程，得到深冷高压储氢气瓶服役工况参数；
步骤二，根据步骤一得到的深冷高压储氢气瓶服役工况参数， 计算层合板低温下所受的目标参数域；根据计算得到的目标参数 域，校核现有层合板力学性能，若满足要求，则用于深冷高压储 氢气瓶制造，若不满足要求，则根据计算结果对层合板的复合材 料性能进行改性以满足要求；对改性后的层合板力学性能进行测 试，测试结果用于构建深冷高压储氢气瓶设计用复合材料性能数 据库；
步骤三，根据构建的深冷高压储氢气瓶设计用复合材料性能数 据库，对缠绕层进行力学分析；根据缠绕层力学分析的结果，进 行缠绕铺层优化设计，得到最优缠绕工艺参数。
与现有技术相比，本发明的有益效果是：
本发明提出的深冷高压储氢气瓶设计方法，综合考虑深冷与 高压双重极端环境对气瓶的影响，以服役工况预测的边界条件指 导气瓶复合材料层设计，对复合材料层进行改性处理和优化设计， 引入低温对复合材料受力的影响，构建了深冷环境下的纤维增强复合材料力学新模型，使得深冷高压储氢气瓶设计步骤更加完善， 设计计算结果更加准确。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细 描述，本申请的其它特征、目的和优点将会变得更明显：
图1为本发明一种深冷高压储氢气瓶设计方法的流程图；
图2为本发明一种深冷高压储氢气瓶设计方法的轴向环向应 力平衡图；
图3为本发明一种深冷高压储氢气瓶设计方法的环向缠绕和 螺旋缠绕线型图；
图4为本发明一种深冷高压储氢气瓶设计方法的三种铺层方 式(A型、B型、C型)下的应力和应变在壁厚上的分布，其中(a) 为三种铺层方式下每层对应的环向应力图；(b)为三种铺层方式下 每层对应的轴向应力图；(c)为三种铺层方式下每层对应的剪切应 力图；(d)为三种铺层方式下每层对应的径向应力与径向位移图；
图5为本发明一种深冷高压储氢气瓶设计方法的失效判断系 数与环向/轴向应力的分布；其中(a)为不同螺旋缠绕角度对应的蔡 吴失效准则系数图；(b)为不同螺旋缠绕角度对应的环向与轴向应 力图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以 理解的是，此处所描述的具体实施例仅仅用于解释相关发明，而 非对该发明的限定。另外还需要说明的是，为了便于描述，附图 中仅示出了与发明相关的部分。
需要说明的是，在不冲突的情况下，本申请中的实施例及实 施例中的特征可以相互组合。下面将参考附图并结合实施例来详 细说明本申请。
请参考图1，一种深冷高压储氢气瓶设计方法，包括四个步骤， 具体是储氢供氢过程理论建模、复合材料设计、纤维缠绕设计以 及建模验证。
第1步：储氢供氢过程理论建模：根据所建的储氢供氢过程 理论模型，得到深冷高压储氢气瓶服役工况参数。
具体步骤包括：
11.考虑节流效应、瓶体结构及绝热层的影响，构建深冷高压 储氢气瓶传热模型；
Figure BDA0003609223180000031
其中，
Figure BDA0003609223180000032
为能量守恒方程中热量输入项；
Figure BDA0003609223180000033
为能量守恒方程中 热量输出项；ms为饱和氢质量；Cs为饱和氢能量密度；ΔT为温差； m为瓶内储氢质量；U为氢的内能；h为焓值；t为时间。
在气瓶传热模型，考虑气瓶与外壳之间的真空绝热层、支撑 结构漏热规律，结合气相供氢过程需要的辅助加热功率，形成热 量输入项
Figure BDA0003609223180000041
根据气瓶保冷要求，优化供氢温度、确定节流制冷 量，形成热量输出项
Figure BDA0003609223180000042
建立系统传热模型。
12.根据气瓶传热模型，构建气瓶供氢过程集总参数模型；
Figure BDA0003609223180000043
Figure BDA0003609223180000044
式中，h为焓值，P为气瓶受到的内部压力，t为时间，T为气瓶 内液氢温度，ρ为质量密度。
以微分思想建立由多个平衡状态构成的热力学瞬态模型，即 集总参数模型。以供氢质量流作为系统质量源项，利用焓值与密 度(质量)变化率计算系统的瞬态温度和压强，通过20-293K、 0.1-20MPa区间内的饱和状态曲线确定气-液-超临界多相流模型 中各相态的质量分布，以质量守恒方程作为迭代计算的判定依据， 确定多相流传质结果，形成描述深冷高压储氢供氢过程的集总参 数热力学模型。
13.根据集总参数模型气瓶复合材料层受力平衡推导出深冷高 压储氢气瓶传质模型；
Figure BDA0003609223180000045
Figure BDA0003609223180000046
式中，P为气瓶受到的内部压力，t为时间，T为气瓶内液氢 温度，σ为复合材料层所受到应力，T′为复合材料层温度，Φ和Ψ为 方程修正系数。
14.根据深冷高压储氢气瓶传质模型，输出深冷高压储氢气瓶 服役工况，所述服役工况包括压力、温度等。
Figure BDA0003609223180000051
Figure BDA0003609223180000052
式中，P为气瓶受到的内部压力，t为时间，T为气瓶内液氢温度， σ为复合材料层所受到应力，T′为复合材料层温度，Φ和Ψ为方程 修正系数。第二步：复合材料设计：根据步骤一得到的深冷高压 储氢气瓶服役工况参数，计算层合板低温下所受应力与产生的变 形等目标参数域；根据计算得到的目标参数域，校核现有层合板 力学性能，若满足要求，则用于深冷高压储氢气瓶制造，若不满 足要求，则根据计算结果对层合板的复合材料性能进行改性以满 足要求；对改性后的层合板力学性能进行测试，测试结果用于构 建深冷高压储氢气瓶设计用复合材料性能数据库。
具体步骤包括：
21.根据步骤一得到的深冷高压储氢气瓶服役工况参数，计算 层合板低温下所受应力与产生变形的目标参数域；
22.根据计算得到的目标参数域，校核现有层合板力学性能。
根据服役工况得到了目标参数域，比如在服役过程中气瓶复 合材料层受到的拉应力为1500-2000MPa，现有层合板的抗拉强度 为1300或者1700MPa，层合板的抗拉强度小于目标参数域最大值 (<2000MPa)，则表示使用过程中会超过材料极限，存在破坏可能 性，不满足使用需求，需要进行改性处理。
改性方法包括纤维改性与树脂改性，纤维改性采用表面氧化、 等离子处理或化学沉积，在纤维表面形成沟痕或突起，增加纤维 表面粗糙度；树脂改性通过在树脂中添加有机改性剂和无机改性 剂，有机改性剂包括聚醚砜、聚氨酯或聚乙二醇等，无机改性剂 包括碳纳米管、石墨烯或二氧化硅。并采用X射线衍射仪对改性 后树脂进行变温测试，得到树脂和纤维不同方向上的低温热膨胀 系数。
24.建立纤维增强树脂基层合板低温力学模型，再根据树脂和 纤维热膨胀系数测试结果，构建静力平衡方程，求解层合板纵向 和横向热膨胀系数，求解公式为：
Figure BDA0003609223180000061
α2＝vf(1+μf)αf+vm(1+μm)αm-(μfvf+μmvm)α1
其中，α1为层合板的纵向热膨胀系数；α2为层合板的横向热膨胀系 数；αf1为层合板中纤维的纵向热膨胀系数；αm为层合板中树脂的热 膨胀系数；αf为层合板中纤维的热膨胀系数；Ef1为层合板中纤维的 纵向弹性模量；vf1为层合板中纤维的纵向体积分数；vm为层合板中 树脂的体积分数；vf为层合板中纤维的体积分数；μm为层合板中树 脂的泊松比；μf为层合板中纤维的泊松比。
25.对改性后层合板力学性能进行测试，测试结果包括层合板 的抗拉强度、弹性模量、压缩强度等力学性能参数，测试结果用 于构建深冷高压储氢气瓶设计用复合材料性能数据库。气瓶工作 过程中工况参数会产生变化，比如内部压力30-35MPa波动，对应 的复合材料内部产生1000-1200Mpa的应力和x-y的应变，这些就 是目标参数域，而通过材料改性就是希望材料能够满足这个区域， 比如这个过程中至少要求抗拉强度大于1200MPa。
第三步：纤维缠绕设计：在层合板设计基础上，将深冷工况 引起的温度载荷项引入纤维缠绕理论，构建各向异性复合材料的 静力学平衡方程，确定深冷高压储氢气瓶的基本尺寸，然后计算 纤维缠绕层各点、各方向的应力应变，并对缠绕工艺进行优化改 进。
31.容器纤维缠绕层及其外绝热结构的热力学建模：复合材料 表面的温度是由多层真空绝缘材料的性能决定的，这部分主要由 热对流和热辐射来计算。靠近碳纤维复合材料层部分的温度实际 上是非常低的，其传热形式主要为热传导。纤维复合层的热传导 计算公式为：
dQc+dQg＝dQ
其中，dQc代表以导热方式传入微元体的净热量，dQg代表微元 体内热源生成的热量，dQ表示微元体内能的增量，其中：
Figure BDA0003609223180000071
其中，x,y,z为空间坐标方向；q为热量，t为时间。
根据傅里叶定律：
Figure BDA0003609223180000072
其中，λ为传热系数；T为环境温度。
可以进一步将写dQc成：
Figure BDA0003609223180000073
同时，单位时间内微元体内热源生成的热量可以表示为：
dQg＝qVdxdydzdt
其中，qV为单位体积内热量，t指时间。
单位时间单位面积的微元体内能的增量为：
Figure BDA0003609223180000074
其中，ρ为常数，u指动量，c为能量密度。
联立式上式，可以得到：
Figure BDA0003609223180000075
令
Figure BDA0003609223180000076
代入上式可得：
Figure BDA0003609223180000077
写成圆柱坐标形式时，该式可以表示为：
Figure BDA0003609223180000081
其中，r为内壁至外壁间任一半径，r0＜r＜ra，
Figure BDA0003609223180000082
为螺旋缠绕角。
所以在圆柱坐标系下深冷高压储氢容器的导热控制方程可以 写为：
Figure BDA0003609223180000083
其中，
Figure BDA0003609223180000084
表示内能产生速率，k表示导热系数，ξ表示热扩散系 数，物理意义是描述物体的热惯性，热扩散系数越大表示热惯性 越小，物体达到与周围环境热平衡的状态越快。
由于管道的长、轴对称和稳态条件，在不产生热量的情况下， 管道内的温度分布仅为半径的函数。导热控制方程可以简化为：
Figure BDA0003609223180000085
表达为积分的形式为：
T＝A+Blnr
式中，A和B为积分常数。
外表面在环境温度T下暴露于自由对流，内表面在内壁温度Tf下暴露于对热流体的强制对流。根据对流换热方程可以列出两个 边界处的传热方程为：
Figure BDA0003609223180000086
其中，T为环境温度，Tf为气瓶内壁温度，T∞为气瓶外部温度，h0,ha分别是复合材料层内外表面的平均对流传热系数，k表示导热系数， r0和ra分别表示深冷高压储氢气瓶的内外半径；利用上面边界条件 可以求出积分常数A和B。
内部温度设为Tf，采用20K温度下的深冷高压氢，外部温度T∞ 为靠近碳纤维复合材料层的最近一层绝热材料的温度为50K。
32.温度载荷项引入计算：
碳纤维在各层中均按特定的方向进行缠绕。储氢容器的内部 温度用Tf表示，参考温度设为Tref，表示容器的最初温度，这里设 置为298K，任意径向位置的工作温度用T(r)表示，所以沿径向的 温差可表示为：
ΔT(r)＝T(r)-Tref
容器外壁产生的应变是由加载的内压和温差导致的热应力共 同作用引起的，可以表示为：
{εi}＝[S]{σi}+αΔT(r)
Figure BDA0003609223180000091
Figure BDA0003609223180000092
其中，
Figure BDA0003609223180000093
表示在热力耦合作用下的柱坐标下的应变矩阵；{εi}表 示沿储氢容器径向方向上第i层纤维应变εi组成的矩阵；[S]表示复合 材料在柱坐标下的柔度矩阵；{σi}表示沿储氢容器径向方向上第i 层纤维应力σi组成的矩阵；α是热膨胀系数；ΔT为温差；d表示微分 符号；u表示位移；
Figure BDA0003609223180000094
表示第i层径向上的应力；
Figure BDA0003609223180000095
表示第i层环 向上的应力；
Figure BDA0003609223180000096
表示第i层径向上的应变；
Figure BDA0003609223180000097
表示第i层环向上的 应变；
Figure BDA0003609223180000098
表示第i层径向上的应变。
使用高阶矩阵快速算法求解上述方程：
Figure BDA0003609223180000101
其中，a,b,d,k为系数；n为纤维缠绕层数；ε0为轴向应变；γ0为 切应变。
深冷高压储氢气瓶的内外半径分别用r0和ra表示。此外，容器 表面温度固定在环境温度Tra。由于氢流体的强制对流作用，假定 内表面温度与流体温度Tr0.一致。根据柱坐标系下的薄壁圆筒的导 热控制方程，储氢容器壁面上的温度沿半径方向的分布如下：
Figure BDA0003609223180000102
其中，h0表示外表面的平均对流换热系数，由所述储氢容器外表 面的绝对压力决定，根据经验结果取为0.1653W/(m2℃)；ha是 内壁面的对流换热系数，取为50W/(m2℃)；k是复合材料的导 热系数。
33.复合材料静力学平衡方程构建：
如图2所示容器受力情况，可以列出轴向和环向内力的受力 平衡方程，轴向平衡表示为：
Figure BDA0003609223180000103
圆周平衡表示为：
Figure BDA0003609223180000104
其中，P表示容器承受的工作压力，d为容器的外半径，t为缠绕 层厚度，Nl是薄膜的轴向内力，Nθ是薄膜的环向内力，
Figure BDA0003609223180000118
为螺旋缠 绕角度，ri为第i层的半径。对上述两式进行化简可以得到：
Nl＝1/2rP，Nθ＝rP
Nθ＝2Nl
其中，r为容器半径。
设计过程中采用环向缠绕和螺旋缠绕相结合的方式，如图3 所示。考虑深冷高压储氢瓶的内衬承担约5％的机械载荷，先计算 出螺旋和环向纤维的厚度：
将环向和螺旋向纤维应力分解为轴向和环向应力并叠加，得 到平衡方程
Figure BDA0003609223180000111
Figure BDA0003609223180000112
其中，α是螺旋缠绕的角度，σf是纤维的轴向拉力，tθ是环向缠绕 纤维厚度，
Figure BDA0003609223180000113
为螺旋缠绕纤维厚度。
考虑在容器的极限压力下的平衡公式，此时纤维达到它的设 计应力。将式Nl＝1/2rP，Nθ＝rP代入得：
Figure BDA0003609223180000114
其中，Pm是容器承受的极限压力，R是容器的外半径，σd是纤维 的设计强度。
Figure BDA0003609223180000115
求解得到纤维厚度：
Figure BDA0003609223180000116
Figure BDA0003609223180000117
34.缠绕工艺优化
采用控制变量和对比实验的方法，使用正交实验法分析纤维 缠绕角度、厚度、张紧力、铺层方式等多种因素对气瓶缠绕层强 度的影响规律，求解热力耦合作用下缠绕角与复合材料层内应力 分布的关系。进一步地，使用Tsai-Wu、Tsai-Hill、Hoffman等强 度失效准则开展气瓶失效预测，分析储氢气瓶在极限工况下的破 坏位置及失效原因，根据失效分析结果，对危险区域进行增强处 理、安全区域进行轻量化设计。
第四步：建模验证
进一步地，对深冷高压储氢气瓶进行附属部件设计、三维建 模及仿真验证。用有限元仿真进行铝合金与碳纤维异种材料之间 受力与变形行为进行仿真，定义不同的单元属性，保证支撑强度 的情况下减少铝合金的用量，进一步降低储氢气瓶的质量。其次， 对铝合金碳纤维双层结构进行界面粘附力和载荷分布测试，结合 两种材料的热膨胀系数确定碳纤维缠绕时需要施加的预应力，减 小二者之间因低温产生的热应力。绝热层采用高真空变密度多层 绝热方法，在辐射热流占主导的高温侧使用较大的层密度来减少 辐射换热，而在低温侧使用较小的层密度来减少固体材料导热， 优化多层绝热材料的整体性能。设计深冷高压储氢瓶专用支撑结 构，支撑结构设计为螺旋形状，增加导热长度，减小两端体积， 从而减少接触面积，支撑结构与复合材料层之间采用柔性接触， 有利于减小因储氢气瓶形变而产生的应力。进一步地，根据轻量 化设计结果和缠绕工艺参数进行深冷高压储氢气瓶瓶身三维建模， 根据服役工况参数在有限元分析软件内进行环境设置，并在所述环境设置下进行样瓶耐压实验、爆破实验、晃动实验等。
实施例：容器容积为141L，容器的工作压力为35MPa，容器 最低工作温度为20K。相应的，爆破压力设为78.75MPa，是工作 压力的2.25倍。采用1/2标准椭球封头，容器的长径比设为4.8。 半径可以由容器的体积推导出来。其长度和直径计算后分别为 1568mm和410mm。
为了得到环向缠绕角，选择容器直径的5％作为碳纤维缠绕带 的宽度：
Figure BDA0003609223180000131
其中，αhoop是环向缠绕角度，dw是碳纤维缠绕带的宽度。
在开始应力分析前，采用网格理论确定所需的纤维层壁厚的厚 度，并将其作为循环计算的初始值。环向和螺旋方向的厚度如下:
Figure BDA0003609223180000132
Figure BDA0003609223180000133
其中，k表示复合材料承受的机械载荷比例(少于95％)，R是容器 的内半径，Pm是最小爆破压力，vf是纤维体积分数，σd是碳纤维 的设计压力(即纤维的拉伸强度)，
Figure BDA0003609223180000134
是螺旋缠绕角度，这里选择5°。
选择Halpin-Tsai方法预测复合材料弹性力学性能，碳纤维的 弹性性能在20-293K的温度范围内是恒定的，环氧树脂基体的弹 性模量从室温(3.6GPa)到20K(8.2GPa)呈线性变化。T300/914 碳纤维/树脂复合材料的纤维体积比为0.57。
代入上述参数，实施例中复合材料层数为：
Figure BDA0003609223180000135
或者20层
Figure BDA0003609223180000136
或者38层
根据失效分析计算公式可知，实施例中的深冷高压储氢容器 在使用58层结构时会发生失效。当层数增加到75层时，才能满 足设计条件。
进一步，对缠绕层和缠绕方向进行分析，共设计三种缠绕结 构。A型，每25层使用相同的缠绕角，每一层包括相对于圆筒的 轴线正负相同的角度的两个缠绕层。B型，每3层构成一个循环， 这两种缠绕方式的螺旋缠绕层与环向缠绕层的比例为1:2。C型， 在A型基础上对层数进行了相应的调整，每26层使用相同的缠绕 角。图4为三种不同铺层顺序的深冷高压储氢容器的应力和位移。 三种铺层顺序构成的储氢容器收到的应力范围相同，但分布不同。 图4(a)可以看出，C型在26-49层受环向应力较小，并且应力 较小部分发生在螺旋缠绕部分。从图4(b)可以看出，C型中间层 承受较大的轴向应力，这也是在螺旋缠绕部分发生的。如果螺旋 缠绕和环向缠绕分别承受更多的轴向应力和环向应力，碳纤维的 强度将得到充分利用。从图4(d)可以看出，C型的径向位移比A 型的径向位移减少了约5.27-6.29％。
在本实施例中，A型缠绕方法始终是首选。但在第30层时， A型Tsai-Wu失效判断系数超过1。因此，在本实施例中，对A型 缠绕方式中的层数进行简单调整，改为C型，此时C型的Tsai-Wu 失效判断系数包络在A和B中间，稳定性较高。选择C型时， Tsai-Wu失效判断系数小于1，可以保证容器的强度。
进一步地，对实施例设计结果进行优化。实施例中螺旋缠绕 的角度是不断变化的，选取各层在各角度下的Tsai-Wu失效判断 系数的最大值，绘制Tsai-Wu失效判断系数随缠绕角的变化趋势 图，如图5所示。图5(a)中，Tsai-Wu失效判断系数随着纵向缠绕 角的增加而增加。该值在角度区间0-14°内小于1，满足容器设计 要求。因此，螺旋缠绕角度应选择在区间0-14°。在后续计算中， 选取螺旋缠绕角为5°，对应的Tsai-Wu失效判断系数为0.984。随 着螺旋缠绕角的增大，Tsai-Wu失效判断系数迅速增大。当完全是 环向缠绕时，容器的轴向应力由垂直于碳纤维的方向承受。此时， 容器的承载能力最弱。图5(b)显示了最外层碳纤维上的轴向应 力和环向应力随角度变化曲线。在0-14°范围内，最外层轴向应力 与环向应力的差异最显著。缠绕角度为89°的环形缠绕时，环向应 力主要由碳纤维方向承受，而轴向应力主要由与碳纤维垂直的方 向承受。因此，在容器设计中应优先选择轴向应力较小的0-14°范 围。最后，本实施例选择的优化缠绕角度为5°。

Claims (10)

1.一种深冷高压储氢气瓶设计方法，其特征在于，包含以下步骤：
步骤一、根据加注终态、供氢控制、瓶体结构及传热模型方程，得到深冷高压储氢气瓶服役工况参数；
步骤二，根据步骤一得到的深冷高压储氢气瓶服役工况参数，计算层合板低温下所受的目标参数域；根据计算得到的目标参数域，校核现有层合板力学性能，若满足要求，则用于深冷高压储氢气瓶制造，若不满足要求，则根据计算结果对层合板的复合材料性能进行改性以满足要求；对改性后的层合板力学性能进行测试，测试结果用于构建深冷高压储氢气瓶设计用复合材料性能数据库；
步骤三，根据构建的深冷高压储氢气瓶设计用复合材料性能数据库，对缠绕层进行力学分析；根据缠绕层力学分析的结果，进行缠绕铺层优化设计，得到最优缠绕工艺参数。
2.根据权利要求1所述的一种深冷高压储氢气瓶设计方法，其特征在于，步骤一，根据节流效应、瓶体结构及能量守恒方程，得到深冷高压储氢气瓶服役工况参数，包括：
根据节流效应、瓶体结构及绝热层，构建深冷高压储氢气瓶传热模型；
以加注终态为初始条件，根据能量守恒方程和构建的传热模型，构建深冷高压储氢气瓶供氢过程集总参数模型；
根据集总参数模型推导出深冷高压储氢气瓶传质模型；
根据深冷高压储氢气瓶传质模型，输出深冷高压储氢气瓶服役工况，所述服役工况包括压力和温度。
3.根据权利要求2所述的一种深冷高压储氢气瓶设计方法，其特征在于，步骤三中，对缠绕层进行力学分析包括：
根据导热控制方程，计算任意径向位置的工作温度，其中导热控制方程为：
T(r)＝A+Blnr
其中，r为气瓶内壁至外壁间任一半径，r0＜r＜ra；T(r)为半径为r时壁面的温度；A和B为积分常数；
确定沿径向的温差：
ΔT(r)＝T(r)-Tref
其中，ΔT(r)为沿径向的温差；Tref为参考温度，一般为室温；
确定容器外壁产生的应变：
{εi}＝[S]{σi}+αΔT(r)
其中，{εi}表示沿储氢容器径向方向上第i层纤维应变εi组成的矩阵；[S]表示复合材料在柱坐标下的柔度矩阵；{σi}表示沿储氢容器径向方向上第i层纤维应力σi组成的矩阵；α是热膨胀系数；
根据缠绕层力学分析的结果，进行缠绕铺层优化设计，得到最优缠绕工艺参数，包括：
根据容器受力情况，计算螺旋和环向纤维的厚度：
Figure FDA0003609223170000021
Figure FDA0003609223170000022
其中，tθ是环向缠绕纤维厚度；tθ为螺旋缠绕纤维厚度；Pm是容器承受的极限压力；R是容器的外半径；σd是纤维的设计强度；θ是螺旋缠绕的角度。
4.根据权利要求3所述的一种深冷高压储氢气瓶设计方法，其特征在于，任意径向位置的工作温度T(r)为：
Figure FDA0003609223170000023
其中，h0表示外表面的平均对流换热系数；ha是内壁面的对流换热系数；k是复合材料的导热系数；r0和ra分别为储氢气瓶的内外半径；Tra和Tr0分别是储氢气瓶的内外温度。
5.根据权利要求1-4任一所述的一种深冷高压储氢气瓶设计方法，其特征在于，
所述步骤三还包括：
根据得到的最优缠绕工艺参数，构建缠绕层模型，对缠绕层模型进行热力耦合失效分析，完成轻量化瓶身设计；
深冷高压储氢气瓶设计方法还包括：
步骤四，根据轻量化设计结果和服役工况参数进行深冷高压储氢气瓶瓶身三维建模，根据服役工况参数在有限元分析软件内进行环境设置，并在设置环境下进行样瓶耐压实验、爆破实验及晃动实验。
6.根据权利要求5所述的一种深冷高压储氢气瓶设计方法，其特征在于，
步骤三中，对缠绕层模型进行热力耦合失效分析，完成轻量化瓶身设计的方法，包括；
对缠绕层模型进行热力耦合失效分析，得到材料-温度-应力共同作用下的气瓶结构损伤机理与失效规律；
再根据所述损伤机理与失效规律进行深冷高压储氢气瓶轻量化设计；
热力耦合失效分析的失效准则为：
Figure FDA0003609223170000031
其中，F11、F12、F22、F33、F1及F2为系数；σ1为x方向上拉压应力；σ2为y方向上拉压应力；τ12为正轴方向上材料的剪切强度；
各系数计算方法为：
Figure FDA0003609223170000032
其中，Xt和Xc表示碳纤维复合材料在纤维方向上的拉伸强度和压缩强度，Yt和Yc表示碳纤维复合材料在垂直纤维方向上的拉伸强度和压缩强度，S为剪切强度。
7.根据权利要求1所述的一种深冷高压储氢气瓶设计方法，其特征在于，步骤一中，所述传热模型方程为：
Figure FDA0003609223170000041
其中，
Figure FDA0003609223170000042
为能量守恒方程中热量输入项；
Figure FDA0003609223170000043
为能量守恒方程中热量输出项；ms为饱和氢质量；Cs为饱和氢能量密度；ΔT为温差；m为瓶内储氢质量；U为氢的内能；h为焓值。
8.根据权利要求1所述的一种深冷高压储氢气瓶设计方法，其特征在于，步骤二中，对复合材料性能进行改性，分为纤维改性与树脂改性；所述纤维改性采用表面氧化、等离子处理或化学沉积，在纤维表面形成沟痕或突起，增加纤维表面粗糙度；所述树脂改性通过在树脂中添加有机改性剂和无机改性剂，所述有机改性剂包括聚醚砜、聚氨酯或聚乙二醇等，所述无机改性剂包括碳纳米管、石墨烯或二氧化硅。
9.根据权利要求8任一所述的一种深冷高压储氢气瓶设计方法，其特征在于，测试所述改性纤维与改性树脂的热膨胀系数，改性后复合材料热膨胀系数计算公式为：
Figure FDA0003609223170000044
α2＝vf(1+μf)αf+vm(1+μm)αm-(μfvf+μmvm)α1
其中，α1为层合板的纵向热膨胀系数；α2为层合板的横向热膨胀系数；αf1为层合板中纤维的纵向热膨胀系数；αm为层合板中树脂的热膨胀系数；αf为层合板中纤维的热膨胀系数；Ef1为层合板中纤维的纵向弹性模量；vf1为层合板中纤维的纵向体积分数；vm为层合板中树脂的体积分数；vf为层合板中纤维的体积分数；μm为层合板中树脂的泊松比；μf为层合板中纤维的泊松比。
10.根据权利要求5所述的一种深冷高压储氢气瓶设计方法，其特征在于，步骤四具体包括：
使用有限元仿真进行铝合金与碳纤维异种材料之间受力与变形行为进行仿真，定义不同的单元属性，保证支撑强度的情况下减少铝合金的用量，进一步降低储氢气瓶的质量；
对所述铝合金碳纤维双层结构进行界面粘附力和载荷分布测试，结合两种材料的热膨胀系数确定碳纤维缠绕时需要施加的预应力，减小二者之间因低温产生的热应力；绝热层采用高真空变密度多层绝热方法，在辐射热流占主导的高温侧使用较大的层密度来减少辐射换热，而在低温侧使用较小的层密度来减少固体材料导热，优化多层绝热材料的整体性能；
设计深冷高压储氢瓶专用支撑结构，支撑结构设计为螺旋形状，增加导热长度，减小两端体积，从而减少接触面积，支撑结构与复合材料层之间采用柔性接触，有利于减小因储氢气瓶形变而产生的应力。
CN202210439056.1A 2022-04-22 2022-04-22 一种深冷高压储氢气瓶设计方法 Pending CN114896719A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210439056.1A CN114896719A (zh) 2022-04-22 2022-04-22 一种深冷高压储氢气瓶设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210439056.1A CN114896719A (zh) 2022-04-22 2022-04-22 一种深冷高压储氢气瓶设计方法

Publications (1)

Publication Number Publication Date
CN114896719A true CN114896719A (zh) 2022-08-12

Family

ID=82717880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210439056.1A Pending CN114896719A (zh) 2022-04-22 2022-04-22 一种深冷高压储氢气瓶设计方法

Country Status (1)

Country Link
CN (1) CN114896719A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116100841A (zh) * 2023-01-13 2023-05-12 江苏集萃复合材料装备研究所有限公司 一种用于高压氢气瓶的制备方法
CN117390937A (zh) * 2023-12-12 2024-01-12 华中科技大学 一种快充条件下含缺陷储氢气瓶损伤预测方法及系统
CN117610385A (zh) * 2024-01-24 2024-02-27 合肥通用机械研究院有限公司 考虑强度和疲劳寿命的ⅳ型储氢气瓶铺层设计方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116100841A (zh) * 2023-01-13 2023-05-12 江苏集萃复合材料装备研究所有限公司 一种用于高压氢气瓶的制备方法
CN116100841B (zh) * 2023-01-13 2024-01-26 江苏集萃复合材料装备研究所有限公司 一种用于高压氢气瓶的制备方法
CN117390937A (zh) * 2023-12-12 2024-01-12 华中科技大学 一种快充条件下含缺陷储氢气瓶损伤预测方法及系统
CN117390937B (zh) * 2023-12-12 2024-03-05 华中科技大学 一种快充条件下含缺陷储氢气瓶损伤预测方法及系统
CN117610385A (zh) * 2024-01-24 2024-02-27 合肥通用机械研究院有限公司 考虑强度和疲劳寿命的ⅳ型储氢气瓶铺层设计方法
CN117610385B (zh) * 2024-01-24 2024-04-09 合肥通用机械研究院有限公司 考虑强度和疲劳寿命的ⅳ型储氢气瓶铺层设计方法

Similar Documents

Publication Publication Date Title
CN114896719A (zh) 一种深冷高压储氢气瓶设计方法
Parnas et al. Design of fiber-reinforced composite pressure vessels under various loading conditions
Kabir Finite element analysis of composite pressure vessels with a load sharing metallic liner
Onder et al. Burst failure load of composite pressure vessels
Chapelle et al. Optimal design of a Type 3 hydrogen vessel: Part I—Analytic modelling of the cylindrical section
Gentilleau et al. Numerical study of influence of temperature and matrix cracking on type IV hydrogen high pressure storage vessel behavior
Hocine et al. Experimental and analytical investigation of the cylindrical part of a metallic vessel reinforced by filament winding while submitted to internal pressure
Zhao et al. Analysis of multilayered carbon fiber winding of cryo-compressed hydrogen storage vessel
Tapeinos et al. Evaluation of the mechanical performance of a composite multi-cell tank for cryogenic storage: Part II–Experimental assessment
Tapeinos et al. Evaluation of the mechanical performance of a composite multi-cell tank for cryogenic storage: Part I-Tank pressure window based on progressive failure analysis
Xia et al. Effects of thickness and winding angle of reinforcement laminates on burst pressure capacity of thermoplastic composite pipes
Shi et al. Partially-plastic theoretical model of thermoplastic composite pipes and comparison of composite failure criteria
Daghighi et al. Design methods for variable-stiffness super-ellipsoidal pressure vessels under thermomechanical loading
Tapeinos et al. Design and analysis of a multi-cell subscale tank for liquid hydrogen storage
Błachut Buckling of multilayered metal domes
Musthak et al. Prediction of transverse directional strains and stresses of filament wound composite pressure vessel by using higher order shear deformation theories
Wang et al. Mechanical analysis of thick-walled filament wound composite pipes under pure torsion load: safety zones and optimal design
Liu et al. Analytical model for the winding process-induced residual stresses of the multilayered filament wound cylindrical composite parts
Saha et al. Gas permeation of composites subjected to impact for cryogenic applications
Zu et al. Minimum weight design of helically and hoop wound toroidal hydrogen storage tanks with variable slippage coefficients
Marchetta et al. Small-scale room-temperature-vulcanizing-655/aerogel cryogenic liquid storage tank for space applications
Tian et al. Progressive damage analysis for multiscale model of linerless composite cryotank and integrated design
Saha et al. Gas permeability and flexural strength of impacted composites for cryogenic propellant tanks
Fu et al. Strength and fatigue study of on-board type III cryo-compressed hydrogen storage cylinder
Adams et al. Integrated Insulation System for Cryogenic Automotive Tanks (iCAT)

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination