FE ANN – An eﬃcient data-driven multiscale approach based on physics-constrained neural networks and automated data mining

Herein, we present a new data-driven multiscale framework called FE ANN which is based on two main keystones: the usage of physics-constrained artiﬁcial neural networks (ANNs) as macroscopic surrogate models and an autonomous data mining process. Our approach allows the eﬃcient simulation of materials with complex underlying microstructures which reveal an overall anisotropic and nonlinear behavior on the macroscale. Thereby, we restrict ourselves to ﬁnite strain hyperelasticity problems for now. By using a set of problem speciﬁc invariants as the input of the ANN and the Helmholtz free energy density as the output, several physical principles, e. g., objectivity, material symmetry, compatibility with the balance of angular momentum and thermodynamic consistency are fulﬁlled a priori. The necessary data for the training of the ANN-based surrogate model, i. e., macroscopic deformations and corresponding stresses, are collected via computational homogenization of representative volume elements (RVEs). Thereby, the core feature of the approach is given by a completely autonomous mining of the required data set within an overall loop. In each iteration of the loop, new data are generated by gathering the macroscopic deformation states from the macroscopic ﬁnite element (FE) simulation and a subsequently sorting by using the anisotropy class of the considered material. Finally, all unknown deformations are prescribed in the RVE simulation to get the corresponding stresses and thus to extend the data set. The proposed framework consequently allows to reduce the number of time-consuming microscale simulations to a minimum. It is exemplarily applied to several descriptive examples, where a ﬁber reinforced composite with a highly nonlinear Ogden-type behavior of the individual components is considered. Thereby, a rather high accuracy could be proved by a validation of the approach.


Introduction
Materials with an underlying meso-or microstructure, e. g., composites, solid foams, dual-phase steels or 3D-printed structures, enable the targeted design of engineering components with respect to their application.However, due to effects such as anisotropy, nonlinear or multiphysics phenomena, the experimental characterization of the effective constitutive behavior of such materials can be very complex.

Multiscale schemes
To avoid such a characterization of the material's effective behavior, computational multiscale schemes can be used.These schemes allow the simulation of engineering components made of materials with underlying microstructure solely based on information about the microstructural arrangement and the properties of the individual components, e. g., matrix and inclusions.Basically, two different types of multiscale schemes exist: the coupled multiscale scheme which is also known as FE 2 approach (finite element square) [1,2,3,4,5,6,7,8] and the decoupled or sequential multiscale scheme [9,10,11,12,13,14].
The FE 2 method allows to completely couple the microscopic and macroscopic scales without the need for the explicit formulation of an effective constitutive model.It is thus universally applicable to arbitrary geometries if the necessary microscopic information are available.However, the decisive disadvantage is the very high computational effort, which results from the solution of the microscopic boundary value problem (BVP) for the homogenization at each integration point of the macroscopic FE mesh.
Within a decoupled multiscale scheme, the material's effective behavior is initially determined from homogenization of representative volume elements (RVEs) and then a suitable constitutive model, the so-called macro or surrogate model, is calibrated by these data.With this model, macroscopic BVPs can now be solved, whereby the influence of the microstructure is implicitly captured by the macro model.Thus, in contrast to the FE 2 method, the explicit solution of the microscopic BVP at each integration point is omitted.The central disadvantage, however, is that the formulation of such a model can be extremely complicated.

Data-based methods in solid mechanics
To circumvent the time consuming task of formulating and calibrating a surrogate model within a decoupled multi-Preprint submitted to Springer arXiv scale scheme, data-based or data-driven techniques are very promising and have the potential to improve or replace traditional constitutive models.These techniques have become increasingly popular in the computational mechanics community during the last years [15,16].In the following, a brief overview of the most common methods and their application to multiscale schemes is given.

Overview on data-based constitutive modeling
A relatively new strategy to substitute classical constitutive equations is the data-driven mechanics approach, initially proposed by Kirchdoerfer and Ortiz [17] and extended to, e. g., noisy data sets [18], finite strains [19], inelasticity [20] or fracture mechanics [21] in the meantime.This method completely avoids to use constitutive equations.Instead, sets of stress-strain tuples which characterize the material's behavior are used.A data-driven solver hence seeks to minimize the distance between the searched solution and the material data set within a proper energy norm, while compatibility and equilibrium have to be satisfied simultaneously.
The construction of so called constitutive manifolds from collected data is an alternative approach which is described for elasticity and inelasticity by Ibañez et al. [22].With this technique, it is possible to strictly fulfill the 2nd law of thermodynamics, i. e., the thermodynamic consistency, by using the GENERIC paradigm (General Equation for Non-Equilibrium Reversible-Irreversible Coupling) during the construction of constitutive manifolds [23].
Besides the previously mentioned techniques, there exist numerous data-based methods originating from the field of machine learning (ML).Probably, the most common technique is the application of artificial neural networks (ANNs), which have already been proposed in the early 90s by the pioneering work of Ghabussi et al. [24].In the last decades, ANNs have been intensively used for mechanical material modeling and simulations by means of the finite element method (FEM), e. g., in [25,26,27,28,29] among others.
However, in general, a large amount of data is required to train ANNs to serve as robust and accurate surrogate models for systems with complex underlying physics, e. g., constitutive models.In this context, a comparatively new branch of ML techniques related to constitutive modeling are approaches classified as physics-informed, physics-constrained [30,31], mechanics-informed [32] or hybrid models [33,34].In these methods, essential physical principles and information are inserted into the ML-based model or parts of classical models are replaced with data-based methods, which leads to an improvement of the extrapolation capability and enables training with sparse data.In the case of ANNs, this can be achieved via the network architecture or by adapted training In the following, the term physics-informed does not directly refer to the PINN-approach (physics-informed neural network) according to Raissi et al. [35].In a PINN, the searched solution field approximated by the ANN, e. g., displacement  (, ), is inserted into the governing partial differential equation (PDE) at collocation points.This expression is then added to the loss, so that the fulfillment of the PDE is enforced.In the context of constitutive modeling, however, the idea to enrich the ML approach with physical knowledge is applied in a similar way.algorithms.By choosing problem-specific invariant sets as the input variables, material symmetries are automatically satisfied for hyperelastic models [36,29,37].Furthermore, the thermodynamic consistency can be fulfilled by choosing the free energy as the output quantity.In order to train the ANN with respect to the stresses, gradients of the output with respect to the input are inserted into the loss [38,36].This technique is also named as Sobolev training in [39].Furthermore, physical knowledge can be inserted via constraint training processes [40].For the consideration of dissipative behavior, Masi et al. [41] proposed an adapted ANN architecture consisting of two feed forward neural networks (FNNs).Thereby, the first network is used to predict the evolution of internal variables and the second for approximating the free energy.Finally, the combination of classical models with data-based techniques is a further method to insert physical knowledge.Thereby, only parts of a model are replaced by data-based techniques, e. g., in plasticity models [33,34,42].

Data-based multiscale modeling and simulation
In the context of multiscale problems, the mentioned databased methods can be used as surrogate models which replace the computationally expensive simulation of RVEs [43].The high flexibility of trained networks with simultaneously excellent prediction qualities is thereby proven in numerous works: In [38] ANNs are used to describe the homogenized response of cubic lattice metamaterials exhibiting large deformations and buckling phenomena.Thereby, several types of hyperelastic ML-based models which fulfill basic physical requirements are used and compared to each other.Based on this, an extension to polyconvex hyperelastic models is given in [37].In [44], FFNs and RNNs (recurrent neural networks) are used to replicate the homogenization of RVEs revealing inelastic behavior of the individual components.Thereby, also unknown paths can be predicted with the trained networks.Similarly, ANNs that replace the inelastic constitutive responses of composite materials sampled by RVE simulations are shown in [45].Thereby, in addition, a deep reinforcement learning combinatorics game is used to automatically find an optimal set of network hyper-parameters from a decision tree.A data-driven multiscale framework called deep material network is shown in [46,47,48].Thereby, the homogenized RVE response is reproduced by a network including a collection of connected mechanistic building blocks with analytical homogenization solutions.With that, a complex effective response could be described without the loss of essential physics.An extension of the deep material network approach to fully coupled thermomechanical multiscale simulations of composite materials is given in [49].
Furthermore, ANNs can be used to link microstructural characteristics and effective properties.For example, in [50], an ML framework for predicting macroscopic yield as a function of crystallographic texture is described.An ML-based multiscale calibration of constitutive models representing the effective response of rate dependent composite materials is applied to brain white matter in [51].Based on a library of calibrated parameters corresponding to a set of microstructural characteristics, an ML model which predicts the constitutive model parameters directly from a new microstructure is trained.
Similarly, an ANN is trained to predict the elastic properties of short fiber reinforced plastics in [52].
Finally, the application of data-based techniques as surrogate models in decoupled multiscale schemes, which enable the simulation of macroscopic engineering components, is promising.In the pioneering work [53], a decoupled multiscale scheme using an ANN-based surrogate model has been shown for elastic composites with cubic microstructures.Another multiscale methodology is presented in [54].Therein, a hybrid macroscopic surrogate model is used, i. e., a traditional constitutive model is combined with a data-driven correction.
An ML-based multiscale framework for the simulation of the elastic response of metals having a polycrystalline microstructure is presented in [39].Therein, the database is generated by using a 3D FFT (fast Fourier transform) solver.
Based on this, [55,42] show an extension to elastoplasticity, whereby ANNs are used for the description of the yield surface and the stress within a hybrid modeling approach.Similarly, the simulation of the elastic-plastic deformation behavior of open-cell foam structures is shown in [33,34].Thereby, a hybrid model including two ANNs is used as the macroscopic surrogate model.Therein, the first network serves for the description of the macroscale yield function and the second one for the prediction of the floating direction.In [49], two scale simulations of thermo-mechanical problems are solved by using deep material networks.Applying the modeling strategy initially proposed in [41], a multiscale scheme for the inelastic behavior of lattice material structures is described in [56].An application of ANNs as surrogate models describing the anisotropic electrical response of graphene/polymer nanocomposites is shown in [57].Furthermore, in [58], the application of recurrent neural networks (RNNs) as surrogate models within multiscale simulations of elastoplastic problems is shown.Furthermore, the RNN-based approach is compared to full FE 2 simulations.An application of RNN surrogate models to viscoplasticity is described in [59].A combination of fully coupled FE 2 simulations with an adaptive switching to ANN-based surrogate models for the complex simulation of RVEs is shown in [60].
In addition, a multiscale framework based on the datadriven mechanics approach [17] is presented in [61] for the application to sand.Thereby, the necessary data set is extracted from lower scale simulations.In [62], a multi-level method is used within a data-driven multiscale scheme which is applied for the simulation of solid foam materials.

Content
Within this contribution, an efficient data-driven multiscale approach called FE ANN , which makes use of physics-constrained ANNs as a macroscopic surrogate model, is presented.The approach allows to consider materials with complex microstructures leading to an overall anisotropic behavior, whereby a restriction to hyperelastic solids is made for now.The ANNbased surrogate model automatically fulfills several physical principles, e. g., objectivity, material symmetry, compatibility with the balance of angular momentum and the thermodynamic consistency.This is done by using a set of problem specific invariants   as the input of the network and the Helmholtz free energy density  as the output, cf.Linka et al. [36], Klein et al. [37] or Linden et al. [63], among others.The data which are used to train the ANN are collected via computational homogenization of an RVE representing the material's microstructure.Thereby, in contrast to most of the datadriven multiscale approaches from the literature, it is not necessary to explore all required macroscopic states of deformation which occur in the macroscopic simulation in advance.Instead, the required data, i. e., effective deformations and corresponding stresses, are determined by homogenization in a fully autonomous way within the framework by collecting the relevant deformations from the macroscopic FE simulation.This procedure is similar to the approach presented by Korzeniowski and Weinberg [62] which is based on the data-driven mechanics approach [17].Here, in addition, the macroscopic deformations are mapped into an invariant space associated to the symmetry group of the considered material.In this space, the selection of relevant states is done so that the number of time-consuming microscale simulations can be reduced to a minimum.Moreover, it is not necessary to perform a rough scan of the relevant deformation area in advance here.The proposed framework is exemplarily applied to several descriptive examples, where a fiber reinforced composite with a highly nonlinear behavior of the individual components is considered.
The organization of the paper is as follows: In Sect.2, the basic equations of the finite strain continuum solid mechanics theory as well as basic principles of hyperelastic models are given.After this, the proposed data-driven multiscale framework is described in Sect.3.This approach is exemplarily applied within several numerical examples in Sect. 4. After a discussion of the results, the paper is closed by concluding remarks and an outlook to necessary future work in Sect. 5.

Continuum solid mechanics
In this section, the basic kinematic and stress quantities as well as general relations of anisotropic hyperelastic constitutive models are summarized shortly.The reader is referred to the textbooks of Haupt [64], Holzapfel [65] or Ogden [66] for a detailed overview.Furthermore, a Hill-type homogenization framework is introduced.For a clear mathematical notation, the space of tensors except for a tensor of rank zero, is used in the following.In Eq. (1), R 3 , N and ⊗ denote the Euclidean vector space, the set of natural numbers and the dyadic product, respectively.Tensors of rank one and two are given by boldface italic symbols in the following, i. e.,  ∈ L 1 or ,  ∈ L 2 .Furthermore, a single or double contraction of two tensors is given by  •  =       ⊗   or  :  =     , respectively.Thereby,   ∈ L 1 denotes a Cartesian basis vector and the Einstein summation convention is used.

Kinematics and stress measures 2.1.1 Kinematics
In the following, a material body C which occupies the reference configuration B 0 ⊂ R 3 at time  0 ∈ R ≥0 and the current configuration B ⊂ R 3 at time  ∈ T := { ∈ R + |  ≥  0 } is considered.The displacement vector  ∈ L 1 of a material point  ∈ C capturing the positions  ∈ B 0 at  0 and  ∈ B at  is given by ( , ) := ( , ) − .Therein,  : B 0 × T → B , ( , ) ↦ →  =: ( , ) denotes a bĳective motion function which is postulated to be continuous in space and time.As further kinematic quantities, the deformation gradient  ∈ L 2 and the Jacobian determinant  ∈ R + are defined by the relations  := (∇  ) T and  := det  > 0 . (2) In the equation above, ∇  is the nabla-operator with respect to reference configuration B 0 .
A deformation measure which is free of rigid body motions is given by the positive definite right Cauchy-Green deformation tensor  :=  T •  ∈  with the space of symmetric second order tensors  :=  ∈ L 2 |  =  T .With that and by using Sylvesters formula, the spectral decomposition of  follows to where 1 ∈ L 2 ,   ∈ R ≥0 and   ∈  denote the identity tensor, the principal stretches and the projection tensors, respectively.The introduced symbol   ∈ {1, 2, 3} indicates the algebraic multiplicity of   .Furthermore,   is defined as the Kronecker delta.

Stress measures
Within the framework of nonlinear continuum solid mechanics, various stress measures can be defined.The symmetric Cauchy stress  ∈ , also known as true stress, is defined with respect to the current configuration B. Furthermore, the 1st and 2nd Piola-Kirchhoff stress tensors  ∈ L 2 and  ∈  follow from the pull back operations Accordingly,  is related to both, B and B 0 , whereas  is completely defined with respect to B 0 .

General properties
The constitutive behavior of the considered solids is restricted to hyperelasticity within this work.Accordingly, a hyperelastic potential which is equal to the Helmholtz free energy density function  : L 2 → R + ,  ↦ → () exists.By applying the procedure of Coleman and Noll [67], the relations ( then follow from the evaluation of the Clausius-Planck inequality [65].With that, the thermodynamic consistency of any hyperelastic model is fulfilled a priori. In addition, there are some further requirements on , which ensure a physically reasonable constitutive behavior [65,66].The normalization condition requires that (1) = 0, i. e., that the free energy vanishes in the undeformed state.Furthermore, the free energy should increase in any case if deformation appears: () > 0 ∀ ≠ 1.The former two conditions imply that  has a global minimum at  = 1.Consequently, the undeformed state is stress-free, i. e.,  (1) = 0 holds.Additionally, the growth condition requires that an infinite amount of energy is needed to infinitely expand the volume or compress it to zero: Finally, the principle of material objectivity states that the free energy is invariant with respect to a superimposed rigid body motion.This statement is expressed by the relation () = ( • ) which holds for all special orthogonal tensors  ∈ Accordingly, the principle of material objectivity is automatically fulfilled if the tensor  is used as the argument of  instead of , i. e.,  :  → R + ,  ↦ → (), which is done in the following.
Finally, polyconvexity of the energy , i. e., convexity with respect to , Cof  :=   −T and det , is a further condition.This condition implies material stability, i. e., Legendre-Hadamard-Ellipticity, see Ebbing [68] for more details.However, polyconvexity is a quite strong requirement on the free energy.

Material symmetry
Besides the previously mentioned requirements, the constitutive equations should also reflect the symmetry of the described material which is expressed by the principle of material symmetry [64,68].According to that, the free energy have to be invariant with respect to all orthogonal transformations belonging to the symmetry group  of the underlying material: In order to describe anisotropic constitutive behavior, the concept of structural tensors can be used [64,68].Depending on the considered material, these tensors are of order two or even higher.They reflect the material's anisotropy and are thus invariant with respect to the symmetry transformations, e. g., The notation  * M  means                  , where the Einstein summation convention is used.If the structural tensors are appended to the list of arguments of , the energy is an isotropic tensor function even if the material is anisotropic which means that holds for all  ∈ .To abbreviate the notation, the sets M 2 := { 1 ,  2 , . . .,   2 } and M 4 := {M 1 , M 2 , . . ., M  4 } have been used in Eq. ( 7).
Finally, a set I := ( 1 ,  2 , . . .,   ) ∈ R ×1 consisting of  ∈ N irreducible scalar valued invariants   ∈ R could be derived by using the Cayley-Hamilton theorem.Consequently, the free energy is expressed by the isotropic tensor function  = ( 1 ,  2 , . . .,   ) which is invariant with respect to all transformations  ∈ .By using the latter definition and applying the chain rule, the 2nd Piola-Kirchhoff stress  follows to where   ∈  denotes tensor generators corresponding to the invariants   [29].

Special anisotropy classes
Within this work, two anisotropy classes are considered, isotropic as well as transversely isotropic materials.Corresponding sets of irreducible invariants are given in the following.
A set of irreducible invariants describing an isotropic hyperelastic solid is given by I The used principal invariants are given by whereby the cofactor of  is defined by Cof  :=  2  −1 .Note that I • is also expressible by the principal stretches   which follows from Eq. ( 3): According to [64,65,68], the structural tensor describing transverse isotropy is given by the second order tensor  =  ⊗ , whereby  ∈ L 1 with | | ≡ 1 is the fiber direction in the undeformed state.A set of irreducible invariants is thus I := ( 1 ,  2 ,  3 ,  4 ,  5 ) ∈ R 5×1 .Therein, the latter two invariants, which capture the material's anisotropy, are given by the expressions  4 :=  :  and  5 :=  :  2 . (11)

Scale transition scheme
In the following, a distinction between two different scales, the micro-and the macroscale is made.The former is characterized by a heterogeneous structure consisting of matrix The Cayley-Hamilton theorem states that a second order tensor fulfills his own eigenvalue equation, e. g.,  3 −  1  2 +  2  −  3 1 = 0, where  1 ,  2 ,  3 denote the principal invariants of , cf.Eq. ( 9).Consequently, any power   with  ≥ 3 as well as the inverse  −1 are expressible in terms of  2 ,  and 1 [64].
Note that the more common expression  2 = 1 2 (tr 2  −tr  2 ) is equivalent to  2 = tr(Cof ).Both expressions can be transformed into each other by using the Cayley-Hamilton theorem.length ℓ ∈ R + , whereas the second considers a macroscopic body B0 ⊂ R 3 of characteristic length l ∈ R + and is assumed to be homogeneous.For the lengths introduced, the relation l ℓ known as scale separation must hold [3].To label macroscopic quantities, they are marked by ( •) in the following.
In order to connect microscopic and macroscopic tensor quantities, an appropriate homogenization scheme is needed.Consequently, each macroscopic point X ∈ B0 gets assigned properties resulting from the behavior of the microstructure.For this purpose, a representative volume element (RVE) of the material is considered on the microscale in the vicinity of X, cf.Fig. 1.An effective macroscopic quantity is then identified from the microscopic field distribution within the RVE by the volume average Using Eq. ( 12), the macroscopic deformation gradient and the 1st Piola-Kirchhoff stress are defined by F :=  and P :=  , respectively [3,4,11].Appropriate boundary conditions for the microscopic BVP, which has to be solved before the volume averaging can be performed, are deducible from the equivalence of the macroscopic and the averaged microscopic energies which is also known as the Hill-Mandel condition.For the considered finite strain setting it is given by the following relation [3,4,11]: where (•) denotes the material time derivative.Regarding purely hyperelastic behavior, the Hill-Mandel condition expresses the equality of the rates of the macroscopic and the averaged microscopic Helmholtz free energies: ψ =  .Consequently, it holds ψ =  .To fulfill Eq. ( 13), several type of boundary conditions (BCs) can be used, whereby periodic BCs given by the spaces Preprint submitted to Springer arXiv are applied within this work [11].In the equation above,  =  •  is the nominal stress vector.Furthermore, (•) + and (•) − denote values on opposing boundaries of the RVE which is supposed to be periodic.The tilde ( •) marks the fluctuation part of a microscopic tensor quantity.

ANN-based multiscale approach with autonomous data generation
Algorithm 1 Procedure of the data-driven multiscale framework FE ANN given as pseudo code.
Based on the summarized continuum theory, the following section introduces the proposed data-driven multiscale scheme FE ANN .The general procedure of this framework is basically subdivided into five steps referred as After the framework is initiated with step (a), the steps (b)-(e) are executed within an overall loop which is given in Algorithm 1. Accordingly, the ANN is trained with the current data set D  of the iteration  ∈ N, the macroscale problem is solved, unknown macroscopic states of deformation F are detected, and the macroscopic stresses P corresponding to the previously unknown deformations are generated via computational homogenization.The loop terminates in an iteration  ≥ 1, after no further deformation states are found and the relevant space of deformation is thus completely sampled by D  .In order to prevent that the macroscopic simulation does not reach the final time step  goal -which could be occur due to a failed convergence of the ANN -but no new deformations are detected within   ∈ { 0 ,  1 , • • • ,   =  end <  goal }, a further inner repeat loop is implemented.Therein, the steps (b)-(d) are performed multiple times, until the final time step  goal is reached or new states of deformation are found.A schematic representation of the framework is given in Fig. 2. The single steps are described in detail in the following.
In order to enable a fully automatized utilization of the framework, a Python wrapper which runs on a high performance cluster (HPC) using the Batch-System SLURM (Simple Linux Utility for Resource Management) has been implemented.This wrapper processes the loop independently by starting the individual jobs and managing the results.

Initial data generation (a)
To start with, an initial data set D 1 := { 1 T , 2 T , . . .,  T } consisting of  ∈ N data tuples  T := (  FRVE ,  PRVE ) ∈ L 2 × L 2 is generated.Basically, a low number of tuples is sufficient to initiate the multiscale scheme.Suitable states of deformation are simple load cases as, e. g., uniaxial tension and equibiaxial tension or simple shear where these load cases could be applied on the RVE in different directions.In the equations above, λ ∈ R + and γ ∈ R denote prescribed effective stretches and shears, respectively.The labeling with (−) means that the corresponding coordinate of the effective 1st Piola-Kirchhoff stress P is prescribed to zero.The stresses  PRVE belonging to the deformations  FRVE are calculated from a computational homogenization according to Subsect.2.3, whereby this is done by using an in house Matlab FE code.The periodic BCs are applied via the master node concept therein [69].

ANN-based macroscopic surrogate model
Within the data-driven multiscale loop, an appropriate macroscopic surrogate model is needed.To this end, a physicsconstrained ANN model, which a priori fulfills several physical conditions, is used, compare Subsect.2.2.
Assuming that the macroscopic material symmetry group of the considered composite is known, and thus the structural tensors M 2 , M 4 corresponding to the material are available, the invariants Ī ( C, M 2 , M 4 ) are determined first.These scalar values are mapped into a normalized domain, i. e.,   : R → [−1, 1] ⊂ R, Ī ↦ → ī ∀ ∈ {1, 2, . . ., } with respect to the training data set, and are arranged in a vector ī = ( ī1 , ī2 , . . ., ī ) ∈ R ×1 .Thereafter, the free energy is predicted by an FNN with the normalized invariants ī serving as input values and ψ as output.Restricting the network architecture to only one hidden layer containing  neurons, it holds  where the monotonously increasing and convex Softplus activation function [37]  : R → (0, ∞),  ↦ → () := ln(1 + exp()) ( 19) is used.The stress prediction is finally done by applying Eq. ( 8): Thus, the surrogate model automatically fulfills several physical principles: thermodynamic consistency, objectivity and material symmetry.A graphical summary of the described surrogate model structure is given in Fig. 3.
In order to fulfill additional physical conditions, it may be necessary to incorporate further non-independent invariants Note that a non-bounded activation function is necessary to fulfill the growth condition, i. e., ψANN ( C) → ∞ as ( J → ∞ ∨ J → 0 + ).Besides the Softplus activation function, further choices are possible, e. g., ELU or ReLu.However, in contrast to the Softplus function, they are not twice continuously differentiable.

Ī *
into the argument list of ψANN .In this case, Eqs. ( 18) and (20) are to be modified according to and where A := { 1 , . . .,   },  = |A|, is a set containing the indices of the additional invariants Ī *  .Accordingly, the growth condition can be additionally guaranteed for J → 0 + by including Ī * 3 := 1/ Ī3 as a further invariant and not only Ī3 itself, where Ī * 3 is defined to be independent, i. e., Ī 3 Ī * 3 = 0. Thus, A = {3} in this case.However, in addition, further constraints have to be satisfied by the weights to ensure that the growth condition is fulfilled.As shown in Appendix A, a possible sufficient condition is given by Therein, the set N := {1, 2, , . . .,  } contains the indices of the hidden layer neurons.In the work [37], the growth condition is fulfilled in a similar way by adding an additional energy term which is not directly included in ψANN .

Training process of the ANN (b)
The ANN-based model is trained with respect to the current data set D  in each iteration  of the multiscale loop, where a random division into training and test data is made.
Within the training process, the weights and bias values ,   ,   ,   and  *  are then determined.In order to enable an adjustment of the ANN to the stress values, gradients of the output with respect to the input are inserted into the loss which is similar to [38,36,39].The training is done by using the SLSQP optimizer (Sequential Least Squares Programming).Thereby, the ANN is trained several times, where the parameters of the best achieved training state are stored at the end [29].An implementation of the described workflow is realized using Python, Tensorflow and SciPy.Within the training the constraint (23) can be switched on optionally.Finally, in order to fulfill the normalization condition, the bias value  is chosen such that ψANN = 0 within the initial (undeformed) state.
Due to local minima within the loss function, the optimization procedure which is applied here depends on the starting values of the weights and biases.Thus, the network is trained several times to overcome this.

Implementation and macroscopic simulation (c)
The model equation (18) has been implemented within the FE toolbox FEniCS [70,71].Therein, the stress relation given by Eq. ( 20) and the material tangent which is required within the solution via a standard Newton-Raphson scheme, are calculated by means of automatic differentiation.
In addition to weights and bias values of the trained ANN, problem specific structural tensors M macro 2 , M macro 4 have to be prescribed within the macroscopic simulation, cf.Subsect.2.2.2.Note that the preferred directions on the macroscale does not necessarily have to match the ones of the RVE.A conversion which is necessary in this case is done in step (e), cf.Subsect.3.3.2.

Autonomous data mining
Besides the physics-constrained surrogate model, the core feature of the data-driven multiscale framework is the data mining process which is done in a fully automatic manner.

Data analysis (d)
Within each iteration of the overall multiscale loop, the local deformation states of the macroscopic body B are collected.To this end, the deformation gradient F is stored at each quadrature point of the FE mesh and at each time increment   .Consequently, the body's full state of deformation is characterized by the set F macro Now, previously unknown deformations have to be detected, whereby it is meaningful to only consider states providing additional information for the material under observation.Since the intrinsic constitutive behavior of the material lives in the space of invariants I ⊂ R (+)×1 , it is useful to perform a transformation into this space at this point [29]: are used, respectively.Thus, by making use of Eq. ( 26), the current data set D  is compared to F macro  within the invariant space.If a state which is contained in F macro  is unique within a given tolerance  tol,1 , it is needed to enrich the data set for the next multiscale iteration step  + 1. Thereby, the set F macro  is searched in an reverse manner, i. e., starting from the last time step  end .If a unique state is identified in the step   , the full time series F( 0 ), F( 1 ), . . ., F(  ) with   ≤  end is added to F new  ⊂ L 2 .Thus, it is possible that deformation states in the set F new  are multiple with respect to F new  itself and/or D  .To avoid that this unnecessarily inflates the training process, a Note that it is useful to save the time series of a new deformation state, i. e., all states F ( 0 ,  1 , . . .,   ) proceeding at a fixed quadrature point.This facilitates the application of the macroscopic deformation within the computational homogenization later on.Furthermore, an extension to path dependent materials requires this mandatory.
Preprint submitted to Springer arXiv further filtering step is performed in the data enrichment step (e).
By using the described technique for the identification of relevant macroscopic deformation states, it is possible to reduce the number of time consuming microscale simulations following in the next step to a minimum.

Data enrichment (e)
After F new  has been determined, it is necessary to identify the stresses belonging to the states  Fmacro .This is done by applying these states within the RVE simulations and calculating the stresses by volume averaging.In order to use only one RVE and anyhow allow different microstructure orientations in the macroscopic body, a transformation of  F is necessary before.To this end, the relation ψ( F) = ψ( • F •  T ) with  ∈  + , which holds for arbitrary material symmetry groups, is used.Thus, the collected state of deformation  Fmacro from the macroscopic sample is processed by the Euclidean transformation The tuples  T = (  FRVE ,  PRVE ) consisting of applied deformations and corresponding stresses are then collected in the set D RVE  .As already mentioned above, it is possible that data tuples in the set D RVE

Examples
In order to demonstrate the ability of the developed datadriven multiscale approach FE ANN described in Sect.3, it is applied to the solution of three numerical examples within this section.Specifically, three macroscopic structures -a cuboid with a circular hole, a torsional sample and the Cook membrane -are considered.All of them consist of a fiber reinforced composite revealing a highly nonlinear behavior of the individual phases.

Microscopic properties of the composite 4.1.1 Constitutive behavior
The constitutive behavior of the composite's individual components, i. e., fibers and matrix, is described by a hyperelastic Ogden model [66].It is given by the free energy density function x 1 x 2 x 3 (a) x 1 -direction x 2 -direction x 3 -direction where,   ,   ∈ R and  ∈ R + are model parameters.The symbol R +  iso  :=  −1/3   denotes the isochoric principal stretches following from the Flory split [72]: i. e., the multiplicative decomposition of  into volumetric  vol and isochoric  iso parts.Note that the introduced material parameters are related to initial shear modulus  init and Poisson's ratio  init via the relations respectively.Furthermore, in order to guarantee a physically meaningful behavior, the parameters   and   have to be restricted by the following constraints: . The stress of the Ogden model is determined by using Eq. ( 5) and follows to The calculation of the projection tensors   related to the right Cauchy-Green deformation tensor  is given in Eq. ( 3).Within the following numerical examples, the material parameters given in Tab. 1 are chosen.With this choice, a highly nonlinear behavior of the matrix phase is achieved, cf.[29].The parameters of the fibers are chosen in such a way, that a neo-Hookean model results.Note that the selected parameters are not related to a real material.

Microstructure and RVE definition
Besides the constitutive behavior of the individual components, the geometric arrangement of the individual phases has to be defined.In order to choose a realistic microstructure, it Preprint submitted to Springer arXiv Table 1: Material parameters for matrix and fiber phases of the composite described by the Ogden model (31).Initial shear modulus  init and Poisson's ratio  init as well as parameter sets   ,   and .The parameters of the matrix phase are chosen according to Kalina et al. [29]. is characterized by a monodisperse stochastic arrangement of fibers with  = 30 % volume fraction here.Furthermore, a minimum distance of  := 0.4 with respect to the fiber radius  = 10 µm is chosen.
According to the homogenization concept, a suitable RVE which is sufficiently large to capture the essential statistic properties of the microstructure is needed.Herein, this condition is checked by applying the  2 -test proposed by Gittman et al. [73].Thereby,  randomly generated RVEs with fixed number of inclusions and volume fraction are generated.The effective response of these RVEs is then determined by a computational homogenization, whereby a specific load case is chosen depending on the property  of interest, e. g., this could be a stress or a stiffness component.Based on these results, a statistical analysis is done by evaluating the scatter of  by means of the quantity If the accuracy of  is sufficient, i. e.,  2 ≤  tol , the tested sample size is the final RVE size.Otherwise, the sample size have to be increased and the analysis will be repeated.
The described  2 -test has been applied to the fiber reinforced material.Thereby, a tolerance of  tol = 2.5 % and a number of  = 10 RVEs for each RVE-class with  inc inclusions have been chosen.As representative load cases, uniaxial tensions according to Eq. ( 16) 1 are applied into the  1 -,  2and the  3 -direction, i. e., perpendicular and parallel to the fiber orientation.Thereby, the stress into the tension direction is evaluated according to Eq. (32).The prescribed tolerance is finally reached for the RVE-class with  inc = 100 fibers.The determined stress-stretch curves and the chosen RVE are depicted in Fig. 4 for the final RVE.Accordingly, a highly nonlinear and anisotropic effective response occurs.
Generation and meshing of the periodic cells have been done by using the python tool gmshModel , whereby the placement of the fibers is realized via the random sequential adsorption (RSA) algorithm [74].The RVEs were meshed by tetrahedron elements with 10 nodes.A total of 94,456 elements is reached for the final RVE with 100 fibers.
It should be noted, that it is not sufficient to represent the microstructure of a fiber reinforced composite with an overall transversely isotropic behavior by a hexagonal unit cell if finite strains are considered.This is due to the loss of the material symmetry which results from the deformation of the RVE cell, cf.Appendix B.

Application of the data-driven multiscale framework for the simulation of macroscopic samples
After the definition of the composite's microscopic properties, the data-driven framework FE ANN is applied for the simulation of several multiscale problems.The macroscopic geometries and applied BCs are depicted in Fig. 5.
In all examples, the following specifications and meta parameters have been used: The effective transversely isotropic behavior, which results from the RVE homogenization, is described by using the set Ī := ( Ī1 , Ī2 , Ī3 , Ī4 , Ī5 , Ī * 3 ) ∈ R 6×1 consisting of six invariants, cf.Subsect.2.2.3.As discussed in Subsect.3.2, the non-independent invariant Ī * 3 = 1/ Ī3 is added to allow the growth condition to be satisfied by construction of the network architecture.As already mentioned, this requires additional constraints for the weights, e. g., Eq. ( 23).However, training under these constraints results in noticeably worse stress predictions.Thus, in order to achieve maximum prediction quality, the constraints were not considered within the training here.The adapted architecture is nevertheless used to allow a better comparability to the same network architecture trained with the constraint, see the study in Appendix A. The ANNs used as the surrogate model consist of only one hidden layer with 15 neurons which has been shown to be sufficiently accurate, where the networks are trained 25 times in each macroscopic iteration.For comparison see the study given in [29].
As shown in Fig. 4(b), the fiber orientation of the RVE is given by  RVE =  3 , where  3 is the Cartesian base vector.In order to transform deformation states for the general case  macro ≠  RVE by using Eq. ( 27), Rodrigues' rotation formula is applied to determine .Therein,  :=  macro ×  RVE is the unit normal vector and  := (  macro ,  RVE ) the angle between the fiber directions.
Finally, the tolerance for the detection of unique deformation states within the invariant space I is set to  tol,1 = 5 %.The tolerance for the filtering step is set to  tol,2 = 1 %, cf.Subsects.3.3.1 and 3.3.2.

Cuboid under tension
As a first example, a cuboid with a circular hole is loaded by tension.To this end, a displacement of û = 0.4 L 1 is prescribed at the top surface, where L 1 is the initial length in  1 -direction.The fiber orientation is chosen to  macro =  1 .The cuboid's geometric dimensions are specified by L 1 × L 2 × L 3 = (100 × 100 × 25) mm.The hole in the center has a radius of 30 mm.

Preprint submitted to Springer arXiv
Table 2: Load cases for the generation of the initial data.The effective deformations FRVE are prescribed according to Eqs. ( 16) and (17).

Load case λmax λmin γmax directions
Multiscale iterations The single iterations within the loop are described for the multiscale simulation of the cuboid in the following.
To start the algorithm, the initial data have to be generated first, cf.Subsect.3.1 step (a).Here, a total number of 18 load cases -six uniaxial tension and compression, six equibiaxial tension and compression as well as six simple shear states -are prescribed to the fiber reinforced RVE.Different loading directions are considered to collect knowledge about the composite's overall anisotropy.The applied stretches and shears as well as the directions are given in Tab. 2. In order to avoid that data providing the same physical information are contained multiple times, a filtering process is applied in the invariant space I .Within I , the collected states cover only a sparse region which is pervaded by a few paths consisting of 193 tuples  T .This is shown in Fig. 6, where the set is exemplarily visualized in the sectional planes Ī1 -Ī3 and Ī1 -Ī5 .However, due to the physical knowledge which is incorporated into the ANN-based macroscopic surrogate model, this sparse data set is sufficient to initiate the multiscale scheme.
the algorithm enters the multiscale loop with iteration 1 and the initial data set is labeled as D 1 .It is used to train the physics-constrained ANN which is afterwards utilized as the RVE surrogate model within the macroscopic simulation, cf.Subsect.3.2 steps (b) and (c).Note that the prediction of the sample's macroscopic fields, i. e., displacement ū( ,   ), deformation F( ,   ), stress P( ,   ), etc., is only an initial guess in this first iteration.This is due to the fact that the ANN does not know the occurring stress-deformation states for the most part at this point.However, it is still sufficient to identify new deformation states.As described in Subsect.3.3.1,this is done by transforming the sample's deformation states, collected into the set F macro 1 , into the invariant space I and to identify unique states by a comparison to D 1 .As already mentioned, the full time series F( 0 ), F( 1 ), . . ., F(  ) with   ≤  end is added to F new  if a unique state is identified in the step   .These deformation paths are prescribed in the RVE Due to the material symmetry (transverse isotropy), several loadings are equivalent from the point of the material, e. g., uniaxial tension in  1 -or  2 -direction.Note that the fiber direction of the RVE is  RVE =  3 .simulation and the corresponding stresses P are determined.
After the subsequent filtering process with the tolerance  tol,2 , a total number of 495 unique tuples is identified.They are given in Fig. 6.
In iteration 2, the new information are used to enrich the data set This set is used to train the ANN which is then applied for the prediction of the sample's macroscopic fields.As shown in Fig. 6, only a low number of new tuples is now detected from F macro 2 .This is due to the fact that the macroscopic simulation terminates already in time step However, using these new data, the prediction of the macroscopic simulation in the subsequent 3rd iteration yields that the sample's deformation field maps to a new region within I which consists of 60 tuples.The described process -consisting of the steps training process, macroscopic simulation, data analysis and data enrichment -is repeated until no further unique data are found which is the case in iteration 5. Note that an overlap of tuples from D  and D new  is possible, since a state is already unique if it deviates only in one of the six relevant invariants.Thus, an intersection in the depicted sectional planes may occur.
In advance, the macroscopic sample's states F macro 5 → D macro 5 mapped to I and the data set D 5 of iteration 5 are depicted within the full invariant space in Fig. 7.It becomes apparent, that D 5 \ D 1 contains tuples which are not included into F macro 5 .Thus, as described above, the predictions of iteration 1-4 are only necessary to gradually approach the correct result of the macroscopic simulation.However, due to the fact that the collected deformations are prescribed within the RVE simulations to determine the corresponding stresses, no defective information could enter the data set.
The stress predictions of the ANN for all tuples in D 5 are compared to the reference values of the RVE simulations in Fig. 8(a).As shown there, an almost perfect prediction occurs for both, training and test data.A single loading path and the corresponding deformed RVE at time  goal are depicted in Fig. 9 exemplarily.
Validation In order to demonstrate the quality of the developed FE ANN approach, the full deformation space of the cuboid is analyzed and the single states are applied to the RVE simulation, respectively.To this end, the symmetric right Cauchy-Green deformation tensors  C are calculated from Preprint submitted to Springer arXiv     The unknown states are prescribed within RVE simulations to get the corresponding stress values.These stresses are then compared to the predictions of the ANN which has been trained by D 5 .As shown in Fig. 8(b), an almost perfect prediction is observed also for these states.Consequently, the proposed FE ANN approach has shown to be highly accurate.
Furthermore, in order to demonstrate that the chosen tolerance  tol,1 = 5 % within the multiscale loop is sufficient, the achieved macroscopic solution is compared to a solution in which this tolerance is prescribed to 2.5 %.As shown in Fig. 11, a relative error of below 2 % with respect to the reference solution occurs for the stress P11 In want of a fully coupled FE 2 implementation, this comparison is used.

Torsional sample
After the detailed description of the multiscale loop within the solution process of the cuboid geometry, a further example is considered.Now, the FE ANN scheme is applied to the solution of a torsional sample with a circular hole given in Fig. 5(b).Accordingly, the fiber direction of the macro sample is chosen as  macro =  2 in this example.Nevertheless, due to the applied transformation given by Eq. ( 27 In order to minimize the computational effort, the scheme is initiated by using the collected data set D 5 from the previous example as the initial data.In this way, a knowledge base is created for a specific material under consideration, which can be used for further simulations and, at the same time, can be continuously expanded.
The multiscale loop now terminates after only 2 iterations which is due to the described initiation with the available data set.As shown exemplarily for the sectional planes Ī1 -Ī3 and Ī1 -Ī5 in Fig. 12, a wide range of relevant states is already covered by the states extracted from the cuboid simulation.Thus, the advantage of a transformation into the invariant space becomes again very clear.Although the deformation of the cuboid tensile specimen and torsional specimen is very different in the C-space, both overlap clearly in invariant space which is in accordance to [29].Caused from the geometry and the anisotropic nonlinear elastic behavior, a complex deformation of the macroscopic sample occurs, cf. the surface plot of C12 on the deformed geometry in Fig. 13.

Cook's membrane
Finally, as a last example, Cook's membrane is simulated by using the developed data-driven multiscale scheme.Thereby, the fiber direction within the membrane is prescribed to To initiate the scheme, the data set from the previous example is used, where this set also contains data collected within the first example.Again, a fast convergence of the multiscale loop is achieved after only 2 overall iterations.As shown exemplarily for the sectional planes Ī1 -Ī3 and Ī1 -Ī5 in Fig. 14, a wide range of relevant states is already covered by the states extracted from the cuboid and the torsional sample simulations.Thus, as already mentioned, the advantage of a transformation into the invariant space is underpinned.Although the deformation of the cuboid tensile specimen, the torsional specimen and the Cook membrane are very different in the C-space, they overlap clearly in invariant space.The deformed macroscopic states with the right Cauchy-Green deformation C11 are depicted in Fig. 15 for the time steps  10 ,  15 and  25 , where the load p is applied linear within the steps   ∈ { 0 ,  1 , . . .,  25 }.Due to the obliquely oriented fibers with respect to the alignment of the Cook membrane in the  1 - 2 -plane, an out of plane    deformation of the Cook membrane occurs.This effect results from the coupling between shear and tension which is a well known effect of fiber reinforced materials.Although a relatively complex response behavior occurs here, the multiscale problem can be solved quickly and without further human supervision even in this example.

Conclusions
In this work, a novel data-driven multiscale approach called FE ANN is presented.It is based on physics-constrained ANNs which are used as highly efficient surrogate models and an unsupervised data mining process.The approach allows the efficient simulation of materials with complex underlying microstructures which reveal an overall anisotropic and nonlinear elastic behavior on the macroscale, e. g., composites, architectured materials with pronounced microstructure or foams.
The framework has been implemented in such a way, that it is usable on a HPC cluster based on the Batch-System SLURM.
Starting from basic kinematics and stress measure definitions, a short revision of anisotropic hyperelastic constitutive models at finite strains is given.Furthermore, a Hill-type homogenization framework is described in brief.Based on this theoretical basis, the developed data-driven multiscale scheme is illustrated in detail.This includes the general procedure and a description of the single steps: initial data generation (a), training process of the ANN (b), macroscopic simulation (c), data analysis (d) and data enrichment (e).Afterwards, the approach is exemplarily applied to the solution of three demonstrative examples, a cuboid under tension, a torsional sample and the Cook membrane.Thereby, the considered macroscopic bodies consist of a fiber reinforced composite revealing a highly nonlinear behavior of the individual phases.Due to the incorporation of physical knowledge into the ANNbased surrogate model, only a small number of computationally expensive RVE simulations was needed to solve the considered macroscopic problems.Furthermore, a rather high accuracy of the surrogate model has been shown within a validation.
Altogether, the presented data-driven approach has shown to be an efficient tool for the solution of complex multiscale problems at finite strains.Due to the implemented unsupervised data mining, it is universally applicable to various macroscopic geometries and BCs.In order to extend the scheme's application area, several extensions have to be made in the future.For instance, an extension to further material symmetry groups [68] have to be made by integrating appropriate invariant sets into the implementation.Finally, an extension to dissipative constitutive behavior [41,56,55] is needed.

A Fulfillment of the growth condition by adjusted physics-constrained ANNs
Within this appended section, it is discussed how the growth condition, i. e., ψANN ( C) → ∞ as ( J → ∞ ∨ J → 0 + ) , can be fulfilled by construction of physics-constrained ANNs.Furthermore, a comparison of such an adapted network with a network neglecting the growth condition is given.

A.1 Network architecture and additional constraints
As discussed in Subsect.3.2, there are several requirements for the ANN to fulfill Eq. ( 35) by construction: First of all, a non-bounded activation function is necessary in any case.This requirement is fulfilled by using the Softplus activation function, cf.Eq. (19).Secondly, the additional invariant Ī * 3 := 1/ Ī3 has to be included in the argument list of ψANN .This is to enable the fulfillment of Eq. ( 35) for J → 0 + .Thirdly, additional constraints have to be satisfied by the weights belonging to Ī3 and Ī * 3 .Accordingly, a possible constraint is given by which is a sufficient condition to guarantee that Eq. ( 35) holds.In the equation above, the set N := {1, 2, , . . .,  } contains the indices of the hidden layer neurons.A proof of Eq. ( 36) is given in the following.Thereby, taking into account the requirements given above, the following network with one hidden layer is considered: Now, the purely volumetric deformation state F = λ1 is analyzed to study the model's behavior for J → 0 + and J → ∞.For this state, the relevant invariants follow to Ī1 = 3 λ2 , Ī2 = 3 λ4 , Ī3 = λ6 , Ī4 = λ2 , Ī5 = λ4 , Ī * 3 = λ−6 .
Additionally, the normalization of the invariants according to ī ( Ī ) := Ī − Īmax  + Īmin represented by a unit cell, whereas the second one consists of 100 fibers to capture for statistical effects.The fiber orientation points in the  3 -direction for both.Exemplarily, a uniaxial tension into the  1 -and the  2 -direction, i. e., perpendicular to the fiber orientation, are considered, where a maximum stretch of λ = 2 is applied.In Fig. 17(a) and (b), the stress-stretch curves and the transverse stretch λ⊥ are depicted for the hexagonal unit cell.Thereby, λ⊥ , which follows due to lateral contraction, is measured in the  2 -or  1 -direction, respectively.Since in the undeformed state all fibers have the same distance in the  1 - 2 -plane, the curves in the initial region are nearly equivalent.However, as the deformation of the RVE increases, the curves deviate more and more from each other.This is due to the fact that the arrangement of the microstructure changes significantly as a result of the deformation which could be termed a deformation induced anisotropy.Regarding the deformed microstructure for tension in the  1 -direction, it no longer corresponds to the arrangement in the case of tension in the  2 -direction (rotated by 90 degrees), cf.Fig. 17(c).Thus, in summary, the material loses the property of transverse isotropy if finite deformations occur.
Compared to this, the same uniaxial loadings are depicted in Fig. 18(a) and (b) for the cell with a random fiber distribution.As one can see there, the curves lie on top of each other over the complete range of stretch λ ∈ [1,2].Thus, in contrast to the hexagonal unit cell, a transversely isotropic effective behavior -which is expected for a fiber reinforced compositeresults even for finite strains, whereby the  1 - 2 -plane is the isotropy plane.In order to simulate the overall behavior of a realistic fiber reinforced composite, the usage of a statistical RVE is thus mandatory for finite strains.

Figure 1 :
Figure 1: Schematic depiction of micro-and macroscale within a multiscale problem.To enable scale separation l ℓ must hold.
(a) initial data generation, (b) training process of the ANN, (c) macroscopic simulation, (d) data analysis, as well as (e) data enrichment.

Figure 2 :
Figure 2: Schematic representation of the data-driven multiscale framework FE ANN : The framework starts with (a) initial data generation by homogenization.Afterwards, the steps (b) training process, (c) macroscopic simulation, (d) data analysis and (e) data enrichment are performed several times until no new data are added.Within the steps (d) and (e), the tolerances  tol,1 and  tol,2 are prescribed, respectively.

Figure 3 :
Figure 3: Structure of the macroscopic surrogate model based on a physics-constrained ANN.Note that   ∈ [−1, 1] only holds for states included in the known region of training data.

𝑖
are multiple with respect to D RVE  itself and/or D  .Thus, a further filtering step is performed, where multiple tuples are sorted out with respect to D  and D RVE  within a given tolerance  tol,2 .To this end, a transformation to the invariant space is used again, see Subsect.3.3.1.Finally, the enriched data set for the next iteration of the multiscale loop follows from D +1 = D  ∪ D new  , where D new  only contains relevant and unique tuples of the added deformations and corresponding stresses.

Figure 4 :
Figure 4: Fiber reinforced composite: (a) effective stress-stretch curves for uniaxial tension in different loading directions and (b) periodic RVE with 100 fibers and  = 30 % volume fraction.The size of the RVE is given by (323.6 × 323.6 × 10) µm and the fiber radius is 10 µm.

Figure 7 :Figure 8 : 3 Figure 9 :
Figure 7: Invariant space I of the cuboid under tension within the last multiscale iteration 5: (a) sectional planes with the full data set D 5 and the final prediction D cuboid 5 of the cuboid's deformation mapped into I as well as (b)-(d) field distribution of the invariants Ī1 , Ī3 and Ī4 on the deformed sample, respectively

Figure 10 :
Figure 10: Deformations occurring in the cuboid simulation: data set D 5 mapped into the C-space (final iteration) and completely unknown data.The fiber direction is given by  RVE =  3 .

|Figure 11 :
Figure 11: Stress field in the cuboid under tension: (a) macroscopic stress field P11 on the deformed configuration and (b) relative error of the determined stress field P11 with respect to Pref 11 obtained from a solution of the multiscale scheme with a reduced detection tolerance of  ref tol,1 = 2.5 %.Within the original simulation, the detection tolerance was  tol,1 = 5 %.

F macro 5 and included into the new set C macro 5 = { 1 Cmacro , 2 5 .
Cmacro , . . .,  Cmacro } ⊂  .(34) This set is then compared to the set D 5 within the space of the tensor C ∈ .Thereby, the contained deformation gradients  FRVE in D 5 are rotated back to  Fmacro by inverting Eq. (27).Subsequently, they are transformed to C-values to enable a comparison with C macro The unique states are then determined within a tolerance of 5 %.The states contained in D 5 and the unique, unknown states in the deformation space, i. e., C macro 5 \ D 5 , are depicted in Fig. 10.Accordingly, a wide region of C macro 5 do not intersect with the training data set which results from the transformation into the invariant space during the multiscale loop.

Figure 12 :
Figure 12: Sampling of the invariant space I within the simulation of the torsional sample.The initial data set contains the load cases given in Tab. 2 and the data collected from the cuboid simulation.The multiscale scheme was terminated in iteration 2. Exemplarily, the sectional planes Ī1 -Ī3 and Ī1 -Ī4 are depicted.

Figure 13 :
Figure 13: Torsional sample: distribution of the macroscopic 1st Piola-Kirchhoff stress P13 and right Cauchy-Green deformation C13 at the final time step.The surface plots are given on the deformed configuration B.
), the same RVE with  RVE =  3 is used.The torsional sample is loaded by specifying a distortion of φ = 45 • around the  1 -axis.The sample's geometric dimensions are specified by L 1 × L 2 × L 3 = (200 × 100 × 100) mm.The hole in the center has a radius of 40 mm.

Figure 14 :
Figure 14: Sampling of the invariant space I within the simulation of the Cook's membrane.The initial data set contains the load cases given in Tab. 2 as well as the data collected from the cuboid and the torsional sample simulation.The multiscale scheme was terminated in iteration 2. Exemplarily, the sectional planes Ī1 -Ī3 and Ī1 -Ī4 are depicted.

Figure 15 :
Figure 15: Cook's membrane: distribution of the macroscopic right Cauchy-Green deformation C11 at the time steps  10 ,  15 and  25 , where  25 is the final time step.The surface plots are given on the deformed configuration B.

Figure 17 :
Figure 17: Uniaxial tension of a hexagonal unit cell with fiber orientation in  3 -direction: (a) and (b) comparison of stress P and transverse stretch λ⊥ for tension in  1 -and  2 -direction and (c) deformed microstructure combined out of 3 × 3 and 3 × 2 periodic unit cells, respectively.The undeformed state is marked with a black dashed line.

Tension in x 1 Figure 18 :
Figure 18: Uniaxial tension of a stochastic RVE with fiber orientation in  3 -direction: (a) and (b) comparison of stress P and transverse stretch λ⊥ for tension in  1 -and  2 -direction and (c) deformed microstructure.The undeformed state is marked with a black dashed line.

Table 3 :
Multiscale iterations of the cuboid under tension: Tuples in the current data set D  and the new states collected in D new  as well as reached time step  end within the macro simulation.Iteration  Tuples in D  Tuples in D new t 10 /t 25 t 15 /t 25 t 25 /t 25