构造各个矩阵
内力f1 = f0 + K * dx
K为tangent stiffnes matrix
dx为当前变形下的一个微小位移
对coratational模型调用testTangentStiffnessMatrix函数，发现原来的wrapped stiffness matrix并不准确。
只有使用准确的tangent stiffness matrix ，误差才为位移的高阶无穷小。（与步长长度平方比例近似恒定）
当位移分量减小到1e-6左右时，由于数值误差，误差不再减小，反而增大。
构造各个矩阵
内力f1 = f0 + K * dx
K为tangent stiffnes matrix
dx为当前变形下的一个微小位移
对coratational模型调用testTangentStiffnessMatrix函数，发现原来的wrapped stiffness matrix并不准确。
只有使用准确的tangent stiffness matrix ，误差才为位移的高阶无穷小。（与步长长度平方比例近似恒定）
当位移分量减小到1e-6左右时，由于数值误差，误差不再减小，反而增大。