Dynamic Slope Stability Back
Geotech

Dynamic Slope Stability — Newmark Method

Enter slope geometry, soil strength parameters, and earthquake input to compute static FS, yield acceleration, and seismic permanent displacement in real time.

Slope Parameters
Slope angle β35°
Cohesion c10 kPa
Friction angle φ30°
Unit weight γ18 kN/m³
Slope height H15 m
Earthquake Parameters
Peak acceleration amax/g0.25
Earthquake duration20 s
Static FS
Yield accel. ky (g)
Permanent disp. D (cm)
Critical radius R (m)
Risk: Low

Key Equations

Yield accel.: $k_y = (FS-1)\sin\beta$

Permanent displacement (Ambraseys-Menu):

$$D = 0.087\frac{v_{max}^2}{a_{max}}\left(\frac{k_y}{a_{max}/g}\right)^{-2.53}$$

FS via infinite slope limit equilibrium

Permanent displacement D vs peak acceleration amax/g (current ●)
Factor of safety FS vs slope angle β (current ●)

What is the Newmark Method for Slope Stability?

🧑‍🎓
What exactly is the Newmark method? I've heard of static slope stability, but what makes this "dynamic"?
🎓
Basically, it's a way to estimate how much an earthquake can make a slope slide. In static analysis, you just check if the slope is stable under gravity. The Newmark method adds the shaking from an earthquake, treating the slope like a rigid block that only slides when the shaking force exceeds a critical threshold. Try moving the "Peak Acceleration" slider in the simulator to see how stronger shaking affects the results.
🧑‍🎓
Wait, really? So the slope isn't failing the whole time during the quake? What's this "yield acceleration" I see in the output?
🎓
Exactly! That's the key insight. The yield acceleration ($k_y$) is the seismic acceleration level needed to just overcome the slope's strength and start it sliding. Below that, the slope holds. Above it, it accumulates permanent displacement. In practice, for a dam or embankment, you'd calculate this to see if shaking from a design earthquake would cause dangerous movement. Change the soil's friction angle ($\phi$) and cohesion ($c$) above to see how they directly influence $k_y$.
🧑‍🎓
So the final "Permanent Displacement" number is the total slide distance? How do you get that from just the acceleration and duration?
🎓
Great question. The simulator uses an empirical formula (like Ambraseys-Menu) that correlates displacement with the *excess* shaking. It integrates the effect of every moment the earthquake acceleration exceeds the yield acceleration. A common case is a highway cut in a seismic zone; engineers use this displacement value to judge risk. Try reducing the slope angle ($\beta$) and notice how both the yield acceleration and the predicted displacement change.

Physical Model & Key Equations

The analysis starts with a static Factor of Safety (FS) against sliding. The dynamic tipping point is the yield acceleration, which scales with how much stronger the slope is than the minimally stable condition.

$$k_y = (FS - 1) \sin \beta$$

Here, $k_y$ is the yield acceleration (in $g$), $FS$ is the static factor of safety, and $\beta$ is the slope angle. A lower $k_y$ means the slope is closer to failure and will start sliding under weaker shaking.

The permanent displacement is estimated using a regression model derived from analyzing records of real earthquake-induced slides. The Ambraseys-Menu equation is one common model.

$$D = 0.087 \frac{v_{max}^2}{a_{max}}\left( \frac{k_y}{a_{max}/g}\right)^{-2.53}$$

Where $D$ is the permanent displacement (cm), $v_{max}$ is the peak ground velocity, $a_{max}$ is the peak ground acceleration, $k_y$ is the yield acceleration, and $g$ is gravity. This shows displacement is highly sensitive to the ratio of yield to peak acceleration.

Real-World Applications

Dam Safety Assessment: Engineers use this method to evaluate the seismic stability of earthfill and tailings dams. For instance, if a predicted displacement exceeds 15 cm for a critical dam, it might trigger a redesign with stronger materials or flatter slopes to reduce risk.

Highway & Railway Embankments: When building transportation corridors in earthquake-prone areas like California or Japan, this analysis ensures that embankments won't suffer excessive settlement or collapse during a quake, blocking vital evacuation and supply routes.

Landslide Hazard Zoning: Geologists and planners apply the Newmark method regionally using GIS to map areas with high displacement potential. A common case is assessing landslide risk for neighborhoods built on steep hillsides after a major fault rupture.

Mine Slope Design: In open-pit mining, the stability of high, steep slopes is critical for worker safety and economic operation. Dynamic analysis helps plan mining sequences and evaluate the need for reinforcement if blasting or regional seismicity could trigger a slide.

Common Misunderstandings and Points to Note

When you start using this tool, there are several pitfalls that engineers, especially those with less field experience, often fall into. A major misunderstanding is thinking that the calculated permanent displacement D directly equals the collapse distance. For example, even if D=0.5m, it does not mean the entire slope will slide 0.5m all at once. The Newmark method provides an estimate of the "average" displacement due to the accumulation of shear strain. Actual failure involves this displacement concentrating locally or developing into a surface slide, so you should treat the D value as a relative indicator for risk comparison.

Next is the setting of input parameters. You must not use the "cohesion c" and "internal friction angle φ" directly from geotechnical investigation reports. During an earthquake, strength degrades (dynamic strength reduction) due to cyclic loading, so it's common to set them at about 70-80% of the static strength. For instance, if static tests yield c=30 kN/m² and φ=30°, for dynamic analysis you would typically use c=24 kN/m² and φ=24°. Be careful not to forget this correction in the tool, as it will lead to an overestimated factor of safety and an underestimated permanent displacement.

Finally, regarding seismic motion input. The tool requires you to input $v_{max}$ and $a_{max}$ as single values, but the period characteristics of the actual seismic wave are crucial. For example, long-period seismic motions affect deeper parts of the slope, increasing the risk of deep-seated slides, not just surface slides. In practice, the basic approach is a "multiple-case analysis," where you input several seismic waves considering the expected earthquake's source characteristics and ground amplification factors, and adopt the most critical result. Remember, this tool is intended for initial screening.

Related Engineering Fields

The concepts behind this dynamic slope stability analysis are actually closely linked to various engineering fields. The first that comes to mind is liquefaction assessment of soils. The "cyclic stress ratio" used in liquefaction evaluation and the "yield acceleration" for slopes share a common principle: comparing the increase in shear stress caused by seismic motion to the material's resistance. The understanding of dynamic strength you gain from slope analysis forms the foundation for comprehending the FL (liquefaction strength) in liquefaction.

Another field is seismic design of structures, particularly foundation and retaining wall design. The permanent displacement generated by a slope during an earthquake must be considered as a "forced displacement" acting on the foundations of retaining walls or bridge piers built on that slope. For instance, if the backfill soil behind a retaining wall is expected to move 10 cm, even if the wall itself doesn't collapse, you need to design joints or bearings that can accommodate that displacement. This is an important part of "soil-structure interaction."

Taking it a step further, there is a connection to block toppling analysis for rock slopes in rock mechanics. While the failure mechanisms differ between soil slope sliding and rock toppling/sliding, the core analytical approach—applying seismic motion as an external force, comparing it to resistance, and evaluating displacement—is the same. Furthermore, advanced numerical simulations like the Discrete Element Method (DEM) used in this field serve as the next step to verify results obtained from simplified methods like the Newmark method with more detailed failure mechanisms.

For Further Learning

If you become interested in the theory behind this tool, we recommend deepening your knowledge in three steps. The first step is to understand "Newmark's original paper and its physical meaning." Newmark initially proposed a sliding block model for dams. The origin is an extremely simple and elegant concept: "the block continues to slide with velocity only for the duration that acceleration exceeds the yield value." Firmly grasping this basic image is the shortcut to understanding more complex equations later.

Next, learn the "assumptions" behind the Ambraseys and Menu empirical equation used in the tool. This equation is a regression formula derived from statistical processing of numerous seismic waves. This means it is not universally applicable to all ground conditions or seismic motions; its accuracy is highest under certain specific conditions (e.g., specific seismic motion characteristics or soil types). Investigating the derivation process and limitations of this equation helps you develop an "engineer's eye" to critically evaluate calculation results rather than accepting them blindly.

As a final step, consider expanding your knowledge to more general methods like the "response displacement method" or "dynamic effective stress analysis using the Finite Element Method (FEM)." The Newmark method treats the slope as a single rigid block, but real ground deforms. Using FEM, you can perform coupled analyses that include the distribution of acceleration and shear stress within the slope, and even the generation and dissipation of excess pore water pressure due to permeability effects (changes in effective stress). Learning FEM after gaining sensitivity to parameters with this tool will make it easier to understand "why those elements and condition settings are necessary."