! 8.6 用实体单元分析圆环振动
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
FINISH
/CLEAR, NOSTART ! 清除数据, 不加载初始化文件
! （1）定义工程选项、分析类型、单元类型、实常数、模态求解选项和材料参数
/FILNAME, EX8.6 ! 定义工程文件名称为EX8.6
/PREP7 ! 加载前处理模块
/TITLE, EX8.6, RADIAL VIBRATIONS OF A CIRCULAR RING FROM AN AXISYMMETRIC MODEL.
ANTYPE, MODAL ! 模态分析
ET, 1, PLANE25 ! 轴对称协调固体单元PLANE25
MXPAND, 1 ! 展开第1阶模态
MP, EX, 1, 207E9 ! 定义第1类材料的弹性模量为EX = 207 GPa
MP, DENS, 1, 7802.2 ! 定义第1类材料的密度DENS = 7802.2 kg/m^3
MP, PRXY, 1, 0 ! 定义泊松比等于PRXY = 0
! （2）定义节点和单元
LOCAL, 11, , 0.254 ! 在圆环内侧254mm处建立局部坐标系
N, 1 ! 在坐标原点（0, 0, 0）定义第1号节点
N, 2, , 0.00127 ! 在坐标（0, 1.27mm, 0）定义第2号节点
! 节点生成命令格式：NGEN, ITIME, INC, NODE1, NODE2, NINC, DX, DY, DZ, SPACE
NGEN, 2, 2, 1, 2, 1, 0.00127 ! 循环2次, 节点号加2, 按照节点1到2的位置, 将X坐标增加0.00127m。
E, 1, 3, 4, 2 ! 过节点1, 3, 4和2定义单元
CP, 1, UX, 1, 2 ! 定义耦合自由度集合1为从1号节点到2号节点的UX自由度
M, 1, UX ! 定义1号节点的UX为主自由度
FINISH ! 退出前处理模块PREP7
! （3）在求解模块中定义位移约束和模态求解选项后, 开始求解
/SOLU ! 加载求解模块
OUTPR, ALL, 1 ! 输出第1阶模态的所有计算结果
D, ALL, UZ, 0 ! 位移约束所有节点的UZ自由度
D, 1, UY, 0 ! 位移约束1号节点的UY自由度
MODE, 0, 1 ! 计算0阶轴对称模态
MODOPT, REDUC ! 使用降阶模态求解器
SOLVE ! 开始进行模态分析的求解
! （4）提取计算结果
*GET, FREQ0, MODE, 1, FREQ ! 提取第1阶频率, 存入变量FREQ0
FINISH ! 退出求解模块SOLUTION
! （5）重新进入求解模块, 改变位移约束后, 重新进行模态分析
/SOLU ! 加载求解模块
OUTPR, ALL, 1 ! 输出第1阶模态的所有计算结果
DDELE, ALL ! 删除所有节点的所有位移约束
M, 1, UX, , , UZ ! 定义1号节点的UX和UZ为主自由度
D, 1, UY ! 位移约束1号节点的UY自由度
MODE, 2, 1 ! 轴对称荷载的第2阶模态
MXPAND, 1, 0, 100 ! 展开频率范围从0到100的第1阶频率
SOLVE ! 开始求解
! （6）提取计算结果
*GET, FREQ2, MODE, 1, FREQ ! 提取第1阶频率, 存入变量FREQ2
*STATUS, PARM ! 显示内存变量的值
FINISH ! 退出求解模块SOLU
