This work focuses on the numerical modeling of fracture and its propagation in heterogeneous materials by means of hierarchical multiscale models based on the FE method, addressing at the same time, the problem of the excessive computational cost through the development, implementation and validation of a set of computational tools based on reduced order modeling techniques.
For fracture problems, a novel multiscale model for propagating fracture has been developed, implemented and validated. This multiscale model is characterized by the following features:
The objectivity of the solution with respect to the failure cell size at the microscale, and the finite element size at the macroscale, was checked. In the same way, its consistency with respect to Direct Numerical Simulations (DNS), was also tested and verified.
For model order reduction purposes, the microscale Boundary Value Problem (VBP), is rephrased using Model Order Reduction techniques. The use of two subsequent reduction techniques, known as: Reduced Order Model (ROM) and HyPer Reduced Order Model (HPROM or HROM), respectively, is proposed.
First, the standard microscale finite element model High Fidelity (HF), is projected and solved in a low-dimensional space via Proper Orthogonal Decomposition (POD). Second, two techniques have been developed and studied for multiscale models, namely: a) interpolation methods, and b) Reduced Order Cubature (ROQ) methods . The reduced bases for the projection of the primal variables, are computed by means of a judiciously training, defining a set of pre-defined training trajectories.
For modeling materials exhibiting hardening behavior, the microscale displacement fluctuations and stresses have been taken as primal variables for the first and second reductions, respectively. In this case, the second reduction was carried out by means of the stress field interpolation. However, it can be shown that the stress projection operator, being computed with numerically converged snapshots, leads to an ill-possed microscale reduced order model. This ill-poseddness is deeply studied and corrected, yielding a robust and consistent solution.
For the model order reduction in fracture problems, the developed multiscale formulation in this work was proposed as point of departure. As in hardening problems, the use of two successive reduced order techniques was preserved.
Taking into account the discontinuous pattern of the strain field in problems exhibiting softening behavior. A domain separation strategy, is proposed. A cohesive domain, which contains the cohesive elements, and the regular domain, composed by the remaining set of finite elements. Each domain has an individual treatment. The microscale Boundary Value Problem (BVP) is rephrased as a saddle-point problem which minimizes the potential of free-energy, subjected to constraints fulfilling the basic hypotheses of multiscale models.
The strain flucuations are proposed as the primal variable for the first reduction, where the high fidelity model is projected and solved into a low-dimensional space via POD. The second reduction is based on integrating the equilibrium equations by means of a Reduced Order Quadrature (ROQ), conformed by a set of integration points considerably smaller than the classical Gauss quadrature used in the high fidelity model.
This methodology had been proven to be more robust and efficient than the interpolation methods, being applicable not only for softening problems, but also for hardening problems.
For the validation of the reduced order models, multiple test have been performed, changing the size of the set of reduced basis functions for both reductions, showing that convergence to the high fidelity model is achieved when the size of reduced basis functions and the set of integration points, are increased. In the same way, it can be concluded that, for admissible errors (lower than ), the reduced order model is times faster than the high fidelity model, considerably higher than the speedups reported by the literature.
BVP Boundary Value Problem
CSDA Continuum Strong Discontinuity Approach
DNS Direct Numerical Simulation
EBA Expanded Basis Approach
EFEM Embedded Finite Element Methodology
EIM Empirical Interpolation Method
FE Finite Element Method
FE FEFE hierarchical multiscale technique
HF High Fidelity
HPROM High-Performance Reduced Order Model
HROM Hyper-Reduced Order Model
MOR Model Order Reduction
POD Proper Orthogonal Decomposition
ROM Reduced Order Model
ROQ Reduced Optimal Quadrature
RUC Repeating Unit Cell
RVE Representative Volume Element
SVD Singular Value Decomposition
During the last decades, a large variety of multiscale strategies focusing on the study and analysis of the mechanical behavior of heterogeneous materials, have been proposed by the computational mechanics community. Based on the work of , these strategies may be divided into three main groups:
The formal treatments were provided by, , [6,7,8,9,10,11]. Two of the most relevant results of this kind of models are, the upper bounds of , and the lower bounds of . Posteriorly, Hill extends those bounds to tensorial entities, particularly, for constitutive tangent tensors .
The hypothesis of these approaches are properly fulfilled if there is a marked scale separation between the phenomena observed at the macroscale, and the ones observed at the microscale. However, nowadays, new approaches have been proposed to overcome this limitation.
Models based on the existence of a RVE can be divided into two main groups:
In virtue of the potential applications in microstructures with complex morphologies, the FE technique is suitable to deal with problems facing material instabilities, like fracture processes. In this sense, some approaches have been proposed [24,25,26,27,28,29,30], among others. One of the main contributions of this work, is an alternative approach with marked differences with respect to the previous ones.
Like the above mentioned models, concurrent models are widely applied. In approaches based on the Finite Element Method, the embedded meshes are not necessarily compatible. However, its computational viability is only for cases with small scale separation, this feature becomes into its main disadvantage. Some concurrent approaches in the field of fracture mechanics have been recently developed, , .
(1) Consistent in the sense that, all possible boundary conditions have to be compatible with the strains obtained at the macroscale.
The study and analysis of fracture in solids has been a topic of research since the last century. The seminal works on this topic were focussed on the Elastic Fracture Mechanics. However, its generalization to nonlinear material behavior is a non trivial task.
Starting from the concept of fracture energy, which has become a central issue in nonlinear fracture mechanics modeling, several techniques have been developed:
The study of heterogeneous materials subjected to softening, and, therefore, to degradation and failure, through multiscale approaches brings additional challenges. The fundamental reason lies in two aspects: (a) it becomes imperative the use of regularized constitutive theories at both scales in order to ensure the well-possednes of the multiscale problem. (b) The size effect, intrinsically related to the fracture energy, and extensively studied by . As a result of this, the homogenized stress tensor, in the post-citrical regime, becomes extremely sensible to the RVE size.
The second issue is the existence of the RVE , and the fulfillment of the basic hypothesis in multiscale modeling.
The necessity to develop specific homogenization techniques, becomes a starting point for obtaining consistent multiscale formulations.  has proposed a methodology that excludes the localization domain in the homogenization process. More recently,  proposed a predetermined size of RVEs.  describes a novel methodology, based on the existence of a macroscopic adhesive interface, which links the macroscale jump of displacements with an equivalent jump at the microscale, imposed by consistent boundary conditions.
Recent works [26,27], describe the material failure by means of nonlocal gradient theories. In this kind of approaches, and, in contrast with other alternatives, the homogenization of the stress field during the post-critical regime, is carried out at the localization zone (this zone corresponds to a subdomain of the RVE). However, other authors claimed some inconsitencies related to this kind of approaches, particularly, about the fact that kinematics at the macroscale is not equivalent to the kinematics modeled at the microscale.
In general, the FE method involving fine space-time discretization and time-dependent homogenization procedures, involves an enormous computational cost, being even larger when facing fracture mechanics problems.
The idea of exploiting the combination of dimensionality reduction and multiscale modeling is certainly not new. A survey of the related literature reveals that, over the last decade, researchers from various scientific disciplines dealing with multiscale problems have begun to consider the model reduction as a potential route –complementary to improvements in software and hardware power –to diminish the often unaffordable cost of multiscale simulations. In the specific context of homogenization-based multiscale methods, the application of model reduction techniques has been addressed by several authors, namely, [53,54,55,56,57]. The strategy adopted in all these works for constructing a cost efficient model of the micro-cell is the standard reduced basis method. The gist of this strategy is to project the governing equations onto a low-order subspace spanned by carefully chosen bases .
Reduced basis methods, in its standard form, suffer from an important limitation when handling nonlinear problems: they reduce notably the number of degrees of freedom –and thus the pertinent equation solving effort–, yet the computational cost associated to the evaluation of the internal forces and jacobians at quadrature points remains the same. Standard reduction methods prove, be effective only when dealing with micro-cells whose constituents obey simple constitutive laws (linear elasticity). In a general inelastic case, the calculation of the stresses at each gauss point is, on its own, a computationally expensive operation and dominates the total cost of the computation. As a consequence, the speed up provided by standard model reduction methods in nonlinear scenarios is practically negligible, and may not compensate the cost associated to the offline construction of the reduced-order bases.
The origin of the first effective proposal on this issue can be traced back to the seminal work of , who suggested to approximate the nonlinear term in the reduced-order equations by a linear combination of a few, carefully chosen basis functions. In the spirit of a offline/online strategy, in the standard reduced basis approach, these spatial bases are computed offline from full-order snapshots of the non-linear term, whereas the corresponding parameter-dependent modal coefficients are determined online by interpolation at a few (as many as basis functions), judiciously pre-selected spatial points. As in classical reduced bases methods, the efficiency of this second or collateral reduction is predicated on the existence of a moderate number ( is the original dimension of the problem) of basis functions whose span accurately approximate the manifold induced by the parametric dependence of the nonlinear contribution. The interpolation method developed by  is known as the EIM; the main ingredients of this method are: a) the use of a greedy algorithm to generate a set of maximally independent bases from the collection of snapshots of the nonlinear term, on the one hand; and b) the recursive selection – also via a greedy algorithm – of spatial locations where the error between the full-order bases and their reconstructed counterparts is maximum1.
In solution methods in which the governing equations are used in its variational form (as in the FE), the reduction of the computational complexity arising from nonlinearities can be, alternatively, achieved by approximating the integrals in which the offending nonlinear function appears, rather than the function itself, as done in the interpolatory and least-square reconstruction techniques discussed above. Based on this observation,  propose a quadrature scheme devised for fast-run integration of the subspace spanned by a representative set of snapshots of the nonlinear integrand.
In what follows, we shall consider as equivalent the appellations HPROM and HROM to refer to reduced basis methods combined with interpolatory or least-square reconstruction schemes.
The development of reduced models for non-homogeneous materials has been tackled in numerous previous contributions, such as , where the proposed reduction techniques are based on Fourier's transforms, or , where a reduced model is applied the homogenization analysis of hyperelastic solids subjected to finite strains. Also, the work in  develops a hyper-reduced model of a monoscale analysis which consider nonlinear material behavior. However, the existing literature barely considers reduced order modeling of non-smooth problems, as is the case of fracture, where discontinuous displacements occur. The multiscale case, when fracture takes place at both scales of the problem, makes the task even much harder. Indeed, only very few contributions have been presented in the literature about this topic, see for example: , which follows an eigendeformation-based methodology, or [63,64] that resort to global–local approaches.
The previous approaches combine projection techniques and, in some cases, empirical criteria to integrate the equilibrium equations in the domain. However, these are ussually ad-hoc techniques, that had been applied to problems with relatively simple crack propagation schemes. Currently, some researchers consider the effective model order reduction of fracture processes, an insolvable problem. This work will reconsider this statements, by developing a robust HPROM formulation, for multiscale fracture problems resulting in high computational speedups.
(1) Maximum in the sense that, the selected points have to be taken from components in which the error between the high fidelity and the HPROM solutions is greater.
The approach adopted in this work, uses a FE method and multiscale hierarchical models. Particularly the FE approach is used, which involves two scales (macroscale and micro/meso-scale) both discretized via finite elements. Infinitesimal strain setting, and first-order homogenization are assumed.
For fracture modeling purposes, the CSDA is adapted to the multiscale setting, and used for modeling propagating fracture at the macroscale level. At the microscale level, the use of predefined cohesive bands, distributed within the components and its interfaces, is proposed. These cohesive bands are endowed with regularized continuum damage models, which induce the crack initiation and propagation.
The Model Order Reduction techniques used in this work, are based on the POD, defining the projection of the full order model into a low-dimension small space, and, on the use of novel interpolation and ROQ schemes to diminish the computational cost generated by the multiscale problem.
The main objectives of this work are:
The remainder of this manuscript is organized in four chapters. Chapter 2 is devoted to the derivation of the multiscale model for propagating fracture, including a brief introduction to the fundamentals of the computational homogenization used in the proposed approach. Chapter 3 deals with the derivation of reduced order models for multiscale, smooth and non-smooth (fracture), problems, as well as some numerical results obtained from the developed models. Chapter 4 provides some concluding remarks and identifies areas for future research. In Appendix A, the participations in national and international conferences, and specialized workshops are listed. In Appendix B, a short summary of the supporting papers is presented. Finally, in Appendix C, the scientific publications supporting this work, and co-authored by the author, are annexed.
In the context of two scale (macroscale – micro/mesoscale) problems, computational homogenization of materials is generally regarded as a way of obtaining point-wise stress–strain constitutive models at the macroscale, accounting for complex micro/mesoscopic material morphology.
The homogenization approach used in this work –commonly know as first-order homogenization– is only valid for materials that display either statistical homogeneity or spatial periodicity.
In consequence, depending on the morphology and random distribution of constituents at the microscale, the definition and existence of a representative sample RVE plays an important role in the material characterization of heterogeneities at the macroscale.
This representative sample, hereafter denoted , is assumed to exhibit several features. One of those corresponds to the size indifference property [70,66,68], which states that if the size of this sample is increased, the response remains identical regardless the admissible boundary conditions on the RVE. The lower size limit for the RVE satisfying the size indifference property is represented by the characteristic length-scale denoted as , giving rise to the existence of the RVE, whereas in microstructures that display periodicity, is known as RUC, or simply unit cell. Furthermore, has to be small enough to be regarded as a point at the macroscale  (, , being the characteristic length of the macroscale , see Fig. 1) this is the so-called scale separation hypothesis.
This section presents a summary of the multiscale variational formulation used in this work. This approach is based on the following fundamental hypotheses:
|Figure 1: Macrosctructure with an embedded local microstructure.|
In the context of the adopted first-order homogenization setting, the microscopic velocity field can be splitted as follows:
where stands for the velocity at the macroscale, stands for the rate of infinitesimal macroscopic strain tensor, the term is a velocity term that varies linearly with , and the velocity fluctuations. The decomposition of the rate of microscopic strain tensor in the finite element framework yields, from the spatial differentiation of Eq. 2.1:
The starting point of multi-scale constitutive settings, is the assumption that the rate of macroscopic strain , at a point of the macro-continuum, is the volume average of the rate of microscopic strain 2.2, over the RVE associated with . This assumption is also interpreted as the fact that the microscale deformations only influence the macroscale behavior through its volume average.
In virtue of 2.2 and 2.3, this condition is equivalent to impose the volume average of the symmetric gradient of the velocity fluctuations to vanish. This condition can be written using the Gauss theorem as a constraint over , involving the whole volume of the RVE, as follows:
where is defined as the space of admissible microscale velocity fluctuations in the RVE, stands for the boundary of the domain , and is the unit normal vector on . Equation 2.4 is also known as the minimum constraint boundary condition.
The actual set of kinematically admissible velocity fields , together with the associated space of virtual kinematically admissible velocities at the microscale, denoted by , play a fundamental role in the variational formulation of the equilibrium problem of the microscale. This space can be defined as follows:
Furthermore, the same arguments can be applied to the total form, and establish that any kinematically admissible displacement fluctuation belongs also to .
The scale bridging equations are completed by introducing the Hill-Mandel Principle of Macro-Homogeneity [20,21]. Based on physical arguments, this Principle states that the macroscopic stress power equates the volume average over the RVE of the microscopic stress power, making both, macroscale and microscale, continuum descriptors energetically equivalent. Thus, departing from:
In particular, taking , and , yields:
where, stands for the macroscopic stress tensor, which turns out to be as the volume average of the microscopic stress . Equation 2.10 is also fulfilled in rate form. In addition to Eq. 2.10, the following condition emerges from the variational equation 2.9 solving for :
Eq. 2.11 defines the variational microscale equilibrium problem (or microscale virtual power principle).
In Computational Fracture Mechanics, hierarchical multiscale methods involve additional issues. In particular:
Along this section, the most important aspects of the proposed multiscale approach are summarized. This multiscale approach is fully detailed in Paper in Sec. 5.2.
|Figure 2: Macroscopic (Structural scale) body (a) subdivision in a non-smooth domain , and a smooth domain (b) h-regularized displacement and strain discontinuity kinematics.|
Considering the body , at the macroscale (see Fig. 2) it is assumed that material points, , of the macroscopic body belong, at the current time , to either one of the two subdomains:
where is the macroscopic displacement field, stands for the time or pseudo-time parameter, and , stands for the symmetric counterpart of .
In Eq. 2.20, stands for the regular (smooth) counterpart of the strain, is a colocation (characteristic) function on (See. Fig. 3), so that the term becomes a h-regularized Dirac's delta function shifted to the center-line, (the macroscopic discontinuity-path at the current time , as shown in Fig. 2-(a)). Thus, in Eq. 2.20, the term corresponds to the non-smooth (discontinuous and h-regularized) localized counterpart of the strains; a space-discontinuous second order tensor for the weak-discontinuity case.
|Figure 3: Colocation function|
Assuming that the fracture at the macroscale has arisen, in turn, by the appearance of failure mechanisms at the microscale level, originated by some type of material failure. The next step is to endow the microscale model with mechanisms to capture the onset and propagation of this material failure. Therefore, without introducing further details, it is considered that the microstructure shall be able to capture some dominant failure mechanisms of the material.
|Figure 4: Outline of the multiscale model for propagating fracture: a) macro and micro scales; b) microcell model accounting for material failure.|
For this purpose, a micro failure cell , of characteristic size , is considered to exist at every material point . It accounts for the material morphology at the microscale (voids, inclusions etc.). In addition, it is endowed with a set of cohesive bands () of very small width , whose position and other geometric properties (typically the normal , see Fig. 4) are predefined. At the current time , the activation (de-cohesion) of a number of those bands, defines the current subset of active bands which constitutes the "active" microscopic failure mechanism, for the considered point .
In principle, there is no intrinsic limitation on the number of the "candidate" cohesive bands to be considered at the failure cell. On one hand, their number and spatial position have to be sufficient to capture the dominant material failure mechanisms at the macroscale. On the other hand, the associated computational cost sets a limitation on the number of such bands. In this context, the following domains at the microscale are considered (see Fig. 4):
where and stand, respectively, for the stress and strain fields at the microscale point, , of the failure cell (corresponding to the macroscale point ), being the microscopic inelastic constitutive tensor derived from the hardening constitutive model.
where stands for a set of internal variables accounting for the inelastic behavior evolution.
where, in continuum damage models, with , and being the rate of the damage internal variable (a scalar for isotropic damage cases).
An advantage of this methodology, in the previous setting, is that the same failure cell morphology is considered to represent the microstructure at every macroscopic point of , both for and . The only difference is the considered constitutive behavior at the cohesive bands , defined in Eqs. 2.21,2.22 and 2.23.
Displacement fluctuations in the CSDA: Considering Eq. 2.1, with a local coordinate system () aligned with the domain (see Fig. 4), and, exhibiting the de-cohesive behavior allocated to the cohesive bands, the smooth part of the microscopic displacement fluctuation field, , can be expressed as:
|Figure 5: Cohesive Band behavior.|
where is the k-regularized Heaviside function shifted to , and is a (smooth) function arbitrarily defined except for the restriction in Eq. 2.25-(c), In Eq. 2.25 , is the apparent jump of across the cohesive band.
Following these statements, the microscale displacement fluctuation is given by (see Fig. 5):
Eq. 2.26 constitutes the displacement counterpart of a k-regularized strong discontinuity kinematics , and proves that the herein proposed cohesive-bands approach, is consistent with a k-regularized strong discontinuity at the cohesive domain . In consequence, the corresponding microscopic strain fluctuation field is given by:
where stands for the k-regularized Dirac delta function, placed at the center line, , of (see Fig. 5-(b)). Thus, the rate of microscopic strain field can be written in terms of the rate of macroscopic strain , and the rate of microscopic displacement fluctuations , as follows:
From Eq. 2.28, it can be concluded, that the second term at the right-hand side becomes unbounded in the limit . In multiscale modeling, this expression is equivalent to the one given, in phenomenological monoscale models, in the context of the Continuum Strong Discontinuity Approach (CSDA) of material failure .
(1) See: Sec.2.1 Paper 2
(2) See: Sec.2.2 Paper 2
One of the most specific features of the proposed multiscale approach, is that the same homogenization setting is used in points of both domains, smooth (), and non-smooth (), coinciding with the approach presented in Sec. 2.1. Other approaches , redefine the failure cell along time, fulfilling conditions of material bifurcation induced by instabilities at the microscale. More complex approaches [66,88,89,90] propose the use of second-order computational homogenization schemes in order to get better accuracy in the prediction of high strain gradients. In this work it is claimed the ability of the proposed approach to induce discrete failure in a first-order homogenization setting, giving rise to objective responses, and proper energy transfer through scales.
An issue appearing in this scenario, widely known in hierarchical multiscale approaches, is its high computational cost. In this context, the proposed model was also conceived to be combined with the use of model order reduction techniques (Paper ) . These techniques have been deeply studied in this work, and their main features are presented in Chapter 3.
In what follows, the consequences of the homogenization procedure based on the Hill-Mandel Principle of Macro-homogeneity are analyzed. The fact that the regularized strong discontinuities appear also at the microscale, being captured by the cohesive bands , is one of the most relevant features of the proposed approach.
|Figure 6: Multiscale model: (a) failure cell with activated failure mode; (b) geometrical characterization of the failure mode.|
For the sake of generality, the RVE is considered composed by several components: a matrix, and randomly distributed inclusions and voids. In addition, a number of cohesive bands are considered defining the set (a sketch is presented in Fig. 6); those cohesive bands allow failure within the matrix, 1 across the aggregates and at the matrix/aggregate interface.
In consonance with the definition of (in particular, the bounded behavior of the microscopic stress field), the second term on the right hand side can be neglected assuming a small enough width of the cohesive bands ().
Finally, after some manipulations of Eq. 2.30, and following the definitions of microscale kinematics in Eq. 2.28, and the lemma in Eq. 23 in Paper , the resulting homogenized constitutive equation fulfills the following:
where, stands for a characteristic length, depending on the activated microscopic failure pattern. The tensorial entities and , are inelastic strains, and play the same role than internal variables in phenomenological models. However, unlike them, here, 2 their evolution is determined, at every macroscopic sampling point , by homogenized values of entities at the corresponding microscopic failure cell . This extends to non-smooth problems, some theoretical results already derived for smooth problems, see [92,93]. In addition, a characteristic length emerges naturally in Eq. 2.36, as the ratio between the measure of the failure cell (area in 2D and volume in 3D), and the measure (length/surface) of the activated microscopic failure mechanism. In consequence this length is of the order of the failure cell size. For a deeper review of the analytical results of this induced homogenized constitutive model, the reader is addressed to Sec. 2.4 in Paper .
The role of the characteristic length, , naturally derived from the present formulation, is not only computational, but it has also other very relevant physical and mechanical implications. Consideration of such a characteristic length, for multiscale based approaches, has been claimed from the material mechanics community , and sometimes introduced in a heuristic way in other approaches . This characteristic length depends on both the specific data of the problem and the local microscopic failure state. Through its consideration, the correct energy transfer between scales and mesh size objectivity can be achieved.
In summary, Eq. 2.35 and Eq. 2.36 retrieve the format of a constitutive model equipped with an internal length and with internal variables whose evolution is described by the microstructure behavior. Although this model will never be used for computational purposes3, it supplies relevant insights on the properties of the resulting homogenized constitutive model.
(1) See: Sec.2.4 Paper 2
(2) See: Sec.2.4 Paper 2
Let us consider, on one hand the fracture energy, corresponding to points , defined as a material property specific for every compound of the heterogeneous RVE, and, on the other hand, the macroscale fracture energy , obtained as an output from the homogenization procedure. According to their definitions, those fracture energies can be computed in terms of fracture energy densities, in terms of the energy dissipation that takes place in bands with bandwidth (at the microscale) and (at the macroscale), respectively.
In virtue of the Hill-Mandel Principle of Macro-Homogeneity, 1 it can be concluded that the macroscopic fracture energy is equivalent to the average of microscopic fracture energy , along the activated failure mechanism at the microscale . Replacing Eq. 2.41 into Eq. 2.9, and after some manipulations, the macroscopic fracture energy is given by the expression :
where is the mean value of the microscopic fracture energy varying along the active failure path. Eq. 2.44 provides the relationship of fracture energies at both scales. In case of an homogeneous fracture energy at the active cohesive bands of the microscale, Eq. 2.40 translates into an exact equivalence of fracture energies along the scales, :
In the light of this result, it can be easily concluded that the fracture energies at the microscale determine, in average, the effective fracture energy at the macroscale. It is stressed the importance of the characteristic length in order to guarantee the proper dissipation between scales. For more details, the reader is addressed to Appendix B in Paper .
(1) See: Sec.2.6 Paper 2
The proposed multiscale formulation has been implemented in a Finite Element model following the setting of a FE strategy. Accordingly, two nested finite element models are used:
In what follows, these two finite element models are described.
Standard quadrilateral finite elements are adopted for the numerical simulation of the cell response. The cohesive bands are also modeled by quadrilateral isoparametric finite elements of very small thickness (high aspect ratio), as shown in Fig. 7-(a), endowed with constitutive models whose behavior is sketched in Fig. 7-(b) and defined through equations 2.21 to 2.23. The remaining finite elements of the cell are endowed with either elastic or inelastic hardening responses. Therefore, only elements on the cohesive bands can exhibit strain localization.
The corresponding nonlinear problem in the failure cell is then solved for the discretized version of the microscale displacement fluctuations, using Eq. 2.11. Dirichlet boundary conditions precluding rigid body motions, and minimal boundary conditions in Eq. 2.4, are also imposed.
|Figure 7: Multiscale model: finite element discretization at the microscale.|
Material failure propagates naturally through the RVE, strain localization takes place at the finite elements defining the cohesive bands. At every time step of the analysis, those finite elements who are in loading state, define the active set of cohesive bands conforming the active failure mechanism.
One of the most critical issues in computational modeling of material failure is the appropriate capture of the crack onset and propagation. When does failure trigger at a given material point? and how does it propagate?, these two questions are the cornerstone of material failure propagation algorithms.
At the microscale, where the morphology and the position of candidate propagation mechanisms are predefined, the two issues are of minor relevance due to the adopted simplified failure-bands model. However, at the macroscale, there is not a predefined failure path, and in principle, any material point may fail and propagate in any direction. To adequately solve the previous questions, the procedure for modeling onset and propagation of discontinuities recently developed for monoscale problems  has been extended to the multiscale setting. The proposed methodology is based on the use of the following specific techniques:
|Figure 8: Evolution of the injection domains for three typical stages () of the discontinuity propagation.|
|Figure 9: Sampling points involved in the numerical integration.|
The injected strain rate at element , with nodes, is the following:
|Figure 10: Weak discontinuity mode. Elemental regularized dipole function .|
In the present multiscale context, the proposed second stage consists of the incremental injection of the following elemental strong discontinuity mode:
in terms of the -regularized Dirac delta function (displayed in Fig. 11), being the direction of the element normal provided by the solution of the discontinuous bifurcation problem presented in Sec. 2.5 in Paper . The resulting variational problem for the injection procedure is summarized in Box A2 - Appendix B in Paper .
|Figure 11: Strong discontinuity mode. Elemental regularized Dirac delta function .|
The resulting procedure is a robust and efficient technique to model propagating material failure in a finite element setting. It is especially appropriate for capturing material failure propagation in coarse meshes, in contraposition of the alternative extra elemental character techniques (phase-field, gradient or non-local damage models), where several elements span the localization band. In addition, its implementation in an existing finite element code has a little intrusive character.
In regards to the space and time integrations, as commented above, injection of weak-discontinuity and strong-discontinuity modes requires, in principle, specific integration rules in space, : a standard four-point Gauss quadrature rule, and two additional sampling points, for injected elements, and so that . Since those domains evolve along time (see Fig. 8), some additional problems on the time-integration of the resulting equilibrium equations are found. To tackle this issue, in  and  is proven that defining some "equivalent" stress entities at the standard Gauss points, the spatial integration can be rephrased as a standard four Gauss points integration rule in the integration domain. This space-time integration rule is fully explained in Appendix B3 in Paper , and the corresponding stress evaluation is also summarized in Box A3.
(1) To switch between stages, a set of control variables are defined, all those detailed in Sec. 3 and Appendix B in Paper
(2) Under the CSDA, the homogenized dissipation is evaluated at the barycenter of the finite element, denoted by
Along this work, some techniques for reducing the unaffordable computational cost inherent to the numerical simulation of multiscale fracture problems have been developed. Those techniques are combined to obtain a hyper-reduced order model HPROM, based on a two-stage strategy:
In what follows, these techniques have been applied to the microscale BVP, while the finite element model at the macroscale remains as the standard one.
The model order reduction concept relies on the premise that, for any input parameter governing the microscale displacement fluctuations , the solution can be approximated by a set of linearly independent basis functions approximately spanning the primal variable1 space.
Following this idea, the off-line stage is devoted to determine via a POD technique, the reduced basis in which the HF solution is projected. Once this basis has been obtained, a subsequent online stage in the reduced-space is considered.
(1) Primal variable is known as the selected variable to perform the reduction process.
Taking as a primal variable the displacement fluctuations, and departing from the problem depicted in Sec. 2.1, a first step consists of determining an approximation1 of the finite element space of kinematically admissible microscale displacement fluctuations . This approximation is obtained as the span of the displacement fluctuation solutions obtained, for a judiciously chosen set of input strain trajectories, every trajectory being discretized into a number of steps . These set of finite element solutions are stored into the snapshot matrix as column vectors:
In consequence, the approximating space for , henceforth called the snapshot space, is then defined as:
where, is the total number of snapshots.
Once the snapshot matrix has been computed, the Elastic-Inelastic decomposition technique is used to determine the reduced basis functions. The reason for it relies on the fact that the SVD applied to the whole matrix , may produce basis with a large number of elements, which makes difficult to retrieve the response of the RVE in some specific cases. Particularly, the elastic response2, might request a much larger number of basis functions, this translating into a significant waste of computational cost.
To eliminate this shortcoming, in this work, it is proposed a time partition of the space of snapshots into elastic (), and inelastic () subspaces.
3obtaining the reduced basis as the combination (spatial sum) of both sub-bases. An orthonormal basis for is determined by taking a low number of elastic snapshots (at a minimum, for 2D problems, for 3D problems), and computing the corresponding orthonormal basis.
Once this set of elastic basis is known, the orthogonal projection of each snapshot onto the orthogonal complement of is computed; with this new set of snapshots, the inelastic basis functions are obtained via SVD. Finally, the assembled basis results the following:
and the reduced-order space , spanned by this base, is:
Placing the elastic modes in the first positions, followed by the essential4 inelastic modes, ensures the reduced-order model to deliver linear elastic solutions with the same accuracy than the HF solutions. For more details, the reader is encouraged to sent to the Appendix B in .
Once the reduced basis is computed, 5 the online stage consists of solving the discrete version of the microscale equilibrium equation (via FE), projected onto the reduced-order space spanned by . To this end, the test and trial functions, and , are approximated by the following linear expansions:
where, and stand for the low-dimensional approximations of trial and test functions, respectively.
PROBLEM A (ROM) (Microscale reduced problem via POD):
Given the macroscale strain, , and the reduced basis for displacement fluctuations , find satisfying:
where denotes the vector containing the coefficients associated to each basis function , being the basic unknowns for the standard reduced-order problem. stands for the reduced strain-displacement matrix “B-matrix” defined as . When using a Gauss quadrature integration scheme, is the total number of Gauss points of the mesh; denotes the weight associated to the g-th Gauss point ; and stand for the reduced B-matrix and the stress vector at Gauss point , respectively .
(1) In general, cannot be precisely determined, such a task will require finite element analyses of the cell under all conceivable strain paths. Rather, one has to be content to construct an approximation of it.
(3) See: App. B Paper 1
(5) See: Sec. 4 Paper 1
Taking advantage of the unbounded character of the microscale strain field typically observed in this kind of problems, the failure cell is splitted into a regular domain (made of elastic matrix and possible inclusions) and a cohesive domain (cohesive bands exhibiting a softening cohesive behavior). Details on this issue can be found in Sec. 3.2.2 in Paper .
In addition to this proposal, the ROM of the failure cell is formulated in an unconventional manner, : in terms of strain fluctuations rather than in terms of conventional displacement fluctuations.
As it will be shown later, it is convenient to rephrase the original problem, posed in terms of displacement fluctuations (PROBLEM-I in Paper ). The primary unknowns of the rephrased problem are now the microscale strain fluctuations instead of its displacement fluctuations, while the constrained original minimum problem (of the standard micro-cell BVP) is rewritten in terms of a Lagrange functional. The resulting formulation is a variationally consistent saddle-point formulation.
Considering the material free energy for the isotropic damage model in , the microscale stress field can be expressed as:
complemented by the evolution laws of the internal variables . Thus, in consonance with the hierarchical multiscale approach, the following parametrized functional can be defined:
where, , is a symmetric second order tensor Lagrange multiplier enforcing condition 2.4 on the microscale strain fluctuations . With this parametrized functional , a saddle-point problem can be stated as:
PROBLEM II (HF) (Microscale saddle-point problem): Given the macroscale strain, , find and satisfying:
where stands for the space of microscale kinematically compatible strain fluctuations and, stands for the evolution equation of the internal variables. After considering that the microscale stress field is given by Eq. 3.9, the following optimality conditions emerge:
Eqs. 3.12 and 3.13 provide the solution of the saddle-point problem stated in Eq. 3.11. It can be proven that Eqs. 3.12 and 3.13 make PROBLEM II equivalent to the original problem in Eq. 2.11, but now rephrased in terms of the microscale strain fluctuations (see PROBLEM I-R in Paper ).
The transition from the high-dimensional finite element space to the reduced-order space, is accomplished by applying the POD technique, now for non-smooth problems. The standard reduced order model is based on the reduction of the strain fluctuation field . The first step consists of generating a collection of solutions (samples) from different trial loading cases, representatives of all possible loading cases.
In each trial case, the microscale strain fluctuation at every Gauss point, is collected and stored in the snapshot matrix as a column vector:
where is the number of snapshots vectors. 1 Therefore, represents a number of sampled solutions obtained with the HF model under different loading conditions. For more details, the reader is addressed to Sec. 4 in .
In order to get a more accurate estimation of the dominant modes of the microscale strain fluctuations, it is convenient to separate the microscale into specific sub-blocks in accordance with the type of material response observed during the load history. The procedure is sketched in Fig. 12:
|Figure 12: (a) All entries of are partitions into two sub-blocks: and , the strain fluctuations from points at regular and cohesive domains, respectively. (b) Snapshots taken from the elastic regime of the failure cell correspond to the sub-block . Snapshots taken during the inelastic regime (at least one Gauss point at is in inelastic state) correspond to the sub-block .|
In addition, the Elastic-Inelastic snapshot decomposition above explained , 2 is also applied to each sub-block. Thus, without loss of generality, the snapshot entries are organized so that the first, entries correspond to Gauss points in , while the remaining entries correspond to Gauss points in .
In accordance with this criterion, the snapshot matrix 3.14 can now be partitioned into sub-matrices as follows:
where and collect the strain fluctuations located outside and inside the cohesive bands, respectively. The right hand side matrix in 3.15 emphasizes the double partition performed in accordance with elastic-inelastic regimes.
After some additional manipulations, the corresponding set of orthonormal basis is obtained as:
The number of reduced basis in is: , where the values and are obtained from the SVD applied to the projected inelastic snapshots. Additional details can be found in Sec. 3.2.2 in Paper .
Once the reduced basis is known, the strain fluctuations are interpolated as a linear combination of the elements of this basis as:
where each element , of the basis , is a microscale strain fluctuation mode and the vector of time dependent coefficients () represents their corresponding amplitudes (the actual unknowns of the problem). 3 In the same way, the variations of the microscale strain fluctuations are expressed similar to 3.18 as:
The problem solved in the online stage is then the following:
PROBLEM III (ROM) (RVE saddle point problem): Given the macroscale strain, , find and satisfying:
The optimality conditions for the problem above yield:
which, expressed in matrix notation, yield:
where is the column vector constituted by piling-up the stress vectors, , evaluated at the integration Gauss points. The column vector is also the pilled-up of repeated values of the same constants vector . 4 The square diagonal matrix , and the rectangular matrix , collect the Gauss weights: , which for plane strain cases () are distributed in sub-block matrices , as follows:
being a matrix with the corresponding -th Gauss weight placed at the diagonal.
The unknowns for the reduced order model are, the modal amplitudes , and the Lagrange multiplier .
Considering the system of equations 3.24 and 3.25 for and , it could be expected that this problem, of equations, should be less computationally costly, than the HF model. However, this is not the case: the numerical simulations with the ROM model do not substantially reduce the computational cost, and little (or even smaller than one) speedups are obtained. This fact highlights that the actual bottleneck for fast online computation is not the solution of the balance equations but, rather, the determination of the stresses, internal forces and stiffness matrices at every integration point of the underlying finite element mesh. Therefore, an additional technique is proposed to reduce the amount of integration points in which the constitutive equation is evaluated.
(1) See: Sec.4 Paper 3
(2) See: Sec.3.2.2 Paper 3
(3) See: Sec.3.1 Paper 3
(4) See: Sec. 3.2 Paper 3
For this purpose, a reduced integration technique has been developed by resorting to a nonconventional method, termed Reduced Optimal Quadrature (ROQ), to integrate the term involving the microscale free energy in 3.20:
Where stands for he ROQ.
The ROQ technique is based on selecting, from the initial set of “Gauss” sampling points, and through an adequate algorithm, an equivalent subset of sampling points , and their new corresponding weights . The success of the reduced integration numerical scheme, in front of the conventional Gauss quadrature, lies on the fact that it is possible to reduce notably the number of involved quadrature points to , being the number of integration points for the Gauss quadrature scheme, keeping under strict control, or even reducing to zero, the numerical error introduced by the reduced quadrature rule. Then, the microscale potential energy in Eq. 3.43, is re-expressed as:
In consequence, the corresponding optimality conditions (equilibrium equations) to be solved during the online stage are:
PROBLEM IV (HPROM) (Microscale reduced saddle-point problem): Given the macroscale strain, , find and satisfying:
A similar procedure could also be used for the integral terms (underlined as “Gauss quadrature”) in Eqs. 3.45 and 3.46. However, this would not produce a substantial computational cost gain due to the fact that those terms are constant (not depending neither on the unknowns of the problem nor on the constitutive internal variables). They are required to be integrated only once, via the standard Gauss quadrature, and the result can be stored, and retrieved when necessary, during the online stage execution.
(1) See: Sec. 5.2 Paper 1
(2) See: Sec. 5.2 Paper 1
(4) See: Sec. 5.3.2 Paper 1
(5) In this context, observed means calulated through the pertinent constitutive equation.
(6) The operator is the so-called selection operator associated to sampling indices , for instance, the restricted matrix of weighted strain modes is defined as
(7) The selection criterion used for the set of sampling points in this approach, is fully explained in Sec. 7 in , details about the optimality criteria, and its corresponding accuracy are deeply detailed in this section.
(8) See: Sec. 6.3 Paper 1
In spite that the goal of the ROQ is to develop a reduced cost interpolation scheme as a general framework for both static and dynamic problems, attention is focussed here on the multiscale quasi-static fracture problems. The minimum number of quadrature points providing an admissible integration error in the free energy integral, 3.43 is based on the optimal linear expansion of in terms of the free energy modes and its corresponding amplitudes . Thus, a similar expression to Eq. 3.18) can be adopted for constructing the reduced microscale strain fluctuations, as follows:
With the previous approximation in hand, the total microscale free energy can be expressed as:
In order to obtain the reduced optimal numerical quadrature rule, the following optimization problem is considered:
OPTIMIZATION PROBLEM : Given the expanded reduced basis , and the set of sampling points , find and satisfying:
Where, and stand for the error committed through the reduced integration of every free energy reduced basis function, and the error in the integration of the volume , respectively. The resulting algorithm (described in the flowchart of Box IV 1 in Sec. 5 of Paper ) returns a sub-set of optimal Gauss points, and the corresponding weights, that integrate exactly the basis and, therefore, the free energy in Eq. 3.47.
Regarding the computation of the microscale energy reduced basis , a SVD-based strategy is used in the off-line stage, similar to that described in Sec. 3.2.1.
The method is again based on the construction of a snapshots matrix, in this case, for the free energy, and the computation of its corresponding reduced basis via SVD. For this purpose, two options appear:
The first method is considered as the algorithmically consistent strategy. However, it is also more expensive than the second one. The reason relies on the fact that, in order to get the reduced basis for the microscale energy , training trajectories have to be computed twice: a) First using the HF model to obtain the strain modes , b) Second, using the ROM model to obtain the corresponding free-energy modes .
Both strategies have been tested and both provide accurate results. However, the later, being the cheaper and simpler one, was adopted as the most convenient.
In summary, both the strains and the free energies of the microscale are sampled simultaneously at the off-line stage, for different sampling trajectories with the HF model, and a series of snapshots of energy, , are evaluated and collected for each Gauss point. Then, the microscale energy snapshot matrix is built as:
In accordance with the position of the Gauss point,9 in the finite element mesh, and following a similar procedure to that adopted in Eq. 3.15, this snapshot matrix is also partitioned into components associated to the domains and as:
and the SVD technique is then separately applied to both partitions of to obtain two distinct (orthogonal) bases, for the elastic regime of both subdomains:
The corresponding inelastic reduced basis functions are also computed via SVD, following a procedure similar to the one described in Sec. 3.2.1. The complete reduced basis for the energy field , is made of the union of and :
The number of basis vectors in is: , where the values of and are obtained from the solution of the SVD applied to the inelastic projected snapshots.
The accuracy of the reduced models, ROM and HPROM, depends on several aspects. In order to assess it, three different sets of tests are done:
Details on this issue can be found in Sec. 5 in Paper .
A squared microscale model, made of a matrix and randomly distributed aggregates, is devised and tested (see figure 14) to simulate the microstructure of a cementitious-like material (concrete). Relevant details about the finite element model are presented in Table 1. To mimic the concrete material response, the failure cell is modeled with three components: aggregates, which are assumed to be elastic, bulk matrix, also assumed elastic, and interfaces (matrix-matrix and matrix-aggregates), which are modeled with cohesive-band equipped with an isotropic damage constitutive law. The properties of the components in the microscale are defined in Table 2.
|Figure 13: Failure cell|
|Number||Number||Number of||Total number of|
|of elements||of D.o.f.||Cohesive Bands||Gauss points ()|
|Cohesive bands of||0.18||2.60||140|
|Cohesive bands of||0.18||–-||–-|
In general terms, Figure 14 can be used as an “abacus” for a-priori selection by the user of the HPROM strategy in a multiscale problem (for a given microstructure at the RVE). For instance, by selecting the admissible error ( 3,5%) in the top figure, the number of strain modes is obtained. Entering in the lower plot, with this result (), one obtains the suitable number of integration points and the expected speedup .
The availability of a catalog (constructed off-line) for a specific RVE microstructure, allows the user's a-priori selection of the appropriated HPROM strategy, by balancing the admissible error vs. the desired speedup.
|Figure 14: HPROM design diagrams. Top: HPROM error in terms of the number of strain modes. Bottom: OQN and obtained speedup in terms of the number of strain modes. By selecting the admissible error (say 3,5%) in the upper diagram, one obtains the requested number of strain modes, Entering with this result in the lower diagram one obtains the suitable number of integration points and the resulting speedup (.|
The test shown in Figure 15 is a benchmark commonly used for testing macroscale propagating fracture models. This concrete-like specimen is considered here to test the qualitative results and convergence properties of the proposed HPROM approach, when utilized in real FE multiscale crack propagation problems.
|Figure 15: L-shaped panel: a) Specimen geometry; b) Finite element mesh|
The geometry of the simulated specimen is depicted in Figure 15-a. As shown in Figure 15-b, the domain of the L-shaped panel is split into two domains: 1) the multiscale domain (with elements) corresponding to the region where the crack is expected to propagate, modeled with the HPROM of the microstructure depicted in Figure 13, and 2) the remaining part of the panel, which is modeled with an elastic monoscale approach (using 1709 elements), where the elasticity tensor is obtained through an elastic homogenization of the micro-structure elastic properties. Even for this (rather coarse) multiscale problem, the high fidelity HF computational solution is extremely costly to handle, until the point that, with the available computational resources1, it was not possible to display the complete action-response curve (in Fig. 16).
However, the remaining structural responses in Figure 16, obtained through a number of HPROM strategies, involve very reasonable computational costs, and they were obtained in advance with no previous knowledge of the HF results. The accuracy is very good, and a response indistinguishable from the HF can be obtained times faster (). A less accurate response, but with a fairly good agreement with the HF can be obtained with .
|Figure 16: L-shaped panel: Structural responses in terms of force P vs. vertical displacement , for different RVE HPROM strategies, and obtained speed-ups.|
In Fig. 17, the evolution of the microscale crack opening is shown. It is worth noting that, both, the microscale failure mechanism and displacement jump vary along the macroscale in agreement with the crack propagation direction observed at the macroscale.
|Figure 17: L-shaped panel: microscale crack activation along the crack-path field, using and .|
This illustrates the new paradigm that is set and the computational possibilities open by the HPROM strategies in computational multiscale modeling explored in this work.
Multiscale modeling is foreseen to become a key approach to enable the next wave of design paradigms for engineering materials and structures. Indeed, it has an excellent potential to account for the physical links between different scales, involving the diverse phenomenologies intervening in the mechanical response of materials (grains, particles, defects, inclusions, etc.).
Quoting from a report by a group of experts to the US National Science Foundation :
". . . . In recent years, a large and growing body of literature in physics, chemistry, biology, and engineering has focused on various methods to fit together simulation models of two or more scales, and this has led to the development of various multi-level modeling approaches. . . .. To date, however, progress on multiscale modeling has been agonizingly slow. Only a series of major breakthroughs will help us establish a general mathematical and computational framework for handling multiscale events and reveal to us the commonalities and limitations of existing methods . . . .".
In this sense, the effort invested in developing and using multiscale models, has been, in many cases, fruitless, due to the involved computational cost in this kind of methodologies. This limitation becomes a bottleneck for multiscale modeling, usually discarded, or, relegated to the availability of supercomputers, and, therefore, not always accessible to the whole computational mechanics community.
In addition, while multiscale models exhibiting material hardening behavior have widely been studied, multiscale models dealing with material softening behavior are in an early stages of development.
Therefore, the development of a reliable, minimally intrusive multiscale fracture models becomes a crucial task, not only in order to have a robust and consistent multiscale fracture numerical tools, but also for developing their related reduced order models that allow their use in complex cases that can be used for industrial purposes, with an affordable cost. These are the fundamental reasons for the research and development about this issues.
A sketch of the overall work carried out in this work is shown in Figure 18. In there, contributions are chronologically numbered and highlighted with a blue arrow. Contributed papers are numbered from P1, corresponding to the the first contribution (Paper ), to P6 (Paper ) the last one; in this context, CB means Chapter in Book.
In what follows, they are specifically commented, and the corresponding conclusions and achievements, are presented.
|Figure 18: Global Flow Chart of the work|
This Article presents the first research developments in this work on MOR techniques applied to multiscale modeling. The scope of this publication is limited to smooth problems and exclude fracture processes. Techniques like interpolation methods via HPROM have been studied.
The concept of a two-stage reduction (ROM-HPROM) is presented. The first reduction, denoted as ROM, is performed via POD, taking the displacement fluctuation field as a primal variable. The second reduction, denoted as HPROM, is performed via interpolation techniques (DEIM) of the microscale stress field.
It is shown that the interpolation-based HPROM obtained in this way, leads to an ill-posed mathematical formulation when the reduction process involves an interpolant constructed using POD modes provided by the primal variable (microscale stress field). This issue has been studied in the paper, and a robust and consistent solution has been proposed.
An additional aspect in this contribution, is the selection of the interpolation points for the stress field. These interpolation points are chosen guided, not only by accuracy requirements, but also by stability considerations.
The method of selection of the interpolation points (Greedy Algorithm) is, at the present, an intensive research field. However, although in the literature there are several alternative algorithms, none of them offers a robust and general treatment to handle this purposes.
Different measures of error have been presented to test the accuracy and the convergence. The work is assessed by the homogenization of a highly complex porous metal material. The results show that, the speed-up factor is about three orders of magnitude, for an error in stress smaller than .
As conclusions of this work, it can be stated that:
In consequence, the numerical results suggest that this HPROM provides accurate solutions to problems exhibiting hardening behavior. However, some questions need to be further analyzed. For example:
These questions motivated the next research work: the development of a reduced order model applied to problems exhibiting discontinuous fields, and in particular, the case of the quasi-brittle fracture.
This work presents a novel approach to two-scale modeling of propagating fracture, based on computational homogenization FE. The specific features of this contribution are:
The approach has been validated and tested using classical benchmarks in fracture mechanics. After validation, some aspects of the proposed approach can be emphasized:
As mentioned, multiscale computational fracture problems and their extension to 3D cases, face a great challenge: the enormous involved computational cost. In consequence, next step is the development of a reduced order model aiming at diminishing the computational burden of the developed multiscale fracture model.
The reduced order model described in Paper was used as a first attempt. However, the results were very unsatisfactory. The conclusions of this interpolation-based approach to multiscale reduced order modeling in fracture cases were:
This suggests additional research and exploration of specific model order reduction techniques for multiscale fracture problems.
(1) The SVD strategy, gives importance to repeated snapshots, and mainly, snapshots which euclidean norm is considerably high.
This article proposes a set of new computational techniques to solve multiscale problems via HPROM techniques. These techniques have been applied to the multiscale model described in 4.1.3, and they are summarized next:
In a first validation stage, in order to test the sensibility of the reduction techniques, a set of three different failure cells have been tested, by increasing the complexity and, consequently, the amount of cohesive bands. A-priori and a-posteriori errors analysis are performed, showing that, increasing the complexity (number of involved operations) at the microscale, the amount of required strain and free energy modes increases only slightly for a given error. This is a clearly promising scenario.
Finally, this reduced multiscale model was also validated and tested with the L-Shape Panel test, comparing the solution with the one given by the HF (obtained with the approach described in the Paper ), and analyzing the impact on the use of different amounts of reduced order basis functions of both, the strain fluctuations and the free energy.
Several aspects of the proposed methodology can be highlighted as new contributions:
At this point it can be argued that only idealized, two-dimensional, problems have been considered. The real interest of many multiscale modeling problems residing on actual three-dimensional problems, the following question arises:
to what extent these techniques can be extended to three-dimensional problems, where the involved RVE complexity and the associated computational cost can be two or three orders of magnitudes larger than in 2D problems?
In Fig. 19, the results obtained from different kind of 2D microscale morphologies are presented. They show a very relevant property: the obtained speedup “scales” linearly with the problem-complexity. Therefore one could think of achievable values of – for the speedup in 3D problems. This fact (in conjunction with, the additional usage of HPC procedures), could turn affordable 3D multiscale fracture modelling.
|Figure 19: Speedup scalability.|
This work has been developed in combination with the reduced order model for non-smooth problems (see 4.1.4). The main objective is to develop the algorithmic procedure in a general setting to be applied to different problems involving integral operators that can be sampled. Not only problems involving multiple scales can be analyzed, but also monoscale (static and dynamic) problems based on the Finite Element method.
It is presented a general framework for the dimensional reduction in terms of numbers of degrees of freedom as well as number of integration points of nonlinear parametrized finite element models.
As in previous cases (see 4.1.2 and 4.1.4), the reduction process is divided into two sequential stages, the first consists of a Galerkin projection of the strain fluctuations, via POD, and the second consists of a novel cubature rule also used in 4.1.4. In this case, this method is deeply studied and analyzed. The distinguish features of the proposed method to be highlighted are:
This model is tested through two structural examples, (quasi-static bending, and resonant vibration of elasto-plastic composite plates). The total amount of integration points is reduced three order of magnitudes, this methodology can be applied to different primary variables, in 4.1.4, attention was focused on use the free energy to determine the reduced integration rule.
Several issues have been improved by this research: firstly, the robustness, one of the most attractive features of the proposed hyper-reduced order model (and in general, of all cubature-based ROMs) is that it preserves the spectral properties of the Jacobian matrix of the finite element motion equations. Secondly, the improved version of the Empirical cubature method, in contrast with other similar techniques proposed in the literature, in which the weights at almost all iterations of the greedy algorithm are calculated with a standard, unconstrained least-squares. In fact, the nonnegative least squares problem is included to filter out small negative weights caused by roundoff errors. And finally, for implementation purposes, the "format" of the finite element method is conserved.
This work presents a brief summary of the two-scale approach for modeling failure propagation, providing details about propagation at the macro and micro levels. This publication is centered in exploring the applicability of the method to structural problems. The four-point bending and the Nooru-Mohamed problems have been chosen as benchmarks, taking the material properties form experimental tests.
In the case of the Nooru-Mohamed test, it has been shown, the influence of the horizontal load (shear force) in the microscale behavior, and the activation of different crack patterns, representing the macroscale changes in the crack propagation scheme. In the four-point bending test, it is displayed the influence on the macroscale propagation scheme, when critical failure mechanisms at the microscale are precluded.
This work has a similar objective than the previous one. A brief summary about the reduced order model based on the two-scale approach for modeling failure propagation, has been presented. This work also presents a summary about the results obtained in the L-Shaped Panel, and the influence of the size of reduced order basis functions (for strain fluctuations and free energy) is presented and analyzed.
This work focuses on exploring different issues of the two-scale approach for modeling failure propagation. Particularly, the total energy dissipation and its relation at both scales is analyzed in some specific fracture problems.
 Oliver, J. and Dias, I.F. and Huespe, A.E. (2014) "Crack-path field and strain-injection techniques in computational modeling of propagating material failure", Volume 274. Computer Methods in Applied Mechanics and Engineering 289–348
 An, S. and Kim, T. and James, D. (2009) "Optimizing cubature for efficient integration of subspace deformations", Volume 27. Cornell University. ACM transactions on graphics 5 165
 Böhm, H. (2013) "Short introduction to basic aspects of continuum micromechanics", Volume 27. Technical Report. Institute of Lightweight Design and Structural Biomechanics (ILSB). Vienna University of Technology 206
 Eshelby, J.D. (1957) "The determination of the field of an ellipsoidal inclusion and related problems", Volume 241. Proc. R. Soc. Lond A 376–396
 Mori, T. and Tanaka, K. (1973) "Average stress in the matrix and average energy of materials with misfitting inclusions", Volume 21. Acta Metall. 571–574
 Nemat-Nasser, S. and Hori, M. (1999) "Micromechanics: overall properties of heterogeneous materials". Elsevier
 Bornert, M. (1999) "Morphologie microstructurale et comportement mécanique; charactérisations expérimentales, approches par bornes et estimations autocohérentes généralisées". PhD. Thesis, Ecole Nationale des Ponts et Chaussées
 Ponte Castañeda, P. and Suquet, P. (1998) "Nonlinear composites", Volume 34. Advances in Applied Mechanics 171–302
 Markov, K. (2000) "Elementary micromechanics of heterogeneous media", Volume. Heterogeneous Media: Micromechanics Modeling Methods and Simulations 1–162
 Milton, G. W. (2002) "The Theory of Composites". Cambridge University Press
 Torquato, S. (2002) "Random Heterogeneous Media". Springer–Verlag
 Voigt, W. (1889) "ber die Beziehung zwischen den beiden Elasticitäts-Constanten isotroper Körper", Volume 38. Ann.Phys. 573–587
 Reuss, A. (1929) "Berechnung der Fliegrenze von Mischkristallen auf Grund der Plastizitätsbe–dingung für Einkristalle", Volume 9. ZAMM 49–58
 Hill, R. (1952) "The elastic behavior of a crystalline aggregate", Volume A65. Proc. Phys. Soc. London 349–354
 Hashin, Z. and Shtrikman, S. (1963) "A variational approach to the theory of the elastic behaviour of multiphase materials", Volume 11. J. Mech. Phys. Solids 127–140
 Ferrer, A. and Oliver, J. and Cante, J. C. and Lloberas-Valls O. (2016) "Vademecum-based approach to multi-scale topological material design", Volume. Advanced Modeling and Simulation in Engineering Sciences 3–23
 Gurson, A.L. (1975) "Plastic flow and fracture behavior of ductile materials incorporating void nucleation, growth and coalescence". Brown University
 Tvergaard, V. (1981) "Influence of voids on shear band instabilities under plane strain conditions", Volume 17. Int. J. Fract. 389–407
 Giusti, S.M. and Blanco, P.J. and Neto, E.A. de Souza and Feijóo, R.A. (2009) "An assessment of the Gurson yield criterion by a computational multi-scale approach", Volume 26(3). Engineering Computations: International Journal for Computer-Aided Engineering and Software 281–301
 Hill, R. (1965) "A self-consistent mechanics of composite materials", Volume 13. J. Mech. Phys. Solids 213–222
 Mandel, J. (1971) "Plasticité Classique et Viscoplasticité". Springer–Verlag. CISM Lecture Notes
 Blanco, Pablo J. and Sánchez, Pablo J. and de Souza Neto, Eduardo A. and Feijóo, Raúl A. (2016) "Variational Foundations and Generalized Unified Theory of RVE-Based Multiscale Models", Volume 23. Archives of Computational Methods in Engineering 2 191–253
 Feyel, F. and Chaboche, J.L. (2000) "FE multiscale approach for modelling the elastoviscoplastic behaviour of long fibre SiC/Ti composite materials", Volume 183. Comput. Meth. App. Mech. Eng. 309–330
 Belytschko, T. and Loehnert, S. and Song, J.H. (2008) "Multiscale aggregating discontinuities: A method for circumventing loss of material stability", Volume 73. Int. J. Numer. Meth. Engng. 869–894
 Song, J.H. and Belytschko, T. (2009) "Multiscale aggregating discontinuities method for micro-macro failure of composites", Volume Part B 40. Composites 417–426
 Nguyen, V.P. and Lloberas-Valls, O. and Stroeven, M. and Sluys, L.J. (2010) "On the existence of representative volumes for softening quasi-brittle materials - A failure zone averaging scheme", Volume 199. Comput. Meth. App. Mech. Eng. 3028–3038
 Nguyen, V.P. and O.Lloberas-Valls and M. Stroeven, L.J. Sluys. (2010) "Homogenization-based multiscale crack modelling: from micro diffusive damage to macro cracks", Volume 200. Comput. Meth. App. Mech. Eng. 1220–1236
 Matous, K. and Kulkarni, M.G. and Geubelle, P.H. (2008) "Multiscale cohesive failure modeling of heterogeneous adhesives", Volume 56. Journal of the Mechanics and Physics of Solids 1511–1533
 Coenen, E.W.C. and Kouznetsova, V.G. and Bosco, E. and Geers, M.G.D. (2012) "A multi-scale approach to bridge microscale damage and macroscale failure: a nested computational homogenization-localization framework", Volume 178. Int. J. Fracture 1-2 157–178
 Toro, S. and Sánchez, P.J. and Huespe, A.E. and Giusti, S.M. and Blanco, P.J. and Feijóo, R.A. (2014) "A two-scale failure model for heterogeneous materials: numerical implementation based on the finite element method", Volume 97. Int. J. Num. Meth. Eng. 5 313–351
 Lloberas-Valls, O. and Rixen, D.J. and Simone, A. and Sluys, L.J. (2012) "Multiscale domain decomposition analysis of quasi-brittle heterogeneous materials", Volume 89. Int. J. Num. Meth. Eng. 11 1337–1366
 Anderson, T.L. (1995) "Fracture Mechanics, Fundamentals and Applications". 2nd, Edition, CRC press
 Pandolfi, A. and Krysl, P. and Ortiz, M. (1999) "Finite element simulation of ring expansion and fragmentation: the capturing of length and time scales through cohesive models of fracture", Volume 95. Int. J. Fracture 1-4 279–297
 Molinari, J.F. and Gazonas, G. and Raghupathy, R. and Rusinek, A. and Zhou, F. (2007) "The cohesive element approach to dynamic fragmentation: the question of energy convergence", Volume 69. Int. J. Num. Meth. Eng. 3 484–503
 Toro, S. and Sánchez, P.J. and Blanco, P.J. and de Souza Neto, E.A. and Huespe, A.E. and Feijóo, R.A. (2016) "Multiscale formulation for material failure accounting for cohesive cracks at the macro and micro scales", Volume 76. Int. J. of Plasticity 75–110
 Rodrigues, E.A. and Manzoli, O.L. and Bitencourt Jr, L.A.G. and Bittencourt, T. (2016) "2D mesoscale model for concrete based on the use of interface element with a high aspect ratio", Volume 94. International Journal of Solids and Structures 112–124
 Pijaudier-Cabot, G. and Bazant, Z. (1987) "Nonlocal damage theory", Volume 113(10). ASCE J. Engrg. Mech. 1512–1533
 Tvergaard, V. and Needleman, A. (1995) "Efects of nonlocal damage in porous plastic solids.", Volume 32(8/9). Int J. Solids Struct. 1063–1077
 Aifantis, E. (1984) "On the microestructural origin of certain inelastic models", Volume 106. ASME J. Engrg Mat. Tech. 326–334
 de Borst, R. and Mühlhaus, H. (1992) "Gradient-dependent plasticity: formulation and algorithmic aspects", Volume 35(3). Int J. Numer. Meth. Engrg. 521–539
 Peerlings, R. and de Borst, R. and Brekelmans, W. and de Vree, J. (1996) "Gradient Enhanced Damage for Quasi-Brittle Materials", Volume 39. Int J. Numer. Meth. Engrg. 3391–3403
 Peerlings, R. and Geers, M. and de Borst, R. and Brekelmans, W. (2001) "A critical comparison of nonlocal and gradient-enhanced softening continua", Volume 38. Int J. Solids Struct. 7723–7746
 Peerlings, R. and de Borst, R. and Brekelmans, W. and Geers, M. (2002) "Localization issues in local and nonlocal continuum approaches to fracture", Volume 21. European Journal of Mechanics A/Solids 175–189
 Steinmmann, P. and Willam, K. (1991) "Localization within the framework of micropolar elasto-plasticity". Advances in continuum mechanics. Springer Verlag 296–313
 Muhlhaus, H. and Vardoulakis, I. (1987) "The thickness of shear bands in granulars materials", Volume 37. Geotechnic 271–283
 Simo, J. and Oliver, J. and Armero, F. (1993) "An analysis of strong discontinuities induced by strain-softening in rate-independent inelastic solids", Volume 12. Comput. Mech. 277–296
 Oliver, J. (2000) "On the discrete constitutive models induced by strong discontinuity kinematics and continuum constitutive equations", Volume 37. Int. J. Solids Struct. 7207–7229
 Oliver, J. and Huespe, A. E. and Pulido, M. D. G. and Chaves, E. (2002) "From continuum mechanics to fracture mechanics: the strong discontinuity approach", Volume 69. Engineering Fracture Mechanics 113–136
 Oliver, J. and Huespe, A.E. (2004) "Theoretical and computational issues in modelling material failure in strong discontinuity scenarios", Volume 193. Comput. Meth. App. Mech. Eng. 2987–3014
 Oliver, J. and Huespe, A. E. (2004) "Continuum approach to material failure in strong discontinuity settings", Volume 193. Comp. Meth. Appl. Mech. in Engrg. 3195–3220
 Bazant, Z.P. and Planas, J. (1998) "Fracture and size efect in concrete and other quasibrittle materials". CRC Press
 Gitman, I.M. and Askes, H. and Sluys, L.J. (2007) "Representative volume: Existence and size determination", Volume 74. Engineering Fracture Mechanics 2518–2534
 Shankar Ganapathysubramanian and Nicholas Zabaras. (2004) "Design across length scales: a reduced-order model of polycrystal plasticity for the control of microstructure-sensitive material properties", Volume 193. Computer Methods in Applied Mechanics and Engineering 45-47 5017 - 5034
 Yvonnet, J. and He, Q-C. (2007) "The reduced model multiscale method (R3M) for the non-linear homogenization of hyperelastic media at finite strains", Volume 223. Journal of Computational Physics 1 341–368
 Boyaval, S. (2007) "Reduced-basis approach for homogenization beyond the periodic setting", Volume. Arxiv preprint math/0702674
 E. Monteiro and J. Yvonnet and Q.C. He. (2008) "Computational homogenization for nonlinear conduction in heterogeneous materials using model reduction", Volume 42. Computational Materials Science 4 704 - 712
 N.C. Nguyen. (2008) "A multiscale reduced-basis method for parametrized elliptic partial differential equations with multiple scales", Volume 227. Journal of Computational Physics 23 9807 - 9822
 Amsallem, David and Cortial, Julien and Carlberg, Kevin and Farhat, Charbel. (2009) "A method for interpolating on manifolds structural dynamics reduced-order models", Volume 80. John Wiley Sons, Ltd. International Journal for Numerical Methods in Engineering 9 1241–1258
 Maxime Barrault and Yvon Maday and Ngoc Cuong Nguyen and Anthony T. Patera. (2004) "An 'empirical interpolation' method: application to efficient reduced-basis discretization of partial differential equations", Volume 339. Comptes Rendus Mathematique 9 667 - 672
 Michel, JC and Moulinec, H and Suquet, P. (2001) "A computational scheme for linear and non-linear composites with arbitrary phase contrast", Volume 52. International Journal for Numerical Methods in Engineering 1-2 139–160
 Ryckelynck, D. (2009) "Hyper-reduction of mechanical models involving internal variables", Volume 77. International Journal for Numerical Methods in Engineering 1 75–89
 Oskay, C. and Fish, J. (2007) "Eigendeformation-based reduced order homogenization for failure analysis of heterogeneous materials", Volume 196. Comput. Meth. App. Mech. Eng. 7 1216–1243
 Zhang, S. and Oskay, C. (2016) "Reduced order variational multiscale enrichment method for elasto-viscoplastic problems", Volume 300. Computer Methods in Applied Mechanics and Engineering 199–224
 Kerfriden, P. and Goury, O. and Rabczuk, T. and Bordas, S.P. (2013) "A partitioned model order reduction approach to rationalise computational expenses in nonlinear fracture mechanics", Volume 256. Computer methods in applied mechanics and engineering 169–188
 Miehe, C. and Schotte, J. and Schröder, J. (1999) "Computational micro-macro transitions and overall moduli in the analysis of polycrystals at large strains", Volume 6. Computational Materials Science 372–382
 Kouznetsova, V.G. (2002) "Computational homogenization for the multi-scale analysis of multi-phase materials". Technische Universiteit Eindhoven
 Neto, E.A. de Souza and Feijóo, R.A. (2006) "Variational foundation on multi-scale constitutive models of solids: small and large strain kinematical formulation". LNCC Research & Development Report No 16
 Drago, A. and Pindera, M. J. (2007) "Micro-macromechanical analysis of heterogeneous materials: Macroscopically homogeneous vs periodic microstructures", Volume 67. Elsevier. Composites science and technology 6 1243–1263
 Weinan, E. (2011) "Principles of Multiscale Modeling". Princeton University, 1st Edition
 Terada, Kenjiro and Hori, Muneo and Kyoya, Takashi and Kikuchi, Noboru. (2000) "Simulation of the multi-scale convergence in computational homogenization approaches", Volume 37. International Journal of Solids and Structures 16 2285–2311
 Gross, D. and Seelig, T. (2011) "Fracture mechanics: with an introduction to micromechanics". Springer–Verlag, 2nd Edition
 Hernández, J.A. and Caicedo, M.A. and Ferrer A. (2017) "Dimensional hyper-reduction of nonlinear finite element models via empirical cubature", Volume 313. Computer Methods in Applied Mechanics and Engineering 687–722
 Hernández, J. A. and Oliver, J. and Huespe, A. and Caicedo, M. (2012) "High-performance model reduction procedures in multiscale simulations". CIMNE, 1nd Edition
 Hughes, T.J.R. (1987) "The Finite Element Method. Linear Static and Dynamic Finite Element Analysis.". Prentice-Hall
 Zienkiewicz, O.C. and Taylor, R.L. (2000) "The Finite Element Method". Butterworth-Heinemann
 Belytshchko, T. and Liu, W. K. and Moran, B. (2001) "Nonlinear Finite Elements for Continua and Structures". J. Wiley & Sons Ltd.
 Bazant, Z. (2010) "Can Multiscale-Multiphysics Methods Predict Softening Damage and Structural Failure?", Volume 8. International Journal for Multiscale Computational Engineering 1 61–67
 Rots, J. (1988) "Computational modelling of concrete fracture". Delft University
 Oliver, J. (1989) "A consistent characteristic length for smeared cracking models", Volume 28. John Wiley Sons, Ltd. International Journal for Numerical Methods in Engineering 2 461–474
 Oliver, J. and Caicedo, M. and Roubin, E. and Huespe, A.E. and HernÃ¡ndez, J.A. (2015) "Continuum approach to computational multiscale modeling of propagating fracture", Volume 294. Computer Methods in Applied Mechanics and Engineering 384–427
 Ghosh, Somnath and Lee, Kyunghoon and Moorthy, Suresh. (1996) "Two scale analysis of heterogeneous elastic-plastic materials with asymptotic homogenization and Voronoi cell finite element model", Volume 132. Computer Methods in Applied Mechanics and Engineering 1 63–116
 Moulinec, H. and Suquet, P. (1998) "A numerical method for computing the overall response of nonlinear composites with complex microstructure", Volume 157. Computer Methods in Applied Mechanics and Engineering 1 69–94
 Michel, J.C. and Moulinec, H. and Suquet, P. (1999) "Effective properties of composite materials with periodic microstructure: a computational approach", Volume 172. Comput. Meth. App. Mech. Eng. 109–143
 Oliver, J. (1995b) "Continuum modelling of strong discontinuities in solids mechanics using damage models", Volume 17(1-2). Comput. Mech. 49–61
 Faria, R. and Oliver, J. and Cervera, M. (1998) "A strain-based plastic viscous-damage model for massive concrete structures", Volume 35. International Journal of Solids and Structures 14 1533–1558
 Oliver, J. and Huespe, A.E. and Blanco, S. and Linero, D.L. (2005) "Stability and robustness issues in numerical modeling of material failure with the strong discontinuity approach", Volume 195. Comput. Meth. App. Mech. Eng. 52 7093–7114
 Oliver, J. (1996a) "Modelling strong discontinuities in solids mechanics via strain softening constitutive equations. Part 1: Fundamentals", Volume 39(21). Int. j. numer. methods eng. 3575–3600
 Geers, M.G.D. and Kouznetsova, V.G. and Brekelmans, W.A.M. (2010) "Multi-scale computational homogenization: Trends and challenges", Volume 234. Journal of Computational and Applied Mathematics 2175–2182
 Otero, F. and Martínez, X. and Oller, S. and Salomón, O. (2015) "An efficient multi-scale method for non-linear analysis of composite structures", Volume 131. Composite Structures 707–719
 Lesicar, T. and Tonkovi, Z. and Sori, J. (2015) "Multi-Scale Modeling of Heterogeneous Materials and the Validation Challenge", Volume 665. Key Engineering Materials 181–184
 Oliver, J. and Caicedo, M. and Huespe, A.E. and HernÃ¡ndez, J.A. and Roubin, E. (2017) "Reduced order modeling strategies for computational multiscale fracture", Volume 313. Computer Methods in Applied Mechanics and Engineering 560–595
 Michel, J.C. and Suquet, P. (2003) "Nonuniform transformation field analysis", Volume 40. International Journal of Solids and Structures 25 6937–6955
 Michel, J.C. and Suquet, P. (2004) "Computational analysis of nonlinear composite structures using the nonuniform transformation field analysis", Volume 193. Computer Methods in Applied Mechanics and Engineering 48–51 5477–5502
 Unger, J. F. (2013) "An FE2-X1 approach for multiscale localization phenomena", Volume 61. Journal of the Mechanics and Physics of Solids 4 928–948
 Simo, J. and Hughes, T. J. R. (1986) "On the variational foundations of assumed strain methods", Volume 53. Journal of Applied Mechanics, ASME 51–54
 Simo, J. and Riffai, M. (1990) "A class of mixed assumed strain methods and the method of incompatible modes", Volume 29. Int. J. Numer. Methods eng. 1595–1638
 Reddy, B. D. and Simo, J. (1995) "Stability and convergence of a class of enhanced strain methods", Volume 32. SIAM, Jornal on Numerical Analysis 1705–1728
 Oliver, J. and Huespe, A. E. and Cante, J. C. and DÃaz, G. (2010) "On the numerical resolution of the discontinuous material bifurcation problem", Volume 83. John Wiley Sons, Ltd. International Journal for Numerical Methods in Engineering 6 786–804
 Hernández, J.A. and Oliver, J. and Huespe, A.E. and Caicedo, M.A. and Cante, J.C. (2014) "High-performance model reduction techniques in computational multiscale homogenization", Volume 276. Computer Methods in Applied Mechanics and Engineering 149–189
 Hernández, J.A. and Oliver, J. and Huespe, A.E. and Caicedo, M.A. and Cante, J.C. (2014) "Computational Homogenization of Inelastic Materials using Model Order Reduction", Volume 141. International Center for Numerical Methods in Engineering (CIMNE) Monograph
 Simo, J. and Hughes, T. (1998) "Computational inelasticity". Springer-Verlag
 Oliver, J. and Caicedo, M. and Huespe, A.E. and Hernández, J.A. and Roubin, E. (2017) "Reduced order modeling strategies for computational multiscale fracture", Volume 313. Computer Methods in Applied Mechanics and Engineering 560–595
 Hoffman, D. (2014) "Numerical Methods for Engineers and Scientists", Volume 2nd. Edition. Marcel Deckker, Inc.
 DeVore, R. and Iserles, A. and Suli, E. (2001) "Foundations of Computational Mathematics", Volume LMS 284. Cambridge University Press
 Ryckelynck, D. (2005) "A priori hyperreduction method: an adaptive approach", Volume 202. Journal of Computational Physics 1 346–366
 Oden, J. T. and Belytschko, T. Fish J. and Hughes, T. J. R. and Johnson, C. and Keyes, L. A and Petzold, L and Srolovitz, L. and Yip, S. (2006) "Simulation-based engineering science". National Science Foundation (USA)
 Farhat, C. and Chapman, T. and Avery, P. (2015) "Structure-preserving, stability, and accuracy properties of the energy-conserving sampling and weighting method for the hyper reduction of nonlinear finite element dynamic models". International Journal for Numerical Methods in Engineering
 A. Arteiro and G. Catalanotti and A.R. Melro and P. Linde and P.P. Camanho. (2014) "Micro-mechanical analysis of the in situ effect in polymer composite laminates", Volume 116. Composite Structures 827–840
 Belytschko, T. and Chen, H. and Xu, J.X. and Zi, G. (2003) "Dynamic crack propagation based on loss of hyperbolicity and a new discontinuous enrichment", Volume 58. Int. J. Numer. Methods Engrg. 1873–1905
 Prabel, B. and Combescure, A. and Gravouil, A. and Marie, S. (2007) "Level set X-FEM non-matching meshes: application to dynamic crack propagation in elastic-plastic media", Volume 69. John Wiley Sons, Ltd. International Journal for Numerical Methods in Engineering 8 1553–1569
 C. Linder and F. Armero. (2009) "Finite elements with embedded branching", Volume 45. Finite Elements in Analysis and Design 4 280-293
 O. Lloberas-Valls and A.E. Huespe and J. Oliver and I.F. Dias. (2016) "Strain injection techniques in dynamic fracture modeling", Volume 308. Computer Methods in Applied Mechanics and Engineering 499–534
Title: High-performance model reduction techniques in computational multiscale homogenization.
Journal of Computer Methods in Applied Mechanics and Engineering
Editors: Thomas J.R. Hughes, J. Tinsley Oden, Manolis Papadrakakis
Title: Continuum approach to computational multiscale modeling of propagting fracture.
Journal of Computer Methods in Applied Mechanics and Engineering
Editors: Thomas J.R. Hughes, J. Tinsley Oden, Manolis Papadrakakis
Title: Reduced Order Modeling strategies for Computational Multiscale Fracture.
Journal of Computer Methods in Applied Mechanics and Engineering
Editors: Thomas J.R. Hughes, J. Tinsley Oden, Manolis Papadrakakis
Title: Dimensional hyper-reduction of nonlinear finite element models via empirical cubature.
Journal of Computer Methods in Applied Mechanics and Engineering
Editors: Thomas J.R. Hughes, J. Tinsley Oden, Manolis Papadrakakis
Title: Continuum Approach to Computational Multi-Scale Modeling of Fracture.
Key Engineering Materials Vol. 627
Advances in Fracture and Damage Mechanics XIII
Editors: J. Alfaiate and M.H. Aliabadi
Trans Tech Publications
Title: Model Order Reduction in computational multiscale fracture mechanics.
Key Engineering Materials Vol. 713
Advances in Fracture and Damage Mechanics XV
Editors: Jesús Toribio, Vladislav Mantic, Andrés Sáez, M.H. Ferri Aliabadi
Trans Tech Publications
Title: Multi-scale (FE) analysis of material failure in cement/aggregate-type composite structures
Computational Modelling of Concrete Structures
Proceedings of EURO–C 2014
Editors: N. Biani; H Mang; Gunther Meschke; Reneé de Borst