Premium
Frictionless geometrically non‐linear contact using quadratic programming
Author(s) -
Johnson A. R.,
Quigley C. J.
Publication year - 1989
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1620280110
Subject(s) - sequential quadratic programming , finite element method , quadratic programming , mathematics , quadratic equation , taylor series , displacement (psychology) , simplex algorithm , mathematical analysis , mathematical optimization , linear programming , geometry , structural engineering , engineering , psychology , psychotherapist
Abstract The successive quadratic programming (SQP) method is used with the finite element method (FEM) to solve frictionless geometrically non‐linear contact problems involving large deformations of the elastica in the presence of flat rigid walls. To formulate the SQP problems, the potential energy (PE) is expanded in a Taylor series of second order in displacement increments about a configuration near a contact solution. The SQP problems consist of minimizing the Taylor expansion of the PE subject to the inequality constraints which represent contact. The quadratic programming (QP) method is made part of a Newton–Raphson (NR) search in which the QP corrections are made when a NR step does not satisfy the constraints. A revised simplex method developed by Rusin is used to solve the QP problems. The elastica is modelled with a total Lagrangian FEM developed by Fried. Solutions are obtained for the end loaded buckled elastica in point contact with a rigid wall and for a uniformly loaded elastica in regional contact with a rigid wall. The problems are also solved using a penalty method. The results obtained for the point contact problem are compared to an analytical solution. Calculations were made to obtain numerical information on maximum load step size and the number of inverse operations required for each load step. Cases in which the elastica stiffened substantially as a result of the initiation of contact are also discussed.