CNG-2型气瓶的数值模拟方法对比.pdfVIP

2016年第 9期 玻 璃 钢 ／复 合 材 料 55 CNG-2型气瓶的数值模拟方法对 比 刘培启 ，柴 森 ，陈祖志 ，古纯霖 (1．大连理工大学化工机械学院，大连 116024；2．中国特种设备检测研究院，北京 100013) 摘要：针对数值模拟时CNG一2型气瓶常用的两种模型——轴对称模型和三维壳单元模型存在的争议 ，对两种模型进行 了详细的对比，并讨论 了纤维缠绕角度对其结构的影响，得到以下结论 ：由于CNG-2型气瓶纤维缠绕角接近 9O。，轴对称模型 和三维壳单元模型的应力分布基本相同；壳单元理论上的缺 陷导致三维壳单元模型纤维层的径向应力为0，是两种模型结果 存在差别的主要原因；对纤维缠绕角接近90。的CNG一2型气瓶进行应力分析时，轴对称模型可以代替三维壳单元模型，且计算 工作量较小；当纤维缠绕角不为90。时，两种模型间的误差随缠绕角变化，且 当缠绕角在90。一70。之间时，误差在 5％以内，在工 程允许范围内仍可使用轴对称模型进行应力分析。相关结论可为CNG-2型气瓶数值分析工作提供有益的指导。 关键词 ：CNG一2型气瓶 ；应力分析；壳单元；轴对称；缠绕角度 中图分类号：TB332 文献标识码 ：A 文章编号：1003—0999(2016)09—0055—06 1 引 言 合气瓶进行应力分析时，普遍采用三维实体单元建 目前 ，纤维缠绕复合材料气瓶已广泛应用于航 立内胆模型，壳单元建立纤维层模型 ]。古海波n0] 空航天、军事和民用领域，对这类由内胆和纤维缠绕 使用该方法分析了复合材料气瓶制造、使用、试验及 层两种不同材料组成的结构进行有效的力学分析具 爆破等各工况下的应力和位移分布，并且优化了预 有重要的工程价值。在以往的众多研究中，纤维缠 紧应力和缠绕层厚度。然而壳单元本身也有缺陷， 绕复合材料气瓶应力分析多采用简化的计算方法， 其本构关系基于一阶剪切变形理论，虽然能够计算 建立了一些缠绕张力的解析计算公式 ¨，，但是这 得到层间应力，但计算精度较低，并且 由于未考虑横 些公式仅适用于气瓶的直筒段，不适用变径的封头 向拉伸刚度，层间的法向应力总会默认为零 ¨。可 部分，具有很大的局限性。目前对该类特殊结构应 见各研究人员对纤维层建模方式的选择上仍存在一 力分析主要采用有限元数值模拟的方法_3]。 些分歧。而且在实际的钢制 内衬环缠绕气瓶 (以下 在现有ANSYS、NASTRAN等大型商业软件中， 简称CNG一2型气瓶)的缠绕工艺中，为便于在封头 对纤维缠绕复合材料气瓶进行应力分析时，纤维缠 处过渡缠绕，防止缠绕层脱落，并使缠绕层周向、轴 绕层的建模通常采用两种单元：三维实体单元和壳 向都能承受一定的应力，气瓶厂家普遍采用 90。环 单元_4J。苏镜元 同时用三维实体单元建立 内胆 缠绕加非90。螺旋缠绕相结合的方式，这种缠绕方 和缠绕层的有限元模型，并且在金属内胆和纤维层 式会破坏CNG一2型气瓶本身的轴对称特性[1’]，进 之间、成型筒体和缠绕带之间创建了两个接触对，但 一 步加剧上述分歧。 是纤维缠绕层层数一般有十几层甚至上百层，这种 本文建立了CNG一2型气瓶的轴对称模型，对比 方法划分的三维实体单元数量密集且普遍呈薄片 三维实体单元和壳单元分析结果的差异，讨论是否 状，可能会造成总体刚度矩阵出现病态，增加求解机 可以使用轴对称模型代替三维壳单元模型，并且分 时，并可能导致求解结果出现误差 ]，同时这种建模 析了纤维缠绕角度对两种模型分析结果差异的影 方式还存在当缠绕角不为90。时难以准确调整实体 响，讨论轴对称模型代替三维壳单