ex4.9.inp

ansys有限元分析程序

! 4.9  用二维实体单元分析等截面悬臂梁的平面弯曲
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH	! 退出以前模块
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
! (1)设置工程选项、单元类型和材料参数
/FILNAME, EX4.9	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX4.9, BENDING OF A SOLID BEAM (PLANE ELEMENTS).
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, PLANE42, , , 3, , , 2	! 定义第1类单元为带厚度二维结构固体单元PLANE42，带1, 3号面输出
MP, EX, 1, 207E3	! 定义第1类材料的弹性模量EX = 207GPa =207 N/mm^2
MP, NUXY, 1, 0.0	! 定义第1类材料的泊松比NUXY ＝ 0.0
R, 1, 25.4	! 定义单元厚度为25.4mm
! (2)定义节点
N, 1	! 在坐标原点(0, 0, 0)定义第1号节点
N, 6, 254	! 在(254mm, 0, 0)定义6号节点
FILL	! 在1号节点和6号节点之间填充生成2号节点到5号节点
! NGEN, ITIME, INC, NODE1, NODE2, NINC, DX, DY, DZ, SPACE
NGEN, 2, 10, 1, 6, 1, , 50.8	! 循环2次, 节点号增量10, 按1号节点到6号节点增量为1的范围, 每
个节点的Y坐标增量为50.8mm, 生成节点11到16
E, 1, 2, 12, 11	! 过节点1, 2, 12和11定义单元
EGEN, 5, 1, 1	! 循环5次, 每次单元两端的节点号增加1, 生成其余的4个单元
! (3)施加第1种情况的边界条件:端部力偶作用情况
D, 1, ALL, , , 11, 10	! 位移约束1号节点的所有位移自由度, 按增量10循环到11号节点.
F, 6, FX, 4448.3	! 在6号节点的X方向施加4448.3N的集中力
F, 16, FX, -4448.3	! 在16号节点的-X方向施加4448.3N的集中力, 和6号节点的力组成力偶
OUTPR, NSOL, 1	! 输出第1荷载步的节点计算结果
OUTPR, ESOL, 1	! 输出第1荷载步的单元计算结果
FINISH	! 退出前处理模块
! (4)进入求解模块求解
/SOLU	! 加载求解模块
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (5)在后处理模块, 提取节点位移计算结果
/POST1	! 进入一般后处理模块
*GET, U1, NODE, 16, U, Y	! 提取16号节点的位移UY, 存入变量U1
PLNSOL, S, X	! 图形显示X方向的正应力
FINISH	! 退出后处理模块
! (6)改变荷载为第二种荷载工况:端部力作用情况，重新计算
/SOLU	! 加载求解模块
F, 6, FX, , , 16, 10	! 在6号节点X方向施加大小等于0的集中力，按增量10循环到16号节点
F, 6, FY, 667.25, , 16, 10	! 在6号节点Y方向施加大小等于667.25N的集中力，按增量10循环到16号节点
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLUTION
! （7）在后处理模块, 提取节点位移计算结果
/POST1	! 进入一般后处理模块
*GET, U2, NODE, 16, U, Y	! 提取16号节点的位移UY, 存入变量U2
PLNSOL, S, X	! 图形显示X方向的正应力
FINISH	! 退出后处理模块