波纹管ANSYS有限元分析

第一章综述

1.1波纹管的发展

波纹管膨胀节是压力容器和管道系统中常用的连接和补偿装置, 其主体是金属波纹管。它的特点是轴向刚度较小, 径向刚度较大, 故能够承受较大的轴向位移,同时还能够承受一定的内、外压, 常用于补偿管道因温度变化、机械位移或振动引起的轴向、横向和角向位移。

运行中的波纹管的波峰波谷部位实际处于弹塑性区域, 但是, 目前波纹管常

用的设计标准是将波纹管视为梁或曲杆, 通过梁理论进行弹性的应力应变分析,

这势必会带来较大的误差。随着数值分析技术和有限元技术的不断发展, 利用非线性有限元技术对波纹管进行较为精确的应力应变分析已成为波纹管研究领域的一

项重要课题。

本文吸取了前人在有限元程序设计上的经验, 根据波纹管自身的非线性特性, 利用Fortran Pow erSta-tion 410与Visual C++ 610混合语言编程技术, 开发出一

套界面友好, 使用方便, 求解精度较高的波纹管有限元应力应变分析专用软件系统。

1.2波纹管非线性特性与有限元程序设计

固体力学问题, 从本质上讲都是非线形的, 线性假设只是实际问题的一种简化。线性分析的基本假设包含三点: 结点位移为无限小量, 材料为线弹性, 加载时边界条件的性质保持不变。如果上述假设中的任一条不能满足, 则属于非线性问题。常见的是不满足前两种假设的非线性行为, 分别称为几何非线性和材料非线性。1.3波纹管的非线性特性

工作状态下的波纹管伸缩变形是一种复杂的非线性变形。在小载荷情况下, 波纹管的应力应变关系呈线性特征, 而当其承受较大的位移补偿时, 结构处于高应力水平, 波峰和波谷的局部处于塑性区, 此时, 一般采用弹性线性强化模型分析波纹管的应力应变关系。另外, 波纹管的单波位移量一般大于壁厚, 属于大位移行为, 同时, 大量的研究表明, 此时的波纹管大部分区域依然处于小应变范围。因此, 运行状况下的波纹管还表现出大位移、小应变的几何非线性特征。

1.3.1非线性有限元程序设计

本软件系统所用的求解程序是一个通用的二维非线性有限元程序, 图 1 是该

程序的流程图, 其主体部分借鉴了文献132中提出的有限元程序设计方法。

图1.1 二维非线性有限元程序流程

1.3.2波纹管有限元计算模型

根据波纹管的几何特征和载荷特征, 系统采用二维轴对称模型进行求解, 单元类型为 8 结点抛物线等参单元, 并开发出相应的前处理程序。以 U 型波纹管为例, 通过输入波纹管的基本波形参数、材料参数和运行状况参数, 便可建立如图 2 所示的计算模型, 自动进行网格划分, 边界约束和等效结点力加载。

图1.2 波纹管二维轴对称有限元计算模型

波纹膨胀节被广泛应用于管路热补偿、隔振、降噪。自 40 年代以来, 世界各工业国家的研究机构和生产厂家对作为其主要组成部分的波纹管进行了广泛的研究。

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

70 年代以后, 随着计算机技术的飞速发展,有限元分析在结构分析中的地位日渐突出, 许多商用有限元软件应运而生。有限元分析以其在解决几何非线性、材料非线性和结构非线性问题方向的独特优势而为愈来愈多的研究人员所接受。非线性有限元计算能够有效地解决传统的解析法无法计算的波纹管弹塑性大变形范围内的载荷) 应力响应问题。但由于以往有限元计算软件方面的局限性, 其分析结果距离能够真正运用到波纹管的强度、刚度以及稳定性设计方面还有一定的差距。

本文建立了波纹管非线性分析模型, 并采用大型有限商用软件ANSYS 进行了计算, 计算结果与试验数据进行了比较, 得到了与波纹管设计计算有关的载荷) 应力响应关系, 并就波纹管非线性有限元计算进行研究。

第二章建模过程及求解

2.1波纹管结构示意图

图2.1 参数含义

2.2第一组数据

表2-1

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.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建立坐标关键点

ANSYS Main Menu:选择Preprocessor →Modeling →Create →Keypoints →In Active CS填写对应坐标

图2.2 坐标点图

2.2.3生成特征点并且连接、倒圆角、布尔操作

选择Preprocessor|Modeling|Create|Lines|Straight Line 弹出Create StraightLine 框，连接各点。如下图所示

选择Preprocessor|Modeling|Operate|Booleans|Add|Lines

图2.3 坐标点连线

2.2.4生成壳体

选择Preprocessor|Modeling|Operate|Booleans|Add|Lines