ex15.2.txt

通过实例详细介绍Ansys在实际工程中的应用

!15.4  利用8节点2D实体单元PLANE82单元计算长柱形天然气罐
!15.4.4  源程序
! 本程序来源于邢静忠等编著的《ANSYS分析实例与工程应用》，机械工业出版社（2004年）
Finish
/Clear
!(1)定义基本参数
h1=298.5         ! 封头的高度坐标298.5 mm
b1=44.5          ! 封头的内半径44.5 mm
b2=b1+23         ! 封头的外半径44.5+23mm厚度
r1=185           ! 罐体内半径185mm
r2=203           ! 罐体外半径203mm
!(2)定义剖面区域边界上的关键点
/Title, Long Cyclinder Pressure Vessel for CNG with PLANE82.
/Prep7           ! 进入前处理模块
k,1,b1,h1        ! 封头起始端内侧的点
k,2,b2,h1        ! 封头起始端外侧的点
k,3,b1,sqrt(r1*r1-b1*b1)     ! 封头内侧剖面线和罐体半球内剖面线的交点
k,4,b2,sqrt(r2*r2-b2*b2)     ! 封头外侧剖面线和罐体半球外剖面线的交点
k,5,r1,0            ! 罐体半球和罐体主体的内侧交点
k,6,r2,0            ! 罐体半球和罐体主体的外侧交点
k,7,0,0             ! 坐标原点
k,8,r1,-1000+h1     ! 罐体主体截断内侧点
k,9,r2,-1000+h1     ! 罐体主体截断外侧点
K,10,0,h1           ! 封头起始端对称线上的点
!(3)定义围成结构剖面区域的线，并定义剖面
L,1,2              ! 封头起始端水平线
L,1,3              ! 封头内侧竖直线
L,2,4              ! 封头外侧竖直线
L,5,8              ! 罐体主体内侧线
L,8,9              ! 罐体主体截断位置水平线
L,6,9              ! 罐体主体外侧线
Larc,5,3,7,r1    ! 过5号，3号关键点，以7号关键点为圆心，r1为半径做半球内圆弧线
Larc,6,4,7,r2    ! 过6号，4号关键点，以7号关键点为圆心，r2为半径做半球外圆弧线
LFILLT,2,7,63, ,     ! 过2号，7号线，以半径63mm做半球内侧位置的倒圆
LFILLT,3,8,40, ,     ! 过3号，8号线，以半径40mm做半球外侧位置的倒圆
AL,all               ! 将所有线围成的区域定义为一个面
!(4)定义单元类型、材料参数和实常数
ANTYPE,STATIC       ! 定义分析类型为静力分析
ET,1,PLANE82        ! 定义第1类单元是8节点二维实体单元(用轴对称选项)
KEYOPT,1,3,1        ! 定义单元的第3个可选项，打开轴对称选项
MP,EX,1,207E3       ! 定义第1类材料的弹性模量为207 GPa
MP,NUXY,1,0.3       ! 定义第1类材料的泊松比NUXY=0.3
!(5)定义位移约束和荷载 
DL,5,,SYMM          ! 给5号线定义对称位移条件
SFL,2,PRES,23, , , ,     ! 给剖面线的内侧线2施加内压荷载23MPa
SFL,4,PRES,23, , , ,     ! 给剖面线的内侧线4施加内压荷载23MPa
SFL,7,PRES,23, , , ,     ! 给剖面线的内侧线7施加内压荷载23MPa
SFL,9,PRES,23, , , ,     ! 给剖面线的内侧线9施加内压荷载23MPa
SFL,1,PRES,-17.681,     ! 给端部封头水平线施加封头传递过来的-17.681MPa
Esize,5                 ! 指定单元尺寸大小为5mm
Amesh,1                 ! 将罐体剖面1号面剖分为面单元
FINISH                  ! 退出前处理模块
!(6)进入求解模块求解
/SOLU              ! 进入求解模块
Solve            ! 开始求解
FINISH           ! 退出求解模块
!(7)进入后处理模块处理计算结果
/POST1          ! 进入后处理模块
Pldisp,2        ! 显示变形图
Plnsol,s,1      ! 显示最大拉应力分布
Plnsol,s,2      ! 显示第2主应力分布
Plnsol,s,3      ! 显示最大压应力分布
Plnsol,s,EQV    ! 显示最大Mises应力分布
Plnsol,s,Z       ! 显示环向应力分布