跪求土木大佬帮看看这个弹塑性本构代码的问题

望舒

之前编写有限元遇到过这样的问题，
NAN的输出结果极有可能是你没有给矩阵（二维数组）所有位置都 赋值。（哪怕是0也得赋值）
你如果不赋值系统就会分配随机数进去，出现大量很大的数字或者NAN。
所以在给刚度矩阵、应力张量、弹性常数等矩阵赋值的时候务必注意不要让系统随机。