The paper describes the derivation and application of a range of numerical algorithms for implementing the Mohr—Coulomb yield criterion in a non‐linear finite element computer program. Emphasis is placed on the difficulties associated with the corners of the yield surface. In contrast to the more conventional forward‐Euler procedures, a backward‐Euler integration technique is adopted. A range of methods, including a ‘consistent approach’ are used to derive the tangent modular matrix. Numerical experiments are presented which involve solution algorithms including the modified and full Newton—Raphson procedures, ‘line‐searches’ and the arc‐length method. It is shown that the introduction of efficient integration and tangency algorithms can lead to very substantial improvements in the convergence characteristics.
CitationDownload as .RIS
MCB UP Ltd
Copyright © 1987, MCB UP Limited