求大神帮解决

linscent

各位大侠好：
      最近使用solid185号单元，建立了一结构模型，初始给定材料参数为各向同性的，然后开始加载，每一次加载后改变材料刚度矩阵，并设定其为各向异性的，给定其21个参数，最开始的几次加载都没有问题，但是后来出现错误信息如下：“The ANEL material 1585(某一单元编号)at temperature=0 is not positive definite。。。”， 但是我查看了输入的材料刚度矩阵并没有出现奇异，是正定的，而且特征值都是挺大的正数，求各位大侠帮帮小弟解决此问题，不胜感激！

PS：当我选择的迭代步长比较大时，比如5000次一循环（此处的步长设定其实就是为了减少或增加计算次数，与非线性之类的步长设定毫无关系），最后可以计算完毕，而不出现错误，但是当我把迭代步长降低到3000时，就出现了上述错误，一时间搞不懂为什么错误会跟我的划分步长有关系，感觉两者是风马牛不相及的啊！

再次感谢各位！

shifang157300

（1）跟步长有关系，非线性迭代（一般是Newton-Rapshon迭代）过程中，每个迭代步越小，线性方程组病态的可能性越小。设想：假如你把每个载荷步的迭代步设置成2个，那么显然，很难由N-R迭代求得正确的节点位移增量，也就无法获得应变及应力，进而无法正确计算内力。
（2）至于为何出现非正定，要看具体参数。建议使用UPFs的usermat.f实现。

linscent

shifang157300 发表于 2013-3-28 09:24
（1）跟步长有关系，非线性迭代（一般是Newton-Rapshon迭代）过程中，每个迭代步越小，线性方程组病态的可 ...

谢谢啊， 不过不大理解

linscent

问题已经解决，请版主删帖！

linscent

问题已经解决，是自己犯了一个低级错误，在进行矩阵特征值求解的时候只是输入了一半信息，另一半对称的没有输入，所以使得自己设定的停止条件没有满足，而ansys检测到了错误，给出提示！