philippe.barabinot@siemens.com Siemens Industry Software SAS Digital Twins, which tend to intervene over the entire life cycle of products from early 92320 Châtillon, France design phase to predictive maintenance through optimization processes, are Full list of author information is available at the end of the article increasingly emerging as an essential component in the future of industries. To reduce the computational time reduced-order modeling (ROM) methods can be useful. However, the spread of ROM methods at an industrial level is currently hampered by the diﬃculty of introducing them into commercial ﬁnite element software, due to the strong intrusiveness of the associated algorithms, preventing from getting robust and reliable tools all integrated in a certiﬁed product. This work tries to circumvent this issue by introducing a weakly-invasive reformulation of the LATIN-PGD method which is TM intended to be directly embedded into Simcenter Samcef ﬁnite element software. The originality of this approach lies in the remarkably general way of doing, allowing PGD method to deal with not only a particular application but with all facilities already included in such softwares—any non-linearities, any element types, any boundary conditions...—and thus providing a new high-performance all-inclusive non-linear solver. Keywords: Digital Twin, Reduced-order model, Non-linear, LATIN-PGD, TM Weakly-invasive, Simcenter Samcef Introduction The digitization of the manufacturing industry is a topic of global interest known as “Industry 4.0”. Companies are prone to generate a digital representation of their entire value chain in order to be more eﬃcient with a more rational use of resources and energy. One talks about “Digital Twins” [1,2] which include two main components: ultra-high ﬁdelity simulations and data collected during the life of the structure [3]. Our approach is more particularly focused on the ﬁrst component—namely performing ultra-high ﬁdelity simulations very quickly—and hence we do not intend to deal with data in this paper. Despite the current strong interest in data, one of the requirements for the Digital Twin ediﬁce is still based on the simulation part. Reduced-order modeling (ROM) methods [4] are a key technology for the simulation component requiring a very fast and robust 3D simulation [5]. They have the advantage of drastically reducing CPU time without requiring any compromise on precision where high-ﬁdelity solutions can always be obtained. The richness of the physical phenom- © The Author(s) 2021. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 0123456789().,–: volV Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 2 of 15 ena remain purely intact, only the way we are representing them—as separate variables format—diﬀers and allows for interesting time savings. Additionally, ROM methods can lead to a better understanding of certain physical phenomena which were previously diﬃcult or even impossible to simulate due to the curse of dimensionality. To be used on a large scale and in a reliable way, ROM methods must be integrated as close as possible into the tools that engineers use on a daily basis, which require the joint use with industrial software. Hence, one requirement is to develop non-intrusive ROM methodologies [6]. Some recent papers are interested in this issue [7–12]. Most of the time, it results from a coupling between industrial software used as a black-box and in-house software that drives the ROM methods for parametrized problems on speciﬁc applications. Another point lies in the management of time-dependent non-linear problems with ROM methods. In this context, the LATIN-PGD method [13,14] is a worthwhile one, which has been the subject of many works highlighting its assets: ostensible reduction of computation time [15,16], construction of virtual charts [7], ability to address all types of non-linearities [17,18] among others. Hence, the LATIN-PGD method is attractive to Siemens Digital Industries Software, as a promising tool in the Digital Twins landscape. As a software company, Siemens Digital Industries Software wants to integrate this method into the TM TM ﬁnite element industrial software Simcenter Samcef part of the Simcenter 3D package. The idea is to get the easiest possible implementation, even if it results in slightly intrusive parts at some places when necessary. This is one of the beneﬁts of having access to the core routines of the software. Three main elements drive the approach: the robustness, the performance and the ease of use by external customers. This work proposes a strategy to use the Proper Generalized Decomposition (PGD) method [19] directly within an industrial software capable of handling any time-dependent non-linear problem in solid mechanics. We develop a weakly-invasive version of the LATIN-PGD method, which diﬀers from the one most commonly used in recent years based on the internal variables description of materials [13,15,20]. Instead, it can be seen as an extension of the ﬁrst functional version [21–24]. Numerous numerical treatments have been made to ensure a maximum generality of the presented approach designed to deal with any non-linearity (material, contact, large deformation,...). In addition to per- formance aspect, robustness remains an essential requirement to be able to deal with any industrial case. While until now the LATIN-PGD method has been used to solve spe- ciﬁc applications (see appendix in [25] for related topics), the presented approach shows the possibilities oﬀered by the method with an industrial software where no particular application is targeted. The novelty is the methodology put in place to use the PGD natively and as seamlessly as possible in an industrial software, thus providing an end-to-end process in one single soft- ware environment. The LATIN-PGD method treats the problem simultaneously in space and time. Regarding parameters of the Digital Twin, they are managed externally via a classical non-intrusive snapshots approach leveraging that the LATIN-PGD algorithm can be initialized from any previous computation [15,26]. In general, the (extra-)parameters cannot be known in advance for all possible particular applications as they vary from one to another. Hence, by robustness concerns, only a space–time PGD is managed by the TM Simcenter Samcef software itself because space and time are always present in non- Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 3 of 15 linear computations. Then, if required, a space–time-parameters PGD can be obtained a posteriori following classical non-intrusive technics [8,10]based on snapshots. TM Simcenter Samcef is organized as diﬀerent modules to be able to simulate diﬀer- ent physics: non-linear mechanical problems, high-speed rotating machines, articulated mechanisms or thermal simulations. As these modules are all part of the same family, it is possible to switch from one analysis to another, combine analyses or conduct co- simulations. In this work, we focus on the module dedicated to non-linear mechanical solver as an example. However, because of its generality, our approach could be used in other modules (thermal one for example) or other ﬁnite element software. The only requirement to set up the LATIN-PGD method is to have a non-linear resolution scheme already in place (Newton–Raphson), which notably solves local non-linearities. This is the strength of the proposed version to reuse without any modiﬁcation the full management of non-linear equations in the local stage. This paper focuses on the implementation—and the diﬃculties associated with—of the LATIN-PGD method in an industrial ﬁnite element software. Fundamentals as conver- gence and robustness studies will be presented in companion paper [27]. At the present time, what has been done ﬁts into the framework of moderate displacements, outside instability zones. In what follows, after a brief reminder of the type of mechanical problems that are con- sidered, the LATIN-PGD strategy and some details on its implementation in Simcenter TM Samcef software are provided. Finally, results on two industrial test-cases are given to highlight the assets. TM Weakly-invasive LATIN-PGD within Simcenter Samcef software A wide variety of problems in many ﬁelds of physics—mechanics, thermics, rotor TM dynamics—can be handled with the Simcenter Samcef software. In this study, we focus on the mechanical solver capable of solving any non-linear evolutionary problem including potentially materials non-linearities, contacts or large deformations. One of the goals is to combine ROM methods (by PGD in this work) with all possible non-linearities that already exist in the software. Types of mechanical problems involved Following a generic manner, any physical system subjected to a whole set of external sol- licitations (displacements, eﬀorts, pressures, temperature...) over a time interval I = [0,T] is considered under quasi-static assumption. Using a classical ﬁnite element formulation, the problem to be solved can be systematically reduced to the following semi-discretized problem: Find q(t) ∈ Q such as F (t) = F t, q(t) ∀t ∈ I (1) d i with F t, q(t) = T t, q(τ t) ∀t ∈ I (2) where Q deﬁnes the space of kinematically admissible ﬁelds subject to suﬃcient regularity, q groups all the displacements unknown, F gathers the external loadings, F represents d i the generalized force associated with internal cohesion eﬀorts and T symbolizes a charac- teristic operator describing all the material non-linearities. This operator T is built from Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 4 of 15 the processing of local elements within the considered industrial software. It will act as a black-box in our approach as explained later. Remark Under the notation T , internal variables can be hidden if some visco-plastic behavioral laws are considered. Let us note in particular that there is an irreversibility in time within the various non-linearities characterized by the operator T i.e. the expression at a given time t depends on the loading history according to a functional formalism. In case of contacts, additional contact forces F have to be taken into account, which requires the introduction of two additional local variables (v, λ)with v the trace of the displacements on the parts in contact and λ the Lagrange multipliers. Hence, the reference problem becomes: Find q(t), λ(t) such as (3) q(t) ∈ Q ∀t ∈ I (4) F (t) = F t, q(t) + F t, v(t), λ(t) ∀t ∈ I (5) i c F t, q(t) = T t, q(τ t) ∀t ∈ I (6) R t, v(t), λ(t) = 0 ∀t ∈ I (7) where in the ﬁnite element formulation, one can deﬁne an operator B such as v = Bq and the operator R represents the non-linear contact equations. This problem is classically solved incrementally in industrial software including Simcenter TM Samcef software. For quasi-static sollicitations, a time implicit scheme is commonly used which consists in solving for each time step a non-linear problem that can be treated by a Newton or a modiﬁed Newton–Raphson algorithm [28–30]. LATIN-PGD methodology The novelty of our approach is to natively set up a space–time PGD method directly TM in an industrial ﬁnite element software (Simcenter Samcef here) in order to best integrate the reduced-order modeling tools into a uniﬁed workﬂow that can be easily handled by customers. One of the requirements is to integrate a new non-linear solver in the software which has to be non-incremental. We choose the non-linear LATIN solver [13,31] designed to work well with the PGD method. This is a non-incremental iterative solver producing at each iteration a better approximation of the solution over the whole space and time. While standard industrial non-linear solvers generally used are incremental by nature, we propose in this work a gateway to eﬀortlessly implement a non-incremental solver directly within incremental one. Schematically, this non-incremental iterative algorithm can be represented as shown in Fig. 1 where two sub-spaces A and are deﬁned. Sub-space A groups the linear equa- d d tions of the problem—equilibrium and kinematic compatibility equations—while sub- space groups the local non-linear equations—behavioral laws, contact equations... The solution of the problem (1–2) then results from the fulﬁllment of this iterative scheme + − with two search directions ϒ and ϒ representing the parameters of the method. More speciﬁcally, Eq. (1) can be identiﬁed as A while Eq. (2)deﬁnes sub-space. The unknowns in this iterative scheme are the generalized displacements q and forces F which diﬀer from classical variables most recently used in the LATIN framework based on internal variables. This choice is guided by the ﬁnite element software. Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 5 of 15 Fig. 1 Convergent iterative scheme for LATIN-PGD algorithm with s = q , F n i,n Usually starting from a linear elastic initialization q (t), the algorithm leads us through a succession of local and global stages by means of the search directions to the terms of a convergent series q (t) . We suppose to be at a given state of the algorithm where k∈N q (t) is considered to be known. Local stage: s (t) −→ s (t) n n Given s (t), a new solution s (t) is sought in sub-space by solving the non-linear equa- n n tions. One has: Constitutive relations F t, q (t) = T t, q (τ t) ∀t ∈ I (8) i,n n n Ascent search direction ϒ q (t) = q (t) ∀t ∈ I (9) n n We use a vertical search direction because this is what is used naturally in Newton– Raphson’s algorithm. Indeed, at each iteration, one provides the displacement and one recovers the associated internal and external forces to be able to solve the equilibrium afterwards. Kinematics being imposed, this step is local in space and operates at each gauss point of each element. By making this choice for the ascent search direction, we can then reuse exactly the same routines as classically used. The only diﬀerence lies in the fact that these routines are no longer used once but several times in order to process all the time steps at each LATIN-PGD iteration. Global stage: s (t) −→ s (t) n n+1 Given s (t) and an operator H that deﬁnes ϒ as a parameter, a new solution s (t)is n n+1 sought in A sub-space by solving the equilibrium equations. One has: Kinematic conditions and equilibrium equation A q (t) ∈ Q ∀t ∈ I (10) n+1 F (t) = F (t) ∀t ∈ I (11) t, q d i,n+1 n+1 Descent search direction ϒ Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 6 of 15 H q (t) − q (t) = F t, q (t) −F t, q (t) ∀t ∈ I (12) i,n+1 i,n n+1 n n+1 n F (t) The choice of search direction operator H is the major parameter of the method. The optimal choice in terms of convergence rate would be to consider the tangent operator at each time steps, while a less costly choice could be to consider only the elastic Hooke operator. In similar fashion, when contacts are involved, one considers s = q , F , v , λ which n i,n n n leads to the following LATIN scheme: Local stage: s (t) −→ s (t) n n Constitutive relations F t, q (t) = T t, q (τ t) ∀t ∈ I (13) i,n n n R t, v(t), λ(t) = 0 ∀t ∈ I (14) Ascent search direction ϒ q (t) = q (t) ∀t ∈ I (15) n n Global stage: s (t) −→ s (t) n n+1 Kinematic conditions and equilibrium equation A q (t) ∈ Q ∀t ∈ I (16) n+1 F (t) = F t, q (t) + F t, v(t), λ(t) ∀t ∈ I (17) i,n+1 c,n+1 d n+1 Descent search direction ϒ H q (t) − q (t) = F t, q (t) − F t, q (t) ∀t ∈ I (18) i,n+1 i,n n+1 n n+1 n λ(t) − λ(t) = h v(t) − v(t) ∀t ∈ I (19) Both search directions H and h are the parameters of the method. One can interpret h as an interface rigidity [18]. The computation of the R operator comes directly from the hosting software as a black-box. No particular modiﬁcation is required for the local stage. In this form, the LATIN-PGD algorithm can be seen as the analogue of the Newton– Raphson algorithm to which convergent loop and time steps loop have been swapped. This is a crucial aspect in terms of ease of implementation as the method could be seamlessly TM integrated into Simcenter Samcef software without breaking anything and without having to redevelop speciﬁc tools at this level. Thus, it is henceforth possible to choose the non-linear solver that will be used for the resolution: either Newton–Raphson or LATIN methods. In the latter case, as only linear equations on the whole space and time are involved in the global stage, one can use PGD technics [4,32]. Focus on the PGD The PGD method [31,33] consists in approximating the correction calculated in the linear global stage by a ﬁnite sum of functions under separate variables: q (t) − q (t) ≈ λ (t) (20) n+1 n k k k=1 Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 7 of 15 This low-rank approximation is the reduced-order basis of the problem (1–2), which is obtained using a greedy algorithm. More speciﬁcally, a ﬁrst step consists in carrying out an update phase of the temporal modes already constituted, the spatial ones being kept orthogonal to each other according to the Hooke operator. Then, a second step consists in enhancing the reduced-order basis by adding some new modes if the update phase were proved to be not suﬃciently eﬀective. Thus, as the LATIN-PGD iterations progress, one comes to improve the solution until convergence. A criterion based on the energy norm averaged over the time interval enables to determine the state of convergence of the algorithm at each iteration and thus set a stop criterion when a certain threshold is reached. Some actions are currently underway to bring out these modes directly within the post- TM processing graphical interface provided by Simcenter 3D platform in order to visual- ized them as easily as possible. The interest is twofold insofar as this data format under separated variables also allows to signiﬁcantly compress the information on the outputs and therefore increases even more the execution speed of such post-processing interface. TM Integration of the LATIN-PGD method into Simcenter Samcef software From a software developer’s point of view, only a few parsimonious modiﬁcations have been necessary to introduced the LATIN-PGD method in pre-existing Simcenter TM Samcef software notably through the use of a dedicated LATIN-PGD fortran module gathering all the new functionalities in an unique place. Everything that makes up the TM signature and the richness of Simcenter Samcef software—the way of treating behav- ioral laws, the way of managing contacts and so on—is fully recycled. In a sense, one of the strength of our approach lies in the privileged place granted to this local stage oper- ating at the level of each element, each gauss point of each time steps, thus guaranteeing robustness and generality. A delicate point that has required particular attention concerns data storage since complete spatial-temporal approximations have to be stored at each LATIN-PGD iteration which is quite uncommon within incremental solvers. Neverthe- less, the judicious use of pre-existing memory spaces, as well as the use of a few dynamic allocations ensure that the methodology can be carried out smoothly, taking advantage in particular of the fact that the solution can be written in a separable variables format. The result is a controlled memory increase. Seen in a macroscopic way, all this boils down to an interface management between diﬀerent mathematical spaces grouping together the equations to solve, therefore to establish an eﬃcient communication between the various parts of the code. Another interesting thing is the notion of having a multi-ﬁdelity solver [34]: typically in design oﬃces, one often only needs to know an idea of the solution, a trend, especially in the upstream phases of the design. This allows a better understanding of the overall functioning of a new system before establishing afterwards which studies should then be carried out in greater depth by calculating high-ﬁdelity solutions, always at low cost. Management of parametrized aspects In an industrial environment, it is not uncommon to perform these non-linear space– time computations several times depending on some sets of parameters. The notion of parameter should be taken in a broad sense. By parameter, we mean for instance material Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 8 of 15 data, boundary conditions, number of cycles, geometry variations (morphing),... In what follows, these parameters are managed in a completely external way taking into account an interesting property of the LATIN-PGD method: the algorithm can be initialized with any previous kinematically admissible ﬁeld. Hence, we do not explicitly refer to the param- eters in the notations although our methodology is designed to deal with models involving parameters as illustrated in the numerical results. There is nothing new here, we use the classical techniques as stated in [15,26]. More precisely, if one notes μ the set of p param- TM eters μ = (μ , ..., μ ), then each LATIN-PGD computation with Simcenter Samcef 1 p (μ) software provides a space–time PGD approximation q for one value of the parameters such as: (μ) (μ) (μ) (μ) q (t) = ξ (t) (21) k k k=1 Next, if desired, one can then compress the information from these diﬀerent reduced- order basis to reconstitute a posteriori a space–time-parameters PGD [8,10]: m p q(t, μ) = ξ (t) ω (μ ) (22) k k k i k=1 i=1 Some industrial illustrations In the following, two industrial test-cases are presented. They both highlight, on the one hand, the possibilities of the LATIN-PGD method incorporating within Simcenter TM Samcef ﬁnite element software to build eﬃciently reduced-order basis and, on the other hand, the gain in computation time compared to classical approach. Combustion engine piston The ﬁrst test-case deals with a combustion engine piston submitted to pressure forces on the upper part. Due to the symmetry of the problem, both in terms of loadings and geometry, only a quarter of the piston is modeled and the boundary conditions are set accordingly. In addition, since the piston can only move vertically in its chamber, zero normal boundary conditions are ﬁxed on the outer cylindrical part. Finally, one leads to the resolution of an evolutionary mechanical problem involving 1 million of degrees of freedom, subjected to cyclic pressure loading in order to simulate the back and forth motion of the piston. The material is non-linear Norton law whose equivalent cumulative plastic strain rate p˙ follows the formula: J (σ) − k p˙ = (23) 1/m Kp where J (σ) represents the equivalent Von Mises stress, k enables to take into account an isotropic hardening and (K, n, m) are three parameters depending on the temperature T of the combustion chamber which is assumed to be uniform for a sake of simplicity. We propose in the following to carry out a parametric study in order to determine the whole space–time solution according to diﬀerent values of the temperature. In particular, we take out the curve showing the evolution of the cumulative plastic strain during cyclic Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 9 of 15 Fig. 2 Cumulative plastic strain output with LATIN-PGD method (left) and deviation from the reference (right) Fig. 3 Temporal evolution of the equivalent cumulative plastic strain according to temperature T for the most loaded element loading resulting. This curve results from 40 computations. As an example, experiments are achieved on 3 loading cycles but more cycles could have been set. Figure 2 gives the spatial distribution of the cumulative plastic strain at the end of the third cycle where the ﬁrst ﬁgure on the left shows the output with LATIN-PGD method and the second ﬁgure on the right displays the deviation from the reference calculated with classical Newton–Raphson scheme. Thus, with only 9 modes, the low-rank approximation of the LATIN-PGD method is able to perfectly capture the high-ﬁdelity solution as the error does not exceed half a percent on all elements. On Figs. 3 and 4 evolution of the cumulative plastic strain along the time steps are exhib- ited. All these calculations have been carried out to achieve a high ﬁdelity solution. While TM almost 22 h are required to achieve all simulations with standard Simcenter Samcef Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 10 of 15 Fig. 4 Equivalent cumulative plastic strain at the end of each cycle according to temperature T on the most loaded element of the structure software, only 7 h are necessary to carry out the entire parametric study with reduced- order model which leads to a speed-up of 3 times faster. Moreover, there is a great disparity in the calculation time: while it remains more or less constant at around 30 min per run with the classical approach, it can vary over a range from 4 to 18 min with new LATIN- PGD method depending on the initialization of the solution. Here, a rather naive path have been pursued consisting in starting from the closest neighbor in the parametric space which proves to be very eﬀective when the gradient evolution remains relatively low (typically the range from 515 to 525 C) and requires a little more time otherwise. In [26] one can ﬁnd strategies to follow for a more judicious selection of the path. However, we note that we always outperform the classical resolution except for the ﬁrst calculation which required 50 min—the time to build a ﬁrst reduced-order basis through PGD from an elastic-solution. Of course, one could probably gain even more by implementing a ﬁrst phase of estimation or extrapolation of the reduced-order basis as an initialization of the algorithm for a new calculation. The better the estimation, the shorter the computation time, but the fact remains that a reliable high-ﬁdelity solution will always be obtained by proceeding again a few iterations with the LATIN-PGD algorithm if necessary. Turbine blade-disk assembly The second test-case concerns a turbine blade-disk assembly involving material non- linearities—plasticity and creep—and contact non-linearities at the junction between the disk and the blade. The assembly is subjected to centrifugal forces as well as pressure ones on the surface of the blade. Additionally, a temperature gradient is applied according to the radial coordinate. Due to the rotating machine loads and geometry, cyclic symmetry conditions are also prescribed on the lateral parts. This geometry consists of 1 million of degrees of freedom, including 40,000 Lagrange multipliers, resulting from diﬀerent types of volumic elements—tetrahedrons or hexahedrons—and speciﬁc elements at the interfaces between the diﬀerent parts of the mesh to write kinematic conditions. Because Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 11 of 15 Fig. 5 Initial conﬁguration (gray) and deformed geometry (color) for the most loading time step of its relatively large size (about 200 mm) and the loadings, the blade unscrews as illustrated on Figure 5, where one can appreciate the deformed geometry in color with respect to the initial conﬁguration in gray. A maximum displacement of 3.86 mm is observed at the end of the structure. In what follows, a two parameters study is undertaken where the parameters are the Young modulus E and the thermal expansion coeﬃcient α of the blade. Values of these parameters are assumed to be insuﬃciently known, mainly related to the manufacturing process. Due to loading, a small variation in the value of these coeﬃcients can aﬀect the maximum Von Mises stress in the contact area where creep occurs over long periods of time, which can damage the structure faster than expected. Although other parameters may be of inﬂuence, we focus on these two parameters for a sake of simplicity. Figure 6 −6 presents the stress map obtained for a given point (E = 220 MPa and α = 14 · 10 −1 K ) in the parametric space. One gives the results of the LATIN-PGD method included TM in Simcenter Samcef software and the remaining error with respect to the reference. The latter is taken as the classical solution from the commercial software (using Newton– Raphson method). One can note a very good concordance of the solutions as the relative error remains less than one percent in the area of interest for 10 PGD modes. Figure 7 shows the response surface of the maximum Von Mises stress in the contact area according to the two parameters. This chart highlights a non-linear variation accord- ing to the Young modulus E and an almost linear one on the second parameter α.In order to achieve this, 70 computations have been carried out by scanning the parametric Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 12 of 15 Fig. 6 Von Mises stress map [MPa] with the LATIN-PGD method (left) and deviation from the reference in percent [%] (right) 1.42 1.41 1.4 −5 1.39 220 ·10 1.38 200 −1 α (K ) E (MPa) Fig. 7 Maximum Von Mises stress in the contact area according to both E and α parameters resulting from 70 parametric points uniformly pulled across the parametric space space sequentially in a uniform manner. Thanks to both ROM methods and LATIN-PGD algorithm allowing to restart the computation from a close previous converged point in the parametric space, only 3 h have been necessary to build the virtual chart while more TM than 21 h would have been necessary using classical Simcenter Samcef software based solely on Newton–Raphson algorithm. This leads to 7 times faster computation. This test-case highlights one of the advantages having the LATIN-PGD method incor- porated natively within an industrial software insofar as one can then exploit all the capabilities of the hosting software at the local stage. As an example, 3 diﬀerent behavior laws (depending on the temperature) have been considered here, as well as diﬀerent types of elements, in addition to the contact behavior. Maximum Von Mises stress (MPa) Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 13 of 15 Conclusions In this work, the LATIN-PGD method have been directly implemented into Simcenter TM Samcef ﬁnite element software, which allows to build and then enrich reduced-order model as desired at the industrial level. Here, a very general way of looking is adopted with PGD by making the best use of all the local behavior pre-existing in industrial softwares. The PGD method becomes in fact integrated as close as possible to the solver, enabling to combine performance, robustness and a wide range of applications that can be carried out. Moreover, a ﬁrst time saving on parametric studies has been demonstrated but a signiﬁcantly additional gain could still be made by exploiting hyper-reduction techniques [35–37] and in particular the Reference Point Method [14,38] well suited for LATIN-PGD method. For each quantity of interest, the latter consists in storing the information in an even more condensed manner at a few cleverly chosen points, an analytical formula being subsequently used to reconstruct the ﬁeld over the whole space and time. These methods are currently being investigated and should enable to go further in the reduction of com- putation time insofar as up to two thirds of the computation time can be spent in the local stage reﬂecting the fact that the current algorithm complexity is not independent of the number of degrees of freedom. Furthermore, one also keep in mind that the gain is directly correlated to the complexity of the problem: thus, considering a dozen of parameters and several hundreds or thousands of loading cycles, the performance would be even more better. Examples on material and contact non-linearities have been presented above, but this very general methodology can deal with any other type of non-linearity already present TM in the industrial Simcenter Samcef software. This is the strength of this approach to be able to treat any type of ﬁnite element, any geometry, any boundary condition, any loading or even any non-linearities all using reduced-order modeling methods in order to circumvent prohibitive computation times. However, at the present time, computations are limited to moderate displacements, outside instability zones. One of our next task will be to overcome this limitation. Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 14 of 15 Acknowledgements This work is funded by Siemens Industry Software SAS and Association Nationale de la Recherche et de la Technologie (ANRT). Authors also acknowledge the support of the CNRS research group 3587 AMORE. Declarations Competing interests The authors declare that they have no competing interests. Author details Siemens Industry Software SAS 92320 Châtillon, France, Université Paris-Saclay, ENS Paris-Saclay, CNRS, LMT - Laboratoire de Mécanique et Technologie, 91190 Gif-sur-Yvette, France. Received: 2 February 2021 Accepted: 24 August 2021 References 1. Chinesta F, Cueto E, Abisset-Chavanne E, Duval JL, Khaldi FE. Virtual, digital and hybrid twins: a new paradigm in data-based engineering and engineered data. Arch Comput Methods Eng. 2020;27(1):105–34. https://doi.org/10. 1007/s11831-018-9301-4. 2. Hartmann D, van der Auweraer H. Digital twins [cs]. 2020. arXiv:2001.09747. 3. Glaessgen E, Stargel D. The digital twin paradigm for future NASA and US Air Force vehicles. In: 53rd AIAA/ASME/ASCE/AHS/ASC structures, structural dynamics and materials conference. American Institute of Aero- nautics and Astronautics. https://doi.org/10.2514/6.2012-1818. 4. Benner P, Grivet-Talocia S, Quarteroni A, Rozza G, Schilders W, Silveira LM, editors. Model order reduction: snapshot- based methods and algorithms, vol. 2. Berlin: De Gruyter; 2021. 5. Chinesta F, Huerta A, Rozza G, Willcox K. Model reduction methods. In: Encyclopedia of computational mechanics, vol. 3. 2nd ed. New York: American Cancer Society; 2017. p. 1–36. https://doi.org/10.1002/9781119176817.ecm2110. 6. Giraldi L, Litvinenko A, Liu D, Matthies HG, Nouy A. To be or not to be intrusive? The solution of parametric and stochastic equations-the “plain vanilla” Galerkin case. SIAM J Sci Comput. 2014;36(6):2720–44. https://doi.org/10. 1137/130942802. 7. Courard A, Néron D, Ladevèze P, Ballere L. Integration of PGD-virtual charts into an engineering design process. Comput Mech. 2016;57(4):637–51. https://doi.org/10.1007/s00466-015-1246-y. 8. Ibáñez R, Abisset-Chavanne E, Ammar A, González D, Cueto E, Huerta A, Duval JL, Chinesta F. A multidimensional data-driven sparse identiﬁcation technique: the sparse proper generalized decomposition. Complexity. 2018. https:// doi.org/10.1155/2018/5608286. 9. Zou X, Conti M, Díez P, Auricchio F. A nonintrusive proper generalized decomposition scheme with application in biomechanics. Int J Numer Methods Eng. 2018;113(2):230–51. https://doi.org/10.1002/nme.5610. 10. Borzacchiello D, Aguado JV, Chinesta F. Non-intrusive sparse subspace learning for parametrized problems. Arch Comput Methods Eng. 2019;26(2):303–26. https://doi.org/10.1007/s11831-017-9241-4. 11. Casenave F, Akkari N, Bordeu F, Rey C, Ryckelynck D. A nonintrusive distributed reduced-order modeling frame- work for nonlinear structural mechanics—application to elastoviscoplastic computations. Int J Numer Methods Eng. 2020;121(1):32–53. https://doi.org/10.1002/nme.6187. 12. Tsiolakis V, Giacomini M, Sevilla R, Othmer C, Huerta A. Nonintrusive proper generalised decomposition for parame- terised incompressible ﬂow problems in OpenFOAM. Comput Phys Commun. 2020;249:107013. https://doi.org/10. 1016/j.cpc.2019.107013. 13. Ladevèze P. Nonlinear computational structural mechanics: new approaches and non-incremental methods of calculation. New York: Springer; 1999. 14. Ladevèze P. On reduced models in nonlinear solid mechanics. Eur J Mech A/Solids. 2016;60:227–37. https://doi.org/ 10.1016/j.euromechsol.2016.08.005. 15. Relun N, Néron D, Boucard PA. A model reduction technique based on the PGD for elastic-viscoplastic computational analysis. Comput Mech. 2013;51(1):83–92. https://doi.org/10.1007/s00466-012-0706-x. 16. Néron D, Boucard P-A, Relun N. Time-space PGD for the rapid solution of 3D nonlinear parametrized problems in the many-query context. Int J Numer Methods Eng. 2015;103(4):275–92. https://doi.org/10.1002/nme.4893. 17. Abdali A, Benkrid K, Bussy P. Simulation of sheet cutting by the large time increment method. J Mater Process Technol. 1996;60(1):255–60. https://doi.org/10.1016/0924-0136(96)02339-4. 18. Oumaziz P, Gosselet P, Boucard P-A, Guinard S. A parallel non-invasive mixed domain decomposition— implementation and applications to mechanical assemblies. Finite Elem Anal Des. 2019;156:24–33. https://doi.org/ 10.1016/j.ﬁnel.2019.01.004. 19. Chinesta F, Ladevèze P. Proper generalized decomposition. In: Model order reduction: snapshot-based methods and algorithms, vol. 2. Berlin: De Gruyter; 2020. p. 97–138. https://doi.org/10.1515/9783110671490-003. 20. Bhattacharyya M, Fau A, Nackenhorst U, Néron D, Ladevèze P. A LATIN-based model reduction approach for the simulation of cycling damage. Comput Mech. 2018;62(4):725–43. https://doi.org/10.1007/s00466-017-1523-z. 21. Ladevèze P. Sur une famille d’algorithmes en mécanique des structures. Comptes-rendus des séances de l’Académie des sciences. Série 2, Mécanique-physique, chimie, sciences de l’univers, sciences de la terre. 1985;300(2):41–4. 22. Boisse P, Bussy P, Ladevèze P. A new approach in non-linear mechanics: the large time increment method. Int J Numer Methods Eng. 1990;29(3):647–63. https://doi.org/10.1002/nme.1620290312. Barabinot et al. Adv. Model. and Simul. in Eng. Sci. (2021) 8:22 Page 15 of 15 23. Bussy P, Rougée P, Vauchez P. The large time increment method for numerical simulation of metal forming processes. In: NUMETA. New York: Elsevier Science Ltd; 1990. p. 102–9. 24. Cognard J-Y, Ladevèze P. The large time increment method applied to cyclic loadings. In: Zyczkowski M, editor. Creep in structures. Berlin: Springer; 1991. p. 555–62. https://doi.org/10.1007/978-3-642-84455-362. 25. Scanﬀ R, Nachar S, Boucard P-A, Néron D. A study on the LATIN-pgd method: analysis of some variants in the light of the latest developments. Arch Comput Methods Eng. 2020. https://doi.org/10.1007/s11831-020-09514-1. 26. Heyberger C, Boucard P-A, Néron D. A rational strategy for the resolution of parametrized problems in the PGD framework. Comput Methods Appl Mech Eng. 2013;259:40–9. https://doi.org/10.1016/j.cma.2013.03.002. 27. Scanﬀ R, Néron D, Ladevèze P, Barabinot P, Cugnon F, Delsemme J-P. Weakly-invasive LATIN-PGD for solving time- dependent non-linear parametrized problems in solid mechanics; 2021. (In preparation). 28. Zienkiewicz OC, Taylor RL. The ﬁnite element method for solid and structural mechanics. New York: Elsevier; 2005. 29. de Borst R, Crisﬁeld MA, Remmers JJC, Verhoosel CV. Nonlinear ﬁnite element analysis of solids and structures. Chichester: Wiley; 2012. 30. Bonnet M, Frangi A, Rey C. The ﬁnite element method in solid mechanics. New York: McGraw Hill Education; 2014. 31. Ladevèze P. PGD in linear and nonlinear computational solid mechanics. In: Chinesta F, Ladevèze P, editors. Separated representations and PGD-based model reduction: fundamentals and applications. Vienna: Springer; 2014. p. 91–152. https://doi.org/10.1007/978-3-7091-1794-1_3. 32. Chinesta F, Ladevèze P, editors. Separated representations and PGD-based model reduction: fundamentals and applications. Courses and Lectures/International Centre for Mechanical Sciences, vol. 554. Wien: Springer; 2014. 33. Chinesta F, Keunings R, Leygue A. The proper generalized decomposition for advanced numerical simulations. SpringerBriefs in applied sciences and technology. Cham: Springer Science & Business Media; 2014. https://doi.org/ 10.1007/978-3-319-02865-1. 34. Nachar S, Boucard P-A, Néron D, Rey C. Multi-ﬁdelity Bayesian optimization using model-order reduction for viscoplas- tic structures. Finite Elem Anal Des. 2020;176:103400. https://doi.org/10.1016/j.ﬁnel.2020.103400. 35. Barrault M, Maday Y, Nguyen NC, Patera AT. An ‘empirical interpolation’ method: application to eﬃcient reduced-basis discretization of partial diﬀerential equations. Comptes Rendus Mathematique. 2004;339(9):667–72. https://doi.org/ 10.1016/j.crma.2004.08.006. 36. Ryckelynck D. A priori hyperreduction method: an adaptive approach. J Comput Phys. 2005;202(1):346–66. https:// doi.org/10.1016/j.jcp.2004.07.015. 37. Maday Y, Mula O. A generalized empirical interpolation method: application of reduced basis techniques to data assimilation. In: Brezzi F, Colli Franzone P, Gianazza U, Gilardi G, editors. Analysis and numerics of partial diﬀerential equations. Springer INdAM Series Springer INdAM Series Springer INdAM Series. Milano: Springer Milan; 2013. p. 221–35. https://doi.org/10.1007/978-88-470-2592-9_13. 38. Capaldo M, Guidault P-A, Néron D, Ladevèze P. The reference point method, a “hyperreduction” technique: application to PGD-based nonlinear model reduction. Comput Methods Appl Mech Eng. 2017;322:483–514. https://doi.org/10. 1016/j.cma.2017.04.033. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional aﬃliations.
"Advanced Modeling and Simulation in Engineering Sciences" – Springer Journals
Published: Oct 4, 2021
Keywords: Digital Twin; Reduced-order model; Non-linear; LATIN-PGD; Weakly-invasive; Simcenter SamcefTM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\hbox {TM}}$$\end{document}
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.