Three Dimensional Modeling of the Wedge Pillar Portion of the WIPP Geomechanical Evaluation (Room G) In Situ Experiment

INTRODUCTION The Waste Isolation Pilot Plant (WIPP) is a research facility located in Southeastern New Mexico. The WIPP is being developed by the U. S. Department of Energy (DOE) to demonstrate the safe storage of defense-related radioactive wastes in bedded salt. Since the WIPP is a research facility. a number of large-scale in situ experiments have been planned and are currently under construction (Munson, 1983). These experiments were designed to study various aspects of nuclear waste storage in bedded salt such as mechanical responses and creep closures of drifts, thermal response due to heated canisters, and thermally-induced fluid migration. One purpose of the WIPP is to develop and demonstrate a general predictive structural analysis capability for a bedded salt repository. One set of experiments, called the Geomechanical Evaluation (Room G), is heavily instrumented to study the creep around several rooms with different configuration. A Layout of the Geomechanical Evaluation room is shown in Figure 1. One portion of it consists of wedge shaped pillar where two drifts intersect at an angle of 7.5 degrees. The pillar can be seen at the bottom of Figure 1. One purpose of the wedge pillar experiment is to study progressive pillar failure in the tip region. Another is to determine how nonuniform pillar thickness affects creep closure of the drifts. An important aspect of the experiments is the correlation between experimental data and corresponding pretest finite element analyses of the site. The finite element simulations serve two purposes. First, the computer simulations aid in understanding the experimental results by providing calculated stress. strain and displacement fields that cannot be measured directly. Second, comparison of calculated and measured drift closures serves as a method for validating or improving the finite element models and the constitutive models employed. The wedge pillar geometry required a 3-D finite element creep calculation. Results from these 3-D calculations will be presented in this paper. A 2-D plane strain double drift model which approximates the wedge pillar geometry at a slice perpendicular to the drift has also been performed and will be compared to the 3-D results in this paper. FINITE ELEMENT COMPUTER PROGRAMS Two finite element computer programs. JAC (Biffle, 1984) and JAC3D (Biffle. 1986), were used. JAC is a 2-D finite element program developed for quasi-static analysis of non-linear solids. It employs the conjugate gradient iterative technique to obtain a solution Spatial integration is performed using a single gauss point in each four node quadrilateral element. An hourglass viscosity technique is used to control the zero energy modes that occur with single point integration. The single point integration combined with the explicit nature of the program and exploitation of CRAY-1 computer architecture results in very efficient execution J4: results were compared to results obtained with eight other 2-D structural creep computer codes in the second WIPP benchmark exercise (Morgan, 1981) JAC3D was derived from JAC to treat 3-D finite element models and has many of the same characteristics including single point integration and hourglass viscosity. These characteristics have made 3-D creep analyses more reasonable by significantly reducing computation time The first exercise of the creep capability