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

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 \ "$$c  
new )<:TpMdUk  
config fluid 5<R m{  
set fluid off \m GY'0  
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 [!-gb+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 vIbM@Y4 '?  
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 y,3ZdY"  
,3y9yJQa*#  
model fl_iso qS FtQ4  
pro perm 0.599e-13 poro 0.0559 range group soil2 ,<r&] eC  
pro perm 0.599e-15 poro 0.04 range group soil3  7-!n-  
pro perm 0.599e-15 poro 0.04 range group soil1 UNff&E-  
ini pp 1.5e6 range group soil1 DQm%=ON7  
ini pp 0 range group soil3 |=m.eU  
ini pp 0 range group soil2 e)g&q'O  
ini fdensity 1.52 }Mt1C~{(  
ini fmod 0.026e9 So'.QWzX  
  7K:V<vX5  
fix pp range z -0.1 0.1 ;底部约束 =4a:)g'  
mpMAhm
:  
model mohr +8T^
q,  
def cs %kjG[C  
;五煤底板 v|o{AL:ei  
b_1=E_1/(3*(1-2*p_1)) !W9:)5^X  
s_1=E_1/(2*(1+p_1)) d%"XsbO  
;五煤 uMS+,dXy  
b_2=E_2/(3*(1-2*p_2)) L
zNfMvh  
s_2=E_2/(2*(1+p_2)) u0 tlf  
;三煤底板 ;ryNfP%  
b_3=E_3/(3*(1-2*p_3)) gJ'pwSA  
s_3=E_3/(2*(1+p_3)) !NkCki"W  
end &Xqxuy ]J  
5$D"uAp<V  
set E_1=1e9 p_1=0.36 mV$ebFco0  
set E_2=1.2e9 p_2=0.37 ElV!C}g  
set E_3=1.3e9 p_3=0.36 4n@lrcq(  
cs 5;UIz@BJ  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 rUj\F9*5#  
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 -6HwGfU  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 ]b!n ;{5  
initial density = 2100 range group soil1 f,KB BBbG  
initial density = 1800 range group soil2 JHt U"  
initial density = 2210 range group soil3 cN8Fn4g
q  
y~@zfJ5/^  
fix x range x -.1 .1 'in%Gii  
fix x range x 9.9 10.1 Kbf(P95+uL  
fix y range y -0.1 0.1 %BP>,E/w  
fix y range y 9.9 10.1 AXW.`~ 4  
fix z range z -.1 .1 k[;)/LfhS  
set gravity 0 0 -10 pB
8D  
step 10 <\u3p3"[4  
save wei.sav Y}N\|*ye-  
ini xdis=0 ydis=0 zdis=0 bYnq,JRA  
;根据体积应变计算渗透率 "4)N]Nj  
config zextra 2 $2?AJ/2r$b  
def xxx L2ydyXIsd  
  whilestepping 0!_?\)X  
  pz=zone_head _y_}/  
  loop while pz # null <ak[`]  
      b=z_vsi(pz) wG3b{0  
      if z_group(pz)='soil2' q!e
E~O;A  
        a=0.0599 YAOfuas]j  
        k=0.599e-13 D3X4@sM  
        z_prop(pz,'poro')=(a+b)/(1+b) _>
Raw  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 +*:mKx@Nw  
      else 49$<:{~  
        a=0.04 
.4)oZ  
        k=0.599e-15 7upko9d/  
        z_prop(pz,'poro')=(a+b)/(1+b) !S#3mT-  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 MK=oGzK  
      endif 7|DG1p9C  
      z_extra(pz,1)=z_prop(pz,'poro') 0lg$zi x(  
      z_extra(pz,2)=z_prop(pz,'per') v{VF>qEP  
      pz=z_next(pz) -}|L<~  
  end_loop u;-&r'J>  
end %+N]$Q  
xxx O {1" I  
;开挖 Pc`d]*BYi  
model null range group soil2 x 6 14 y 6 14 z 4 4.5 Cp6S2v I  
     >|E]
??v  
save jiang1.sav T8x)i\<  
   5M0Q'"`F:  
  ir_XU/ve  
  |Ev|A9J!