Finite-Element Analysis of Inclined Piezocone Penetration Test in Clays

A three-dimensional finite-element analysis was performed to analyze the effect of soil anisotropy on the inclined piezocone penetration test in normally consolidated clay. The piezocone penetration was numerically simulated based on a large strain formulation using the commercial finite-element code ABAQUS, and the anisotropic modified cam clay model (AMCCM) was chosen and implemented into ABAQUS through the user subroutine UMAT. For verification purposes, numerical simulations were first performed on previously conducted calibration chamber tests, and the predicted results were compared with the measured values. For different initial stress conditions and different penetration angles, the cone tip resistance profile; excess pore pressure profile at the cone tip; typical stress, strain and excess pore pressure distributions around the cone; and excess pore pressure dissipation at the cone tip are provided. This study shows that when the initial stress state is anisotropic, the soil behavior is different under different angles of penetration.     

Determination of Pile Base Resistance in Sands

Advances in the design of axially loaded piles are desirable because significant cost savings may result. Well-designed piles settle by amounts that are well tolerated by the superstructure and induce strains around the pile base that are far removed from failure. To investigate the development of base resistance for a given soil condition and increasing settlements, piles embedded in sand are modeled using the finite-element method with a nonlinear elastic-plastic model. Based on the load-settlement response obtained from the finite-element analysis and cone penetration resistance obtained from cavity expansion and stress rotation analyses, values of normalized base resistance, defined as base resistance divided by cone penetration resistance, are obtained. The relationship between base resistance and cone resistance is useful in the design of deep foundation using cone penetration test results. The effect of the initial coefficient of earth pressure at rest K0 on normalized base resistance values is also investigated. Several case histories, including both nondisplacement and displacement piles, are used for comparison with the theoretical results.     

Immediate Settlement of Shallow Foundations Bearing on Clay

Immediate and long-term settlement checks are an integral part of foundation design. Therefore, reasonably accurate estimates of the immediate settlement of shallow foundations bearing on clay are necessary, particularly for highly plastic clays or organic soils, for which the immediate settlement may be significant. This immediate settlement is due entirely to the distortion of the clay underneath the shallow foundations because, in the short term, there is no opportunity for change in the clay volume. Since soil stress-strain response is nonlinear even at small strains, design procedures based on linear elasticity cannot accurately predict soil deformations. Hence, an immediate settlement analysis that takes soil nonlinearity into account is needed. In this paper, finite-element analysis is used to develop design charts that can be used to estimate the immediate settlement of axially loaded square, rectangular, and strip footings bearing on clay. The clay is modeled with a simple nonlinear constitutive relationship. A design example is included to illustrate how the proposed procedure can be readily applied in practice with the knowledge of the undrained shear strength and the initial shear modulus of the clay.     

Evaluation of Remolded Shear Strength and Sensitivity of Soft Clay Using Full-Flow Penetrometers

The undrained remolded shear strength of soft clays is of importance in geosystem design, particularly for offshore structures. Common methods to estimate remolded shear strength, such as correlations with cone penetration data, direct measurement with an in situ field vane shear device, and laboratory measurements, produce varied results and can be particularly costly and time consuming. Full-flow penetrometers (T-bar and Ball) provide an alternative rapid method to estimate remolded shear strength and soil sensitivity through remolding soil by repeated cycling of the penetrometer up and down over a given depth interval. The cyclic penetration resistance degradation curve inherently contains information regarding remolded strength and sensitivity. The objective of this paper is to assess the ability of full-flow penetrometers to predict remolded strength and soil sensitivity, and to develop a suite of predictive correlations in which these properties can be estimated in the absence of complementary laboratory or in situ test data. To accomplish this, full-flow penetration profiles and cyclic tests were performed at five well characterized soft clay sites, which together represent the broad range of soils in which the penetrometers will be often used. A previously developed model for the reduction in penetration resistance with cycling is modified to predict the entire degradation curve, including the remolded penetration resistance using only measurements obtained during initial penetrometer penetration and extraction. Using field vane shear strength as the reference measurement, correlations are developed to predict soil sensitivity and remolded shear strength based solely on full-flow penetrometer data, which is particularly useful in site investigation programs where site specific data are not yet available or are sparse. Finally, the usefulness of these relationships is demonstrated by implementing them for two additional soft clay sites.     

Determination of Hydraulic Conductivity Using Piezocone Penetration Test

An elastoplastic, finite-strain, coupled theory of mixtures in an updated Lagrangian reference frame is applied to the piezocone penetration test to estimate the hydraulic conductivity of the soil via analysis of the steady-state excess pore pressure generated during piezocone penetration. The results of this approach were compared with piezocone penetration test data. It showed that reliable hydraulic conductivities can be estimated conveniently without performing pore pressure dissipation tests. This study also shows that the change in the dimensionless excess pore pressure (excess pore pressure is normalized by the effective overburden pressure) at the cone tip is almost constant when the dimensionless hydraulic conductivity (hydraulic conductivity is normalized by the penetration speed and cone radius, hereafter called DLHC) is less than 10?7 or greater than 10?4. It is also shown that the drainage condition around the cone tip is close to a fully undrained condition when the DLHC of the soil is less than 10?7, while it is close to a fully drained condition when the DLHC of the soil is greater than 10?4.     

Multiple-Porosity Contaminant Transport by Finite-Element Method

An exponential finite-element model for multiple-porosity contaminant transport in soils is proposed. The model combines three compartments for dissolved contaminants: a primary compartment of diffusion–advection transport with nonequilibrium sorption, a secondary compartment with diffusion in rectangular or spherical soil blocks, and a tertiary compartment for immobile solutions within the primary compartment. Hence the finite-element model can be used to solve four types of mass-transfer problems which include: (1) intact soils, (2) intact soils with multiple sources of nonequilibrium partitioning, (3) soils with a network of regularly spaced fissures, and (4) structured soils. Hitherto, mobile/immobile compartments, fissured soils, and nonequilibrium sorption have been treated separately or in pairs. A Laplace transform is applied to the governing equations to remove the time derivative. A Galerkin residual statement is written and a finite-element method is developed. Both polynomial and exponential finite elements are implemented. The solution is inverted to the time domain numerically. The method is validated by comparison to analytical and boundary element predictions. Exponential elements perform particularly well, speeding up convergence significantly. The scope of the method is illustrated by analyzing contamination from a set of four waste repositories buried in fissured clay.     

黏土中静压沉桩离心模型

采用西澳大学室内鼓轮式离心机,在预先固结的高岭黏土中开展不同离心力场(50g,125g及250g,g为重力加速度)条件下的模型压桩试验、T-bar试验和静力触探试验,分析了模型桩在贯入过程、静置稳定过程中桩身径向应力(σ_r)的变化规律,并对后期桩体拉伸载荷阶段的径向应力变化值(Δσ_r)及桩侧摩阻力变化情况行了探讨,揭示了在不同超固结比(OCRs)黏土中静压桩侧摩阻力的演变特性.在此基础上,通过两种经验公式方法对桩侧摩承载力进行了预测计算和对比分析.研究结果表明:沉桩过程中桩端相对高度(h/B)对桩身径向应力的发展变化有很大的影响,桩身不同位置(h/B)的总径向应力对同一贯入深度而言,存在桩侧径向应力退化现象;基于静力触探试验提出的经验方法,能有效考虑静力触探锥端阻力(q_t)和桩端相对高度(h/B)因素的影响,将其应用于黏土沉桩时桩侧摩阻力的预测,可取得与试验实测结果较吻合的结果.研究成果对软土地区静压桩施工与承载力设计具有一定的工程指导意义.      

Evolution of the Pore-Pressure Field around a Moving Conical Penetrometer of Finite Size

A solution is developed for the evolution of buildup, steady, and postarrest dissipative pore-fluid pressure fields that develop around a finite-radius conical penetrometer advanced in a saturated linearly elastic porous medium. The analog with cone penetrometer testing is direct and is used to enable continuous distributions of permeability and diffusivity to be determined with depth. This analysis reveals the direct dependence of penetration rate on the induced fluid pressure field magnitudes, and predicts that a penetration rate threshold limit exists with respect to pore-pressure generation. This represents the essence of a partially drained system. The developed pore-pressure field is determined to be a function of the dissipation rate of the material, the penetration rate, and the storage effects of the advecting medium. Analysis of the pore-pressure field under start-up conditions reveals that the time required to reach steady state is strongly influenced by the penetration rate and the pressure-dissipation properties of the material. Analysis of the developed stable pressure fields illustrates the inversely proportional relationship that exists between penetration rate and pore-pressure magnitudes at the cone surface; representing the influence of storage in the medium on stable pore-pressure magnitudes. Stable pressure fields below the penetration threshold limit, UD ? 10?1, form a spherical response around the cone tip transitioning to an elongated radial response for penetration rates above this limit. Postarrest analysis indicates that the prearrest penetration rate strongly influences the dissipation rate and pattern of dissipation. The developed analysis can be correlated with CPTu-recovered data to independently evaluate permeability magnitudes during steady penetration.     

Three-Dimensional Large Deformation FE Analysis of Square Footings in Two-Layered Clays

In strong over soft two-layered clays, there is a potential for the footing to experience a punch-through failure, where the footing penetrates a large distance at a short time after the initial peak resistance is reached. Three-dimensional (3D) large deformation finite-element analyses using 3D RITSS method were conducted to simulate the penetration responses of square footings in strong over soft clays. The effects of surface soil heave and soil layer interface deformation during footing penetration were studied in weightless soils. Fitted equations were proposed to express the footing capacity response against the penetration depth. Based on the fitted equations, formulas to calculate footing peak bearing factor and the corresponding penetration depth were developed. The peak footing capacity factor and the corresponding penetration depth increases with the increasing of soil layer strength ratio, relative top soil layer thickness and soil weight factor, thus the potential of punch-through failure was reduced accordingly. It was also found that the soil weight effect can be a simple surcharge based on the formula developed in the weightless soil. Design charts for the peak footing capacity factor and the corresponding penetration depth were developed.     

Determination of Strength and Index Properties of Fine-Grained Soils Using a Soil Minipenetrometer

This paper describes the use of a soil minipenetrometer (SMP) to determine the strength and index properties of fine-grained soils. The SMP has been developed to allow both fall cone and quasi-static penetration tests to be performed. Displacement controlled quasi-static penetration tests can be used for the direct measurement of undrained shear strength, both for remolded and undisturbed samples. In addition the quasi-static penetration test can be used to define an additional lower plastic limit parameter, the PL100, which represents the moisture content of a fine-grained soil with an undrained strength 100 times that at the liquid limit. This approach offers the advantage that removal of the coarse fraction is not required to estimate the PL100.     

Liquefaction Resistance of Soils from Shear-Wave Velocity

A simplified procedure using shear-wave velocity measurements for evaluating the liquefaction resistance of soils is presented. The procedure was developed in cooperation with industry, researchers, and practitioners and evolved from workshops in 1996 and 1998. It follows the general format of the Seed-Idriss simplified procedure based on standard penetration test blow count and was developed using case history data from 26 earthquakes and >70 measurement sites in soils ranging from fine sand to sandy gravel with cobbles to profiles including silty clay layers. Liquefaction resistance curves were established by applying a modified relationship between the shear-wave velocity and cyclic stress ratio for the constant average cyclic shear strain suggested by R. Dobry. These curves correctly predicted moderate to high liquefaction potential for >95% of the liquefaction case histories and are shown to be consistent with the standard penetration test based curves in sandy soils. A case study is provided to illustrate application of the procedure. Additional data are needed, particularly from denser soil deposits shaken by stronger ground motions, to further validate the simplified procedure.     

Finite Strain, Anisotropic Modified Cam Clay Model with Plastic Spin.?II: Application to Piezocone Test

Micromechanical and anisotropic behavior of soils is analyzed using the plastic spin and anisotropic modified Cam clay model (AMCCM) in an updated Lagrangian reference frame. The micromechanical behavior of soils is especially important in large strain problems. The anisotropic behavior inherently exists in soils. This study uses AMCCM with the plastic spin to evaluate the pore pressure response in cone penetration tests. This test is a typical example for large strain problems. The analytical results indicate that the predicted pore pressure is closer to the measured one when the plastic spin is incorporated in AMCCM. It also shows that the micromechanical behavior is paramount in the vicinity of the cone tip where large strains and rotations occur.     

Liquefaction and Soil Failure During 1994 Northridge Earthquake

The 1994 Northridge, Calif., earthquake caused widespread permanent ground deformation on the gently sloping alluvial fan surface of the San Fernando Valley. The ground cracks and distributed deformation damaged both pipelines and surface structures. To evaluate the mechanism of soil failure, detailed subsurface investigations were conducted at four sites. Three sites are underlain by saturated sandy silts with low standard penetration test and cone penetration test values. These soils are similar to those that liquefied during the 1971 San Fernando earthquake, and are shown by widely used empirical relationships to be susceptible to liquefaction. The remaining site is underlain by saturated clay whose undrained shear strength is approximately half the value of the earthquake-induced shear stress at this location. This study demonstrates that the heterogeneous nature of alluvial fan sediments in combination with variations in the ground-water table can be responsible for complex patterns of permanent ground deformation. It may also help to explain some of the spatial variability of strong ground motion observed during the 1994 earthquake.     

Deformation and Pore-Water Pressure Responses of Elastic Viscoplastic Soil

The deformation and pore-water pressure responses of clayey soils are of great interest to civil engineers. In this paper, displacements and pore-water pressures of a clay subjected to the loading of a strip footing are simulated using a fully coupled finite-element (FE) consolidation method incorporated with a newly developed 3D elastic viscoplastic (EVP) model for the clay. A brief introduction to the 3D EVP model and its implementation in the FE analysis is presented. The 3D EVP model can describe the time-dependent stress-strain behavior of clayey soils, including volumetric creep. The main objective of this paper is to examine how the viscosity (or creep parameter ψ/V) of the clay affects the deformation and pore-water pressure responses of the clay. For this, the value of the creep parameter ψ/V is varied in the FE analysis. When viscous nature is taken into account, the pore-water pressure in the soil is higher than that without consideration of the viscous nature. The phenomenon of pore-water pressure increase due to creep is studied in this paper. It is found that larger creep parameter ψ/V results in higher pore-water pressure and larger deformation in the soil. The difference of the pore-water pressure due to the Mandel-Cryer effects and the creep is investigated using the FE model and discussed in this paper. In addition, a few other parameters (Poisson's ratio ν, permeability k, clay layer thickness h, and thickness h to a half footing width a ratio h∕a) are also varied to investigate their influence on deformation and pore-water pressure of the soil with creep. It is found that, the lower the permeability of soil, the higher is the pore-water pressure and the larger is the local deformation. The thickness of the soil layer also has a great influence on the pore-water pressure induced by the viscous effect. All these increased pore-water pressures result from a balance of the pore-water pressures induced by creep (and the Mandel-Cryer effects or both) and dissipated because of drainage.     

Effect of Penetration Rate on Cone Penetration Resistance in Saturated Clayey Soils

In this paper, the effects of penetration rate on cone resistance in saturated clayey soils are investigated. Shear strength rate effects in clayey soils are related to two physical processes: the increase of shear strength with increasing rate of loading and the increase of shear strength as the process transitions from undrained to drained. Special focus is placed on this second effect. Cone penetration tests were performed at various penetration rates both in the field and in a calibration chamber, and the resulting data were analyzed. The field cone penetration tests were performed at two test sites with fairly homogeneous clayey silt and silty clay layers located below the groundwater table. Additionally, tests with both cone and flat-tip penetrometers in sand-clay mixtures were performed in a calibration chamber to investigate the change in drainage conditions from undrained to partially drained and from partially drained to fully drained. A series of flexible-wall permeameter tests were conducted in the laboratory for various clayey sand mixtures prepared at various mixing ratios in order to obtain values of the coefficient of consolidation, which is required to estimate the penetration rates below which penetration is drained and above which penetration is undrained. A correlation between cone resistance and drainage conditions was established based on the results of the calibration chamber and field penetration tests.     

Relationship between the Undrained Shear Strength, Water Content, and Mineralogical Properties of Fine-Grained Soils

The relationship between the undrained shear strength of fine-grained soils and the water content can be described with a nonlinear function in which the type of soil is determined by two parameters. It is well known that these parameters depend mainly on the mineral compositions of soils; these relationships, however, have not yet been investigated. The findings described in this paper define those mineralogical properties of soils which determine the values of both parameters. Experimentally obtained results suggest that the parameters primarily depend on the size of the clay minerals, their quantity in soil composition, and the interlayer water quantity in the expanding clay minerals. As this dependence is well defined, the parameters, and thus the undrained shear strength at different water content, can be defined from knowledge of these mineralogical soil properties.     

Numerical Parametric Study of Piezocone Penetration Test in Clays

This paper presents a numerical model for the simulation of the piezocone penetration test that is used to carry out a parametric study of the piezocone penetration test in cohesive soils. The piezocone penetration is numerically simulated using an axisymmetric elasto-plastic large deformation finite-element analysis code. The numerical simulation is accomplished in two stages. First, the piezocone is expanded radially from an initial small radius (0.1ro) to the piezocone radius, ro, at the specified depth. Second, the continuous penetration of the piezocone is simulated by applying incremental vertical displacements of the nodes representing the piezocone boundary. The constraint approach is used to model the soil-piezocone interface friction. The Mohr-Coulomb frictional criterion is used to define the sliding potential of the nodes. The main objective of this paper is to present the numerical model and to investigate the effect of the lateral and vertical stresses and the overconsolidation ratio on the cone tip resistance and the developed excess pore pressure around the piezocone. The variation of the horizontal and vertical hydraulic conductivity coefficients on the developed spatial excess pore pressure and its dissipation are also investigated. The results of the numerical study are also compared with the miniature piezocone penetration tests in cohesive soil specimens conducted at the Louisiana State University Calibration Chamber. Results of this study are in good agreement with the measured values.     

Analysis of Factors Influencing Soil Classification Using Normalized Piezocone Tip Resistance and Pore Pressure Parameters

This paper discusses the development of a framework for classifying soil using normalized piezocone test (CPTU) data from the corrected tip resistance (qt) and penetration pore-water pressure at the shoulder (u2). Parametric studies for normalized cone tip resistance (Q = qcnet/σv0′) and normalized excess pressures (Δu2/σv0′) as a function of overconsolidation ratio (OCR = σvy′/σv0′) during undrained penetration are combined with piezocone data from clay sites, as well as results from relatively uniform thick deposits of sands, silts, and varietal clays from around the globe. The study focuses on separating the influence of yield stress ratio from that of partial consolidation on normalized CPTU parameters, which both tend to increase Q and decrease the pore pressure parameter (Bq = Δu2/qcnet). The resulting recommended classification chart is significantly different from existing charts, and implies that assessment of data in Q–Δu2/σv0′ space is superior to Q–Bq space when evaluating piezocone data for a range of soil types. Still, there are zones of overlap for silty soils and heavily overconsolidated clays, thus requiring that supplementary information to Q and Δu2/σv0′ be obtained in unfamiliar geologies, including variable rate penetration tests, dissipation tests, CPT friction ratio, or soil sampling.     

Estimation of Footing Settlement in Sand

The settlement of foundations under working load conditions is an important design consideration. Well‐designed foundations induce stress‐strain states in the soil that are neither in the linear elastic range nor in the range usually associated with perfect plasticity. Thus, in order to accurately predict working settlements, analyses that are more realistic than simple elastic analyses are required. The settlements of footings in sand are often estimated based on the results of in situ tests, particularly the standard penetration test (SPT) and the cone penetration test (CPT). In this article, we analyze the load‐settlement response of vertically loaded footings placed in sands using both the finite element method with a nonlinear stress‐strain model and the conventional elastic approach. Calculations are made for both normally consolidated and heavily overconsolidated sands with various relative densities. For each case, the cone penetration resistance qc is calculated using CONPOINT, a widely tested program that allows computation of qc based on cavity expansion analysis. Based on these analyses, we propose a procedure for the estimation of footing settlement in sands based on CPT results.     

Effect of Surface Heave on Response of Partially Embedded Pipelines on Clay

The as-laid embedment of an on-bottom pipeline strongly influences the resulting thermal insulation, and the resistance to subsequent axial and lateral movement of the pipeline. Reliable assessment of these parameters is essential for the design of offshore pipelines. Static vertical penetration of a pipe into a soft clay seabed—which can be modeled as an undrained process—causes heave of soil on each side of the pipeline. The heaved soil contributes to the vertical penetration resistance and the horizontal capacity. This paper describes a series of large deformation finite-element analyses of pipe penetration, supported by a simple analytical assessment of the heave process. The conventional bearing capacity approach to the analysis of pipe penetration is reviewed, and modifications for the effects of soil weight and heave are presented. It is shown that in soft soil conditions—which are typical for deep water—the soil self-weight contributes a significant portion of the vertical penetration resistance and horizontal capacity. If heave is neglected, the soil weight leads to a vertical force due to buoyancy, based on Archimedes’ principle. When heave is considered, the soil weight contributes an additional component of vertical load, exceeding simple buoyancy, due to the distorted geometry of the soil surface. Archimedes’ principle does not apply. The finite-element analyses, benchmarked against rigorous plasticity solutions, are used to calibrate simple expressions for predicting static vertical pipe penetration, and the resulting horizontal capacity. These simple solutions allow the conventional bearing capacity approach to be used in a manner which correctly accounts for the effects of soil self-weight and heave. An approximate solution for predicting the “local” pipe embedment—relative to the raised soil level immediately adjacent to the pipe—is derived. The local embedment significantly exceeds the nominal embedment relative to the original soil surface. This effect counteracts the tendency for heave to reduce the embedment by raising the penetration resistance.