Fault rupture propagation through stratified sand–clay deposits and engineered earth structures: a meshfree and critical-state modeling approach

Permanent deformation and uplift caused by fault rupture is one of the most significant hazards posed by earthquakes on the built environment. In this paper, we use smoothed particle hydrodynamics (SPH) to explore the effects of soil layering or stratification on the trajectories and deformation patterns caused by rupturing reverse faults in bedrock, as well as in the foundations of engineered earth structures. SPH is a continuum meshfree numerical method highly adept at modeling large deformation problems in geotechnics. Through the use of constitutive models involving softening behavior as well as critical state type models, we isolate the effects of rigid body rotation from critical state behavior of soil in helping explain the frequently observed rotation of shear bands emanating from the bedrock fault. This analysis is facilitated by the fact that the SPH method allows us to track the propagation of shear bands over substantial amounts of vertical uplift (more than 50% of the total height of the soil deposit), far beyond many previous computational studies employing the finite element method (FEM). We observe and characterize various emergent features including fault bifurcations, stunted faults, and tension cracking, while providing insights into practical guidelines regarding the potential surface distortion width, and the critical amount of fault displacement required for surface rupture depending on the multilayered constitution of the soil deposit. Finally, we predict the expected amount of surface distortion and internal damage to earthen embankments depending on varying fault location and soil makeup.