This paper describes the results of a 2D and 3D numerical modeling carried out on an unstable natural rock slope (volume = 110,000–335,000 m3) located in north-eastern Italy (Friuli Venezia Giulia Region). It was aimed at investigating the mechanical behavior and internal rock mass damage of the unstable slope before the collapse. The numerical simulations testify that the stability condition of the slope is close to the limit equilibrium (strength reduction factor, SRF = 1.03–1.13), as demonstrated by the considerable rock mass damage observed on the field. The overall mechanical behavior of the slope is mainly governed by the kinematic conditions of the secondary internal blocks, which, in turn, depend on the geometry and mechanical properties of the major discontinuities that delimit the adjacent blocks. Slope failure is achieved through internal rock mass damage represented by internal shear and tensile ruptures localized in correspondence with over-stressed zones. The failure mechanism is characterized by sliding along pre-existing discontinuities and inner damage in the form of the enucleation of shear bands that originate internal secondary failure surfaces and/or damaged rock mass zones. The stress–strain modeling predicts intense slope deformations in zones where rock mass damage actually occurred. This paper emphasizes the decisive connection between the geomechanical field survey and numerical modeling. The comparison of surface geological data acquired on the field with the mechanical indicators obtained from the numerical analyses can significantly improve knowledge of the rock mass damage process that involves unstable slopes approaching failure condition.
