本压缩文件中包含求解二维标量波动方程初边值问题的线性元三角形质量集中有限元方法MATLAB源程序。这里只解释下质量集中:一般二维波方程的常规有限元离散格式为MU_tt+KU=F,其中M称为质量矩阵,K为刚度矩阵,F为荷载向量;M中各元素都是通过一个内积积分的表达式计算而来,质量集中的意思就是通过特定的数值积分公式替代内积分使得M为对角矩阵,即除主对角线外其余元素全为0。关于质量集中的论文可参考“精读--Cohen_2001_Higher order triangular FE with mass lumping for the wave equation”。需要说明的是,本压缩包中还包含常规的线性元三角形有限元法,可在程序内自己手动选择使用质量集中或常规有限元。
1
正则性假设或拟一致条件是进行有限元理论分析的基本假设,本文在不满足该假设的条件下,采用一种更简便的方法给出了双线性元在求解二阶问题时的超逼近结果;在此基础上,又得到了其在中心点的点态超收敛结果;最后的数值试验验证了理论分析的正确性。
2021-05-12 10:02:54 948KB 自然科学 论文
1