本人模拟深基坑开挖，初始应力场塑性求解在不平衡力数量级e3过后一直不收敛，一直上下波动，按esc终止后，plot发现zdis 大的惊人，往后面开挖算了几步，地面沉降都好几米，找了很久找不到原因。是否地层参数的问题，求大神指教，命令如下：
def inisoil
array soil_z(20)
array exca_z(20)
array mp_z(20)
array soi_depth_seg(20)
array mp_depth_seg(20)
array exca_depth_seg(20)
global_group_id=1
end
;initialize the arrays needed
inisoil
def gen_out_ground
loop n (1,soil_num)
p0x=0
p0y=0
p0z=soil_z(n)
p1x=r_max
p1y=0
p1z=soil_z(n)
p2x=0
p2y=0
p2z=soil_z(n+1)
;to assert that p2z<p0z
if p2z>p0z
oo=out('p2z<p0z,gen zone error!')
endif
p3x=0
p3y=r_max
p3z=soil_z(n)
p4x=r_max
p4y=0
p4z=soil_z(n+1)
p5x=0
p5y=r_max
p5z=soil_z(n+1)
p6x=r_max
p6y=r_max
p6z=soil_z(n)
p7x=r_max
p7y=r_max
p7z=soil_z(n+1)
p8x=exca_x_max
p8y=0
p8z=soil_z(n)
;other parameters needed
d1=exca_x_max
d2=exca_y_max
d3=d1
d4=d2
n1=inner_seg
n3=arc_seg
n2=soi_depth_seg(n)
n4=out_seg
r1=1
r2=1
r3=1
r4=1.5
group_id = 'soil'+string(global_group_id)
command
gen zone radtunnel p0 p0x,p0y,p0z p1 p1x,p1y,p1z p2 p2x,p2y,p2z p3 p3x,p3y,p3z &
dim d1,d2,d3,d4 size n1 n2 n3 n4 ratio r1 r2 r3 r4 group group_id;fill
endcommand
global_group_id=global_group_id+1
endloop
end
def gen_exca_region
loop n (1,exca_num)
p0x=0
p0y=0
p0z=exca_z(n)
p1x=exca_x_max
p1y=0
p1z=exca_z(n)
p2x=0
p2y=0
p2z=exca_z(n+1)
p3x=0
p3y=exca_y_max
p3z=exca_z(n)
n1=exca_seg
n2=exca_depth_seg(n)
n3=exca_arc_seg
r1=1
r2=1
r3=1
group_id = 'exca'+string(global_group_id)
command
gen zone brick p0 p0x,p0y,p0z p1 p1x,p1y,p1z p2 p2x,p2y,p2z p3 p3x,p3y,p3z &
size n1 n2 n3 ratio r1 r2 r3 group group_id
endcommand
global_group_id=global_group_id+1
endloop
end
def gen_full_model
command
gen zone reflect ori 0,0,0 dip 0 dd 90 norm 0,-1,0
gen zone reflect ori 0,0,0 dip 0 dd 90 norm -1,0,0
endcommand
end
;initial the settings
set grav 0,0,-10
set mech force 500 ratio 1.0e-5
set log on
;***************************************************
;the following code gen the geometry of the ground
;***************************************************
call genzone.txt
def set_soil
soil_z(1)=0
soil_z(2)=-6
soil_z(3)=-8.9
soil_z(4)=-12.7
soil_z(5)=-14
soil_z(6)=-16.4
soil_z(7)=-17
soil_z(8)=-19.6
soil_z(9)=-20.6
soil_z(10)=-36.6 exca_z(1)=0
exca_z(2)=-2
exca_z(3)=-5
exca_z(4)=-12
exca_z(5)=-17
exca_z(6)=-23
exca_z(7)=-27.836
exca_z(8)=-36.6
soi_depth_seg(1)=2
soi_depth_seg(2)=1
soi_depth_seg(3)=1
soi_depth_seg(4)=1
soi_depth_seg(5)=1
soi_depth_seg(6)=1
soi_depth_seg(7)=1
soi_depth_seg(8)=1
soi_depth_seg(9)=3 exca_depth_seg(1)=1
exca_depth_seg(2)=1
exca_depth_seg(3)=1
exca_depth_seg(4)=1
exca_depth_seg(5)=1
exca_depth_seg(6)=1
exca_depth_seg(7)=3
end
set_soil
set soil_num 9
set r_max 100
set inner_seg 10
set out_seg 10
set arc_seg 10
set exca_x_max 10.3
set exca_y_max 16.6gen_out_ground
set exca_num 7
set exca_x_max 10.3
set exca_y_max 16.6
set exca_arc_seg 20
set exca_seg 20gen_exca_region
gen_full_model
;************************************************************
;set the materail of the soil
;************************************************************
def trans
p_ratio=K_mod/(1+K_mod)
E_mod=y_mod*3
s_mod = E_mod / (2.0 * (1.0 + p_ratio))
b_mod = E_mod / (3.0 * (1.0 - 2.0 * p_ratio))
end
;人工填筑土
mod mohr rang z -6 0
set y_mod 2.5e6
set K_mod 0.4
trans
pro bulk b_mod shear s_mod coh 10.0e3 fri 8.0 dilation 0.0 dens 1.9e3 rang z -6 0
;松散卵石
mod mohr rang z -8.9 -6
set y_mod 25e6
set K_mod 0.3
trans
pro bulk b_mod shear s_mod coh 0 fri 35.0 dilation 0.0 dens 1.9e3 rang z -8.9 -6 ;稍密卵石
mod mohr rang z -12.7 -8.9
set y_mod 30e6
set K_mod 0.3
trans
pro bulk b_mod shear s_mod coh 0 fri 40 dilation 0.0 dens 2e3 rang z -12.7 -8.9
;松散卵石 洪积层
mod mohr rang z -14 -12.7
set y_mod 25e6
set K_mod 0.30
trans
pro bulk b_mod shear s_mod coh 0 fri 25 dilation 0.0 dens 1.9e3 rang z -14 -12.7;稍密卵石洪积层
mod mohr rang z -16.4 -14
set y_mod 30e6
set K_mod 0.3
trans
pro bulk b_mod shear s_mod coh 0 fri 45 dilation 0.0 dens 2e3 rang z -16.4 -14
;细砂（松散）
mod mohr rang z -17 -16.4
set y_mod 6e6
set p_ratio 0.4
trans
pro bulk b_mod shear s_mod coh 0 fri 27 dilation 0.0 dens 1.8e3 rang z -17 -16.4
;稍密卵石洪积层
mod mohr rang z -19.6 -17
set y_mod 30e6
set K_mod 0.3
trans
pro bulk b_mod shear s_mod coh 0 fri 45 dilation 0.0 dens 2e3 rang z -19.6 -17 ;细砂（松散）
mod mohr rang z -20.6 -19.6
set y_mod 6e6
set p_ratio 0.4
trans
pro bulk b_mod shear s_mod coh 0 fri 27 dilation 0.0 dens 1.8e3 rang z -20.6 -19.6
;密实卵石洪积层
mod mohr rang z -36.6 -20.6
set y_mod 55e6
set p_ratio 0.2
trans
pro bulk b_mod shear s_mod coh 0 fri 52 dilation 0.0 dens 2.25e3 rang z -36.6 -20.6
;************************************************************
;set the fixty
;************************************************************
;fix the bottom
fix x y z rang z -36.61 -36.59
fix x rang x 100.1 99.9
fix x rang x -100.1 -99.9fix y rang y 100.1 99.9
fix y rang y -100.1 -99.9
hist id=1 unbal
solve elas
