5-7.txt

是用apdl语言编的

/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