Print Email Facebook Twitter Two solution strategies to improve the computational performance of Sequentially Linear Analysis for quasi-brittle structures Title Two solution strategies to improve the computational performance of Sequentially Linear Analysis for quasi-brittle structures Author Pari, M. (TU Delft Applied Mechanics) Swart, W. (TU Delft Electrical Engineering, Mathematics and Computer Science) van Gijzen, M.B. (TU Delft Numerical Analysis) Hendriks, M.A.N. (TU Delft Applied Mechanics) Rots, J.G. (TU Delft Applied Mechanics) Faculty Electrical Engineering, Mathematics and Computer Science Date 2019 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. Subject direct linear solveriterative linear solverlow-rank matrix correctionnonlinear finite element analysissequentially linear analysis To reference this document use: http://resolver.tudelft.nl/uuid:f2d72f14-6c68-4058-a825-12aaf78d0a21 DOI https://doi.org/10.1002/nme.6302 ISSN 0029-5981 Source International Journal for Numerical Methods in Engineering, 121 (2020) (10), 2128-2146 Part of collection Institutional Repository Document type journal article Rights © 2019 M. Pari, W. Swart, M.B. van Gijzen, M.A.N. Hendriks, J.G. Rots Files PDF nme.6302.pdf 2.78 MB Close viewer /islandora/object/uuid:f2d72f14-6c68-4058-a825-12aaf78d0a21/datastream/OBJ/view