隧道地震响应命令

finish
/clear
/prep7

!创建几何模型
k,1,0,2.7388794
k,2,-2.0711036,1.7033277
k,3,-2.6388794,0
k,4,-2.3359835,-1.9455903
k,5,-1.6736199,-2.6205221
k,6,0,-2.8972106
k,1000,,,
larc,1,2,1000,2.589
larc,2,3,1000,2.839
larc,3,4,1000,6.4
larc,4,5,1000,1.05
larc,5,6,1000,5.2
lsymm,x,all
nummrg,all
numcmp,all
lcomb,1,2
lcomb,6,7
nummrg,all
numcmp,all
lcomb,3,4
lcomb,7,8
nummrg,all
numcmp,all
lcomb,2,3
lcomb,5,6
nummrg,all
numcmp,all

al,all
blc4,-36,-69,72,84
asba,2,1,,,keep
numcmp,all

!切分
asel,s,,,2
wprota,,,90
asbw,all
wpoffs,,,6 !沿Z轴正方向移动6
asbw,all
wpoffs,,,-12 !沿Z轴负方向移动6
asbw,all
wpcsys,-1 !将工作面恢复到初始状态
wprota,,90
asbw,all
wpoffs,,,6
asbw,all
wpoffs,,,-12
asbw,all
wpcsys,-1
alls
nummrg,all
numcmp,all

!定义材料属性

*set,th1,0.3
et,1,beam3
et,2,plane42
keyopt,2,3,2

r,1,th1,th1*th1*th1/12,th1,,

mp,ex,1,31e9 !!二衬参数
mp,prxy,1,0.3
mp,dens,1,2500

mp,ex,2,1.40e8 !!北京地层
mp,prxy,2,0.3
mp,dens,2,2049
tb,dp,2
tbdata,1,3.4e4,30

!赋予属性
lsel,s,,,1,4
latt,1,1,1
asel,all
aatt,2,,2

!隧道周围划分网格
lsel,s,,,1,3,2
lesize,all,,,8
lsel,s,,,2,4,2
lesize,all,,,8
lsel,s,,,1,4
lmesh,all

lsel,s,,,24,25
lsel,a,,,35
lsel,a,,,44
lesize,all,,,8,3
lsel,s,,,29,30
lesize,all,,,4
lsel,a,,,32,34,2
lesize,all,,,4
lsel,s,,,38,39
lsel,a,,,41,43,2
lesize,all,,,4

asel,s,,,1
amesh,all
alls
amap,16,1,29,19,2
amap,17,4,17,29,1
amap,12,3,24,19,2
amap,13,4,17,24,3

!横线
lsel,s,,,15
lsel,a,,,18
lesize,all,,,4
lsel,a,,,16
lsel,a,,,19
lesize,all,,,4
lsel,s,,,12
lsel,a,,,13
lsel,a,,,23
lsel,a,,,27
lsel,a,,,36
lesize,all,,,20
lsel,s,,,5
lsel,a,,,7
lsel,a,,,11
lsel,a,,,14
lsel,a,,,20
lesize,all,,,20

!竖线
lsel,s,,,9,10
lsel,a,,,37
lsel,a,,,21,22
lesize,all,,,6

lsel,s,,,42
lsel,a,,,40
lesize,all,,,4

lsel,s,,,33
lsel,a,,,31
lesize,all,,,4

lsel,s,,,17
lsel,a,,,26
lsel,a,,,28
lsel,a,,,8
lsel,a,,,6
lesize,all,,,42

alls
asel,u,,,1
asel,u,,,12,13
asel,u,,,16,17
amesh,all
alls

acel,,9.8

asel,s,,,1 !隧道部分
esla,s
ekill,all
alls

nsel,s,loc,x,-36
d,all,uy,0
nsel,s,loc,x,36
d,all,uy,0
nsel,s,loc,y,-69
d,all,all
alls

!读入加速度时程
alls
*DIM,shuiping,ARRAY,2,1000,0, , ,
!*
*CREATE,ansuitmp
*VREAD,shuiping(1,1),'shuiping','txt','E:\',1000, , , ,
(e9.3,e11.3)
*END
/INPUT,ansuitmp

!!!!!计算
/solu
antype,4
trnopt,full
lumpm,0
!autots,on !使用自动时间步
!pred,on !打开时间步长预测器
!lnsrch,on !打开线性搜索
!nlgeom,on !打开大位移效果
!nropt,full !设定牛顿-拉普森选项
!cnvtol,f,,0.02,2,0.5 !设置收敛项
!CUTCONTROL,DSPLIMIT,1000000000
!NCNV,2, !计算不收敛时不退出
!NEQIT,40 !每一子步中方程的迭代次数(默认是25)
!SOLCONTROL,OFF, , !

开关 不考虑位移收敛

*do,t,1,75,1
time,0.02*t
kbc,1
nsub,1
acel,shuiping(2,t)
alphad,0.069 !对结果影响不大
betad,0.035
solve
*enddo

*do,t,76,150,1
time,0.02*t
kbc,1
nsub,50
acel,shuiping(2,t)
alphad,0.069 !对结果影响不大
betad,0.035
solve
*enddo
*do,t,151,750,1
time,0.02*t
kbc,1
nsub,600
acel,shuiping(2,t)
alphad,0.069 !对结果影响不大
betad,0.035
solve
*enddo
save
!finish




!竖向
/post26
/axlab,y,uy !定义坐标轴!y轴为uy
nsol,2,2811,u,y,uy_2 !显示点的时间历程变化曲线plvar,2
prvar,2 ！数据文本
plvar,2 ！图象

/axlab,y,vy
deriv,3,2,1,,,,,1
plvar,3
prvar,3
/axlab,y,ay
deriv,4,3,1,,,,,1
plvar,4
prvar,4
plvar,2,3,4
eep,all



!水平
/post26
/axlab,y,ux
nsol,2,2811,u,x,ux_2
plvar,2 !图象
prvar,2 !数据
/axlab,y,vx
deriv,3,2,1,,,,,1
plvar,3
prvar,3
/axlab,y,ax
deriv,4,3,1,,,,,1
plvar,4
prvar,4
plvar,2,3,4
!eep,all