ex2.1.inp

ansys有限元分析程序

! 2.1  铰接杆在外力作用下的变形
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除已有的数据
! （1）设置工程选项, 定义几何参数
/FilName, EX2.1	! 定义工程文件名称
/TITLE, EX2.1, DEFLECTION OF A HINGED SUPPORT
*AFUN, DEG 	! 将三角函数的默认单位由弧度改为角度
THETA=30  	! 定义变量Theta为杆件与水平面夹角30度
A=10000	! 定义变量A为节点3的水平X坐标10000mm
B=A/2*TAN(THETA)	! 节点2的铅直Y坐标
! （2）进入前处理模块, 定义节点
/PREP7	! 进入前处理模块
N, 1	! 在坐标原点(0, 0, 0)定义第1号节点
N, 2, A/2, -B	! 在 (A/2, -B) 定义第2号节点
N, 3, A	! 在 (A, 0) 定义第3号节点
! （3）在前处理模块中, 定义单元类型、材料参数和各个单元
ET, 1, LINK1	! 定义第1类单元为二维杆单元LINK1
R, 1, 1000	! 定义第1类实常数, 杆件截面面积A=1000mm^2
MP, EX, 1, 210E3	! 定义第1类材料的弹性模量EX＝210000N/mm^2
E, 1, 2	! 过节点1, 2定义单元1
E, 2, 3	! 过节点2, 3定义单元2
! （4）在前处理模块中, 定义位移约束
D, 1, ALL, , , 3, 2	! 约束1号节点的所有位移自由度, 按增量2循环到3来约束3号节点
F, 2, FY, -1000	! 在2号节点施加沿着-Y方向大小为1000N的集中力
OUTPR, , 1	! 输出第1个荷载步的基本计算结果
FINISH	! 退出前处理模块PREP7
! （5）进入求解模块, 开始求解
/SOLU	! 进入求解模块
SOLVE	! 开始静力分析求解
FINISH	! 退出求解模块SOLU
! （6）进入后处理模块, 显示节点位移和杆件内力
/POST1	! 进入后处理模块
PLDISP, 2	! 显示结构变形图, 参数“2”表示用虚线绘制出原来结构轮廓
MID_NODE = NODE (A/2, -B, 0 )	! 寻找距离 (A/2, -B, 0 )最近的节点, 存入MID_NODE
*GET, DISP, NODE, MID_NODE, U, Y	! 提取节点MID_NODE上的位移UY, 存入DISP
LEFT_EL = ENEARN (MID_NODE)	! 寻找距离节点MID_NODE最近的单元, 存入LEFT_EL
ETABLE, STRS, LS, 1	! 用轴向应力SAXL的编号”LS, 1”定义单元表STRS
*GET, STRSS, ELEM, LEFT_EL, ETAB, STRS	! 从单元表STRS中提取LEFT_EL单元的应力, 存入STRSS
*STATUS, PARM	! 显示所有内存变量的值
FINISH	! 退出后处理模块