请教vumat中如何求解J3不变量

si13

初次接触vumat，程序也编的不多，在此请教几个问题
1、子程序中不是有提供 SPRINC 和 VSPRINC 函数吗？这两个函数看帮助描述的应该是
求主应力的函数，为什么在论坛里搜的结果里面似乎都表示子程序里面不能直接求主应力，
需要自己使用“盛金法”、“jacobi法”等编程序求出，难到是我对这两个函数的作用理解有误？
2、在程序中如何求解J3不变量？
这个我原本打算先求偏应力张量的三个主应力，相乘得到
但是在探索求解主应力的方法上发现有大牛指明可使用J2、J3、平均应力等编程求出主应力
这就陷入了循环的怪圈

这让我困惑不已，我理论基础不深，上述可能在理论上理解有误，
请各位路过高手指点迷津！

shawn2008

我常用下面这篇论文的方法求解主应力

Linear transfomation-based anisotropic yield functions  Original Research Article
International Journal of Plasticity, Volume 21, Issue 5, May 2005, Pages 1009-1039
F. Barlat, H. Aretz, J.W. Yoon, M.E. Karabin, J.C. Brem, R.E. Dick
Open Show preview  |   PDF (718 K)   |   Related articles  |  Related reference work articles

http://www.sciencedirect.com/sci ... i/S0749641904001160

si13

非常感谢！
J3不变量我也知道如何计算了，直接根据定义来就成
编程太少  经验不足~~