《毕业论文（设计）基于ISO 11439 标准及ANSYS 平台的CNG 气瓶有限元应力分析08956.doc》由会员分享，可在线阅读，更多相关《毕业论文（设计）基于ISO 11439 标准及ANSYS 平台的CNG 气瓶有限元应力分析08956.doc（12页珍藏版）》请在三一文库上搜索。
1、基于ISO 11439标准及ANSYS 平台的CNG气瓶有限元应力分析岳中第 马 宁 刘 彬（北京航空制造工程研究所，100024）摘要：本文针对铝胆复合材料气瓶设计与制造过程中提出的有限元应力分析问题，根据ISO 11439标准，对若干关键技术进行了综述和研究，突出的应用成果说明了这些技术的可靠性与实用性。1 CNG气瓶设计及其应用背景气瓶作为一种没有固定使用地点和环境、无专职使用和操作人员的特殊压力容器，已经在国内外得到了广泛的研究和应用。突出的范例有医护领域的自救呼吸装置及数以百万计的各类公交汽车用的压缩天然气气瓶。其中，复合材料缠绕气瓶是用钢胆或铝胆作为内衬，外面用高强度纤维复合材料缠
2、绕制成，具有良好的抗疲劳、抗晶间腐蚀性能，且有较轻、便于搬运及应用等优点。目前，发达国家已有300多万只复合材料缠绕气瓶在安全的使用着。传统工程意义的“薄膜理论”，“网格理论”，对复合材料气瓶的设计与分析计算，已不能满足要求。气瓶作为一种特殊的压力容器，在安全性、稳定性及可靠性等性能评定方面，根据ISO 11439国际标准的要求，有严格的设计、制造、实验及性能验证的规范。本文涉及的CNG-3型气瓶，是ISO 11439标准中其制造技术最为复杂的一种，内胆材料为铝，纤维全缠绕。目前，铝内胆的设计与制造工艺，国内外已经形成十分先进与稳定的专利技术。作者针对全缠绕的CNG-3型气瓶，进行了材料非线性
3、有限元分析方面的探索，并充分利用ANSYS的二次开发平台，开发了适合CNG-3等类型气瓶在各种压力工况下的应力分析与性能评定的软件，从而能够准确、快速、方便地确定复合材料气瓶是否符合工程应用环节的要求。2 对CNG-3型气瓶进行应力分析的若干关键技术 根据设计制造CNG-3型气瓶的基本实践，我们必需面对并真正解决若干关键技术问题。2.1 ISO 11439对复合气瓶的应力分析要求根据ISO 11439国际标准，此类气瓶总体上被分为四类，即CNG-1，CNG-2，CNG-3和CNG-4。其中， CNG-1指全金属型气瓶；CNG-2指金属内衬纤维缠绕浸渍树脂型气瓶，其纤维缠绕方向只是周向单向缠绕；
4、CNG-3也指金属内衬纤维缠绕浸渍树脂型气瓶，但其纤维缠绕方向是全方向的；CNG-4指非金属内衬纤维缠绕型气瓶，即全复合材料气瓶（如内衬为塑料材料）。根据我们的理解，对于CNG-3，ISO设计规范的基本要求为： 建立可用的材料非线性分析方法与流程； 建立正确的材料非线性的弹-塑性、应力-应变模型； 建立全缠绕状态的复合材料机械性能参数计算方法； 建立多纤维缠绕层或不同纤维材料构成的材料本构方程与平衡方程； 分析纤维缠绕张力作用下的预应力； 不同压力工况的应力分析；必须进行紧缩压力、零压力、工作压力、水压试验压力、最小爆炸压力等工况下的材料非线性应力计算； 最大的缺陷尺寸，应能确保在循环压力工况
5、下的安全可靠性，同时可保证LBB性能（爆炸前先泄露）。2.2 铝胆与复合材料缠绕层的组合模型对于铝胆复合材料气瓶，其外部纤维缠绕层是由高强度增强纤维及其基体，如玻璃纤维/环氧树脂、有机纤维/环氧树脂及碳纤维/环氧树脂等缠绕而成的。因此，它们具有比强度大、比刚度高等诸多的优点。在高压下，铝胆很快进入塑性，而纤维增强复合材料还在弹性状态。实验表明，这种组合模型至少在一个材料方向上呈现出应力-应变关系的非线性。同时，石墨纤维/环氧树脂、Kevlar纤维/环氧树脂及硼纤维/环氧树脂具有相当高的剪切非线性；硼纤维/铝复合材料具有相当高的横向或剪切非线性；碳/碳复合材料在各材料主轴方向均存在着较高的非线性
6、等。这种非线性主要是由于基体性能的非线性和在受载过程中所产生的损伤缺陷引起的，它与基体的种类、纤维的形态、复合体积比率、纤维缠绕方向及承受载荷形式有关。显然，这是一个复合材料本构关系问题，通过应力-应变的物理非线性表现出来。如何准确地描述本构关系，对于我们能否较精确地进行力学分析和结构设计是非常重要的。材料非线性是指材料本构关系的非线性。若不考虑变形对材料的影响，平衡方程简化成线性，同时应变和位移的关系也是线性的，这就是小弹塑性问题。在塑性区域内，应力和应变之间不存在弹性范围的广义虎克定律。塑性变形的过程是不可逆的，变形功部分转化为热能，最后状态的应力与变形途径有关。2.3预应力结构与材料非线
7、性 正确应用预应力技术，是纤维缠绕复合材料气瓶的一个显著特点，也是气瓶应力分析中的一个重要环节。预应力结构形式是在60年代由预应力钢丝缠绕结构发展起来的。随着航空航天工业的发展，以及许多部门如原子能工业中的高速离心机和导弹武器系统等对机械结构高强度、低重量和高可靠性方面提出了更高的要求，碳纤维等非金属材料缠绕式预应力结构得到快速的发展和应用。 预应力容器的基本力学特点是，内壁是其强度最薄弱之处，起关键作用的是切向拉应力，而径向压缩应力恒等于内压之负值，这是无法改变的。预应力容器的基本设计思想，就是在内壁预先建立一种切向压应力预应力，和拉伸的工作切向应力叠加，以大大减少甚至于完全消除切向拉应力，
8、在合成状态建立起切向压应力或压应变。这样，就带来一系列静强度和疲劳强度上的优点，它们是： 材料的不同力学性能得到充分发挥； 承载能力大； 抗疲劳强度高； 安全可靠，无爆炸危险。2.4 金属材料双线性各向同性强化模式 s A A B sy 拉伸 O e 2sy 压缩 B 图1 铝胆的Bauschinger效应 铝胆复合材料气瓶应力分析的的关键性技术之一，是实现铝胆复合材料气瓶的材料非线性有限元分析。铝胆复合材料气瓶在设计中充分利用了结构的材料非线性特点，通过自紧压力（AUTOFRETTAGE）的拉伸, 使金属铝胆达到屈服状态，产生塑性效应（即Bauschinger 效应），但外层缠绕的复合材料由
9、于高强度而处在线性变形状态。当结构进行卸载至零压力，铝胆因塑性变形不能恢复原始状态，受到外层材料的紧缩压力，产生了自紧应力。 图1表明了铝胆的Bauschinger效应，拉伸方向的屈服应力增加，导致后继的压缩屈服应力在数量上降低，因此屈服应力高低之间存在2sy的差异。这种材料双线性各向同性强化模式代表了应力应变曲线以及Von Mises的屈服准则。2.5 ANSYS与数值模拟技术许多工程分析问题，都归结为在给定边界条件下求解它们的控制方程（常微分方程或偏微分方程）的问题。但是，能用解析方法求出精确解的只有方程性质比较简单，且几何边界相当规则的少数问题。对于大多数的工程技术问题，由于物体的几何形
10、状复杂或者问题的某些特征是非线性的，则很少有解析解。因此，人们在广泛吸收现代数学、力学理论的基础上，借助于现代高速计算机来获得满足工程要求的数值解，这种数值模拟技术是现代工程学形成和发展的最重要推动力之一。有限元法作为工程应用中的一种最重要的数值模拟技术的地位越来越得到人们的公认。ANSYS系统是国际上最具知名度的有限元分析软件之一，也已作为我国压力容器标准委员会推荐的基础软件之一，用于压力容器性能分析与评定。它不仅是融结构、热、流体、电磁、声学于一体的大型通用有限元分析软件，而且还为用户提供了一个不断改进功能的二次开发平台，具体包括：结构高度非线性分析、电磁分析、计算流体动力学分析、优化设计
11、、接触分析、自适应网格划分、大变形/有限转动功能以及利用ANSYS参数设计语言（APDL）的扩展宏命令功能。FEACAC（The Finite Element Analysis for the Cylinder of Aluminum-liner with Composites），便是在 ANSYS二次开发平台发展的，它使CNG的设计分析与性能评定完全实现了三“化”（通用化、系列化、参数化），成为CNG快速设计与性能评定的重要基础手段。3 基于ISO 11439及ANSYS 的FEACAC 开发国际上先进的CNG设计制造商提供的EPDRP技术报告（Engineering Preliminary
12、 Design Review Package）表明，对复合材料气瓶进行有限元应力分析计算是气瓶设计最重要的基础技术。其中，七种工况（自紧压力、零压力、工作压力、实验压力、最小爆炸压力等）的结构非线性应力分析，是设计的核心关键技术，也是CNG-3国际标准认证所要求的基本技术环节。我们在 ANSYS二次开发平台上发展的FEACAC， 包括了三个宏程序包。它们是： （1） ALUM-专用于金属内胆的参数化FEA宏程序；（2） COMP_M-专用于金属胆复合气瓶方案设计的参数化FEA宏程序；（3） COMP-专用于金属胆CNG-3型复合气瓶整体强度校核的参数化FEA宏程序。FEACAC这些宏程序包分别
13、为铝胆设计及旋压加工、复合气瓶方案设计与优化、性能评定与应力校核开发的。图2给出了FEACAC的运行处理流程。1 FEACAC程序文件 用户参数库文件 实用函数库文件2 调用ANSYS 线性与非线性分析器，建立结构非线性分析流程3 调用ANSYS 后置处理器，给出应力分析结果与强度校核报告图2 FEACAC处理流程其中， 框1是用APDL语言开发的FEACAC宏程序库，自动产生FEA几何模型，并建立各压力工况的载荷模型与边界约束条件。框2建立了有效的结构非线性流程，及选择了合适的求解大型稀疏平衡方程的算法。它利用ANSYS功能强大的结构非线性分析器，进行线性与非线性的FEA分析。框3调用ANS
14、YS 后置处理器，直接产生应力分析结果与强度校核报告，给出金属胆及各复合材料纤维缠绕层纤维方向的应力结果，并自动给出与标准应力分析结果的校核误差报告。3.1 ALUM图3是FEACAC的金属胆（如，ALUM.）FEA的软件树结构, 包括4个程序。其中，*.le表示线性分析，模型由实验曲线给出；*.lt表示线性分析，模型由椭圆曲线理论给出。2d/3d表明了2维与3维两类不同模型。ALUMAlum_3d.ltAlum_3d.leAlum_2d.ltAlum_2d.leA_ser.lib，Cylmac.lib图3. ALUM(铝胆) FEA软件的树结构3.1 COMP_M图4是FEACAC的方案设计
15、软件树结构。此时仅取中段圆柱体，按2d模型程序计算。其中，comp_2d.ns表示单一工况的非线性分析；comp_2d.na表示七种工况的非线性分析。Verify.ns，verify.na是对上述两分析程序的结果应力进行处理的宏程序，它将给出方案设计的非线性应力分析报告，以及与某些标准应力对比的校核报告。实践表明，COMP_M是铝胆复合气瓶设计使用频度最大、最快速的方案选择与优化的重要工具。COMP_MVerify.naVerify.nsComp_2d.naComp_2d.nsM_User.lib，Cylmac.lib图4. FEACAC初步设计阶段的软件树结构3.2 COMP COMPCom
16、p_2d.enaComp_2d.ensComp_2d.tnaComp_2d.tnsC_user.lib，Cylmac.lib图5. FEACAC详细设计阶段的软件树结构 图5是FEACAC详细设计阶段的软件树结构。这是在ANSYS平台上开发的专用软件，它将建立整体结构强度校核的模型，包括4个2d的模型产生与非线性分析程序。其中，comp_2d.tns表示单一工况、具有理论模线的非线性分析；comp_2d.tna表示七种工况、具有理论模线的非线性分析；comp_2d.ens表示单一工况、具有实验外型模线的非线性分析；comp_2d.ena表示七种工况、具有实验外型模线的非线性分析。Verify.
17、ns，verify.na是对上述两分析程序的结果应力进行校核的程序。4 工程应用及结论研制某产品的实践表明，基于ANSYS平台及国际标准ISO 11439开发的FEACAC软件与美国知名专业技术资询公司(HEI)所用的FEAS软件相比较，针对同一产品的两组不同的设计与制造参数，进行七种载荷工况（自紧压力、零压力、工作压力、循环压力、试验压力、最小爆破压力、平均爆破压力）的工程分析，能够得到相同工程精度的数值分析结果，附录1表明其中一种产品设计的应力校核误差分析。值得注意的是，由于结构工作压力与循环压力工况下，其内胆环向应力接近于零或压缩状态，因而相应的绝对误差较小，但相对误差较大。图6、7分别
18、是2d加载模型与零压力下自紧应力云图。图8、9分别是封头及过渡段结构的局部应力云图。我们认为，对于高压容器的设计与制造来说，基于ISO 11439标准与ANSYS平台的二次开发是十分重要的，它能使用户以较低的成本和较短的时间与国际设计与认证标准接轨，对市场的新产品设计与开发形成快速有效的反应能力。5参考文献1. INTERNATIONAL STANDARD ISO 11493，Gas cylinders-High pressure cylinders for the on-board storage of natural gas as a fuel for automotive vehicle
19、s. 20002. 王勖成 邵敏。有限元法基本原理和数值方法。清华大学出版社，19973. 张志民。复合材料结构力学。北京航空航天大学出版社，19934. 黄再满 蒋鞠慧 薛忠民 黄毓圣。复合材料天然气气瓶预紧压力的研究5. 段登平，刘正兴，罗海安。纤维缠绕壳体材料非线性及大变形分析计算。复合材料学报，第16卷，第1期：142-1486. 刘永纯 董雨达。复合材料压缩天然气（CNG）气瓶应用研究。 纤维复 合材料。 2000，17（3），50-527. 王国强。 实用工程数值模拟技术及其在ANSYS上的实践。 西北工业大学出版社，1999 图6 FEA的2D整体加载模型 图7 零压力自紧应力云
20、图图8 封头的局部应力分布 图9 结构过渡段的应力分布云图附录1. 应力校核误差分析表（一） Stress Verification I On Sidewall Cylinder From BAMTRI and HyPerComp 1. STRESS NAME LIST OF RESULTS ITEMS (1). PRESSURE = LINER INTERNAL PRESSURE待添加的隐藏文字内容3 (2). LOGITU = LINER LOGITUDINAL STRESS (3). HOOP = LINER HOOPS STRESS (4). V.MISES = VON MISES S
21、TRESS (5). C.HOOP = CARBON HOOPS FIBER STRESS (6). C.HELIC = CARBON FULL HELICALS FIBER STRESS (7). C.HIGH = CARBON HIGH ANGLE HELICALS FIBER (8). G.HELIC = GLASS HELICALS FIBER STRESS (9). G.HOOP = GLASS HOOPS FIBER STRESS 2. RELATIVE TOLERANCE LIST OF STRESS RESULTS COMPARISONS PRESSURE LOGITU HOO
22、P V.MISES C.HOOP C.HELIC C.HIGH G.HELIC G.HOOP ( Note: PRESSURE in MPa, STRESS in MPa) (1). AUTOFRETTAGE PRESSURE 0.0% -4.0% -2.7% 0.1% 2.3% -1.5% -0.7% -2.0% 2.7% (2). ZERO PRESSURE 0.0% 0.2% 2.2% -2.5% -1.4% -1.9% -5.1% -3.0% -1.4% (3).WORKING PRESSURE 0.0% -19.5% -19.0% -9.7% 1.4% -1.7% -1.6% -2.
23、3% 1.7% (4). CYLCLE PRESSURE 0.0% -10.0% -6.6% -3.0% 1.7% -1.6% -1.4% -2.0% 2.0% (5). TEST PRESSURE 0.0% -6.2% -4.2% -1.4% 2.0% -1.6% -1.1% -2.2% 2.3% (6). MIN. BURST PRESSURE 0.4% -4.4% -2.6% 0.8% 3.1% -0.7% 0.2% -0.9% 3.7% (7). AVERGE BURST PRESSURE 0.4% -6.7% -0.7% 2.4% 3.2% -0.1% 0.9% -0.3% 4.5%
24、 * Stress Data Was Finished 本文作者简介岳中第： 研究员, 政府特殊津贴专家，长期从事飞机结构强度分析软件开发及CAD/CAE集成研究，参与了对德、美的国际技术合作项目。曾负责航空领域和北京市多个工程项目，获得近10项国家及部级科学技术进步奖。马 宁：工程师，从事气瓶设计与制造。刘 彬：硕士研究生，从事气瓶设计与性能评定方面的研究。Editors note: Judson Jones is a meteorologist, journalist and photographer. He has freelanced with CNN for four years,
25、covering severe weather from tornadoes to typhoons. Follow him on Twitter: jnjonesjr (CNN) - I will always wonder what it was like to huddle around a shortwave radio and through the crackling static from space hear the faint beeps of the worlds first satellite - Sputnik. I also missed watching Neil
26、Armstrong step foot on the moon and the first space shuttle take off for the stars. Those events were way before my time.As a kid, I was fascinated with what goes on in the sky, and when NASA pulled the plug on the shuttle program I was heartbroken. Yet the privatized space race has renewed my child
27、hood dreams to reach for the stars.As a meteorologist, Ive still seen many important weather and space events, but right now, if you were sitting next to me, youd hear my foot tapping rapidly under my desk. Im anxious for the next one: a space capsule hanging from a crane in the New Mexico desert.It
28、s like the set for a George Lucas movie floating to the edge of space.You and I will have the chance to watch a man take a leap into an unimaginable free fall from the edge of space - live.The (lack of) air up there Watch man jump from 96,000 feet Tuesday, I sat at work glued to the live stream of t
29、he Red Bull Stratos Mission. I watched the balloons positioned at different altitudes in the sky to test the winds, knowing that if they would just line up in a vertical straight line we would be go for launch.I feel this mission was created for me because I am also a journalist and a photographer,
30、but above all I live for taking a leap of faith - the feeling of pushing the envelope into uncharted territory.The guy who is going to do this, Felix Baumgartner, must have that same feeling, at a level I will never reach. However, it did not stop me from feeling his pain when a gust of swirling win
31、d kicked up and twisted the partially filled balloon that would take him to the upper end of our atmosphere. As soon as the 40-acre balloon, with skin no thicker than a dry cleaning bag, scraped the ground I knew it was over.How claustrophobia almost grounded supersonic skydiverWith each twist, you
32、could see the wrinkles of disappointment on the face of the current record holder and capcom (capsule communications), Col. Joe Kittinger. He hung his head low in mission control as he told Baumgartner the disappointing news: Mission aborted.The supersonic descent could happen as early as Sunday.The
33、 weather plays an important role in this mission. Starting at the ground, conditions have to be very calm - winds less than 2 mph, with no precipitation or humidity and limited cloud cover. The balloon, with capsule attached, will move through the lower level of the atmosphere (the troposphere) wher
34、e our day-to-day weather lives. It will climb higher than the tip of Mount Everest (5.5 miles/8.85 kilometers), drifting even higher than the cruising altitude of commercial airliners (5.6 miles/9.17 kilometers) and into the stratosphere. As he crosses the boundary layer (called the tropopause), he
35、can expect a lot of turbulence.The balloon will slowly drift to the edge of space at 120,000 feet (22.7 miles/36.53 kilometers). Here, Fearless Felix will unclip. He will roll back the door.Then, I would assume, he will slowly step out onto something resembling an Olympic diving platform.Below, the
36、Earth becomes the concrete bottom of a swimming pool that he wants to land on, but not too hard. Still, hell be traveling fast, so despite the distance, it will not be like diving into the deep end of a pool. It will be like he is diving into the shallow end.Skydiver preps for the big jumpWhen he ju
37、mps, he is expected to reach the speed of sound - 690 mph (1,110 kph) - in less than 40 seconds. Like hitting the top of the water, he will begin to slow as he approaches the more dense air closer to Earth. But this will not be enough to stop him completely.If he goes too fast or spins out of contro
38、l, he has a stabilization parachute that can be deployed to slow him down. His team hopes its not needed. Instead, he plans to deploy his 270-square-foot (25-square-meter) main chute at an altitude of around 5,000 feet (1,524 meters).In order to deploy this chute successfully, he will have to slow t
39、o 172 mph (277 kph). He will have a reserve parachute that will open automatically if he loses consciousness at mach speeds.Even if everything goes as planned, it wont. Baumgartner still will free fall at a speed that would cause you and me to pass out, and no parachute is guaranteed to work higher
40、than 25,000 feet (7,620 meters).It might not be the moon, but Kittinger free fell from 102,800 feet in 1960 - at the dawn of an infamous space race that captured the hearts of many. Baumgartner will attempt to break that record, a feat that boggles the mind. This is one of those monumental moments I will always remember, because there is no way Id miss this.