Para explorar la distribución y la evolución de la roca circundante de las cavernas subterráneas de las centrales hidroeléctricas, se llevaron a cabo las pruebas de compresión dinámica para la muestra de granito y se estableció la relación constitutiva dinámica adecuada para la programación informática en términos de relación constitutiva de daño elástico-plástico estática basada en la deformación. Introduciendo cinco terremotos de diversos espectros de frecuencia, se realizó el análisis dinámico totalmente no lineal para algunas grandes cavernas subterráneas de la región de desfiladeros y montañas de la provincia china de Sichuan con el método de la diferencia explícita. Los resultados de la simulación muestran que, a medida que las ondas sísmicas de entrada se hacen fuertes gradualmente, la zona de daños se extiende generalmente en profundidad desde el centro de las paredes laterales. La zona de daños alcanza el valor máximo en las zonas inmediatamente posteriores a la llegada del pico de aceleración sísmica y se mantiene constante a partir de entonces. Según el análisis inicial, si la frecuencia dominante de la onda de entrada se acerca a la frecuencia natural de las cavernas subterráneas, el área máxima de la zona de daños aumentará. Estos resultados pueden proporcionar datos prácticos para el diseño y la construcción antisísmica.