! 5.2 圆形薄膜大变形静力分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART ! 清除数据，不加载初始化文件
! （1）设置工程选项, 单元类型、实常数和材料参数
/FILNAME, EX5.2
/PREP7 ! 加载前处理模块
/TITLE, EX5.2, LARGE DEFLECTION OF A CIRCULAR MEMBRANE.
ET, 1, SHELL51, , , , 1 ! 定义第1类单元为轴对称结构壳体单元SHELL51
R, 1, 2.54e-3 ! 定义厚度为0.00254mm
MP, EX, 1, 207E3 ! 定义第1类材料的弹性模量EX＝207 GPa
MP, NUXY, 1, 0.25 ! 定义第1类材料的泊松比NUXY＝0.25
MP, ALPX, 1, 0.556E-5 ! 定义第1类材料的热膨胀系数ALPX＝0.556E-5
! (2)定义节点和单元
N, 1 ! 在坐标原点定义第1号节点
N, 11, 254 ! 在坐标(254mm, 0)定义第11号节点
FILL ! 在1号节点和11号节点之间填充生成其余节点
E, 1, 2 ! 用1号节点和2号节点定义单元
EGEN, 10, 1, 1 ! 按照前面单元循环10次, 生成其余的9个单元
FINISH ! 退出前处理模块PREP7
! (3)在求解模块, 设置求解选项, 施加位移约束后求解
/SOLU ! 加载求解模块
ANTYPE, STATIC ! 定义分析类型为静力分析
NLGEOM, ON ! 打开大变形选项
SSTIF, ON ! 打开应力强化选项
CNVTOL, F, , , , 2 ! 只使用力的收敛准则, 最小参考值为2N
OUTPR, , LAST ! 只输出最后一个荷载步的基本结果
TUNIF, -50 ! 设置温度变化, 形成温度应力
NSEL, S, LOC, X, 0 ! 选择X=0位置的节点
DSYM, SYMM, X ! 施加关于X轴(YZ平面)的对称边界条件
NSEL, ALL ! 选择所有节点
D, 11, ALL ! 固定11号节点的所有自由度
SOLVE ! 开始求解
! (4)设置递增步长荷载, 循环求解
KBC, 1 ! 设置固定步长荷载增量
SF, ALL, PRES, 689.5e-6 ! 在所有单元上施加压力荷载689.5Pa
SOLVE ! 开始求解
! (5)移去温度应力, 在求解
TUNIF, 0 ! 移去温度变化, 以删除温度应力
SOLVE ! 开始求解
! (6)在后处理模块, 定义单元表, 提取单元应力计算结果
/POST1 ! 进入一般后处理模块
SET, 3 ! 读入第3个荷载步的结果
*GET, UY, NODE, 1, U, Y ! 提取1号节点的位移UY值, 存入变量UY
ESEL, S, ELEM, , 1, 1 ! 选择1号单元
ETABLE, CENT, LS, 5 ! 以壳体中面薄膜应力“LS, 5”为内容定义单元表CENT
ESORT, CENT ! 对单元表CENT按大小排序
*GET, PRSCNT, SORT, , MAX ! 提取排序后的最大值, 存入变量PRSCNT
ESEL, S, ELEM, , 10, 10 ! 选择10号单元
ETABLE, CEN, LS, 5 ! 以壳体中面薄膜应力“LS, 5”为内容定义单元表CEN
ESORT, CEN ! 对单元表CEN按大小排序
*GET, PRSOUT, SORT, , MAX ! 提取排序后的最大值, 存入变量PRSOUT
*STATUS, PARM ! 显示保存在变量中的结果
FINISH
