Low intrusive coupling of implicit and explicit time integration schemes for structural dynamics

Simulation of low energy impacts on composite structures is a key feature in aeronautics. Unfortunately it involves very expensive numerical simulations: on the one side, the structures of interest have large dimensions and need fine volumic meshes (at least locally) in order to properly capture damage. On the other side, explicit simulations are commonly used to lead this kind of simulations (Lopes et al., 2009 1; Bouvet, 2009 2), which results in very small time steps to ensure the CFL condition (Courant et al., 1967 3). Implicit algorithms are actually more difficult to use in this situation because of the lack of smoothness of the solution that can lead to prohibitive number of time steps or even to non-convergence of Newton-like iterative processes. It is also observed that non-smooth phenomena are localized in space and time (near the impacted zone). It may therefore be advantageous to adopt a multiscale space/time approach by splitting the structure into several substructures with their own space/time discretization and their own integration scheme. The purpose of this decomposition is to take advantage of the specificities of both algorithms families: explicit scheme focuses on non-smooth areas while smoother parts (actually linear in this work) of the solutions are computed with larger time steps with an implicit scheme. We propose here an implementation of the Gravouil-Combescure method (GC) (Combescure and Gravouil, 2002 4) by the mean of low intrusive coupling between the implicit finite element analysis (FEA) code Zset/Zebulon (Z-set official website, 2013 5) and the explicit FEA code Europlexus (Europlexus official website, 2013 6). Simulations of low energy impacts on composite stiffened panels are presented. It is shown on this application that large time step ratios can be reached, thus saving computation time. HighlightsNon-intrusive GC implementation.Industrial code-coupling.Simulation of low energy impact on laminated stiffened panel.Python interface to finite element codes.

[1]  Alain Combescure,et al.  Damage model with delay effect Analytical and numerical studies of the evolution of the characteristic damage length , 2003 .

[2]  O. Zienkiewicz,et al.  A note on mass lumping and related processes in the finite element method , 1976 .

[3]  Olivier Allix,et al.  A two‐scale approximation of the Schur complement and its use for non‐intrusive coupling , 2011 .

[4]  Charbel Farhat,et al.  A transient FETI methodology for large‐scale parallel implicit computations in structural mechanics , 1994 .

[5]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

[6]  Christophe Bouvet,et al.  Validation of low velocity impact modelling on different stacking sequences of CFRP laminates and influence of fibre failure , 2013 .

[7]  Alain Combescure,et al.  A numerical scheme to couple subdomains with different time-steps for predominantly linear transient analysis , 2002 .

[8]  S. Krenk,et al.  Multi‐time‐step domain coupling method with energy control , 2010 .

[9]  O. Allix,et al.  Delayed-Damage Modelling for Fracture Prediction of Laminated Composites under Dynamic Loading , 1997 .

[10]  Zafer Gürdal,et al.  LOW-VELOCITY IMPACT DAMAGE ON DISPERSED STACKING SEQUENCE LAMINATES. PART I: EXPERIMENTS , 2009 .

[11]  Sylvain Gavoille Stratégies de couplage de modèles discret-continu en dynamique explicite , 2009 .

[12]  Thomas J. R. Hughes,et al.  Improved numerical dissipation for time integration algorithms in structural dynamics , 1977 .

[13]  J.-F. Maire,et al.  A multiscale progressive failure approach for composite laminates based on thermodynamical viscoelastic and damage models , 2007 .

[14]  Alain Combescure,et al.  A monolithic energy conserving method to couple heterogeneous time integrators with incompatible time steps in structural dynamics , 2011 .

[15]  Vincent Faucher,et al.  A time and space mortar method for coupling linear modal subdomains and non-linear subdomains in explicit structural dynamics , 2003 .

[16]  R. Courant,et al.  On the Partial Difference Equations, of Mathematical Physics , 2015 .

[17]  O. C. Zienkiewicz,et al.  An alpha modification of Newmark's method , 1980 .

[18]  A. Prakash,et al.  A FETI‐based multi‐time‐step coupling method for Newmark schemes in structural dynamics , 2004 .

[19]  Zafer Gürdal,et al.  Low-velocity impact damage on dispersed stacking sequence laminates. Part II: Numerical simulations , 2009 .

[20]  Jintai Chung,et al.  Explicit time integration algorithms for structural dynamics with optimal numerical dissipation , 1996 .

[21]  C. Farhat,et al.  A method of finite element tearing and interconnecting and its parallel solution algorithm , 1991 .

[22]  Olivier Allix,et al.  A 3D damage analysis of low-velocity impacts on laminated composites , 2002 .

[23]  Steen Krenk,et al.  Extended state‐space time integration with high‐frequency energy dissipation , 2008 .

[24]  Thomas Hauser,et al.  Recent Updates to the CFD General Notation System (CGNS) , 2012 .

[25]  Chloé Dupleix-Couderc Sur une stratégie multi-échelle d'analyse des grands délaminages en dynamique transitoire , 2011 .

[26]  P. Gosselet,et al.  Non-overlapping domain decomposition methods in structural mechanics , 2006, 1208.4209.

[27]  Christophe Bouvet,et al.  Low velocity impact modeling in composite laminates capturing permanent indentation , 2012 .

[28]  A. Peirce Computer Methods in Applied Mechanics and Engineering , 2010 .

[29]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[30]  Olivier Allix,et al.  A coupling strategy for adaptive local refinement in space and time with a fixed global model in explicit dynamics , 2014 .

[31]  Alain Combescure,et al.  Arlequin framework for multi-model, multi-time scale and heterogeneous time integrators for structural transient dynamics , 2013 .

[32]  O. Allix,et al.  Non-intrusive and exact global/local techniques for structural problems with local plasticity , 2009 .