Using Piezocone Dissipation Test to Estimate the Undrained Shear Strength (su), Nkt and NΔu Factors in Cohesive Soils

Show more

1. Introduction

Prediction of soil properties from piezocone test (CPTu) data in cohesive soils is routinely carried out in geotechnical design. This is possible because there is a general recognition that analytical and numerical analysis techniques and constitutive models of soil behaviour are now sufficiently developed to produce good agreement between cone parameters and independently measured soil properties (e.g. Lunne et al., 1997 ; Yu, 2004 ; Schnaid, 2009 ). A theoretical frame is essential, because the CPTu cannot measure the undrained shear strength directly and therefore CPTu assessment of s_{u} relies on a combination of theory and empirical correlations.

Penetration tests in clay are generally undrained, and therefore the cone tip resistance q_{c} can be related to s_{u} as follows:

${q}_{t}={N}_{kt}\cdot {s}_{u}+{\sigma}_{0}$ (1)

where N_{kt} is a theoretical cone factor and σ_{0} is the in situ total stress. The cone factor may be determined using simple bearing capacity formulations, cavity expansion or strain path method (e.g. Terzaghi, 1943 ; Meyerhof, 1956 ; Caquot & Kerisel, 1953 ; Baligh, 1985 ; Teh & Houlsby, 1991 ; Yu & Whittle, 1999 ). Yu (2004) pointed out that while each theory may be used alone for cone penetration analysis, better predictions of cone penetration mechanisms may be achieved if some of the methods are used in combination. A combination of strain path analysis and finite element calculations was used by Teh & Houlsby (1991) to model cone penetration in a Von Mises soil. Yu & Whittle (1999) proposed a cone factor estimated from both strain path analysis and cavity expansion methods. In this approach, the strain path solution developed by Baligh (1985) was used to estimate the size of the plastic zone produced by penetration. Once the plastic zone is established, spherical cavity expansion was used to determine the stress distribution and therefore cone resistance. Burns & Mayne (1998) use cavity expansion-critical state framework to model monotonic and dilatory response with regard to time.

Whereas theoretical solutions have been contributing in the understanding of the fundamental mechanics of cone penetration, empirical correlations are still widely used in practice to estimate s_{u} from cone resistance. Values of cone factor (N_{kt}) often fall in the range from 10 to 20 and are influenced by soil plasticity, overconsolidation ratio, sample disturbance, strain rate and scale effects, as well as the reference test from which s_{u} has been established (e.g. Aas et al., 1986 ; Mesri, 1975 ; Lunne et al., 1997 ). However in overconsolidated clays the values of N_{k} often fall outside the predicted range and there are no constitutive modes to support empirical evidences.

A potential alternative to overcome the existing uncertainty related to N_{kt} is to use the excess pore pressure to estimate s_{u}. Several relationships have been proposed based on theoretical or semi-theoretical approaches using cavity expansion theory ( Battaglio et al., 1986 ; Campanella et al., 1985; Randolph & Wroth, 1979; Vesic, 1972):

${s}_{u}=\frac{{u}_{2}-{u}_{o}}{{N}_{\Delta u}}$ (2)

where N_{Δu} is shown from cavity expansion to vary in the range of 2 to 20. Lunne et al. (1997) recommend using a value of N_{Δu} of between 7 and 10. It has been advocated that these methods have the advantage of increased accuracy in the measurement of Δu, mainly in soft clays where Δu can be very large (e.g. Campanella et al., 1985 ).

The rationality in using Equation (2) is that cone resistance and excess pore pressures generated during cone penetration into fine grained soils will be dependent on the same parameters—stress state, stress history, soil stiffness—and can therefore be associated in CPTu predictions. Proposed semi-empirical solutions ( Massarch & Broms, 1981

A new method proposed by Mántaras et al. (2014) is applied here. This analysis advocates a different approach using dissipation tests and linking the measured piezocone maximum dissipation excess pore-water pressure (Δu_{max}) to the undrained shear strength (s_{u}). The mathematical solution proposed by Burns & Mayne (1998) is used as reference given the fact that excess pore water pressures are computed through a combination of the octahedral and the shear-induced components, allowing both normally and overconsolidated clays to be modeled from pore pressures measured immediately behind the cone shoulder (u_{2}).

The verification of the aforementioned method requires well-documented cases, which are not always available in engineering practice. This work aims to expand the database on which the method is based, and to verify the consistency of results within the empirical framework that relates cone factors and normalized excess pore pressure.

2. Mathematical Formulation ( Mántaras et al., 2014 )

The mathematical solution proposed by Burns & Mayne (1998) is based on the cavity expansion-critical state framework for the monotonic and dilatory response with regard to time. The excess pore water pressures, Δu, at any time can be compared with the initial values during penetration, $\Delta {u}_{i}={u}_{2i}-{u}_{0}$ , represented as:

$\Delta {u}_{i}={\left(\Delta {u}_{oc}\right)}_{i}+{\left(\Delta {u}_{\tau}\right)}_{i}$

where: ${\left(\Delta {u}_{oc}\right)}_{i}=\left(2/3\right)M{{\sigma}^{\prime}}_{vo}{\left(OCR/2\right)}^{\Lambda}\mathrm{ln}{I}_{r}$ is the octahedral component during penetration and ${\left(\Delta {u}_{\tau}\right)}_{i}={{\sigma}^{\prime}}_{vo}{\left[1-\left(OCR/2\right)\right]}^{\Lambda}$ is the shear-induced component during penetration, with OCR (the overconsolidation ratio), ${\varphi}^{\prime}$ the effective friction angle, $M=\left(6\mathrm{sin}{\varphi}^{\prime}\right)/\left(3-\mathrm{sin}{\varphi}^{\prime}\right)$ and L the compressibility ratio (1 − CC/CR). Departing from this concept, Burns & Mayne (1998) derived the following equation for the normalized excess pore pressure:

$\frac{\Delta u}{{{\sigma}^{\prime}}_{v}}=\frac{\frac{2}{3}M\cdot \mathrm{ln}\left({I}_{r}\right){\left[\frac{OCR}{2}\right]}^{\Lambda}}{1+50\cdot T}+\frac{1-{\left[\frac{OCR}{2}\right]}^{\Lambda}}{1+5000\cdot T}$ (3)

where: Ir is the rigidity index and T a dimensionless time factor.

Being Tmax the value of T where the excess pore pressures reaches its maximum value, it is possible to found it analitically from the equation:

$\partial \left(\Delta u/{{\sigma}^{\prime}}_{v}\right)/\partial T=0$ resulting the expression:

The maximum value of T could be obtained

${T}_{\mathrm{max}}=-\frac{1}{50}\frac{\left[\alpha {\left(\frac{OCR}{2}\right)}^{\Lambda}+1-{\left(\frac{OCR}{2}\right)}^{\Lambda}-9.9\sqrt{\alpha \cdot {\left(\frac{OCR}{2}\right)}^{\Lambda}\cdot \left[{\left(\frac{OCR}{2}\right)}^{\Lambda}-1\right]}\right]}{100\cdot \alpha \cdot {\left(\frac{OCR}{2}\right)}^{\Lambda}+1-{\left(\frac{OCR}{2}\right)}^{\Lambda}}$ (4)

Note that Tmax is the value of T where the excess pore pressure reaches its maximum value taking the first derivative of Equation (3) and setting the gradient of the function equal to zero, it is possible to determine the maximum value of the normalized excess pore pressure in relation to time factor T.

(Tmax) combined with Δumax and related to the normalized undrained shear strength from direct shear undrained test (Wroth, 1984) yields the undrained pore pressure ratio ${N}_{\Delta {u}_{\mathrm{max}}}\left(=\Delta {u}_{\mathrm{max}}/{s}_{u}\right)$ defined as the ratio of the maximum excess pore pressure and the undrained shear strength:

${N}_{\Delta {u}_{\mathrm{max}}}=\left(\frac{\Delta {u}_{\mathrm{max}}}{{s}_{u}}\right)=\frac{\frac{\frac{2}{3}M\cdot \mathrm{ln}\left({I}_{r}\right){\left[\frac{OCR}{2}\right]}^{\Lambda}}{1+50\cdot {T}_{\mathrm{max}}}+\frac{1-{\left[\frac{OCR}{2}\right]}^{\Lambda}}{1+5000\cdot {T}_{\mathrm{max}}}}{\frac{1}{2}\mathrm{sin}\left(\varphi \right)\cdot OC{R}^{\Lambda}\cdot \mathrm{log}\left({I}_{r}\right)}$ (5)

Based on Equation (5), the derived formulation shows little sensitivity to variations on OCR and ${\varphi}^{\prime}$ , and for typical soil parameters can be reduced to a relatively simple expression of a constant times the logarithm of the rigidity index, Ir (as illustrated on Figure 1 for OCR).

Figure 1. Undrained pore pressure ratio ${N}_{\Delta {u}_{\mathrm{max}}}$ (Mántaras et al.., 2014).

Advantages of interpreting pore pressures measurements compared to the more conventional piezocone penetration data are: (a) the cone factor ${N}_{\Delta {u}_{\mathrm{max}}}$ is less affected by soil rigidity and stress history because measurements result from pore pressure flow around the probe rather than of complete soil displacement caused by the cone penetration test and (b) well-defined failure mechanisms of flow around the probe allow for sound theoretical analysis of pore pressure dissipation and shear strength.

From the mathematical viewpoint the above expression is rigorously consistent; expressing the ratio between values of pore pressure measured during piezocone dissipation and the undrained shear strength in Direct Simple Shear conditions.

3. Case Studied

To illustrate the possibilities offered by the approach proposed by Mántaras et al. (2014), it will be applied to data from a well-documented site research program of a soft clay deposit in the metropolitan area of Porto Alegre, Rio Grande do Sul, Brazil. The complete research program consists of eleven study points using different types of tests in each one.

In order to visualize the results of the different tests and the application of the method, two points of very different characteristics were chosen to be presented in this work. In both cases, undrained shear strength values predicted from classical approaches are consistent between them (su obtained from Vane test, N_{kt} and N_{Δu}).

The first boring is illustrated in Figure 2, showing a homogeneous clay profile deposit. Undrained strength, su, obtained from the dissipation test method proposed by Mántaras et al. (2014) using ${N}_{\Delta {u}_{\mathrm{max}}}$ =8 ( ${N}_{\Delta {u}_{\mathrm{max}}}=4\ast \mathrm{log}{I}_{r}$ , with Ir equal to 100 = 8) are of the same order of magnitude as those calculated from other methods.

Figure 2. Metropolitan Region Porto Alegre, Brazil. (a) Homogeneous geotechnical profile. (b) Dissipation tests.

The second boring is illustrated in Figure 3, showing a layered deposit. The different grain size of layers could also be estimated from assessment of horizontal coefficient of consolidation, Ch, from dissipation tests. Once again, dissipation test results yielded su values of the same order of magnitude as those calculated from other methods.

Table 1 summarized the results for the whole campaign. Aiming to verify the consistency of assessed undrained strength with this novel approach, the cones factor N_{kt} and N_{Δu} values were back calculated in Table 1 and presented in Figure 4 and Figure 5 in the same framework proposed by Lunne et al. (1997).

Despite the dispersion of the results, it is clear the trends of N_{kt} vs. Bq and N_{Δu} vs. Bq values and are consistent with those found in international literature.

Figure 3. Metropolitan Region Porto Alegre, Brazil. (a) Erratic geotechnical profile. (b) Dissipation tests.

Figure 4. . Back calculated cone factors N_{kt} vs. B_{q} using s_{u} (
${\Delta {u}_{\mathrm{max}}}_{}$ method).

Figure 5. Back calculated cone factors N_{kt} vs. B_{q} using s_{u} (
${\Delta {u}_{\mathrm{max}}}_{}$ method).

Table 1. All the results obtained in the complete study.

4. Closing Remarks

The paper applies a new method proposed by Mántaras et al. (2014) to estimate the undrained shear strength from piezocone dissipation tests. Stress history, shear strength and compressibility are the critical factors affecting the accuracy of predictions and are properly taken into account.

The method yields an undrained pore pressure ratio ${N}_{\Delta {u}_{\mathrm{max}}}$ of about 8 which is within the range of early recommended values.

Also the N_{kt} and N_{Δu} factors calculated from the results of undrained resistance obtained from dissipation tests applying this new mathematical formulation are consistent with those referenced in the international literature and encourage its use in the engineering practice.

References

[1] Aas, G., Lacase, M., Lunne, T., & Hoeg, K. (1986). Use of in Situ Test for Foundation Design on Clay. In Proceedings ASCE Specialty Conference in Situ 86: Use of in Situ Tests in Geotechnical Engineering (pp. 1-30). New York.

[2] Baligh, M. M. (1985). Strain Path Method. Journal of the Soil Mechanics and Foundations Division, 11, 1108-1136.

https://doi.org/10.1061/(ASCE)0733-9410(1985)111:9(1108)

[3] Battaglio, M., Bruzzi, D., Jamiolkowski, M., & Lancellotta, R. (1986). Interpretation of CPTs and CPTUs. In Proceedings 4th International Geotechnical Seminar (pp. 129-143). Singapore.

[4] Burns, S. E., & Mayne, P. W. (1998). Monotonic and Dilatory Pore-Pressure Decay during Piezocone Tests in Clay. Canadian Geotechnical Journal, 35, 1063-1073.

https://doi.org/10.1139/cgj-35-6-1063

[5] Campanella, R. G., Robertson, P. K., Gillespie, D., & Greig, J. (1985). Recent Developments in Situ Testing of Soils. In Proceedings of XI ICSMFE (Vol. 2, pp. 849-854). San Francisco, California.

[6] Caquot, A., & Kérisel, J. (1953). Sur le Terme de Surfacedans Le Calcul de Fondations en Milieu Puérulent. In Proceedings of the Third International Conference on Soil Mechanics and Foundation Engineering (Vol. 1). Switzerland

[7] Lunne, T., Robertson, P. K., & Powell, J. J. M. (1997). Cone Penetration Testing in Geotechnical Practice (312 p). Boca Raton, FL: CRC Press.

[8] Mántaras, F. M., Odebrecht, E., & Schnaid, F. (2014). Using Piezocone Dissipation Test to Estimate the Undrained Shear Strength in Cohesive Soil. Canadian Geotechnical Journal, 52, 318-325.

https://doi.org/10.1139/cgj-2014-0176

[9] Massarch, K. R., & Broms, B. B. (1981). Pile Driving in Clay Sloples. In Proceedings International Conference Soil Mechanics and Foundation (Vol. 3, pp. 469-474). Stockholm.

[10] Mesri, G. (1975). Discussion on “New Design Procedure for Stability of Soft Clays”. Journal of the Geotechnical Engineering Division, ASCE, 101, 409-412.

[11] Meyerhof, G. G. (1956). Penetration Tests and Bearing Capacity of Cohesionless Soils. Journal of the Soil Mechanics and Foundations Division, ASCE, 82, 1-19.

[12] Randolph, M. F., & Wroth, C. P. (1979). An Analytical Solution for the Consolidation around a Driven Pile. International Journal for Numerical and Analytical Methods in Geomechanics, 3, 217-229.

https://doi.org/10.1002/nag.1610030302

[13] Schnaid, F. (2009). In Situ Testing in Geomechanics (Vol. 1. 329 p.) Oxon: Taylor & Francis.

[14] Teh, C. I., & Houlsby, G. T. (1991). An Analytical Study of the Cone Penetration Test in Clay. Géotechnique, 41, 17-34.

https://doi.org/10.1680/geot.1991.41.1.17

[15] Terzaghi, R. (1943). Theoretical Soil Mechanics (510 p.). New York: John Wiley & Sons.

https://doi.org/10.1002/9780470172766

[16] Vesic, A. S. (1972). Expansion of Cavities in Infinite Soil Mass. Journal of the Soil Mechanics and Foundations Division, ASCE, 98, 265-290.

[17] Wroth, C. P. (1984). The Interpretation of in Situ Soil Test. 24th Rankine Lecture. Géotechnique, 34, 449-489.

https://doi.org/10.1680/geot.1984.34.4.449

[18] Yu, H. S. (2004). In Situ Soil Testing: From Mechanics to Interpretation, James K. Mitchell Lecture. In Proceedings on Geotechnical and Geophysical Characterization (Vol. 1, pp. 3-38). Porto: Mill Press.

[19] Yu, H. S., & Whittle, A. J. (1999). Combining Strain Path Analysis and Cavity Expansion Theory to Estimate Cone Resistance in Clay. Unpublished Notes.