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

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 x~,?Zj)n?C  
new 4>eg@sN  
config fluid @)B5^[4(;  
set fluid off !Dn1pjxc  
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 24k}~"We  
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 #Yb9w3N  
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 ?fwr:aP~  
ep1Ajz.l  
model fl_iso WED7]2>  
pro perm 0.599e-13 poro 0.0559 range group soil2 =7Gi4X%  
pro perm 0.599e-15 poro 0.04 range group soil3 fC:\Gh5  
pro perm 0.599e-15 poro 0.04 range group soil1 B ~bU7.Cd  
ini pp 1.5e6 range group soil1  ]@ 0V  
ini pp 0 range group soil3 Z?[J_[ZtR3  
ini pp 0 range group soil2 d`*vJ#$>2  
ini fdensity 1.52 %ieAY-<"  
ini fmod 0.026e9 \}
(-9dr  
  8#3cmpx4  
fix pp range z -0.1 0.1 ;底部约束 (.t:sn"P  
\shoLp   
model mohr Eu_0n6J  
def cs c=mFYsSv  
;五煤底板 nr%^:u  
b_1=E_1/(3*(1-2*p_1)) +n]Knfi  
s_1=E_1/(2*(1+p_1)) #!<s& f|O  
;五煤 PltPIu)F  
b_2=E_2/(3*(1-2*p_2)) Xc<9[@  
s_2=E_2/(2*(1+p_2)) n m$G4Q  
;三煤底板 ?AH<y/i<Y  
b_3=E_3/(3*(1-2*p_3)) Yhdt8[ 2  
s_3=E_3/(2*(1+p_3)) SaScP  
end Wsz9X;  
ZACn_gd[5  
set E_1=1e9 p_1=0.36 .j,&/y&  
set E_2=1.2e9 p_2=0.37 F=lj$?4{  
set E_3=1.3e9 p_3=0.36 "A$Y)j<#G  
cs X*;p;N  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 Fsdn2{g8U  
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 2^j9m}`  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 Lg53 Ms%  
initial density = 2100 range group soil1 `OK }q  
initial density = 1800 range group soil2 1&WFs6  
initial density = 2210 range group soil3 5rhdm?Ls0  
&Vm[5X
W  
fix x range x -.1 .1 =swcmab;  
fix x range x 9.9 10.1 CF4y$aC#  
fix y range y -0.1 0.1 mQ)l`
wGh  
fix y range y 9.9 10.1 &0blHDMj{#  
fix z range z -.1 .1 vdM\scO:  
set gravity 0 0 -10 {uG_)GFr0  
step 10 g9Ll>d)tE3  
save wei.sav Fs~-exY1  
ini xdis=0 ydis=0 zdis=0 &}?e:PEy  
;根据体积应变计算渗透率 YRXe j  
config zextra 2 z[Qv}pv  
def xxx q3:' 69  
  whilestepping Mqy`j9FbL  
  pz=zone_head !pXz-hxKT  
  loop while pz # null b s*Z{R  
      b=z_vsi(pz) qdZYaS ~  
      if z_group(pz)='soil2' [!$>:_Vq/  
        a=0.0599 \/GY0s  
        k=0.599e-13 [O'p&j@  
        z_prop(pz,'poro')=(a+b)/(1+b) 1 ErYob.p  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3  )BB a  
      else >TJ$Z3
 
        a=0.04 1+#8} z:  
        k=0.599e-15 F3'G9Xf8Q=  
        z_prop(pz,'poro')=(a+b)/(1+b) NoIdO/vy"  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 YEj U3^@  
      endif -"H9W:  
      z_extra(pz,1)=z_prop(pz,'poro') b /ySt<  
      z_extra(pz,2)=z_prop(pz,'per') 2y,wN"qH*  
      pz=z_next(pz) .42OSV  
  end_loop Sy55w={  
end YG|T;/-  
xxx $<#sCrNX  
;开挖 `w}"
0+V  
model null range group soil2 x 6 14 y 6 14 z 4 4.5 >?V->7QLP  
     D.}b<kDD  
save jiang1.sav t~M_NEPxV  
   u,q#-d0g;  
  }mK_d9dx  
  D%btlw?{