ex3.4.inp

ansys有限元分析程序

! 3.4  工字形截面外伸梁的平面弯曲
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
Finish
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
/FILNAME, EX3.4	! 定义文件名称
/PREP7	! 加载前处理模块
/TITLE, EX3.4, BEAM STRESSES AND DEFLECTIONS.	! 设置标题
! （1）定义求解类型，实常数和材料参数
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, BEAM3	! 定义第1类单元为二维弹性梁单元BEAM3
KEYOPT, 1, 9, 9	! 在9个中点位置上输出计算结果
R, 1, 35e-4, 4.5E-5, 0.20	! 定义第1类实常数、截面积35cm^2、惯性矩4500cm^4和高度20cm
MP, EX, 1, 210E9	! 定义第1类材料的弹性模量EX＝210Gpa
MP, PRXY, 1, 0.3	! 定义第1类材料的泊松比PRXY=0.3
! （2）定义节点位置
N, 1	! 在坐标原点(0，0)定义1号节点
N, 5, 4	! 在(4m，0)位置定义5号节点
FILL	! 利用节点填充FILL生成中间3个节点
! （3）定义单元连接关系
E, 1, 2	! 过1，2节点定义第1个单元
EGEN, 4, 1, 1	! 按照前面的单元模式，生成4次，每次两端节点编号加1，生成其余3个单元
! （4）定义位移约束和荷载条件
D, 2, UX, , , , , UY	! 定义2号节点的位移边界条件，UX=0, UY=0
D, 4, UY	! 定义4号节点的位移边界条件，UY=0
SFBEAM, 1, 1, PRES, 5000	! 定义第1个单元的表面荷载5000N/m
SFBEAM, 4, 1, PRES, 5000	! 定义第4个单元的表面荷载5000N/m
FINISH	! 退出前处理模块PREP7
! （5）进入求解模块，定义求解选项并求解
/SOLU	! 加载求解模块
OUTPR, BASIC, 1	! 只输出第1个荷载步的基本计算结果
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! （6）进入后处理模块，列表显示节点位移计算结果和图形显示变形图
/POST1	! 进入一般后处理模块
SET, 1, 1	! 读入第1个荷载步的计算结果
PRNSOL, U, COMP	! 列表显示节点线位移和节点线位移的矢量和
PRNSOL, ROT, COMP	! 列表显示节点转角和节点角位移的矢量和
PLDISP, 1	! 带原来模型图，显示变形后的模型
! （7）提取指定位置的节点和单元，定义单元表并获取弯曲应力
MID_NODE = NODE (2, , , )	! 选择距离(2, 0, 0)位置最近的节点，命名为MID_NODE
*GET, DISP, NODE, MID_NODE, U, Y	! 提取名称为MID_NODE的节点的竖向位移UY，存入变量DISP
MID_ELM = ENEARN (MID_NODE)	! 选择距离MID_NODE节点最近的单元, 命名为MID_ELM
ETABLE, STRS, LS, 3	! 以LS, 3（-Y一侧的弯曲应力SBYB）为内容定义单元表STRS
*GET, STRSS, ELEM, MID_ELM, ETAB, STRS	! 提取单元表STRS中单元MID_ELM的应力，存入STRSS
*STATUS, PARM
FINISH