! 7.9  用殼體單元計(jì)算懸臂等強(qiáng)度梁的自由振動(dòng)
! 本程序來源于邢靜忠編著的《ANSYS應(yīng)用實(shí)例與分析》，科學(xué)出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數(shù)據(jù), 不加載初始化文件
! （1）定義工程選項(xiàng)、分析類型、單元類型、實(shí)常數(shù)和材料參數(shù)
/FILNAME, EX7.9	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關(guān)閉智能化網(wǎng)格功能
/TITLE, EX7.9, VIBRATION OF A WEDGE (SHELL63 WITH BENDING STIFFNESS).
ANTYPE, MODAL  	! 指定分析類型為模態(tài)分析
MXPAND, 1     	! 只展開第1階模態(tài)
MODOPT, REDUC, 2, , , 2	! 用降階方法求解前2階振型
ET, 1, SHELL63, 2	! 使用只考慮彎曲效應(yīng)的厚殼單元SHELL63
R, 1, 25.4E-3	! 定義殼體厚度為25.4mm
MP, EX, 1, 207E9	! 定義第1類材料的彈性模量EX＝207GPa
MP, DENS, 1, 7781	! 定義第1類材料的密度DENS＝7781kg/m^3
MP, NUXY, 1, 0	! 定義第1類材料的泊松比NUXY＝0
! （2）通過關(guān)鍵點(diǎn)定義線, 通過關(guān)鍵點(diǎn)定義面, 對(duì)面進(jìn)行面單元剖分得到有限元模型
K, 1	! 在坐標(biāo)原點(diǎn)（0, 0, 0）定義第1個(gè)關(guān)鍵點(diǎn)
K, 2, 406.4E-3	! 在坐標(biāo)(406.4mm, 0)定義第2個(gè)關(guān)鍵點(diǎn)
K, 3, 406.4E-3, 50.8E-3	! 在坐標(biāo)(406.4mm, 50.8mm)定義第3個(gè)關(guān)鍵點(diǎn)
L, 2, 3	! 過2, 3關(guān)鍵點(diǎn)定義直線
LESIZE, 1, , , 1	! 指定定義的1號(hào)線, 單元剖分段數(shù)為1
A, 1, 2, 3, 3	! 過關(guān)鍵點(diǎn)1, 2, 3, 3定義1號(hào)面
ARSYM, Y, 1	! 關(guān)于Y坐標(biāo)用1號(hào)面, 對(duì)稱生成2號(hào)面
NUMMRG, KP	! 合并重合的關(guān)鍵點(diǎn)
ESIZE, , 4	! 指定單元剖分段數(shù)為4段
AMESH, 1, 2	! 對(duì)定義好的1號(hào)面到2號(hào)面進(jìn)行面單元剖分
! （3）選擇某些位置的節(jié)點(diǎn), 施加邊界條件, 并指定輸出選項(xiàng)
/PNUM, NODE, 1	! 打開節(jié)點(diǎn)顯示
NSEL, S, LOC, Y, 0	! 選擇所有在Y=0位置上的節(jié)點(diǎn)
M, ALL, UZ	! 將所有選擇集中的節(jié)點(diǎn)的UZ自由度設(shè)定為主自由度
NSEL, S, LOC, X, 406.4E-3	! 選擇所有在X=406 mm位置上的節(jié)點(diǎn)
D, ALL, UZ, , , , , ROTX, ROTY, UX, UY	! 約束選擇集中節(jié)點(diǎn)的UZ, ROTX, ROTY, UX, UY自由度
NSEL, ALL	! 選擇所有結(jié)點(diǎn)
OUTPR, ALL, 1	! 指定輸出選項(xiàng)為第1階頻率的所有計(jì)算輸出參數(shù)
OUTRES, ALL, 0	! 指定所有輸出參數(shù)不用保存到數(shù)據(jù)庫
FINISH	! 退出前處理模塊PREP7
! （4）進(jìn)入求解模塊, 開始求解并提取頻率計(jì)算結(jié)果
/SOLU	! 加載求解模塊
SOLVE	! 開始模態(tài)分析的求解
*GET, FREQ1, MODE, 1, FREQ	! 獲取第1階頻率, 存入變量FREQ1
FINISH	! 退出求解模塊SOLU
! （5）重新進(jìn)入前處理模塊, 定義使用帶殼體效應(yīng)的殼單元
/PREP7	! 加載前處理模塊
SMRT, OFF	! 關(guān)閉智能化功能
C***      USING SHELL63 ELEMENTS WITH SHELL OPTION
ET, 1, SHELL63	! 使用彎曲效應(yīng)和薄膜效應(yīng)的單元SHELL63
FINISH 	! 退出前處理模塊PREP7
! （6）用殼體單元重新計(jì)算
/SOLU	! 加載求解模塊
SOLVE	! 開始模態(tài)分析的求解
*GET, FREQ2, MODE, 1, FREQ	! 獲取第1階頻率, 存入變量FREQ2
*STATUS, PARM  	! 顯示內(nèi)存變量的值
FINISH	! 退出后處理模塊