在我学习有限元方法的第一门课程之后,我了解到元素和全局刚度矩阵是对称的。刚度矩阵对称的原因是什么?是因为弱形式的双线性项的对称属性吗?
FEM:刚度矩阵的对称性
计算科学
有限元
2021-12-04 21:34:26
2个回答
如果双线性形式是对称的(包括您选择施加边界条件的任何方式),并且您选择使用相同类型的试验和测试函数,则矩阵将是对称的。您可以通过写出问题的形式来看到这一点:
如果算子,刚度矩阵是对称的的 PDE 是自伴随的,即如果你有对于任何一对函数在合适的函数空间中,其中表示两个函数之间的内积, 例如(内积)。
刚度矩阵是通过投影 PDE 获得的在有限元函数基础上, 形式如下:(Galerkin 公式),然后在表格中寻找解决方案,通过将第二个方程注入第一个方程,您可以得到离散方程,形式为, 在哪里是未知数的向量(函数基础中解的坐标),是投影到函数基础上的 RHS (),是质量矩阵 () 和矩阵的系数. 因此,如果运营商是自伴随的,你有和矩阵是对称的。
这是最常见的偏微分方程的情况,但是当物理变得更复杂时也有例外......
其它你可能感兴趣的问题