A solver combining reduced basis and convergence acceleration with applications to non‐linear elasticity

An iterative solver is proposed to solve the family of linear equations arising from the numerical computation of non-linear problems. This solver relies on two quantities coming from previous steps of the computations: the preconditioning matrix is a matrix that has been factorized at an earlier step and previously computed vectors yield a reduced basis. The principle is to define an increment in two sub-steps. In the first sub-step, only the projection of the unknown on a reduced subspace is incremented and the projection of the equation on the reduced subspace is satisfied exactly. In the second sub-step, the full equation is solved approximately with the help of the preconditioner. Last, the convergence of the sequences is accelerated by a well-known method, the modified minimal polynomial extrapolation. This algorithm assessed by classical benchmarks coming from shell buckling analysis. Finally, its insertion in path following techniques is discussed. This leads to non-linear solvers with few matrix factorizations and few iterations. Copyright © 2009 John Wiley & Sons, Ltd.

[1]  Claude Brezinski,et al.  Convergence acceleration during the 20th century , 2000 .

[2]  Khalide Jbilou,et al.  On block minimal residual methods , 2007, Appl. Math. Lett..

[3]  L. Duigou,et al.  Convergence acceleration of iterative algorithms. Applications to thin shell analysis and Navier–Stokes equations , 2009 .

[4]  H. Sadok,et al.  Vector extrapolation methods. applications and numerical comparison , 2000 .

[5]  J. Awrejcewicz,et al.  Asymptotic Solution of the Theory of Shells Boundary Value Problem , 2007 .

[6]  S. Liao,et al.  Beyond Perturbation: Introduction to the Homotopy Analysis Method , 2003 .

[7]  B. Cochelin,et al.  A critical review of asymptotic numerical methods , 1998 .

[8]  Michel Potier-Ferry,et al.  Mathematical and numerical connections between polynomial extrapolation and Padé approximants: applications in structural mechanics , 2004 .

[9]  G. Dhatt,et al.  Modélisation des structures par éléments finis , 1990 .

[10]  Michel Potier-Ferry,et al.  An iterative method based upon Padé approximants , 1999 .

[11]  E. Riks.  Some computational aspects of the stability analysis of nonlinear structures , 1984 .

[12]  Ji-Huan He,et al.  Comparison of homotopy perturbation method and homotopy analysis method , 2004, Appl. Math. Comput..

[13]  Michel Potier-Ferry,et al.  An iterative process based on homotopy and perturbation techniques , 2000 .

[14]  Michel Potier-Ferry,et al.  Projection techniques to improve high-order iterative correctors , 2004 .

[15]  D. O’Leary.  The block conjugate gradient algorithm and related methods , 1980 .

[16]  T. Weiland,et al.  Subspace projection extrapolation scheme for transient field simulations , 2004, IEEE Transactions on Magnetics.

[17]  T. Steinmetz,et al.  Multiple right-hand side techniques for the numerical simulation of quasistatic electric and magnetic fields , 2008 .

[18]  Hamid Zahrouni,et al.  Computing finite rotations of shells by an asymptotic-numerical method , 1999 .

[19]  Frédéric Guyomarc'h,et al.  A Deflated Version of the Conjugate Gradient Algorithm , 1999, SIAM J. Sci. Comput..

[20]  K. Y. Sze,et al.  Popular benchmark problems for geometric nonlinear analysis of shells , 2004 .

[21]  Bruno Cochelin,et al.  Asymptotic-numerical methods and pade approximants for non-linear elastic structures , 1994 .

[22]  Faezeh Toutounian,et al.  Global least squares method (Gl-LSQR) for solving general linear systems with several right-hand sides , 2006, Appl. Math. Comput..

[23]  E. Ramm,et al.  Shell theory versus degeneration—a comparison in large rotation finite element analysis , 1992 .