CN112460464A - 一种基于纤维缠绕工艺的高压储氢容器的设计方法 - Google Patents

一种基于纤维缠绕工艺的高压储氢容器的设计方法 Download PDF

Info

Publication number
CN112460464A
CN112460464A CN202011162433.9A CN202011162433A CN112460464A CN 112460464 A CN112460464 A CN 112460464A CN 202011162433 A CN202011162433 A CN 202011162433A CN 112460464 A CN112460464 A CN 112460464A
Authority
CN
China
Prior art keywords
hydrogen storage
pressure hydrogen
storage container
fiber
pressure
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.)
Granted
Application number
CN202011162433.9A
Other languages
English (en)
Other versions
CN112460464B (zh
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.)
Jiangsu Guofu Hydrogen Energy Technology Equipment Co Ltd
Original Assignee
Jiangsu Guofu Hydrogen Energy Technology Equipment Co Ltd
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 Jiangsu Guofu Hydrogen Energy Technology Equipment Co Ltd filed Critical Jiangsu Guofu Hydrogen Energy Technology Equipment Co Ltd
Priority to CN202011162433.9A priority Critical patent/CN112460464B/zh
Publication of CN112460464A publication Critical patent/CN112460464A/zh
Application granted granted Critical
Publication of CN112460464B publication Critical patent/CN112460464B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C1/00—Pressure vessels, e.g. gas cylinder, gas tank, replaceable cartridge
    • F17C1/005—Storage of gas or gaseous mixture at high pressure and at high density condition, e.g. in the single state phase
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C1/00—Pressure vessels, e.g. gas cylinder, gas tank, replaceable cartridge
    • F17C1/02—Pressure vessels, e.g. gas cylinder, gas tank, replaceable cartridge involving reinforcing arrangements
    • F17C1/04—Protecting sheathings
    • F17C1/06—Protecting sheathings built-up from wound-on bands or filamentary material, e.g. wires
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2203/00—Vessel construction, in particular walls or details thereof
    • F17C2203/01—Reinforcing or suspension means
    • F17C2203/011—Reinforcing means
    • F17C2203/012—Reinforcing means on or in the wall, e.g. ribs
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2203/00—Vessel construction, in particular walls or details thereof
    • F17C2203/06—Materials for walls or layers thereof; Properties or structures of walls or their materials
    • F17C2203/0634—Materials for walls or layers thereof
    • F17C2203/0658—Synthetics
    • F17C2203/0663—Synthetics in form of fibers or filaments
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2221/00—Handled fluid, in particular type of fluid
    • F17C2221/01—Pure fluids
    • F17C2221/012—Hydrogen
    • 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/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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/30—Hydrogen technology
    • Y02E60/32—Hydrogen storage

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Filling Or Discharging Of Gas Storage Vessels (AREA)

Abstract

本发明公开了一种基于纤维缠绕工艺的高压储氢容器的设计方法，具体包括如下步骤：步骤1：确定高压储氢容器的基本参数；步骤2：获取纤维复合层及内衬层的材料性能参数；步骤3：计算纤维复合层的缠绕厚度值和内衬层的内衬厚度值：（1）根据步骤1和步骤2确定设计准则和设计参数
Figure DEST_PATH_IMAGE001
；（2）根据高压储氢容器在极限压力下的平衡公式和Ⅲ型高压储氢容器的承压特点列平衡方程式；（3）根据平衡方程式求解：
Figure 976398DEST_PATH_IMAGE002
、
Figure DEST_PATH_IMAGE003
、
Figure 319917DEST_PATH_IMAGE004
；（4）求解
Figure DEST_PATH_IMAGE005
；步骤4：利用步骤2，验证步骤3中计算结果：
Figure 511864DEST_PATH_IMAGE002
、
Figure 122974DEST_PATH_IMAGE003
、
Figure 476595DEST_PATH_IMAGE004
、
Figure 274787DEST_PATH_IMAGE005
的准确性。设计中对计算结果进行校核和修正，可以提升各计算结果的准确度，为高压储氢容器的设计制造、结构优化提供了理论依据，降低制造成本。

Description

一种基于纤维缠绕工艺的高压储氢容器的设计方法
技术领域
本发明涉及高压储氢技术领域，尤其涉及一种基于纤维缠绕工艺的高压储氢容器的设计方法。
背景技术
高压储氢技术具有充装方便、成本低等优点，是目前最常见、最成熟的储氢技术。在高压储氢技术中，纤维增强复合材料以其优异性能常被用来作为高压储氢容器的主要承载层，采用纤维增强复合材料制成的纤维缠绕式高压储氢容器因具有较高的比强度、轻质、耐腐蚀性、耐疲劳等特点，被广泛应用于火箭、卫星、新能源汽车、医疗设备等产品上。
目前，在采用纤维增强复合材料制成的纤维缠绕式高压储氢容器的传统设计中，缺乏对计算结果：纤维复合层的螺旋缠绕纤维厚度、纤维复合层的环向缠绕纤维厚度、内衬层的内衬厚度、纤维复合层的纤维预应力的校核和修正，而传统设计中使用网格理论来计算上述结果时，计算结果往往与工程实际应用相差很大，造成制造成本增加等缺陷。
发明内容
本发明所需解决的技术问题是：提供一种能对计算结果进行校核和修正的基于纤维缠绕工艺的高压储氢容器的设计方法。
为解决上述问题，本发明采用的技术方案是：所述的一种基于纤维缠绕工艺的高压储氢容器的设计方法，具体包括如下步骤：
步骤1：确定高压储氢容器的基本参数：高压储氢容器承受的工作压力Pw、高压储氢容器承受的极限压力Pm、高压储氢容器的容积V、高压储氢容器的筒身长L、高压储氢容器的内半径R、纤维复合层螺旋缠绕的缠绕角α；
步骤2：获取纤维复合层及内衬层的材料性能参数：纤维复合材料的抗拉强度σd、纤维复合材料的抗拉模量EF、内衬材料的弯曲屈服强度σip、内衬材料的弹性模量E0、内衬材料的泊松比μ；
步骤3：计算纤维复合层的缠绕厚度值和内衬层的内衬厚度值；
(1)根据步骤1中高压储氢容器的基本参数和步骤2中纤维复合层及内衬层的材料性能参数确定设计准则和设计参数λ；其中，λ为III型高压储氢容器的纤维强度利用系数；
(2)根据高压储氢容器在极限压力下的平衡公式和III型高压储氢容器的承压特点，列如下平衡方程式：
a、轴向和环向内力平衡方程式为：
Figure BDA0002744660840000021
Figure BDA0002744660840000022
b、极限压力下平衡受力平衡方程式为：
Figure BDA0002744660840000023
σdtαsin2α+σdtθ＝λRPm；
c、预应力下在极限压力的受力平衡方程式为：
Figure BDA0002744660840000024
σdtαsin2α+σdtθ+σipti＝RPm；
以上各方程式中参数的含义如下所示：
d-高压储氢容器的内直径；
R-高压储氢容器的内半径；
P-高压储氢容器承受的任意压力(P＜Pm)；
Pm-高压储氢容器承受的极限压力；
Nl-薄膜的轴向内力；
Nθ-薄膜的环向内力；
σd-纤维复合材料的设计强度；
σip-内衬材料的弹性极限；
tα-纤维复合层的螺旋缠绕纤维厚度；
tθ-纤维复合层的环向缠绕纤维厚度；
ti-内衬层的内衬厚度；
α-纤维复合层螺旋缠绕的缠绕角；
λ-III型高压储氢容器的纤维强度利用系数；
t-t＝tα+tθ+ti；
ri-高压储氢容器的任意半径(R＜ri＜R+t)；
(3)根据轴向和环向内力平衡的平衡方程式、极限压力下平衡受力平衡的平衡方程式、预应力下在极限压力的受力平衡方程式求解得到：纤维复合层的螺旋缠绕纤维厚度tα、纤维复合层的环向缠绕纤维厚度tθ、内衬层的内衬厚度ti；
Figure BDA0002744660840000031
Figure BDA0002744660840000032
Figure BDA0002744660840000033
(4)求解纤维复合层的纤维预应力值σp：根据变形协调条件εθ＝εiθ，其中，εθ是纤维的环向应变，εiθ是内衬的环向应变；联立步骤3中(2)和(3)的公式，求得纤维复合层的纤维预应力值σp；
Figure BDA0002744660840000041
以上各方程式中参数的含义如下所示：
E0-内衬材料的弹性模量；
EG-纤维复合材料的弹性模量；
μ-内衬材料的泊松比；
σip-内衬材料的弹性极限；
σiθ-内衬材料的环向应力；
Pw-高压储氢容器承受的工作压力；
K-与材料中纤维和树脂比重相关的量；
步骤4：利用步骤2中纤维复合层及内衬层的材料性能参数，验证步骤3中计算结果：tα、tθ、ti、σp的准确性；
(1)高压储氢容器筒身段的应变矩阵为：
{ε}＝[Sij]{σ1}；
(2)根据高压储氢容器筒身段的应变矩阵和步骤2中纤维复合层及内衬层的材料性能参数，计算高压储氢容器筒身段的整体刚度矩阵；
应力应变的转换矩阵
Figure BDA0002744660840000042
其中
Figure BDA0002744660840000043
那么转换坐标系后的应变矩阵为：
{εx}＝[Tε]-1[Sij][Tσ]{σx}；
求出柔度矩阵：
Figure BDA0002744660840000051
得到整体刚度矩阵：
Figure BDA0002744660840000052
考虑三项应力状态，应力与应变的关系：
Figure BDA0002744660840000053
(3)计算应力分布；
a、列出平衡方程、几何方程和物理方程：
Figure BDA0002744660840000054
Figure BDA0002744660840000055
{εr}＝[Sij]{σr}；
b、列出关于应力的常微分方程
Figure BDA0002744660840000056
其中：
Figure BDA0002744660840000057
Figure BDA0002744660840000058
Figure BDA0002744660840000059
Figure BDA00027446608400000510
Figure BDA00027446608400000511
Figure BDA0002744660840000061
c、根据上述应力的常微分方程的通解并利用边界条件：
σr|r＝b＝-P0，σr|r＝d＝0，求解待定系数，最终得到纤维复合层的环向应力
和纤维复合层的径向应力关于半径r的计算公式：
Figure BDA0002744660840000062
Figure BDA0002744660840000063
其中：
Figure BDA0002744660840000064
Figure BDA0002744660840000065
G3＝ζk2；
Figure BDA0002744660840000066
Figure BDA0002744660840000067
Figure BDA0002744660840000068
Figure BDA0002744660840000069
同理，内衬层在高压储氢容器筒身段的环向应力和径向应力关于r的计算公式为：
σir＝[Pw(R2d2-R2b2)+P0(R2b2-b2r2)]·[r2(R2-b2)]-1；
σiθ＝[Pw(R2r2-R2b2)+P0(R2b2-b2r2)]·[r2(R2-b2)]-1；
以上各方程式中参数的含义如下所示：
Pw-高压储氢容器承受的工作压力；
P0-外套对内衬的约束压力；
R-高压储氢容器的内半径；
b-内衬层的外半径；
r-高压储氢容器的任意半径(a＜r＜d)；
σθ-纤维复合层的环向应力；
σr-纤维复合层的径向应力；
σir-内衬层的环向应力；
σiθ-内衬层的径向应力；
(4)根据Tsai-Wu强度准则判断高压储氢容器筒身段是否失效：若根据Tsai-Wu强度准则判断为不失效，则步骤3的(1)中确定的设计参数λ准确，求解得到的纤维复合层的螺旋缠绕纤维厚度tα、纤维复合层的环向缠绕纤维厚度tθ、内衬层的内衬厚度ti、纤维复合层的纤维预应力值σp数值准确；若根据Tsai-Wu强度准则判断为失效，则重新修正确定设计参数λ，重复步骤3和步骤4直至根据Tsai-Wu强度准则判断为不失效；
Tsai-Wu强度准则：
Fiσi+Fijσiσj≤1，(i，j＝1，2，3)；
根据实验可以确定各项强度系数，其中：
Figure BDA0002744660840000071
Figure BDA0002744660840000072
Figure BDA0002744660840000073
Figure BDA0002744660840000074
Figure BDA0002744660840000081
当满足Fiσi+Fijσiσj＞1时，就判断为失效；
以上各方程式中参数的含义如下所示：
Xt-纤维复合材料在X方向上的抗拉强度；
Xc-纤维复合材料在X方向上的抗压强度；
Yt-纤维复合材料在Y方向上的抗拉强度；
Yc-纤维复合材料在Y方向上的抗拉强度。
进一步地，前述的一种基于纤维缠绕工艺的高压储氢容器的设计方法，其中，步骤3的(1)中的设计准则为：当高压储氢容器达到极限压力时纤维达到其强度极限；高压储氢容器达到工作压力时内衬材料达到弹性极限。
进一步地，前述的一种基于纤维缠绕工艺的高压储氢容器的设计方法，其中，步骤4中转换坐标系后的应变矩阵：
{εx}＝[Tε]-1[Sij][Tσ]{σx}；转换后的坐标系指以高压储氢容器的筒身轴线方向为x轴，高压储氢容器的圆筒径向方向为z轴，y为xz所成平面的法向；转换前x轴沿纤维方向，其他两轴与转换后描述相同。
本发明的有益效果是：将纤维的应力校核与材料失效判断引入到高压储氢容器的设计之中，对计算出的纤维复合层的缠绕厚度值等计算结果进行校核和修正，可以提升各计算结果的准确度，为高压储氢容器的设计制造、结构优化提供了理论依据，进一步降低采用纤维增强复合材料制成的纤维缠绕式高压储氢容器的制造成本。
附图说明
图1本发明所述的一种基于纤维缠绕工艺的高压储氢容器的设计方法的原理图。
图2是使用MATLAB对环向应力分布的模拟结果图。
图3是使用MATLAB对径向应力分布的模拟结果图。
图4是高压储氢容器的内部结构示意图。
图5是图4的局部结构示意图。
图6是图4中A-A剖视方向的局部结构示意图。
具体实施方式
下面结合附图及优选实施例对本发明所述的技术方案作进一步详细的说明。
如图1所示，本发明所述的一种基于纤维缠绕工艺的高压储氢容器的设计方法，所述的设计方法具体包括如下步骤：
步骤1：如图4、图5和图6所示，确定高压储氢容器的基本参数。
本发明采用50MPa储氢瓶，破裂压力按照SAE J2579的2.25倍安全系数，为112.5MPa。
综上，高压储氢容器的基本参数为：高压储氢容器承受的工作压力Pw为50MPa、高压储氢容器承受的极限压力Pm为112.5MPa、高压储氢容器的容积V为1000L、高压储氢容器的筒身长L为2.11m、高压储氢容器的内半径R为0.37m、纤维复合层螺旋缠绕的缠绕角α设定为40度。
步骤2：如图4、图5和图6所示，获取纤维复合层2及内衬层1的材料性能参数。
本发明以东丽公司T700S为例，其抗拉强度4900MPa，抗拉模量230GPa。内衬选用6061铝合金，弯曲屈服强度255MPa，弹性模量69.7GPa。泊松比0.33。环氧树脂基碳纤维复合材料的参数为σb＝2678.8Mpa，E1＝181Gpa，E2＝E3＝10.3Gpa，G12＝7.17Gpa，v＝0.28。
综上，纤维复合材料的抗拉强度σd为4900MPa、纤维复合材料的抗拉模量EF为230GPa、内衬材料的弯曲屈服强度σip为255MPa、内衬材料的弹性模量E0为69.7GPa、内衬材料的泊松比μ为0.33。
步骤3：计算纤维复合层2的缠绕厚度值和内衬层1的内衬厚度值：
(1)根据步骤1中高压储氢容器的基本参数和步骤2中纤维复合层及内衬层的材料性能参数确定设计准则和设计参数λ；其中，λ为III型高压储氢容器的纤维强度利用系数；这里设计参数λ设定为0.95。
其中，设计准则为：当高压储氢容器达到极限压力时纤维达到其强度极限；高压储氢容器达到工作压力时内衬材料达到弹性极限。
(2)根据高压储氢容器在极限压力下的平衡公式和III型高压储氢容器的承压特点，列如下平衡方程式：
a、轴向和环向内力平衡方程式为：
Figure BDA0002744660840000101
Figure BDA0002744660840000102
b、极限压力下平衡受力平衡方程式为：
Figure BDA0002744660840000103
σdtαsin2α+σdtθ＝λRPm；
c、预应力下在极限压力的受力平衡方程式为：
Figure BDA0002744660840000104
σdtαsin2α+σdtθ+σipti＝RPm；
以上各方程式中参数的含义如下所示：
d-高压储氢容器的内直径；
R-高压储氢容器的内半径；
P-高压储氢容器承受的任意压力(P＜Pm)；
Pm-高压储氢容器承受的极限压力；
Nl-薄膜的轴向内力；
Nθ-薄膜的环向内力；
σd-纤维复合材料的设计强度；
σip-内衬材料的弹性极限；
tα-纤维复合层的螺旋缠绕纤维厚度；
tθ-纤维复合层的环向缠绕纤维厚度；
ti-内衬层的内衬厚度；
α-纤维复合层螺旋缠绕的缠绕角；
λ-III型高压储氢容器的纤维强度利用系数；
t-t＝tα+tθ+ti；
ri-高压储氢容器的任意半径(R＜ri＜R+t)。
(3)根据轴向和环向内力平衡的平衡方程式、极限压力下平衡受力平衡的平衡方程式、预应力下在极限压力的受力平衡方程式求解得到：纤维复合层的螺旋缠绕纤维厚度tα、纤维复合层的环向缠绕纤维厚度tθ、内衬层的内衬厚度ti；
Figure BDA0002744660840000111
Figure BDA0002744660840000112
Figure BDA0002744660840000121
将各已知参数代入上述公式中，最终得到如下计算结果：
tα＝6.89mm；
tθ＝5.24mm；
ti＝8.08mm。
(4)求解纤维复合层的纤维预应力值σp：根据变形协调条件εθ＝εiθ，其中，εθ是纤维的环向应变，εiθ是内衬的环向应变；联立步骤3中(2)和(3)的公式，求得纤维复合层的纤维预应力值σp；
Figure BDA0002744660840000122
以上各方程式中参数的含义如下所示：
E0-内衬材料的弹性模量；
EG-纤维复合材料的弹性模量；
μ-内衬材料的泊松比；
σip-内衬材料的弹性极限；
σiθ-内衬材料的环向应力；
Pw-高压储氢容器承受的工作压力；
K-与材料中纤维和树脂比重相关的量；
将各已知参数代入上述公式中，最终得到如下计算结果：σp＝342.57MPa。
步骤4：利用步骤2中纤维复合层及内衬层的材料性能参数，验证步骤3中计算结果：tα、tθ、ti、σp的准确性：
(1)高压储氢容器筒身段的应变矩阵为：
{ε}＝[Sij]{σ1}。
(2)根据高压储氢容器筒身段的应变矩阵和步骤2中纤维复合层及内衬层的材料性能参数，计算高压储氢容器筒身段的整体刚度矩阵；
应力应变的转换矩阵
Figure BDA0002744660840000131
其中
Figure BDA0002744660840000132
那么转换坐标系后的应变矩阵为：
{εx}＝[Tε]-1[Sij][Tσ]{σx}；
转换后的坐标系指以高压储氢容器的筒身轴线方向为x轴，高压储氢容器的圆筒径向方向为z轴，y为xz所成平面的法向；转换前x轴沿纤维方向，其他两轴与转换后描述相同。
求出柔度矩阵：
Figure BDA0002744660840000133
得到整体刚度矩阵：
Figure BDA0002744660840000134
考虑三项应力状态，应力与应变的关系：
Figure BDA0002744660840000135
(3)计算应力分布：
a、列出平衡方程、几何方程和物理方程：
Figure BDA0002744660840000136
Figure BDA0002744660840000137
{εr}＝[Sij]{σr}；
b、列出关于应力的常微分方程
Figure BDA0002744660840000141
其中：
Figure BDA0002744660840000142
Figure BDA0002744660840000143
Figure BDA0002744660840000144
Figure BDA0002744660840000145
Figure BDA0002744660840000146
Figure BDA0002744660840000147
c、根据上述应力的常微分方程的通解并利用边界条件σr|r＝b＝-P0，σr|r＝d＝0，求解待定系数，最终得到纤维复合层的环向应力和纤维复合层的径向应力关于半径r的计算公式：
Figure BDA0002744660840000148
Figure BDA0002744660840000149
其中：
Figure BDA00027446608400001410
Figure BDA0002744660840000151
G3＝ζk2；
Figure BDA0002744660840000152
Figure BDA0002744660840000153
Figure BDA0002744660840000154
Figure BDA0002744660840000155
求解待定系数主要是数值解法，本领域人员可自行推导，属于本领域常规解法，因而这里不再赘述。
同理，内衬层在高压储氢容器筒身段的环向应力和径向应力关于r的计算公式为：
σir＝[Pw(R2d2-R2b2)+P0(R2b2-b2r2)]·[r2(R2-b2)]-1；
σiθ＝[Pw(R2r2-R2b2)+P0(R2b2-b2r2)]·[r2(R2-b2)]-1；
以上各方程式中参数的含义如下所示：
Pw-高压储氢容器承受的工作压力；
P0-外套对内衬的约束压力；
R-高压储氢容器的内半径；
b-内衬层的外半径；
r-高压储氢容器的任意半径(a＜r＜d)；
σθ-纤维复合层的环向应力；
σr-纤维复合层的径向应力；
σir-内衬层的环向应力；
σiθ-内衬层的径向应力；
图2是使用MATLAB对环向应力分布的模拟结果图，图3是使用MATLAB对径向应力分布的模拟结果图。从图2和图3中可以看出环向应力和径向应力的变化趋势：纤维复合层都是随径向位移的增大而减小，而内衬层的环向应力σir的值则在400MPa附近。且在内衬和纤维接触面处均发生应力突变，从图2和图3中可以看出纤维复合层承受更大的环向应力和径向应力。
(4)根据Tsai-Wu强度准则判断高压储氢容器筒身段是否失效：若根据Tsai-Wu强度准则判断为不失效，则步骤3的(1)中确定的设计参数λ准确，求解得到的纤维复合层的螺旋缠绕纤维厚度tα、纤维复合层的环向缠绕纤维厚度tθ、内衬层的内衬厚度ti、纤维复合层的纤维预应力值σp数值准确；若根据Tsai-Wu强度准则判断为失效，则重新修正确定设计参数λ，重复步骤3和步骤4直至根据Tsai-Wu强度准则判断为不失效；
Tsai-Wu失效判据式：
Fiσi+Fijσiσj≤1，(i，j＝1，2，3)；
Xt＝1760Mpa；Xc＝1016.63Mpa；Yt＝25Mpa；Yc＝117Mpa；根据实验可以确定各项强度系数，其中：
Figure BDA0002744660840000161
Figure BDA0002744660840000162
Figure BDA0002744660840000163
Figure BDA0002744660840000164
Figure BDA0002744660840000165
当满足Fiσi+Fijσiσj＞1时，就判断为失效；
以上各方程式中参数的含义如下所示：
Xt-纤维复合材料在X方向上的抗拉强度；
Xc-纤维复合材料在X方向上的抗压强度；
Yt-纤维复合材料在Y方向上的抗拉强度；
Yc-纤维复合材料在Y方向上的抗拉强度；
将各已知参数代入上述公式中，最终得到如下计算结果：代入Tsai-Wu后的系数是0.974，满足Tsai-Wu失效判据式，气瓶结构安全。
以上所述仅是本发明的较佳实施例，并非是对本发明作任何其他形式的限制，而依据本发明的技术实质所作的任何修改或等同变化，仍属于本发明要求保护的范围。
本发明的优点是：将纤维的应力校核与材料失效判断引入到高压储氢容器的设计之中，对计算出的纤维复合层的缠绕厚度值等计算结果进行校核和修正，可以提升各计算结果的准确度，为高压储氢容器的设计制造、结构优化提供了理论依据，进一步降低采用纤维增强复合材料制成的纤维缠绕式高压储氢容器的制造成本。

Claims (3)

1.一种基于纤维缠绕工艺的高压储氢容器的设计方法，其特征在于：所述的设计方法具体包括如下步骤：
步骤1：确定高压储氢容器的基本参数：高压储氢容器承受的工作压力Pw、高压储氢容器承受的极限压力Pm、高压储氢容器的容积V、高压储氢容器的筒身长L、高压储氢容器的内半径R、纤维复合层螺旋缠绕的缠绕角α；
步骤2：获取纤维复合层及内衬层的材料性能参数：纤维复合材料的抗拉强度σd、纤维复合材料的抗拉模量EF、内衬材料的弯曲屈服强度σip、内衬材料的弹性模量E0、内衬材料的泊松比μ；
步骤3：计算纤维复合层的缠绕厚度值和内衬层的内衬厚度值；
(1)根据步骤1中高压储氢容器的基本参数和步骤2中纤维复合层及内衬层的材料性能参数确定设计准则和设计参数λ；其中，λ为Ⅲ型高压储氢容器的纤维强度利用系数；
(2)根据高压储氢容器在极限压力下的平衡公式和Ⅲ型高压储氢容器的承压特点，列如下平衡方程式：
a、轴向和环向内力平衡方程式为：
Figure FDA0002744660830000011
Figure FDA0002744660830000012
b、极限压力下平衡受力平衡方程式为：
Figure FDA0002744660830000013
σdtαsin2α+σdtθ＝λRPm；
c、预应力下在极限压力的受力平衡方程式为：
Figure FDA0002744660830000014
σdtαsin2α+σdtθ+σipti＝RPm；
以上各方程式中参数的含义如下所示：
d-高压储氢容器的内直径；
R-高压储氢容器的内半径；
P-高压储氢容器承受的任意压力(P<Pm)；
Pm-高压储氢容器承受的极限压力；
Nl-薄膜的轴向内力；
Nθ-薄膜的环向内力；
σd-纤维复合材料的设计强度；
σip-内衬材料的弹性极限；
tα-纤维复合层的螺旋缠绕纤维厚度；
tθ-纤维复合层的环向缠绕纤维厚度；
ti-内衬层的内衬厚度；
α-纤维复合层螺旋缠绕的缠绕角；
λ-Ⅲ型高压储氢容器的纤维强度利用系数；
t-t＝tα+tθ+ti；
ri-高压储氢容器的任意半径(R＜ri＜R+t)；
(3)根据轴向和环向内力平衡的平衡方程式、极限压力下平衡受力平衡的平衡方程式、预应力下在极限压力的受力平衡方程式求解得到：纤维复合层的螺旋缠绕纤维厚度tα、纤维复合层的环向缠绕纤维厚度tθ、内衬层的内衬厚度ti；
Figure FDA0002744660830000021
Figure FDA0002744660830000022
Figure FDA0002744660830000023
(4)求解纤维复合层的纤维预应力值σp：根据变形协调条件εθ＝εiθ，其中，εθ是纤维的环向应变，εiθ是内衬的环向应变；联立步骤3中(2)和(3)的公式，求得纤维复合层的纤维预应力值σp；
Figure FDA0002744660830000031
以上各方程式中参数的含义如下所示：
E0-内衬材料的弹性模量；
EG-纤维复合材料的弹性模量；
μ-内衬材料的泊松比；
σip-内衬材料的弹性极限；
σiθ-内衬材料的环向应力；
Pw-高压储氢容器承受的工作压力；
K-与材料中纤维和树脂比重相关的量；
步骤4：利用步骤2中纤维复合层及内衬层的材料性能参数，验证步骤3中计算结果：tα、tθ、ti、σp的准确性；
(1)高压储氢容器筒身段的应变矩阵为：
{ε}＝]Sij]{σ1}；
(2)根据高压储氢容器筒身段的应变矩阵和步骤2中纤维复合层及内衬层的材料性能参数，计算高压储氢容器筒身段的整体刚度矩阵；
应力应变的转换矩阵
Figure FDA0002744660830000041
其中
Figure FDA0002744660830000042
那么转换坐标系后的应变矩阵为：
{εx}＝[Tε]-1[Sij][Tσ]{σx}；
求出柔度矩阵：
Figure FDA0002744660830000043
得到整体刚度矩阵：
Figure FDA0002744660830000044
考虑三项应力状态，应力与应变的关系：
Figure FDA0002744660830000045
(3)计算应力分布；
a、列出平衡方程、几何方程和物理方程:
Figure FDA0002744660830000046
Figure FDA0002744660830000047
{εr}＝[Sij[{σr}；
b、列出关于应力的常微分方程
Figure FDA0002744660830000048
其中：
Figure FDA0002744660830000051
Figure FDA0002744660830000052
Figure FDA0002744660830000053
Figure FDA0002744660830000054
Figure FDA0002744660830000055
Figure FDA0002744660830000056
c、根据上述应力的常微分方程的通解并利用边界条件：
σr|r＝b＝-P0,σr|r＝d＝0，求解待定系数，最终得到纤维复合层的环向应力和纤维复合层的径向应力关于半径r的计算公式：
Figure FDA0002744660830000057
Figure FDA0002744660830000058
其中：
Figure FDA0002744660830000059
Figure FDA00027446608300000510
G3＝ζk2；
Figure FDA00027446608300000511
Figure FDA00027446608300000512
Figure FDA00027446608300000513
Figure FDA0002744660830000061
同理，内衬层在高压储氢容器筒身段的环向应力和径向应力关于r的计算公式为：
σir＝[Pw(R2d2-R2b2)+P0(R2b2-b2r2)]·[r2(R2-b2)]-1；
σiθ＝[Pw(R2r2-R2b2)+P0(R2b2-b2r2)]·[r2(R2-b2)]-1；
以上各方程式中参数的含义如下所示：
Pw-高压储氢容器承受的工作压力；
P0-外套对内衬的约束压力；
R-高压储氢容器的内半径；
b-内衬层的外半径；
r-高压储氢容器容器的任意半径(a<r<d)；
σθ-纤维复合层的环向应力；
σr-纤维复合层的径向应力；
σir-内衬层的环向应力；
σiθ-内衬层的径向应力；
(4)根据Tsai-Wu强度准则判断高压储氢容器筒身段是否失效：若根据Tsai-Wu强度准则判断为不失效，则步骤3的(1)中确定的设计参数λ准确,求解得到的纤维复合层的螺旋缠绕纤维厚度tα、纤维复合层的环向缠绕纤维厚度tθ、内衬层的内衬厚度ti、纤维复合层的纤维预应力值σp数值准确；若根据Tsai-Wu强度准则判断为失效，则重新修正确定设计参数λ，重复步骤3和步骤4直至根据Tsai-Wu强度准则判断为不失效；
Tsai-Wu强度准则：
Fiσi+Fijσiσj≤1，(i，j＝1，2，3)；
根据实验可以确定各项强度系数，其中：
Figure FDA0002744660830000071
Figure FDA0002744660830000072
Figure FDA0002744660830000073
Figure FDA0002744660830000074
Figure FDA0002744660830000075
当满足Fiσi+Fijσiσj＞1时，就判断为失效；
以上各方程式中参数的含义如下所示：
Xt-纤维复合材料在X方向上的抗拉强度；
Xc-纤维复合材料在X方向上的抗压强度；
Yt-纤维复合材料在Y方向上的抗拉强度；
Yc-纤维复合材料在Y方向上的抗拉强度。
2.根据权利要求1所述的一种基于纤维缠绕工艺的高压储氢容器的设计方法，其特征在于：步骤3的(1)中的设计准则为：当高压储氢容器达到极限压力时纤维达到其强度极限；高压储氢容器达到工作压力时内衬材料达到弹性极限。
3.根据权利要求1或2所述的一种基于纤维缠绕工艺的高压储氢容器的设计方法，其特征在于：步骤4中转换坐标系后的应变矩阵：
{εx}＝[Tε]-1[Sij][Tσ]{σx}；转换后的坐标系指以高压储氢容器的筒身轴线方向为x轴，高压储氢容器的圆筒径向方向为z轴，y为xz所成平面的法向；转换前x轴沿纤维方向，其他两轴与转换后描述相同。
CN202011162433.9A 2020-10-27 2020-10-27 一种基于纤维缠绕工艺的高压储氢容器的设计方法 Active CN112460464B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011162433.9A CN112460464B (zh) 2020-10-27 2020-10-27 一种基于纤维缠绕工艺的高压储氢容器的设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011162433.9A CN112460464B (zh) 2020-10-27 2020-10-27 一种基于纤维缠绕工艺的高压储氢容器的设计方法

Publications (2)

Publication Number Publication Date
CN112460464A true CN112460464A (zh) 2021-03-09
CN112460464B CN112460464B (zh) 2021-12-21

Family

ID=74835523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011162433.9A Active CN112460464B (zh) 2020-10-27 2020-10-27 一种基于纤维缠绕工艺的高压储氢容器的设计方法

Country Status (1)

Country Link
CN (1) CN112460464B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116100841A (zh) * 2023-01-13 2023-05-12 江苏集萃复合材料装备研究所有限公司 一种用于高压氢气瓶的制备方法
CN114183684B (zh) * 2021-11-26 2023-06-09 江苏国富氢能技术装备股份有限公司 一种碳纤维全缠绕气瓶的增材补强及碳纤维缠绕方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5012950A (en) * 1988-04-30 1991-05-07 Holger Knappe Plastic container for liquids or gases
CN101029707A (zh) * 2007-03-21 2007-09-05 重庆汽车研究所 车用钢内衬纤维材料环向缠绕层增强cng气瓶和制造方法
CN102135178A (zh) * 2010-12-30 2011-07-27 西安航天复合材料研究所 一种用于压力容器的干纱缠绕成型方法
JP2013173345A (ja) * 2012-01-24 2013-09-05 Toyota Motor Corp 高圧ガスタンクの製造方法および製造装置
CN103672388A (zh) * 2013-11-29 2014-03-26 北京宇航系统工程研究所 一种封头-筒身一体化纤维缠绕复合材料气瓶的设计方法
CN105443973A (zh) * 2015-11-09 2016-03-30 大连理工大学 薄壁金属内衬纤维复合材料缠绕压力容器及其制造工艺
CN106017554A (zh) * 2016-05-18 2016-10-12 浙江大学 同时测量纤维缠绕金属内衬压力容器的温度和压力的装置
US20180045371A1 (en) * 2016-08-09 2018-02-15 Hyundai Motor Company High pressure tank
DE102018005571A1 (de) * 2018-07-13 2020-01-16 Daimler Ag Verfahren zur Herstellung eines Druckgasbehälters

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5012950A (en) * 1988-04-30 1991-05-07 Holger Knappe Plastic container for liquids or gases
CN101029707A (zh) * 2007-03-21 2007-09-05 重庆汽车研究所 车用钢内衬纤维材料环向缠绕层增强cng气瓶和制造方法
CN102135178A (zh) * 2010-12-30 2011-07-27 西安航天复合材料研究所 一种用于压力容器的干纱缠绕成型方法
JP2013173345A (ja) * 2012-01-24 2013-09-05 Toyota Motor Corp 高圧ガスタンクの製造方法および製造装置
CN103672388A (zh) * 2013-11-29 2014-03-26 北京宇航系统工程研究所 一种封头-筒身一体化纤维缠绕复合材料气瓶的设计方法
CN105443973A (zh) * 2015-11-09 2016-03-30 大连理工大学 薄壁金属内衬纤维复合材料缠绕压力容器及其制造工艺
CN106017554A (zh) * 2016-05-18 2016-10-12 浙江大学 同时测量纤维缠绕金属内衬压力容器的温度和压力的装置
US20180045371A1 (en) * 2016-08-09 2018-02-15 Hyundai Motor Company High pressure tank
DE102018005571A1 (de) * 2018-07-13 2020-01-16 Daimler Ag Verfahren zur Herstellung eines Druckgasbehälters

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114183684B (zh) * 2021-11-26 2023-06-09 江苏国富氢能技术装备股份有限公司 一种碳纤维全缠绕气瓶的增材补强及碳纤维缠绕方法
CN116100841A (zh) * 2023-01-13 2023-05-12 江苏集萃复合材料装备研究所有限公司 一种用于高压氢气瓶的制备方法
CN116100841B (zh) * 2023-01-13 2024-01-26 江苏集萃复合材料装备研究所有限公司 一种用于高压氢气瓶的制备方法

Also Published As

Publication number Publication date
CN112460464B (zh) 2021-12-21

Similar Documents

Publication Publication Date Title
CN112460464B (zh) 一种基于纤维缠绕工艺的高压储氢容器的设计方法
Zu et al. Design and analysis of filament-wound composite pressure vessels based on non-geodesic winding
CN113111517A (zh) 一种非等极孔纤维缠绕压力容器的设计方法
CN113154247B (zh) 一种高压储氢气瓶的制备方法
CN111368439A (zh) 一种基于缠绕成型工艺的压力容器的设计方法
CN103672388B (zh) 一种封头-筒身一体化纤维缠绕复合材料气瓶的设计方法
CN112297464B (zh) 一种基于纤维缠绕工艺的深冷高压储氢容器的设计方法
CN112963722A (zh) 一种金属内衬碳纤维缠绕复合材料气瓶及其制造方法
CN112329301B (zh) 一种适用于金属内衬复合材料气瓶的最优自紧压力确定方法
CN111998220A (zh) 一种高压复合轻量化储氢瓶
US6325108B1 (en) Prestressed composite cryogenic piping
Kang et al. Parametric optimization for multi-layered filament-wound cylinder based on hybrid method of GA-PSO coupled with local sensitivity analysis
Zu et al. Application of variable slippage coefficients to the design of filament wound toroidal pressure vessels
CN210088449U (zh) 一种高压复合轻量化储氢瓶
Coskun et al. Comparison of mechanical properties of the Type 1 and Type 2 composite hydraulic cylinder designs: A numerical study
Akcay et al. Analysis of multilayered composite cylinders under thermal loading
Özaslan et al. Design and validation of a filament wound composite rocket motor case
Karimi et al. Designing and optimizing of composite and hybrid drive shafts based on the bees algorithm
CN113962096B (zh) 一种导弹发射筒三向应力确定方法及系统
Razin The Problem of Optimum Design of Composite Housings of Solid Propellant Rocket Engines
Bae et al. Optimal design for compressed natural gas composite vessel by using coupled model with liner and composite layer
Srividya et al. Failure analysis of GFRP and CFRP composite laminated pressure vessel
CN206399283U (zh) 含金属内衬碳纤维复合材料身管
Zu et al. Integrated Design and Production of Filament-Wound Composite Structures: Compromise between Strength and Manufacturability
Elgohary et al. Stress analysis of multi-layered composite cylinders subjected to various loadings

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
GR01 Patent grant
GR01 Patent grant