ex5.9.inp

来自「ansys有限元分析程序」· INP 代码 · 共 61 行

! 5.9  两端简支开口柱壳在自重作用下的静力分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH	! 退出以前模块
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! (1)设置工程选项、分析类型、实常数、材料参数
/FILNAME, EX5.9	! 定义工程文件名称
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化网格功能
/TITLE, EX5.9, BARREL VAULT ROOF UNDER SELF-WEIGHT.
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, SHELL93	! 定义第1类单元是8节点壳单元SHELL93
R, 1, 0.25	! 定义壳体厚度为0.25m
MP, EX, 1, 4.32E8	! 定义第1类材料的弹性模量EX ＝ 432 MPa
MP, NUXY, 1, 0.0	! 定义第1类材料的泊松比NUXY＝0.0
MP, DENS, 1, 36.7347	! 定义第1类材料的密度DENS ＝ 36.7347 kg/m^3
! (2)在柱坐标系下定义关键点、柱面并对其进行单元剖分
CSYS, 1	! 切换坐标系到柱坐标系
K, 1, 25, 50	! 在坐标(25m, 50度)定义第1个关键点
K, 2, 25, 50, 25	! 在坐标(25m, 50度, 25m)定义第2个关键点
KGEN, 2, 1, 2, 1, , 40	! 循环2次, 按照1号节点到2号节点范围, 角度坐标增量40度, 生成3, 4节点
A, 1, 3, 4, 2	! 过关键点1, 3, 4和2定义柱面
ESIZE, , 4	! 定义单元剖分尺寸为分为4段
/View, 1, 1, 1, 1	! 将窗口1的视点切换到等轴测方向
AMESH, 1	! 对1号面执行面单元剖分
! (3)选择对称边界上的点, 施加对称边界条件, 对端部边界点施加墙体约束
CSYS, 0	! 切换坐标系到整体直角坐标系
NSEL, S, LOC, X	! 选择位于X＝0位置的节点
DSYM, SYMM, X	! 对选择集中的所有节点施加关于X轴(YZ平面)的对称位移条件
NSEL, S, LOC, Z	! 选择位于Z＝0位置的节点
DSYM, SYMM, Z	! 对选择集中的所有节点施加关于Z轴(XY平面)的对称位移条件
NSEL, S, LOC, Z, 25	! 选择位于Z＝25m位置的节点
D, ALL, UX, 0, , , , UY, ROTZ	! 约束Z＝25m位置的节点的UX, UY和ROTZ自由度
NSEL, ALL	! 选择所有节点
ACEL, , 9.8	! 定义Y方向的重力加速度ACEL＝9.8 m/s^2
FINISH	! 退出前处理模块PREP7
! (4)进入求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始静力分析的求解
FINISH	! 退出求解模块SOLU
! (5)进入后处理模块提取计算结果
/POST1	! 进入一般后处理模块
NSEL, S, NODE, , 1, 2, 1	! 选择上面A点和B点位置的节点1和2
ESLN, S    	! 选择包含1号节点和2号节点的单元
PRNSOL, U, COMP	! 列表显示位移计算结果及其合成位移结果
*GET, UYA, NODE, 1, U, Y	! 提取1号节点的Y方向位移UY, 并保存到变量UYA
*GET, UXA, NODE, 1, U, X	! 提取1号节点的X方向位移UX, 并保存到变量UXA
! (6)在柱坐标系下显示, 提取计算结果
RSYS, 1	! 将结果坐标系切换到柱坐标系下
SHELL, TOP	! 选择壳体的上表面作为计算结果输出对象
PRNSOL, S, COMP	! 列表显示应力分量计算结果
*GET, SIGZ_TOP, NODE, 1, S, Z	! 提取1号节点的Z方向应力SZ, 并保存到变量SIGZ_TOP
*GET, SIGY_TOP, NODE, 2, S, Y	! 提取2号节点的Y方向应力SY, 并保存到变量SIGY_TOP
! (7)提取壳体下表面的应力计算结果
SHELL, BOT	! 选择壳体的下表面作为计算结果输出对象
PRNSOL, S, COMP	! 列表显示应力分量计算结果
*GET, SIGZ_BOT, NODE, 1, S, Z	! 提取1号节点的Z方向应力SZ, 并保存到变量SIGZ_BOT
*GET, SIGY_BOT, NODE, 2, S, Y	! 提取2号节点的Y方向应力SY, 并保存到变量SIGY_BOT
*STATUS, PARM	! 显示保存在变量中的结果
FINISH