ex7.8.inp

ansys有限元分析程序

! 7.8  有轴向压力作用的简支梁的自由振动
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! （1）定义工程选项、分析类型、单元类型、实常数和材料参数
/FILNAME, EX7.8	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX7.8, LATERAL VIBRATION OF AN AXIALLY LOADED BAR.
ANTYPE, STATIC	! 定义分析类型为静力分析
PSTRES, ON	! 打开预应力开关(对于含有预应力的模态分析, 必须打开该选项)
ET, 1, BEAM4	! 定义第1类单元为空间梁单元BEAM4
R, 1, 25.81E-4, 55.493E-8, 55.493E-8, 50.8E-3, 50.8E-3	! 定义梁单元的截面几何参数：截面积25.81E-4 m^2, 两个抗弯惯性矩55.493E-8 m^4和两个截面高度50.8E-3 m.
MP, EX, 1, 207E9	! 定义第1类材料的弹性模量EX = 207GPa
MP, DENS, 1, 7781	! 定义第1类材料的密度DENS = 7781kg/m^3
MP, PRXY, , 0.3	! 定义第1类材料的泊松比PRXY = 0.3
! （2）定义节点位置和单元连接
N, 1	! 在坐标原点(0, 0, 0)定义第1号结点
N, 14, 2.032	! 在最右端(2.032m, 0, 0)处定义节点14
FILL	! 在1号节点和14号节点中间, 填充生成其余的节点
E, 1, 2	! 定义过1, 2节点的第1个单元
EGEN, 13, 1, 1	! 按照前一个单元的样式, 循环13次生成其余的12个单元
! （3）施加位移约束和荷载
D, 1, UY, , , 14, , ROTX, ROTZ	! 约束节点1到14的UY自由度和绕X轴Z轴的转动自由度
D, 1, UX, , , , , UZ	! 约束1号节点的UX和UZ自由度, 1号节点只剩下绕Y轴的转动自由度
D, 14, UZ	! 约束14号节点的UZ自由度
F, 14, FX, -177932	! 在14号节点上施加沿着-X方向大小等于177932N的集中力
M, 2, UZ, 13	! 选择2号到13号节点的UZ为主自由度
FINISH	! 退出前处理模块PREP7
! （4）进入求解模块, 指定求解选项并求解
/SOLU	! 进入求解模块
OUTPR, BASIC, 1	! 指定输出选项为第1阶频率的基本输出参数
SOLVE	! 开始静力分析的求解
FINISH	! 退出求解模块SOLUTION
! （5）进入一般后处理模块, 提取频率计算结果
/POST1	! 进入一般后处理模块
ETABLE, STRS, LS, 1	! 以杆件轴向应力为内容, 定义单元表STRS
*GET, STRSS, ELEM, 13, ETAB, STRS	! 从单元表STRS中提取13号单元的应力, 存入STRSS
*GET, DEF, NODE, 14, U, X	! 提取14号节点的水平位移UX, 存入变量DEF
! （6）重新进入求解模块, 用模态求解方法求解
/SOLU  	! 加载求解模块
ANTYPE, MODAL	! 指定求解类型为模态分析
PSTRES, ON	! 打开预应力开关
MXPAND, 3 	! 指定展开前3阶模态
MODOPT, LANB, 3	! 选择分块LANCZOS特征值求解算法
SOLVE 	! 开始模态分析的求解
! （7）提取频率计算结果
*GET, FREQ1, MODE, 1, FREQ	! 获取第1阶频率, 存入变量FREQ1
*GET, FREQ2, MODE, 2, FREQ	! 获取第2阶频率, 存入变量FREQ2
*GET, FREQ3, MODE, 3, FREQ	! 获取第3阶频率, 存入变量FREQ3
*STATUS, PARM  	! 显示内存变量的值
FINISH	! 退出后处理模块