ex5.5.inp

ansys有限元分析程序

! 5.5  简支和固支圆板的在不同荷载作用下的弯曲
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
! （1）设置工程参数、分析类型、单元类型、实常数和材料参数
/FILNAME, EX5.5
/PREP7	! 加载前处理模块
/TITLE, EX5.5, BENDING OF A CIRCULAR PLATE USING AXISYMMETRIC SHELL ELEMENTS.
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, SHELL51	! 定义第1类单元为轴对称壳体单元SHELL51
R, 1, 25.4	! 定义单元厚度为25.4mm
MP, EX, 1, 207E3	! 定义第1类材料的弹性模量为EX=207 GPa
MP, NUXY, 1, 0.3	! 定义第1类材料的泊松比NUXY=0.3
! （2）定义关键点、线和面, 并剖分为面单元
K, 1	! 在坐标原点定义第1个关键点
K, 2, 1016	! 在坐标(1016mm, 0)定义第2个关键点
K, 3, 508	! 在坐标(508mm, 0)定义第3个关键点
L, 1, 3	! 过关键点1, 3定义直线
L, 3, 2	! 过关键点3, 2定义直线
LESIZE, 1, , , 5, 20	! 将1号线剖分为5段, 间隔比例按20递增
LESIZE, 2, , , 5, .05	! 将2号线剖分为5段, 间隔比例按0.05递减
LMESH, ALL	! 将所有线剖分为线单元
/PNUM, ELEM, 1	! 打开单元编号显示开关
EPLOT	! 图形显示单元布置
FINISH	! 退出前处理模块PREP7
! (3)按第1种情况计算:均布荷载、固支端
/SOLU	! 加载求解模块
OUTPR, , 1	! 输出第1荷载步的基本结果
DK, 1, UX, , , , ROTZ	! 约束第1个关键点的X方向自由度和Z方向转动自由度
DK, 2, ALL  	! 约束第2个关键点的所有自由度
SFE, ALL, 1, PRES, , 0.041370	! 在所有线表面施加0.041370 MPa的压力荷载
SOLVE	! 开始求解
! (4)按第2种情况计算:中心集中荷载、固支边
FK, 1, FY, 33539	! 在第1个关键点上施加沿着-Y方向的集中力33539 N
SFE, ALL, 1, PRES, , 0	! 删除表面压力荷载
SOLVE	! 开始求解
! （5）按第3种情况计算:均布荷载、简支边
DKDELE, 2, ROTZ	! 删除夹紧边界条件约束ROTZ=0
FK, 1, FY	! 定义第1个关键点上的力为0，这样就删去了前面在这个关键点上的定义
SFE, ALL, 1, PRES, , 0.010343	! 在所有线表面施加0.010343 MPa的压力荷载
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (6)在后处理模块中显示变形情况
/POST1	! 进入一般后处理模块
/DSCALE, 1, 35	! 定义位移放大系数35
/WINDOW, 1, -1, 1, -1, -.333	! 定义1号窗口
SET, 1, 1	! 读入第1个荷载步的计算结果
PLDISP, 1	! 带原来模型图, 显示变形后的模型
MID_NODE = NODE (0, 0, 0)	! 选择距离(0, 0, 0)最近的节点, 存入变量MID_NODE
*GET, DEF_C1, NODE, MID_NODE, U, Y	! 提取节点MID_NODE的位移UY, 存入变量DEF_C1
ETABLE, STRS, LS, 9	! 以壳体底部的拉应力”LS, 9”为内容, 定义单元表STRS
*GET, STRSS_C1, ELEM, 10, ETAB, STRS	! 从单元表STRS提取单元10的应力, 存入变量STRSS_C1
! (7)读入第2荷载步的计算结果, 显示变形图
/WINDOW, 1, OFF	! 关闭1号窗口的显示
/WINDOW, 2, -1, 1, -.333, .333, 1	! 定义2号窗口
/NOERASE	! 不擦除已有窗口的内容
SET, 2, 1	! 读入第2荷载步的计算结果
PLDISP	! 显示第2荷载步的变形图
*GET, DEF_C2, NODE, MID_NODE, U, Y	! 提取节点MID_NODE的位移UY, 存入变量DEF_C2
ETABLE, STRS, LS, 9	! 以壳体底部的拉应力”LS, 9”为内容, 定义单元表STRS
*GET, STRSS_C2, ELEM, 10, ETAB, STRS	! 从单元表STRS提取单元10的应力, 存入变量STRSS_C2
! （8）读入第3荷载步的计算结果, 显示变形图
/WINDOW, 2, OFF	! 关闭2号窗口
/WINDOW, 3, -1, 1, 0.333, 1, 1	! 定义3号窗口
SET, 3, 1	! 读入第3荷载步的计算结果
PLDISP	! 显示第3荷载步的变形图
*GET, DEF_C3, NODE, MID_NODE, U, Y	! 提取节点MID_NODE的位移UY, 存入变量DEF_C3
ETABLE, STRS, LS, 9	! 以壳体底部的拉应力”LS, 9”为内容, 定义单元表STRS
*GET, STRSS_C3, ELEM, 1, ETAB, STRS	! 从单元表STRS提取单元1的应力, 存入变量STRSS_C3
*STATUS, PARM	! 显示保存在变量中的结果
FINISH