! 8.1  圓形張緊薄膜的自由振動
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數據，不加載初始化文件
! （1）定義工程選項、分析類型、單元類型、材料參數和實常數
/FILNAME, EX8.1	! 定義工程文件名稱為EX8.1
/PREP7	! 加載前處理模塊/PREP7
/TITLE, EX8.1, VIBRATION OF A STRETCHED CIRCULAR MEMBRANE.
ANTYPE, STATIC	! 定義分析類型為靜力分析
PSTRES, ON	! 打開預應力選項（該選項對于有預應力的結構振動非常重要）
ET, 1, SHELL51	! 定義第1類單元為軸對稱殼單元SHELL51
MP, EX, 1, 207E9	! 定義第1類材料的彈性模量EX=207GPa
MP, PRXY, 1, 0.3	! 定義第1類材料的泊松比PRXY＝0.3
MP, DENS, 1, 7805	! 定義第1類材料的密度DENS＝7805kg/m^3
R, 1, 0.254e-3	! 定義第1類實常數定義殼體厚度0.254mm
! （2）通過節點定義建立模型
N, 1  	! 在坐標原點（0, 0, 0）定義第1號節點
N, 10, 0.381	! 在坐標（0, 0.381mm, 0）定義第10號節點
FILL	! 在1號節點和10號節點中間填充生成其余8個節點
E, 1, 2	! 過節點1和節點2定義軸對稱殼單元
EGEN, 9, 1, 1	! 用前面定義的單元, 循環9次, 每次單元的節點號增加1, 生成其余8個單元
! （3）施加位移約束和荷載
D, 1, UX, , , , , ROTZ	! 約束1號節點的沿著X方向的線位移UX和繞Z軸的角位移ROTZ
D, 10, UY, , , , , UZ	! 約束10號節點的沿著Y方向的線位移UY和沿著Z軸的線位移UZ
F, 10, FX, 41924	! 在節點10的X方向施加集中力, 大小為41924N
FINISH	! 退出前處理模塊PREP7
! （4）進入求解模塊, 指定求解輸出選項后求解
/SOLU 	! 加載求解模塊
OUTPR, BASIC, 1	! 指定只輸出第1階模態的基本計算結果
SOLVE  	! 開始求解
FINISH 	! 退出求解模塊SOLU
! （5）在一般后處理模塊中, 定義單元表, 提取單元的應力值
/POST1 	! 進入一般后處理模塊
ETABLE, STRS, LS, 1	! 以殼應力最大拉應力”LS, 1”為內容, 定義單元表STRS
*GET, STRSS, ELEM, 9, ETAB, STRS	! 從單元表STRS提取單元9的最大拉應力, 存入STRSS
FINISH	! 退出一般后處理模塊
! （6）重新進入求解模塊, 開始模態分析
/SOLU 	! 加載求解模塊
ANTYPE, MODAL	! 指定分析類型是模態分析
PSTRES, ON	! 指定包含預應力效應
MXPAND, 3	! 展開前3階模態
MODOPT, REDUC, , , , 9	! 用降階方法求解前9階模態
M, 1, UY, 9	! 定義節1到節點9的UY方向為主自由度方向
SOLVE	! 開始進行模態分析的求解
! （7）提取頻率計算結果
*GET, FREQ1, MODE, 1, FREQ	! 提取第1階頻率計算結果, 存入變量FREQ1
*GET, FREQ2, MODE, 2, FREQ	! 提取第2階頻率計算結果, 存入變量FREQ2
*GET, FREQ3, MODE, 3, FREQ	! 提取第3階頻率計算結果, 存入變量FREQ3
*Status, parameters	! 列出內存變量的值
FINISH	! 退出后處理模塊