On a multiscale strategy and its optimization for the simulation of combined delamination and buckling

Abstract

International audienceThis paper investigates a computational strategy for studying the interactions between multiple through-the-width delaminations and global or local buckling in composite laminates taking into account possible contact between the delaminated surfaces. In order to achieve an accurate prediction of the quasi-static response, a very refined discretization of the structure is required, leading to the resolution of very large and highly nonlinear numerical problems. In this paper, a nonlinear finite element formulation along with a parallel iterative scheme based on a multiscale domain decomposition are used for the computation of 3D mesoscale models. Previous works by the authors already dealt with the simulation of multiscale delamination assuming small perturbations. This paper presents the formulation used to include geometric nonlinearities into this existing multiscale framework and discusses the adaptations that need to be made to the iterative process in order to ensure the rapid convergence and the scalability of the method in the presence of buckling and delamination. These various adaptations are illustrated by simulations involving large numbers of DOFs

A parallel noninvasive multiscale strategy for a mixed domain decomposition method with frictional contact

“…Remark The choice of the search direction is always a complex question since a good search direction should be able to provide information on the stiffness of the whole structure (even if the multiscale approach enables to correct a bad initial guess), but this information may be dependent on local information and in, particular, on the contact state. In the work of Saavedra et al, it was proposed to regularly update the search direction according to the current contact status. Another work by the aforementioned authors illustrated the many parameters that affected the convergence (geometry, heterogeneity, load case, partitioning, search direction, macrobasis).…”

“…In the work of Saavedra et al, 40 it was proposed to regularly update the search direction according to the current contact status. Another work by the aforementioned authors 34 illustrated the many parameters that affected the convergence (geometry, heterogeneity, load case, partitioning, search direction, macrobasis).…”

Substructured formulations of nonlinear structure problems - influence of the interface condition

“…Studies have proved that they needed to reproduce short-range effects (like local heterogeneity) but also structural effects (like the anisotropy induced by the slenderness of plate structures [25]). When one wishes to choose an invariant scalar (or tensor in case of anisotropy) for each interface, it can be beneficial to use a coarse model for its estimation [26].…”

Mixed nonlinear localization strategy for the buckling and post-buckling analyses of structures: parameter optimization and path following implementation

“…A dual version has been applied to damage analysis in Pebrel, Rey and Gosselet (2008). A mixed version has been applied to the delamination of composites structures in Allix, Kerfriden and Gosselet (2010a) and to delamination-buckling interaction in Saavedra, Allix and Gosselet (2012).…”