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

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 u^HC1r|%  
new 0wM2v[^YO  
config fluid {kCCpU  
set fluid off U($bR|%D  
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 Nm;y
L  
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 w; [ndZCY7  
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 >Q5et1c  
GRM:o)4;#  
model fl_iso %b8ig1  
pro perm 0.599e-13 poro 0.0559 range group soil2 "DN0|%`M/  
pro perm 0.599e-15 poro 0.04 range group soil3 U"nk AW  
pro perm 0.599e-15 poro 0.04 range group soil1 &8p]yo2zO  
ini pp 1.5e6 range group soil1 {)V!wSi  
ini pp 0 range group soil3 oT76)O  
ini pp 0 range group soil2 ~ _ ogeD  
ini fdensity 1.52 H(L.k;B  
ini fmod 0.026e9 in-|",O`Z  
  lZ5LHUzP  
fix pp range z -0.1 0.1 ;底部约束 Qt-7jmZw1  
PHQ7  
model mohr E}?n^Zf  
def cs WP5cC@x  
;五煤底板 x*H,eY3  
b_1=E_1/(3*(1-2*p_1)) S2y_5XJ<D  
s_1=E_1/(2*(1+p_1)) 'e&L53n  
;五煤 i}u,_ }  
b_2=E_2/(3*(1-2*p_2)) #mT\B[4h  
s_2=E_2/(2*(1+p_2)) -!o*A>N  
;三煤底板 @477|LO
  
b_3=E_3/(3*(1-2*p_3)) voX4A pl  
s_3=E_3/(2*(1+p_3)) eILdq*  
end 4S%s
=vw  
` nd/N#  
set E_1=1e9 p_1=0.36 ,> Ya%;h2k  
set E_2=1.2e9 p_2=0.37 `$PdI4~J  
set E_3=1.3e9 p_3=0.36 {ef9ov Xk  
cs o,r72>|  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 F
IMM\W  
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 >W?i+,g  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 91f{qq=#J{  
initial density = 2100 range group soil1 6!39t  
initial density = 1800 range group soil2 ^LI\W'K  
initial density = 2210 range group soil3 PHAM(iC&D  
\wZ 4enm  
fix x range x -.1 .1 |WkWZZ^  
fix x range x 9.9 10.1 scc+r  
fix y range y -0.1 0.1 PJfADB7Y  
fix y range y 9.9 10.1
_cc37[  
fix z range z -.1 .1 v(0IQ  
set gravity 0 0 -10 ^;
8dl.;  
step 10 Ez1-Nx  
save wei.sav @23?II$=@  
ini xdis=0 ydis=0 zdis=0 ahh&h1q7|  
;根据体积应变计算渗透率 L`2(u!i J  
config zextra 2 a!@(bb z>  
def xxx 2v$\mL  
  whilestepping &0b\E73  
  pz=zone_head } yb"/
jp  
  loop while pz # null 'G1~ A +  
      b=z_vsi(pz) wC>}9
OM  
      if z_group(pz)='soil2' @vc9L  
        a=0.0599 (X?HuWTm  
        k=0.599e-13
C1`fJhy  
        z_prop(pz,'poro')=(a+b)/(1+b) 7+NBcZuG9  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 w,T-vf  
      else R@df~  
        a=0.04 BT^Im=A  
        k=0.599e-15 Kj7 ?_o{  
        z_prop(pz,'poro')=(a+b)/(1+b) %!\iII  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 ([ jm=[E^  
      endif 9oj#5Hq  
      z_extra(pz,1)=z_prop(pz,'poro') <B Vx%  
      z_extra(pz,2)=z_prop(pz,'per') `_iK`^(-  
      pz=z_next(pz) u<U8LR=)V5  
  end_loop #Zg pm"MW  
end 1Y:JGon  
xxx iBPIj;,  
;开挖 |4C5;"Pc  
model null range group soil2 x 6 14 y 6 14 z 4 4.5 6vfut$)[{  
     Y1wH_!%b  
save jiang1.sav 1rIL[(r4  
   5(>ux@[qI:  
  79z/(T
+  
  Dnl|B\