This study examines the seismic stability of an underground powerhouse cavern located in the Lesser Himalayan region of Nepal. Both static and seismic loading conditions are analyzed using the finite element method (FEM) and the distinct element method (DEM). Rock mass properties are derived from field investigations and laboratory testing, while empirical correlations are applied to estimate rock mass strength and deformation modulus. Pseudo-static analyses are performed using the FEM-based software Rock and Soil-2-Dimensionsl (RS2) Version 11.027, and dynamic analyses are conducted using the DEM-based software Universal Distinct Element Code (UDEC) Version 5.0 to evaluate deformation and stress redistribution around the cavern. Seismic fragility curves are developed to quantify the probability of damage under varying seismic intensities. Results indicate that a peak ground acceleration (PGA) of 0.25 g increases cavern wall deformation by approximately 15–20 mm compared to static conditions. Fragility analysis shows a probability exceeding 68% for slight damage, while the probability of collapse remains low at approximately 1.7%. Seismic loading also significantly alters stress redistribution along the cavern boundary. Overall, the combined use of numerical modeling and fragility analysis provides a probabilistic framework for assessing seismic risk in underground caverns, offering valuable insights for the design and safety evaluation of hydropower projects in seismically active Himalayan regions.
Adhikari et al. (Thu,) studied this question.