波纹管ansys有限元分析.docxVIP

第一章综述 波纹管的发展 波纹管膨胀节是压力容器和管道系统屮常用的连接和补偿装置，其主体是金 属波纹管。它的特点是轴向刚度较小，径向刚度较人，故能够承受较人的轴向位移, 同时还能够承受一定的内、外压，常用于补偿管道因温度变化、机械位移或振动引 起的轴向、横向和角向位移。 运行中的波纹管的波峰波谷部位实际处于弹塑性区域，但是，目前波纹管常 用的设计标准是将波纹管视为梁或曲杆，通过梁理论进行弹性的应力应变分析， 这势必会带来较人的误差。随着数值分析技术和有限元技术的不断发展，利用非线 性有限元技术对波纹管进行较为精确的应力应变分析己成为波纹管研究领域的一 项重要课题。 本文吸取了前人在有限元程序设计上的经验，根据波纹管自身的非线性特性， 利用Fortran Pow erSta-tion 410与 Visual C++ 610混合语言编程技术，开发出一 套界面友好，使用方便，求解精度较高的波纹管有限元应力应变分析专用软件系 统。 1.2波纹管非线性特性与有限元程序设计 固体力学问题，从本质上讲都是非线形的，线性假设只是实际问题的一种简 化。线性分析的基木假设包含三点：结点位移为无限小量，材料为线弹性，加载时 边界条件的性质保持不变。如果上述假设中的任一条不能满足，则属于非线性问题。 常见的是不满足前两种假设的非线性行为，分别称为几何非线性和材料非线性。 1.3波纹管的非线性特性 工作状态下的波纹管伸缩变形是一种复杂的非线性变形。在小载荷情况下，波 纹管的应力应变关系呈线性特征，而当其承受较大的位移补偿时，结构处于高应力 水平，波峰和波谷的局部处于塑性区，此时，一般采用弹性线性强化模型分析波纹 管的应力应变关系。另外，波纹管的单波位移量一般大于壁厚，属于大位移行为， 同时，大量的硏究表明，此时的波纹管大部分区域依然处于小应变范围。因此，运 行状况下的波纹管还表现出大位移、小应变的儿何非线性特征。 戟荷增%循环ALGOR设置变量以识别解法类熨INCREM按给定的載荷系数增盘加載STIFFP计算单元刚度矩年迭代循环FRONT波前法求解联立方程组RESIDU计算残余力矢量CONVER 戟荷增%循环 ALGOR设置变量以识别解法类熨 INCREM按给定的載荷系数增盘加載 STIFFP计算单元刚度矩年 迭代循环 FRONT波前法求解联立方程组 RESIDU计算残余力矢量 CONVER检査求解过程是否收敛 OUTPUT输出本次栽荷增精下的结果 图1.1二维非线性有限元程序流程 1. 3. 2波纹管有限元计算模型 根据波纹管的几何特征和载荷特征，系统采用二维轴对称模型进行求解，单元 类型为8结点抛物线等参单元，并开发出相应的前处理程序。以U型波纹管为例, 通过输入波纹管的基本波形参数、材料参数和运行状况参数，便可建立如图2所示 的计算模型，自动进行网格划分，边界约朿和等效结点力加载。 图1.2波纹管二维轴对称有限元计算模型 1. 3.1非线性有限元程序设计 木软件系统所用的求解程序是一个通用的二维非线性有限元程序，图1是该 程序的流程图,其主体部分借鉴了文献132中提出的有限元程序设计方法。 程序的流程图, 波纹膨胀节被广泛应用于管路热补偿、隔振、降噪。自40年代以来，批界各 工业国家的研究机构和生产厂家对作为其主要组成部分的波纹管进行了广泛的研 究。 在不同工况下对波纹管进行应力）应变的分析是评估其对各种工况承受能力的 基础和重要组成部分。70年代之前的结构分析主要采用解析法，根据近似简单梁、 近似圆柱体、近似壳体的假设，依据弹性理论得到波纹管壳体中的应力）应变值 ［1］。但由于波纹管木身是一种较为复杂的轴对称薄壁壳体，且在绝大多数工况下材 料处于弹塑性大变形范围内，因而在解析解与波纹管材料的实际响应之间存在着较 大的误差。为使两者间达到近似的统一，学者们在解析法得到的近似解中引入了由图 表形式给出的修正系数。这种方法FI前正广泛应用于设计以及生产部门，如EJMA 规范中的应力计算式［刀。 70年代以后，随着计算机技术的飞速发展，有限元分析在结构分析中的地位日 渐突出，许多商用有限元软件应运而生。有限元分析以其在解决儿何非线性、材料 非线性和结构井线性问题方向的独特优势而为愈来愈多的研究人员所接受。非线性 有限元计算能够有效地解决传统的解析法无法计算的波纹管弹塑性大变形范围内的 载荷）应力响应问题。但由于以往有限元计算软件方面的局限性，其分析结果距离 能够真正运用到波纹管的强度、刚度以及稳定性设计方面还有一定的差距。 本文建立了波纹管非线性分析模型，并釆用大型有限商用软件ANSYS进行了计 算，计算结果与试验数据进行了比较，得到了与波纹管设计计算有关的载荷）应力 响应关系，并就波纹管非线性有限元计算进行研究。 第二章建模过程及求解 2.1波纹管结