弹性波正演模型的快速有限元算法

用有限元法做弹性波正演模型计算已进入实用阶段,但有限元法所要求的大存储量和计算速度慢的缺点依然存在。本文在深入分析有限元算法的基础上,注意到刚度矩阵,耦合矩阵在特殊单元中的对称性、耦合矩阵具有的转置性及某些特殊特征,采用正方形单元并将与节点相乘的元素归并到时间离散中运算,使计算速度提高一倍,使整体矩阵的存放量减少一半。文中最后给出了一个计算实例,证实该快速算法的正确性。