我用的是flac2d5.0版本的,软件提示渗透系数为0,没有成功赋值,貌似是wk11(i,p_z)这个表述有问题,但是查fish手册又只找到wk11,没有更详细的说明,查不到正确的表述方式...
希望有大神可以帮帮我
def k_modify
p_z=zone_head
d_kv0=1e-9
d_pa=101325.0
d_A=6
d_B=6
loopwhile p_z#null
sigma_1=-1.0*z_sig3(p_z)
d_wk11=d_kv0*(1+sigma_1/d_pa)^(-d_B)
d_wk22=d_A*wk11
wk11(i,p_z)=d_wk11
wk22(i,p_z)=d_wk22
p_z=z_next(p_z)
endloop
endk_modify