ex3.5(1).inp

ansys有限元分析程序

! 3.5  矩形截面梁的纵横弯曲分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据，不加载初始化文件
/FILNAME, EX3.5(1)	! 定义文件名称
/PREP7	! 加载前处理模块
/TITLE, EX3.5(1), TIE ROD WITH LATERAL LOADING, NO STRESS STIFFENING.
! （1）定义求解类型，实常数和材料参数
ANTYPE, STATIC	! 定义分析类型为静力分析
ET, 1, BEAM4, , , , , , 1	! 定义第1类单元为三维弹性梁单元BEAM4
R, 1, 0.403225e-2, 0.135492e-5, 0.135492e-5, 0.0635, 0.0635	! 定义第1类实常数、 梁的截面积0.403225e-2m^2, 、惯性矩IZZ=IYY=0.135492e-5 m^4、截面高度B=H=6.35cm
MP, EX, 1, 210E9	! 定义第1类材料的弹性模量EX＝210GPa
MP, PRXY, , 0.3	! 定义材料的泊松比为PRXY＝0.3
! （2）定义节点位置
N, 1	! 在(0, 0)定义第1个节点
N, 5, 2.54	! 在(2.54, 0)定义第5个节点
FILL	! 利用节点填充FILL命令生成中间3个节点
! （3）定义单元连接关系
E, 1, 2	! 过1, 2节点定义第1个单元
EGEN, 4, 1, 1	! 按照前面的单元模式，生成4次，每次两端节点号加1，生成其余3个单元
! （4）定义位移约束和荷载条件
D, ALL, UY, , , , , ROTX, ROTZ	! 约束所有节点的UY、ROTX、ROTZ自由度
D, 1, UZ	! 约束1号节点的UZ自由度
NSEL, S, , , 5	! 选择编号不小于5的所有节点
DSYM, SYMM, X	! 定义所有节点(5号节点)关于x轴的对称位移边界条件
NSEL, ALL	! 选择所有节点
! （5）施加荷载
F, 1, FX, -97740	! 在1号节点施加沿着X方向的集中力-97740N
SFBEAM, ALL, 1, PRES, 314	! 在梁表面施加均布荷载314N/m
FINISH	! 退出前处理模块PREP7
! （6）进入求解模块，开始求解
/SOLU	! 加载求解模块
SOLVE	! 开始静力分析求解
FINISH	! 退出求解模块SOLU
! （7）进入后处理模块，显示位移计算结果，提取左边支座节点的转角和跨中挠度
/POST1	! 加载后处理模块
NSEL, S, , , 1, 5, 4	! 从1到5按照增量4循环，选中1号节点和5号节点
PRNSOL, U, Z	! 列表显示选择集中所有节点的Z方向的位移UZ
PRNSOL, ROT, Y	! 列表显示选择集中所有节点的绕Y轴的转角ROTY
NSEL, ALL	! 选择所有节点
PRRSOL	! 列表显示被约束节点的支反力（Print Reaction SOLution）
RGHT_END = NODE (2.54, 0, 0)	! 选择距离(2.54, 0, 0)最近的节点(只有5号节点), 命名为RGHT_END
LFT_END  = NODE (0, 0, 0)	! 选择距离(0, 0, 0)最近的节点，命名为LFT_END
*GET, UZ_MX_C2, NODE, RGHT_END, U, Z	! 获取节点RGHT_END的线位移UZ
*GET, SLOPE_C2, NODE, LFT_END, ROT, Y	! 获取节点LFT_END的转角ROTY
FINISH	! 退出后处理模块Post1
! （8）在时间历程后处理器中，叠加处理5号节点位置的弯矩
/POST26	! 进入时间历程后处理器Post26
RFORCE, 2, RGHT_END, M, Y	! 叠加RGHT_END节点的弯矩MY，存入第2个Post26变量
STORE	! 保存变量
*GET, M_MX_C2, VARI, 2, EXTREM, VMAX	! 获取第2个Post26变量的最大值
FINISH	! 退出后处理模块
! （9）重新进入前处理模块，重新计算考虑应力强化效应的压杆
/PREP7	! 加载前处理模块
/TITLE, EX3.5(2), TIE ROD WITH LATERAL LOADING, STRESS STIFFENING PRESENT.
SSTIF, ON	! 打开应力强化选项
NSUBST, 5	! 指定这次分析所要使用的荷载子步为5步
AUTOTS, ON	! 激活自动时间步长
FINISH	! 退出前处理模块PREP7
! (10) 进入求解模块，指定最小收敛误差后，开始非线性求解
/SOLU	! 加载求解模块
CNVTOL, F, , .0001, , 1	! 最小收敛误差
SOLVE	! 开始求解
FINISH	! 退出求解模块SOLU
! (11) 进入后处理模块，显示位移计算结果，提取左边支座节点的转角和跨中挠度
/POST1	! 进入一般后处理模块
NSEL, S, , , 1, 5, 4	! 从1到5按增量4循环，选中1号节点和5号节点
PRNSOL, U, Z	! 列表显示选择集中所有节点（1号和5号节点）的Z方向的位移UZ
PRNSOL, ROT, Y	! 列表显示选择集中所有节点的绕Y轴的转角ROTY
PRRSOL	! 列表显示被约束节点的支反力
*GET, UZ_MX_C1, NODE, RGHT_END, U, Z	! 获取节点RGHT_END的线位移UZ
*GET, SLOPE_C1, NODE, LFT_END, ROT, Y	! 获取节点LFT_END的转角ROTY
FINISH	! 退出求解模块SOLU
! (12) 在时间历程后处理器中，叠加处理跨中（5号节点）位置的弯矩
/POST26	! 进入时间历程后处理器Post26
RFORCE, 2, RGHT_END, M, Y	! 叠加RGHT_END节点的弯矩MY，存入第2个Post26变量
STORE	! 保存变量
*GET, M_MX_C1, VARI, 2, EXTREM, VMAX	! 获取第2个Post26变量的最大值
*STATUS, PARM
FINISH