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

大家帮我看看为什么会出现'property 'poro' not found in zone 1,在这里谢谢了。 K#bdb  
new RAMkTS  
config fluid nR)/
k,3W  
set fluid off Ed[ tmaEuV  
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 cIqk=_]  
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 Wyu$J  
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 PI{;3X}9$,  
6<sd6SM  
model fl_iso utIR\e#:B  
pro perm 0.599e-13 poro 0.0559 range group soil2 ConXP\M-  
pro perm 0.599e-15 poro 0.04 range group soil3 f49pIcAq  
pro perm 0.599e-15 poro 0.04 range group soil1 dQ6:c7hp>D  
ini pp 1.5e6 range group soil1 C/Vs+aW n  
ini pp 0 range group soil3 2e_ssBbb  
ini pp 0 range group soil2 61|uvTX  
ini fdensity 1.52 5W/!o&x~7  
ini fmod 0.026e9 hE>ux"_2/  
  }{>)2S  
fix pp range z -0.1 0.1 ;底部约束 +_Nr a  
<a/TDW  
model mohr <<#-IsT  
def cs `"CIy_m  
;五煤底板 i#/,Q1yEn  
b_1=E_1/(3*(1-2*p_1)) 4]L5%=atn  
s_1=E_1/(2*(1+p_1)) ACQc 0:q  
;五煤 Rh!B4oB4  
b_2=E_2/(3*(1-2*p_2)) r*#ApM"L  
s_2=E_2/(2*(1+p_2)) EBjSK
/  
;三煤底板 eQh@.U*S)  
b_3=E_3/(3*(1-2*p_3)) /Qh  
s_3=E_3/(2*(1+p_3)) Oax*3TD  
end 7_Yxz$m  
nY(>|!  
set E_1=1e9 p_1=0.36 wQR>S>p  
set E_2=1.2e9 p_2=0.37 .FUws  
set E_3=1.3e9 p_3=0.36 `#O%ZZ+  
cs @tQu3Rq@  
property bulk = b_1 shear = s_1 friction = 25 cohesion = 0.2e6 ten = 0.04e6 range group soil1 nx#0*r}5  
property bulk = b_2 shear = s_2 friction = 23 cohesion = 0.1e6 ten = 0.04e6 range group soil2 \Eq,4-q  
property bulk = b_3 shear = s_3 friction = 21 cohesion = 0.26e6 ten = 0.06e6 range group soil3 ~{Gbu
oH  
initial density = 2100 range group soil1 4.7YIM  
initial density = 1800 range group soil2 Dqss/vw
V  
initial density = 2210 range group soil3 on?<3eED  
zrt\]h+  
fix x range x -.1 .1 N ]
KS\  
fix x range x 9.9 10.1 Dep.Qfv{-  
fix y range y -0.1 0.1 NfG<!  
fix y range y 9.9 10.1 !cP2,l'f  
fix z range z -.1 .1 cj ?aCVa  
set gravity 0 0 -10 r)Q/YzXx*  
step 10 oNW.-gNT  
save wei.sav bI~ R6o  
ini xdis=0 ydis=0 zdis=0 }^|g|xl!  
;根据体积应变计算渗透率 0=N,y  
config zextra 2 Lhg4fuos@)  
def xxx RmF,x9  
  whilestepping ~s.~X5  
  pz=zone_head !*2cK>`  
  loop while pz # null wOW#A}m'vj  
      b=z_vsi(pz) ^L}ICm_#  
      if z_group(pz)='soil2' Z) nB  
        a=0.0599 ITcgpK6k  
        k=0.599e-13 ;U7o)A;  
        z_prop(pz,'poro')=(a+b)/(1+b) ze* =7  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 \o-9~C\c*  
      else HEqTlnxUu  
        a=0.04 % zP]z  
        k=0.599e-15 7Mx F? I  
        z_prop(pz,'poro')=(a+b)/(1+b) bv %Bo4s  
        z_prop(pz,'perm')=(k/(1+b))*(1+b/a)^3 xJCMxt2Y  
      endif W 7xh  
      z_extra(pz,1)=z_prop(pz,'poro') %6A."sePO  
      z_extra(pz,2)=z_prop(pz,'per') =;9
%Q{  
      pz=z_next(pz) ug?gVK  
  end_loop :bBLP7eyV  
end Ew`(x30E  
xxx @9,=|kxK  
;开挖 D~~
"wos  
model null range group soil2 x 6 14 y 6 14 z 4 4.5 Ck`-<)uN  
     'gBns  
save jiang1.sav >"LHr&;m&h  
   6I)[6R  
  12{F  
  Sq#AnD6To