ex7.9.inp

ansys有限元分析程序

! 7.9  用壳体单元计算悬臂等强度梁的自由振动
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除数据, 不加载初始化文件
! （1）定义工程选项、分析类型、单元类型、实常数和材料参数
/FILNAME, EX7.9	! 定义工程文件名称
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化网格功能
/TITLE, EX7.9, VIBRATION OF A WEDGE (SHELL63 WITH BENDING STIFFNESS).
ANTYPE, MODAL  	! 指定分析类型为模态分析
MXPAND, 1     	! 只展开第1阶模态
MODOPT, REDUC, 2, , , 2	! 用降阶方法求解前2阶振型
ET, 1, SHELL63, 2	! 使用只考虑弯曲效应的厚壳单元SHELL63
R, 1, 25.4E-3	! 定义壳体厚度为25.4mm
MP, EX, 1, 207E9	! 定义第1类材料的弹性模量EX＝207GPa
MP, DENS, 1, 7781	! 定义第1类材料的密度DENS＝7781kg/m^3
MP, NUXY, 1, 0	! 定义第1类材料的泊松比NUXY＝0
! （2）通过关键点定义线, 通过关键点定义面, 对面进行面单元剖分得到有限元模型
K, 1	! 在坐标原点（0, 0, 0）定义第1个关键点
K, 2, 406.4E-3	! 在坐标(406.4mm, 0)定义第2个关键点
K, 3, 406.4E-3, 50.8E-3	! 在坐标(406.4mm, 50.8mm)定义第3个关键点
L, 2, 3	! 过2, 3关键点定义直线
LESIZE, 1, , , 1	! 指定定义的1号线, 单元剖分段数为1
A, 1, 2, 3, 3	! 过关键点1, 2, 3, 3定义1号面
ARSYM, Y, 1	! 关于Y坐标用1号面, 对称生成2号面
NUMMRG, KP	! 合并重合的关键点
ESIZE, , 4	! 指定单元剖分段数为4段
AMESH, 1, 2	! 对定义好的1号面到2号面进行面单元剖分
! （3）选择某些位置的节点, 施加边界条件, 并指定输出选项
/PNUM, NODE, 1	! 打开节点显示
NSEL, S, LOC, Y, 0	! 选择所有在Y=0位置上的节点
M, ALL, UZ	! 将所有选择集中的节点的UZ自由度设定为主自由度
NSEL, S, LOC, X, 406.4E-3	! 选择所有在X=406 mm位置上的节点
D, ALL, UZ, , , , , ROTX, ROTY, UX, UY	! 约束选择集中节点的UZ, ROTX, ROTY, UX, UY自由度
NSEL, ALL	! 选择所有结点
OUTPR, ALL, 1	! 指定输出选项为第1阶频率的所有计算输出参数
OUTRES, ALL, 0	! 指定所有输出参数不用保存到数据库
FINISH	! 退出前处理模块PREP7
! （4）进入求解模块, 开始求解并提取频率计算结果
/SOLU	! 加载求解模块
SOLVE	! 开始模态分析的求解
*GET, FREQ1, MODE, 1, FREQ	! 获取第1阶频率, 存入变量FREQ1
FINISH	! 退出求解模块SOLU
! （5）重新进入前处理模块, 定义使用带壳体效应的壳单元
/PREP7	! 加载前处理模块
SMRT, OFF	! 关闭智能化功能
C***      USING SHELL63 ELEMENTS WITH SHELL OPTION
ET, 1, SHELL63	! 使用弯曲效应和薄膜效应的单元SHELL63
FINISH 	! 退出前处理模块PREP7
! （6）用壳体单元重新计算
/SOLU	! 加载求解模块
SOLVE	! 开始模态分析的求解
*GET, FREQ2, MODE, 1, FREQ	! 获取第1阶频率, 存入变量FREQ2
*STATUS, PARM  	! 显示内存变量的值
FINISH	! 退出后处理模块