! 8.4 悬臂宽板的模态分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH ! 退出以前模块
/CLEAR, NOSTART ! 清除内存中的所有数据, 不读入初始化文件
! （1）定义工程选项、单元类型、材料参数和实常数
/FILNAME, EX8.4 ! 定义文件名称EX8.4
/UNITS, SI ! 使用国际制单位
/TITLE, EX8.4, Vibration Analysis of Cantilever Plate.
/PREP7 ! 进入前处理模块PREP7
ET, 1, SHELL63 ! 定义第1类单元为壳单元SHELL63
MP, EX, 1, 2.07e11 ! 定义第1类材料的弹性模量EX=207 GPa
MP, DENS, 1, 7800 ! 定义第1类材料的密度DENS=7800 kg/m^3
MP, NUXY, 1, 0.33 ! 定义第1类材料的泊松比NUXY=0.33
R, 1, 0.03 ! 定义单元的第1类实常数:厚度为30mm
! （2）定义节点位置
N, 1, 0, -0.06 ! 在（0, -60mm）定义1号节点
N, 11, 0.5, -0.06 ! 在（500mm, -60mm）定义11号节点
FILL, 1, 11 ! 填充生成1到11号节点中间的所有节点
N, 21, 0, 0 ! 在（0, 0）定义21号节点
N, 31, 0.5, 0 ! 在（500mm, 0）定义31号节点
FILL, 21, 31 ! 填充生成21到31号节点中间的所有节点
N, 41, 0, 0.06 ! 在（0, 60mm）定义41号节点
N, 51, 0.5, 0.06 ! 在（500mm, 60mm）定义51号节点
FILL, 41, 51 ! 填充生成41到51号节点中间的所有节点
! （3）定义单元
E, 1, 2, 22, 21 ! 用节点1, 2, 22, 21定义第1个单元
EGEN, 10, 1, 1 ! 按第1个单元的编号规则生成10次, 得到1到10单元
EGEN, 2, 20, 1, 10 ! 按第1单元到第10个单元的编号模式生成2次, 增量为20, 得到11到20单元
FINISH ! 退出前处理模块/PREP7
! （4）进入求解模块, 指定求解选项和位移约束后求解
/SOLU ! 进入求解模块SOLUTION
ANTYPE, MODAL ! 申明求解类型是模态分析
MODOPT, LANB, 5, 0, 0, , OFF, , 2 ! 使用分块Lanczos方法求解前5阶振型和频率
D, 1, ALL, 0, 0, 11, 1 ! 固定1号节点, 用增量1循环到11, 生成1到11号共11个固定点
SOLVE ! 开始模态分析的求解
FINISH ! 退出后模块
! （5）在后处理模块中, 图形显示和动画显示模态
/POST1 ! 进入后处理模块POST1
SET, LIST ! 列出各阶频率计算结果
/VIEW, 1, 1, 1, 1 ! 切换到等轴测（1, 1, 1）方向
SET, FIRST ! 读入第1阶模态数据
PLDISP, 1 ! 在图形窗口显示结构振型
SET, Next ! 读入下一阶振型数据, 使用多次可以依次得到各阶振型和频率
PLDISP, 1 ! 在图形窗口显示结构振型
ANMODE, 10, 0.05 ! 用10帧每隔0.05秒钟的动画显示振型
FINISH
