A generic and efficient Taylor series based continuation method using a quadratic recast of smooth nonlinear systems

International audienceThis paper is concerned with a Taylor series based continuation algorithm, ie, the so-called Asymptotic Numerical Method (ANM). It describes a generic continuation procedure that apply the ANM principle at best, that is to say, that presents a high level of genericity without paying the price of this genericity by low computing performances. The way to quadratically recast a system of equation is now part of the method itself, and the way to handle elementary transcendental function is detailed with great attention. A sparse tensorial formalism is introduced for the internal representation of the system, which, when combines with a block condensation technique, provides a good computational efficiency of the ANM. Three examples are developed to show the performance and the versatility of the implementation of the continuation tool. Its robustness and its accuracy are explored. Finally, the potentiality of this method for complex non linear finite element analysis is enlightened by treating 2D elasticity problem with geometrical nonlinearities