A modified Lemke Algorithm for dynamic rigid plastic response of skeletal structures. (November 2021)

• Record Type: Journal Article
• Authors: Khan, Azam; Ahmad, Irshad; Khan, Wajid; Khan, Niaz Bahadur; Aslam, Muhammad
• Is Part Of: Computers & structures. Volume 256(2021)
• Publication Date: 2021-11
• Subjects: Rigid-plastic response -- Lemke Algorithm -- Non-holonomic behavior -- Linear complementarity problem -- Quadratic program -- Degeneracy; Structural engineering -- Data processing -- Periodicals; Electronic data processing -- Structures, Theory of -- Periodicals; 624.171
• DOI: 10.1016/j.compstruc.2021.106638
• Languages: English
• ISSNs: 0045-7949

Abstract

Highlights: Modified Lemke Algorithm is used for dynamic rigid plastic response. Finite element formulation is in the form of Linear Complementarity Problem. The algorithm is adapted to handle strongly nonlinear rigid-plastic systems. The algorithm output for a pulse loaded simply supported beam is investigated. Abstract: This paper proposes a modified Lemke algorithm to determine the non-holonomic response of rigid plastic skeletal structures subjected to extreme dynamic loading. The basic formulation for the dynamic rigid-plastic response has a mathematical form of linear complementarity problem (LCP), and, equivalently, a pair of dual quadratic programs (QP). Previous attempts using the Wolfe type LCP solver exhibited very pronounced sensitivity with semidefinite LCPs, which appear in this class of dynamics. Therefore, the current study offers a numerically robust Lemke algorithm, which is considerably adapted to deal effectively with the presence of unrestricted variables and the semidefinite characteristics of structural matrices. In addition, degeneracy in the LCP solution of rigid plastic structures is efficiently handled by using the Lexicographic Minimum Ratio test. To validate the formulation and the algorithm developed, bench tests are conducted, which involve a simply supported beam and a portal frame subjected to a uniformly distributed rectangular pulse loading.