! 8.5  固支圓板的自由振動
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數據, 不加載初始化文件
! （1）設置工程選項、單元類型、材料參數
/FILNAME, EX8.5	! 定義文件名稱
/PREP7  	! 加載前處理模塊
SMRT, OFF  	! 關閉智能化網格功能
/TITLE, EX8.5, NATURAL FREQUENCY OF A FLAT CIRCULAR PLATE WITH CLAMPED EDGE.
ET, 1, PLANE2, , , 1	! 使用軸對稱平面單元PLANE2
MP, EX, 1, 207E9	! 定義第1類材料的彈性模量EX ＝ 207 GPa
MP, PRXY, , 0.3	! 定義材料的泊松比為PRXY ＝ 0.3
MP, DENS, 1, 7801.6	! 定義第1類材料的密度DENS = 7801.6 kg/m^3
! （2）定義關鍵點、線和面, 并將面剖分為面單元
K, 1    	! 在坐標原點(0, 0, 0)定義第1個關鍵點
K, 2, 0.4318 	! 在坐標(0.4318m, 0)定義第2個關鍵點
KGEN, 2, 1, 2, 1, , 0.0127	! 循環2次, 按照節點1到節點2的坐標, Y坐標增加12.7mm
L, 1, 2      	! 過關鍵點1, 2定義直線
L, 3, 4     	! 過關鍵點3, 4定義直線
LESIZE, ALL, , , 10	! 所有直線在剖分單元時, 剖分為10段
ESIZE, , 1 	! 單元尺寸為剖分為1段
A, 1, 2, 4, 3	! 過節點1, 2, 3和4定義面
AMESH, 1	! 對1號面進行面單元剖分
! （3）選擇邊界節點, 施加位移約束
NSEL, R, LOC, X, 0	! 重新選擇位于X＝0的節點
D, ALL, UX	! 對所有節點的X方向位移施加位移約束UX＝0
NSEL, ALL	! 選擇所有節點
DK, 2, ALL, , , 1	! 約束2號關鍵點的所有自由度, 并將其標識為約束展開
DK, 4, UX, , , 1	! 約束4號關鍵點的UX自由度, 并將其標識為約束展開
FINISH	! 退出前處理模塊PREP7
! （4）指定求解類型和求解方法后, 開始模態分析的求解
/SOLU  	! 加載求解模塊
ANTYPE, MODAL	! 指定分析類型為模態分析
MODOPT, SUBSP, 9	! 指定使用子空間迭代法求解器, 展開前9解模態
EQSLVE, PCG	! 指定方程求解器為PCG（預條件的共軛梯度循環求解器）
SOLVE	! 開始模態分析的求解
! （5）提取頻率計算結果
*GET, F1, MODE, 1, FREQ	! 提取第1階頻率, 存入變量F1
*GET, F2, MODE, 2, FREQ	! 提取第2階頻率, 存入變量F2
*GET, F3, MODE, 3, FREQ	! 提取第3階頻率, 存入變量F3
*STATUS, PARM	! 查看保存到變量中的計算結果
FINISH	! 退出求解模塊SOLU