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

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 +7 j/.R  
new %gSmOW2.c^  
config fluid Vj8-[ww!  
set fluid off dNqj|Vu  
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 Buq(L6P9r  
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 0jxXUWO  
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 ]-a
/)8  
9PG{>W$M  
model fl_iso []yIz1P=j  
pro perm 0.599e-13 poro 0.0559 range group soil2 I3sf
O
U  
pro perm 0.599e-15 poro 0.04 range group soil3 -*ZQ=nomN  
pro perm 0.599e-15 poro 0.04 range group soil1 0Ny +NE:6M  
ini pp 1.5e6 range group soil1 'IVC!uL,%  
ini pp 0 range group soil3 on5\rY<I:@  
ini pp 0 range group soil2 2\|sXC  
ini fdensity 1.52 2,Z@<  
ini fmod 0.026e9 X
>o*eN  
  #u\~AO?h  
fix pp range z -0.1 0.1 ;底部约束 )S`Yl;oL  
/!6 VP |  
model mohr (sH4T>  
def cs @<4U &  
;五煤底板 pa?AKj]  
b_1=E_1/(3*(1-2*p_1)) **AkpV)  
s_1=E_1/(2*(1+p_1)) | "M1+(k7  
;五煤 -y3[\zNe  
b_2=E_2/(3*(1-2*p_2)) a%6=sqxE  
s_2=E_2/(2*(1+p_2)) *J':U>p  
;三煤底板 3|)cT1ej  
b_3=E_3/(3*(1-2*p_3)) |\| v%`r2  
s_3=E_3/(2*(1+p_3)) |m*l/@1  
end $ER$|9)KD  
_ogN  
set E_1=1e9 p_1=0.36 m@L>6;*  
set E_2=1.2e9 p_2=0.37 zA$ f$J7\^  
set E_3=1.3e9 p_3=0.36 :iQJ
9Hdz  
cs Q4S:/"*v8  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 /Zx"BSu
 
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 [{&OcEf  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 "^ dMCS@  
initial density = 2100 range group soil1 M-+!z5q~d  
initial density = 1800 range group soil2 V@gG x  
initial density = 2210 range group soil3 !Zc#E,  
d:cOdm>,  
fix x range x -.1 .1 A5%Now;.cf  
fix x range x 9.9 10.1 Zm6jF  
fix y range y -0.1 0.1 NQiu>Sg  
fix y range y 9.9 10.1 U!BZsVx  
fix z range z -.1 .1 J5Q.v;  
set gravity 0 0 -10 Zk|PQfi+  
step 10 |A0U3$S=  
save wei.sav Hxc>?  
ini xdis=0 ydis=0 zdis=0 qzZ;{>_f  
;根据体积应变计算渗透率 u'T>Y1I  
config zextra 2  : ?Z9  
def xxx h{jm  
  whilestepping o{hZjn-  
  pz=zone_head |GPYbxzc  
  loop while pz # null i`}9VaUG  
      b=z_vsi(pz) OP_\V8=  
      if z_group(pz)='soil2' 0dD.xuor  
        a=0.0599 @|OGxQo
C  
        k=0.599e-13 'fzJw  
        z_prop(pz,'poro')=(a+b)/(1+b) dww4o~hO  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 5;XU6Rz!  
      else ds{)p<LpT  
        a=0.04 Y10  
        k=0.599e-15 C&EA@U5X^  
        z_prop(pz,'poro')=(a+b)/(1+b) N!\1O,  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 rV-Xsf7Z  
      endif Aaz:C5dtU  
      z_extra(pz,1)=z_prop(pz,'poro') )WvOa] :  
      z_extra(pz,2)=z_prop(pz,'per') !SLfAFcS  
      pz=z_next(pz) 4>B=k  
  end_loop *W 04$N  
end DwTVoCC  
xxx M3JV^{O/DV  
;开挖 `:0Auw9h  
model null range group soil2 x 6 14 y 6 14 z 4 4.5 ;vk>k0S  
      2q9$5   
save jiang1.sav n%QWs1 b  
   X@A8~kj1  
  oXVx9dZ  
  |a'$v4dCF