Numerical modeling of skin friction and penetration problems in geotechnical engineering

Abstract

Numerical modeling using finite element method (FEM) is well-recognized as a powerful method for both engineers and researchers to solve boundary value problems. In the modeling of geotechnical problems, the analyses are often limited to simple static problems with either steady-state effective or total stress approach while the transient response (development and dissipation of excess pore water pressure, uex) is seldom considered. Besides, infinitesimal small soil deformation is usually assumed. The simulation is further complicated when the soil-structure interaction problems involve significant soil displacements; like a pile subject to negative skin friction (NSF) and a cone/pile penetration. However, conventional FEM analysis prematurely terminates due primarily to excessive mesh distortion. One could see that simulating a transient problem with large deformation and distortion remains a great challenge. In this study, advanced FE simulations are performed to give new insights into the problems of (1) a pile subject to NSF; and (2) a cone penetration. The transient response of the NSF problem is modeled with the fluid-coupled consolidation technique and geometric nonlinearity. The fluid-coupled cone penetration problem is modeled with a newly developed adaptive approach.

The NSF and cone penetration simulations involve complex soil-structure interface modeling. Two types of modified interface responses are developed and verified which consider fluid coupling. The developed algorithm is applied to back analyze a case history of a pile subject to NSF induced by surcharge loading. Promising results were shown. Development of dragload and neutral plane (NP) with time is studied. NP locates at 75% of the pile embedded length (D) in long-term. Next, a parametric study is performed to investigate the influences of pile geometries, ground compressibility and loading conditions towards the pile responses. The long-term NP locates at around 0.55D to 0.65D in the studied engineering scenarios. The maximum downdrag can be up to 10% of the pile diameter. NP shifts upward when the head load increases. A simple design chart is proposed which helps engineers to estimate the long-term axial load distribution. An illustrative example is given to demonstrate the application and performance of the chart.

The study is extended to investigate the cone penetration problem. An advanced adaptive method is developed and implemented into the FE package ABAQUS to resolve the problems of numerical instability, excessive mesh distortion and premature termination. The proposed method is verified by modeling a ground consolidation problem. Next, total stress back analysis of cone penetration is conducted with the proposed method. The development of cone factor predicted by the proposed method gives a better match with the laboratory result when comparing with the built-in ALE method. Next, the development and dissipation of uex during cone advancing with the proposed method and fluid-coupled technique is investigated. uex develops dramatically around the cone tip. The soil permeability is back calculated from the dissipation test and agrees well with the input value. It is believed that the construction effects of a press-in pile and the subsequence NSF on that pile can be modeled by utilizing the finding of this study.