通過對懸臂梁進行模態分析及提取剛度矩陣及質量矩陣完整程序。

finish

/clear

/config,nres,20000

/prep7

ee=6.96e10

b=0.5

h=0.05

lcd=5

aa=b*h

iz=b*h*h*h/12

iy=h*b*b*b/12

et,1,beam4

r,1,aa,iz,iy,h,b

mp,ex,1,6.96e10

mp,dens,1,2730

mp,prxy,1,0.33

mp,alpx,1,1e-5

k,1

k,2,5

l,1,2

lesize,all,,,20

numoff,node,1

lmesh,all

!!!節點重新編號

n,22,5,0,0

nummrg,node,,,,high

finish

alls

/solu

dk,1,all

!模態分析

/SOL 

ANTYPE,2   

MODOPT,LANB,10  

EQSLV,SPAR  

MODOPT,LANB,20,0,99999999, ,OFF   

SOLVE

finish

/post1

set,list

finish

!!!!創建子結構part1

/filnam,part1  

/solu

antype,substr      !分析類型 子結構

seopt,part1,2       !子結構一

!創建part1

nsel,s,node,,1,8

esln,r,1,all

cm,part1,elem

!創建interface

nsel,r,node,,8

cm,interface,node