! 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