Posted in Geotechnical Engineering, Soil Mechanics

Getting to the Bottom of Terzaghi and Peck’s Lateral Earth Pressures for Braced Cuts

One of those “things” in geotechnical engineering that looks like “settled science” but may not be is the whole business of lateral earth pressures for braced cuts. (An example of one is shown at the right.) Textbooks of all kinds (including Soils in Construction) show pressure profiles and solution techniques that are “definitive.” Or are they? This article is more about asking questions than delivering another round of “definitive” answers, but hopefully it will at least spark some thought and perhaps make practitioners more careful in their application of these methods.

The Basic Problem

Based on experience, first in Berlin and later in Chicago, Terzaghi (and later also Ralph Peck) developed a set of pressure distributions as shown at the left. These are at variance with those usually associated with retaining walls in general and sheet pile walls in particular. The theory behind these (certainly for the clays) had its genesis in Rankine theory adapted for cohesive soils, but the distribution is rather different. These distributions have been reproduced in many textbooks and reference books, including Soils in Construction and Sheet Pile Design by Pile Buck.

It’s worth noting that there are other pressure distributions that have been formulated other than the ones shown above, as outlined by Boone and Westland (2005).

Along with the distributions came the method of using them: the “hinged method” of analysing the sheet pile wall. Strictly speaking the sheet pile wall is a continuous beam with multiple supports. Since there are usually more than two struts and supports, to use a continuous beam requires a statically indeterminate beam. Applying hinges (as shown at the right) can make the beam statically determinate. Although methods for solving statically indeterminate beams existed in Terzaghi’s day (remember it was geotechnical methods which then and now lag the rest of civil engineering in advancement,) converting the problem to a statically determinate one was convenient for computational purposes.

But then comes the kicker: as generally presented, if the distributions above are used, they must be used with the hinged method, even when analysing a braced cut wall using a braced cut method with a continuous beam is nearly trivial now. Why is this?

How It Came About

Like so many things in civil engineering, the investigation of pressures on braced cuts came about as a result of tragedy. As noted by Rogers (2013):

On the evening of December 1, 1938 Terzaghi delivered a terse lecture titled “The danger of excavating subways in soft clays beneath large cities.” The lecture focused on his recent experiences with construction of the Berlin Subway, which was hampered by a high water table in running sands. These conditions had contributed to the sudden failure of a shored excavation which killed 20 workers in August 1935. He made a convincing case for proper geotechnical oversight during construction if similar tragedies were to be avoided in Chicago.

The lecture with its graphic images of the dead bodies beneath the collapsed bulkhead along the Hermann Goring Strasse succeeded in scaring his audience to death, and promptly found the State Street Property Owners’ Association and City of Chicago bidding for Terzaghi’s services. The City wanted him to advise them on how best to monitor progress of excavations and ground settlement, differentiating what structural or architectural damage was caused by subway construction.

Both Terzaghi and Ralph Peck ended up doing the monitoring. The soils in Chicago were predominantly soft clays, so the earth pressures were different. Much of the theory and application behind this is documented in Peck (1943). (Interesting side note: one of the lines of the subway ran past the site of Vulcan’s old Milwaukee Avenue plant where the first Warrington-Vulcan steam hammers were designed and built.)

One might ask, “How did they come up with the earth pressure distributions?” They did so–and this is the key to the problem–by measuring the reactions on the braces. They did this in the face of the fact that, as is usual with braced cuts, the braces were put in successively with excavations, and that much of the movement of the wall–and thus the mobilisation of the earth pressures–was in place before the braces were installed. (It is easy to forget the importance of that mobilisation, but both Terzaghi and Peck were well aware of it and its effects.) To translate the loads on the braces into a pressure distribution, they adopted Terzaghi’s procedure from the Berlin subway as follows (Peck (1943)):

The vertical members of the sheeting are assumed to be hinged at each strut except the uppermost one, and a hinge is assumed to exist at the bottom of the cut. The abscissas of the pressure diagram “A” represent the intensity of horizontal pressure required to produce the measured strut loads. A study of such diagrams for all of the measured profiles disclosed that the maximum abscissa never exceeded the value KA Ya H. Every measured set of strut loads resulted in a different pressure diagram “A,” all of which were found to lie within the boundaries of the trapezoid indicated by the dotted lines. Thus, if strut loads are computed on the basis of this trapezoid, they will most probably be on the safe side.

This, therefore, is the origin of the requirement to use a hinged wall where there were no actual hinges. At the time it was a reasonable solution. As noted earlier, solutions for continuous beams existed but back-figuring the pressures using them would have been a formidable “inverse problem” given the computing power of the day. Doing this, however, raises as many questions as it answers, such as the following:

  • If a continuous beam had been used, would the pressure distribution have been different?
  • What is the relationship between the pressure distribution computed by Terzaghi’s method and what is actually experienced by the wall? Put another way, did Terzaghi’s simplification of the structural situation compromise his distribution? (No doubt some conservatism in the pressure distributions offset that problem.)
  • If the pressure distributions are right for engineering purposes, is it still necessary to use a hinged solution? Especially with beam software, a continuous beam is much simpler to analyse and structurally more representative of the actual sheeting and bracing.

Peck himself was well aware of the limitations of the method; he made the following admission:

It is apparent, therefore, that it is useless to attempt to compute the real distribution of lateral pressure over the sheeting. Of far greater practical importance is the statistical investigation of the variation in strut loads actually measured, in order to determine the maximum loads that may be expected under ordinary construction procedures.

This too raises another question: in developing distributions primarily to determine brace/strut loads, do we compromise the accuracy of determining the maximum moment in the sheeting itself?

Moving Forward

It’s difficult to really know how to answer many of the questions this problem raises. Some suggestions are as follows:

  • It is hoped that there is enough conservatism in these earth pressure distributions to accommodate either method. That’s likely, as inspection of some of Peck (1943) curves will attest. That likeliness is buttressed by the fact that these methods came out of a deadly accident.
  • More comparisons of hinged and continuous beam models are needed. There is one in Sheet Pile Design by Pile Buck and another in the post on this site A Simple Example of Braced Cut Analysis. These are simply not enough to establish a trend one way or another, although the results are interesting and hold promise.
  • A “hand” solution based on parametric studies using FEA or another numerical method would move things forward considerably. Obviously these are limited by the accuracy of the soil modelling but they can be applied to a wider variety of cases.
  • Field tests should include measurement of actual lateral earth pressures on the sheeting at various points. The use of strut loads, although easier to measure with the technology of the 1930’s and 1940’s, is still indirect. Another interesting approach is to use an inverse method and a continuous beam with existing data, although this is not as satisfactory as direct measurement of earth pressures.

References

  • Boone, S.J. & Westland, J. (2005) “Design of excavation support using apparent earth pressure diagrams: consistent design or consistent problem?” Fifth International Symposium on Geotechnical Aspects of Underground Construction in Soft Ground, International Conference on Soil Mechanics and Geotechnical Engineering, 809 – 816.
  • Peck, R.B. (1943) “Earth-Pressure Measurements In Open Cuts, Chicago (Ill.) Subway.” Transactions of the American Society of Civil Engineers, Vol 108, pp 1008-1036.
  • Rogers, D.A. (2013) “Ralph Peck’s Circuitous Path to Professor of Foundation Engineering (1930-48)” Presented at the Seventh International Conference on Case Histories in Geotechnical Engineering.
Posted in Academic Issues, Geotechnical Engineering

Seepage and Bottom Heave Calculations for Sheet Pile Braced Cut Trenches

Part of Soils in Construction‘s presentation of dewatering is this topic. I cover it in my treatment of flow nets in Soil Mechanics: Groundwater and Permeability II; however, Soils in Construction uses a less computationally intensive approach. In this piece I’ll explain that, give some better graphics than the book had available at the time of publication, and compare them with the flow net/FEA results I discussed in my Soil Mechanics course.

Let us consider the problem of a braced cut, which is commonly used for “cut and cover” construction. Such a construction is shown at the right. One of the challenges of temporary works such as this is to insure that a) there is sufficient pumping capacity to keep the “steel trench” dry, and b) the hydraulic gradient of the water coming up into the trench is sufficiently low to avoid soil boiling and bottom heave due to that soil boiling. (Bottom heave can take place due to other factors as well.) An example of that (showing a flow net) is below to the left.

Because the cut is internally braced, it is usually possible to make the sheeting walls simply penetrate to just below the bottom of the trench. However, in order to mitigate the effects of water flowing from around the cut into the bottom, the sheeting can be extended. The idea is that, the longer the extensions, the longer distance the water has to flow, the increased resistance of the soil to flow, and the lower the hydraulic gradients, which both reduce the flow overall and the possibility that the soil at the bottom of the trench will boil, i.e., enter into a quick condition.

Soils in Construction shows two methods of dealing with this problem:

  1. A “rule of thumb” for sheeting extension; and
  2. Charts to determine the minimum extension of the sheeting. These were developed by Marsland (1953). In the book they were taken from NAVFAC DM 7.01, but since the book’s publication NAVFAC DM 7.1 has redone the graphics, and you can see that below (and click on it to download). It is important to note that there is an error in the lower part of this figure; I checked it against Marsland (1953) and have modified it a bit to restore it to Marsland’s original formulation, the following example will show how it is supposed to be used.

Example Problem

I have used this example for many years and feature it in Soil Mechanics: Groundwater and Permeability II. Consider the braced cut shown below.

The parameters for this problem are as follows:

  • Braced Sheet Piling Excavation
    • Depth of excavation = 38’
    • Width of excavation = 32’
    • Impervious layer 20’ below the toe of the sheeting
    • Length of Sheet Piling = 66’
  • Water table at the excavation level on the excavation side
  • Water table 6’ below the top of the sheeting on the soil side
  • Variables for chart above
    • Hw = 38 – 6 = 32′
    • D = 66 – 38 = 28′
    • Hl = 20′
    • H = D + Hl = 28 + 20 = 48′
    • W = 16′
    • Even though the graphic for the problem doesn’t show it, the problem statement indicates an impervious layer below the excavation; thus, we will use the lower chart for this problem.
  • Soil Conditions
    • Uniform medium sand, k = 0.0003 ft/sec
    • Saturated Unit Weight = 115 pcf

We need to, one way or another, insure that the design does not experience soil boiling.

The “rule of thumb” in Soils in Construction states that “the depth of penetration of the cut-off wall below the bottom of the excavation should be a third of the “length” computed. For this wall, the ratio is 28/66 = 42%, which means the rule of thumb is achieved.

Turning to the chart above, we must compute three quantities:

  • The x-axis ratio of half width of excavation to net hydrostatic head, or W/Hw = 16/32 = 0.5;
  • the y-axis ratio of penetration required to net hydrostatic head, or D/Hw = 28/32 = 0.875; and
  • the ratio of the depth between the bottom of the excavation and the impervious layer to the net hydrostatic head, or H/Hw = 48/32 = 1.5.

In this case we only have two values of H/Hw to work with: 1 and 2. For H/Hw = 1, FS ~ 2.5 (this is a very rough extrapolation.) For H/Hw = 2, FS = 1.75. Since H/Hw = 1.5 is in the middle between the two, the FS = (2.5+1.75)/2 = 2.125.

If we want to check our results by neglecting the impervious layer, we use the upper chart, and assuming the soil is closer to being a loose sand, FS ~ 1.4.

So how does this all compare to a flow net/FEA result? Same is given below.

At the right is a flow net generated by the finite element program SEEP-W. At the left is a chart showing the direction of the water flow (arrows) and the hydraulic head (coloured bands.) An in-depth explanation of these can be found at Soil Mechanics: Groundwater and Permeability II. The results are as follows:

  • Gradient at bottom of excavation = 0.48
  • Factor of safety = 1.77
  • Total flow = 1.02 gpm/ft of wall

The critical hydraulic gradient, computed by the methods shown in Computing Pore Water Pressure and Effective Stress in Upward (and Downward) Flow in Soil, is 115/62.4 – 1 = 0.843, which checks with the factor of safety for the actual gradient.

The factor of safety from the chart ignoring the presence of the impervious layer is the closest to the FEA/flow net result. For the chart with the impervious layer, the extrapolation for the lower factor of safety should perhaps be ignored.

A couple of other charts (based on Marsland (1953)) are shown below.

References

  • Marsland, A.R. 1953. “Model Experiments to Study the Influence of Seepage on the Stability of a Sheeted Excavation in Sand.” Geotechnique, 3(6), 223-241.
Posted in Academic Issues, Geotechnical Engineering

The 2:1 Method for Estimating Stresses Under Foundations

Readers of my post Analytical Boussinesq Solutions for Strip, Square and Rectangular Loads and those related to it know that the math related to these methods can get complicated, and in any case the idea of a “purely elastic” soil response to load is purely theoretical. So is there a simpler way? The most common simplification used is the 2:1 Method, shown at the right.

The method basically assumes the following:

  • Uniformly loaded foundation
  • Only stresses of interest are under the centroid of the foundation, as shown, and only vertical stresses are considered
  • Stresses decrease as if there is a truncated pyramid below the foundation with a slope of 1H:2V.

As shown, to determine the additional vertical stress \Delta \sigma_v at a distance z below the centre of the foundation, the equation for a rectangular foundation of width B and length L (can be interchanged, but conventionally B is the smaller of the two) for a load Q on the surface is given by the equation

\Delta \sigma_v = \frac{Q}{(B+z)(L+Z)} (1)

If the equation is written in terms of the unit load q, it becomes

\Delta \sigma_v = \frac{q}{(1+\frac{z}{B})(1+\frac{z}{L})} (2)

Obviously for square foundations B=L but the solution is the same.

As an example, let us consider a foundation where Q=100 kN,\,B=5m,\,L=8m,\,z=3m. By substitution into Equation (1), \sigma_v = 1.14\,kPa. If we compute the unit load to be q = \frac{100}{(5)(8)}=2.5\,kPa, substitution into Equation (2) yields the same result.

Although it is possible to apply this method to circular foundations, it is just as easy to use Boussinesq theory under the centre of the foundation. The equation for a circle of radius r and all other variables the same is

\Delta \sigma_v =  q(1 - \frac{z^3}{(z^2+r^2)^\frac{3}{2}}) (3)

The 2:1 method is the only method taught in Soils in Construction. It is generally used with shallow foundations, but can also be applied to pile groups in clay, as shown at the right. It’s simplicity and reasonable accuracy has brought it acceptance, and further information on it can be found in the Soils and Foundations Reference Manual. The concept has also been applied to pile toes, and you can see this in STADYN Wave Equation Program 10: Effective Hyperbolic Strain-Softened Shear Modulus for Driven Piles in Clay.

Posted in Academic Issues, Geotechnical Engineering

Specific Surfaces of Soil Materials

One concept that appears at the end of Chapter 1 of Soils in Construction is that of the specific surfaces of soil materials, including data from Lambe and Whitman’s classic text Soil Mechanics. Additional information on this topic can be found in Tsytovich, which is reproduced below. (Tsytovich refers to this as “unit surface area,” but the concept is the same.)

A factor of importance in the evaluation of the properties of solid soil particles is their mineral composition. Thus, some minerals, such as quartz and feldspar, interact only slightly with the surrounding water, whereas other minerals, for instance, montmorillonite, can interact substantially more actively and in a different way. The smaller the particles of a soil, the greater their unit surface area (per cm3 or per gram) and the larger the number of centres of interaction with the surrounding water and in contacts between solid particles proper. For instance, particles of kaolin (a clay mineral) have a unit surface area of 10 m2/g, whereas those of montmorillo­nite have a very large unit surface area of 800 m2/g, which inevitably affects the properties of natural soils containing particles of montmorillonite. The presence of particles of mica (which are very slip­pery and have only a negligible shear resistance) has an essential effect on the physical properties of the soils containing such particles; this circumstance must always be taken into account.

Posted in Academic Issues, Geotechnical Engineering

The “Why” and “What” of Soils in Construction

Anyone familiar with the history of geotechnical engineering is aware that its development can, to some extent, be tracked with the development of its textbooks.  Early textbooks tended to be vague and empirical in nature.  With new books more theory is found, especially after the works of Terzaghi, Peck, Tschebotarioff and Taylor.  By the early 1970’s there was a large selection of textbooks for the undergraduate instructor to choose from in the topics of soil mechanics, foundations or books that featured both.

This selection, which peaked with the end of the “Golden Age” of geotechnical engineering, has thinned considerably, as it has with most textbooks.  Soils in Construction made its debut in 1974, right in the middle of the last burst of activity in the field.  So why has this textbook endured when so many others have fallen by the wayside?

The answer is simple: it wasn’t aimed at undergraduate civil engineering students but at construction management ones, and in that respect it was far enough ahead of its time to endure but no so far as to die before its time would come.

When I started my career in the 1970’s, people with geotechnical knowledge working for contractors were a rare breed.  This is not to say that such people had not taken their place; the likes of Lazarus White comes to mind.  But most contractors–especially small and medium size firms, firms which frequently specialised in one or more aspects of geotechnical construction–did not have on staff people with a working knowledge of applied geotechnical theory.

Contractors were not alone in this lack.  State DOT’s were likewise short of people with this type of understanding, and they had large sums of public money entrusted to them.  The FHWA saw the need to address this issue, and the Soils and Foundations Reference Manual was the result of that effort.  Although it can be used in a college setting (I have done this in my Soil Mechanics and Foundation Design and Analysis courses) it takes a lot of work, more work than most academics are used to putting into an undergraduate course.

Soils in Construction is the answer to this dilemma.  It is geared towards construction management students whose mathematical level may not be up to that of their engineering counterparts.  (But…I always told my students that the only calculus they’d get in my courses is if they didn’t brush their teeth; even for them it is the nature of basic geotechnical engineering.)  It enables them to grasp the basics of the application of soil mechanics to practical problems, including temporary works, whose engineering is frequently overlooked but which is often vital for the successful completion of the permanent works to follow.

With this the authors commend this work to our readers, hoping that it will result in more successful geotechnical projects for contractor, owner and engineer alike.