ex3.11.inp

ansys有限元分析程序

! 3.11  偏心受压杆的大变形分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
/FILNAME, EX3.11	! 定义文件名称
! （1）定义工程选项、分析类型、实常数和材料参数等
/PREP7	! 加载前处理模块
/TITLE, EX3.11, LARGE DEFLECTION ECCENTRIC COMPRESSION OF SLENDER COLUMN.
ANTYPE, STATIC	! 定义分析类型为静力分析
NLGEOM, ON	! 激活大变形选择项
ET, 1, BEAM54	! 定义第1类单元为二维锥形非对称梁单元BEAM54
R, 1, 2167.7, 54.106e4, 14.732, 42.671	! 定义第1类实常数：截面积2167.7mm^2，抗弯惯性矩54.106cm^4，中性轴到上下边的距离分别为14.732mm和42.671mm
RMORE, , , , -14.732, , -14.732	! 两端偏心距e=14.732mm（HYT2, HYB2, DX1, DY1, DX2, DY2）
MP, EX, 1, 207E3	! 定义第1类材料的弹性模量EX = 207e3 N/mm^2
MP, NUXY, 1, 0.3	! 定义第1类材料的泊松比NUXY = 0.3
! （2）定义节点
N, 1	! 在坐标原点(0, 0, 0)定义第1号节点
N, 5, , 1524	! 在坐标(0, 1524mm, 0)定义第5号节点
FILL	! 在第1号节点和第5号节点中间填充生成其余节点
! （3）定义单元
E, 1, 2	! 过节点1和节点2定义单元
EGEN, 4, 1, 1	! 用前面一个单元的模式，循环4次，生成其余的3个单元
! （4）施加位移约束和荷载
D, 1, ALL	! 完全固定1号节点
F, 5, FY, -17793	! 定义沿着-Y轴的大小为17793N的集中力
FINISH	! 退出前处理模块PREP7
! （5）进入求解模块求解
/SOLU	! 加载求解模块
CNVTOL, F, , 1E-4	! 设置非线性分析力的收敛误差1E-4 N
CNVTOL, M, , 1E-1	! 设置非线性分析力偶的收敛误差为1E-1 N.m
OUTPR, , LAST	! 设置只输出最后一个荷载步的结果
SOLVE	! 开始大变形分析的非线性迭代求解
FINISH	! 退出求解模块SOLUTION
! （6）在后处理模块，提取计算结果
/POST1	! 进入一般后处理模块
END_NODE = NODE (0, 3048/2, 0)	! 选择距离(0, 1524mm, 0)最近的节点, 命名为END_NODE
*GET, DEF, NODE, END_NODE, U, X	! 提取节点END_NODE的位移UX, 存入变量DEF
ETABLE, TENS, NMISC, 1	! 以BEAM54的最大拉应力(NMISC, 1)为内容定义单元表TENS
ETABLE, COMP, NMISC, 2	! 以BEAM54的最大压应力(NMISC, 2)为内容定义单元表COMP
*GET, STS_TENS, ELEM, 1, ETAB, TENS	! 提取第1号单元的最大拉应力, 存入变量STS_TENS
*GET, STS_COMP, ELEM, 1, ETAB, COMP	! 提取第1号单元的最大压应力, 存入变量STS_COMP
*STATUS, PARM
FINISH