ET,1,PIPE59,2
ET,2,link10
ET,3,PIPE59,2
ET,4,COMBIN14 !竖向弹簧
ET,5,COMBIN40 !横向弹簧
ET,6,COMBIN40
ET,7,COMBIN40
MP,DENS,1,2968
MP,DENS,2,2300
MP,DENS,3,100
MP,EX,1,2.6E+10
MP,EX,2,2.1E+11
MP,EX,3,2.0E+13
MP,EX,4,2.1E+11
mp,damp,4,0.05
KEYOPT,4,2,3 ! 定义方向
KEYOPT,5,3,2
KEYOPT,6,3,3
KEYOPT,7,3,2
MP,PRXY,1,0.3
MP,PRXY,2,0.3
MP,PRXY,3,0.3
TB,WATE,1,,,0
TBMODIF,1,3,400
TBMODIF,1,4,1080
TB,WATE,2,,,0
TBMODIF,1,3,400
TBMODIF,1,4,1080
R,1,15.95,1.32,1,2,0,0,0,
R,2,0.12, !张拉预应力
!R,2,0.62,0.000
R,3,3,1,1,2
!R,4,6E8,0.1 !竖向弹簧刚度
R,4,3*2*1.16e8,0.1
R,5,0.9*2*2*7.8e7,0.15,10,,0.2*1.7e7,0.1*2*2*7.8e7 !双线性弹簧相关参数combin 40
R,6,3.0e8,0.15,10,,3203000,2.373E7
R,7,3.0E8,0.15,10,,3203000,2.373E7
建模
n,1,0,0,-30
n,651,4680,0,-30
n,1201,0,0,-30
n,1202,4680,0,-30
n,701,72,7.3,-30
n,764,4608,7.3,-30
n,801,72,100,-180
n,864,4608,100,-180
n,901,72,-7.3,-30
n,964,4608,-7.3,-30
n,1001,72,-100,-180
n,1064,4608,-100,-180
fill,1,651
fill,701,764
fill,801,864
fill,901,964
fill,1001,1064
finish
/prep7
type,1
real,1
mat,1
*do,i,1,650
en,i,i,i+1
*enddo
finish
!建立悬索模型
/prep7
type,2
mat,2
real,2
*do,i,1,64
en,700+i,700+i,800+i
en,900+i,900+i,1000+i
*enddo
finish
/prep7
type,3
mat,3
real,3
*do,i,1,64
e,700+i,10*i+1
e,900+i,10*i+1
*enddo
finish
!建立弹簧模型
/prep7
type,4
real,4
mat,4
en,1201,1201,1
en,1202,1202,651
/prep7
type,5
real,5
mat,4
en,1203,1201,1
en,1204,1202,651
时辰分析
!分析 地震 载荷
finish
/config,nres,5000
FREQ1=0.18
/solu
antype,trans
pstres,on
nlgeom,on
!sstif,on
DAMPRATIO=0.02
ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926
ALPHAD,0
BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)! 结构阻尼
*do,i,1,1000
*do,k,1,2
fdele,all,all
kbc,0
kk=i*2+k-1
time,0.02*kk
convtol,f,,0.01
NSUBST,10,0,0
OUTRES,ERASE
OUTRES,BASI,LAST
acel,,acelt(kk,1)*0.001,acelv(kk,1)*0.001+10 !施结构加惯性力
f,2,fz,fz1(kk,1),,650 ! 施加水动力惯性
f,2,fz,fz1(i,2),,650 !施加水拈性力
solve
*enddo
*enddo |