之前我们在做图形学有限元中,使用的是stvk的方法,通过应变计算能量密度
而更加在工程学中通用的方法需要建立起应变strain到应力stress的映射
1.各向同性材料
所谓各向同性材料,就是对一种材料来说,从各个方向对他所施加的力,所造成的形变效果是一样的
我们可以认为P是形变矩阵的一个函数,(P是第一皮耶尔应力,而在之前的模拟中P是通过Green strain计算出来的,而Green strain是一个关于F形变矩阵的一个函数)将F做一个SVD可以得到两个旋转和一个缩放的形变,旋转—拉伸—旋转,将旋转拉出去,应力主要基于拉伸这部分,也就是拉伸矩阵,基于此根据一些数学力学公式,我们将P用三个量表示
那么使用的对于能量关于应变的公式就有以下这些,其中在工程学中应用较多的是新胡克定律
其中绿色部分是用来抵抗拉伸的,蓝色部分是为了阻止体积的改变
其中fung model一般用于模拟人体组织
根据这些定义可以推导出应力矩阵
这个时候的伪代码
stvk和新胡可的比较