Computing Limit Loads by Minimizing a Sum of Norms

This paper treats the problem of computing the collapse state in limit analysis for a solid with a quadratic yield condition, such as, for example, the von Mises condition. After discretization with the finite element method, using divergence-free elements for the plastic flow, the kinematic formulation reduces to the problem of minimizing a sum of Euclidean vector norms, subject to a single linear constraint. This is a nonsmooth minimization problem, since many of the norms in the sum may vanish at the optimal point. Recently an efficient solution algorithm has been developed for this particular convex optimization problem in large sparse form. The approach is applied to test problems in limit analysis in two different plane models: plane strain and plates. In the first case more than 80% of the terms in the objective function are zero in the optimal solution, causing extreme ill conditioning. In the second case all terms are nonzero. In both cases the method works very well, and problems are solved which are larger by at least an order of magnitude than previously reported. The relative accuracy for the solution of the discrete problems, measured by duality gap and feasibility, is typically of the order 10-8.

[1]  Michael L. Overton,et al.  A quadratically convergent method for minimizing a sum of euclidean norms , 1983, Math. Program..

[2]  E. Anderheggen,et al.  Finite element limit analysis using linear programming , 1972 .

[3]  E. Christiansen.  Computation of limit loads , 1981 .

[4]  A. Conn,et al.  A second-order method for solving the continuous multifacility location problem , 1982 .

[5]  E. Christiansen.  Limit analysis of collapse states , 1996 .

[6]  E. Christiansen,et al.  Limit analysis in plasticity as a mathematical programming problem , 1980 .

[7]  Michael A. Saunders,et al.  A projected Lagrangian algorithm and its implementation for sparse nonlinear constraints , 1982 .

[8]  R. Temam.  Navier-Stokes Equations , 1977 .

[9]  P. G. Ciarlet,et al.  Basic error estimates for elliptic problems , 1991 .

[10]  A. Conn,et al.  A Stable Algorithm for Solving the Multifacility Location Problem Involving Euclidean Distances , 1980 .

[11]  Knud D. Andersen.  An Efficient Newton Barrier Method for Minimizing a Sum of Euclidean Norms , 1993, SIAM J. Optim..

[12]  Paul H. Calamai,et al.  A projected newton method forlp norm location problems , 1987, Math. Program..

[13]  H. G. Petersen,et al.  Estimation of convergence orders in repeated richardson extrapolation , 1989 .

[14]  K. O. Kortanek,et al.  Computing material collapse displacement fields on a cray X-MP/48 by the LP primal affine scaling algorithm , 1990 .

[15]  I. F. Collins,et al.  Limit analysis of plates and isoperimetric inequalities , 1994, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[16]  D. Goldfarb.  Extension of Davidon’s Variable Metric Method to Maximization Under Linear Inequality and Equality Constraints , 1969 .

[17]  E. Christiansen.  Limit Analysis for Plastic Plates , 1980 .

[18]  Knud D. Andersen.  An infeasible dual affine scaling method for linear programming , 1993 .

[19]  E. Christiansen,et al.  Computation of the collapse state in limit analysis using the LP primal affine scaling algorithm , 1991 .

[20]  Edmund Christiansen,et al.  Computations in limit analysis for plastic plates , 1983 .

[21]  Knud D. Andersen,et al.  Limit Analysis with the Dual Affine Scaling Algorithm , 1995 .

[22]  Edmund Christiansen,et al.  On the collapse solution in limit analysis , 1986 .

[23]  Michele Capurso.  Limit analysis of continuous media with piecewise linear yeld condition , 1971 .