Posted in Uncategorized

Analysing a Gravity Wall Using Vector Analysis

One of the reasons I was interested in teaching Statics at Lee University was because I was continually disappointed at my students’ memory of their statics. Statics is crucial in the design and analysis of geotechnical structures, and most of the problems–at the undergraduate level at least–aren’t that involved, or at least I thought they weren’t. A great deal of the problem is that geotechnical statics usually involves converting distributed loads into resultants, which Statics–and Mechanics of Materials for that matter–generally associate this with beam problems, not always the case with geotechnical problems.

Another culprit is that Statics, in the U.S. at least, is a vector proposition from the start. At the University of Tennessee at Chattanooga where I taught, it was called “Vector Statics,” which gives the game away early. (At Lee we use the same book and teach the same material, but simply title it “Statics.”)

But what if we applied a vector approach to a simple geotechnical problem? That’s what we’re going to do here with a concrete gravity wall. I will use the method outlined in my post A Simplified Method to Design Cantilever Gravity Walls. You can refer to the theory there, I will try to keep it to a minimum. The wall is pictured at the top of the post, I will reproduce it below.

Analysing Overturning

We have three forces acting on the wall:

  1. The weight of the gravity wall itself, Wconc
  2. The weight of the soil trapped by the heel of the wall, Wsoil
  3. The lateral force of the soil on the wall, Fh

Forces 1 and 2 are determined by computing the cross-sectional area of the concrete and soil and multiplying each by the unit weight as shown above, and then converting the result to a vector force and placing it at the centroid of the area (another Statics topic.) Instead of the “manual” approach in A Simplified Method to Design Cantilever Gravity Walls, the was was drawn in CAD and both the areas and centroids were determined automatically. You can see the magnitudes and locations of those resultants above.

The lateral force of the soil is computed using Rankine’s theory. The first thing is to determine the working internal friction angle of the soil by applying the Shear Mobilisation Factor SMF. Assuming an SMF of 2/3, that friction angle changes from the 30 degree one shown above to a 21.05 degree one, which is applied to the formula for Rankine active pressures for level backfill,

k_h = \frac{1-sin\phi}{1+sin\phi} (1)

Doing that results in a kh = .471. The force on the wall is then determined by the formula

F_h = \frac{k_h\gamma_{soil}H^2}{2} (2)

The division by two reflects the fact that soil effective stress (and thus lateral earth pressure) increases linearly with depth (like a fluid,) creating a triangular distribution (yet another concept from Statics.)

At this point there the resisting forces R and T are not defined. The forces themselves are easily computed by summing forces in the x and y directions. Doing this, we have

R = -W_{conc}-W_{soil} (3a)
T = -F_{h} (3b)

The location of F–along the surface of the footing–is evident. The location of R is not; it is some distance x from the toe (Point “A”) of the footing. We can obtain x by summing moments around Point “A,” and with a vector method that means taking cross products of the moment arms with the forces.

Converting both the moment arms r and the forces to vector notation yields the following:

  • Concrete Weight: r = 2.358 i + 4.049 j, Wcon = −4.35 j
  • Soil Weight: r = 5.576 i + 8.264 j, Wsoil = −6.30 j
  • Lateral Earth Pressure: r = 8 i + 4 j, Fh = −4.073261616 i
  • Vertical Footing Force: r = x i, R = 10.65 j (Equation (3a))

The force T does not enter into this because its line of action runs through Point “A,” thus its moment is zero as its moment arm is zero.

The cross product moments around the toe (Point “A” in the drawing) are as follows:

  • Concrete Weight: \left [\begin {array}{ccc} i&j&k\\{\medskip} 2.358& 4.049&0 \\{\medskip}0&- 4.35&0\end {array}\right ]  = -10.257 k
  • Soil Weight: \left [\begin {array}{ccc} i&j&k\\{\medskip} 5.576& 8.264&0\\{\medskip}0&- 6.3&0\end {array}\right ] = -35.129k
  • Lateral Earth Pressure: \left [\begin {array}{ccc} i&j&k\\{\medskip}8&4&0\\{\medskip}- 4.073&0&0\end {array}\right ] = 16.293k
  • Vertical Footing Force: \left [\begin {array}{ccc} i&j&k\\{\medskip}x&0&0\\{\medskip}0& 10.65&0\end {array}\right ] = 10.65 x k

Summing these moments,

29.093\,k+ 10.65\,xk=0 (4)

Solving yields x = 2.732′. At this point we need to determine whether this is an acceptable location or not for the force. The goal is for the pressure to be positive (downward) along the entire surface of the footing. There are two ways of determining this:

We will do the latter. The middle third of this foundation falls between 2.67′ < x < 5.33′, so the vertical footing force is within the middle third (barely.) As I noted in A Simplified Method to Design Cantilever Gravity Walls, “In this case we make a common assumption that, as long as the resultant force of the wall is within the kern and there are no negative pressures on the base, overturning will not be experienced. It is certainly possible to do an explicit overturning analysis to check this result.”

Analysing Sliding

With the lack of keys or deep foundations, the only lateral resistance to sliding is the friction force T. We computed that force based on Equation (3b,) but in reality that force cannot exceed–and there should be a factor of safety in that inequality–the frictional force possible, which is defined by the equation

T_{max} = \mu R (5)

in which case

FS_{Sliding} = \frac{T_{max}}{T} (6)

Equation (5) is written in “mechanical engineers format.” Geotechnical engineers understand all too well the concept of a friction angle. In my post Explaining the Relationship Between the Coefficient and the Angle of Friction I relate the two from a non-geotech standpoint; we can turn Equation (5) into a more “geotech-friendly” form by noting that

\mu = \tan\phi (7)

Let us assume that the value of \phi is the same under the wall as next to the wall, and let us also assume that the friction angle between the base and the soil is the same as the friction angle of the soil overall, as was done in A Simplified Method to Design Cantilever Gravity Walls. That being the case, \mu = \tan(30) = 0.577 . Substituting into Equation (5,) T_{max} = 0.577 \times 10.65 = 6.15\,kips . The factor of safety from Equation (6) is thus FS_{Sliding} = \frac{6.15}{4.073} = 1.51, which is barely over the minimum criterion for usual loads given in A Simplified Method to Design Cantilever Gravity Walls.

Observations

  • The use of vectors for this problem is overkill from a computational standpoint. It also requires locating the centroid/CG of the two regions in both the x- and y-directions, although with using CAD this is trivial. On the other hand doing it using vectors is more “bullet proof” in that the student is not required to “think” but just “plug and chug” without having to identify lines of action and perpendicular moment arms.
  • The fact that the word “barely” appears in both analyses should inspire some additional conservatism in the design. The simplest way to improve the situation would be to move the heel to the right, which would shift the resisting forces away from the toe (and thus increase their resisting moment) and also put the footing force resultant deeper into the middle third.
  • Both bearing capacity and settlement of the wall’s foundation, the methodology for which are discussed in A Simplified Method to Design Cantilever Gravity Walls, are beyond the scope of this post. Also beyond the scope of this post is the structural design of the wall and of course the global stability of the wall as well.

Posted in Academic Issues, Geotechnical Engineering

Closing the Loop (or at least trying to) on Hough’s Settlement Method

My post Getting to the Legacy of B.K. Hough and his Settlement Method got a good amount of interest. Unfortunately it left several “loose ends,” some of which is because we don’t have full information on some of Hough’s methodology, and others because we don’t have full information on the FHWA’s thinking in their resurrection of the method. The objective of this post is to clear up some of these loose ends, especially the latter.

What we needed was some additional information on the background for both Hough’s work and the FHWA’s interpretation of that work. That information came from an unexpected source: the FHWA’s Design and Construction of Driven Pile Foundations, 2016 Edition. This has been available on vulcanhammer.info just about since it was first published. Now some of you are asking, “Why is this guy waiting until now to come out with stuff he’s had for years? Doesn’t he read the material he offers for download?” In my defence, I didn’t expect the solution to come from the driven pile manual, since it’s primarily a shallow foundation issue. It’s there because it’s proposed for use with groups of driven piles in cohesionless soils.

With all that said let’s look at some of the issues this “new” source tackles.

The SPT Correction Issue

The basic equation for Hough’s method is similar to that used for one-dimensional consolidation settlement and is

S=\frac{H_{o}}{C'} \log_{10} \frac{\sigma'_f}{\sigma'_o} (1)

where

  • S = settlement of cohesionless soil layer
  • H_o = layer thickness
  • C' = compression coefficient
  • \sigma'_o = effective stress at the centre of the layer
  • \sigma'_f = effective stress at the centre of the layer plus the induced stress from the surface at the centre of the layer

The tricky part has always been C' . My criticisms of the FHWA’s implementation of the method in the post Getting to the Legacy of B.K. Hough and his Settlement Method were as follows:

  • Hough’s Method may or may not have used a 60% (N60) efficient hammer to develop the method. It’s entirely possible but Hough (1959) doesn’t say what type of hammer he used to develop the method.
  • The FHWA implementation of Hough’s method included an SPT correction for overburden, which was absent in Hough (1959.)
  • An absence of any mention of Hough (1969,) where he modified the C' values.
  • The lack of an attempt to convert the chart to equations. My original post demonstrates how this can be done, using the method of Hough (1969.)

Design and Construction of Driven Pile Foundations, 2016 Edition addresses all of this as follows:

Cheney and Chassie (2002) report that FHWA experience with this method indicates the method is usually conservative and can overestimate settlements by a factor of 2. This conservatism is attributed to the use of the original bearing capacity index chart from Hough (1959) which was based upon SPT donut hammer data. Based upon average energy variations between SPT donut, safety, and automatic hammers reported in technical literature, Figure 7-43 now includes a correlation between SPT N values from safety and automatic hammers and bearing capacity index. The safety hammer values are considered N60 values. This modification should improve the accuracy of settlement estimates with this method.

The following should be noted about these changes:

  • The assumption that Hough used a “donut” hammer is a reasonable one based on the technology of the time, but it’s still an assumption. Hough doesn’t tell us the kind of SPT hammer he used, or even how he came up with the C’ values shown above.
  • The chart above (Figure 7-43) shows the correlation for three types of hammers: donut, safety, and automatic hammers (which now “rule the roost” in testing.) However, it still insist that the safety hammer values should be corrected for overburden, something else absent from Hough’s study.
  • There is still no awareness of Hough (1969.)

So we can say that we have, perhaps, made some progress toward a solution, but at this point we are not quite where we would like to be.

My Thinking on the Way Forward

We usually think in terms of progress in this field in terms of peer-reviewed articles. But we’ve had the peer-reviewed articles, the experts weigh in on what they mean, and another set of experts try to make things better. My own solution to this problem would run like this:

  • Let’s accept the FHWA’s idea that Hough used a donut hammer for his original work. Donut hammers have a “standard” efficiency of 45%. To get the SPT blow counts to an N60 value (60% efficiency) we need to divide the donut SPT values by 60/45 = 4/3. (That’s what’s going on in the graph above.)
  • Let’s use Hough’s method as he presented it in Hough (1969) and assume that those values too came from donut hammers.
  • Let’s lose the overburden correction; it wasn’t in the original and I don’t see how one can justify putting it into this method.

If we implement all of this, the chart now looks like this:

If formulae are preferred (and that’s normally the case these days) they will again be in the form

C' = A \exp^{BN} (2)

and the coefficients will be as follows:

Soil TypeAB
Organic Silt, Little Clay7.220.0305
Inorganic Sandy Silt18.270.0279
Very Well Graded Fine to Coarse Sand22.860.0270
Well Graded Clean Fine to Coarse Sand28.220.0289
Well Graded Silty Sand and Gravel32.850.0289
Uniform Clean Inorganic Silt37.020.0294
Very Uniform Clean Medium Sand (Similar to Std. Ottawa)58.660.0299

The “B” coefficients run in a fairly narrow range and have an average of 0.0289. It’s possible using this or a more sophisticated method to apply the same value of B to all of the soils.

Any comments from those of you who have used Hough’s Method, or have research on the topic, would be greatly appreciated.

References

  • Hough, B.K. (1959). “Compressibilty as the Basis for Soil Bearing Value,” Journal of the Soil Mechanics and Foundations Division, ASCE, Vol. 85, Part 2.
  • Hough, B.K. (1969). Basic Soils Engineering. Second Edition. New York: Ronald Press Company.