/TITLE,cable-stayed bridge
/prep7
et,1,beam4
et,2,link10
r,1,25.6,20,2000,16,1.6,,,21.8
r,2,16,39.7,77.7,4.7,3.4,,,61.6
r,3,40,213.3,83.3,5,8,,,333.3
r,4,54,450,200,6,9,,,648
r,5,1,0.08333333,0.08333333,1,1,,,1/3
r,6,0.12,0
mp,ex,1,3.5e10
mp,prxy,1,0.17
mp,dens,1,2500
mp,ex,2,10e15
mp,prxy,2,0
mp,dens,2,0
mp,ex,3,1.9e10
mp,prxy,3,0.25
mp,dens,3,1200
*do,i,1,59
x=-174*2+(i-1)*6
y1=-14
y2=14
n,3*(i-1)+1,x   ! 1,4,7,...
n,3*(i-1)+2,x,y1 ! 2,5,8,...
n,3*i,x,y2     ! 3,6,9,...
*end do
type,1
real,1
mat,1
*do,i,1,58
j=3*(i-1)+1
e,j,j+3
*end do
type,1
real,5
mat,2
*do,i,1,59
j=3*(i-1)+1
j1=3*(i-1)+2
j2=3*i
e,j,j1
e,j,j2
*end do
i=59*3
n,i+1,-174,-10,-30
n,i+2,-174,10,-30
n,i+3,-174,-15
n,i+4,-174,15
*do,j,1,5
k=i+4+j
n,k,-174,0,60+(j-1)*18
*end do
type,1
real,4
mat,1
e,i+1,i+3
e,i+2,i+4
type,1
real,3
mat,1
e,i+3,i+5
e,i+4,i+5
type,1
real,2
mat,1
*do,j,1,4
k=i+4+j
e,k,k+1
*end do
!>>>>>>>>>>>>>>>>>>>>>>>>>>
type,2
real,6
mat,3
e,i+5,89
e,i+5,90
type,2
real,6
mat,3
*do,j,1,8
e,i+6,89+3*j
e,i+6,89-3*j
e,i+6,90+3*j
e,i+6,90-3*j
*end do
*do,k,1,3
*do,j,1,7
e,i+6+k,113+(k-1)*21+3*j
e,i+6+k,65-(k-1)*21-3*j
e,i+6+k,114+(k-1)*21+3*j
e,i+6+k,66-(k-1)*21-3*j
*end do
*end do
i=i+9
nsym,x,i,all
esym,,i,all
nummrg,all
type,1
real,5
mat,2
j=i
i=i-9
n,1000,-174
e,1000,i+3
e,1000,i+4
n,2000,174
e,2000,i+3+j
e,2000,i+4+j
cpintf,uy
cpintf,uz
cpintf,rotx
cpintf,rotz
nsel,,loc,z,-30
d,all,all
allsel
nsel,,loc,x,-348
nsel,r,loc,y,0
d,all,uy
d,all,uz
allsel
nsel,,loc,x,348
nsel,r,loc,y,0
d,all,uy
allsel
NUMCMP,all
/solu
solve
/prep7
r,6,0.012,0.012
/solu
lumpm,on
pstres,on
solve
finish
upcoord,1.0,on
save
/solu
antype,modal
modopt,lanb,20
mxpand,10,0,100
lumpm,on
pstres,on
solve
finish

!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
/resume
!将生成的地震数据文件放到工作目录下
*dim,AAZ1,array,2,50
*dim,AAY1,array,2,50
*CREATE,ansuitmp
*vread,AAZ1(1,1),ACELZ,txt,,50
(e9.3 , e11.3)
*END
/INPUT,ansuitmp
*CREATE,ansuitmp
*vread,AAY1(1,1),ACELY,txt,,50
(e9.3 , e11.3)
*END
/INPUT,ansuitmp
N=117
*dim,fcar,array,N
*do,i,1,N
tt=(i-1)*0.1
fcar(i)=1000*cos(10*tt)
*enddo
*dim,fwind,array,N,N
*do,i,1,N
tt=(i-1)*0.1
p=50*sin(1.5*tt)
fwind(1,i)=3*3*p
*do,j,1,N-2
  fwind(j+1,i)=6*3*p
*enddo
fwind(N,i)=fwind(1,i)
*enddo
FINISH
/config,nres,2000
/solu
antype,trans
M=20
*do,i,1,M
time,i*0.1
auto,off
deltim,0.05
kbc,0
outres,basic
cnvtol,f,,0.01
*if,i,le,50,then
  acel,,-AAY1(2,i),10-AAZ1(2,i)
*ELSE
  ACEL,0,0,10
*endif
*if,i,eq,1,then
  nsel,s,loc,x,-174*2
  nsel,r,loc,y,0
  nsel,r,loc,Z,0
  f,all,fz,-fcar(j)
  allsel
*else
  *do,j,2,M
  nsel,s,loc,x,-174*2+(j-1)*6
  nsel,r,loc,y,0
  nsel,r,loc,Z,0
  f,all,fz,-fcar(j)
  allsel
  nsel,s,loc,x,-174*2+(j-2)*6
  nsel,r,loc,y,0
  nsel,r,loc,Z,0
  fdele,all,fz
  allsel  
  *enddo
*endif
*do,j,1,N
  nsel,s,loc,y,-14
  nsel,s,loc,x,-174*2+(j-1)*6
  f,all,fy,fwind(j,i)
  allsel
*enddo
lswrite,i
*enddo
lssolve,1,M
finish