Through field measurements, it was shown that the shear stresses along a fully grouted rock blot were proportional to the differences between the displacement of a bolt rod and that of surrounding rock. Accordingly the theoretical equation for the displacement along a rock bolt, R(x), can be proposed as differential form; $ \frac{d^2R(x)}{dx^2}-aR(x)=-a^2U(x)$ where U(x) is the displacement of surrounding rock and α is a constant. In this study above equation was applied to the boundary element analysis which can simulate the supporting behavior of each rock bolt around a tunnel and the results of the implementation of boundary element technique were shown.

The effects of geologic structures such as rock joins and bedding planes on the thermal conductivity of a discontinuous rock mass are studied. The expressions for the equivalent thermal conductivities of jointed rock masses are derived and found to be anisotropic. The degree of anisotropy depends primarily on the thermal properties contrast between the joint phase and surrounding intact rock, the joint density expressed as volume fraction and the inclination angle of the joint. Within the context of 2-dimensional finite element heat transfer scheme, the isotherms around a circular hole are analyzed for both the isotropic and anisotropic rock masses in 3 different thermal boundary conditions. i.e. temperature, heat flux and convection boundary conditions. The temperature in the stratified anisotripic rock mass is greatly influenced by the thermal properties of the rock formation in contact with the heat source. Using the excavation-temperature coupled elastic plastic finite element method, analyzed is the thermo-mechanical stability of a circular opening subjected to 10$0^{\circ}C$ at a depth of 527m. It is found that the thermal stress concentration was enough to deteriorate the stability and form a plastic yield zone around the opening, in contrast to the safety factor greater than 2 resulted form the excavation-only analysis.

In this study, the elastic modulus and the initial stresses of the rock were calculated through back analysis of in-situ displacements measured during excavation of the underground caverns. Results from back analysis were employed to determine the redistributed stresses the displacements and relaxed zone in the rock around the caverns, which supplement the geological characterization results. To verify the reliability of the back analysis program the elastic modulus and the initial stresses were obtained from inputting the displacements calculated by FEM. These were compared with the assumed normalized stresses in FEM and were in a reasonable agreement with an error of more or less than 3%.

Natural stress measurements have been made at two sites at the depth of 280m from surface by means of stress relief overcoring methods using three directonal deformation gage. Attempts have been made to determine the state of natural stress in the rock and provide useful basic data to investigate the stress distribution and the determination of yield zone around powerhouse cavern. The magnitude and the direction of the miximum principal stress obstained from in-situ stress measurements is -96.1kgf/$\textrm{cm}^2$ and N38$^{\circ}$W, N35$^{\circ}$W respectively. Vertical stresses are in approximately agreement with the theoretical value. The ratio of measured to theoretical stresses are 85% at two sites. The ratio of average horizontal to vertical stresses(k=($\sigma$h)ave/$\sigma$v is 1.07.