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)

我在设置变量kq时，就算添加了if(S>0,Kq,1）但是运算还是没有用，而且我的初始值设置为0.75，S本身的范围就应该在0~1之间，怎么会出现负值呢？

我在计算的时候也会遇见这种情况，我通过调整边界条件和初值解决的。