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

Nitsche’s method for two and three dimensional NURBS patch coupling ; ; et al in Computational Mechanics (in press) We present a Nitche’s method 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 linear ... [more ▼] We present a Nitche’s method 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 linear elastostatics in two and and three-dimensions. The method can deal with surface-surface or volume-volume coupling, and we show how it can be used to handle heterogeneities such as inclusions. We also present preliminary results on modal analysis. This simple coupling method has the potential to increase the applicability of NURBS-based isogeometric analysis for practical applications. [less ▲] Detailed reference viewed: 1117 (75 UL)Finite strain poro-hyperelasticity: an asymptotic multi-scale ALE-FSI approach supported by ANNs Dehghani, Hamidreza ; Zilian, Andreas in Computational Mechanics (2023) This contribution introduces and discusses a formulation of poro-hyperelasticity at finite strains. The prediction of the time-dependent response of such media requires consideration of their ... [more ▼] This contribution introduces and discusses a formulation of poro-hyperelasticity at finite strains. The prediction of the time-dependent response of such media requires consideration of their characteristic multi-scale and multi-physics parameters. In the present work this is achieved by formulating a non-dimensionalised fluid–solid interaction problem (FSI) at the pore level using an arbitrary Lagrange–Euler description (ALE). The resulting coupled systems of PDEs on the reference configuration are expanded and analysed using the asymptotic homogenisation technique. This approach yields three partially novel systems of PDEs: the macroscopic/effective problem and two supplementary microscale problems (fluid and solid). The latter two provide the microscopic response fields whose average value is required in real-time/online form to determine the macroscale response (a concurrent multi-scale approach). In order to overcome the computational challenges related to the above multi-scale closure, this work introduces a surrogate approach for replacing the direct numerical simulation with an artificial neural network. This methodology allows for solving finite strain (multi-scale) porohyperelastic problems accurately using direct automated differentiation through the strain energy. Optimal and reliable training data sets are produced from direct numerical simulations of the fully-resolved problem by including a simple real-time output density check for adaptive sampling step refinement. The data-driven approach is complemented by a sensitivity analysis of the RVE response. The significance of the presented approach for finite strain poro-elasticity/poro-hyperelasticity is shown in the numerical benchmark of a multi-scale confined consolidation problem. Finally, to show the robustness of the method, the system response is dimensionalised using characteristic values of soil and brain mechanics scenarios. [less ▲] Detailed reference viewed: 32 (3 UL)ANN-aided incremental multiscale-remodelling-based finite strain poroelasticity Dehghani, Hamidreza ; Zilian, Andreas in Computational Mechanics (2021) Mechanical modelling of poroelastic media under finite strain is usually carried out via phenomenological models neglecting complex micro-macro scales interdependency. One reason is that the mathematical ... [more ▼] Mechanical modelling of poroelastic media under finite strain is usually carried out via phenomenological models neglecting complex micro-macro scales interdependency. One reason is that the mathematical two-scale analysis is only straightforward assuming infinitesimal strain theory. Exploiting the potential of ANNs for fast and reliable upscaling and localisation procedures, we propose an incremental numerical approach that considers rearrangement of the cell properties based on its current deformation, which leads to the remodelling of the macroscopic model after each time increment. This computational framework is valid for finite strain and large deformation problems while it ensures infinitesimal strain increments within time steps. The full effects of the interdependency between the properties and response of macro and micro scales are considered for the first time providing a more accurate predictive analysis of fluid-saturated porous media which is studied via a numerical consolidation example. Furthermore, the (nonlinear) deviation from Darcy’s law is captured in fluid filtration numerical analyses. Finally, the brain tissue mechanical response under the uniaxial cyclic test is simulated and studied. [less ▲] Detailed reference viewed: 76 (6 UL)Poroelastic model parameter identification using artificial neural networks: on the effects of heterogeneous porosity and solid matrix Poisson ratio Dehghani, Hamidreza ; Zilian, Andreas in Computational Mechanics (2020), 66 Predictive analysis of poroelastic materials typically require expensive and time-consuming multiscale and multiphysics approaches, which demand either several simplifications or costly experimental tests ... [more ▼] Predictive analysis of poroelastic materials typically require expensive and time-consuming multiscale and multiphysics approaches, which demand either several simplifications or costly experimental tests for model parameter identification. This problem motivates us to develop a more efficient approach to address complex problems with an acceptable computational cost. In particular, we employ artificial neural network (ANN) for reliable and fast computation of poroelastic model parameters. Based on the strong-form governing equations for the poroelastic problem derived from asymptotic homogenisation, the weighted residuals formulation of the cell problem is obtained. Approximate solution of the resulting linear variational boundary value problem is achieved by means of the finite element method. The advantages and downsides of macroscale properties identification via asymptotic homogenisation and the application of ANN to overcome parameter characterisation challenges caused by the costly solution of cell problems are presented. Numerical examples, in this study, include spatially dependent porosity and solid matrix Poisson ratio for a generic model problem, application in tumour modelling, and utilisation in soil mechanics context which demonstrate the feasibility of the presented framework. [less ▲] Detailed reference viewed: 136 (5 UL)Isogeometric analysis of thin Reissner-Mindlin shells: locking phenomena and B-bar method ; ; et al in Computational Mechanics (2020), 65(5), 1323-1341 We propose a local type of B-bar formulation, addressing locking in degenerated Reissner–Mindlin shell formulation in the context of isogeometric analysis. Parasitic strain components are projected onto ... [more ▼] We propose a local type of B-bar formulation, addressing locking in degenerated Reissner–Mindlin shell formulation in the context of isogeometric analysis. Parasitic strain components are projected onto the physical space locally, i.e. at the element level, using a least-squares approach. The formulation allows the flexible utilization of basis functions of different orders as the projection bases. The introduced formulation is much cheaper computationally than the classical $$\bar{B}$$B¯ method. We show the numerical consistency of the scheme through numerical examples, moreover they show that the proposed formulation alleviates locking and yields good accuracy even for slenderness ratios of $$10^5$$105, and has the ability to capture deformations of thin shells using relatively coarse meshes. In addition it can be opined that the proposed method is less sensitive to locking with irregular meshes. [less ▲] Detailed reference viewed: 148 (4 UL)Non-localised contact between beams with circular and elliptical cross-sections Magliulo, Marco ; Lengiewicz, Jakub ; Zilian, Andreas et al in Computational Mechanics (2020), 65 The key novelty of this contribution is a dedicated technique to e fficiently determine the distance (gap) function between parallel or almost parallel beams with circular and elliptical cross-sections ... [more ▼] The key novelty of this contribution is a dedicated technique to e fficiently determine the distance (gap) function between parallel or almost parallel beams with circular and elliptical cross-sections. The technique consists of parametrizing the surfaces of the two beams in contact, fixing a point on the centroid line of one of the beams and searching for a constrained minimum distance between the surfaces (two variants are investigated). The resulting unilateral (frictionless) contact condition is then enforced with the Penalty method, which introduces compliance to the, otherwise rigid, beams' cross-sections. Two contact integration schemes are considered: the conventional slave-master approach (which is biased as the contact virtual work is only integrated over the slave surface) and the so-called two-half-pass approach (which is unbiased as the contact virtual work is integrated over the two contacting surfaces). Details of the finite element formulation which is suitably implemented using Automatic Di fferentiation techniques are presented. A set of numerical experiments shows the overall performance of the framework and allows a quantitative comparison of the investigated variants. [less ▲] Detailed reference viewed: 181 (50 UL)Multiple crack detection in 3D using a stable XFEM and global optimization ; ; Bordas, Stéphane in Computational Mechanics (2018) A numerical scheme is proposed for the detection of multiple cracks in three dimensional (3D) structures. The scheme is based on a variant of the extended finite element method (XFEM) and a hybrid ... [more ▼] A numerical scheme is proposed for the detection of multiple cracks in three dimensional (3D) structures. The scheme is based on a variant of the extended finite element method (XFEM) and a hybrid optimizer solution. The proposed XFEM variant is particularly well-suited for the simulation of 3D fracture problems, and as such serves as an efficient solution to the so-called forward problem. A set of heuristic optimization algorithms are recombined into a multiscale optimization scheme. The introduced approach proves effective in tackling the complex inverse problem involved, where identification of multiple flaws is sought on the basis of sparse measurements collected near the structural boundary. The potential of the scheme is demonstrated through a set of numerical case studies of varying complexity. [less ▲] Detailed reference viewed: 209 (7 UL)Gradient plasticity crack tip characterization by means of the extended finite element method ; ; Bordas, Stéphane in Computational Mechanics (2017) Strain gradient plasticity theories are being widely used for fracture assessment, as they provide a richer description of crack tip fields by incorporating the influence of geometrically necessary ... [more ▼] Strain gradient plasticity theories are being widely used for fracture assessment, as they provide a richer description of crack tip fields by incorporating the influence of geometrically necessary dislocations. Characterizing the behavior at the small scales involved in crack tip deformation requires, however, the use of a very refined mesh within microns to the crack. In this work a novel and efficient gradient-enhanced numerical framework is developed by means of the extended finite element method (X-FEM). A mechanism-based gradient plasticity model is employed and the approximation of the displacement field is enriched with the stress singularity of the gradient-dominated solution. Results reveal that the proposed numerical methodology largely outperforms the standard finite element approach. The present work could have important implications on the use of microstructurally-motivated models in large scale applications. The non-linear X-FEM code developed in MATLAB can be downloaded from www.empaneda.com/codes. [less ▲] Detailed reference viewed: 241 (2 UL)Multi-phase-field modeling of anisotropic crack propagation for polycrystalline materials Nguyen, Thanh Tung ; ; et al in Computational Mechanics (2017), 60 Detailed reference viewed: 107 (0 UL)Automatised selection of load paths to construct reduced-order models in computational damage micromechanics: from dissipation-driven random selection to Bayesian optimization ; ; Bordas, Stéphane et al in Computational Mechanics (2016) In this paper, we present new reliable model order reduction strategies for computational micromechanics. The difficulties rely mainly upon the high dimensionality of the parameter space represented by ... [more ▼] In this paper, we present new reliable model order reduction strategies for computational micromechanics. The difficulties rely mainly upon the high dimensionality of the parameter space represented by any load path applied onto the representative volume element. We take special care of the challenge of selecting an exhaustive snapshot set. This is treated by first using a random sampling of energy dissipating load paths and then in a more advanced way using Bayesian optimization associated with an interlocked division of the parameter space. Results show that we can insure the selection of an exhaustive snapshot set from which a reliable reduced-order model can be built. [less ▲] Detailed reference viewed: 417 (32 UL)Modelling interfacial cracking with non-matching cohesive interface elements Nguyen, Viet Ha ; ; Bordas, Stéphane et al in Computational Mechanics (2016), 58(5), 731-746 Interfacial cracking occurs in many engineering problems such as delamination in composite laminates, matrix/interface debonding in fibre reinforced composites etc. Computational modelling of these ... [more ▼] Interfacial cracking occurs in many engineering problems such as delamination in composite laminates, matrix/interface debonding in fibre reinforced composites etc. Computational modelling of these interfacial cracks usually employs compatible or matching cohesive interface elements. In this paper, incompatible or non-matching cohesive interface elements are proposed for interfacial fracture mechanics problems. They allow non-matching finite element discretisations of the opposite crack faces thus lifting the constraint on the compatible discretisation of the domains sharing the interface. The formulation is based on a discontinuous Galerkin method and works with both initially elastic and rigid cohesive laws. The proposed formulation has the following advantages compared to classical interface elements: (i) non-matching discretisations of the domains and (ii) no high dummy stiffness. Two and three dimensional quasi-static fracture simulations are conducted to demonstrate the method. Our method not only simplifies the meshing process but also it requires less computational demands, compared with standard interface elements, for problems that involve materials/solids having a large mismatch in stiffnesses. [less ▲] Detailed reference viewed: 141 (4 UL)Nitsche’s method for two and three dimensional NURBS patch coupling ; ; et al in Computational Mechanics (2014), 53(6), 1163-1182 We present a Nitche’s method 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 linear ... [more ▼] We present a Nitche’s method 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 linear elastostatics in two and and three-dimensions. The method can deal with surface-surface or volume-volume coupling, and we show how it can be used to handle heterogeneities such as inclusions. We also present preliminary results on modal analysis. This simple coupling method has the potential to increase the applicability of NURBS-based isogeometric analysis for practical applications. [less ▲] Detailed reference viewed: 822 (18 UL)On efficient and reliable stochastic generation of RVEs for analysis of composites within the framework of homogenization Salnikov, Vladimir ; ; in Computational Mechanics (2014) Detailed reference viewed: 101 (0 UL)A hybrid smoothed extended finite element/level set method for modeling equilibrium shapes of nano-inhomogeneities ; Bordas, Stéphane ; in Computational Mechanics (2013) Interfacial energy plays an important role in equilibrium morphologies of nanosized microstructures of solid materials due to the high interface-to-volume ratio, and can no longer be neglected as it does ... [more ▼] Interfacial energy plays an important role in equilibrium morphologies of nanosized microstructures of solid materials due to the high interface-to-volume ratio, and can no longer be neglected as it does in conventional mechanics analysis. When designing nanodevices and to understand the behavior of materials at the nano-scale, this interfacial energy must therefore be taken into account. The present work develops an effective numerical approach by means of a hybrid smoothed extended finite element/level set method to model nanoscale inhomogeneities with interfacial energy effect, in which the finite element mesh can be completely independent of the interface geometry. The Gurtin-Murdoch surface elasticity model is used to account for the interface stress effect and the Wachspress interpolants are used for the first time to construct the shape functions in the smoothed extended finite element method. Selected numerical results are presented to study the accuracy and efficiency of the proposed method as well as the equilibrium shapes of misfit particles in elastic solids. The presented results compare very well with those obtained from theoretical solutions and experimental observations, and the computational efficiency of the method is shown to be superior to that of its most advanced competitor. © 2013 Springer-Verlag Berlin Heidelberg. [less ▲] Detailed reference viewed: 406 (7 UL)Efficient recovery-based error estimation for the smoothed finite element method for smooth and singular linear elasticity ; ; et al in Computational Mechanics (2013), 52(1), 37-52 An error control technique aimed to assess the quality of smoothed finite element approximations is presented in this paper. Finite element techniques based on strain smoothing appeared in 2007 were shown ... [more ▼] An error control technique aimed to assess the quality of smoothed finite element approximations is presented in this paper. Finite element techniques based on strain smoothing appeared in 2007 were shown to provide significant advantages compared to conventional finite element approximations. In particular, a widely cited strength of such methods is improved accuracy for the same computational cost. Yet, few attempts have been made to directly assess the quality of the results obtained during the simulation by evaluating an estimate of the discretization error. Here we propose a recovery type error estimator based on an enhanced recovery technique. The salient features of the recovery are: enforcement of local equilibrium and, for singular problems a "smooth + singular" decomposition of the recovered stress. We evaluate the proposed estimator on a number of test cases from linear elastic structural mechanics and obtain efficient error estimations whose effectivities, both at local and global levels, are improved compared to recovery procedures not implementing these features. © 2012 Springer-Verlag. [less ▲] Detailed reference viewed: 358 (4 UL)An Adaptive Multiscale Method for Quasi-static Crack Growth ; ; Bordas, Stéphane et al in Computational Mechanics (2013) This paper proposes an adaptive atomistic-continuum numerical method for quasi-static <br />crack growth. The phantom node method is used to model the crack in the continuum region <br />and a molecular ... [more ▼] This paper proposes an adaptive atomistic-continuum numerical method for quasi-static <br />crack growth. The phantom node method is used to model the crack in the continuum region <br />and a molecular statics model is used near the crack tip. To ensure self-consistency in the bulk, <br />a virtual atom cluster is used to model the material of the coarse scale. The coupling between <br />the coarse scale and ne scale is realized through ghost atoms. The ghost atom positions are <br />interpolated from the coarse scale solution and enforced as boundary conditions on the ne <br />scale. The ne scale region is adaptively enlarged as the crack propagates and the region behind <br />the crack tip is adaptively coarsened. An energy criterion is used to detect the crack tip loca- <br />tion. The triangular lattice in the ne scale region corresponds to the lattice structure of the <br />(111) plane of an FCC crystal. The Lennard-Jones potential is used to model the atom-atom <br />interactions. The method is implemented in two dimensions. The results are compared to pure <br />atomistic simulations; they show excellent agreement. [less ▲] Detailed reference viewed: 480 (12 UL)XLME interpolants, a seamless bridge between XFEM and enriched meshless methods ; ; et al in Computational Mechanics (2013) In this paper, we develop a method based on local maximum entropy shape functions together with enrichment functions used in partition of unity methods to discretize problems in linear elastic fracture ... [more ▼] In this paper, we develop a method based on local maximum entropy shape functions together with enrichment functions used in partition of unity methods to discretize problems in linear elastic fracture mechanics. We obtain improved accuracy relative to the standard extended finite element method at a comparable computational cost. In addition, we keep the advantages of the LME shape functions, such as smoothness and non-negativity. We show numerically that optimal convergence (same as in FEM) for energy norm and stress intensity factors can be obtained through the use of geometric (fixed area) enrichment with no special treatment of the nodes near the crack such as blending or shifting. © 2013 Springer-Verlag Berlin Heidelberg. [less ▲] Detailed reference viewed: 385 (1 UL)Mesh adaptivity driven by goal-oriented locally equilibrated superconvergent patch recovery ; ; et al in Computational Mechanics (2013) Goal-oriented error estimates (GOEE) have become popular tools to quantify and control the local error in quantities of interest (QoI), which are often more pertinent than local errors in energy for ... [more ▼] Goal-oriented error estimates (GOEE) have become popular tools to quantify and control the local error in quantities of interest (QoI), which are often more pertinent than local errors in energy for design purposes (e.g. the mean stress or mean displacement in a particular area, the stress intensity factor for fracture problems). These GOEE are one of the key unsolved problems of advanced engineering applications in, for example, the aerospace industry. This work presents a simple recovery-based error estimation technique for QoIs whose main characteristic is the use of an enhanced version of the Superconvergent Patch Recovery (SPR) technique previously used for error estimation in the energy norm. This enhanced SPR technique is used to recover both the primal and dual solutions. It provides a nearly statically admissible stress field that results in accurate estimations of the local contributions to the discretisation error in the QoI and, therefore, in an accurate estimation of this magnitude. This approach leads to a technique with a reasonable computational cost that could easily be implemented into already available finite element codes, or as an independent postprocessing tool. © 2013 Springer-Verlag Berlin Heidelberg. [less ▲] Detailed reference viewed: 234 (2 UL)A computational library for multiscale modeling of material failure ; ; Bordas, Stéphane et al in Computational Mechanics (2013) We present an open-source software framework called PERMIX for multiscale modeling and simulation of fracture in solids. The framework is an object oriented open-source effort written primarily in Fortran ... [more ▼] We present an open-source software framework called PERMIX for multiscale modeling and simulation of fracture in solids. The framework is an object oriented open-source effort written primarily in Fortran 2003 standard with Fortran/C++ interfaces to a number of other libraries such as LAMMPS, ABAQUS, LS-DYNA and GMSH. Fracture on the continuum level is modeled by the extended finite element method (XFEM). Using several novel or state of the art methods, the piece software handles semi-concurrent multiscale methods as well as concurrent multiscale methods for fracture, coupling two continuum domains or atomistic domains to continuum domains, respectively. The efficiency of our open-source software is shown through several simulations including a 3D crack modeling in clay nanocomposites, a semi-concurrent FE-FE coupling, a 3D Arlequin multiscale example and an MD-XFEM coupling for dynamic crack propagation. © 2013 Springer-Verlag Berlin Heidelberg. [less ▲] Detailed reference viewed: 668 (4 UL)Enriched residual free bubbles for semiconductor device simulation ; Bordas, Stéphane ; et al in Computational Mechanics (2012), 50(1), 119-133 This article outlines a method for stabilising the current continuity equations which are used for semiconductor device simulation. Residual-free bubble functions (RfBF) are incorporated into a finite ... [more ▼] This article outlines a method for stabilising the current continuity equations which are used for semiconductor device simulation. Residual-free bubble functions (RfBF) are incorporated into a finite element (FE) implementation that are able to prevent oscillations which are seen when using the conventional Bubnov-Galerkin FE implementation. In addition, it is shown that the RfBF are able to provide stabilisation with very distorted meshes and curved interface boundaries. Comparison with the commonly used SUPG scheme is made throughout, showing that in the case of 2D problems the RfBF allow faster convergence of the coupled semiconductor device equations, especially in the case of distorted meshes. [less ▲] Detailed reference viewed: 126 (0 UL) |
||