/filename,ex6-1                        	        !趨
/title, The Dynamic Analysis of the Bridge   	!趨
/prep7                               	        !ǰģ
!**********************************
!嵥ԪͼԪʵ
!**********************************
Et,1,beam4                            	!嵥Ԫ
Et,2,shell63
Et,3,link10
!Ԫʵ
*dim,BBS,,43                          	!ʼϲ漸 
*dim,HHS,,43
*dim,SSS,,43
*dim,IIYYS,,43
*dim,IIZZS,,43
!
*do,i,1,43                          	!ͨDOѭ鸳ֵ
  BBS(i)=1.5*(i-1)/42+3
  HHS(i)=(i-1)/21+3
  SSS(i)=BBS(i)*HHS(i)  
  IIzzS(i)=BBS(i)*HHS(i)**3/12
  IIyyS(i)=HHS(i)*BBS(i)**3/12
*enddo                            	!DOѭֵ
!
*dim,BS,,42                        	!ʼϲƽ
*dim,HS,,42
*dim,SS,,42
*dim,IYYS,,42
*dim,IZZS,,42
*do,i,1,42                          	!ͨDOѭ鸳ֵ
  hS(i)=(hhS(i)+hhS(i+1))/2
  bS(i)=(bbS(i)+bbS(i+1))/2
  SS(i)=(SSS(i)+SSS(i+1))/2
  IYYS(i)=(IIYYS(i)+IIYYS(i+1))/2
  IZZS(i)=(IIZZS(i)+IIZZS(i+1))/2
*enddo                            	!DOѭֵ
!
*dim,BBX,,80                      	!ʼ²漸
*dim,HHX,,80
*dim,SSX,,80
*dim,IIYYX,,80
*dim,IIZZX,,80
!
*do,i,1,80                         		!ͨDOѭ鸳ֵ
  BBX(i)=2.5
  HHX(i)=(i-1)*4.7164/79+5
  SSX(i)=BBX(i)*HHX(i)  
  IIzzX(i)=BBX(i)*HHX(i)**3/12
  IIyyX(i)=HHX(i)*BBX(i)**3/12
*enddo                           	!DOѭֵ
!
*dim,BX,,79                      	!ʼ²ƽ
*dim,HX,,79
*dim,SX,,79
*dim,IYYX,,79
*dim,IZZX,,79
*do,i,1,79                        	!ͨDOѭ鸳ֵ
  BX(i)=(BBX(i)+BBX(i+1))/2
  HX(i)=(HHX(i)+HHX(i+1))/2
  SX(i)=(SSX(i)+SSX(i+1))/2
  IYYX(i)=(IIYYX(i)+IIYYX(i+1))/2
  IZZX(i)=(IIZZX(i)+IIZZX(i+1))/2
*enddo                          	!DOѭֵ
!
*do,i,1,42                        	!ϲԪʵ
  r,i,SS(i),IZZS(i),IYYS(i),bS(i),hS(i)
*enddo
!
*do,i,43,114                      	!²Ԫʵ
  r,i,SX(i-42),IZZX(i-42),IYYX(i-42),bX(i-42),hX(i-42)
*enddo
!ǵԪʵ
r,115,1,,,,,                        	!-1Ԫʵ
r,116,4,,,,,                        	!-2Ԫʵ
r,117,1,,,,,                        	!-3Ԫʵ
!Ԫʵ
*dim,cable_strain_bs,,2              	!屳Ԫʵ
cable_strain_bs(1)=3.978873577e-3    	!ʼӦ䣨ԤӦ2500
cable_strain_bs(2)=3.978873577e-3    	!ʼӦ䣨ԤӦ2500
!
R,118,0.031416,cable_strain_bs(1)     	!屳Ԫʵ
R,119,0.031416,cable_strain_bs(2)
! 
*dim,cable_strain_zs,,11               !Ԫʵ
cable_strain_zs(1)=2.419155135e-3      !ʼӦ䣨ԤӦ380
cable_strain_zs(2)=2.419155135e-3      !ʼӦ䣨ԤӦ380
cable_strain_zs(3)=2.2918311805e-3     !ʼӦ䣨ԤӦ360
cable_strain_zs(4)=2.2918311805e-3     !ʼӦ䣨ԤӦ360
cable_strain_zs(5)=2.228169203e-3      !ʼӦ䣨ԤӦ350
cable_strain_zs(6)=2.228169203e-3      !ʼӦ䣨ԤӦ350
cable_strain_zs(7)=1.909859317e-3      !ʼӦ䣨ԤӦ300
cable_strain_zs(8)=1.65521141e-3       !ʼӦ䣨ԤӦ260
cable_strain_zs(9)=1.65521141e-3       !ʼӦ䣨ԤӦ260
cable_strain_zs(10)=1.40056340e-3      !ʼӦ䣨ԤӦ220
cable_strain_zs(11)=1.273239545e-3     !ʼӦ䣨ԤӦ200
!
*do,i,120,130                        !ʵ
  R,i,0.007853982,cable_strain_zs(i-119)  
*enddo
!ӵԪʵ
r,131,16,21.3333,21.3333,4,4           !Ԫʵ
!ϲ
mp,ex,1,3.4e10                      !ϲ
mp,prxy,1,0.3
mp,dens,1,2.5e3
!
MP,EX,2,3.5e10                     !
mp,prxy,2,0.3
MP,DENS,2,2.5e3
!
MP,EX,3,2e11                       !屳
mp,prxy,3,0.3
MP,DENS,3,7.85e3
!**********************************
!ɹؼ
!**********************************
!Ĺؼ
*do,i,54.5,33.5,-0.5                 	  !ϲؼ
  k,110-i*2,-(82.1495+i)/7.403322,i,0
*enddo
!
*do,i,33,0,-0.5                      !²ؼ
  k,110-i*2,-15.6213*i/33.5,i,7.4839*i/33.5-7.4839
*enddo
*do,i,-0.5,-2.5,-0.5
  k,110-i*2,-2.7978*i/6,i,-i/2-7.4839
*enddo
*do,i,33,0,-0.5
  k,182-i*2,-15.6213*i/33.5,i,-7.4839*i/33.5+7.4839
*enddo
*do,i,-0.5,-2.5,-0.5
  k,182-i*2,-2.7978*i/6,i,i/2+7.4839
*enddo
!ؼ(Ϊ)
k,192,-2,0,-6                        !ؼ
k,193,0,0,-6
k,194,2,0,-6
k,195,9,0,-6
!
*do,i,196,205                       !ؼ
  k,i,(i-195)*9+9,0,-6
*enddo
k,206,105,0,-6
!
k,207,-2,0,6                        !ؼ
k,208,0,0,6
k,209,2,0,6
k,210,9,0,6
!
*do,i,211,220                       !ؼ
  k,i,(i-210)*9+9,0,6
*enddo
k,221,105,0,6
!屳êֹؼ
k,222,-31,0,14
k,223,-31,0,-14
!**********************************
!ɵԪ
!********************************** 
*do,i,1,42                       !ϲ
  l,i,i+1                        !ؼ֮  
  lsel,u,type,,1,131             !ѡѶԵ
  latt,1,i,1,,,,                 !ָߵĲϱšʵšԪͺ
  lesize,all,,,1,,1              !趨ߵķָ
  lmesh,all,all                  !߻
  alls                           !ѡȫͼԪ
*enddo                           !DOѭ־ 
!
*do,i,43,114                     !²
  l,i,i+1
  lsel,u,type,,1,131
  latt,1,i,1,,,,
  lesize,all,,,1,,1
  lmesh,all,all
  alls
*enddo
!
l,43,116                         	!²Զ
lsel,u,type,,1,131
latt,1,43,1,,,,
lesize,all,,,1,,1
lmesh,all,all
alls
!
*do,i,116,186
  l,i,i+1
  lsel,u,type,,1,131
  latt,1,i-72,1,,,,
  lesize,all,,,1,,1
  lmesh,all,all
  alls
*enddo
!
l,9,222                             !ֱ
lsel,u,type,,1,131
latt,3,118,3,,,,
lesize,all,,,1,,1
lmesh,all,all
alls
!
l,9,223
lsel,u,type,,1,131
latt,3,119,3,,,,
lesize,all,,,1,,1
lmesh,all,all
alls
!
*do,i,1,11                           !
  l,i*4+1,209-i
  l,i*4+1,221-i
  lsel,u,type,,1,131
  latt,3,i+119,3,,,,
  lesize,all,,,1,,1
  lmesh,all,all
  alls
*enddo
!
l,110,193                           !Ӵ
l,182,208
lsel,u,type,,1,131
latt,3,131,1,,,,
lesize,all,,,1,,1
lmesh,all,all
alls
!
wproat,0,-90,0                                !
cyl4,-5.6695,-70,76,90,64,118.09989655,0
wproat,0,90,0
!
a,188,192,207,191                             !
lsel,u,type,,1,131
aatt,2,115,2,,,,
lesize,all,2,,,,2
lesize,213,2,,,,2
lesize,214,2,,,,2
lesize,215,2,,,,2
lesize,216,2,,,,2
lesize,217,,,,4,2
lesize,218,2,,,,2
lesize,219,,,,4,2
alls
amesh,all
!
a,192,193,208,207                              !潻紦
a,193,194,209,208
lsel,u,type,,1,131
aatt,2,116,2,,,,
lesize,all,2,,,,1
alls
amesh,all
!
*do,i,1,12                                      !
  a,i+193,i+194,i+209,i+208
*enddo
lsel,u,type,,1,131
aatt,2,117,2,,,,
lesize,all,2,,,,1
alls
amesh,all
save
finish
!**********************************
!ء
!**********************************
/SOLU
d,115,all     			            !ڽڵʩλԼ
d,187,all
d,188,all
d,189,all
d,235,uy
d,239,uy
d,280,uy
d,292,uy
d,735,uy
d,739,uy
acel,,10      			            !ʩٶ
!**********************************
!ģ̬ģ
!**********************************
antype,modal               	   !趨ģ̬
modopt,subsp,20                    !趨ģ̬ѡ
solve
save,rst11,db
finish
!**********************************
!Ӧ׷ģ
!**********************************
/solu         			
Antype,spectr                       !趨׷
Spopt,sprs,20,yes                   !趨׷ѡ
Svtype,2                            !趨Ӧ
Sed,1,1,                            !趨
Freq,0.1181,0.1228,0.1436,0.1576,0.1629,0.1639,0.1761,0.2197,0.2273,0.2358,0.2588,0.2674,0.3209,0.3478,0.3739,0.4002,0.4938,0.7503,1.3831,1.5370   !趨Ƶֵ
Sv,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.07996,0.0662,0.0454,0.0262,0.0238                                  !趨Ƶֵ
Solve
Save,rst2,db
Finish
!******ģ̬չ******
/solu      					    
Antype,modal                      !趨ģ̬
Expass,on                         !ģ̬չ
Mxpand,20,,,yes,0.005             !趨ģ̬չѡ
Solve
Save,rst3,db
Finish
!*******ϲģ̬******
/solu     					     
Antype,spectr                        !趨׷
Srss,0.02,disp                       !趨ģ̬ϲѡ
Solve
Save,rst4,db
Finish
/post26                           	    !ʱ̺鿴
Nsol,2,2,u,x,ux_2          		    !ڵ2 XλΪ2
Nsol,3,2,u,y,uy_3          		    !ڵ2 YλΪ3
NSOL,4,400,U,X, UX_4    		    !ڵ400 XλΪ4
NSOL,5,400,U,Y, UY_5    		    !ڵ400 YλΪ5
STORE,MERGE                   	   !洢ı
PLVAR,2,3,4,5,                     !ͼʾ2345
Finish
!**********************************
!˲̬ģ
!**********************************
*DIM, horizon, ARRAY,2,50,0, , ,        !ʼhorizon
*DIM, vertical, ARRAY,2,50,0, , ,       !ʼvertical
*CREATE, ansuitmp                       !꿪ʼʶ
*VREAD,horizon(1,1),'horizon','txt',' ',50, , , , , , !horizon.txtļݵhorizon
(e9.3,e11.3)                             !趨ݶʽ
*END                                     !ʶ 
/INPUT,ansuitmp   
*CREATE,ansuitmp 
*VREAD,vertical(1,1),'vertical','txt',' ',50, , , , , , !vertical.txtļݵvertical
(e9.3,e11.3)
*END
/SOL                       	                !˲̬ģ
ANTYPE,4					!÷Ϊ˲̬
TRNOPT,FULL				        !趨ⷽΪFull
*do,t,1,50,1					!ѭ⿪ʼ־ʱΪ5s
time,0.1*t					!ÿزʱ䣬ʱΪ0.1
kbc,0						!üزļʽΪramp
nsub,10						!ÿزӲ
acel,horizon (2,t),vertical (2,t)		!ʩˮƽٶ
allsel,all                      	        !ѡȫͼԪ 
solve						!
*enddo                     !ѭ־
FINISH                     !ģ־
/POST26                    !ʱ̺鿴
Nsol,2,2,u,x,ux_2          !ڵ2 XλΪ2
Nsol,3,2,u,y,uy_3          !ڵ2 YλΪ3
NSOL,4,400,U,X, UX_4       !ڵ400 XλΪ4
NSOL,5,400,U,Y, UY_5       !ڵ400 YλΪ5
STORE,MERGE                !涨ı
PLVAR,2,4                  !ʾڵ2ڵ400XλӦ
PLVAR,3,5                  !ʾڵ2ڵ400YλӦ
Finish                     !ģ 
/EXIT,ALL                  !˳
