Error：the solution has a non-positive definite stiffness matrix

fusizhide101

热力场耦合分析，计算刚开始不就就出现这个错误，求指教啊！！

Cob

I guess you use the implicit method. if you did use implicit method, you may want to specify BFGS as solution method. That may solve your problem, because BFGS method will form a positive definite matrix artifically.
And I think something may be wrong in your model too, try to double check your model carefully.

The information you provided is too limit, what I suggest just based on my guess. Hope it can help a little bit.

adriano198707

Cob 发表于 2011-5-10 01:15
I guess you use the implicit method. if you did use implicit method, you may want to specify BFGS as ...

你好，我想请问你一下，我在做一个正交切削的仿真模拟，目前力学模拟已经有结果了，按照某些论文所述添加热固耦合参数后提交计算发现计算到一半出现non-positive definite stiffness matrix in the first time step.
想请问一下是不是热固耦合分析的时间步长过大引起的？
例外，我试着将热固耦合的时间不长都设置为0（默认值），发现计算没多久温度输出文件tprint变得很大（几个G）于是就中断了计算。

这是我的原帖地址，我附上了K文件，请您抽空帮我看一下行吗？谢谢了
http://forum.simwe.com/forum.php ... D154%26typeid%3D154

adriano198707

Cob 发表于 2011-5-10 01:15
I guess you use the implicit method. if you did use implicit method, you may want to specify BFGS as ...

你好，我又试了一下，将材料的比热容提高了1000倍，计算时进行下去了，但是感觉温度的升高不够，才6°C。

能否告知一下做热固耦合时需要注意些什么参数么？
标准单位又是怎样的呢？

我用的是ISO单位，m-kg-s

Cob

adriano198707 发表于 2012-3-26 17:27
你好，我又试了一下，将材料的比热容提高了1000倍，计算时进行下去了，但是感觉温度的升高不够，才6°C。 ...

non-positive definite stiffness matrix 一般都是因为模型有些问题。用BFGS 求解法能部分缓解问题，但是不一定能彻底。
你把比热容提高了1000倍，据我猜测是以前的错误+错误的参数反而使计算稳定了。你还是在仔细查查模型。另外，Q=cmDeltaT, 你切削过程中转化的热能是不变的（Q），你把比热容提高了1000倍（C）, 温度自然升得少了

adriano198707

Cob 发表于 2012-3-28 22:46
non-positive definite stiffness matrix 一般都是因为模型有些问题。用BFGS 求解法能部分缓解问题，但是 ...

你好，感谢你百忙之中回复我的帖子。

我先用错误的材料参数进行模拟，现在计算还是没能完全完成，但是会出现负温度。
想请问一下模型一般来说会有些什么样的问题，我是新手，希望能得到指点。

谢谢了

adriano198707

Cob 发表于 2012-3-28 22:46
non-positive definite stiffness matrix 一般都是因为模型有些问题。用BFGS 求解法能部分缓解问题，但是 ...

另外，我在另外一个帖子上看到说当单元删除后隐式计算会出现错误，有这回事么？
原帖地址：http://forum.simwe.com/thread-865874-1-1.html

Cob

adriano198707 发表于 2012-3-28 23:01
你好，感谢你百忙之中回复我的帖子。

我先用错误的材料参数进行模拟，现在计算还是没能完全完成，但是会 ...

模型错误真是各种各样，我真回答不了这个问题。把你作过的，都检查检查吧。

adriano198707

Cob 发表于 2012-3-28 23:16
模型错误真是各种各样，我真回答不了这个问题。把你作过的，都检查检查吧。 ...

你好，因为我之前先做了一个力学分析，分析结果还是可以接受的，当我增加了耦合分析之后出现了错误。
我的耦合分析关键字参考了各种网上的资料的关键字，用不同的关键字试了很多遍。

不知能否帮我看下我K文件耦合部分的关键字？
谢谢了

zmh0611

mark，收藏了！！！