A practical numerical approach for large deformation problems in soil

A practical method is presented for numerical analysis of problems in solid (in particular soil) mechanics which involve large strains or deformations. The method is similar to what is referred to as arbitrary Lagrangian-Eulerian, with simple infinitesimal strain incremental analysis combined with regular updating of co-ordinates, remeshing of the domain and interpolation of material and stress parameters. The technique thus differs from the Lagrangian or Eulerian methods more commonly used. Remeshing is accomplished using a fully automatic remeshing technique based on normal offsetting, Delaunay triangulation and Laplacian smoothing. This technique is efficient and robust. It ensures good quality shape and distribution of elements for boundary regions of irregular shape, and is very quick computationally. With remeshing and interpolation, small fluctuations appeared initially in the load-deformation results. In order to minimize these, different increment sizes and remeshing frequencies were explored. Also, various planar linear interpolation techniques were compared, and the unique element method found to work best.Application of the technique is focused on the widespread problem of penetration of surface foundations into soft soil, including deep penetration of foundations where soil flows back over the upper surface of the foundation. Numerical results are presented for a plane strain footing and an axisymmetric jack-up (spudcan) foundation, penetrating deeply into soil which has been modelled as a simple Tresca or Von Mises material, but allowing for increase of the soil strength with depth. The computed results are compared with plasticity solutions for bearing capacity. The numerical method is shown to work extremely well, with potential application to a wide range of soil-structure interaction problems.