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

A hyper-reduction method using adaptivity to cut the assembly costs of reduced order models Hale, Jack ; ; Baroli, Davide et al in Computer Methods in Applied Mechanics and Engineering (in press) At every iteration or timestep of the online phase of some reduced-order modelling schemes, large linear systems must be assembled and then projected onto a reduced order basis of small dimension. The ... [more ▼] At every iteration or timestep of the online phase of some reduced-order modelling schemes, large linear systems must be assembled and then projected onto a reduced order basis of small dimension. The projected small linear systems are cheap to solve, but assembly and projection are now the dominant computational cost. In this paper we introduce a new hyper-reduction strategy called reduced assembly (RA) that drastically cuts these costs. RA consists of a triangulation adaptation algorithm that uses a local error indicator to con- struct a reduced assembly triangulation specially suited to the reduced order basis. Crucially, this reduced assembly triangulation has fewer cells than the original one, resulting in lower assembly and projection costs. We demonstrate the efficacy of RA on a Galerkin-POD type reduced order model (RAPOD). We show performance increases of up to five times over the baseline Galerkin-POD method on a non-linear reaction-diffusion problem solved with a semi-implicit time-stepping scheme and up to seven times for a 3D hyperelasticity problem solved with a continuation Newton-Raphson algorithm. The examples are implemented in the DOLFIN finite element solver using PETSc and SLEPc for linear algebra. Full code and data files to produce the results in this paper are provided as supplementary material. [less ▲] Detailed reference viewed: 259 (36 UL)A Bayesian framework to identify random parameter fields based on the copula theorem and Gaussian fields: Application to polycrystalline materials Rappel, Hussein ; ; et al in Journal of Applied Mechanics (in press) For many models of solids, we frequently assume that the material parameters do not vary in space, nor that they vary from one product realization to another. If the length scale of the application ... [more ▼] For many models of solids, we frequently assume that the material parameters do not vary in space, nor that they vary from one product realization to another. If the length scale of the application approaches the length scale of the micro-structure however, spatially fluctuating parameter fi elds (which vary from one realization of the fi eld to another) can be incorporated to make the model capture the stochasticity of the underlying micro-structure. Randomly fluctuating parameter fields are often described as Gaussian fields. Gaussian fi elds however assume that the probability density function of a material parameter at a given location is a univariate Gaussian distribution. This entails for instance that negative parameter values can be realized, whereas most material parameters have physical bounds (e.g. the Young's modulus cannot be negative). In this contribution, randomly fluctuating parameter fi elds are therefore described using the copula theorem and Gaussian fi elds, which allow di fferent types of univariate marginal distributions to be incorporated, but with the same correlation structure as Gaussian fields. It is convenient to keep the Gaussian correlation structure, as it allows us to draw samples from Gaussian fi elds and transform them into the new random fields. The bene fit of this approach is that any type of univariate marginal distribution can be incorporated. If the selected univariate marginal distribution has bounds, unphysical material parameter values will never be realized. We then use Bayesian inference to identify the distribution parameters (which govern the random fi eld). Bayesian inference regards the parameters that are to be identi fied as random variables and requires a user-defi ned prior distribution of the parameters to which the observations are inferred. For the homogenized Young's modulus of a columnar polycrystalline material of interest in this study, the results show that with a relatively wide prior (i.e. a prior distribution without strong assumptions), a single specimen is su ciffient to accurately recover the distribution parameter values. [less ▲] Detailed reference viewed: 114 (7 UL)Frictional interactions for non-localised beam-to-beam and beam-inside-beam contact ; Lengiewicz, Jakub ; Zilian, Andreas et al in International Journal for Numerical Methods in Engineering (in press) This contribution presents the extensions of beam-to-beam and beam-inside-beam contact schemes of the same authors towards frictional interactions. Since the schemes are based on the beams’ true surfaces ... [more ▼] This contribution presents the extensions of beam-to-beam and beam-inside-beam contact schemes of the same authors towards frictional interactions. Since the schemes are based on the beams’ true surfaces (instead of surfaces implicitly deduced from the beams’ centroid lines), the presented enhancements are not only able to account for frictional sliding in the beams’ axial directions, but also in the circumferential directions. Both the frictional beam-to-beam approach as well as the frictional beam-inside-beam approach are applicable to shear-deformable and shear-undeformable beams, as well as to beams with both circular and elliptical cross-sections (although the cross-sections must be rigid). A penalty formulation is used to treat unilateral and frictional contact constraints. FE implementation details are discussed, where automatic differentiation techniques are used to derive the implementations. Simulations involving large sliding displacements and large deformations are presented for both beam-to-beam and beam-inside-beam schemes. All simulation results are compared to those of the frictionless schemes. [less ▲] Detailed reference viewed: 36 (1 UL)Beam-inside-beam contact: Mechanical simulations of slender medical instruments inside the human body Magliulo, Marco ; Lengiewicz, Jakub ; Zilian, Andreas et al in Computer Methods and Programs in Biomedicine (2020), 196 Background and Objective This contribution presents a rapid computational framework to mechanically simulate the insertion of a slender medical instrument in a tubular structure such as an artery, the ... [more ▼] Background and Objective This contribution presents a rapid computational framework to mechanically simulate the insertion of a slender medical instrument in a tubular structure such as an artery, the cochlea or another slender instrument. Methods Beams are employed to rapidly simulate the mechanical behaviour of the medical instrument and the tubular structure. However, the framework’s novelty is its capability to handle the mechanical contact between an inner beam (representing the medical instrument) embedded in a hollow outer beam (representing the tubular structure). This “beam-inside-beam” contact framework, which forces two beams to remain embedded, is the first of its kind since existing contact frameworks for beams are “beam-to-beam” approaches, i.e. they repel beams from each other. Furthermore, we propose contact kinematics such that not only instruments and tubes with circular cross-sections can be considered, but also those with elliptical cross-sections. This provides flexibility for the optimization of patient-specific instruments. Results The results demonstrate that the framework’s robustness is substantial, because only a few increments per simulation and a few iterations per increment are required, even though large deformations, large rotations and large curvature changes of both the instrument and tubular structure occur. The stability of the framework remains high even if the modulus of the inner tube is thousand times larger than that of the outer tube. A mesh convergence study furthermore exposes that a relatively small number of elements is required to accurately approach the reference solution. Conclusions The framework’s high simulation speed originates from the exploitation of the rigidity of the beams’ cross-sections to quantify the exclusion between the inner and the hollow outer beam. This rigidity limits the accuracy of the framework at the same time, but this is unavoidable since simulation accuracy and simulation speed are two competing interests. Hence, the framework is particularly attractive if simulation speed is preferred over accuracy. [less ▲] Detailed reference viewed: 185 (25 UL)Accelerating fatigue simulations of a phase-field damage model for rubber Loew, Pascal Juergen ; ; Peters, Bernhard et al in Computer Methods in Applied Mechanics and Engineering (2020), 370(113247), Phase-field damage models are able to describe crack nucleation as well as crack propagation and coalescence without additional technicalities, because cracks are treated in a continuous, spatially finite ... [more ▼] Phase-field damage models are able to describe crack nucleation as well as crack propagation and coalescence without additional technicalities, because cracks are treated in a continuous, spatially finite manner. Previously, we have developed a phase-field model to capture the rate-dependent failure of rubber, and we have further enhanced it to describe failure due to cyclic loading. Although the model accurately describes fatigue failure, the associated cyclic simulations are slow. Therefore, this contribution presents an acceleration scheme for cyclic simulations of our previously introduced phase-field damage model so that the simulation speed is increased to facilitate large-scale simulations of industrially relevant problems. We formulate an explicit and an implicit cycle jump method, which, depending on the selected jump size, reduce the calculation time up to 99.5%. To circumvent the manual tuning of the jump size, we also present an adaptive jump size selection procedure. Thanks to the implicit adaptive scheme, all material parameters are identified from experiments, which include fatigue crack nucleation and crack growth. Finally, the model and its parameters are validated with additional measurements of the fatigue crack growth rate. [less ▲] Detailed reference viewed: 74 (0 UL)Generalized quasicontinuum modeling of metallic lattices with geometrical and material nonlinearity and variability Chen, Li ; Beex, Lars ; et al in Computer Methods in Applied Mechanics and Engineering (2020), 366(112878), We propose a generalized quasicontinuum method to model the mechanical response of 3D lattice structures. The method relies on the spatial coupling of fully-resolved domains and coarse-grained domains. In ... [more ▼] We propose a generalized quasicontinuum method to model the mechanical response of 3D lattice structures. The method relies on the spatial coupling of fully-resolved domains and coarse-grained domains. In the fully-resolved domain, the full micro-structure is taken into account. In the coarse-grained domain, the kinematics of the micro-structure are individually interpolated based on their connectivity. On top of that, the contributions of the microstructure to the governing equations in the coarse-grained domain are sampled using only a few unit cells. In both domains, geometrical and material variability along the strut can be naturally taken into account using a 3D co-rotational beam finite element with embedded plastic hinges. We verify the approach for BCC lattices, demonstrating that the new method can capture both material and geometrical non-linearities of single struts at a fraction of the cost of a direct numerical simulation. [less ▲] Detailed reference viewed: 48 (3 UL)DATA DRIVEN SURGICAL SIMULATIONS Deshpande, Saurabh ; Bordas, Stéphane ; Beex, Lars et al Scientific Conference (2020, July) Detailed reference viewed: 35 (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: 116 (44 UL)Contact between shear-deformable beams with elliptical cross-sections Magliulo, Marco ; Zilian, Andreas ; Beex, Lars in Acta Mechanica (2020), 231 Slender constituents are present in many structures and materials. In associated mechanical models, each slender constituent is often described with a beam. Contact between beams is essential to ... [more ▼] Slender constituents are present in many structures and materials. In associated mechanical models, each slender constituent is often described with a beam. Contact between beams is essential to incorporate in mechanical models, but associated contact frameworks are only demonstrated to work for beams with circular cross-sections. Only two studies have shown the ability to treat contact between beams with elliptical cross-sections, but those frameworks are limited to point-wise contact, which narrows their applicability. This contribution presents initial results of a framework for shear-deformable beams with elliptical cross-sections if contact occurs along a line or at an area (instead of at a point). This is achieved by integrating a penalty potential over one of the beams’ surfaces. Simo-Reissner Geometrically Exact Beam (GEB) elements are employed to discretise each beam. As the surface of an assembly of such beam elements is discontinuous, a smoothed surface is introduced to formulate the contact kinematics. This enables the treatment of contact for large sliding displacements and substantial deformations. [less ▲] Detailed reference viewed: 113 (6 UL)Fatigue phase-field damage modeling of rubber using viscous dissipation: Crack nucleation and propagation Loew, Pascal Juergen ; Peters, Bernhard ; Beex, Lars in Mechanics of Materials (2020), 142 By regularizing sharp cracks within a pure continuum setting, phase-damage models offer the ability to capture crack nucleation as well as crack propagation. Crack branching and coalescence can ... [more ▼] By regularizing sharp cracks within a pure continuum setting, phase-damage models offer the ability to capture crack nucleation as well as crack propagation. Crack branching and coalescence can furthermore be described without any additional efforts, as geometrical descriptions of the cracks are not required. In this contribution, we extend our previous phase-field model for rate-dependent fracture of rubbers in a finite strain setting (Loew et al., 2019) to describe damage under cyclic loading. The model is derived from the balance of mechanical energy and introduces a fatigue damage source as a function of the accumulated viscous dissipation under cyclic loading. We use uniaxial cyclic tension to present the influence of the fatigue material parameters and to confirm the model’s energy balance. The parameters are subsequently identified using monotonic and cyclic experiments of a plane stress nature. Finally, the model is validated by separate experiments, which demonstrate that the model accurately predicts (fatigue) crack nucleation as well as propagation. [less ▲] Detailed reference viewed: 131 (7 UL)Adaptive equation-free multiscale modeling of metallic lattices with geometrical nonlinearity and variability Chen, Li ; ; Beex, Lars et al Scientific Conference (2019, September 12) An equation-free concurrent multiscale framework is proposed to model 3D metallic lattice structures. The proposed equation-free multiscale method (EFMM) is effectively a generalization of the ... [more ▼] An equation-free concurrent multiscale framework is proposed to model 3D metallic lattice structures. The proposed equation-free multiscale method (EFMM) is effectively a generalization of the quasicontinuum method [2] and relies on the use of fully-resolved domains (FRD) in which all details of the lattice micro-structure are captured, and of coarse-grained domains (CGD) in which a model reduction is performed by interpolation and summation steps. The particularity of the lattice geometrical description is that cross section variations along the lattice struts (caused by the manufacturing process) are explicitly represented by their discretization in several beam finite elements, both in the FRDs and CGDs. The interpolation step of the EFMM refers to a kinematic approximation of the lattice deformation within CGDs based on the displacement of a reduced number of material points. One of the originalities of this work is the consideration of a separate interpolation of each type of kinematic variables within the CGDs, as a function of the connectivity of the lattice beam nodes (i.e. taking the location of different cross sections into account) and their kinematical pattern. This, together with accounting for geometric nonlinearity, by the development and implementation of a 3D co-rotational beam finite element [1], are innovative contributions. Choosing the appropriate sizes of the FRDs and the CGDs for a lattice to be simulated is a trade-off because larger FRDs prevail the accuracy but compromise the efficiency while larger CGDs do the opposite. Since the required sizes of the FRDs and CGDs are generally not known a priori for specific applications, an adaptive coarse-graining strategy is developed. To be specific, the whole lattice is initially considered as a CGD. Two kinds of error indicator are proposed (e.g. the Zienkiewicz-Zhu error indicator [4, 3] and the error indicator based on the discrepancy of strain energy). The error indicator guides on: 1) introducing more material points and rearranging the interpolation for the CGDs; 2) changing the localization-prone parts of the lattice into FRDs. The adaptive EFMM is applied to metallic BCC lattices with various sizes and loading conditions. By comparing to the results of those of the direct numerical simulation (DNS), it is shown that geometrical non-linearities can be captured at a fraction of the DNS cost. [less ▲] Detailed reference viewed: 30 (2 UL)Equation-free multiscale modeling of metallic lattices with geometrical and material nonlinearity and variability Chen, Li ; ; Beex, Lars et al Scientific Conference (2019, September 05) An nonlinear equation-free concurrent multiscale numerical framework, being the generalization of the quasicontinuum method [2] is proposed in this contribution to model 3D metallic lattice structures ... [more ▼] An nonlinear equation-free concurrent multiscale numerical framework, being the generalization of the quasicontinuum method [2] is proposed in this contribution to model 3D metallic lattice structures. The proposed equation-free multiscale method (EFMM) relies on the use of fully-resolved domains (FRD) in which all of the details of the lattice micro-structure are captured, and of coarse-grained domains (CGD) in which a model reduction is performed by interpolation and summation steps. The particularity of the lattice geometry description is that cross section variations along the lattice struts (that are experimentally observed as a result of the manufacturing process) are explicitly represented by their discretization in several beam finite elements, both in the FRDs and CGDs. The interpolation step of the EFMM refers to a kinematic approximation of the lattice deformation within CGDs based on the movement of a reduced number of material points at the CGD corners. One of the originalities of this work is the consideration of a separate interpolation of each type of degrees of freedom within the CGDs, as a function of the connectivity of the lattice beam nodes (i.e. taking the location of different cross sections into account) and their kinematical pattern. This, together with accounting for plasticity, by the development and implementation of a 3D co-rotational beam finite element [1] with embedded plastic hinges [3], are unprecedented and original contributions. The EFMM is applied to metallic BCC lattices with various sizes and loading conditions. By comparing to direct numerical simulation (DNS), it is shown that both material and geometrical non-linearities can be captured at a fraction of the DNS cost (the computational time is reduced by 97.27% while introducing an error of only 3.76%). [less ▲] Detailed reference viewed: 37 (5 UL)Fusing the Seth-Hill strain tensors to fit compressible elastic material responses in the nonlinear regime Beex, Lars in International Journal of Mechanical Sciences (2019), 163 Strain energy densities based on the Seth-Hill strain tensors are often used to describe the hyperelastic mechanical behaviours of isotropic, transversely isotropic and orthotropic materials for ... [more ▼] Strain energy densities based on the Seth-Hill strain tensors are often used to describe the hyperelastic mechanical behaviours of isotropic, transversely isotropic and orthotropic materials for relatively large deformations. Since one parameter distinguishes which strain tensor of the Seth-Hill family is used, one has in theory the possibility to t the material response in the nonlinear regime. Most often for compressible deformations however, this parameter is selected such that the Hencky strain tensor is recovered, because it yields rather physical stress-strain responses. Hence, the response in the nonlinear regime is in practise not often tailored to match experimental data. To ensure that elastic responses in the nonlinear regime can more accurately be controlled, this contribution proposes three generalisations that combine several Seth-Hill strain tensors. The generalisations are formulated such that the stress-strain responses for in finitesimal deformations remain unchanged. Consequently, the identifi cation of the Young's moduli, Poisson's ratios and shear moduli is not a ffected. 3D fi nite element simulations are performed for isotropy and orthotropy, with an emphasis on the identifi cation of the new material parameters. [less ▲] Detailed reference viewed: 129 (4 UL)Clustering Based Model Order Reduction For Hyper Elastoplastic Material Models Vijayaraghavan, Soumianarayanan ; Beex, Lars ; et al Presentation (2019, July 29) Detailed reference viewed: 25 (5 UL)Bayesian Identification of Mean-Field Homogenization model parameters and uncertain matrix behavior in non-aligned short fiber composites ; ; et al in Composite Structures (2019), 220 We present a stochastic approach combining Bayesian Inference (BI) with homogenization theories in order to identify, on the one hand, the parameters inherent to the model assumptions and, on the other ... [more ▼] We present a stochastic approach combining Bayesian Inference (BI) with homogenization theories in order to identify, on the one hand, the parameters inherent to the model assumptions and, on the other hand, the composite material constituents behaviors, including their variability. In particular, we characterize the model parameters of a Mean-Field Homogenization (MFH) model and the elastic matrix behavior, including the inherent dispersion in its Young's modulus, of non-aligned Short Fibers Reinforced Polymer (SFRP) composites. The inference is achieved by considering as observations experimental tests conducted at the SFRP composite coupons level. The inferred model and material law parameters can in turn be used in Mean-Field Homogenization (MFH)-based multi-scale simulations and can predict the confidence range of the composite material responses. [less ▲] Detailed reference viewed: 150 (12 UL)A stochastic Galerkin cell-based smoothed finite element method (SGCS-FEM) ; Beex, Lars ; Bordas, Stéphane et al in International Journal of Computational Methods (2019), 17(8), In this paper, the cell based smoothed finite element method is extended to solve stochastic partial diff erential equations with uncertain input parameters. The spatial field of Young's moduli and the ... [more ▼] In this paper, the cell based smoothed finite element method is extended to solve stochastic partial diff erential equations with uncertain input parameters. The spatial field of Young's moduli and the corresponding stochastic results are represented by Karhunen-Lo eve expansion and polynomial chaos expansion, respectively. The Young's Modulus of structure is considered to be random for stochastic static as well as free vibration problems. Mathematical expressions and the solution procedure are articulated in detail to evaluate the statistical characteristics of responses in terms of static displacements and free vibration frequencies. The feasibility and eff ectiveness of the proposed SGCS-FEM method in terms of accuracy and lower requirement on the mesh size in the solution domain over that of conventional FEM for stochastic problems are demonstrated by carefully chosen numerical examples. From the numerical study, it is inferred that the proposed framework is computationally less demanding without compromising accuracy. [less ▲] Detailed reference viewed: 146 (2 UL)Estimating fibres' material parameter distributions from limited data with the help of Bayesian inference Rappel, Hussein ; Beex, Lars in European Journal of Mechanics. A, Solids (2019), 75 Numerous materials are essentially structures of discrete fibres, yarns or struts. Considering these materials at their discrete scale, one may distinguish two types of intrinsic randomness that affect ... [more ▼] Numerous materials are essentially structures of discrete fibres, yarns or struts. Considering these materials at their discrete scale, one may distinguish two types of intrinsic randomness that affect the structural behaviours of these discrete structures: geometrical randomness and material randomness. Identifying the material randomness is an experimentally demanding task, because many small fibres, yarns or struts need to be tested, which are not easy to handle. To avoid the testing of hundreds of constituents, this contribution proposes an identification approach that only requires a few dozen of constituents to be tested (we use twenty to be exact). The identification approach is applied to articially generated measurements, so that the identified values can be compared to the true values. Another question this contribution aims to answer is how precise the material randomness needs to be identified, if the geometrical randomness will also influence the macroscale behaviour of these discrete networks. We therefore also study the effect of the identified material randomness to that of the actual material randomness for three types of structures; each with an increasing level of geometrical randomness. [less ▲] Detailed reference viewed: 189 (32 UL)Multiscale fracture: a natural connection between reduced order models and homogenisation Bordas, Stéphane ; Beex, Lars ; Chen, Li et al Scientific Conference (2019, May 13) Detailed reference viewed: 115 (9 UL)A Tutorial on Bayesian Inference to Identify Material Parameters in Solid Mechanics Rappel, Hussein ; Beex, Lars ; Hale, Jack et al in Archives of Computational Methods in Engineering (2019) The aim of this contribution is to explain in a straightforward manner how Bayesian inference can be used to identify material parameters of material models for solids. Bayesian approaches have already ... [more ▼] The aim of this contribution is to explain in a straightforward manner how Bayesian inference can be used to identify material parameters of material models for solids. Bayesian approaches have already been used for this purpose, but most of the literature is not necessarily easy to understand for those new to the field. The reason for this is that most literature focuses either on complex statistical and machine learning concepts and/or on relatively complex mechanical models. In order to introduce the approach as gently as possible, we only focus on stress–strain measurements coming from uniaxial tensile tests and we only treat elastic and elastoplastic material models. Furthermore, the stress–strain measurements are created artificially in order to allow a one-to-one comparison between the true parameter values and the identified parameter distributions. [less ▲] Detailed reference viewed: 596 (101 UL)Identifying elastoplastic parameters with Bayes' theorem considering double error sources and model uncertainty Rappel, Hussein ; Beex, Lars ; et al in Probabilistic Engineering Mechanics (2019), 55 We discuss Bayesian inference for the identi cation of elastoplastic material parameters. In addition to errors in the stress measurements, which are commonly considered, we furthermore consider errors in ... [more ▼] We discuss Bayesian inference for the identi cation of elastoplastic material parameters. In addition to errors in the stress measurements, which are commonly considered, we furthermore consider errors in the strain measurements. Since a difference between the model and the experimental data may still be present if the data is not contaminated by noise, we also incorporate the possible error of the model itself. The three formulations to describe model uncertainty in this contribution are: (1) a random variable which is taken from a normal distribution with constant parameters, (2) a random variable which is taken from a normal distribution with an input-dependent mean, and (3) a Gaussian random process with a stationary covariance function. Our results show that incorporating model uncertainty often, but not always, improves the results. If the error in the strain is considered as well, the results improve even more. [less ▲] Detailed reference viewed: 335 (67 UL) |
||