Stress hessian computation in FEM
介绍
在进行软体模拟时,如果使用牛顿法计算最优的下降方向,需要计算能量密度函数关于位置 的Hessian矩阵,即。其中是一个四面体的四个顶点的位置。
在进行软体模拟时,如果使用牛顿法计算最优的下降方向,需要计算能量密度函数关于位置 的Hessian矩阵,即。其中是一个四面体的四个顶点的位置。
首先,已知
这里把后三列写成一个3x3的矩阵。的分量表示为,的分量表示为,则能量密度函数关于位置 的梯度为:
Projective Dynamics是一种用于软体模拟的方法,算法分为local step和global step两个部分。其中local step可以对于每个四面体约束并行计算,global step只需要求解一个线性方程组,而他的矩阵非常特殊,是一个Gram矩阵,因此可以预先用Cholesky分解。本文主要介绍local step的实现的推导过程。用Corotated strain model的CUDA的实现作为例子。