基坑开挖_实体单元模拟

new
gen zone brick p0 0 0 0 p1 10.35 0 0 p2 0 4.5 0 p3 0 0 60 size 11 4 60
gen zone brick p0 10.35 0 0 p1 11.35 0 0 p2 10.35 4.5 0 p3 10.35 0 60 size 1 4 60
gen zone brick p0 11.35 0 0 p1 150 0 0 p2 11.35 4.5 0 p3 11.35 0 60 size 45 4 60 rat 1.05 1 1
group wall ran x 10.35 11.35 z 24 60
group soil ran gr wall not
m m
prop sh 1.85e7 bu 4.0e7 fric 33 cohesion 1e3 dil 3 ten 1e-3 density 1960 ran z 45 60
prop sh 1.73e7 bu 3.75e7 fric 33 cohesion 1e3 dil 3 ten 1e-3 density 1950 ran z 36 45
prop sh 1.31e7 bu 2.83e7 fric 33 cohesion 1e3 dil 3 ten 1e-3 density 1970 ran z 17 36
prop sh 1.88e7 bu 4.08e7 fric 33 cohesion 1e3 dil 3 ten 1e-3 density 1990 ran z 0 17

def ini_szz
_grad =1960*10
_szzl=-_grad*15
_szz0=-_grad*60
command
ini szz add _szz0 grad 0 0 _grad ran z 45 60
ini szz add _szzl ran z 0 45
endcommand
_grad =1950*10
_szzl=-_grad*9
_szz0=-_grad*45
command
ini szz add _szz0 grad 0 0 _grad ran z 36 45
ini szz add _szzl ran z 0 36
endcommand
_grad =1970*10
_szzl=-_grad*19
_szz0=-_grad*36
command
ini szz add _szz0 grad 0 0 _grad ran z 17 36
ini szz add _szzl ran z 0 17
endcommand
_grad =1990*10
_szzl=-_grad*17
_szz0=-_grad*17
command
ini szz add _szz0 grad 0 0 _grad ran z 0 17
endcommand
end
ini_szz
def ini_sxx_syy
pnt=zone_head
loop while pnt#null
val=0.455*z_szz(pnt)
z_sxx(pnt)=val
z_syy(pnt)=val
pnt=z_next(pnt)
endloop
end
ini_sxx_syy

fix x z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 149.9 150.1
fix y
set gravity 10
so

ini xdis 0 ydis 0 zdis 0
ini xve 0 yve 0 zve 0

m elas ran gr wall
prop bu 8.95e9 sh 8.17e9 ran gr wall
gen separate wall
int 1 face ran group wall x 10.3 10.4
int 1 maxedge 0.5
int 1 prop kn=6.47e7 ks=6.47e7 tens=1.0e-3 coh=0.8e3 fric=22
int 2 face ran group wall x 11.3 11.4
int 2 maxedge 0.5
int 2 prop kn=6.47e7 ks=6.47e7 tens=1.0e-3 coh=0.8e3 fric=22
int 3 face ran group wall z 23.9 24.1
int 3 maxedge 0.5
int 3 prop kn=4.58e7 ks=4.58e7 tens=1.0e-3 coh=0.8e3 fric=22

;;;;;;;;step1;;;;;;;;;;;;
;;;;;;;;lower GWL;;;;;;;;
;apply nstress add -57e4 grad 0 0 10e3 ran x 11.3 11.4 z 56 57 gr wall
;apply nstress add -1.0e4 ran x 11.3 11.4 z 24 56 gr wall
;so
;;;;;;;;;excavation;;;;;;
model null range x -.1 10.4 z 57 60
so
;;;;;;;;;preload;;;;;;;;;
sel beam id 1 begin (0,2.25,57.5) end (10.35,2.25,57.5) nseg 10
sel beam id 1 prop dens=7.9e3 emod=1.26e11 nu=0.3 xcarea=0.01739 xcj=5.31e-4 xciy=3.95e-4 &
xciz=1.36e-4 ydirection=(0 1 0)
sel delete link ran id 1
sel link id 1 2 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid ran id 1
plo sk
plo add sel geom black red cid on scale 0.02
sel node local xdir=(0,1,0) ydir=(0,0,1) range id 1
sel node fix lsys range id 1
sel node fix x xr yr zr range id 1
sel node apply force 0.65e6 0 0 ran id 1
hist id 1 sel beamsel force fx end1 cid 10
so

sel node fix x z xr yr zr range id 1

;;;;;;;;step2;;;;;;;;;;;;
;apply nstress add -56e4 grad 0 0 10e3 ran x 11.3 11.4 z 52 56 gr wall
;apply nstress add -4.0e4 ran x 11.3 11.4 z 24 52 gr wall
;so
model null range x -.1 10.4 z 53 57
so
sel beam id 2 begin (0,2.25,54.1) end (10.35,2.25,54.1) nseg 10
sel beam id 2 prop dens=7.8e3 emod=1.26e11 nu=0.3 xcarea=0.02195 xcj=8.78e-4 xciy=6.54e-4 &
xciz=2.24e-4 ydirection=(0 1 0)
sel delete link ran id 2
sel link id 2 13 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid ran id 2
plo sk
plo add sel geom black red cid on scale 0.02
sel node local xdir=(0,1,0) ydir=(0,0,1) range id 12
sel node fix lsys range id 12
sel node fix x xr yr zr range id 12
sel node apply force 1.25e6 0 0 ran id 12
hist id 2 sel beamsel force fx end1 cid 20
so
sel node fix x z xr yr zr range id 12

;;;;;;;;step3;;;;;;;;;;;;
;apply nstress add -52e4 grad 0 0 10e3 ran x 11.3 11.4 z 49 52 gr wall
;apply nstress add -3.0e4 ran x 11.3 11.4 z 24 49 gr wall
;so
model null range x -.1 10.4 z 50 53
so
sel beam id 3 begin (0,2.25,50.9) end (10.35,2.25,50.9) nseg 10
sel beam id 3 prop dens=1.6e4 emod=1.26e11 nu=0.3 xcarea=0.05924 xcj=2.45e-3 xciy=1.83e-3 &
xciz=6.2e-4 ydirection=(0 1 0)
sel delete link ran id 3
sel link id 3 24 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid ran id 3
plo sk
plo add sel geom black red cid on scale 0.02
sel node local xdir=(0,1,0) ydir=(0,0,1) range id 23
sel node fix lsys range id 23
sel node fix x xr yr zr range id 23
sel node apply force 1.85e6 0 0 ran id 23
hist id 3 sel beamsel force fx end1 cid 30
so
sel node fix x z xr yr zr range id 23

;;;;;;;;step4;;;;;;;;;;;;
;apply nstress add -49e4 grad 0 0 10e3 ran x 11.3 11.4 z 45 49 gr wall
;apply nstress add -4.0e4 ran x 11.3 11.4 z 24 45 gr wall
;so
model null range x -.1 10.4 z 46 50
so
sel beam id 4 begin (0,2.25,47.4) end (10.35,2.25,47.4) nseg 10
sel beam id 4 prop dens=1.6e4 emod=1.26e11 nu=0.3 xcarea=0.05924 xcj=2.45e-3 xciy=1.83e-3 &
xciz=6.2e-4 ydirection=(0 1 0)
sel delete link ran id 4
sel link id 4 35 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid ran id 4
plo sk
plo add sel geom black red cid on scale 0.02
sel node local xdir=(0,1,0) ydir=(0,0,1) range id 34
sel node fix lsys range id 34
sel node fix x xr yr zr range id 34
sel node apply force 1.55e6 0 0 ran id 34
hist id 4 sel beamsel force fx end1 cid 40
so
sel node fix x z xr yr zr range id 34

;;;;;;;;step5;;;;;;;;;;;;
;apply nstress add -45e4 grad 0 0 10e3 ran x 11.3 11.4 z 42 45 gr wall
;apply nstress add -3.0e4 ran x 11.3 11.4 z 24 42 gr wall
;so
model null range x -.1 10.4 z 43 46
so
sel beam id 5 begin (0,2.25,44.2) end (10.35,2.25,44.2) nseg 10
sel beam id 5 prop dens=1.6e4 emod=1.26e11 nu=0.3 xcarea=0.05924 xcj=2.4

5e-3 xciy=1.83e-3 &
xciz=6.2e-4 ydirection=(0 1 0)
sel delete link ran id 5
sel link id 5 46 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid ran id 5
plo sk
plo add sel geom black red cid on scale 0.02
sel node local xdir=(0,1,0) ydir=(0,0,1) range id 45
sel node fix lsys range id 45
sel node fix x xr yr zr range id 45
sel node apply force 1.9e6 0 0 ran id 45
hist id 5 sel beamsel force fx end1 cid 50
so
sel node fix x z xr yr zr range id 45

;;;;;;;;step6;;;;;;;;;;;;
;apply nstress add -42e4 grad 0 0 10e3 ran x 11.3 11.4 z 39 42 gr wall
;apply nstress add -3.0e4 ran x 11.3 11.4 z 24 39 gr wall
;so
model null range x -.1 10.4 z 40 43
so
plo sk
plo add sel geom black red cid on scale 0.02