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

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 fqp7a1qQl  
new cqU/Y_%l'  
config fluid Dqo:X`<bT  
set fluid off 9?mOLDu}Q0  
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 YUE1 '}  
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 Q,LWZw~"  
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 ,;-cz-,  
l9Cy30O6  
model fl_iso uC$4TnoQx.  
pro perm 0.599e-13 poro 0.0559 range group soil2 :r|P?;
t(  
pro perm 0.599e-15 poro 0.04 range group soil3 W SeRV?+T  
pro perm 0.599e-15 poro 0.04 range group soil1 @eD~FNf-]  
ini pp 1.5e6 range group soil1 [}g5Z=l  
ini pp 0 range group soil3 -T="Ml&  
ini pp 0 range group soil2 #JT%]!  
ini fdensity 1.52 &T4Cn@  
ini fmod 0.026e9
'2B0D|r"a  
 
uX5B>32  
fix pp range z -0.1 0.1 ;底部约束 K {N;k-  
FU3IK3}  
model mohr ,Y/>*,J  
def cs jC }u>AB  
;五煤底板 <PA$hTYM  
b_1=E_1/(3*(1-2*p_1)) O8}s*}]  
s_1=E_1/(2*(1+p_1)) Y&Nv>o_}5  
;五煤 Bf}0'MK8zQ  
b_2=E_2/(3*(1-2*p_2)) hFF&(t2{^  
s_2=E_2/(2*(1+p_2)) *g_>eNpXD  
;三煤底板 ~d
1RD  
b_3=E_3/(3*(1-2*p_3)) gQzF C&g  
s_3=E_3/(2*(1+p_3)) bsS|!KT  
end vf'jz
`Z  
^K@GK  
set E_1=1e9 p_1=0.36 (6_/n&mF  
set E_2=1.2e9 p_2=0.37 p<.!::*%(  
set E_3=1.3e9 p_3=0.36 FH@e:-*=  
cs k/f_@
8  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 l~>rpG  
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 oFA$X Y  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 wB0WR  
initial density = 2100 range group soil1 iO#xIl<  
initial density = 1800 range group soil2 ,kuFTWB  
initial density = 2210 range group soil3 J'oz P^N  
m3ZOq B-  
fix x range x -.1 .1 Z|7I }i  
fix x range x 9.9 10.1 PzG:M7  
fix y range y -0.1 0.1 %-zH]"Q$  
fix y range y 9.9 10.1 <L[)P{jn?p  
fix z range z -.1 .1 6I'VXdeN  
set gravity 0 0 -10 2Uw}'J_N  
step 10 NxRiEe#m  
save wei.sav 1JY90l$ME  
ini xdis=0 ydis=0 zdis=0 C]@B~X1H^  
;根据体积应变计算渗透率 [+2^n7R  
config zextra 2 = ~R3*GN  
def xxx A=>6$L];'  
  whilestepping p5 PON0dS  
  pz=zone_head lJpv
 
  loop while pz # null %JU23c*  
      b=z_vsi(pz) A^jm<~  
      if z_group(pz)='soil2' nFOG=>c}  
        a=0.0599 [&59n,R`  
        k=0.599e-13 aQ3vG08L>  
        z_prop(pz,'poro')=(a+b)/(1+b) m=b+V#4i(  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 (W6\%H2u  
      else W;*vcbP  
        a=0.04 JC-yiORVr  
        k=0.599e-15 #k]0[;1os  
        z_prop(pz,'poro')=(a+b)/(1+b) -CNv=vj 3  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 {twf7.eY  
      endif {+59YO  
      z_extra(pz,1)=z_prop(pz,'poro') &~6O;}\  
      z_extra(pz,2)=z_prop(pz,'per') )4O`%9=M&  
      pz=z_next(pz) +2enz!z#k  
  end_loop ^,r;/c9A8  
end ^hTJp{  
xxx X%qR6mMfT7  
;开挖 {nWtNyJpS  
model null range group soil2 x 6 14 y 6 14 z 4 4.5   r3K:  
     `o-<,  
save jiang1.sav jWjK-q@Y  
   w9}IM149  
  W_lNvzag  
  AZP>\Dq