level 4
金慧之父
楼主
new
gen zon radcyl p0 0 0 0 p1 12.55 0 0 p2 0 48 0 p3 0 0 12.55 &
size 5 96 10 4 dim 5.25 5.25 5.25 5.25 rat 1 1 1 1.2 group outsiderocksy
gen zon cshell p0 0 0 0 p1 5.25 0 0 p2 0 48 0 p3 0 0 5.25 &
size 1 96 10 4 dim 4.75 4.75 4.75 4.75 rat 1 1 1 1 group concretlinersy fill group insiderocksy
gen zon radcyl p0 0 0 0 p1 0 0 -10 p2 0 48 0 p3 12.55 0 0 &
size 5 96 10 4 dim 4.25 5.25 4.25 5.25 rat 1 1 1 1.2 group outsiderockxy
gen zon cshell p0 0 0 0 p1 0 0 -4.25 p2 0 48 0 p3 5.25 0 0 &
size 1 96 10 4 dim 3.75 4.75 3.75 4.75 rat 1 1 1 1.2 group concretlinerxy fill group insiderockxy
gen zon brick p0 0 0 12.55 p1 12.55 0 12.55 p2 0 48 12.55 p3 0 0 35.25 &
size 5 96 8 group outsiderock1
gen zon brick p0 0 0 -32.75 p1 12.55 0 -32.75 p2 0 48 -32.75 p3 0 0 -10 &
size 5 96 12 group outsiderock2
gen zon brick p0 12.55 0 0 p1 62.25 0 0 p2 12.55 48 0 p3 12.55 0 12.55 &
size 20 96 5 group outsiderock3
gen zon brick p0 12.55 0 -10 p1 62.25 0 -10 p2 12.55 48 -10 p3 12.55 0 0 &
size 20 96 8 group outsiderock3
gen zon brick p0 12.55 0 12.55 p1 62.25 0 12.55 p2 12.55 48 12.55 p3 12.55 0 35.25 size 20 96 8 group outsiderock4
gen zon brick p0 12.55 0 -32.75 p1 62.25 0 -32.75 p2 12.55 48 -32.75 p3 12.55 0 -10 size 20 96 12 group outsiderock5
gen zon reflect dip 90 dd 90 orig 0 0 0
pl bl g
set gravity 0 0 -10
fix z range z -32.76,-32.74
fix x range x -62.26,-62.24
fix x range x 62.24,62.26
fix y range y -0.01,0.01
fix y range y 47.99,48.01
model mohr
ini density 1800
prop bulk=7.01e8 shear=2.8e8 fric=23 coh=3e5 tension=0.8e5
step 3000
save Gsol.sav
def abc
loop i(1,14)
if i>0 then
if i<13 then
x0=4*i-4
x1=4*i
command
model null range group insiderocksy y x0 x1
model null range group concretlinersy y x0 x1
sel shell id=1 range x -5.26 5.26 z -0.01 5.26 cyl end1 0 x0 5.25 end2 0 x1 5.25 rad 5.25
sel shell id=1 prop iso=(23e9 0.2) thick=0.25
plot block group
end_command
end_if
end_if
if i>1 then
if i<15 then
x2=4*i-8
x3=4*i-4
command
mo null range group insiderockxy y x2 x3
mo null range group concretlinerxy y x2 x3
sel shell id=1 range x -5.26 5.26 z -4.26 0.01 cyl end1 0 x2 -4.25 end2 0 x3 4.25
sel shell id=1 prop iso=(23e9 0.2) thick=0.25
plot block group
end_command
end_if
end_if
command
model elas range group concretlinersy y x2 x3
model elas range group concretlinerxy y x2 x3
ini density 2600 range group concretlinersy y x2 x3
ini density 2600 range group concretlinerxy y x2 x3
prop bulk=12.8e9 shear=16.67e9 range group concretlinersy y x2 x3
prop bulk=12.8e9 shear=16.67e9 range group concretlinerxy y x2 x3
plot block group
end_command
endloop
end
abc
step 3000
save step1.sav
plot add axes red
plot cont zdisp
plot cont xdisp
plot cont szz
plot cont sxx
plot cont szz range group concretlinersy any group concretlinerxy any
plot cont sxx range group concretlinersy any group concretlinerxy any
plot add axes red
plot block state
save step2.sav
2017年06月12日 01点06分
1
gen zon radcyl p0 0 0 0 p1 12.55 0 0 p2 0 48 0 p3 0 0 12.55 &
size 5 96 10 4 dim 5.25 5.25 5.25 5.25 rat 1 1 1 1.2 group outsiderocksy
gen zon cshell p0 0 0 0 p1 5.25 0 0 p2 0 48 0 p3 0 0 5.25 &
size 1 96 10 4 dim 4.75 4.75 4.75 4.75 rat 1 1 1 1 group concretlinersy fill group insiderocksy
gen zon radcyl p0 0 0 0 p1 0 0 -10 p2 0 48 0 p3 12.55 0 0 &
size 5 96 10 4 dim 4.25 5.25 4.25 5.25 rat 1 1 1 1.2 group outsiderockxy
gen zon cshell p0 0 0 0 p1 0 0 -4.25 p2 0 48 0 p3 5.25 0 0 &
size 1 96 10 4 dim 3.75 4.75 3.75 4.75 rat 1 1 1 1.2 group concretlinerxy fill group insiderockxy
gen zon brick p0 0 0 12.55 p1 12.55 0 12.55 p2 0 48 12.55 p3 0 0 35.25 &
size 5 96 8 group outsiderock1
gen zon brick p0 0 0 -32.75 p1 12.55 0 -32.75 p2 0 48 -32.75 p3 0 0 -10 &
size 5 96 12 group outsiderock2
gen zon brick p0 12.55 0 0 p1 62.25 0 0 p2 12.55 48 0 p3 12.55 0 12.55 &
size 20 96 5 group outsiderock3
gen zon brick p0 12.55 0 -10 p1 62.25 0 -10 p2 12.55 48 -10 p3 12.55 0 0 &
size 20 96 8 group outsiderock3
gen zon brick p0 12.55 0 12.55 p1 62.25 0 12.55 p2 12.55 48 12.55 p3 12.55 0 35.25 size 20 96 8 group outsiderock4
gen zon brick p0 12.55 0 -32.75 p1 62.25 0 -32.75 p2 12.55 48 -32.75 p3 12.55 0 -10 size 20 96 12 group outsiderock5
gen zon reflect dip 90 dd 90 orig 0 0 0
pl bl g
set gravity 0 0 -10
fix z range z -32.76,-32.74
fix x range x -62.26,-62.24
fix x range x 62.24,62.26
fix y range y -0.01,0.01
fix y range y 47.99,48.01
model mohr
ini density 1800
prop bulk=7.01e8 shear=2.8e8 fric=23 coh=3e5 tension=0.8e5
step 3000
save Gsol.sav
def abc
loop i(1,14)
if i>0 then
if i<13 then
x0=4*i-4
x1=4*i
command
model null range group insiderocksy y x0 x1
model null range group concretlinersy y x0 x1
sel shell id=1 range x -5.26 5.26 z -0.01 5.26 cyl end1 0 x0 5.25 end2 0 x1 5.25 rad 5.25
sel shell id=1 prop iso=(23e9 0.2) thick=0.25
plot block group
end_command
end_if
end_if
if i>1 then
if i<15 then
x2=4*i-8
x3=4*i-4
command
mo null range group insiderockxy y x2 x3
mo null range group concretlinerxy y x2 x3
sel shell id=1 range x -5.26 5.26 z -4.26 0.01 cyl end1 0 x2 -4.25 end2 0 x3 4.25
sel shell id=1 prop iso=(23e9 0.2) thick=0.25
plot block group
end_command
end_if
end_if
command
model elas range group concretlinersy y x2 x3
model elas range group concretlinerxy y x2 x3
ini density 2600 range group concretlinersy y x2 x3
ini density 2600 range group concretlinerxy y x2 x3
prop bulk=12.8e9 shear=16.67e9 range group concretlinersy y x2 x3
prop bulk=12.8e9 shear=16.67e9 range group concretlinerxy y x2 x3
plot block group
end_command
endloop
end
abc
step 3000
save step1.sav
plot add axes red
plot cont zdisp
plot cont xdisp
plot cont szz
plot cont sxx
plot cont szz range group concretlinersy any group concretlinerxy any
plot cont sxx range group concretlinersy any group concretlinerxy any
plot add axes red
plot block state
save step2.sav
