Computational & Technology Resources
an online resource for computational,
engineering & technology publications
Computational Science, Engineering & Technology Series
TRENDS IN ENGINEERING COMPUTATIONAL TECHNOLOGY
Edited by: M. Papadrakakis, B.H.V. Topping
Three-Dimensional Numerical Analysis of a Dynamic Structure, Saturated Soil and Pore Fluid Interaction Problem
A.H.C. Chan1 and J. Ou2
1Department of Civil Engineering, University of Birmingham, United Kingdom
A.H.C. Chan, J. Ou, "Three-Dimensional Numerical Analysis of a Dynamic Structure, Saturated Soil and Pore Fluid Interaction Problem", in M. Papadrakakis, B.H.V. Topping, (Editors), "Trends in Engineering Computational Technology", Saxe-Coburg Publications, Stirlingshire, UK, Chapter 17, pp 335-353, 2008. doi:10.4203/csets.20.17
Keywords: three-dimensional, soil-structure interaction, earthquake, Biot formulation, soil dynamics.
A numerical method for the calculation of the overall behaviour of a non-trivial three-dimensional saturated soil-structure system is presented in this paper. The numerical modelling of saturated soil is discussed in detail together with comparisons with soil-structure physical experiments performed on the centrifuge facilities at the Rensselaer Polytechnic Institute for the VELACS project.
In order to model the saturated soil, three main components are needed for the numerical simulation and they are as follows:
The numerical procedure employs the Biot formulation as the basic mathematical framework, the Finite Element method for spatial discretisation and the Generalised Newmark scheme for time stepping . The u-p formulation is used with displacement (u) and pore pressure (p) as the primary unknowns and the constitutive relationship used is the Pastor-Zienkiewicz mark III model .
The numerical procedure is validated using Test No.12 from the VELACS  experiments, is a dynamic three-dimensional soil-structure interaction test with a model of a structure in the form of a block being placed on a saturated sand layer overlaid by a thin silt layer and subjected to earthquake-like base motion. The tests were performed at 100g centrifugal acceleration at Princeton University, Rensselaer Polytechnic Institute and University of California, Davis .
Very good comparisons with the results from the physical centrifuges test had been obtained. All the accelerations predicted were in good agreement with the measured values. Very good comparisons were also obtained for most of the pore water pressure transducers though in some locations, the level of excess pore pressure rise was slightly lower than that measured. Around 10cm settlement was predicted by the numerical analysis which is between the value obtained at UC Davis (9cm) and RPI (13cm). But a very wide spread of the measured vertical settlement was observed therefore the measured final settlement value may not be too trustworthy for comparison purposes however the numerical prediction correctly modelled the trend at which the settlement developed.
These results confirmed that the assumptions and the approximation used in the numerical procedure are sound and further analyses on various three-dimensional examples and further investigations on the use of iterative solvers instead of the direct solver used in this study are underway.
purchase the full-text of this chapter (price £20)