波纹管ANSYS有限元分析

1、第一章综述1.1 波纹管的发展波纹管膨胀节是压力容器和管道系统中常用的连接和补偿装置，其主体是金属波纹管。它的特点是轴向刚度较小，径向刚度较大，故能够承受较大的轴向位移，同时还能够承受一定的内、外压，常用于补偿管道因温度变化、机械位移或振动引起的轴向、横向和角向位移。运行中的波纹管的波峰波谷部位实际处于弹塑性区域，但是，目前波纹管常用的设计标准是将波纹管视为梁或曲杆，通过梁理论进行弹性的应力应变分析，这势必会带来较大的误差。随着数值分析技术和有限元技术的不断发展，利用非线性有限元技术对波纹管进行较为精确的应力应变分析已成为波纹管研究领域的一项重要课题。本文吸取了前人在有限元程序设计上的经验，根

2、据波纹管自身的非线性特性，利用FortranPowerSta-tion410与VisualC+610混合语言编程技术,开发出一套界面友好，使用方便，求解精度较高的波纹管有限元应力应变分析专用软件系统。1.2 波纹管非线性特性与有限元程序设计固体力学问题，从本质上讲都是非线形的，线性假设只是实际问题的一种简化。线性分析的基本假设包含三点：结点位移为无限小量，材料为线弹性，加载时边界条件的性质保持不变。如果上述假设中的任一条不能满足，则属于非线性问题。常见的是不满足前两种假设的非线性行为，分别称为几何非线性和材料非线性。1.3 波纹管的非线性特性工作状态下的波纹管伸缩变形是一种复杂的非线性变形。在

3、小载荷情况下，波纹管的应力应变关系呈线性特征，而当其承受较大的位移补偿时，结构处于高应力水平，波峰和波谷的局部处于塑性区，此时，一般采用弹性线性强化模型分析波纹管的应力应变关系。另外，波纹管的单波位移量一般大于壁厚,属于大位移行为，同时，大量的研究表明，此时的波纹管大部分区域依然处于小应变范围。因此，运行状况下的波纹管还表现出大位移、小应变的几何非线性特征。1.3.1 非线性有限元程序设计本软件系统所用的求解程序是一个通用的二维非线性有限元程序，图1是该程序的流程图，其主体部分借鉴了文献132中提出的有限元程序设计方法。ZERO把需要累加的数组置零图1.1二维非线性有限元程序流程RF顼DU计算

4、残余力矢最FRONT潴前法求解联立方程组EEFFP计算单元刚窿矩阵ALGOR设置变量以识别解皆觉型COLVER检杳求解过程是否收放|迭代结拜期布增量循环.1.3.2 波纹管有限元计算模型根据波纹管的几何特征和载荷特征，系统采用二维轴对称模型进行求解，单元类型为8结点抛物线等参单元，并开发出相应的前处理程序。以U型波纹管为例，通过输入波纹管的基本波形参数、材料参数和运行状况参数，便可建立如图2所示的计算模型，自动进行网格划分，边界约束和等效结点力加载。图1.2波纹管二维轴对称有限元计算模型波纹膨胀节被广泛应用于管路热补偿、隔振、降噪。自40年代以来，世界各工业国家的研究机构和生产厂家对作为其主要

5、组成部分的波纹管进行了广泛的研究。在不同工况下对波纹管进行应力）应变的分析是评估其对各种工况承受能力的基础和重要组成部分。70年代之前的结构分析主要采用解析法，根据近似简单梁、近似圆柱体、近似壳体的假设，依据弹性理论得到波纹管壳体中的应力）应变值1。但由于波纹管本身是一种较为复杂的轴对称薄壁壳体，且在绝大多数工况下材料处于弹塑性大变形范围内，因而在解析解与波纹管材料的实际响应之间存在着较大的误差。为使两者间达到近似的统一，学者们在解析法得到的近似解中引入了由图表形式给出的修正系数。这种方法目前正广泛应用于设计以及生产部门，如EJMA规范中的应力计算式2o70年代以后，随着计算机技术的飞速发展，

6、有限元分析在结构分析中的地位日渐突出，许多商用有限元软件应运而生。有限元分析以其在解决几何非线性、材料非线性和结构非线性问题方向的独特优势而为愈来愈多的研究人员所接受。非线性有限元计算能够有效地解决传统的解析法无法计算的波纹管弹塑性大变形范围内的载荷）应力响应问题。但由于以往有限元计算软件方面的局限性，其分析结果距离能够真正运用到波纹管的强度、刚度以及稳定性设计方面还有一定的差距。本文建立了波纹管非线性分析模型,并采用大型有限商用软件ANSYS!行了计算，计算结果与试验数据进行了比较，得到了与波纹管设计计算有关的载荷）应力响应关系，并就波纹管非线性有限元计算进行研究。第二章建模过程及求解2.1

7、 波纹管结构示意图图2.1参数含义2.2 第一组数据表2-1关外径内径壁厚波距波高波纹管波峰半波谷半层波参工DD2STB长度L径R1径R2数数(mm)(mm)(mm)(mm)(mm)(mm)(mm)(mm)(m)(m)数二7736981907544023.721.314值2.2.1 具体坐标1(0.0.0)2(61.3.0.0)3(61.3.70.0)4(108.7.70.0)5(108.7.0.0)6(151.3.0.0)7(151.3.70.0)8(198.7.70.0)9(198.7.0.0)10(241.3.0.0)11(241.3.0.0)12(288.7.70.0)13(288.7

8、.0.0)14(331.3.0.0)15(331.3.70.0)16(378.7.70.0)17(378.7.0.0)18(440.0.0)19(0.-349.0)20(440.-349.0)2.2.2 建立坐标关键点ANSYSMainMenu:选择Preprocessor-Modeling-Create-Keypoints一InActiveCS填写对应坐标图2.2坐标点图2.2.3 生成特征点并且连接、倒圆角、布尔操作选择Preprocessor|Modeling|Create|Lines|StraightLineW出CreateStraightLine框，连接各点。如下图所示选择Prepr

9、ocessor|Modeling|Operate|Booleans|Add|Lines图2.3坐标点连线2.2.4 生成壳体选择Preprocessor|Modeling|Operate|Booleans|Add|Lines图2.4螺线管壳体2.2.5 划分网格ANSYSMainMenu:Preprocessor-Meshing-MeshTooH(SizeControls)lines图2.5网格划分2.2.6 分析计算选择Solution|Solve|CurrentLS得出结果2.2.7 偏移量,-w;=:±:'5;一框三直LIJILL.i-.HnJ_.:.1-M.!:.:!

10、:r二二JI副圈制园Ei-Fm":='fl一滴£±一.-F!T-1rt一一1卜MHm二-k吗昌¥泪三受1|督13一1*IWEHEtBUDI：3rUfEM3nI9W7工UtiETS口"才图2.6偏移图2.2.8 分析计算图2.7应力云图这个是应力云图，从图2.7中可以看出当最大位移是45mmi勺时候最大应应力是491mpa最小应力是3.3281mpa.最大应力接近屈服极限480mpa,因此可以认为波纹管最大补偿是45mm图2.8位移云图这个是位移图2.8,当位移是45mm的时候达到304材料的最大屈服极限2.3 第二组数据表2-2参数外径

11、D1(mm)内径D2(mm)壁厚S(mm)波距T(mm)波高B(mm)波纹管长度L(mm)波峰半径R1(mm)波谷半径R2(mm)层数(m)波数(m)数值7236481907544023.721.3142.3.1 具体坐标1(0.0.0)2(61.3.0.0)3(61.3.70.0)4(108.7.70.0)5(108.7.0.0)6(151.3.0.0)7(151.3.70.0)8(198.7.70.0)9(198.7.0.0)10(241.3.0.0)11(241.3.0.0)12(288.7.70.0)13(288.7.0.0)14(331.3.0.0)15(331.3.70.0)16(

12、378.7.70.0)17(378.7.0.0)18(440.0.0)19(0.-324.0)20(440.-324.0)2.3.2 定义材料参数因为我们选的材料是304，上网查的304钢所对应的弹性模量是1.9e11pa,泊松比是0.3052.3.3 建立坐标关键点ANSYSMainMenu:选择Preprocessor-Modeling-Create-Keypoints一InActiveCS填写对应坐标图2.9坐标点图2.3.4 生成特征点并且连接、倒圆角、布尔操作选择Preprocessor|Modeling|Create|Lines|StraightLine单出CreateStraig

13、htLine框，连接各点。如下图所示选择Preprocessor|Modeling|Operate|Booleans|Add|Lines图2.10坐标连线2.3.5 生成壳体选择Preprocessor|Modeling|Operate|Booleans|Add|Lines图2.11螺纹管壳体2.3.6 划分网格ANSYSMainMenu:Preprocessor-Meshing-MeshTooH(SizeControls)linesAN：P.SI图2.12网格划分2.2.7分析计算选择Solution|Solve|CurrentLS得出结果2.3.8 偏移量图2.13偏移图JTH16-,3Q

14、LC>2.3.9 应力及波纹处受力分析图2.14应力云图这个是应力云图，从图2.14中可以看出当最大位移是43mm的时候最大应应力是436.04mpa,最小应力是9.1958mpa.最大应力接近屈服极限480mpa,因此可以认为波纹管最大补偿是43mm图2.15位移云图这个是位移云图2.15,当位移是43mm的时候达到304材料的最大屈服极限此时的位移补偿量为43mm第三章结论说明3.1 通过上面两组数据的对比分析，得由以下结论(1)内径对波纹管的补偿量是很大的，第一组数据内外径D=773mmd=698波纹管的补偿量是45mm勺时候它的应力接近最大屈服极限480mpa;第二组数据内外径D

15、=723mmd=648波纹管的补偿量是43mm勺时候它的应力接近最大屈服极限480mpa.内外径相差50mm波纹管的补偿量就相差2mm由此可以看出内外径对波纹管的影响是不可忽略的(2)波纹管通常要求有尽可能大的补偿量，因此波纹管在满足承压的前提下，尽可能减小内径，以获得较高的柔性。参考文献1叶青.微机原理与接口技术.清华大学出版社，2011:9-1.2薛华成.管理信息系统M.北京：清华大学出版社，1993:99-103.3李晖.高性+能计算机若干关键问题研究D.中国科学技术大学，2009:125-138.沈阳化工大学科亚学院工程软件实践应用专业：机械设计制造及其自动化班级:学生姓名：指导教师:

16、工程软件应用实践任务书机械制造及其自动化专业班学生:论文题目：内外径对U形波纹管膨胀节补偿能力的影响具体参数:参数1参数2外径Di(mm)773723内径D2(mm)698648壁厚(mm)11波距T(mm)9090波高B(mm)7070波纹管长度L(mm)440440波峰半径R1(mm)23.723.7波谷半径R2(mm)21.321.3层数m11波数n44材料304304指导教师：201年月日目录第一章综述11.1 波纹管的发展11.2 波纹管非线性特性与有限元程序设计11.3 波纹管的非线性特性11.3.1 非线性有限元程序设计21.3.2 波纹管有限元计算模型2第二章建模过程及求解42.1 波纹管结构示意图42.2 第一组数据42.2.1 具体坐标42.2.2 建立坐标关键点42