型高压储罐优化设计方案方法.docx

1、型高压储罐优化设计方案方法3型高压储罐优化设计方法P. Francescato，A. Gillet，D. Leh，P. Saffr波尔多大学，法国国家科学研究中心摘要本研究承受内压下比较不同的优化设计方法生产的两个角缠绕结构。近些年，研究集中在材料和他们破坏行为以及绕组优化技术在船舶几何关系和制造过程，但这些研究没有建立明确的设计规则。在本文中，两种方法，一种是基于一个基本的CPT和渐进破坏和枚举优化，另一个更复杂的有限元优化相比最近结果基于清洁发展机制的理论和一个3 d有限元技术将一个改进遗传算法。结果被用来比较的优点每种方法的局限，并推导出简单的设计规则和计算策略，将可靠地获得一个很好的C

2、PU成本/效益比。关键词： 储氢罐； 缠绕结构； 高压储罐.1.介绍因为他们的高剪裁功率，和他们的非常有趣的与常规材料相比，复合材料的特定的属性材料越来越多地用在高性能的域，例如为航空航天，汽车和船舶等行业。特别地，在过去的十年中，大量的研究已经做了各种方法设计纤维缠绕储氢罐。但是尚未找到明确的解决方案进行氢车辆在绝对安全的条件下，在可接受的成本。最流行的存储解决方案，包括高压的筒形罐压力高达35兆帕存储）与制造纤维缠绕技术，采用高强度UD碳纤维/环氧树脂复合材料的金属3型船）或塑料类型4艘）眼线。长丝缠绕是一个过程，得到最佳的纤维/基体的比例，因此，最佳的刚性/重量比在长丝卷绕过程中，纤维束

3、被放置在一个旋转的可拆卸的芯棒和纤维的取向范围内是纤维滑移限制。高压储罐和一般轴对称结构是典型的类型的复合部件它可以是高度根据他们必须承受的负荷进行优化，幸好可能性的纤维取向的纤维缠绕的过程提供。优化的重量复合贮氢罐在爆破压力下是一个很大的挑战，因为最佳的复合物取决于物理结构和机械性能的材料，形状，制造，处理以及经济和环境的制约。鉴于高的参数数目，定义一个最优配置是复杂的，一般需要计算机辅助设计工具。已经做了相当多的研究，以了解和发展更精确和有效的计算预测的最终强度的复合材料结构的方法，包括具体的故障损害的法律和优化技术。因此，任何介绍，这些作品是一定的限制，所以，我们只保留了最近的研究中的特

4、定复合材料压力容器的优化设计问题。第一个研究了专门对最后的关闭和圆顶进行了优化设计这是主要的和最复杂的部分压力容器的设计。一般理论假设是简单的方法来考虑的结构特性，只从一个容器设计与测地线的圆顶连续的安排，只是体现出纯张力。霍贾提出 1首先显示的属性树脂的形状的设计中具有相当的重要性2使用增量的轴对称的方法在圆顶计算绕组的模式，生产复合材料铺层的厚度不变3分析方程和图形分析方法，设计层压复合材料压力容器下的应变和强度的限制。利用经典层合板理论和广义平面应力模型，Parnas和Katirci开发了一个分析程序设计和复合材料压力预测的行为4在各种负载条件下的船舶。这项研究是能够建立对于一个给定的机

5、械加载组合的最优卷绕角。瓦西里耶夫等。 5提出了一种优化的圆顶形状与最大的可靠性和合理的质量纤维胶带纵向和横向钢筋。 Kim等人。6建立了纤维缠绕的优化内部压力下的结构的设计方法。一个半测地路径算法被用来确定可能的绕组图案，以滑移考虑。 Kim等人。73型储罐和优化工作谐振边缘设计，使用半测地路径算法，修改遗传算法和渐进失效分析，找到最好的铝 - 碳纤维/环氧树脂组合制造对于一个给定的内部压力的最轻的储罐。在这项工作中，进行渐进失效分析的方法。8由刘和郑提出了能源为基础的损伤模型损伤力学理论。9预测一个简单的铝 - 碳纤维/环氧渐进破坏复合材料层合板3型储罐内部压力增加。模拟的三维有限元模型的

6、建立，考虑到失败的层压板的损伤过程模式：纤维断裂，基体开裂和剪切破坏。后来，同一作者和他的同事使用了不同的顶部方法消耗臭氧层物质，以确定最佳的设计相同的复合储氢容器：自适应遗传算法10人工免疫系统的方法11。本文的目的是检查用不同的方法第一层的故障FPF）损害的进展和菲纳故障FF）为类型3储罐进行单向形式，单向增加内部压力和，以确定最佳的设计解决方案方面不同损伤的结果。优化设计是考虑一个下达到最小重量钻压力和材料的限制。通常的方法的实现这个目标，简要描述的，主要区别在于他们的方式和检测损坏的组成层数，采取有分的行为的板块数。众所周知，大部分退化模型分为三个主要类别Padhi等12，Spotts

7、wood和Palazotto13）一）总伤害的方法失败的帘布层的刚度和强度的减小到零，除了性状的纤维的方向，这种方法往往会导致两个步骤的过程等同于的抵销的理论，）逐步破坏方法：减少的刚性依赖于失败模式的帘布层和管理经验退化因素当层发生故障时，有一个应力再分配在机智层和其余的重复增量的分析通常在哪一层失败，在纤维方向的负荷，称为FF负载），三）连续损伤力学CDM）理论与数值计算程序有限元FE）方法消耗臭氧层物质），导致的损伤张量的精确定义和其演变规律，并已被证明是最现代的方法，获得损害开始和摄取的复合结构，尽管需要的大量的实验努力识别损伤参数和大量的CPU时间数值实现。在这篇文章中，不同的计算策

8、略的确定FPF和FF，所以3型储罐的优化设计，相比彼此之间并与最近的出版物中提出的结果，尤其是6,7两种方法。2.参考解决方案Kim等人。6,7，以执行两个版本的优化设计纤维缠绕型3储罐内部压力下。他们的优化设计方法，可能包括最好的程序此刻，那就是：一个半测地路径算法，一个进步故障分析CDM理论）和遗传算法GA）一个三维有限元方法的数字，尽管与高数值和实验成本。优化设计算法应用两种有代表性的高压氢存储坦克，在本文章中简称为T1和T2。这两个高压容器3型储罐组成的薄铝板由碳纤维增强塑料CFRP）的衬垫包裹由于高压氢内部承受的应力。储罐的结构的层压材料组成的三层图1a）：铝衬套，沿其整个长度的螺旋

9、层罐在其纵向轴线卷绕角度h，从和箍层仅在筒状部的储罐轴的槽）。的材料性质和几何和基本设计条件的这两个的Al-CFRP坦克总结分别在表1和表2中。 限制的研究领域应用的设计变量优化和遗传算法的变量在表3中。我们认为相同纤维体积分数和厚度每个回圈层，固定在0.55毫M，T1，T2和0.8毫M的。应用程序的目标是确定的四个最佳设计变量nh个螺旋状的层的数目，nh个箍层的数目，H的卷绕角度和el的衬套的厚度），满足于表2中的基本设计条件。所获得的结果为两个的问题，提出进一步一段与我们自己的结果。 第一种方法，M1：工程方法3.1.原理和假设 目的是提出解决一个简单的方法，以最小的努力和可接受的合理的精

10、度，设计问题对均匀内承受高压储存容器压力。所考虑的问题是，减少尽可能在突发压力复合容器的重量。由基姆等人得到的结果。【6】这些标准坦克证明在螺旋状层位于纵向纤维应力在圆柱部分大于位于除了在穹顶，穹顶的底部附近，被排除在外从我们的研究中图1A）。同时，它似乎是合理的假设高压罐的优化设计可以实现经典板理论CPT） 1 爱基尔霍夫，应用它的圆柱部分见 10 ）。本发明的方法进行下列因素和假设：研究只侧重于容器筒节受到内部压力P与近端的影响加载径向，环向和轴向坐标分别为用R，H和Z衬垫的厚度EL内半径R0是恒定的在圆柱体的整个长度容器的总厚度等于：e=e1 +e+eke是螺旋层的厚度，是从Z纵轴的缠绕

11、角，eh是环向层的厚度，两者都是常数的圆柱体的整个长度图1b） 厚度比/R0是小于1.1；在这种情况下，薄壁的分析可以满足 15 使圆柱形压力容器进行到处平面H，Z）纵向NZ 和横向NH合成图1C）：NH= PR0，NZ = N / 2尽管没有弯矩应用层压板的非对称性，和横向剪切变形忽略CPT）合金被认为是各向同性和理想弹塑性；它服从vonM塞斯准则UD的碳纤维/环氧复合材料被认为是线弹性、横观各向同性和符合故障类型的标准：最大应力，TsaiHill或Hashin。这些标准是详情在12的平面应力 第弗朗西斯卡托等人。/复合结构，942018）20872096图1。a）纵向纤维应力 6 ，b）从

12、6中研究了一个3型槽结构和c）圆筒状部分的有代表性的板在纵向的NZ和横向NH膜负载。表1一种铝衬和单向碳纤维复合材料本构特性 161.3 GPa 149.12 GPaE2, E3 8.82 GPa 10.558 GPaG12, G13 5.331 GPa 4.138 GPaG23 2.744 GPa 3.311 GPav12, v13 0.33 0.253v23 0.45 0.421Xt (在纤维方向的拉伸强度 2300 MPa 2408 MPaYt (在transversedirection抗拉强度30 MPa 37 MPaS (平面剪切强度 50 MPaa55 MPaa密度 1.58x10

13、-6kg/mm3 1.608x10-6kg/mm3 一个标准的值，因为在文献中没有提出。 6,7 。绕组和交织的影响是不考虑；只有UD平面堆叠序列被认为是在这些条件下，圆柱壳的一小部分可以分析，膜下加载板(LT是k层的破坏准则和S是安全的因素.3.2.方法损伤规律和优化在层压板失败的第一层，平面应力在不同层的罐的圆柱形部分可以准确计算程序1）也可用于检测负载的第一层出现故障FPF）。FPF一般是独立的失效准则见 12 和进一步的结果）。一般来说，FPF涉及某些类型的基质降解，下面将要看到的，可以被认为是可以容忍的层压直至最终失效FF）的结构，这通常是由于纤维断裂。该基质降解导致在横向剪切刚度E

14、T和GLT的帘布层但不大大影响的纵向刚度和强度的ZL和XT显著减少。当一个层失败，有剩余的完整的层恒载应力重分布图2）。基质的降解和其发展是在复合材料失效的一个非常重要的研究领域，它是为工程师了解和分析层压复合行为FPF FF为了评估损伤容限，因此，例如基础，复合材料结构的实际安全系数。对损坏的面内弹性法，横观各向同性k-ply，表达的是经典的 16 ：第弗朗西斯卡托等人。/复合结构，942018）20872096表2槽T1 T2几何描述和设计基础条件来自基姆等人。 6,7 ）。槽T1槽T2asteb对称槽工业型对称的3箱在气缸半径的几何127 mm 186.44 mm描述圆顶内半径45.8

15、mm船舶的总长度454.4 mm 823.37 mm基本设计1.最大工作压力13.79兆帕PW1.内部的最大工作压力为35 MPa PW条件2.衬垫的产量是禁止的2.纤维定向设计强度为878.2兆帕3.复合材料的安全系数3 3.衬垫的设计强度为123.2兆帕4.衬垫厚度被认为是恒定的4.实际爆破压力的压力的比值2.355.减轻重量是最重要的设计目标5.衬垫厚度被认为是恒定的 6.减轻重量是最重要的设计目标表3T1和T2对称3型槽的设计变量来自基姆等人 6,7 ）。槽T1槽T2asteb对称槽工业型对称的3箱螺旋层设计变量的最小值：2最大值：7最小值：2最大值：5个数箍层数最小值：1最大值：16

16、 最小值：1最大值：16螺旋和环向层1.10.55 1.60.8厚度mm）衬砌厚度最小值：1.0最大值：2.5 步骤0.1最小值：0.9最大值：4.0 步骤0.1mm）绕组可行12.0, 13.0, 14.0, 15.5, 16.5, 17.0, 6.5, 7.5, 8.0, 8.5, 9.0, 11.0, 12.0, 的角度18.0, 21.0, 21.5,22.5, 23.5, 24.0, 13.0, 13.5, 15.0, 16.0, 17.0,17.5,24.5, 25.0, 25.5,27.5, 28.0, 28.5,18.5, 19.0, 19.5, 20.0, 20.5, 21.

17、5, 22.0 29.5, 30.5, 31.5, 32.5, 33.0, 33.5,34.0, 34.5图2.对不同的损伤模型理论的内在压力应变曲线。PFPF：第一层破坏压力,dT(k和dLT(k，依赖于复杂的模型CDM型）和复杂的实验程序加载/卸载，看到 16,17 ），我们采用第一种方法，在实证退化因子的增量荷载分析是用于减少的k-ply根据确定的Hashin失效准则检测失效模式的材料特性请参阅巴布罗 14 ，阿克拉斯和李 18 ）.这种退化的因素通常是裂纹密度相关，是指国际光纤故障在帘布层K除以实验区的长度，乘以厚度诺普斯和BGLE 19 ）。裂纹密度的识别，因此降解因子，需要基于光学

18、或声学技术，大量的实验工作。在目前情况下，一个单一的任意矩阵退化的因素是用来减少退化层选择的刚度值如下0指出原完整）的财产，相当于=0和1是一个完全损坏的财产。为了简单而有效的计算，本模型是基于一个固定值的损伤模型，这是用来按一定关系减少退化层的刚度特性3）。一个渐进的过程，是选择在降解的因素是改变迭代直到预测FF负载对应，在目前的情况下，到第一纤维断裂的检测。这种方法会导致一个分段线性应力应变曲线黑色虚线见图2）和会，没有实验，允许的值在01）的影响的测定和关于CDM模型及其精度在图2中的黑色曲线）和后来的实验曲线。FF负载通过此方法计算可能会低于实际的断裂载荷，从而高估了铺层设计保守的方法

19、）；但是，在任何情况下，破坏过程的仔细的解释和产生的FF负载是必要的。它也有可能找到一个FF负载近似由一个简单的两步方法，其刚度和失败的铺层的强度降低到接近于零的值；本方法等效于网格理论，因此具有有限的适用性，虽然通过该方法得到的FF负载一般用渐进破坏法 14 一致图2）。最后，我们注意到，在图2中提出的不同的曲线代表的最终破坏类型对应的层合板强度控制由纤维(LT可以在一个自己写的分析程序，计算实现，plyby层，在平面应变和应力基于CPT。对于一个固定的目标函数和固定的基本设计条件如表2型），该程序可以很容易地从一个定义域的研究分类数控组合，找到最佳的解决方案。这可以在几秒钟内完成相应的研究

20、表3域的组合，一个简单的FORTRAN和MATLAB 程序。本文首次运用常规叫做FPF优化。这个程序更有趣的使用是在FF找到最佳的设计，需要考虑的损伤萌生和在加载增加其进展的帐户。在5节中描述的进一步应用，我们选择了Hashin失效准则具有可识别的故障模式和高精度平面随着刚度折减系数，从0.7到0.99不等的逐步退化的方法。这第二次使用的程序被称为FF的优化。在本文的其余部分，我们将这一工程设计方法为M1。 第弗朗西斯卡托等人。/复合结构，942018）20872096图3.一）变化的螺旋层衬板和纵向应力等效VonM塞斯应力；b）破坏准则的水平在厚度和沿T1的横截面的最优解W = 11.43公

21、斤/平方M问题P2）m2-fpf。4.第二种方法，平方M：有限元法和渐进损伤和遗传算法的优化4.1.原理和假设为了实现对约束的设计问题的表2中所给出的最优解，轴对称缠绕结构的优化设计程序的开发：一个半测地路径算法，Reddy的渐进失效分析和有限元进行了遗传算法分析。这个程序是由基姆等人使用相似的。在参考文献6,7。的M2的方法具体假设如下：内部的圆顶形状是固定的椭圆的一部分，尊重坦克T1和T2施加的几何条件。由于层的特殊结构，采用了一些二维轴对称四边形元素与2 X 2集成点，允许有限的CPU成本的应力状态很好的近似。网格是由基姆等人使用相同。：在圆柱部分20元素和60元素的圆顶区域内端径向位移

22、是固定的，(ii一个统一的轴面力，相应的在水箱底部压力的影响，应用于衬,(iii的平面对称条件。蜿蜒的路径是用半测地路径算法计算 20 ：在和是R相对于z第一和第二衍生物。是切向和法向力之间的比率，它必须保持低于滑移系数。k = 0，一个在clairault方程发现同样的趋势：rsina = CSTE。方程=0如果弹性，I(n=1的损坏或可塑性的实例。优化了遗传算法已成功地应用在最近的优化问题求解复合基础板和结构 21 。我们采用整数编码的个体，一个统一的均匀变异算子。约束条件考虑的原则为基础的可行个体优越：一个可行的个人而不是一个系统保留不可行。在随后的一代一代人口置换和程序的建立。这 个选

23、择过程是用来保存大量的母体高性能个人。通常情况下，东北最好的个人的母体是再抄写在儿童人口。失踪的个体是由再生产。 22 表明，该技术可以实现在不同的极值的几个亚群。这是与许多局部最优组合结构优化问题是有价值的。循环停止时进行的分析的数量达到一个预 定义的阈值。由于设计变量数目低，我们10000点分析停止。作为一个线性有限元计算具有约0.15的平均持续时间，一个优化的时间约为25分钟。一个非 线性有限元计算的平均时间为10秒，所以我们选择了只使用5000个评价，以获取不到一天的CPU时间。在本文的其余部分，我们将本次设计方法为M2。第弗朗西斯卡托等人。/复合结构，942018）20872096表

24、4为T1槽优化设计结果得到的层压材料的FPF M1 M2方法和相应的最大应力在每个物料在圆筒部分PW = 13.79兆帕。a 相应的最优的每单位面积的重量在罐筒节。b 价值近似CPT。c 最大应力准则。5.结果与讨论5.1.设计罐T1的结果 根据表2中描述的设计条件下，为T1的优化问题可以表述为：当和分别为螺旋和箍选择之间的最大应力破坏准则。TsaiHill 或者 Hashin标准。在某些条件下7）中，问题P1）仅减少只是M1的方法的第一步骤中，即一个简单的FPF计算。在表4中的所有对应于最小的最佳层压解决方案质量和不同的标准，包括的值中的应力每层，列出和相比，Kim等人的最优解。7。特别是，

25、强调Kim的在不同的帘布层在工作压力下，呃溶液测定由一个简单的线性CLT的方法，因为它们不显示在参考文献。 6,7。这些值确认，Kim的解决方案可被视为FPF的解决方案，工作压力PW=13.79 MPa条件下，因为没有的应力超过在罐的圆筒形部分的各种材料的强度。在这些条件下的问题P1）解决的方法M1在其FPF形式导致的高估设计，因为计算不需要任何变化的由于安全因素的损伤本构层衬套无屈服）和为CFRP安全系数3.0）。我们可能期望FPF设计，最佳的解决方案是独立的选择失败的标准，因为横向的矩阵失败螺旋帘布层中的模式是显性的，和考虑到相同的方式，由不同标准。对于P1），我们的解决方案比Kim重要，

26、这是可以理解的，针对不同的应力得到的每个方法。这种差异可能结果，我们完全理解Kim的设计条件。但这些条件都没有明确界定的文献。6,7，我们提出要解决两个新的问题，P2）和P3），定义分别由关系式8）和9）。 像P1），问题P2）是一个“FPF的问题”，其中只有为Hashin纤维和基体模式的安全因素被区分。对于这个问题，每单位面积的最小权重W大幅度降低，确认的横向的高冲击强度的最佳FPF溶液的CFRP。最佳的的解决方案和相应的权重W获得的M1和M2都非常接近，并轻于Kim的解决方案。导致2螺旋和4相同的圆筒部的层叠体的箍层非常相似的缠绕角度32和28），唯一的不同的是衬套的厚度。 冯M塞斯和纵向

27、应力变化分别在衬套和螺旋层）获得从最佳的P2-M2-FPF溶液表4）表示在图3。为M2的方法，储罐的面积附近的圆顶从优化过程中排除，因为由于应力集中效应边界条件的有限元模型图3a）。在M2方法结果是通过以下方式获得的3-D层状的那些非常相似固体的商业软件ANSYS有限元分析。图3a）。此过程的细节和进一步的结果节中介绍。5.2， M2的方法证明，任何局部弯曲或剪切效应可以证明任何加固的圆周罐的圆筒部的刚性，是清楚本章Kim的解决方案而不是我们的解决方案M1和M2层）。尽管它的简单性，P2-M1-FPF方法似乎是提供一个良好的，但没考虑到本设计主要是罐的近似占主导地位的圆柱部分，即使我们确实观察

28、到一个适度增加纵向应力层螺旋气缸穹顶结附近见表4，图。3aand）第5.2节）。可能是最佳的高压罐设计非常敏感的底部施加的边界条件圆顶和客观需要比较我们的方法和Kim等人的方法。这些规定需要精确的边界条件参考文献的作者。 6。表 最小化的每单位区域W受的重量：最大工作内压PW13:79兆帕不变的几何条件：P1von Mises失效的条件下的线性无屈服：螺旋层的Hashin矩阵失效模式为1:0，安全系数等于 标准：圈层的Hashin标准，安全系数等于3点光纤故障模式 最后，问题P3）中，所描述的关系式9），是一个“FF问题在其中的安全系数的工作直接施加到压力的优化设计，必须取得爆破压力PB =3PW 41.37兆帕。方法M1和M2与Reddy的损害法）被用来解决这个问题，尽管M1-M2和方法之间可能是不同的假设Kim的方法上面提到的。所有类型3高压储罐带有或不带有衬套产量超过一般设计的工作压力，分别设计P3a）和P3b）在关系式9）中。如果产量是允许的，储罐的重量有效地减少，如表5中可以看出，在示出M1和M2的方法得到的最好的解决方案相比Kim的解决方案。对于这个特定的问题，以及相应的M1-最佳的解决方案