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

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 sB
'Z9  
new wh)Ujgd  
config fluid Dvz}sQZ  
set fluid off z2DjYTm[~  
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 e!.r- v9  
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  SED_^  
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 a+#Aitd  
Mr u  
model fl_iso GK\`8xWE  
pro perm 0.599e-13 poro 0.0559 range group soil2
q:X&)f  
pro perm 0.599e-15 poro 0.04 range group soil3 t?6_^ 08  
pro perm 0.599e-15 poro 0.04 range group soil1 XX;MoE~MM  
ini pp 1.5e6 range group soil1 U5pg<xI  
ini pp 0 range group soil3 hB 36o9|9  
ini pp 0 range group soil2 j%@wQVxq  
ini fdensity 1.52 RY9h^q*  
ini fmod 0.026e9 `6BjNV  
  NAt; r  
fix pp range z -0.1 0.1 ;底部约束 A@qwD300Vo  
)EptyH  
model mohr "Z70 jkW[  
def cs dlB?/J<  
;五煤底板 7PG|e#  
b_1=E_1/(3*(1-2*p_1)) zv[pfD7a  
s_1=E_1/(2*(1+p_1)) TOvpv@?-  
;五煤 :3$-Qv X  
b_2=E_2/(3*(1-2*p_2)) +Oa+G.;)o4  
s_2=E_2/(2*(1+p_2)) Mw9;O6  
;三煤底板 >UDb:N[  
b_3=E_3/(3*(1-2*p_3)) t!JD]j>q
 
s_3=E_3/(2*(1+p_3)) ,8/Con|o  
end u{dI[?@  
hm+,o_+  
set E_1=1e9 p_1=0.36 PShluhY  
set E_2=1.2e9 p_2=0.37 o
A'LQ  
set E_3=1.3e9 p_3=0.36 mw,\try  
cs 0CFON2I  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 EXD Qr'"  
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 %L;;W,l$`)  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 xY,W[?3CY  
initial density = 2100 range group soil1 xKW`m  
initial density = 1800 range group soil2 hfQx$cv6  
initial density = 2210 range group soil3 uiEAi  
0^IHBN?9  
fix x range x -.1 .1 sn6:\X<[  
fix x range x 9.9 10.1 rW<sQ0  
fix y range y -0.1 0.1 Q@"mL  
fix y range y 9.9 10.1 E`aAPk_y  
fix z range z -.1 .1 pg:1AAhT[  
set gravity 0 0 -10 U#{^29ik=o  
step 10 k,UezuV  
save wei.sav B^C!UWN>%X  
ini xdis=0 ydis=0 zdis=0 X5`#da  
;根据体积应变计算渗透率 ?}D|]i34  
config zextra 2 IS9}@5`'  
def xxx +o7Np|Ou  
  whilestepping ;[}<xw3):  
  pz=zone_head tbr1mw'G  
  loop while pz # null WUnmUW[/  
      b=z_vsi(pz) Jfs$VGZP;  
      if z_group(pz)='soil2' 9=wt9` ?  
        a=0.0599 EX.`6,:+2  
        k=0.599e-13 ]J aV +b'O  
        z_prop(pz,'poro')=(a+b)/(1+b) a&*fk?o  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3
wf[B-2q)  
      else @=kDaPme92  
        a=0.04 UNB'Xjp}@  
        k=0.599e-15 z%cpV{Nu  
        z_prop(pz,'poro')=(a+b)/(1+b) }VUrn2@-4  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 b9(_bsc  
      endif B-`,h pp  
      z_extra(pz,1)=z_prop(pz,'poro') A^9RGz4=  
      z_extra(pz,2)=z_prop(pz,'per') yS)73s/MrY  
      pz=z_next(pz) @! gJOy  
  end_loop zYW+Goz/C  
end B9glPcy}SS  
xxx
1dgy-$H~  
;开挖 (4WAoye|  
model null range group soil2 x 6 14 y 6 14 z 4 4.5 G?t<4MTv  
     snW=9b)m  
save jiang1.sav
(6*CORE   
   5[py{Gq  
  uN)o|7  
  bZB7t`C5