/prep7
!注意：簡支板的模型必須用殼單元，不能用實體建模，它是無厚度的，建成一塊有厚度的板是錯誤的，會產生很大的偏差
!目前四個角上的節點只能通過編號選擇，不過如果單元等分數發生變化，則節點號會隨著變化；所以改變等分數時，程序相應部分要改變。
lengthx=1.5     ! x方向長度
lengthy=0.3     !y方向長度
lengthz=0.005
div=30          !單元等分數
nmodes=15       !分析的模態階數
!
ET,1,SHELL63
R,1,lengthz, , , , , ,
!
DENS,1,2770         ! (Rho) density  
EX,1,71e9
NUXY,1,0.33
!
K,1,0,0,0,  
K,2,lengthx,lengthy,0,  
BLC4,0,0,lengthx,lengthy
!
LESIZE,ALL, , ,div,1,1   
CM,_Y,AREA                !以下8行皆為對面積的網格化，批處理，無需理解。
ASEL, , , ,       1 
CM,_Y1,AREA 
CHKMSH,'AREA'   
CMSEL,S,_Y    
AMESH,_Y1   
CMDEL,_Y
CMDEL,_Y1   
CMDEL,_Y2  
!
LSEL,S, , ,       2             !area 4 faces to area 6
LSEL,A, , ,       4 
NSLL,S,1
D,ALL,UX,,,,,UY,UZ,rotx,rotz    ! Constrain all displacements to zero at the walls  
!
LSEL,S, , ,       1             !area 5 faces to area 7
LSEL,A, , ,       3 
NSLL,S,1
D,ALL,UX,,,,,UY,UZ,roty,rotz    ! Constrain all displacements to zero at the walls  
!
LSEL,S, , ,       1             ! select crossecting lines and get the points of the corners 
NSLL,S,1
LSEL,S, , ,       2
LSEL,A, , ,       4 
NSLL,R,1
!DDELE,all,all,,
DDELE,all,rotx
DDELE,all,roty
!
LSEL,S, , ,       3 
NSLL,S,1
LSEL,S, , ,       2
LSEL,A, , ,       4 
NSLL,R,1
!DDELE,all,all,,
DDELE,all,rotx
DDELE,all,roty
!
nsel,all
esel,all
!
FINISH                ! begin to solve  
/SOLU   
!*  
ANTYPE,2
!*  
MODOPT,LANB,nmodes    !nmodes為賦值變量，模態階數
EQSLV,SPAR  
MXPAND,nmodes, , ,0 
LUMPM,0 
PSTRES,0
!*  
MODOPT,LANB,nmodes,0,100000000, ,OFF, ,0
/STAT,SOLU  
SOLVE