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

This 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 stable and optimally convergent LaTIn-CutFEM algorithm for multiple unilateral contact problems

In this paper, we propose a novel unfitted finite element method for the simulation of multiple body contact. The computational mesh is generated independently of the geometry of the interacting solids, which can be arbitrarily complex. The key novelty of the approach is the combination of elements of the CutFEM technology, namely the enrichment of the solution field via the definition of overlapping fictitious domains with a dedicated penalty-type regularisation of discrete operators, and the LaTIn hybrid-mixed formulation of complex interface conditions. Furthermore, the novel P1-P1 discretisation scheme that we propose for the unfitted LaTIn solver is shown to be stable, robust and optimally convergent with mesh refinement. Finally, the paper introduces a high-performance 3D level-set/CutFEM framework for the versatile and robust solution of contact problems involving multiple bodies of complex geometries, with more than two bodies interacting at a single point

Regularized regressions for parametric models based on separated representations

Regressions created from experimental or simulated data enable the construction of metamodels, widely used in a variety of engineering applications. Many engineering problems involve multi-parametric physics whose corresponding multi-parametric solutions can be viewed as a sort of computational vademecum that, once computed offline, can be then used in a variety of real-time engineering applications including optimization, inverse analysis, uncertainty propagation or simulation based control. Sometimes, these multi-parametric problems can be solved by using advanced model order reduction—MOR-techniques. However, solving these multi-parametric problems can be very costly. In that case, one possibility consists in solving the problem for a sample of the parametric values and creating a regression from all the computed solutions. The solution for any choice of the parameters is then inferred from the prediction of the regression model. However, addressing high-dimensionality at the low data limit, ensuring accuracy and avoiding overfitting constitutes a difficult challenge. The present paper aims at proposing and discussing different advanced regressions based on the proper generalized decomposition (PGD) enabling the just referred features. In particular, new PGD strategies are developed adding different regularizations to the s-PGD method. In addition, the ANOVA-based PGD is proposed to ally them

Méthode des éléments finis : Développements récents adaptés au calcul de structures complexes

National audienceLe développement et la commercialisation de nombreux codes de calcul par la méthode des éléments finis rendent la simulation de systèmes mécaniques accessible à tous, même aux plus novices d'entre nous. Les systèmes de modélisation et de CAO permettent aisément d'obtenir un modèle éléments finis, même pour des géométries complexes. De plus, les systèmes de post-traitement graphiques tracent des cartes d'isovaleurs (de contraintes, de déformations, …) avec une qualité et un nombre de couleurs impressionnants. Néanmoins il convient de rester lucide et d'être conscients de la qualité des calculs que nous menons. Même si elle se banalise, la méthode des éléments finis reste un outil complexe et savoir bien la maîtriser n'est pas chose aisée. Nous présentons ici quelques notions sur le contrôle des calculs et sur certains outils en cours de développement pour donner à l'utilisateur l'assurance d'une qualité de son modèle. Un exemple de maîtrise de modèle élastique 3D est détaillé. De plus, nous essayons de donner quelques notions sur les difficultés rencontrées lors de l'étude de structures complexes qui n'est pas encore, de nos jours, à la portée de tous. Par "complexe", nous entendons complexité du système mécanique (géométrie, liaisons, …), complexité des sollicitations appliquées (chargements divers, dynamiques, cycliques, …, grands déplacements, …) et complexité du modèle de matériau utilisé (plasticité, viscoplasticité, endommagement, …).Nous présentons les problèmes numériques rencontrés lors de l'étude, par la méthode des éléments finis, de structures complexes et nous donnons quelques pistes pour faciliter le traitement de certains types de problèmes. Un exemple de méthode adaptée au calcul d'assemblages complexes de structures élastiques est proposé. Après avoir donné quelques idées sur les problèmes de taille et de coût de calcul, nous présentons les perspectives offertes par les calculateurs à architecture parallèle modernes. Un exemple illustre ces possibilités

Sur la prise en compte des variations de chargement et de géométrie dans les assemblages

L'objet de ce travail est de développer une méthode spécialisée pour l'étude paramétrique d'assemblages incluant de multiples zones de contact frottant. Cette méthode tire parti du fait que les non-linéarités et/ou les incertitudes sont localisées dans les liaisons. On utilise une représentation unifiée de ces zones de liaison sous forme d'interfaces en introduisant une décomposition de l'assemblage. L'algorithme de résolution itératif utilisé (méthode LATIN) permet un découplage du traitement des non-linéarités et/ou incertitudes locales du traitement des problèmes linéaires globaux. Chaque configuration d'un assemblage correspond à un choix de valeurs des paramètres décrivant les chargements extérieurs et les géométries des zones de contact. Plutôt que de réaliser un calcul complet pour chaque jeu de paramètres, nous utilisons la capacité de la méthode LATIN à réutiliser simplement la solution d'un problème donné (associé à un jeu de paramètres) pour en résoudre d'autres. Les exemples proposés concernent des assemblages de structures 3D avec des liaisons par contact frottant

Physics-informed and graph neural networks for enhanced inverse analysis

Purpose : This paper presents an original approach for learning models, partially known, of particular interest when performing source identification or structural health monitoring. The proposed procedures employ some amount of knowledge on the system under scrutiny as well as a limited amount of data efficiently assimilated. Design/methodology/approach : Two different formulations are explored. The first, based on the use of informed neural networks, leverages data collected at specific locations and times to determine the unknown source term of a parabolic partial differential equation. The second procedure, more challenging, involves learning the unknown model from a single measured field history, enabling the localization of a region where material properties differ. Findings: Both procedures assume some kind of sparsity, either in the source distribution or in the region where physical properties differ. This paper proposed two different neural approaches able to learn models in order to perform efficient inverse analyses.Originality/valueTwo original methodologies are explored to identify hidden property that can be recovered with the right usage of data. Both methodologies are based on neural network architecture. Originality/value : Two original methodologies are explored to identify hidden property that can be recovered with the right usage of data. Both methodologies are based on neural network architecture

Port-metriplectic neural networks: thermodynamics-informed machine learning of complex physical systems

We develop inductive biases for the machine learning of complex physical systems based on the port-Hamiltonian formalism. To satisfy by construction the principles of thermodynamics in the learned physics (conservation of energy, non-negative entropy production), we modify accordingly the port-Hamiltonian formalism so as to achieve a port-metriplectic one. We show that the constructed networks are able to learn the physics of complex systems by parts, thus alleviating the burden associated to the experimental characterization and posterior learning process of this kind of systems. Predictions can be done, however, at the scale of the complete system. Examples are shown on the performance of the proposed technique.Comment: 9 pages, 5 figure

Recalage de structures légères par une approximation polynomiale en vue de leur contrôle actif

Le recalage de structures légères en vue du contrôle de leurs vibrations est une problématique importante. Un des outils mécaniques performant dans ce domaine est le recalage basé sur l'erreur en relation de comportement. Cet outil est adapté ici au contrôle actif en contexte incertain de telles structures grâce une description polynomiale de l'algorithme de calcul : les inconnues des polynômes sont la variabilité des paramètres (matériau ou conditions aux limites) qu'on souhaite recaler

Modular parametric PGD enabling online solution of partial differential equations

In the present work, a new methodology is proposed for building surrogate parametric models of engineering systems based on modular assembly of pre-solved modules. Each module is a generic parametric solution considering parametric geometry, material and boundary conditions. By assembling these modules and satisfying continuity constraints at the interfaces, a parametric surrogate model of the full problem can be obtained. In the present paper, the PGD technique in connection with NURBS geometry representation is used to create a parametric model for each module. In this technique, the NURBS objects allow to map the governing boundary value problem from a parametric non-regular domain into a regular reference domain and the PGD is used to create a reduced model in the reference domain. In the assembly stage, an optimization problem is solved to satisfy the continuity constraints at the interfaces. The proposed procedure is based on the offline--online paradigm: the offline stage consists of creating multiple pre-solved modules which can be afterwards assembled in almost real-time during the online stage, enabling quick evaluations of the full system response. To show the potential of the proposed approach some numerical examples in heat conduction and structural plates under bending are presented

Learning data-driven reduced elastic and inelastic models of spot-welded patches

Solving mechanical problems in large structures with rich localized behaviors remains a challenging issue despite the enormous advances in numerical procedures and computational performance. In particular, these localized behaviors need for extremely fine descriptions, and this has an associated impact in the number of degrees of freedom from one side, and the decrease of the time step employed in usual explicit time integrations, whose stability scales with the size of the smallest element involved in the mesh. In the present work we propose a data-driven technique for learning the rich behavior of a local patch and integrate it into a standard coarser description at the structure level. Thus, localized behaviors impact the global structural response without needing an explicit description of that fine scale behaviors