Posted in Uncategorized

Featured in Print: Materials Testing

Materials Testing CoverMost Soil Mechanics courses feature a lab course to go with them, either integral or separate.  The textbook manufacturers are quick to take advantage of that with expensive books and software.  ASTM specs are just about standard in the U.S., but trying to teach out of an ASTM spec is a trying business.

An alternative to this is Materials Testing, the print version of the U.S. Army’s FM 5-472.  In addition to the soils tests, it features many other items, some of which are not geotechnical tests at all:

  1. Materials Testing Overview
  2. Soils Testing
  3. Bituminous Mixtures Testing
  4. Concrete Testing
  5. Soil Stabilization

The soils tests in particular are as follows:

  • Moisture-Content Determination
  • Specific-Gravity-of-Solids Determination (ASTM D 854-92)
  • Grain-Size Analysis and Distribution (ASTM D 422-63 and ASTM 2217-85)
  • Liquid Limit, Plastic Limit, and Plasticity Index Determination (ASTM D 4318-95a)
  • Laboratory Compaction Characteristics of Soil Using Modified Effort (Compaction Test) (ASTM D 1557-91)
  • In-Place Density Determination
  • CBR Tests

You can view the table of contents below.

Pages from FM5-472_Page_1Pages from FM5-472_Page_2Pages from FM5-472_Page_3Pages from FM5-472_Page_4Pages from FM5-472_Page_5Pages from FM5-472_Page_6Pages from FM5-472_Page_7Pages from FM5-472_Page_8

Materials Testing is an excellent lab manual, supplementary textbook, or handy resource for the testing and improvement of soils, asphalt or concrete.


Posted in STADYN

Improved Methods for Forward and Inverse Solution of the Wave Equation for Piles

This slideshow requires JavaScript.

This doctoral research project at the University of Tennessee at Chattanooga is now complete, and we are pleased to present the following:


This dissertation discusses the development of an improved method for the static and dynamic analysis of driven piles for both forward and inverse solutions. Wave propagation in piles, which is the result of pile head (or toe) impact and the distributed mass and elasticity of the pile, was analysed in two ways: forward (the hammer is modelled and the pile response and capacity for a certain blow count is estimated) or inverse (the force-time and velocity-or displacement-time history from driving data is used to estimate the pile capacity.) The finite element routine developed was a three dimensional model of the hammer, pile and soil system using the Mohr-Coulomb failure criterion, Newmark’s method for the dynamic solution and a modified Newton method for the static solution. Soil properties were aggregated to simplify data entry and analysis. The three-dimensional model allowed for more accurate modelling of the various parts of the system and phenomena that are not well addressed with current one-dimensional methods, including bending effects in the cap and shaft response of tapered piles. Soil layering was flexible and could either follow the grid generation or be manually input. The forward method could either model the hammer explicitly or use a given force-time history, analysing the pile response. The inverse method used an optimization technique to determine the aggregated soil properties of a given layering scheme. In both cases the static axial capacity of the pile was estimated using the same finite element model as the dynamic method and incrementally loaded. The results were then analysed using accepted load test interpretation criteria. The model was run in test cases against current methods to verify its features, one of which was based on actual field data using current techniques for both data acquisition and analysis, with reasonable correlation of the results. The routine was standalone and did not require additional code to use.

Dynamic Pile Testing Results
Crescent Foundation Demonstration Test Pile – Vulcan SC9 Hammer
Kenner, Louisiana

Brian Mondello and Sean Killingsworth
May 2014

This report presents the results from dynamic pile testing, and related data analysis, performed during the initial drive testing of the subject Test Pile on April 30, 2014, at the above referenced jobsite location in Kenner, Louisiana. The primary test objective was the monitoring of the hammer/driving system performance. Additionally, the testing objectives included the monitoring of dynamic pile driving stresses, pile structural integrity, and pile static bearing capacity. These objectives were met by means of a Pile Driving Analyzer® (PDA), Model PAX, which uses the Case Method for numerical computations. An additional analysis was performed on a selected test record using the CAPWAP® computer program. Discussions on the testing equipment, analytical procedures, theory, application, and limitations are presented in Appendix A. Testing and analysis results are presented in Appendix B.

A video of the SC-9 hammer featured in Mondello and Killingsworth:

Dynamic Response of Footings and Piles

Wing Tai Peter To
PhD Dissertation
University of Manchester
February 1985

Dynamic response analyses can be regarded as stress wave propagation problems. The solution of such by the finite element method entails more consideration than static problems, since sources of inaccuracies such as dispersion, spurious oscillations due to mesh gradation, wave reflection at transmitting boundaries, as well as instability or inaccuracy due to temporal operators and discretisation can arise. The criteria for formulating a finite element model for dynamic response analysis have been investigated. Using the relatively simple von Mises soil model (satisfactory for undrained saturated clay) three categories of problems have been investigated:

  1. The dynamic response analyses of surface footings subjected to periodic and impact loading have been performed in order to evaluate the finite element model design criteria. An approximate analysis is also performed in reducing a three-dimensional indirect impact problem to a two-dimensional analysis.
  2. Vibratory pile driving is a relatively new but somewhat unreliable technique of pile installation. Penetration is instantaneous if conditions are right, but with the high hire charges and uncertainty in success the technique is unpopular, especially in clays. In the work presented it is shown that vibratory installation is possible in cohesive soils at the fundamental frequency for vertical pile translation, if a high enough dynamic oscillatory force is provided. Penetration mechanisms have also been exploited.
  3. On the other hand, impact pile driving is reliable and widely adopted in terrestrial as well as offshore construction. Experience in one dimensional wave equation analysis is discussed, and further numerical evaluation of the parameters involved has been carried out by a more elaborate axisymmetric finite element model. In cohesive soils a closed-ended pile may be driven more easily than an equivalent open ended pile, depending on the level of the internal soil column and the soil properties. In the light of the growing popularity of nondestructive determination of the axial load-carrying capacity of piles by dynamic methods, the possibility of correlating the soil resistance mobilised in dynamic conditions to the ultimate static capacity is queried. The semi-empirical Case method has been assessed in detail.
Posted in STADYN

Lessons Learned From Writing Finite Element Code for Pile Dynamics

Back in 1936, Lazarus White (of Spencer, White and Prentis fame) wrote the following about the dynamic formulae:

I read some papers last night where some of these pile driving formulas were derived, and the result was that my sleep was very much disturbed.

While there are many reasons to lose sleep while writing a Ph.D. dissertation, dynamic formulae aren’t among them. However, in the course of researching, coding, debugging, and writing what ended up as Warrington (2016), it became clear that there are things that induce sleep loss in geotechnical engineering in general and pile dynamics in particular.

As the title indicates, this is intended to be a “lessons learned” type of presentation. Given the manifold uncertainties that accompany any geotechnical engineering topic, the dual temptations to either completely discard any past experience or to dogmatically defend the current “state of the art” should be avoided. Given the environment we work in, coming to “definitive” solutions is not an easy task. The industry’s/profession’s innate conservatism, especially in codifying methods and materials, is well founded. On the other hand, there comes a point where serious issues “out there” need to be discussed and new (or at least different) solutions need to be considered.

The current state of the technology of pile dynamics, which in practical terms begins with Smith (1960) and Rausche (1970), has advanced not only the design and installation of driven piles but, though high-speed dynamic testing applied to bored piles and pile integrity testing, deep foundations in general. Since both (and the improvements that have come since then) are numerical methods with computer implementation, to say that computer power has advanced in the intervening years states the obvious. Numerical methods that existed even then, such as those documented in Carnahan, Luther and Wilkes (1969) can be brought to bear on the problem in new ways. Plasticity in finite element analysis is well established, documented in works such as Owen and Hinton (1980). Optimization techniques such as those described in Gill, Murray and Wright (1981) can be used for signal matching techniques, which are crucial in determining the soil resistance to driving (SRD) from dynamic pile data.

So without further introduction, on to the lessons learned.

We Have Gone As Far with the Smith Model as We Can Go

Given the state of soil dynamics when Smith (1960) was published, it is remarkable that his model has survived and continued (with modifications) as the “standard” model for soil response in the one-dimensional wave equation for piles. It has been tested (although really meaningful testing is difficult) and has been used in many applications, both in forward (predictive) and inverse (interpretive) application of the wave equation to piles. It is relatively simple to implement and the quantities of quake and damping have been extensively correlated.

However, any numerical model of a physical phenomenon has weaknesses, and the Smith model is no exception. Some of these are as follows:

There is no really accepted correlation between the soil properties as determined using standard testing procedures and the Smith parameters. This is especially critical with the soil damping, which has been shown to be one of the most critical properties to properly quantify in a wave equation analysis (Meseck (1985)).

The Smith model in the forward mode virtually requires the use of static methods to determine the SRD, which have numerous problems of their own. (Their proliferation is reminiscent of that of the dynamic formulae.) An additional complication is the diffuse state of static load test interpretation. While most dynamic methods in the U.S. presume the use of Davisson’s Method, this method is not universal in application.

The nature of Smith damping has never been completely enunciated. When velocity-based damping is used in the modelling of a physical system, generally speaking it refers to some kind of viscous friction, such as in a “dashpot.” Most of the “damping” in a dynamic pile system is the radiation of energy into the distributed mass and elasticity (and plasticity) of the soil surrounding the pile.

There is no good way to integrate many phenomena experienced in pile driving into the Smith model. Perhaps the most important of these is pile set-up, which is currently the “hot” research topic. In most cases the goal of these studies is to determine the ratio of the SRD (used in the wave equation analysis) to the ultimate static capacity of the pile. What is needed is a method to predict pile set-up, principally through the computation of elevated pore water pressures.

The Smith model’s neat division of soil resistance into shaft and toe does not reflect the reality of many driven piles, especially tapered piles, piles with plugging and piles with non-planar toes.

There are many “add-ons” that we use in the application of the Smith model (and other purely 1D wave equation analyses. But many of these are just that: add-ons to a model that in itself does not take into account many of the phenomena experienced in pile driving. Like Vergil in the Divine Comedy, the Smith model has been a faithful guide, but now we must move on if we want to really advance the science of pile dynamics.

Mohr-Coulomb Isn’t Perfect, But For The Moment It’s The Best We Have

If you realize that I’m advocating for a “3D” model of the hammer-pile-soil system, you’re right. (If we stick with round, axisymmetric piles, we can use a 2D model.) But that raises as many questions as it answers. Some of these relate to the hammer, including the gas pressures in diesel hammers, the modelling of intextensible places such as the hammer-cushion, hammer-cap or cap-pile interfaces, and cushion plasticity. But the key question, as always, relates to the soil. What soil model will we use? How will we take into consideration the inherently non-linear nature of the soil? How will we model the inevitable plasticity? What soil properties will we use?

At this point, in spite of its limitations, the Mohr-Coulomb model, familiar to just about everyone in the geotechnical industry, is probably the best overall soil model to use in pile dynamics. It certainly has its limitations; it is, in some ways, an attempt to render an elasto-plastic simplification of a hyperbolic (Duncan and Chang (1970)) soil response. In that respect it justifies the characterization of Massarsch (1983) as “crude.” But it remains overall the best model to use (Abbo (2011); McCarron (2013)) for the following reasons:

It is suitable for a wide variety of soils, and piles are driven into diverse stratigraphies, both from one site to the next and along the length of the pile.

The current standard testing methodology is “designed” around Mohr-Coulomb. If we are serious about addressing the lack of correlation between properties used in pile dynamics and soil properties, we need to be able to test for these properties. Until we revise our standard testing methodologies, this limitation will stand.

It is well represented in finite element codes that are used for geotechnical analysis. Warrington (2016) represented the first application of a standalone finite element code for pile dynamics in over twenty years. There are advantages and disadvantages in both approaches, but implementation of Mohr-Coulomb in finite element code is well established.

For these and other reasons, Mohr-Coulomb remains the best model for failure and plasticity that we have available. Like everything else, this is not written in stone; but until other things advance in the industry, this is the best we have.

We Need To Understand Why Reducing Pile Head Data Is So Difficult

The ability to estimate static pile response from high-velocity dynamic data is the cornerstone of pile dynamics. Although it is “intuitively obvious” (a favourite phrase of academics,) quantifying it has been something of a challenge. Today methods in use for doing this estimation are basically signal-matching techniques. The concept is to have the model return the same force-time and velocity-time histories as come from the field. In processing the data, we compare the two “point by point;” if the sum of the differences (or the square of the differences) goes below some kind of tolerance criterion, we say the signals match and the model replicates the physical system.

Although the results have been reasonable problems remain, requiring a great deal of human intervention. So why is this problem so difficult? The answer, in terms of optimization techniques, is the presence of multiple local minima of the aforementioned “sums.” This is a “fancy” term that needs some explanation.

Optimization techniques attempt to find a minimum (or maximum) value for a function of several variables. Consider the example of the famous Rosenbrock Function, given (in a typical form) by the equation

R(x,y) = (1-x)2 + 100(y-x2)2

A plot of this is shown below.


This function is often used to test optimization techniques, by starting at a point which away from the known minimum and checking to see how well the technique converges on the actual minimum, which in this case is . Although the function is difficult, it is reasonably “regular.”

Now consider a topographical map of a hilly region, or better the soil profile for a Karst topography. Now there are many local minima; it is easy for an optimization technique looking for the “bottom” of the valley without realizing that there is yet another valley over the hill that is actually lower! Using an automated technique, especially one based on Newton’s Method, can easily get trapped in such a minimum, thinking it has found the solution to the problem when in fact it has not.

Getting to the real minimum has been a central frustration in the development of signal matching techniques for pile dynamics, because the inverse problem in pile dynamics, like the ground the piles are driven into, has many local minima. This problem has impeded the automation of the search process more than any other. The increase in computer power available to engineers is an opportunity to “get to the bottom” of the problem, although in doing so a consistently successful methodology is necessary.

We Need to Instrument Piles in More Places than the Pile Head

The earliest instrumented pile study was that of Glanville (1938), and the piles were instrumented at the head, middle and toe. Jorg Osterberg showed us how to test the piles with a load cell at the toe, and Fellenius (2014) has advocated instrumenting piles at locations other than the head for static testing. And yet today, eighty years after it was first done, we generally still instrument the pile head only for both static and dynamic testing. The need to do so for dynamic testing is, if anything, more urgent than with static testing, and with the expendable instrumentation we have available it should not be that big of a problem. Doing so would give us a more complete picture of what’s really going on during driving and make the search for an SRD much simpler.

The Advance of Numerical Methods Should Occasion Some Changes in Geotechnical Education

At this point we turn away somewhat from the problem at hand and look at geotechnical education in general. If we take a serious look at our basic geotechnical courses, much of the material we have is little different than that which was published in Terzaghi and Peck (1948). While we cannot denigrate those who went before us just because we have found out something new, nevertheless the environment we operate in is vastly different from that which those two giants strode in a half century and more ago. The biggest change is the introduction and proliferation of numerical methods.

Many geotechnical engineers are wary of the advance of these methods. Some will tell you that the results obtained from hand calculations based on experience are no worse (and many cases better) than those from numerical methods. Often this is correct. Part of the problem is that the numerical methods are based on the same theory as the closed form solutions. The advantage of numerical methods is not necessarily in underlying theory but in the ability to better simulate the physical system through geometric replication. But another part of the problem is that most numerical methods, presented in packages, are “black boxes” to many practitioners, who frequently do not have a complete grasp of their inner workings. The temptation is thus great to receive the results uncritically, without support from more conventional methods or experience. It is thus incumbent upon geotechnical educators, starting at the undergraduate level, to give students some understanding of how these methods work, and also in their strengths and weaknesses, so that they can form an intelligent opinion of the results.

The biggest difference between the way geotechnical professionals have employed theory in the past and the way numerical methods handle it is that, in the first case, elasticity and plasticity are handled separately, while in the second the two are simulated in the same model, which alternates from one to the other depending upon the system modelled and the loads placed upon it. For example, in settlement problems we frequently consider elastic settlement separate from consolidation settlement (primary and secondary,) which in turn is separate from bearing capacity. In physical reality all of these are taking place at the same time, and the model that can accurately simulate this is of value.

It is for this reason that this educator has shifted towards a text such as Verruijt and van Bars (2007). Such an approach emphasizes more of the “continuum mechanics” of soil mechanics than is customarily done. This is arguably a more “theoretical” approach to soil mechanics. It is not a complete solution to the problem, but it is a start.


These are some observations and proposals derived from a long (three year) study on the subject of pile dynamics. Geotechnical engineering tends to be a conservative field of study and practice, but we must not let that impede real progress towards foundations that are more economical and reliable than before. Ultimately numerical methods in general and finite element analysis in particular will advance in usage; if we implement them intelligently and educate our new engineers in their inner workings, the profession will advance.

If we can achieve that goal, we can all sleep better at night.


  • Abbo, A.J., Lyamin, A.V., Sloan, S.W. and Hambleton, J.P. (2011). “A C2 continuous approximation to the Mohr-Coulomb yield surface.” International Journal of Solids and Structures, 48(21), 3001-3010.
  • Carnahan, B., Luther, H.A. and Wilkes, J.O. (1969). Applied Numerical Methods, John Wiley & Sons, Inc., New York, NY
  • Duncan, J.M., and Chang, C.Y. (1970). “Nonlinear Analysis of Stress and Strain in Soils.” Journal of the Soil Mechanics and Foundations Division, American Society of Civil Engineers, 96(5), 1629-1653.
  • Fellenius, B.H. (2014) Basics of Foundation Design. (March 13, 2015)
  • Gill, P.E., Murray, W. and Wright, M.H. (1981). Practical Optimization, Academic Press, Inc., London, England.
  • Glanville, W.H., Grime, G., Fox, E.N, and Davies, W.W (1938). An Investigation of the Stresses in Reinforced Concrete Piles During Driving, Department of Scientific and Industrial Research, British Building Research Board, London, England.
  • Massarsch, K. R. (1983). “Vibration Problems in Soft Soils.” Proceedings of the Symposium on Recent Developments in Laboratory and Field Tests and Analysis of Geotechnical Problems, Asian Institute of Technology, A. A. Balkema, Rotterdam, 539-549.
  • McCarron, W.O. (2013). “Numerical Modeling Applications in Offshore Petroleum Developments.” Geo-Strata, 17(3), 34-39.
  • Meseck, H. (1985). “Application of a Wave Equation Programme to Establish the Bearing Capacity of Driven Piles.” Proceedings of the International Symposium on Penetrability and Drivability of Piles, Volume 2, Japanese Society of Soil Mechanics and Foundation Engineering, 84-90.
  • Owen, D.R.J., and Hinton, E. (1980). Finite Elements in Plasticity: Theory and Practice. Pineridge Press, Swansea, Wales.
  • Rausche, F. (1970). Soil Response from Dynamic Analysis and Measurements on Piles. Ph.D. Dissertation, Case Western Reserve University, Cleveland, OH.
  • Smith, E.A.L. (1960). “Pile-Driving Analysis by the Wave Equation.” Journal of the Soil Mechanics and Foundations Division, 127(1) 35-61.
  • Terzaghi, K., and Peck, R.B. (1948) Soil Mechanics in Engineering Practice, McGraw-Hill Company, New York, NY.
  • Verruijt, A., and van Bars, S. (2007). Soil Mechanics. VSSD, Delft, the Netherlands.
  • Warrington, D.C. (2016). Improved Methods For Forward And Inverse Solution Of The Wave Equation For Piles. Ph.D. Dissertation, University of Tennessee at Chattanooga, Chattanooga, TN.
Posted in Uncategorized

Two Deep Foundations Programs–BENT1 and PX4C3–Now Available

This week we feature two programs for the analysis of deep foundations: BENT1 and PX4C3.  Although these names may not be familiar, if you’re in the deep foundation engineering field you’ll probably recognise some of their successors.

Let’s start with BENT1, the program description is as follows:

BENT1 is a program designed to analyse two dimensional problems involving pile supported foundations subjected to inclined & eccentric loadings. It consists of an interactive solution for the three equilibrium equations. The purpose of the iterative procedure is to find deflected position of the structure so that equilibrium and compatibility are satisfied. The pile cap is assumed to be rigid in the analysis. Subroutines MAKE and COM62 are used in this program. A complete description of the program can be found in the document “Background Theory and Documentation of Five University of Texas Soil-Structure Interaction Computer Programs,” Miscellaneous Paper K-75-2, by N. Radhakrishnan and F. Parker.

COM62 is the ancestor of all of the programs with that in the name, and ultimately programs such as LPILE.

Then we consider PX4C3:

PX4C3 is a finite difference program used to compute load settlement characteristics on an axially loaded pile of constant outside diameter. A set of load transfer curves along the pile (i.e., skin friction developed on the side of the pile relative to the absolute axial displacement of the pile section) & four point resistance curve at the pile tip (i.e., relationship between the total axial soil resistance on the base of the pile tip & the pile tip movement) are used in program to obtain non-linear soil-pile relationships. Finite difference equations are used to achieve compatibility between pile displacement & load transfer along the pile & between soil resistance & load transfer along the tip of the pile. A complete description of the program can be found in the document “Background Theory and Documentation of Five University of Texas Soil-Structure Interaction Computer Programs,” Miscellaneous Paper K-75-2, by N. Radhakrishnan and F. Parker.

This is the ancestor of programs such as APILE.

Complete documentation for both programs is included, and can also be found here.

Download BENT1 and PX4C3

Posted in Uncategorized

Soils in Construction, Now at Waveland

It’s official: the construction management textbook Soils in Construction, Fifth Edition by W.L. Schroeder, S.E. Dickenson and D.C. Warrington is now at Waveland Press.  As someone who has dealt with contractors for his entire working career, I know that an understanding of the essentials of soil mechanics and foundations is crucial for successful–and profitable–completion of just about any civil engineering project.  Most soil mechanics and foundations textbooks are aimed at design professionals or engineering students; Soils in Construction, Fifth Edition is aimed at the contractor and the construction management student. It presents soil mechanics concepts in an easy to understand fashion. It includes many worked examples and illustrations to make understanding of the basic concepts simple.

And the really good news is that, with Waveland as publisher, the price has been reduced considerably.

Chapters in the book are as follows:

  1. Physical Character of Soil Constituents
  2. Natural Soil Deposits
  3. Soil Index Properties
  4. Soil Classification
  5. Stress Analysis and Engineering Properties
  6. The Contract and Contract Documents
  7. Interpretation of Soils Reports
  8. Embankment and Construction Control
  9. Dewatering
  10. Excavations and Excavation Supports
  11. Foundation Construction
  12. Construction Access and Haul Roads

The book also has two appendices:

  • Laboratory Testing Exercises
  • Pile Hammer Specifications

You can order a copy of the book by clicking on the graphic above or here.  Academics who would like more information or order a desk copy can do that here.