关于z_prop(pz,'poro')的错误 [复制链接]

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 E\$W_Lmr  
new eDB;cN  
config fluid i6N',&jFU  
set fluid off {>;R?TG]$  
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 QSj]ZA  
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 9!tW.pK5  
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 :Qq#Z  
P|> ~_$W  
model fl_iso A)KZa
"EX  
pro perm 0.599e-13 poro 0.0559 range group soil2 .p$(ZH =~  
pro perm 0.599e-15 poro 0.04 range group soil3 S[QrS7  
pro perm 0.599e-15 poro 0.04 range group soil1 I2DpRMy  
ini pp 1.5e6 range group soil1 mBC+6(5V  
ini pp 0 range group soil3 'f|o{  
ini pp 0 range group soil2 _XT pU  
ini fdensity 1.52 B1Oq!k  
ini fmod 0.026e9 Q=:|R
3U/  
  hzC>~Ub5  
fix pp range z -0.1 0.1 ;底部约束 r_.S>]  
^}C\zW  
model mohr jqkqZF  
def cs O9p|a%o  
;五煤底板 "ITIhnE  
b_1=E_1/(3*(1-2*p_1)) lRdChoL$2  
s_1=E_1/(2*(1+p_1)) 6zn5UW#q  
;五煤 4/)k)gLI  
b_2=E_2/(3*(1-2*p_2)) F<w/PMb  
s_2=E_2/(2*(1+p_2)) b>$S<td  
;三煤底板 1nOCQ\$l  
b_3=E_3/(3*(1-2*p_3)) /Q )\+  
s_3=E_3/(2*(1+p_3)) s(8W_4&'  
end ,2)6s\]/b  
!VK
|u8i
 
set E_1=1e9 p_1=0.36 )_NO4`ejs/  
set E_2=1.2e9 p_2=0.37 *owU)  
set E_3=1.3e9 p_3=0.36 yppo6HGD  
cs Xm2'6f,  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 Cls%M5
MH  
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 kNL\m[W8$  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 {y;n:^  
initial density = 2100 range group soil1 4`R(?  
initial density = 1800 range group soil2 W'.m'3#z  
initial density = 2210 range group soil3 w*MpX U<  
V1JIht>Opo  
fix x range x -.1 .1 PxE3K-S)G  
fix x range x 9.9 10.1 \|ao`MMaD<  
fix y range y -0.1 0.1 9k=3u;$v  
fix y range y 9.9 10.1 v9UD%@tZ  
fix z range z -.1 .1 D sWSGb  
set gravity 0 0 -10 D,ln
)["xm  
step 10 FCn_^l)EA  
save wei.sav fxIf|9Qi`  
ini xdis=0 ydis=0 zdis=0 ,?XCyHSgWW  
;根据体积应变计算渗透率 bYPKh  
config zextra 2 YAmb`CP  
def xxx >"<Wjr8W!$  
  whilestepping &t-kpA|EG  
  pz=zone_head _-Fs#f8  
  loop while pz # null VD\=`r)nT  
      b=z_vsi(pz) +'w3 =2Bo  
      if z_group(pz)='soil2' xP,hTE  
        a=0.0599 d`6 'Z  
        k=0.599e-13  rXU\  
        z_prop(pz,'poro')=(a+b)/(1+b) K_-MYs.  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 j8`BdKg  
      else 5 u0HI  
        a=0.04 E+JqWR5  
        k=0.599e-15 :/Qq@]O>  
        z_prop(pz,'poro')=(a+b)/(1+b) 4a]P7fx-  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 40<mrVl  
      endif y^%y<~f
 
      z_extra(pz,1)=z_prop(pz,'poro') >dG[G>  
      z_extra(pz,2)=z_prop(pz,'per') C>w|a  
      pz=z_next(pz) &8 x-o,  
  end_loop BVO<e \>3  
end K96<M);:g  
xxx r>U@3%0&  
;开挖 O8.5}>gDn.  
model null range group soil2 x 6 14 y 6 14 z 4 4.5 #1G:lhkC  
     IN G@B#Cl  
save jiang1.sav ?3xzd P  
   N8FF3}> g  
  ;d9QAN&0}  
  D5HZ2cz|a