编写的fish 程序有问题,该怎么解决
下面的程序运行后出现“cannot raise neg. number to real power”这是因为什么参数设置错了吗?我理解,请高手指点new
restore 初始孔压.sav
;位移和速度归零
ini xdisp 0
ini ydisp 0
ini zdisp 0
ini xvel 0
ini yvel 0
ini zvel 0
config fl
;...fluid flow model...
model fl_iso
prop poros 0.5range group 1 any group 2 any
model fl_null range group 3
ini fden 1000
ini fmod 6.7e3 ftens -1e-3
water density 1000
def ana_sol
d_density=1450
w_density=1000
oumiga=0.15
aerfa=0.006
c_a=1/aerfa
c_n=1.45
c_m=1-1/c_n
sita_r=0.15
sita_s=0.5
sita=oumiga*d_density/w_density
sita1=sita-sita_r
sita2=sita_s-sita_r
h_xili=c_a*[(sita1/sita2)^(-1/c_m)-1]^(1-c_m)
perm_s=1e-11
k_perm=perm_s*{1-(aerfa*h_xili)^(c_n-1)*^(-c_m)}^2/^(c_m/2)
pnt=gp_head
loop while pnt # null
command
prop perm k_perm group 1 any group 2 any
endcommand
pnt=gp_next(pnt)
endloop
end
ana_sol 幂指数函数底数不能为负数,给底数加个绝对值abs函数即可解决
页:
[1]