ex7.6.inp

ansys有限元分析程序

! 7.6  简支梁的自振频率计算
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH	! 结束前面的模块
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! （1）定义工程选项、单元类型、材料参数和实常数
/FILNAME, EX7.6(1)	! 定义工程文件名称
/PREP7	! 加载前处理模块
/TITLE, EX7.6(1), FUNDAMENTAL FREQUENCY OF A SIMPLY SUPPORTED BEAM.
ET, 1, BEAM3	! 定义第1类单元为二维弹性梁单元BEAM3
MP, EX, 1, 207E9	! 定义第1类材料的弹性模量EX＝207 GPa
MP, DENS, 1, 7781	! 定义第1类材料的密度DENS＝7781 kg/m^3
MP, GXY, 1, 207E9/2.6	! 定义第1类材料的剪切弹性模量GXY＝207E9/2.6
R, 1, 25.81e-4, 55.493e-8, 50.8e-3	! 定义二维梁的截面积25.81e-4 m^2, 抗弯惯性矩55.493e-8 m^4和截面高度50.8e-3 m.
! （2）通过关键点的定义建立模型
K, 1	! 在坐标原点(0, 0, 0)定义第1个关键点
K, 2, 2.032	! 在最右端支座位置(2.032m, 0)定义第2个关键点
L, 1, 2	! 过关键点1, 2定义直线
LESIZE, ALL, , , 4	! 定义所有直线在剖分单元时, 剖分为4份
LMESH, 1	! 用梁单元BEAM3剖分第1条线段
/PNUM, NODE, 1	! 打开节点编号显示开关
EPLOT 	! 图形显示单元布置情况
FINISH	! 退出前处理模块PREP7
! （3）在求解模块定义模态分析选项, 施加位移约束后开始模态分析的求解
/SOLU	! 加载求解模块
ANTYPE, MODAL	! 指定分析类型为模态分析
MODOPT, REDUC, 3, , , 3	! 指定用降阶方法求解前3阶模态
MXPAND, 3	! 展开到第3阶模态
M, 3, UY, 5	! 将第3个节点到第5个节点的UY自由度设置为主自由度
OUTPR, ALL, 1	! 输出第1阶频率的所有内容
DK, ALL, UX	! 约束所有关键点在X方向的自由度UX
DK, ALL, UY	! 约束所有关键点在Y方向的自由度UY
! （4）在求解模块, 开始模态分析的求解
PSOLVE, ELFORM	! 开始生成单元刚度矩阵
PSOLVE, TRIANG	! 三角化总体刚度矩阵
PSOLVE, EIGREDUC	! 用HOUSEHOLDER方法计算特征值和特征向量 
PSOLVE, EIGEXP	! 展开特征向量的解答
FINISH	! 退出求解模块SOLUTION
! （5）在时间历程后处理器获取频率计算结果
/POST26	! 进入时间历程后处理器Post26
*GET, FREQ1, MODE, 1, FREQ	! 获取第1阶频率计算结果, 存入变量FREQ1
*GET, FREQ2, MODE, 2, FREQ	! 获取第2阶频率计算结果, 存入变量FREQ2
*GET, FREQ3, MODE, 3, FREQ	! 获取第3阶频率计算结果, 存入变量FREQ3
*STATUS, PARM  	! 显示内存变量的值
! （6）进入一般后处理模块/POST1, 图形显示振型
/POST1	! 进入一般后处理模块/POST1
set, 1, 1	! 读入第1阶模态的计算结果
PLDISP, 2	! 图形显示第1阶振型
set, 1, 2	! 读入第2阶模态的计算结果
PLDISP, 2	! 图形显示第2阶振型
set, 1, 3	! 读入第3阶模态的计算结果
PLDISP, 2	! 图形显示第3阶振型
FINISH	! 退出后处理模块