波纹管ANSYS有限元分析.docx
《波纹管ANSYS有限元分析.docx》由会员分享，可在线阅读，更多相关《波纹管ANSYS有限元分析.docx（14页珍藏版）》请在冰豆网上搜索。
波纹管ANSYS有限元分析
第一章综述
1.1波纹管的发展
波纹管膨胀节是压力容器和管道系统中常用的连接和补偿装置，其主体是金属波纹管。
它的特点是轴向刚度较小，径向刚度较大，故能够承受较大的轴向位移，同时还能够承受一定的内、外压，常用于补偿管道因温度变化、机械位移或振动引起的轴向、横向和角向位移。
运行中的波纹管的波峰波谷部位实际处于弹塑性区域，但是，目前波纹管常
用的设计标准是将波纹管视为梁或曲杆，通过梁理论进行弹性的应力应变分析，这势必会带来较大的误差。
随着数值分析技术和有限元技术的不断发展，利用非线
性有限元技术对波纹管进行较为精确的应力应变分析已成为波纹管研究领域的一项重要课题。
本文吸取了前人在有限元程序设计上的经验，根据波纹管自身的非线性特性，利用FortranPowerStation410与VisualC++610混合语言编程技术,开发出一套界面友好，使用方便，求解精度较高的波纹管有限元应力应变分析专用软件系统。
1.2波纹管非线性特性与有限元程序设计
固体力学问题，从本质上讲都是非线形的，线性假设只是实际问题的一种简化。
线性分析的基本假设包含三点：
结点位移为无限小量，材料为线弹性，加载时边界条件的性质保持不变。
如果上述假设中的任一条不能满足，则属于非线性问题。
常见的是不满足前两种假设的非线性行为，分别称为几何非线性和材料非线性。
1.3波纹管的非线性特性
工作状态下的波纹管伸缩变形是一种复杂的非线性变形。
在小载荷情况下，波
纹管的应力应变关系呈线性特征，而当其承受较大的位移补偿时，结构处于高应力水平，波峰和波谷的局部处于塑性区，此时，一般采用弹性线性强化模型分析波纹管的应力应变关系。
另外，波纹管的单波位移量一般大于壁厚,属于大位移行为，同时，大量的研究表明，此时的波纹管大部分区域依然处于小应变范围。
因此，运行状况下的波纹管还表现出大位移、小应变的几何非线性特征。
1.3.1非线性有限元程序设计
本软件系统所用的求解程序是一个通用的二维非线性有限元程序，图1是该
程序的流程图，其主体部分借鉴了文献
132中提出的有限元程序设计方法。
AUSOR强股变量以以别将法英5：
ISCREM投给定的赖荷系数塘量加载
STIFTPHI）用元刚度矩阵
CXJTPLT给出本次莪黄增量下的结果
图1.1二维非线性有限元程序流程
1.3.2波纹管有限元计算模型
根据波纹管的几何特征和载荷特征，系统采用二维轴对称模型进行求解，单元类型为8结点抛物线等参单元，并开发出相应的前处理程序。
以U型波纹管为例，通过输入波纹管的基本波形参数、材料参数和运行状况参数，便可建立如图2所示的计算模型，自动进行网格划分，边界约束和等效结点力加载。
图1.2波纹管二维轴对称有限元计算模型
波纹膨胀节被广泛应用于管路热补偿、隔振、降噪。
自40年代以来，世界各
工业国家的研究机构和生产厂家对作为其主要组成部分的波纹管进行了广泛的研
究。
在不同工况下对波纹管进行应力）应变的分析是评估其对各种工况承受能力的基础和重要组成部分。
70年代之前的结构分析主要采用解析法，根据近似简单梁、近似圆柱体、近似壳体的假设，依据弹性理论得到波纹管壳体中的应力）应变值[1]。
但由于波纹管本身是一种较为复杂的轴对称薄壁壳体，且在绝大多数工况下材料处于弹塑性大变形范围内，因而在解析解与波纹管材料的实际响应之间存在着较大的误差。
为使两者间达到近似的统一，学者们在解析法得到的近似解中引入了由图表形式给出的修正系数。
这种方法目前正广泛应用于设计以及生产部门，如EJMA
规范中的应力计算式[2]o
70年代以后，随着计算机技术的飞速发展，有限元分析在结构分析中的地位日渐突出，许多商用有限元软件应运而生。
有限元分析以其在解决几何非线性、材料非线性和结构非线性问题方向的独特优势而为愈来愈多的研究人员所接受。
非线性有限元计算能够有效地解决传统的解析法无法计算的波纹管弹塑性大变形范围内的载荷）应力响应问题。
但由于以往有限元计算软件方面的局限性，其分析结果距离
能够真正运用到波纹管的强度、刚度以及稳定性设计方面还有一定的差距。
本文建立了波纹管非线性分析模型,并采用大型有限商用软件ANSYS!
行了计算，计算结果与试验数据进行了比较，得到了与波纹管设计计算有关的载荷）应力响应关系，并就波纹管非线性有限元计算进行研究。
第二章建模过程及求解
3.1波纹管结构示意图
JRL
一.-h—疝s.—■aa
rr\gf—
I_uuu
图2.1参数含义
3.2第一组数据
表2-1
外外径内径壁厚波距波高波纹管波峰半波谷半层波
参
聂DD2STB长度L径R1径R2数数
数
（mm）（mm）（mm）（mm）（mm）（mm）（mm）（mm）（m）（m）
数
二7736981907544023.721.314
值
3具体坐标
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）
3建立坐标关键点
ANSYSMainMenu:
选择Preprocessor-Modeling-Create-Keypoints一In
ActiveCS填写对应坐标
KinjuslTf«MUH
图2.2坐标点图
3生成特征点并且连接、倒圆角、布尔操作
选择Preprocessor|Modeling|Create|Lines|StraightLineW出CreateStraightLine框，连接各点。
如下图所示
选择Preprocessor|Modeling|Operate|Booleans|Add|Lines
ATJS\S|
图2.3坐标点连线
3生成壳体
选择Preprocessor|Modeling|Operate|Booleans|Add|Lines
图2.4螺线管壳体
3划分网格
ANSYSMainMenu:
Preprocessor-Meshing-MeshTooH（SizeControls）lines
图2.5网格划分
3分析计算
选择Solution|Solve|CurrentLS得出结果
3偏移量
图2.6偏移图
3分析计算
图2.7应力云图
这个是应力云图，从图2.7中可以看出当最大位移是45mmi勺时候最大应应力是491mpa最小应力是3.3281mpa.最大应力接近屈服极限480mpa,因此可以认为波纹管最大补偿是45mm
这个是位移图2.8,当位移是45mm的时候达到304材料的最大屈服极限
3.3第二组数据
表2-2
参
数
外径
D
（mm）
内径
D2
（mm）
壁厚
S
（mm）
波距
T
（mm）
波高
B（mm）
波纹管长度L（mm）
波峰半径R1（mm）
波谷半
径R2（mm）
层数（m）
波数（m）
数值
723
648
1
90
75
440
23.7
21.3
1
4
3
具体坐标
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.-324.0）20（440.-324.0）
3定义材料参数
因为我们选的材料是304，上网查的304钢所对应的弹性模量是1.9e11pa,泊松比是0.305
3建立坐标关键点
ANSYSMainMenu:
选择Preprocessor-Modeling-Create-Keypoints一In
ActiveCS填写对应坐标
图2.9坐标点图
3生成特征点并且连接、倒圆角、布尔操作
选择Preprocessor|Modeling|Create|Lines|StraightLine单出CreateStraightLine框，连接各点。
如下图所示
选择Preprocessor|Modeling|Operate|Booleans|Add|Lines
图2.10坐标连线
3生成壳体
选择Preprocessor|Modeling|Operate|Booleans|Add|Lines
图2.11螺纹管壳体
3划分网格
ANSYSMainMenu:
Preprocessor-Meshing-MeshTooH（SizeControls）lines
图2.12网格划分
2.2.7分析计算
选择Solution|Solve|CurrentLS得出结果
偏移量
图2.13偏移图
应力及波纹处受力分析
图2.14应力云图
这个是应力云图，从图
2.14中可以看出当最大位移是43mm的时候最大应应力
是436.04mpa,最小应力是9.1958mpa.最大应力接近屈服极限480mpa,因此可以认为波纹管最大补偿是43mm
图2.15位移云图
这个是位移云图2.15,当位移是43mm的时候达到304材料的最大屈服极限此时的位移补偿量为43mm
第三章结论说明
通过上面两组数据的对比分析，得由以下结论
（1）内径对波纹管的补偿量是很大的，
第一组数据内外径D=773mmd=698波纹管的补偿量是45mm勺时候它的应力接近最大屈服极限480mpa;
第二组数据内外径D=723mmd=648波纹管的补偿量是43mm勺时候它的应力接近最大屈服极限480mpa.
内外径相差50mm波纹管的补偿量就相差2mm由此可以看出内外径对波纹管的影响是不可忽略的
（2）波纹管通常要求有尽可能大的补偿量，因此波纹管在满足承压的前提下，尽可能减小内径，以获得较高的柔性。
参考文献
[1]叶青.微机原理与接口技术.清华大学出版社，2011:
9-1.
[2]薛华成.管理信息系统[M].北京：
清华大学出版社，1993:
99-103.
[3]李晖.高性+能计算机若干关键问题研究[D].中国科学技术大学，2009:
125-138.
沈阳化工大学科亚学院
工程软件实践应用
专业：
机械设计制造及其自动化班级:
学生姓名：
指导教师:
工程软件应用实践任务书
机械制造及其自动化专业班学生:
论文题目：
内外径对U形波纹管膨胀节补偿能力的影响
具体参数:
参数1
参数2
外径Di（mm）
773
723
内径D2（mm）
698
648
壁厚（mm）
1
1
波距T（mm）
90
90
波高B（mm）
70
70
波纹管长度L（mm）
440
440
波峰半径R1（mm）
23.7
23.7
波谷半径R2（mm）
21.3
21.3
层数m
1
1
波数n
4
4
材料
304
304
指导教师：
201年月日
第一章综述1
波纹管的发展1
波纹管非线性特性与有限元程序设计1
波纹管的非线性特性1
非线性有限元程序设计2
波纹管有限元计算模型2
第二章建模过程及求解4
波纹管结构示意图4
第一组数据4
具体坐标4
建立坐标关键点4
生成特征点并且连接、倒圆角、布尔操作5
生成壳体5
划分网格6
分析计算6
偏移量6
分析计算7
第二组数据8
具体坐标8
定义材料参数8
建立坐标关键点8
生成特征点并且连接、倒圆角、布尔操作9
生成壳体9
划分网格10
分析计算10
偏移量10
应力及波纹处受力分析11
第三章结论说明13
3.1通过上面两组数据的对比分析，得出以下结论13
参考文献14
升级会员 