/filename,ex6-4                      	!趨
/title,The Dynamic Analysis of the Dam  !趨
/prep7                             	!ǰģ
*set,h,150               		!-Ӹ߶
Et,1,plane42             		!õԪ
!ϲ
Mp,ex,1,2.9e10          	        !趨135m»ϲ
Mp,prxy,1,0.17
Mp,dens,1,2500
Mp,ex,2,2.5e10          		!趨135mϻϲ
Mp,prxy,2,0.17
Mp,dens,2,2500
Mp,ex,3,2.9e10         		     	!趨ʯϲ
Mp,prxy,3,0.3
Mp,dens,3,2600
!ؼ
K,1,
K,2,0.9*0.7*h+0.1*h
K,3,0.1*h,0.9*h
K,4,,0.9*h
!
L,1,2
L,2,3
L,3,4
L,4,1
!
Al,1,2,3,4
Rectng,0,0.1*h,0.9*h,h
!ദؼ㼰
K,9,,0.2*h
K,10,-0.1*h
A,1,9,10
!漰ؼ
K,11,-1.5*h
K,12,-1.5*h,-2*h
K,13,0.9*0.7*h+2.1*h,-2*h
K,14,0.9*0.7*h+2.1*h
A,11,12,13,14
K,15,-0.1*h,-2*h
K,16,,-2*h
K,17,0.9*0.7*h+0.1*h,-2*h
A,1,10,15,16
A,16,17,2,1
Aplot
aovlap,all
aglue,all
/pnum,area,1
/pnum,line,1
NUMMRG,ALL,all                   		!ϲظԪ
NUMCMP,ALL                      		!ѹб
APLOT 
!񻮷
Lsel,s,,,15,16,1               !4񻮷֣ӦڰĶ
Lesize,all,,,10                !趨߷ָΪ10
lsel,s,,,3,4,1
lesize,all,,,10
mat,2
amesh,4
!
lsel,s,,,1                     !7񻮷֣Ӧڰ²
lesize,all,,,10
lsel,s,,,2
lesize,all,,,30
allsel
lsel,s,,,5,21,16
lccat, all                !521Ӳӳ񻮷
Mat,1                     !Ϻ1
MSHAPE,0,2D               !趨״Ϊ2Dı
MSHKEY, 1                 !趨񻮷ַʽΪӳ񻮷
amesh,7
ldele,22                  !ɾ22
amesh,1                   !1񻮷֣Ӧڰ
allsel
lsel,s,,,8,14,2           !2񻮷֣Ӧڰ²Ļ
lsel,a,,,9
lesize,all,,,20
lsel,s,,,17,18,1
lesize,all,,,15
lsel,s,,,19,20,1
lesize,all,,,20
Mat,3
MSHAPE,0,2D 
MSHKEY, 1
amesh,2
amesh,3                  !3񻮷֣Ӧڰ²Ļ
amesh,5                  !5񻮷֣Ӧڰ²Ļ
amesh,6                  !6񻮷֣Ӧڰ²Ļ
save
!ģ
/solu
Allsel
Lsel,s,,,8,9,1                      	!ѡı߽
Dl,all,,ux,0                       	!ʩXԼ
Allsel
Lsel,s,,,18,20,2                    	!ѡױͼԪ
Lsel,a,,,11,13,2                    	!ѡױͼԪ
Dl,all,,all                         	!ԵױʩȫԼ
Acel,,10                          	!ʩٶ
!ģ̬ʼӿռ䷨ģ̬ȡǰ10׵ģ̬
antype,modal                      	!趨ģ̬
modopt,subsp,10                   	!趨ģ̬ѡ 
solve
save,dam_dynamic_rst1,db
finish
!Ӧ׷ʼ
/solu         
Antype,spectr              	!趨Ϊ׷
Spopt,sprs,10,yes           	!ָΪӦ
Svtype,2                  	!趨Ӧ
Sed,1,1,                   	!趨
Freq,0.1603,0.1671,0.1864,0.1932,0.2191,0.2914,0.3223,0.3301,0.3669,0.4781
Sv,2,2,2,2,2,2,2,2,2,0.9128    	!趨Ƶʼֵ
Solve
Save,rst2,db
Finish
!ģ̬չ
/solu     
Antype,modal              	!趨ģ̬
Expass,on                 	!ģ̬չѡ
Mxpand,10,,,yes,0.005       	!ģ̬չ
Solve
Save,rst3,db
Finish
!ϲģ̬
/solu    
Antype,spectr              	!趨Ϊ׷
Srss,0.02,disp              	!趨ģ̬ϲѡ
Solve
Save,rst4,db
Finish
!ͨú鿴
/post1
/input,,mcom              		!ϲģ̬
Esel,s,mat,,1
Esel,a,mat,,2
Eplot
Set,first                  		!һغɲ
Plnsol,u,sum,1,1           		!λƷֵֵͼ
Set,next                  		!һغɲ
Plnsol,u,sum,1,1
Set,next
Plnsol,u,sum,1,1
Set,first
Plnsol,s,1,0,1
Plnsol,epto,1,0,1
Set,next
Plnsol,s,1,0,1
!ʱ̺鿴
/post26
NSOL,2,22,U,X, UX_2     		!趨ڵ22xλΪ2
NSOL,3,22,U,Y, UY_3     		!趨ڵ22yλΪ3
XVAR,1                 		        !趨xΪʱ
PLVAR,2,3,              		!ͼʾ23 
FINISH                  		!ͨúģ
/EXIT,ALL               		!ȫݲ˳
