Source author record

Pierre Kerfriden

Pierre Kerfriden appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

A Bayesian multiscale CNN framework to predict local stress fields in structures with microscale features

Multiscale computational modelling is challenging due to the high computational cost of direct numerical simulation by finite elements. To address this issue, concurrent multiscale methods use the solution of cheaper macroscale surrogates as boundary conditions to microscale sliding windows. The microscale problems remain a numerically challenging operation both in terms of implementation and cost. In this work we propose to replace the local microscale solution by an Encoder-Decoder Convolutional Neural Network that will generate fine-scale stress corrections to coarse predictions around unresolved microscale features, without prior parametrisation of local microscale problems. We deploy a Bayesian approach providing credible intervals to evaluate the uncertainty of the predictions, which is then used to investigate the merits of a selective learning framework. We will demonstrate the capability of the approach to predict equivalent stress fields in porous structures using linearised and finite strain elasticity theories.

preprint2022arXiv

Concurrent multiscale analysis without meshing: Microscale representation with CutFEM and micro/macro model blending

In this paper, we develop a novel unfitted multiscale framework that combines two separate scales represented by only one single computational mesh. Our framework relies on a mixed zooming technique where we zoom at regions of interest to capture microscale properties and then mix the micro and macroscale properties in a transition region. Furthermore, we use homogenization techniques to derive macro model material properties. The microscale features are discretized using CutFEM. The transition region between the micro and macroscale is represented by a smooth blending function. To address the issues with ill-conditioning of the multiscale system matrix due to the arbitrary intersections in cut elements and the transition region, we add stabilization terms acting on the jumps of the normal gradient (ghost-penalty stabilization). We show that our multiscale framework is stable and is capable to reproduce mechanical responses for heterogeneous structures in a mesh-independent manner. The efficiency of our methodology is exemplified by 2D and 3D numerical simulations of linear elasticity problems.

preprint2014arXiv

Certification of projection-based reduced order modelling in computational homogenisation by the Constitutive Relation Error

In this paper, we propose upper and lower error bounding techniques for reduced order modelling applied to the computational homogenisation of random composites. The upper bound relies on the construction of a reduced model for the stress field. Upon ensuring that the reduced stress satisfies the equilibrium in the finite element sense, the desired bounding property is obtained. The lower bound is obtained by defining a hierarchical enriched reduced model for the displacement. We show that the sharpness of both error estimates can be seamlessly controlled by adapting the parameters of the corresponding reduced order model.

preprint2013arXiv

A goal-oriented reduced basis method for the wave equation in inverse analysis

In this paper, we extend the reduced-basis methods developed earlier for wave equations to goal-oriented wave equations with affine parameter dependence. The essential new ingredient is the dual (or adjoint) problem and the use of its solution in a sampling procedure to pick up "goal-orientedly" parameter samples. First, we introduce the reduced-basis recipe --- Galerkin projection onto a space $Y_N$ spanned by the reduced basis functions which are constructed from the solutions of the governing partial differential equation at several selected points in parameter space. Second, we propose a new "goal-oriented" Proper Orthogonal Decomposition (POD)--Greedy sampling procedure to construct these associated basis functions. Third, based on the assumption of affine parameter dependence, we use the offline-online computational procedures developed earlier to split the computational procedure into offline and online stages. We verify the proposed computational procedure by applying it to a three-dimensional simulation dental implant problem. The good numerical results show that our proposed procedure performs better than the standard POD--Greedy procedure in terms of the accuracy of output functionals.

preprint2013arXiv

Isogeometric cohesive elements for two and three dimensional composite delamination analysis

Isogeometric cohesive elements are presented for modeling two and three dimensional delaminated composite structures. We exploit the knot insertion algorithm offered by NURBS (Non Uniform Rational B-splines) to generate cohesive elements along delamination planes in an automatic fashion. A complete computational framework is presented including pre-processing, processing and post-processing. They are explained in details and implemented in MIGFEM--an open source Matlab Isogemetric Analysis code developed by the authors. The composite laminates are modeled using both NURBS solid and shell elements. Several two and three dimensional examples ranging from standard delamination tests (the mixed mode bending test), the L-shaped specimen with a fillet, three dimensional (3D) double cantilever beam and a 3D singly curved thick-walled laminate are provided. To the authors' knowledge, it is the first time that NURBS-based isogeometric analysis for two/three dimensional delamination modeling is presented. For all examples considered, the proposed framework outperforms conventional Lagrange finite elements.

preprint2013arXiv

Locally equilibrated stress recovery for goal oriented error estimation in the extended finite element method

Goal oriented error estimation and adaptive procedures are essential for the accurate and efficient evaluation of numerical simulations that involve complex domains. By locally improving the approximation quality we can solve expensive problems which could result intractable otherwise. Here, we present an error estimation technique for enriched finite element approximations that is based on an equilibrated recovery technique, which considers the stress intensity factor as the quantity of interest. The locally equilibrated superconvergent patch recovery is used to obtain enhanced stress fields for the primal and dual problems defined to evaluate the error estimate.

preprint2013arXiv

Mesh adaptivity driven by goal-oriented locally equilibrated superconvergent patch recovery

During the last decade there has been an increase on the use of goal-oriented error estimates which help to quantify and control the local error on a quantity of interest (QoI) that might result relevant for design purposes (e.g. the mean stress or mean displacement in a particular area, the stress intensity factor for fracture problems,...). Residual-based error estimators have been used to estimate the error in quantities of interest for finite element approximations. This work presents a recovery-based error estimation technique for QoI whose main characteristic is the use of an enhanced version of the Superconvergent Patch Recovery (SPR) technique developed by Zienkiewicz and Zhu. This enhanced version of the SPR technique, used to recover the primal and dual solutions, provides a nearly statically admissible stress field that results in accurate estimations of the error in the QoI.

preprint2013arXiv

Nitsche's method for two and three dimensional NURBS patch coupling

A Nitche's method is presented to couple non-conforming two and three dimensional NURBS (Non Uniform Rational B-splines) patches in the context of isogeometric analysis (IGA). We present results for elastic stress analyses under the static condition of two and three dimensional NURBS geometries. The contribution fills the gap in the literature and enlarges the applicability of NURBS-based isogeometric analysis.

preprint2013arXiv

Space-time goal-oriented reduced basis approximation for linear wave equation

In this paper, we study numerically the linear damped second-order hyperbolic partial differential equation (PDE) with affine parameter dependence using a goal-oriented approach by finite element (FE) and reduced basis (RB) methods. The main contribution of this paper is the "goal-oriented" proper orthogonal decomposition (POD)-Greedy sampling procedure within the RB approximation context. First, we introduce the RB recipe: Galerkin projection onto a space $Y_N$ spanned by solutions of the governing PDE at $N$ selected points in parameter space. This set of $N$ parameter points is constructed by the standard POD-Greedy sampling procedure already developed. Second, based on the affine parameter dependence, we make use of the offline-online computational procedures: in the offline stage, we generate the RB space; in the online stage, given a new parameter value, we calculate rapidly and accurately the space-time RB output of interest and its associated asymptotic error. The proposed goal-oriented POD-Greedy sampling procedure can now be implemented and will look for the parameter points such that it minimizes this (asymptotic) output error rather than the solution error (or, error indicator which is the dual norm of residual) as in the standard POD-Greedy procedure. Numerical results show that the new goal-oriented POD-Greedy sampling procedure improves significantly the accuracy of the space-time output computation in comparison with the standard POD-Greedy one. The method is thus ideally suited for repeated, rapid and reliable evaluation of input-output relationships within the space-time setting.

preprint2013arXiv

Virtual Delamination Testing through Non-Linear Multi-Scale Computational Methods: Some Recent Progress

This paper deals with the parallel simulation of delamination problems at the meso-scale by means of multi-scale methods, the aim being the Virtual Delamination Testing of Composite parts. In the non-linear context, Domain Decomposition Methods are mainly used as a solver for the tangent problem to be solved at each iteration of a Newton-Raphson algorithm. In case of strongly nonlinear and heterogeneous problems, this procedure may lead to severe difficulties. The paper focuses on methods to circumvent these problems, which can now be expressed using a relatively general framework, even though the different ingredients of the strategy have emerged separately. We rely here on the micro-macro framework proposed in (Ladevèze, Loiseau, and Dureisseix, 2001). The method proposed in this paper introduces three additional features: (i) the adaptation of the macro-basis to situations where classical homogenization does not provide a good preconditioner, (ii) the use of non-linear relocalization to decrease the number of global problems to be solved in the case of unevenly distributed non-linearities, (iii) the adaptation of the approximation of the local Schur complement which governs the convergence of the proposed iterative technique. Computations of delamination and delamination-buckling interaction with contact on potentially large delaminated areas are used to illustrate those aspects.

preprint2012arXiv

A partitioned model order reduction approach to rationalise computational expenses in multiscale fracture mechanics

We propose in this paper an adaptive reduced order modelling technique based on domain partitioning for parametric problems of fracture. We show that coupling domain decomposition and projection-based model order reduction permits to focus the numerical effort where it is most needed: around the zones where damage propagates. No \textit{a priori} knowledge of the damage pattern is required, the extraction of the corresponding spatial regions being based solely on algebra. The efficiency of the proposed approach is demonstrated numerically with an example relevant to engineering fracture.

preprint2012arXiv

Statistical extraction of process zones and representative subspaces in fracture of random composite

We propose to identify process zones in heterogeneous materials by tailored statistical tools. The process zone is redefined as the part of the structure where the random process cannot be correctly approximated in a low-dimensional deterministic space. Such a low-dimensional space is obtained by a spectral analysis performed on pre-computed solution samples. A greedy algorithm is proposed to identify both process zone and low-dimensional representative subspace for the solution in the complementary region. In addition to the novelty of the tools proposed in this paper for the analysis of localised phenomena, we show that the reduced space generated by the method is a valid basis for the construction of a reduced order model.

preprint2011arXiv

A relocalization technique for the multiscale computation of delamination in composite structures

In this article, we present numerical enhancements of a multiscale domain decomposition strategy based on a LaTIn solver and dedicated to the computation of the debounding in laminates. We show that the classical scale separation is irrelevant in the process zones, which results in a drop in the convergence rate of the strategy. The scalability is restored by performing nonlinear subresolutions in the vincinity of the front of the crack, at each prediction stage of the iterative solver.

preprint2011arXiv

A three-scale domain decomposition method for the 3D analysis of debonding in laminates

The prediction of the quasi-static response of industrial laminate structures requires to use fine descriptions of the material, especially when debonding is involved. Even when modeled at the mesoscale, the computation of these structures results in very large numerical problems. In this paper, the exact mesoscale solution is sought using parallel iterative solvers. The LaTIn-based mixed domain decomposition method makes it very easy to handle the complex description of the structure; moreover the provided multiscale features enable us to deal with numerical difficulties at their natural scale; we present the various enhancements we developed to ensure the scalability of the method. An extension of the method designed to handle instabilities is also presented.

preprint2011arXiv

Bridging Proper Orthogonal Decomposition methods and augmented Newton-Krylov algorithms: an adaptive model order reduction for highly nonlinear mechanical problems

This article describes a bridge between POD-based model order reduction techniques and the classical Newton/Krylov solvers. This bridge is used to derive an efficient algorithm to correct, "on-the-fly", the reduced order modelling of highly nonlinear problems undergoing strong topological changes. Damage initiation problems are addressed and tackle via a corrected hyperreduction method. It is shown that the relevancy of reduced order model can be significantly improved with reasonable additional costs when using this algorithm, even when strong topological changes are involved.

preprint2011arXiv

Improved multiscale computational strategies for delamination

This paper presents a three-scale computational strategy for the study of composite modeled at the mesoscale so that delamination can be reliably simulated. The solver is based on a LaTIn approach so that nonlinearities can be tackled at the local scale. We show how search directions which are parameters of the method need to be updated according to the delamination status of interfaces. We also present the parallelization of the macro problem which transmits long-range phenomena and warranty the scalabilty of the method; this parallelization is based on the balancing domain decomposition method, it leads to the introduction of a third (super-macro) scale.

preprint2011arXiv

Local/global model order reduction strategy for the simulation of quasi-brittle fracture

This paper proposes a novel technique to reduce the computational burden associated with the simulation of localised failure. The proposed methodology affords the simulation of damage initiation and propagation whilst concentrating the computational effort where it is most needed, i.e. in the localisation zones. To do so, a local/global technique is devised where the global (slave) problem (far from the zones undergoing severe damage and cracking) is solved for in a reduced space computed by the classical Proper Orthogonal Decomposition, while the local (master) degrees of freedom (associated with the part of the structure where most of the damage is taking place) are fully resolved. Both domains are coupled through a local/global technique. This method circumvents the difficulties associated with model order reduction for the simulation of highly non-linear mechanical failure and offers an alternative or complementary approach to the development of multiscale fracture simulators.

preprint2011arXiv

On the control of the load increments for a proper description of multiple delamination in a domain decomposition framework

In quasi-static nonlinear time-dependent analysis, the choice of the time discretization is a complex issue. The most basic strategy consists in determining a value of the load increment that ensures the convergence of the solution with respect to time on the base of preliminary simulations. In more advanced applications, the load increments can be controlled for instance by prescribing the number of iterations of the nonlinear resolution procedure, or by using an arc-length algorithm. These techniques usually introduce a parameter whose correct value is not easy to obtain. In this paper, an alternative procedure is proposed. It is based on the continuous control of the residual of the reference problem over time, whose measure is easy to interpret. This idea is applied in the framework of a multiscale domain decomposition strategy in order to perform 3D delamination analysis.

preprint2011arXiv

On the role of enrichment and statical admissibility of recovered fields in a-posteriori error estimation for enriched finite element methods

Purpose: This paper aims at assessing the effect of (1) the statical admissibility of the recovered solution; (2) the ability of the recovered solution to represent the singular solution; on the accuracy, local and global effectivity of recovery-based error estimators for enriched finite element methods (e.g. the extended finite element method, XFEM). Design/methodology/approach: We study the performance of two recovery techniques. The first is a recently developed superconvergent patch recovery procedure with equilibration and enrichment (SPR-CX). The second is known as the extended moving least squares recovery (XMLS), which enriches the recovered solutions but does not enforce equilibrium constraints. Both are extended recovery techniques as the polynomial basis used in the recovery process is enriched with singular terms for a better description of the singular nature of the solution. Findings: Numerical results comparing the convergence and the effectivity index of both techniques with those obtained without the enrichment enhancement clearly show the need for the use of extended recovery techniques in Zienkiewicz-Zhu type error estimators for this class of problems. The results also reveal significant improvements in the effectivities yielded by statically admissible recovered solutions. Originality/value: This work shows that both extended recovery procedures and statical admissibility are key to an accurate assessment of the quality of enriched finite element approximations.