FEM:刚度矩阵的对称性

计算科学 有限元
2021-12-04 21:34:26

在我学习有限元方法的第一门课程之后,我了解到元素和全局刚度矩阵是对称的。刚度矩阵对称的原因是什么?是因为弱形式的双线性项的对称属性吗?

2个回答

如果双线性形式是对称的(包括您选择施加边界条件的任何方式),并且您选择使用相同类型的试验和测试函数,则矩阵将是对称的。您可以通过写出问题的形式来看到这一点:

Kij=b(ϕi,ϕj)=b(ϕj,ϕi)=Kji

如果算子,刚度矩阵是对称的L的 PDE 是自伴随的,即如果你有Lf,g=f,Lg对于任何一对函数(f,g)在合适的函数空间中,其中u,v表示两个函数之间的内积u,v, 例如Ωuvdx(L2内积)。

刚度矩阵是通过投影 PDE 获得的Lf=g在有限元函数基础上(Φi), 形式如下:i,Lf,Φi=g,Φi(Galerkin 公式),然后在表格中寻找解决方案f=αjΦj,通过将第二个方程注入第一个方程,您可以得到离散方程,形式为Kα=Bβ, 在哪里α是未知数的向量(函数基础中解的坐标(Φi)),β是投影到函数基础上的 RHS (βi=g,Φi),B是质量矩阵 (bi,j=Φi,Φj) 和矩阵的系数ki,j=LΦi,Φj. 因此,如果运营商L是自伴随的,你有LΦi,Φj=LΦj,Φi和矩阵K是对称的。

这是最常见的偏微分方程的情况,但是当物理变得更复杂时也有例外......