三维弹性波有限元模拟在微机上的实现

The finite element modelling of 3-D elastic wave implemented in microcomputer

1. 中国地质科学院岩石圈研究中心, 北京100037;2. 长春地质学院地球物理系, 长春130062

Researeh Center of Lithosphere, China Academy of Geological Sciences, Beijin8 100037

针对三维弹性波有限元模拟要求计算机内存大和计算时间长的困难,本文采用一种特殊剖分法和刚度矩阵的上半带等带宽压缩存贮法,并采用集中质量矩阵和中心差分时间积分显格式的结合,极大地节省了内存,提高了运算速度,在微机上实现了三维弹性波有限元模拟的计算。理论模型的试算验证了算法的正确性和可行性。

In order to deal with large computer storage and long time of calculation in the finite element modelling of 3-D elastic wave, we used a particular segment method and the bandwidth compact storage method of the upper semi-band of stiffness matrix. Simultaneously, we make use of the central mass matrix and the explicit formula of center time difference. Using these methods to save computer storage and increase calculation speed, the finite clement modelling of 3-D elastic wave is implemented in micro-computer. The correctness and efficiency of this algorithm is proved by the results of theoretical modelling.