TY - GEN

T1 - Coordinate perturbation method for upper bound limit analysis

AU - Hambleton, J. P.

AU - Sloan, S. W.

N1 - Copyright:
Copyright 2011 Elsevier B.V., All rights reserved.

PY - 2011

Y1 - 2011

N2 - Conventional methods for numerical upper bound limit analysis utilize a predefined mesh of elements with nodes at fixed coordinates, and there the objective is to find the kinematically admissible velocity field associated with the mesh that minimizes the limit load. In this paper, a numerical method for optimizing the applied load with respect to both the velocity field and nodal coordinates is presented. The formulation is for the classical case of rigid elements (blocks) separated by velocity discontinuities, considering plane strain and weightless material obeying the Mohr-Coulomb yield condition. The central idea of the approach is to successively perturb velocities and nodal coordinates starting from an initial solution corresponding to a predefined mesh, where the equations governing the perturbed velocities and nodal coordinates are derived through linearization of exact expressions of the limit load and the conditions imposed on velocity jumps by the flow rule. Despite approximations made in perturbation steps, the method furnishes a final solution which satisfies the flow rule and energy balance exactly (within numerical tolerances) and therefore provides a rigorous upper bound on the true collapse load. The effectiveness of the method is demonstrated by comparing numerical results with the analytical solution for a benchmark problem in soil mechanics.

AB - Conventional methods for numerical upper bound limit analysis utilize a predefined mesh of elements with nodes at fixed coordinates, and there the objective is to find the kinematically admissible velocity field associated with the mesh that minimizes the limit load. In this paper, a numerical method for optimizing the applied load with respect to both the velocity field and nodal coordinates is presented. The formulation is for the classical case of rigid elements (blocks) separated by velocity discontinuities, considering plane strain and weightless material obeying the Mohr-Coulomb yield condition. The central idea of the approach is to successively perturb velocities and nodal coordinates starting from an initial solution corresponding to a predefined mesh, where the equations governing the perturbed velocities and nodal coordinates are derived through linearization of exact expressions of the limit load and the conditions imposed on velocity jumps by the flow rule. Despite approximations made in perturbation steps, the method furnishes a final solution which satisfies the flow rule and energy balance exactly (within numerical tolerances) and therefore provides a rigorous upper bound on the true collapse load. The effectiveness of the method is demonstrated by comparing numerical results with the analytical solution for a benchmark problem in soil mechanics.

UR - http://www.scopus.com/inward/record.url?scp=80052106525&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=80052106525&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:80052106525

SN - 9789609875011

T3 - Computational Geomechanics, COMGEO II - Proceedings of the 2nd International Symposium on Computational Geomechanics

SP - 373

EP - 384

BT - Computational Geomechanics, COMGEO II - Proceedings of the 2nd International Symposium on Computational Geomechanics

T2 - 2nd International Symposium on Computational Geomechanics, COMGEO II

Y2 - 27 April 2011 through 29 April 2011

ER -