请教:关于consistent tangent matrix

willnessless

通常,consistent tangent matrix可以这样计算

Cct=R-RaaR/aRa                                  (1)

R={I+lamda*C*[Par(a)/Par(xigma)]} 的逆矩阵*C.  (2)

而传统的弹塑性刚度矩阵Cep就是把 (1)中的R 换成弹性矩阵 C.
------------------------------------------------------------------------------------------

问题是:

对于一些特殊的屈服面, 例如抗拉强度面, F=P+常数=0, 它的法向a就是常数矢量, Par(a)/Par(xigma) 就是0. 代入(2)---->R=C

也就是说,consistent tangent matrix就是传统的弹塑性刚度矩阵Cep. 那如果对于D-P模型,设置了抗拉强度, 其实用不用consistent tangent matrix意义并不大啊. 因为系统的收敛速度收到那个抗拉强度面的限制, 收敛速度相当于采用传统的弹塑性刚度矩阵Cep的速度.(当然是指有积分点的应力进入拉屈服状态时).

不知道我的理解对不对? 因为对于D-P, 如果收敛速度拉屈服面控制,那就不用去求那个R矩阵,这样更省事啊.

[ 本帖最后由 willnessless 于 2008-11-6 10:30 编辑 ]

willnessless

Is there anybody interested in this?

willnessless

点名啦, aogudusi版主请指点, 发表高见....

版主是做过D-P的umat的,而且也采用了设置抗拉强度的办法来处理顶点,如果我没记错的话.

willnessless

顶下.

顺便请版主帮个忙, 把这个贴划归到"用户子程序与二次开发"那个分栏. 要不老是不容易找. 多谢啦...

willnessless

看来理解错了,对于特殊的屈服面,consistent tangent matrix确实就是弹塑性刚度矩阵Cep, 但是收敛速度应该还是跟用consistent tangent matrix一样快. 所以在抗拉面, 并不会影响收敛的速度.