Underground excavations in hard rock subject to dynamic loading caused by nuclear detonation may experience violent collapse (rockbursting) similar to what is observed in deep mining operations. This proposal aims at developing a methodology to determine the circumstances for which rockburst will occur, based on the assumption that collapse of the cavity is triggered by a surface buckling of the walls of the excavations. The methodology will rely on a finite element discretization to determine the conditions of instability, by searching the conditions for which one of the values of the incremental stiffness matrix vanishes. The objectives of this phase I proposal are to assemble l: finite element code, implement constitutive equations to model pressure-sensitive dilatant rock, and validate the numerical model against theoretical solutions on the surface instability of a half-plane subject to lateral loads.