[求助]关于屈服准则编程程序 [复制链接]

小弟现在在编写一个Mohr-Coulomb的屈服准备，但是小弟不是很清楚应该按怎么一个程序去编写。 oNiS"\t  
目前只清楚，先提取出应力分量，带入屈服函数，计算是否大于0，若大于则代表屈服，然后按塑性流动法则计算出体应变增量和剪应变增量，但是然后呢？ XY%8yII6  
希望有懂的大侠们能指导下，小弟跪谢啊。小弟的QQ：148305732 "~FXmKcX  
加QQ跟在线回答都可以

不知你是用隐式求解还是显式求解，如用隐式求解，编写mohr-coulomb比较困难，因为不规则六面体在棱上容易发生歧义，需要定义许多判断函数，否则很容易造成结果发散。 VUx~Y'b  
非线性问题建议用显式求解方式。 j*6!7u.,K  
Y ,}p  
无论你用有限体积法或者fem方法，均可以获得单元或高斯点上的应力，然后计算主应力及主方向，将主应力代入屈服方程（包括拉伸屈服及剪切屈服），选择其中一个屈服方程，修正主应力，修正主应力的方式就是将屈服面外的应力拉回到屈服面上，具体可参考大部分有限元书籍或者flac3d的手册。接着是将修正后的主应力计算为全局坐标系下的应力分量，并返回到节点应力上。

你说的这个大致的流程我知道，我主要就是想知道下细致的步骤，修正角点我用Koiter方法，我现在在看一些塑性力学的书籍！