!9.4  懸臂長板的自由振動頻率
! 本程序來源于邢靜忠等編著的《ANSYS分析實例與工程應用》，機械工業出版社（2004年）
FINISH
/CLEAR, NOSTART
! (1) 定義工程選項,分析類型,單元類型,實常數,模態求解選項和材料參數
/FILENAME, EX9.4
/PREP7
SMRT,OFF
/TITLE, EX9.4 , VIBRATION OF A FLAT PLATE
C***   VIBRATION PROBS. IN ENGR., TIMOSHENKO, 3RD. ED., PAGE 338, ART. 53
C***     USING SHELL63 ELEMENTS WITH SHELL OPTION
ANTYPE,MODAL
ET,1,SHELL63
R,1,25.4e-3
MXPAND,1
MODOPT,REDUC,2,,,2
MP,EX,1,207e9
MP,DENS,1,7780.81
MP,NUXY,1,0
! (2) 定義關鍵點,線,通過對線的單元劃分得到有限元模型
K,1,,-50.8e-3
K,2,406.4e-3,-50.8e-3
KGEN,2,1,2,1,,101.6e-3
L,1,3
L,2,4
LESIZE,ALL,,,2
A,1,2,4,3
ESIZE,,4
AMESH,1
! (3) 選擇邊界上的節點,施加位移邊界條件
/PNUM, NODE, 1
NSEL,S,LOC,Y,0
M,ALL,UZ
NSEL,S,LOC,X,406.4e-3
D,ALL,ALL
NSEL,ALL
OUTPR,ALL,1
/REPLOT
FINISH
! (4) 在求解模塊求解模態分析
/SOLU
SOLVE
! (5) 提取計算結果
*GET,FREQ,MODE,1,FREQ
! (6) 定義數組變量,將計算結果及其與理論解的比值輸出到文本文件
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = '     f, '
LABEL(1,2) = 'Hz      '
*VFILL,VALUE(1,1),DATA,128.12
*VFILL,VALUE(1,2),DATA,FREQ
*VFILL,VALUE(1,3),DATA,ABS(FREQ/128.12)
/COM
/OUT,  EX9_4 , out
/COM,-------------------   EX9.4  RESULTS COMPARISON --------------
/COM,
/COM,                 |   TARGET   |   ANSYS   |   RATIO
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,'   ',F10.2,'  ',F10.2,'   ',1F5.3)
/COM,----------------------------------------------------------
/OUT
FINISH
*LIST,  EX9_4 , out