ex5.2.inp

ansys有限元分析程序

! 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