Browse ORBi

- What it is and what it isn't
- Green Road / Gold Road?
- Ready to Publish. Now What?
- How can I support the OA movement?
- Where can I learn more?

ORBi

Isogeometric boundary element methods for three dimensional fatigue crack growth ; ; et al Report (2015) The isogeometric boundary element method (IGABEM) based on NURBS is adopted to model fracture problem in 3D. The NURBS basis functions are used in both crack representation and physical quantity ... [more ▼] The isogeometric boundary element method (IGABEM) based on NURBS is adopted to model fracture problem in 3D. The NURBS basis functions are used in both crack representation and physical quantity approximation. A stable quadrature scheme for singular integration is proposed to enhance the robustness of the method in dealing with highly distorted element. The convergence study in crack opening displacement is performed for penny-shaped crack and elliptical crack. Two ways to extract stress intensity factors (SIFs), the contour $M$ integral and virtual crack closure integral, are implemented based on the framework of dual integral equations. An algorithm is outlined and validated to be stable for fatigue crack growth, thanks to the smoothness not only in crack geometry but also in stress/SIFs solution brought by IGABEM. [less ▲] Detailed reference viewed: 536 (14 UL)Scale selection in nonlinear fracture mechanics of heterogeneous materials ; ; Bordas, Stéphane in Philosophical Magazine (2015), 95(28-30), 3328-3347 A new adaptive multiscale method for the non-linear fracture simulation of heterogeneous materials is proposed. The two major sources of errors in the finite element simulation are discretization and ... [more ▼] A new adaptive multiscale method for the non-linear fracture simulation of heterogeneous materials is proposed. The two major sources of errors in the finite element simulation are discretization and modelling errors. In the failure problems, the discretization error increases due to the strain localization which is also a source for the error in the homogenization of the underlying microstructure. In this paper, the discretization error is controlled by an adaptive mesh refinement procedure following the Zienkiewicz–Zhu technique, and the modelling error, which is the resultant of homogenization of microstructure, is controlled by replacing the macroscopic model with the underlying heterogeneous microstructure. The scale adaptation criterion which is based on an error indicator for homogenization is employed for our non-linear fracture problem. The control of both discretization and homogenization errors is the main feature of the proposed multiscale method. [less ▲] Detailed reference viewed: 204 (12 UL)Advances in error estimation for homogenisation ; ; et al in 13th U.S. National Congress on Computational Mechanics (2015, July) In this paper, the concept of modeling error is extended to the homogenisation of elliptic PDEs. The main difficulty is the lack of a full description of the diffusion coefficients. We overcome this ... [more ▼] In this paper, the concept of modeling error is extended to the homogenisation of elliptic PDEs. The main difficulty is the lack of a full description of the diffusion coefficients. We overcome this obstacle by representing them as a random a field. Under this framework, it is possible to quantify the accuracy of the surrogate model (the homogenised model) in terms of first moments of the energy norm and quantities of interest. This work builds on the seminal work of [1]. The methodology here presented rely on the Constitutive Relation Error (CRE) which states that a certain measures of the primal and dual surrogate model upper bound the exact error. The surrogate model, in agreement with homogenisation, is deterministic. This property exploited to obtain bounds whose computation is also deterministic. It is also shown that minimising the CRE in the set of homogenisation schemes leads us to an optimal surrogate that is closely related to the classical Voigt and Reuss models. Numerical examples demonstrate that the bounds are easy and affordable to compute, and useful as long as the mismatch between he diffusion coefficients of the microstructure remain small. In the case of high mismatch, extensions are proposed, through the introduction of stochastic surrogate models.. [1]Romkes, Albert, J. Tinsley Oden, and Kumar Vemaganti."Multi-scale goal-oriented adaptive modeling of random heterogeneous materials." Mechanics of materials 38.8(2006): 859-872. [less ▲] Detailed reference viewed: 120 (5 UL)Quasicontinuum methods for planar beam lattices (abstract) Beex, Lars ; ; et al Scientific Conference (2015, July) Detailed reference viewed: 312 (1 UL)Multi-scale fracture, model reduction, CAD and image as a model Bordas, Stéphane ; ; Beex, Lars et al Scientific Conference (2015, June 24) Detailed reference viewed: 304 (7 UL)The stable GFEM. Convergence, accuracy and Diffpack implementation ; ; Bordas, Stéphane et al Presentation (2015, May 12) Detailed reference viewed: 95 (1 UL)Error estimation in homogenisation ; ; et al Presentation (2015, January 30) Detailed reference viewed: 128 (3 UL)Adaptive methods for multiscale fracture Bordas, Stéphane ; ; et al in International Journal of Engineering Science (2015, January 01) Adaptive methods for multiscale fracture In this work, we discuss two classes of methods to reduce the complexity of (multi scale) fracture simulations. In a first part, we discuss algebraic model ... [more ▼] Adaptive methods for multiscale fracture In this work, we discuss two classes of methods to reduce the complexity of (multi scale) fracture simulations. In a first part, we discuss algebraic model reduction. We show that algebraic model reduction such as the proper orthogonal decomposition cannot be used directly because of the lack of corelation introduced by the damage or cracks. We demonstrate the use of proper orthogonal decompositions by subdomains as a candidate to reduce computational expenses in non-linear fracture simulations whilst controlling the error level. We then consider algebraic model reduction, namely the proper orthogonal decomposition(POD) to drastically reduce the computational time associated with computing the response of representative volume elements (RVEs) used in homogenization, e.g. by the FE2 method. The snapshots are obtained by solving the RVE boundary value problem for various loading paths. To speed-up the computations, system approximation through the discrete empirical interpolation (DEIM) is used and allows the evaluation of the internal forces for only a small subset of the elements making the RVE structure. In a second part, we propose an adaptive hybrid multiscale method for modelling fracture in a heterogeneous material that is composed of orthotropic grains with cohesive interfaces between the grains. Instead of a direct solver, FE2 method [1] based on homogenisation is employed in order to compute the effective behaviour of the heterogeneous microscopic material on the coarser scale. At this scale the modelling error due to the homogenisation is still low [3]. The coarse scale is discretized with unstructured triangular finite elements, and adaptive mesh refinement is used to control the discretizsation error. While the mesh refinement keeps the discretisation error with in a certain range, the modelling error increases due to the fact that by refining the coarse elements, the scale separation assumption which is a key issue for homogenisation may no longer be fulfilled [4]. Whereas the modelling error is inversely proportional to the size of the coarse elements, a critical element size can be found that corresponds to the critical value of the modelling error. A critical zone emerges when the size of a coarse element reaches the critical size, or if the underlying representative volume element of the microstructure loses stability due to localisation (lack of scale separation). Thereafter, a zoom-in process is triggered that replaces the corresponding coarse elements of the critical zone with high resolution microscale mesh to which it glues the coarse scale mesh through a strong coupling technique using Lagrange multipliers [5]. The high resolution region can gradually be extended to include the newly emerging critical zones. A local arc-length technique is adopted to trace the highly non-linear curve of the global load-displacement by controlling the opening of microscopic cohesive cracks in the fully resolved regions. The proposed adaptive multiscale method allows us to introduce progressive discrete micro cracks at the macroscale. The unstructured mesh enables us to model problems with non-regular shapes, and the arc-length method, defined over multiple scales, allows the regularisation of softening problems that are treated in quasi-statics. We exercise this method on the simulation of polycrystalline fracture, where each grain is considered orthotropic and compare results to direct numerical simulation. [less ▲] Detailed reference viewed: 664 (22 UL)Interfacial shear stress optimization in sandwich beams with polymeric core using non-uniform distribution of reinforcing ingredients ; ; et al in Composite Structures (2015) Core shearing and core/face debonding are two common failure states of sandwich beams which are mainly the result of excessive shear stresses in the core. Generally, the core made of homogeneous Fiber ... [more ▼] Core shearing and core/face debonding are two common failure states of sandwich beams which are mainly the result of excessive shear stresses in the core. Generally, the core made of homogeneous Fiber Reinforced Polymer (FRP) shows better shear resistance in comparison with that made of pure polymer. Usually, this enhancement is however somewhat limited. This paper proposes a methodology to decrease interfacial stresses by presenting the optimal distribution of reinforcing ingredients in the polymeric matrix. For this purpose, a Non-Uniform Rational Bspline (NURBS) based reinforcement distribution optimizer is developed. This technique aims at the local stress minimization within any arbitrary zone of the design domain. In our methodology, optimization and model analysis (calculation of the objective function and the design constraints) have common data sets. The quadratic NURBS basis functions smoothly define the reinforcement distribution function as a NURBS surface. The core and face sheets are modeled as multi-patches and compatibility in the displacement field is enforced by the penalty method. An adjoint sensitivity method is devised to minimize the objective function within areas of interest defined over arbitrary regions in the design domain. It is also used for efficient updating of design variables through optimization iterations. The method is verified by several examples. [less ▲] Detailed reference viewed: 225 (12 UL)A fast, certified and "tuning-free" two-field reduced basis method for the metamodelling of parametrised elasticity problems ; ; Bordas, Stéphane in Computer Methods in Applied Mechanics and Engineering (2015) This paper proposes a new reduced basis algorithm for the metamodelling of parametrised elliptic problems. The developments rely on the Constitutive Relation Error (CRE), and the construction of separate ... [more ▼] This paper proposes a new reduced basis algorithm for the metamodelling of parametrised elliptic problems. The developments rely on the Constitutive Relation Error (CRE), and the construction of separate reduced order models for the primal variable (displacement) and flux (stress) fields. A two-field Greedy sampling strategy is proposed to construct these two fields simultaneously and efficient manner: at each iteration, one of the two fields is enriched by increasing the dimension of its reduced space in such a way that the CRE is minimised. This sampling strategy is then used as a basis to construct goal-oriented reduced order modelling. The resulting algorithm is certified and "tuning-free": the only requirement from the engineer is the level of accuracy that is desired for each of the outputs of the surrogate. It is also one order of magnitude more efficient in terms of computational expenses than competing methodologies. [less ▲] Detailed reference viewed: 398 (12 UL)Gradient Smoothing in Finite Elasticity: near-incompressibility ; Bordas, Stéphane ; Report (2015) Detailed reference viewed: 142 (10 UL)Locally equilibrated stress recovery for goal oriented error estimation in the extended finite element method Bordas, Stéphane ; ; et al in Computers and Structures (2015) Goal oriented error estimation and adaptive procedures are essential for the accurate and efficient evaluation of finite element numerical simulations that involve complex domains. By locally improving ... [more ▼] Goal oriented error estimation and adaptive procedures are essential for the accurate and efficient evaluation of finite element numerical simulations that involve complex domains. By locally improving the approximation qual- ity, for example, by using the extended finite element method (XFEM), we can solve expensive problems which could result intractable otherwise. Here, we present an error estimation technique for enriched finite element approxi- mations 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. [less ▲] Detailed reference viewed: 209 (11 UL)Probabilistic multiconstraints optimization of cooling channels in ceramic matrix composites ; ; Bordas, Stéphane et al in Composites. Part B, Engineering (2015), 81 This paper presents a computational reliable optimization approach for internal cooling channels in Ceramic Matrix Composite (CMC) under thermal and mechanical loadings. The algorithm finds the optimal ... [more ▼] This paper presents a computational reliable optimization approach for internal cooling channels in Ceramic Matrix Composite (CMC) under thermal and mechanical loadings. The algorithm finds the optimal cooling capacity of all channels (which directly minimizes the amount of coolant needed). In the first step, available uncertainties in the constituent material properties, the applied mechanical load, the heat flux and the heat convection coefficient are considered. Using the Reliability Based Design Optimization (RBDO) approach, the probabilistic constraints ensure the failure due to excessive temperature and deflection will not happen. The deterministic constraints restrict the capacity of any arbitrary cooling channel between two extreme limits. A “series system” reliability concept is adopted as a union of mechanical and thermal failure subsets. Having the results of the first step for CMC with uniformly distributed carbon (C-) fibers, the algorithm presents the optimal layout for distribution of the C-fibers inside the ceramic matrix in order to enhance the target reliability of the component. A sequential approach and B-spline finite elements have overcome the cumbersome computational burden. Numerical results demonstrate that if the mechanical loading dominates the thermal loading, C-fibers distribution can play a considerable role towards increasing the reliability of the design. [less ▲] Detailed reference viewed: 129 (4 UL)Shape optimisation with isogeometric boundary element methods ; Bordas, Stéphane ; Presentation (2014, December) Detailed reference viewed: 135 (3 UL)Multiscale computational mechanics: industrial applications Bordas, Stéphane ; ; Beex, Lars et al Presentation (2014, November 25) Detailed reference viewed: 190 (7 UL)Multiscale fracture across scales and time Bordas, Stéphane ; Scientific Conference (2014, November 11) Multi-scale Computational Mechanics in Aerospace Engineering Flying is today one of the safest ways to spend our time. In the United Kingdom, for example, it is 33,000 times more likely to die from a ... [more ▼] Multi-scale Computational Mechanics in Aerospace Engineering Flying is today one of the safest ways to spend our time. In the United Kingdom, for example, it is 33,000 times more likely to die from a clinical error than from an air crash. This is probably the consequence of over a century of experience building, starting with the Wright brothers at the beginning of the 20th century to the most recent aerospace developments culminating in technological giants such as the Airbus A380 and the Boeing Dreamliner, through the enlightening catastrophic events of the "Comet Aircraft”, ``Liberty Ships'' and many others. Yet, with the increasing urge to increase flight efficiency, decrease costs and Carbon emissions, airlines have been pushed to drive down the weight of aircraft, whilst guaranteeing their safety. This push for lighter aircraft has progressively seen a reduction in the use of metallic components which have been slowly replaced by composite materials. Such composite materials are made up of two or more phases of which they exploit the mechanical complementarity. For some applications, such as thermal barrier coatings, thermal complementarity is also leveraged. Yet, these novel materials, and especially their failure mechanisms and durability have proven difficult to understand, both through physical and virtual, in silico, experiments. One of the reasons for this is the large ratio between the size of the smallest constituent relevant in the description of failure mechanisms (e.g. 5-10 micron diameter carbon fibres) and the size of the structure (79m wingspan A380). In this presentation, we will briefly review advances in modeling and simulation of failure across the scales. We will discuss non exhaustively some of the recent advances in this field, ranging from adaptive atomistic modeling of fracture to algebraic model reduction methods for severely non-linear problems, including homogenization. We will also discuss the relevance of such simulations in daily engineering practice and claim that devising interactive simula- tors able to let engineers interact with the composite structure of interest and thus develop intuition about these advanced and complex materials. We will conclude by making a parallel between the difficulties encountered in modeling complex aerospace components and those met in personalized medicine, by discussing briefly the concept of Digital Twin. [less ▲] Detailed reference viewed: 426 (14 UL)Multiscale Quasicontinuum Methods for Dissipative Truss Models and Beam Networks Beex, Lars ; ; et al Presentation (2014, November 05) Detailed reference viewed: 138 (4 UL)Cardiff/Luxembourg Computational Mechanics Research Group Bordas, Stéphane ; ; Hale, Jack et al Poster (2014, November) Detailed reference viewed: 187 (7 UL)Shape optimisation directly from CAD: an isogeometric boundary element approach ; Bordas, Stéphane ; Report (2014) Detailed reference viewed: 115 (4 UL)Error controlled adaptive multiscale method for fracture in polycrystalline materials ; ; Bordas, Stéphane Report (2014) A lack of separation of scales is the major hurdle hampering predictive and computationally tractable simulations of fracture over multiple scales. In this thesis an adaptive multiscale method is ... [more ▼] A lack of separation of scales is the major hurdle hampering predictive and computationally tractable simulations of fracture over multiple scales. In this thesis an adaptive multiscale method is presented in an attempt to address this challenge. This method is set in the context of FE2 Feyel and Chaboche [2000] for which computational homogenisation breaks down upon loss of material stability (softening). The lack of scale separation due to the coalescence of microscopic cracks in a certain zone is tackled by a full discretisation of the microstructure in this zone. Polycrystalline materials are considered with cohesive cracks along the grain boundaries as a model problem. Adaptive mesh re nement of the coarse region and adaptive initiation and growth of fully resolved regions are performed based on discretisation error and homogenisation error criteria, respectively. In order to follow sharp snap-backs in load-displacement paths, a local arc-length technique is developed for the adaptive multiscale method. The results are validated against direct numerical simulation [less ▲] Detailed reference viewed: 103 (1 UL) |
||