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

Manimaran Pari, Wouter Swart, Martin van Gijzen, Max Hendriks, Jan Rots

Research output: Contribution to journalArticleScientificpeer-review

7 Citations (Scopus)
58 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, 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 minimised. Comparison and validation against a traditional parallel direct sparse solver, with regard to a two-dimensional (2D) and three-dimensional (3D) benchmark study, illustrates the improved performance of the Woodbury-based direct solver over its counterparts, especially for large 3D problems.

Original languageEnglish
Pages (from-to)2128-2146
Number of pages19
JournalInternational Journal for Numerical Methods in Engineering
Volume121 (2020)
Issue number10
DOIs
Publication statusPublished - 2019

Keywords

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

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