A discrete element method for evaluating the seismic performance of concrete gravity dam-reservoir systems under main shock-aftershock events

Abstract

Dams are crucial for water supply, flood prevention, and hydroelectric power generation. Often located in seismically active regions, they are vulnerable to main shock-aftershock (MS-AS) sequences, which can compromise structural integrity and hydraulic safety. Critical aspects of dam response to MS–AS events remain unclear, particularly the required rest time between successive events and threshold AS-to-MS intensity measure ratios that could serve as predictors of additional damage. This study addresses these gaps by analyzing concrete gravity dam–reservoir systems of three heights (50 m, 100 m, and 150 m) using the developed discrete element–based approach coupled with displacement/pressure-based mixed finite elements for the reservoir. Empirical rest time equations were derived from 124 as-recorded ground motions, while seismic performance under varying intensity levels was evaluated using 14 as-recorded MS–AS sequences. Damage was quantified using discrete indices of base crack length, maximum base crack width, and maximum total upstream crack width. Results indicate that AS primarily propagate existing cracks at lower intensities, whereas higher intensities generate new cracks along the upstream face, increasing crack widths by 25–30% on average. The 50 m high dam remained within the mild damage category, while taller dams occasionally reached moderate levels, posing potential seepage risks. Threshold AS-to-MS ratios for four different intensity measures were identified. These findings provide mechanistic insight into crack propagation under MS-AS events, providing practical guidance for post-earthquake dam safety assessment, inspection prioritization, and incorporating sequential seismic effects into design and emergency planning.