ABSTRACT

The paper describes an application of parallel high performance computing for the finite element modelling of large scale 3D problems of geomechanics. This application concerns computations of stress fields induced by mining and requires solution of large linear systems with nearly four million of unknowns. Moreover, these systems may be singular and slightly inconsistent. The solution of these systems is performed by the conjugate gradient methods which are stabilized by projections when solving singular systems and preconditioned by the displacement decomposition and incomplete factorization. These methods, with either fixed or variable preconditioning, are implemented on a powerful workstation and a distributed memory parallel multicomputer. The performance of these solvers is shown.