A minimally-intrusive fully 3D separated plate formulation in computational structural mechanics

Most of mechanical systems and complex structures exhibit plate and shell components. Therefore, 2D simulation, based on plate and shell theory, appears as an appealing choice in structural analysis as it allows reducing the computational complexity. Nevertheless, this 2D framework fails for capturing rich physics compromising the usual hypotheses considered when deriving standard plate and shell theories. To circumvent, or at least alleviate this issue, authors proposed in their former works an in-plane–out-of-plane separated representation able to capture rich 3D behaviors while keeping the computational complexity the one of 2D simulations. In the present paper we propose an efficient integration of fully 3D descriptions into existing plate software

Efficient simulation tools for real-time monitoring and control using model order reduction and data-driven techniques

Numerical simulation, the use of computers to run a program which implements a mathematical model for a physical system, is an important part of today technological world. It is required in many scientific and engineering fields to study the behaviour of systems whose mathematical models are too complex to provide analytical solutions and it makes virtual evaluation of systems responses possible (virtual twins). This drastically reduces the number of experimental tests for accurate designs of the real system that the numerical model represents. However these virtual twins, based on classical methods which make use of a rich representations of the system (ex. finite element method), rarely allows real-time feedback, even when considering high performance computing, operating on powerful platforms. In these circumstances, the real-time performance required in some applications are compromised. Indeed the virtual twins are static, that is, they are used in the design of complex systems and their components, but they are not expected to accommodate or assimilate data so as to define dynamic data-driven application systems. Moreover significant deviations between the observed response and the one predicted by the model are usually noticed due to inaccuracy in the employed models, in the determination of the model parameters or in their time evolution. In this thesis we propose different methods to solve these handicaps in order to perform real-time monitoring and control. In the first part Model Order Reduction (MOR) techniques are used to accommodate real-time constraints; they compute a good approximation of the solution by simplifying the solution procedure instead of the model. The accuracy of the predicted solution is not compromised and efficient simulations can be performed (digital twins). In the second part Data-Driven modelling are employed to fill the gap between the parametric solution computed by using non-intrusive MOR techniques and the measured fields, in order to make dynamic data-driven application systems, DDDAS, possible (Hybrid Twins).La simulación numérica, el uso de ordenadores para ejecutar un programa que implementa un modelo matemático de un sistema físico, es una parte importante del mundo tecnológico actual. En muchos campos de la ciencia y la ingeniería es necesario estudiar el comportamiento de sistemas cuyos modelos matemáticos son demasiado complejos para proporcionar soluciones analíticas, haciendo posible la evaluación virtual de las respuestas de los sistemas (gemelos virtuales). Esto reduce drásticamente el número de pruebas experimentales para los diseños precisos del sistema real que el modelo numérico representa. Sin embargo, estos gemelos virtuales, basados en métodos clásicos que hacen uso de una rica representación del sistema (por ejemplo, el método de elementos finitos), rara vez permiten la retroalimentación en tiempo real, incluso cuando se considera la computación en plataformas de alto rendimiento. En estas circunstancias, el rendimiento en tiempo real requerido en algunas aplicaciones se ve comprometido. En efecto, los gemelos virtuales son estáticos, es decir, se utilizan en el diseño de sistemas complejos y sus componentes, pero no se espera que acomoden o asimilen los datos para definir sistemas de aplicación dinámicos basados en datos. Además, se suelen apreciar desviaciones significativas entre la respuesta observada y la predicha por el modelo, debido a inexactitudes en los modelos empleados, en la determinación de los parámetros del modelo o en su evolución temporal. En esta tesis se proponen diferentes métodos para resolver estas limitaciones con el fin de realizar un seguimiento y un control en tiempo real. En la primera parte se utilizan técnicas de Reducción de Modelos para satisfacer las restricciones en tiempo real; estas técnicas calculan una buena aproximación de la solución simplificando el procedimiento de resolución en lugar del modelo. La precisión de la solución no se ve comprometida y se pueden realizar simulaciones efficientes (gemelos digitales). En la segunda parte se emplea la modelización basada en datos para llenar el vacío entre la solución paramétrica, calculada utilizando técnicas de reducción de modelos no intrusivas, y los campos medidos, con el fin de hacer posibles los sistemas de aplicación dinámicos basados en datos (gemelos híbridos).La simulation numérique, c'est-à-dire l'utilisation des ordinateurs pour exécuter un programme qui met en oeuvre un modèle mathématique d'un système physique, est une partie importante du monde technologique actuel. Elle est nécessaire dans de nombreux domaines scientifiques et techniques pour étudier le comportement de systèmes dont les modèles mathématiques sont trop complexes pour fournir des solutions analytiques et elle rend possible l'évaluation virtuelle des réponses des systèmes (jumeaux virtuels). Cela réduit considérablement le nombre de tests expérimentaux nécessaires à la conception précise du système réel que le modèle numérique représente. Cependant, ces jumeaux virtuels, basés sur des méthodes classiques qui utilisent une représentation fine du système (ex. méthode des éléments finis), permettent rarement une rétroaction en temps réel, même dans un contexte de calcul haute performance, fonctionnant sur des plates-formes puissantes. Dans ces circonstances, les performances en temps réel requises dans certaines applications sont compromises. En effet, les jumeaux virtuels sont statiques, c'est-à-dire qu'ils sont utilisés dans la conception de systèmes complexes et de leurs composants, mais on ne s'attend pas à ce qu'ils prennent en compte ou assimilent des données afin de définir des systèmes d'application dynamiques pilotés par les données. De plus, des écarts significatifs entre la réponse observée et celle prévue par le modèle sont généralement constatés en raison de l'imprécision des modèles employés, de la détermination des paramètres du modèle ou de leur évolution dans le temps. Dans cette thèse, nous proposons di érentes méthodes pour résoudre ces handicaps afin d'effectuer une surveillance et un contrôle en temps réel. Dans la première partie, les techniques de Réduction de Modèles sont utilisées pour tenir compte des contraintes en temps réel ; elles calculent une bonne approximation de la solution en simplifiant la procédure de résolution plutôt que le modèle. La précision de la solution n'est pas compromise et des simulations e caces peuvent être réalisées (jumeaux numériquex). Dans la deuxième partie, la modélisation pilotée par les données est utilisée pour combler l'écart entre la solution paramétrique calculée, en utilisant des techniques de réduction de modèles non intrusives, et les champs mesurés, afin de rendre possibles des systèmes d'application dynamiques basés sur les données (jumeaux hybrides)

Efficient simulation tools for real-time monitoring and control using model order reduction and data-driven techniques

Cotutela: Universitat Politècnica de Catalunya i École Centrale de Nantes.Numerical simulation, the use of computers to run a program which implements a mathematical model for a physical system, is an important part of today technological world. It is required in many scientific and engineering fields to study the behaviour of systems whose mathematical models are too complex to provide analytical solutions and it makes virtual evaluation of systems responses possible (virtual twins). This drastically reduces the number of experimental tests for accurate designs of the real system that the numerical model represents. However these virtual twins, based on classical methods which make use of a rich representations of the system (ex. finite element method), rarely allows real-time feedback, even when considering high performance computing, operating on powerful platforms. In these circumstances, the real-time performance required in some applications are compromised. Indeed the virtual twins are static, that is, they are used in the design of complex systems and their components, but they are not expected to accommodate or assimilate data so as to define dynamic data-driven application systems. Moreover significant deviations between the observed response and the one predicted by the model are usually noticed due to inaccuracy in the employed models, in the determination of the model parameters or in their time evolution. In this thesis we propose different methods to solve these handicaps in order to perform real-time monitoring and control. In the first part Model Order Reduction (MOR) techniques are used to accommodate real-time constraints; they compute a good approximation of the solution by simplifying the solution procedure instead of the model. The accuracy of the predicted solution is not compromised and efficient simulations can be performed (digital twins). In the second part Data-Driven modelling are employed to fill the gap between the parametric solution computed by using non-intrusive MOR techniques and the measured fields, in order to make dynamic data-driven application systems, DDDAS, possible (Hybrid Twins).La simulación numérica, el uso de ordenadores para ejecutar un programa que implementa un modelo matemático de un sistema físico, es una parte importante del mundo tecnológico actual. En muchos campos de la ciencia y la ingeniería es necesario estudiar el comportamiento de sistemas cuyos modelos matemáticos son demasiado complejos para proporcionar soluciones analíticas, haciendo posible la evaluación virtual de las respuestas de los sistemas (gemelos virtuales). Esto reduce drásticamente el número de pruebas experimentales para los diseños precisos del sistema real que el modelo numérico representa. Sin embargo, estos gemelos virtuales, basados en métodos clásicos que hacen uso de una rica representación del sistema (por ejemplo, el método de elementos finitos), rara vez permiten la retroalimentación en tiempo real, incluso cuando se considera la computación en plataformas de alto rendimiento. En estas circunstancias, el rendimiento en tiempo real requerido en algunas aplicaciones se ve comprometido. En efecto, los gemelos virtuales son estáticos, es decir, se utilizan en el diseño de sistemas complejos y sus componentes, pero no se espera que acomoden o asimilen los datos para definir sistemas de aplicación dinámicos basados en datos. Además, se suelen apreciar desviaciones significativas entre la respuesta observada y la predicha por el modelo, debido a inexactitudes en los modelos empleados, en la determinación de los parámetros del modelo o en su evolución temporal. En esta tesis se proponen diferentes métodos para resolver estas limitaciones con el fin de realizar un seguimiento y un control en tiempo real. En la primera parte se utilizan técnicas de Reducción de Modelos para satisfacer las restricciones en tiempo real; estas técnicas calculan una buena aproximación de la solución simplificando el procedimiento de resolución en lugar del modelo. La precisión de la solución no se ve comprometida y se pueden realizar simulaciones efficientes (gemelos digitales). En la segunda parte se emplea la modelización basada en datos para llenar el vacío entre la solución paramétrica, calculada utilizando técnicas de reducción de modelos no intrusivas, y los campos medidos, con el fin de hacer posibles los sistemas de aplicación dinámicos basados en datos (gemelos híbridos).La simulation numérique, c'est-à-dire l'utilisation des ordinateurs pour exécuter un programme qui met en oeuvre un modèle mathématique d'un système physique, est une partie importante du monde technologique actuel. Elle est nécessaire dans de nombreux domaines scientifiques et techniques pour étudier le comportement de systèmes dont les modèles mathématiques sont trop complexes pour fournir des solutions analytiques et elle rend possible l'évaluation virtuelle des réponses des systèmes (jumeaux virtuels). Cela réduit considérablement le nombre de tests expérimentaux nécessaires à la conception précise du système réel que le modèle numérique représente. Cependant, ces jumeaux virtuels, basés sur des méthodes classiques qui utilisent une représentation fine du système (ex. méthode des éléments finis), permettent rarement une rétroaction en temps réel, même dans un contexte de calcul haute performance, fonctionnant sur des plates-formes puissantes. Dans ces circonstances, les performances en temps réel requises dans certaines applications sont compromises. En effet, les jumeaux virtuels sont statiques, c'est-à-dire qu'ils sont utilisés dans la conception de systèmes complexes et de leurs composants, mais on ne s'attend pas à ce qu'ils prennent en compte ou assimilent des données afin de définir des systèmes d'application dynamiques pilotés par les données. De plus, des écarts significatifs entre la réponse observée et celle prévue par le modèle sont généralement constatés en raison de l'imprécision des modèles employés, de la détermination des paramètres du modèle ou de leur évolution dans le temps. Dans cette thèse, nous proposons di érentes méthodes pour résoudre ces handicaps afin d'effectuer une surveillance et un contrôle en temps réel. Dans la première partie, les techniques de Réduction de Modèles sont utilisées pour tenir compte des contraintes en temps réel ; elles calculent une bonne approximation de la solution en simplifiant la procédure de résolution plutôt que le modèle. La précision de la solution n'est pas compromise et des simulations e caces peuvent être réalisées (jumeaux numériquex). Dans la deuxième partie, la modélisation pilotée par les données est utilisée pour combler l'écart entre la solution paramétrique calculée, en utilisant des techniques de réduction de modèles non intrusives, et les champs mesurés, afin de rendre possibles des systèmes d'application dynamiques basés sur les données (jumeaux hybrides).Postprint (published version

Simulation tools for biomechanical applications with PGD-based reduced order models

Cotutela Universitat Politècnica de Catalunya i Università degli Studi di PaviaNumerical simulation tools are generally used in all modern engineering fields, especially those having difficulties in performing large number of practical experiments, such as biomechanics. Among the computational methods, Finite Element (FE) is an essential tool. Nowadays, the fast-growing computational techniques, from the upgrading hardware to the emerging of novel algorithm, have already enabled extensive applications in biomechanics. For applications that require fast response and/or multiple queries, Reduced Order Modelling (ROM) methods have been developed based on existing methods such as FE, and have eventually enabled real-time numerical simulation for a large variety of engineering problems. In this thesis, several novel computational techniques are developed to explore the capability of Proper Generalised Decomposition (PGD), which is an important approach of ROM. To assess the usability of the PGD-based ROM for biomechanical applications, a real human femur bone is chosen to study its mechanical behaviour as an example. Standard image-based modelling procedure in biomechanics is performed to create an FE model which is then validated with in vitro experimental results. As a basis of this work, the medical image processing has to be performed, in order to generate an available FE model. This model is validated according to data collected from a previously performed \textit{in vitro} experimental test. The full procedure of image-based model generation and the validation of generated model is described in Chapter 2. As a major objective of this thesis, a non-intrusive scheme for the PGD framework is developed in Chapter 3. It is implemented using in-house developed Matlab (Mathworks, USA) code to conduct the PGD work flow, and calling Abaqus as an external solver for devised fictitious mechanical problems. The transformation of data from computed tomography (CT) image set to FE model including inhomogeneous material properties is subjected to some physical constraints, and when applying the load, there are also geometric constraints limiting the locations where load could be applied. These constraints will lead to a constrained parameter space, which possibly has difficulty to be separated in a Cartesian fashion. Therefore, a novel strategy to separate the parameters in a collective manner is proposed in Chapter 4. Chapter 5 details a comprehensive application in biomechanics, the methodologies proposed in Chapter 3 and 4 are applied on the practical model generated in Chapter 2. As a typical application of the PGD vademecum, a material property identification problem is discussed. Further PGD vademecum is generated using the identified material properties with variable loading locations, and with this vademecum, real-time mechanical response of the femur is available. In addition, for the purpose of extending the methodologies to orthotropic materials, which is commonly used in biomechanics, in Chapter 6 another linear elastic model is investigated with the non-intrusive PGD scheme. Nowadays, isogeometric analysis (IGA) is a very popular tool in computational mechanics. It is appealing to take advantage of non-uniform rational B-splines (NURBS) to discretise the model. For PGD, using B-splines for the discretisation of the parameter space could improve the quality of vademecum, especially for problems involving sensitivities with respect to the parameters during the online computations. It is important and necessary to extend the PGD framework to nonlinear solid mechanics, because most biological soft tissues have been observed nonlinear mechanical behaviours. Consequently, in Chapter 7 we have developed a PGD framework for the St.Venant-Kirchhoff constitutive model using the Picard linearisation which is consistent with the fixed-point iteration algorithm commonly used in PGD. In Chapter 8, conclusive remarks are addressed as well as forecasts of possible future works.Postprint (published version

On the coupling of local 3D solutions and global 2D shell theory in structural mechanics

Most of mechanical systems and complex structures exhibit plate and shell components. Therefore, 2D simulation, based on plate and shell theory, appears as an appealing choice in structural analysis as it allows reducing the computational complexity. Nevertheless, this 2D framework fails for capturing rich physics compromising the usual hypotheses considered when deriving standard plate and shell theories. To circumvent, or at least alleviate this issue, authors proposed in their former works an in-plane-out-of-plane separated representation able to capture rich 3D behaviors while keeping the computational complexity of 2D simulations. However, that procedure it was revealed to be too intrusive for being introduced into existing commercial softwares. Moreover, experience indicated that such enriched descriptions are only compulsory locally, in some regions or structure components. In the present paper we propose an enrichment procedure able to address 3D local behaviors, preserving the direct minimally-invasive coupling with existing plate and shell discretizations. The proposed strategy will be extended to inelastic behaviors and structural dynamics

Spurious-free interpolations for non-intrusive PGD-based parametric solutions: Application to composites forming processes

Non-intrusive approaches for the construction of computational vademecums face different challenges, especially when a parameter variation affects the physics of the problem considerably. In these situations, classical interpolation becomes inaccurate. Therefore, classical approaches for the construction of an offline computational vademecum, typically by using model reduction techniques, are no longer valid. Such problems are faced in different physical simulations, for example welding path problems, resin transfer molding, or sheet compression molding, among others. In such situations, the interpolation of precomputed solutions at prescribed parameter values (built using either intrusive or non intrusive techniques) generates spurious numerical artifacts. In this work, we propose an alternative interpolation and simulation strategy by using physically-based morphing of spaces. The morphing will transform the uncompatibe physical domains of the problem’s solution into a compatible one, where an interpolation free of artifacts can be performed. Later on, an inverse transformation can be used to push-back the solution. Different relevant examples are illustrated in this work to motivate the use of the proposed metho

Numerical Investigations of a High Frequency Pulsed Gaseous Fuel Jet Injection into a Supersonic Crossflow

The investigation of fuel delivery mechanisms is a critical design point in the development of supersonic combustion ramjet (scramjet) technology. Primary challenges include proper penetration of the jet in the supersonic cross-flow while keeping total pressure losses and wall drag to a minimum. To reduce drag and heat loads especially at high burner entry Mach numbers it is desirable to use a minimally intrusive means of fuel delivery. Pulsation of gaseous jets has been shown to increase penetration and mixing in subsonic flows. A limited number of experimental studies and even fewer numerical studies have suggested that when applied to supersonic crossflows, gaseous jets pulsed in the kilohertz range of frequencies improve jet penetration and mixing. To improve on the limited number of numerical studies of pulsed jets in supersonic crossflows (PJISF), 2D and 3D computational fluid dynamics (CFD) simulation models of non-excited (steady) and sinusoidally excited (pulsed) jets were constructed using ANSYS FLUENT 15.0. The 2D results included pulsation at 8, 16, 32 and 48 kHz. These simulation results showed that pulsation at 16 kHz provided the best penetration improvement in the jet near field and far field among all frequencies sampled. A 3D wall-modeled Large Eddy Simulation (WMLES) was constructed with the goals resolving large scale turbulent flow structure and time evolution of a jet pulsed in a supersonic crossflow, as well as to compare the effects of sinusoidal pulsation at 16 kHz with steady injection for the same flow conditions as the 2D case. A comparison of the jet trajectories between the steady and pulsed injection cases demonstrated that for sinusoidal pulsation of a jet at 16 kHz over the equivalent cycle averaged injection total pressure and momentum flux ratio, pressure, jet penetration is improved over the steady jet, up to 50% in the near field of the jet. Furthermore, improved mass concentration decay associated with jet-crossflow mixing and far field total pressure recovery has been demonstrated as a result of pulsation of the jet

Sub-minute in situ Fracture Test in a Lab CT-scanner

International audienceThe present study aims at demonstrating the feasibility of performing a fracture test in less than one minute in a lab CT scanner despite the severe time constraints of tomography acquisition. After introducing the basic concepts of Projection-based Digital Volume Correlation (P-DVC), the specific implementation of this methodology to a wedge splitting test on a refractory material is presented. The kinematics of the test is described over a mesh tailored to the sample geometry, and the elastic behavior of the sample is exploited through finite element computations to provide sensitivity fields of experimental boundary conditions to allow for their "measurements." Enhancing the simulation to account for crack advance with extended finite element analyses allows the sensitivity of the procedure to the crack position to be assessed. A confidence interval for the refractory toughness is finally obtained

Proper Generalized Decomposition for Parametric Study and Material Distribution Design of Multi-Directional Functionally Graded Plates Based on 3D Elasticity Solution

The use of mesh-based numerical methods for a 3D elasticity solution of thick plates involves high computational costs. This particularly limits parametric studies and material distribution design problems because they need a large number of independent simulations to evaluate the effects of material distribution and optimization. In this context, in the current work, the Proper Generalized Decomposition (PGD) technique is adopted to overcome this difficulty and solve the 3D elasticity problems in a high-dimensional parametric space. PGD is an a priori model order reduction technique that reduces the solution of 3D partial differential equations into a set of 1D ordinary differential equations, which can be solved easily. Moreover, PGD makes it possible to perform parametric solutions in a unified and efficient manner. In the present work, some examples of a parametric elasticity solution and material distribution design of multi-directional FGM composite thick plates are presented after some validation case studies to show the applicability of PGD in such problems

Application of PGD separation of space to create a reduced-order model of a lithium-ion cell structure

Lithium-ion cells can be considered a laminate of thin plies comprising the anode, separator, and cathode. Lithium-ion cells are vulnerable toward out-of-plane loading. When simulating such structures under out-of-plane mechanical loads, subordinate approaches such as shells or plates are sub-optimal because they are blind toward out-of-plane strains and stresses. On the other hand, the use of solid elements leads to limitations in terms of computational efficiency independent of the time integration method. In this paper, the bottlenecks of both (implicit and explicit) methods are discussed, and an alternative approach is shown. Proper generalized decomposition (PGD) is used for this purpose. This computational method makes it possible to divide the problem into the characteristic in-plane and out-of-plane behaviors. The separation of space achieved with this method is demonstrated on a static linearized problem of a lithium-ion cell structure. The results are compared with conventional solution approaches. Moreover, an in-plane/out-of-plane separated representation is also built using proper orthogonal decomposition (POD). This simply serves to compare the in-plane and out-of-plane behaviors estimated by the PGD and does not allow computational advantages relative to conventional techniques. Finally, the time savings and the resulting deviations are discussed