Two solution strategies to improve the computational performance of Sequentially Linear Analysis for quasi-brittle structures

Research output: Contribution to journalArticleScientificpeer-review

22 Downloads (Pure)

Abstract

Sequentially Linear Analysis (SLA), an event-by-event procedure for finite element (FE) simulation of quasi-brittle materials, is based on sequentially identifying a critical integration point in the FE model, to reduce its strength and stiffness, and the corresponding critical load multiplier (λcrit ), to scale the linear analysis results. In this article, two strategies are proposed to efficiently reuse previous stiffness matrix factorisations and their corresponding solutions in subsequent linear analyses, since the global system of linear equations representing the FE model changes only locally. The first is based on a direct solution method in combination with the Woodbury matrix identity, to compute the inverse of a low-rank corrected stiffness matrix relatively cheaply. The second is a variation of the traditional incomplete LU preconditioned conjugate gradient method (CG), wherein the preconditioner is the complete factorisation of a previous analysis step’s stiffness matrix. For both the approaches, optimal points at which the factorisation is recomputed are determined such that the total analysis time is minimized. Comparison and validation against a traditional parallel direct sparse solver, with regard to a 2D and 3D benchmark study, illustrates the improved performance of the Woodbury based direct solver over its counterparts, especially for large 3D problems.
Original languageEnglish
Number of pages19
JournalInternational Journal for Numerical Methods in Engineering
DOIs
Publication statusPublished - 2019

Keywords

  • nonlinear finite element analysis
  • sequentially linear analysis
  • direct linear solver
  • iterative linear solver
  • low‐rank matrix correction

Fingerprint Dive into the research topics of 'Two solution strategies to improve the computational performance of Sequentially Linear Analysis for quasi-brittle structures'. Together they form a unique fingerprint.

Cite this