One-step and linear multistep methods for nonlinear consolidation

Abstract The behavior of numerical ODE methods for the solution of nonlinear consolidation is investigated. One-step algorithms based on the generalized trapezoidal method and linear multistep methods based on the backward differentiation formulas (BDF) are considered. Mechanical deformation is coupled with fluid flow via a formulation based on Biot's three-dimensional consolidation theory. The resulting nonlinear FE matrix equations are cast in such a form that the one-step and the BDF methods can be coded in a similar fashion via the use of a displacement difference operator. Accuracy and stability analyses are presented for both the one-step and BDF methods. In particular, the effect of variable step size on the stability of the BDF methods is discussed in the context of a zero-stability criterion so that bounds on a uniformly growing step size can be established. The generalized trapezoidal methods and the backward differentiation formulas of order up to three are shown to converge with the expected order, save for round-off errors, for nonlinear consolidation problems where the soil is modeled as an elasto-plastic material of the Cam-Clay type.

[1]  Per Lötstedt,et al.  Numerical solution of nonlinear differential equations with algebraic constraints I: Convergence results for backward differentiation formulas , 1986 .

[2]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[3]  Ronaldo I. Borja,et al.  Cam-Clay plasticity, part I: implicit integration of elastoplastic constitutive relations , 1990 .

[4]  W. Spaans.  The finite element methods , 1975 .

[5]  Jean H. Prevost,et al.  Implicit-explicit schemes for nonlinear consolidation , 1983 .

[6]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[7]  J. Dieudonne.  Foundations of Modern Analysis , 1969 .

[8]  G. Pinder,et al.  Computational Methods in Subsurface Flow , 1983 .

[9]  C. W. Gear,et al.  Simultaneous Numerical Solution of Differential-Algebraic Equations , 1971 .

[10]  C. W. Geart,et al.  THE EFFECT OF VARIABLE MESH SIZE ON THE STABILITY OF MULTISTEP METHODS , 1974 .

[11]  James M. Ortega,et al.  Iterative solution of nonlinear equations in several variables , 2014, Computer science and applied mathematics.

[12]  J. C. Small,et al.  Elasto-plastic consolidation of soil , 1976 .

[13]  J. C. Small,et al.  An investigation of the stability of numerical solutions of Biot's equations of consolidation , 1975 .

[14]  J. C. Small,et al.  The analysis of finite elasto-plastic consolidation , 1979 .

[15]  Ronaldo I. Borja,et al.  Linearization of elasto‐plastic consolidation equations , 1989 .

[16]  C. W. Gear,et al.  Automatic integration of Euler-Lagrange equations with constraints , 1985 .

[17]  Chandrakant S. Desai,et al.  Two numerical schemes for nonlinear consolidation , 1981 .

[18]  Ronaldo I. Borja,et al.  Cam-Clay plasticity, Part II: implicit integration of constitutive equation based a nonlinear elastic stress predictor , 1991 .

[19]  J. C. Simo,et al.  Non‐smooth multisurface plasticity and viscoplasticity. Loading/unloading conditions and numerical algorithms , 1988 .

[20]  K. Park.  An Improved Stiffly Stable Method for Direct Integration of Nonlinear Structural Dynamic Equations , 1975 .

[21]  K. Roscoe,et al.  ON THE GENERALIZED STRESS-STRAIN BEHAVIOUR OF WET CLAY , 1968 .

[22]  Thomas J. R. Hughes,et al.  Implicit-explicit finite elements in nonlinear transient analysis , 1979 .

[23]  Edward Kavazanjian,et al.  A constitutive model for the stress–strain–time behaviour of ‘wet’ clays , 1985 .

[24]  M. Biot.  THEORY OF ELASTICITY AND CONSOLIDATION FOR A POROUS ANISOTROPIC SOLID , 1955 .

[25]  A. Nordsieck,et al.  On Numerical Integration of Ordinary Differential Equations , 1953 .

[26]  Roland W. Lewis,et al.  A comparison of time marching schemes for the transient heat conduction equation , 1975 .

[27]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[28]  G. Dahlquist.  Convergence and stability in the numerical integration of ordinary differential equations , 1956 .