本压缩文件中包含求解二维标量波动方程初边值问题的线性元三角形质量集中有限元方法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