• Geomechanics and Engineering
• /
• v.8 no.4
• /
• pp.559-594
• /
• 2015

A nonlinear creep function incorporated into the elastic visco-plastic model may describe the long-term soil deformation more accurately. However, by applying the conventional procedure, there are challenges to determine the model parameters due to limitation of suitable data points. This paper presents a numerical solution to obtain several parameters simultaneously for a nonlinear elastic visco-plastic (EVP) model using the available consolidation data. The finite difference scheme using the Crank-Nicolson procedure is applied to solve a set of coupled partial differential equations of the time dependent strain and pore water pressure dissipation. The model parameters are determined by applying the algorithm of trust-region reflective optimisation in conjunction with the finite difference solution. The proposed method utilises all available consolidation data during dissipation of the excess pore water pressure to determine the required model parameters. Moreover, the reference time in the elastic visco-plastic model can readily be adopted as a unit of time; denoting creep is included in the numerical predictions explicitly from the very first time steps. In this paper, the settlement predictions of thick soft clay layers are presented and discussed to evaluate and compare the accuracy and reliability of the proposed method against the graphical procedure to obtain the model parameters. In addition, comparison of the available experimental results to the numerical predictions confirms the accuracy of the numerical procedure.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.