Posted in Geotechnical Engineering

Rankine and Coulomb Earth Pressure Coefficients

For retaining walls, computing these is important; but the textbooks and reference books are frequently confusing and sometimes wrong. These formulae, derived using Maple, should clear up a few things, although they’re a) not always in the format you’re used to and b) subject to the Terms and Conditions of this site.

Let’s start with a diagram and the basic Coulomb formulae, from NAVFAC DM 7.02.

It’s important for to have a nomenclature chart for any lateral earth pressure coefficient formulae.

Coulomb Active Coefficients

All angles nonzero:

K_{{a}}=\left (\cos(\phi-\theta)\right )^{2}\left (\cos(\theta)\right)^{-2}\left (\cos(\theta+\delta)\right )^{-1}\left (1+\sqrt {{\frac {\sin(\phi+\delta)\sin(\phi-\beta)}{\cos(\theta+\delta)\cos(\theta-\beta)}}}\right )^{-2}

Vertical Wall:

K_{{a}}=\left(\cos(\phi)\right)^{2}\cos(\beta)\left(\cos(\delta)\cos(\beta)+2\,\cos(\delta)\sqrt{{\frac{\sin(\phi+\delta)\sin(\phi-\beta)}{\cos(\delta)\cos(\beta)}}}\cos(\beta)+\sin(\phi+\delta)\sin(\phi-\beta)\right)^{-1}

Level Backfill:

K_{{a}}=\left (\cos(\phi-\theta)\right )^{2}\left (\cos(\theta)\right)^{-1}\left (\cos(\theta+\delta)\cos(\theta)+2\,\cos(\theta)\cos(\theta+\delta)\sqrt {{\frac {\sin(\phi+\delta)\sin(\phi)}{\cos(\theta+\delta)\cos(\theta)}}}+\sin(\phi+\delta)\sin(\phi)\right )^{-1}

Vertical Wall and Level Backfill:

K_{{a}}=\left(\cos(\phi)\right)^{2}\left(\cos(\delta)+2\,\cos(\delta)\sqrt{{\frac{\sin(\phi+\delta)\sin(\phi)}{\cos(\delta)}}}+\sin(\phi+\delta)\sin(\phi)\right)^{-1}

Coulomb Passive Coefficients

All angles nonzero:

K_{{p}}=\left (\cos(\phi+\theta)\right )^{2}\left (\cos(\theta)\right)^{-2}\left (\cos(\theta-\delta)\right )^{-1}\left (1-\sqrt {{\frac {\sin(\phi+\delta)\sin(\phi+\beta)}{\cos(\theta-\delta)\cos(\theta-\beta)}}}\right )^{-2}

Vertical Wall:

K_{{p}}=\left(\cos(\phi)\right)^{2}\cos(\beta)\left(\cos(\delta)\cos(\beta)-2\,\cos(\delta)\sqrt{{\frac{\sin(\phi+\delta)\sin(\phi+\beta)}{\cos(\delta)\cos(\beta)}}}\cos(\beta)+\sin(\phi+\delta)\sin(\phi+\beta)\right)^{-1}

Level Backfill:

K_{{p}}=-\left (\cos(\phi+\theta)\right )^{2}\left (\cos(\theta)\right )^{-1}\left (-\cos(\theta-\delta)\cos(\theta)+2\,\cos(\theta)\cos(\theta-\delta)\sqrt {{\frac {\sin(\phi+\delta)\sin(\phi)}{\cos(\theta-\delta)\cos(\theta)}}}-\sin(\phi+\delta)\sin(\phi)\right )^{-1}

Vertical Wall and Level Backfill:

K_{{p}}=\left(\cos(\phi)\right)^{2}\left(\cos(\delta)-2\,\cos(\delta)\sqrt{{\frac{\sin(\phi+\delta)\sin(\phi)}{\cos(\delta)}}}+\sin(\phi+\delta)\sin(\phi)\right)^{-1}

Rankine Active Coefficients

All angles nonzero (except obviously \delta = 0 ):

K_{{a}}=\left (\cos(\phi-\theta)\right )^{2}\left (\cos(\theta)\right)^{-3}\left (1+\sqrt {{\frac {\sin(\phi)\sin(\phi-\beta)}{\cos(\theta)\cos(\theta-\beta)}}}\right )^{-2}

Vertical Wall:

K_{{a}}=\left (\cos(\phi)\right )^{2}\cos(\beta)\left (\cos(\beta)+2\,\sqrt {{\frac {\sin(\phi)\sin(\phi-\beta)}{\cos(\beta)}}}\cos(\beta)+\sin(\phi)\sin(\phi-\beta)\right )^{-1}

Level Backfill:

K_{{a}}=\left (\cos(\phi-\theta)\right )^{2}\left (\cos(\theta)\right)^{-1}\left (\left (\cos(\theta)\right )^{2}+2\,\left (\cos(\theta)\right )^{2}\sqrt {{\frac {\left (\sin(\phi)\right )^{2}}{\left (\cos(\theta)\right )^{2}}}}+1-\left (\cos(\phi)\right )^{2}\right )^{-1}

Vertical Wall and Level Backfill:

K_{{a}}={\frac {\left (\cos(\phi)\right )^{2}}{\left (1+\sqrt {\left (\sin(\phi)\right )^{2}}\right )^{2}}}

Rankine Passive Coefficients

All angles nonzero (except obviously \delta = 0 ):

K_{{p}}=\left (\cos(\phi+\theta)\right )^{2}\left (\cos(\theta)\right)^{-3}\left (1-\sqrt {{\frac {\sin(\phi)\sin(\phi+\beta)}{\cos(\theta)\cos(\theta-\beta)}}}\right )^{-2}

Vertical Wall:

K_{{p}}=-\left (\cos(\phi)\right )^{2}\cos(\beta)\left (-\cos(\beta)+2\,\sqrt {{\frac {\sin(\phi)\sin(\phi+\beta)}{\cos(\beta)}}}\cos(\beta)-\sin(\phi)\sin(\phi+\beta)\right )^{-1}

Level Backfill:

K_{{p}}=-\left (\cos(\phi+\theta)\right )^{2}\left (\cos(\theta)\right )^{-1}\left (-\left (\cos(\theta)\right )^{2}+2\,\left (\cos(\theta)\right )^{2}\sqrt {{\frac {\left (\sin(\phi)\right )^{2}}{\left (\cos(\theta)\right )^{2}}}}-1+\left (\cos(\phi)\right )^{2}\right )^{-1}

Vertical Wall and Level Backfill:

K_{{p}}={\frac {\left (\cos(\phi)\right )^{2}}{\left (1-\sqrt {\left (\sin(\phi)\right )^{2}}\right )^{2}}}

Posted in Geotechnical Engineering, Pile Driving Equipment

Vulcan and Sheet Piling — vulcanhammer.info

The development of sheet piling and the Warrington-Vulcan hammers was about the same time (along with concrete piles, other types of steel piles, and the Engineering News Formula.) The years leading up to World War I were ones of rapid development in the marine construction and deep foundation industries, and Vulcan was in the middle […]

Vulcan and Sheet Piling — vulcanhammer.info
Posted in Geotechnical Engineering

The Difference Between Flamant’s and Terzaghi’s Solution for Line Loads Behind Retaining Walls

Generally speaking, in Soil Mechanics courses elastic solutions on semi-infinite half spaces are presented to allow the geotechincal engineer to estimate the stresses induced in a soil by a load at the surface. Also presented in Soil Mechanics courses are charts and equations to estimate the lateral pressures on retaining walls induced by vertical loads on the surface behind the retaining walls. How either of these came into being is generally not explained; additionally, the fact that they are related is not explained either. The purpose of this article is to explain that relationship and, in the case of the retaining walls, how the original equations have been modified to suit experimental data. For simplicity’s sake, the discussed will be restricted to line loads.

The Original Equations: Flamant

As pointed out by Verruijt, the original equation for the stresses (vertical, horizontal and shear) were first set forth by Flamant in 1892. The equations are shown below.

Flamant Solution for a Line Load on a Semi-Infinite Elastic Medium, from Verruijt

The derivation of this equation can be found here or in Verruijt. Note carefully that the vertical load induces stresses in the horizontal direction (and in shear) as well as the vertical direction.

Verruijt then shows the following:

Mirror Flamant Loads on either side of an arbitrary dividing line and the computation of stress against a rigid wall.

Basically, by static equilibrium, if you replace the left half of the two mirrored loads (on the left) with a rigid wall, the horizontal stresses would be the same on the centre axis/wall as induced by two line loads. The resulting stresses and the resultant for the stress distribution are given below.

Flamant’s Equation for Horizontal Stresses Against a Rigid Wall, Including the resultants

The Modified Equations: Terzaghi

Generally, however, these equations are not presented in books such as DM 7, Sheet Pile Design by Pile Buck or others in this exact form. The following chart (frequently copied) comes from the Soils and Foundations Manual:

The line load is in the upper right hand corner. For values of m (the ratio of the distance from the line load to the wall over the height of the wall) greater than or equal to 0.4, the two results are the same. For those less than that, they are different. (In Verruijt’s notation, m = a/z.) The difference is because most books use the formulation of Terzaghi (1954). He explains the difference as follows:

However, the application of the line load tends to produce a lateral deflection of the vertical section, and the flexural rigidity of the bulkhead interferes with that deflection…However, for values smaller than (m=)0.4, the discrepancy between observed and computed values increases with decreasing values values of m…

From Terzaghi, K. (1954) “Anchored Bulkheads.” Transactions of the American Society of Civil Engineers, Vol. 119, Issue 1.

The whole issue of the flexibility of the retaining wall has been the chief complicating factor in this discussion, going back to Spangler’s tests in the 1930’s.

Comparing the Two Solutions

As an illustration, consider the pressure distribution situation when m=0.3:

Comparison of Flamant/Verruijt and Terzaghi Solutions for Line Load Pressures, m = 0.3

The pressures have been made dimensionless for generalisation. The Flamant solution comes to a higher peak nearer to the surface but falls off more rapidly down the wall. Terzaghi’s solution is more evenly distributed.

Now consider the situation at m=0.5:

Comparison of Flamant/Verruijt and Terzaghi Solutions for Line Load Pressures, m = 0.5

The two are identical in this range.

We can also consider the resultants as well:

Comparison of Flamant/Verruijt and Terzaghi Solutions for Line Load Resultants

The y-axis is made dimensionless by dividing the resultant by the vertical line load. For values of m less than or equal to 0.4, the results are different; for greater, they are the same.

In general, we can say that Flamant’s original formulation is more conservative. In the event that a deeper understanding of the interaction of surface loads with a retaining wall is desired, a finite element analysis needs to be done.

Posted in Geotechnical Engineering

G-I Geo Legends Series Interviews Harry Poulos

Most geotechnical engineers will recognize the name Harry Poulos. The Geo-Institute’s Geo-Legends series recently posted an interview with Professor Poulos of Coffey Engineering and the University of Sydney. He has worked on the foundations of some of the most well-known skyscrapers in the world in Dubai and elsewhere.

G-I Geo Legends Series Interviews Harry Poulos