z_prop出现的错误

jwj15tty

帮我看看为什么会出现property 'poro' not found in zone 1,我在上面都付了值啊。求高手指点指点
new
config fluid
set fluid off
gen zone brick p0 0 0 0 p1 20 0 0 p2 0 20 0 p3 0 0 4
size 20 20 4 group soil1
gen zone brick p0 0 0 4 p1 20 0 4 p2 0 20 4 p3 0 0 6
size 20 20 4 group soil2
gen zone brick p0 0 0 6 p1 20 0 6 p2 0 20 6 p3 0 0 10
size 20 20 4 group soil3

model fl_iso
pro perm 0.599e-13 poro 0.0559 range group soil2
pro perm 0.599e-15 poro 0.04 range group soil3
pro perm 0.599e-15 poro 0.04 range group soil1
ini pp 1.5e6 range group soil1
ini pp 0 range group soil3
ini pp 0 range group soil2
ini fdensity 1.52
ini fmod 0.026e9
  
fix pp range z -0.1 0.1 ;底部约束

model mohr
def cs

b_1=E_1/(3*(1-2*p_1))
s_1=E_1/(2*(1+p_1))

b_2=E_2/(3*(1-2*p_2))
s_2=E_2/(2*(1+p_2))

b_3=E_3/(3*(1-2*p_3))
s_3=E_3/(2*(1+p_3))

end

set E_1=1e9 p_1=0.36
set E_2=1.2e9 p_2=0.37
set E_3=1.3e9 p_3=0.36

cs
property bulk = b_1 shear = s_1 friction = 25 cohesion
= 0.2e6 ten = 0.04e6 range group soil1
property bulk = b_2 shear = s_2 friction = 23 cohesion
= 0.1e6 ten = 0.04e6 range group soil2
property bulk = b_3 shear = s_3 friction = 21 cohesion
= 0.26e6 ten = 0.06e6 range group soil3
initial density = 2100 range group soil1
initial density = 1800 range group soil2
initial density = 2210 range group soil3

fix x range x -.1 .1
fix x range x 9.9 10.1
fix y range y -0.1 0.1
fix y range y 9.9 10.1
fix z range z -.1 .1
set gravity 0 0 -10
step 10
save wei.sav
ini xdis=0 ydis=0 zdis=0
config zextra 2
def xxx
  whilestepping
  pz=zone_head
  loop while pz # null
      b=z_vsi(pz)
      if z_group(pz)='soil2'
        a=0.0599
        k=0.599e-13
        z_prop(pz,'poro')=(a+b)/(1+b)
        z_prop(pz,'per')=(k/(1+b))*(1+b/a)^3
      else
        a=0.04
        k=0.599e-15
        z_prop(pz,'poro')=(a+b)/(1+b)
        z_prop(pz,'per')=(k/(1+b))*(1+b/a)^3
      endif
      z_extra(pz,1)=z_prop(pz,'poro')
      z_extra(pz,2)=z_prop(pz,'per')
      pz=z_next(pz)
  end_loop
end
xxx
;开挖
model null range group soil2 x 6 14 y 6 14 z 4 4.5
solve     
save jiang1.sav

doctorluo

z_prop(pz,'poro')改为z_prop(pz,'porosity')
z_prop(pz,'per')改为 z_prop(pz,'permeablity‘)