kq=1e-7*S^0.5*((1-(1-S^(1/0.26))^0.26)^2,表示的渗透率
S是自变量,表示的是冻土相对饱和度,初始值:0.75
当瞬态设置为:range(0,0.1,1) range(1.2,0.5,6),运算没有问题,当运算时间变长:range(0,0.1,1) range(1.2,0.5,6) range(6.2,1,30),就会出现以下错误,
非线性求解器不收敛。
时间:4.672。
试图计算负值的非整数次幂。
函数: ^
计算变量失败。
Variable: comp1.kq, Defined as: (ks*(comp1.S^l))*(((comp1.S^(1/m))^m)^2)
|
|