! 5.4 圆柱形薄壳在均匀内压作用下的静力分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART ! 清除数据，不加载初始化文件
/FILNAME, EX5.4 ! 工程工件名称
/TITLE, EX5.4, CYLINDRICAL MEMBRANE UNDER PRESSURE.
! （1）定义工程选项、分析类型、单元类型、实常数和材料参数
/PREP7 ! 加载前处理模块
ANTYPE, STATIC ! 定义分析类型为静力分析
ET, 1, SHELL41 ! 定义第1类单元薄壳单元SHELL41
R, 1, 25.4 ! 定义厚度25.4 mm
MP, EX, 1, 207E3 ! 定义第1类材料的弹性模量EX=207 GPa
MP, NUXY, 1, 0.3 ! 定义第1类材料的泊松比NUXY=0.3
! （2）在柱坐标系下, 定义节点和单元
CSYS, 1 ! 切换为柱坐标系
/View, 1, 1, 1, 1 ! 将窗口1切换到等轴侧视点方向
N, 1, 1524 ! 在(1524mm, 0, 0)位置定义节点1
N, 2, 1524, , 254 ! 在(1524mm, 0, 254mm)位置定义节点
NGEN, 2, 2, 1, 2, 1, , 10 ! 按照1, 2号节点, 循环2次, 节点号加2, Y坐标加10度, 生成3, 4节点
NROTAT, ALL ! 旋转节点坐标系到柱坐标系方向
E, 1, 2, 4, 3 ! 过节点1, 2, 4, 3定义单元
! （3）定义耦合自由度、位移约束和荷载力
CP, 1, UX, 1, 2, 3, 4 ! 将1, 2, 3, 4号节点的UX径向自由度定义为第1个耦合自由度
CP, 2, UZ, 2, 4 ! 将2, 4号节点的UZ轴向自由度定义为第2个耦合自由度
D, 1, UZ, , , 3, 2 ! 位移约束1号节点的UZ自由度, 按增量2循环到3号节点
D, ALL, UY ! 位移约束所有节点的UY自由度
SFE, 1, 4, PRES, , -103.43*25.4 ! 定义1号单元的4号面上的轴向拉应力103.43 MPa
SFE, 1, 1, PRES, , -3.447575 ! 定义内压3.447575 MPa
FINISH ! 退出前处理模块PREP7
! (4)在求解模块, 设置求解选项后求解
/SOLU ! 加载求解模块
OUTPR, , 1 ! 输出第1荷载步的基本结果
SOLVE ! 开始求解
FINISH ! 退出求解模块SOLUTION
! (5)在后处理模块中提取应力计算结果
/POST1 ! 进入一般后处理模块
*GET, STRS_HOP, NODE, 1, S, 2 ! 提取1号节点的第2主应力, 存入变量STRS_HOP
*GET, STRS_AX, NODE, 1, S, 1 ! 提取1号节点的第1主应力, 存入变量STRS_AX
FINISH
