new
restore model.sav;初始模型中结构单元节点与土体间无link
set echo on
set plot jpg
set plot jpg si 1600 1000
plot set ba wh
plot set cap right
plot set cap si 25
plot set rot 30 0 30
sel group shot1 range group shot_concrete z=(-3.25,0.01)
sel group shot2 range group shot_concrete z=(-4.88,-3.25)
sel group shot3 range group shot_concrete z=(-6.4,-4.88)
group excavation range group 12 any group 13 any group 14 any group 15 any
group soil_out range group 16
group exca1 range group excavation z=(-3.25,0.01)
group exca2 range group excavation z=(-4.88,-3.25)
group exca3 range group excavation z=(-6.4,-4.88)
;---------土层分组---------
range name s1 z=(-6.5,0)
;---------initial state---------
;model eals
;prop y=40e6 p=0.3 density=1.8e3
model mohr
prop y=40e6 p=0.35 c=8e3 f=10.0 t=1e3 density=1.78e3 range s1
;----------------------------------
fix x range x=(-150.1,-149.9) any x=(149.9,150.1) any
fix y range y=(-0.1 0.1) any y=(149.9 150.1) any
fix z range z=(-60.1 -59.9)
set gravity=(0.0,0.0,-10.0)
ini szz 0.0 gradient=(0.0,0.0,18e3)
ini sxx 0.0 gradient=(0.0,0.0,9e3)
ini syy 0.0 gradient=(0.0,0.0,9e3)
;------------施加地面荷载15kpa超载---------------
sel shell apply pres -15e3 range group shell_load
sel shell apply pres -15e3 range group shot1 z=(-0.01,0.01)
;----------------------------------------
cyc 1
sel node fix 1sys
;-------单元节点坐标确定顺序liner,geogrid,pile,cable,shell and beam
;------------重新建立link---------------
sel node ini xpos add 0.0
;---------结构边界----------
sel node fix y xr zr range y=(-0.01,0.01)
;--------结构单元赋初参数--------
sel shell prop iso=(1e-10.0.2) density=0.0
sel cable prop emod=1e-10 density=0.0 gr coh=1e20 gr_fri=14.0 gr_k=24e9&
gr_per=0.314 slide=on slide_tol=1e6 ycomp=1e20 ytens=1e20
;-----------------------------------
set mech damp combined
sel set damp combined
hist unbal
p1 hist 1
step 3000
save initial.sav
;-------------开挖计算-------------
new
res initial.sav
ini xdisp=0.0 ydisp=0.0 zdisp=0.0
ini xvel=0.0 yvel=0.0 zvel=0.0
sel node ini xdisp=0.0 ydisp=0.0 zdisp=0.0
sel node ini xvel=0.0 yvel=0.0 zvel=0.0
ini state=0
;-------------阳角坡顶第4网格点位移-------------
hist id 10 sel node ydisp id 26950
hist id 11 sel node ydisp id 26934
;-------------A-A剖面第1排土钉第一单元轴力-------------
hist id 101 sel cablesel force cid 178710
hist id 102 sel cablesel force cid 179324
hist id 103 sel cablesel force cid 179872
hist id 104 sel cablesel force cid 180493
hist id 105 sel cablesel force cid 180931
hist id 106 sel cablesel force cid 181312
hist id 107 sel cablesel force cid 181690
hist id 108 sel cablesel force cid 182025
hist id 109 sel cablesel force cid 182331
;----------stage1----------
model null range group exca1
step 2000
sel cable prop emod=210e9 density=7.8e3 gr_coh=3.14e4 gr_fri=10.0 gr_k=24e9 gr_per=0.314&
slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6 range group nail1
sel shell prop iso=(10e9,0.2) density=2.5e3 thick=0.1 ran group shot1
save stage1.sav
;----------stage2----------
model null range group exca2
step 2000
sel cable prop emod=210e9 density=7.8e3 gr_coh=3.14e4 gr_fri=10.0 gr_k=24e9 gr_per=0.314&
slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6 range group nail2
sel shell prop iso=(10e9,0.2) density=2.5e3 thick=0.1 ran group shot2
save stage2.sav
;----------stage3----------
model null range group exca3
step 2000
sel shell prop iso=(10e9,0.2) density=2.5e3 thick=0.1 ran group shot3
save stage3.sav
