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

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 ZdlZ,vK^.  
new hrNri$  
config fluid TWQf2  
set fluid off Z%N{Y x(  
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 un6grvxr  
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 z~A||@4'  
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 SvAz9>N4  
]3NH[&+  
model fl_iso G! zV=p  
pro perm 0.599e-13 poro 0.0559 range group soil2 VDx=Tsu-  
pro perm 0.599e-15 poro 0.04 range group soil3 Dvo.yn|kB  
pro perm 0.599e-15 poro 0.04 range group soil1 -G(z!ed  
ini pp 1.5e6 range group soil1 3 {\b/NL$  
ini pp 0 range group soil3 ul?'kuYk  
ini pp 0 range group soil2  |UZ#2  
ini fdensity 1.52 {uL<$;#i  
ini fmod 0.026e9 E&#AX:  
  H"WkyvqXb  
fix pp range z -0.1 0.1 ;底部约束 -F/"W  
9-Ikd>9  
model mohr ^p|MkB?uM  
def cs f%an<>j^w  
;五煤底板 G>M# BuU  
b_1=E_1/(3*(1-2*p_1)) u!k\W{  
s_1=E_1/(2*(1+p_1)) 7F>5<Gv:-  
;五煤 3:PBVt=  
b_2=E_2/(3*(1-2*p_2)) c;&m}ImLe.  
s_2=E_2/(2*(1+p_2)) bDeHU$  
;三煤底板 ?>1AT==wI  
b_3=E_3/(3*(1-2*p_3)) mSY;hJi  
s_3=E_3/(2*(1+p_3)) a^~T-;_V  
end KX"?3#U#Fm  
Lf3Ri/@ p  
set E_1=1e9 p_1=0.36 PM=
Q\0  
set E_2=1.2e9 p_2=0.37 k*J}/HO  
set E_3=1.3e9 p_3=0.36 (W
VN*OR?  
cs Cj _Q9/  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 2ELw}9  
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 ?uc=(J+6  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 Ih;6(5z  
initial density = 2100 range group soil1 <?yAIhgN*  
initial density = 1800 range group soil2 K,&)\r kzD  
initial density = 2210 range group soil3 Y1rU  
$LRvPan`  
fix x range x -.1 .1 8+k\0fmy  
fix x range x 9.9 10.1 Dq)V]
Zx  
fix y range y -0.1 0.1 54_CewL1P]  
fix y range y 9.9 10.1 R61.!ql%w  
fix z range z -.1 .1 JX7_/P  
set gravity 0 0 -10 ZK[4n5}  
step 10 Lp5LRw  
save wei.sav %Nwap~=H;  
ini xdis=0 ydis=0 zdis=0 x>A(016:C  
;根据体积应变计算渗透率 )-xx$0mL-  
config zextra 2 f>s3Q\+  
def xxx H1yl88K  
  whilestepping |g;hXr#~  
  pz=zone_head MRxzOs  
  loop while pz # null 4`o_r%   
      b=z_vsi(pz) w,SOvbAxX2  
      if z_group(pz)='soil2' U; JZN  
        a=0.0599 EoD;'+d  
        k=0.599e-13 ?[hIv6c  
        z_prop(pz,'poro')=(a+b)/(1+b) -8-  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3
R+Dx#Wn I  
      else Jt\?
,~,
  
        a=0.04 G#t!{Q}8  
        k=0.599e-15 <.HX_z3l  
        z_prop(pz,'poro')=(a+b)/(1+b) Z `sM(?m  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 ?x7zYE,6  
      endif gXZC%S  
      z_extra(pz,1)=z_prop(pz,'poro') fW~r%u .y  
      z_extra(pz,2)=z_prop(pz,'per') mY*JNx  
      pz=z_next(pz) 3,
DUT{2  
  end_loop a6wPkf7-H  
end y2;uG2IS_g  
xxx ~2@+#1[g8z  
;开挖 3.Z}2F]  
model null range group soil2 x 6 14 y 6 14 z 4 4.5 pUutI|mt/  
     qPDe;$J)  
save jiang1.sav m/Ou$  
   H:Y?("k  
  1v)ur\>R  
  %ZF47P%6