y0=-1.5
y5=0
def excav1
y0=y0+1.5
y1=y0+1.5
command
model null range group 'kaiwa' y @y0 @y1
model null range group 'erchen' y @y0 @y1
model null range group 'chuzhi' y @y0 @y1
end_command
command
solve
end_command
end
@excav1
def excav2
cut_i=2
loop cut (cut_i,60)
cut_cur=cut
step_name='excavation step'+string(cut)
ii=out(step_name)
y0=1.5*(cut-1)
y1=y0+1.5
y2=y0-1.5
y3=y1-1.5
h=y0-y5
command
model null range group 'kaiwa' y @y0 @y1
model null range group 'erchen' y @y0 @y1
model null range group 'chuzhi' y @y0 @y1
end_command
command
model e range group 'chuzhi' y @y2 @y3
pro density 2200 bulk 12.78e9 shear 9.58e9 range group 'chuzhi' y @y2 @y3
end_command
if h>=60 then
y5=y5+9
y4=y5-9
command
model e range group 'erchen' y @y4 @y5
pro density 2300 bulk 17.58e9 shear 13.125e9 range group 'erchen' y @y4 @y5
end_command
end_if
command
solve
end_command
end_loop
end
@excav2