model restore '8-1.sav'
fish define E_modify
p_z=zone.head
d_k=704
d_n=0.38
d_pa=101325.0
loop while p_z # null
sigma_3=-1.0*zone.stress.max(p_z)
E_new=d_k*d_pa*(sigma_3/d_pa)^d_n
zone.prop(p_z,'young')=E_new
p_z=zone.next(p_z)
endloop
end
@E_modify
model save '8-2.sav'
出错原因:Fish: Cannot raise negative number to real power.
自己想的是:幂函数的底数不能为负数,
但是为什么会出现幂函数为负数的情况,怎么解释?求助!!!