new
gen zon cshell p0 0 0 0 p1 0.5 0 0 p2 0 60 0 p3 0 0 0.5 size 4 50 8 3 dim 0.045 0.045 0.045
0.045 rat 1 1 1 1 fill
gen zon radcylinder p0 0 0 0 p1 0.9 0 0 p2 0 60 0 p3 0 0 0.9 size 4 50 8 3 dim 0.5 0.5 0.5
0.5 rat 1 1 1 1
gen zone reflect normal (-1,0,0) origin 0 0 0
gen zone reflect normal (0,0,-1) origin 0 0 0
gen zone reflect normal (-1,0,0) origin -0.9 0 0
gen zone brick p0 -2.7 0 -0.9 p1 0.9 0 -0.9 p2 -2.7 60 -0.9 p3 -2.7 0 -1.4 size 6 30 4
gen zone brick p0 -2.7 0 0.9 p1 0.9 0 0.9 p2 -2.7 60 0.9 p3 -2.7 0 2.9 size 6 30 10
gen zone brick p0 -2.7 0 -1.4 p1 0.9 0 -1.4 p2 -2.7 60 -1.4 p3 -2.7 0 -2.4 size 6 30 5
gen zone brick p0 -2.7 0 2.9 p1 0.9 0 2.9 p2 -2.7 60 2.9 p3 -2.7 0 3.9 size 6 30 3
gen zon brick p0 -4.1 0 2.9 p1 -2.7 0 2.9 p2 -4.1 60 2.9 p3 -4.1 0 3.9 size 2 30 3
gen zon brick p0 -4.1 0 0.9 p1 -2.7 0 0.9 p2 -4.1 60 0.9 p3 -4.1 0 2.9 size 2 30 10
gen zon brick p0 -4.1 0 -0.9 p1 -2.7 0 -0.9 p2 -4.1 60 -0.9 p3 -4.1 0 0.9 size 2 30 8
gen zon brick p0 -4.1 0 -1.4 p1 -2.7 0 -1.4 p2 -4.1 60 -1.4 p3 -4.1 0 -0.9 size 2 30 4
gen zon brick p0 -4.1 0 -2.4 p1 -2.7 0 -2.4 p2 -4.1 60 -2.4 p3 -4.1 0 -1.4 size 2 30 5
gen zone reflect normal (1,0,0) origin 0.9 0 0.9
group zone ding range x -7.7 2.3 y 0 60 z 2.9 3.9
group zone shang range x -7.7 2.3 y 0 60 z 0.9 2.9
group zone jia range x -7.7 2.3 y 0 60 z 0.9 -0.9
group zone 1 range cyl end1 0 0 0 end2 0 60 0 rad 0.045
group zone 2 range cyl end1 -1.8 0 0 end2 -1.8 50 0 rad 0.045
group zone 3 range cyl end1 -3.6 0 0 end2 -3.6 50 0 rad 0.045
group zone 4 range cyl end1 -5.4 0 0 end2 -5.4 50 0 rad 0.045
group zone xia range x -7.7 2.3 y 0 60 z -0.9 -1.4
group zone di range x -7.7 2.3 y 0 60 z -1.4 -2.4
model ela
prop bulk 2.23e9 shear 1.46e9 den 2450 range group ding
prop bulk 1.25e9 shear 0.72e9 den 1340 range group shang
prop bulk 2.57e9 shear 1.7e9 den 2300 range group jia
prop bulk 2.57e9 shear 1.7e9 den 2300 range group 1
prop bulk 2.57e9 shear 1.7e9 den 2300 range group 2
prop bulk 2.57e9 shear 1.7e9 den 2300 range group 3
prop bulk 2.57e9 shear 1.7e9 den 2300 range group 4
prop bulk 1.25e9 shear 0.72e9 den 1340 range group xia
prop bulk 2.85e9 shear 1.88e9 den 2450 range group di
fix x y z range z -2.41 -2.39
fix x range x -7.71 -7.69
fix x range x 2.29 2.31
fix y range y -.1 .1
fix y range y 59.9 60.1
set gravity 0 0 -9.81
apply szz -6e6 range z 3.89 3.91
ini sxx = -12e6
ini syy = -8e6
hist unbal
solve
gen zon cshell p0 0 0 0 p1 0.5 0 0 p2 0 60 0 p3 0 0 0.5 size 4 50 8 3 dim 0.045 0.045 0.045
0.045 rat 1 1 1 1 fill
gen zon radcylinder p0 0 0 0 p1 0.9 0 0 p2 0 60 0 p3 0 0 0.9 size 4 50 8 3 dim 0.5 0.5 0.5
0.5 rat 1 1 1 1
gen zone reflect normal (-1,0,0) origin 0 0 0
gen zone reflect normal (0,0,-1) origin 0 0 0
gen zone reflect normal (-1,0,0) origin -0.9 0 0
gen zone brick p0 -2.7 0 -0.9 p1 0.9 0 -0.9 p2 -2.7 60 -0.9 p3 -2.7 0 -1.4 size 6 30 4
gen zone brick p0 -2.7 0 0.9 p1 0.9 0 0.9 p2 -2.7 60 0.9 p3 -2.7 0 2.9 size 6 30 10
gen zone brick p0 -2.7 0 -1.4 p1 0.9 0 -1.4 p2 -2.7 60 -1.4 p3 -2.7 0 -2.4 size 6 30 5
gen zone brick p0 -2.7 0 2.9 p1 0.9 0 2.9 p2 -2.7 60 2.9 p3 -2.7 0 3.9 size 6 30 3
gen zon brick p0 -4.1 0 2.9 p1 -2.7 0 2.9 p2 -4.1 60 2.9 p3 -4.1 0 3.9 size 2 30 3
gen zon brick p0 -4.1 0 0.9 p1 -2.7 0 0.9 p2 -4.1 60 0.9 p3 -4.1 0 2.9 size 2 30 10
gen zon brick p0 -4.1 0 -0.9 p1 -2.7 0 -0.9 p2 -4.1 60 -0.9 p3 -4.1 0 0.9 size 2 30 8
gen zon brick p0 -4.1 0 -1.4 p1 -2.7 0 -1.4 p2 -4.1 60 -1.4 p3 -4.1 0 -0.9 size 2 30 4
gen zon brick p0 -4.1 0 -2.4 p1 -2.7 0 -2.4 p2 -4.1 60 -2.4 p3 -4.1 0 -1.4 size 2 30 5
gen zone reflect normal (1,0,0) origin 0.9 0 0.9
group zone ding range x -7.7 2.3 y 0 60 z 2.9 3.9
group zone shang range x -7.7 2.3 y 0 60 z 0.9 2.9
group zone jia range x -7.7 2.3 y 0 60 z 0.9 -0.9
group zone 1 range cyl end1 0 0 0 end2 0 60 0 rad 0.045
group zone 2 range cyl end1 -1.8 0 0 end2 -1.8 50 0 rad 0.045
group zone 3 range cyl end1 -3.6 0 0 end2 -3.6 50 0 rad 0.045
group zone 4 range cyl end1 -5.4 0 0 end2 -5.4 50 0 rad 0.045
group zone xia range x -7.7 2.3 y 0 60 z -0.9 -1.4
group zone di range x -7.7 2.3 y 0 60 z -1.4 -2.4
model ela
prop bulk 2.23e9 shear 1.46e9 den 2450 range group ding
prop bulk 1.25e9 shear 0.72e9 den 1340 range group shang
prop bulk 2.57e9 shear 1.7e9 den 2300 range group jia
prop bulk 2.57e9 shear 1.7e9 den 2300 range group 1
prop bulk 2.57e9 shear 1.7e9 den 2300 range group 2
prop bulk 2.57e9 shear 1.7e9 den 2300 range group 3
prop bulk 2.57e9 shear 1.7e9 den 2300 range group 4
prop bulk 1.25e9 shear 0.72e9 den 1340 range group xia
prop bulk 2.85e9 shear 1.88e9 den 2450 range group di
fix x y z range z -2.41 -2.39
fix x range x -7.71 -7.69
fix x range x 2.29 2.31
fix y range y -.1 .1
fix y range y 59.9 60.1
set gravity 0 0 -9.81
apply szz -6e6 range z 3.89 3.91
ini sxx = -12e6
ini syy = -8e6
hist unbal
solve
