Implementation Of Finite Element-- Boundary Integral Linkage Algorithms For Rock Mechanics Applications

INTRODUCTION Methods of linking boundary integral (BI) solutions with finite element (FE) solutions have been well described theoretically in other publications (e.g., Zienkiewicz, et al., 1977). The purpose of this paper is to point out some practical aspects of implementing these linkage algorithms for rock mechanics applications and to describe an alternative linkage algorithm. A not so obvious motivation for using FE-BI linkages is illustrated in Figure 1. The cost of modeling the excavation of the crusher station by the BI method can equal or exceed the cost of an equivalent FE model. This is due to the large surface area which must be modeled by the BI method. In such cases, a more economical approach might be to model the region immediately adjacent to the crusher station by finite elements and link the solution at the periphery of the FE model with a solution for the exterior region. Further economies can be realized if some flexibility in the type of exterior solution is possible. For example, it would be desirable to use a coarse discretization of the FE-BI interface for a BI solution, which would be less demanding of computer storage and time, or to use an analytical solution for the exterior region. The alternative linkage algorithm to be discussed allows such flexibility. BOUNDARY INTEGRAL METHODS There are basically two formulations of BI methods: the direct and the indirect. FE-BI linkage algorithms depend on which formulation is used. The purpose of this section is to briefly describe each method.