当选择static solve的时候,求解的流程如下:1.获得内力
2.qresidual = 外力-内力,qdelta = qresidual, qdelta的非约束元素赋给bufferConstrained
3.用tangentStiffnessMatrix给systemMatrix A赋值
4.求解A * buffer = bufferConstrained
5.把buffer指定为qdelta不受约束的部分
由此可见,实质上是求解 tangentStiffnessMat * dx = 外力 - 内力
而tangentStiffnessMatrix*dx等于内力增量,于是式子就是 内力 + 内力增量 = 外力
实际上就是牛顿法求平衡位置
看ImplicitBackwardEulerSparse关于static solve的代码
原文:http://www.cnblogs.com/dydx/p/4243844.html