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: 798 (61 UL)An Implicit boundary approach for viscous compressible high Reynolds flows using hybrid remeshed particle hydrodynamics method Obeidat, Anas ; Bordas, Stéphane in Journal of Computational Physics (in press) Detailed reference viewed: 134 (21 UL)The edge-based strain smoothing method for compressible and nearly incompressible non-linear elasticity for solid mechanics ; ; et al E-print/Working paper (in press) Detailed reference viewed: 454 (36 UL)ADVANCES IN GEOMETRY INDEPENDENT APPROXIMATIONS ; ; Bordas, Stéphane et al Scientific Conference (2019, April 11) We present recent advances in geometry independent field approximations. The GIFT approach is a generalisation of isogeometric analysis where the approximation used to describe the field variables no ... [more ▼] We present recent advances in geometry independent field approximations. The GIFT approach is a generalisation of isogeometric analysis where the approximation used to describe the field variables no-longer has to be identical to the approximation used to describe the geometry of the domain. As such, the geometry can be described using usual CAD representations, e.g. NURBS, which are the most common in the CAD area, whilst local refinement and meshes approximations can be used to describe the field variables, enabling local adaptivity. We show in which cases the approach passes the patch test and present applications to various mechanics, fracture and multi-physics problems. [less ▲] Detailed reference viewed: 126 (3 UL)Introduction to Isogeometric Analysis Bordas, Stéphane ; Lian, Haojie ; Ding, Chensen Report (2019) Detailed reference viewed: 31 (4 UL)Modelling Complex Systems: a primer - agent-based models, equation-based models, statistical models and Bayesian inference, digital twins Bordas, Stéphane Learning material (2019) Modelling Complex Systems: a primer - agent-based models, equation-based models, statistical models and Bayesian inference, digital twins Detailed reference viewed: 446 (22 UL)Meshing or not meshing - Iso/sub/super-geometric analysis (adaptive unfitted methods for real-time simulations) Immersed collocation methods
Bordas, Stéphane Scientific Conference (2019, January 04) Detailed reference viewed: 156 (14 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: 293 (54 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: 156 (39 UL)A hyper-reduction method using adaptivity to cut the assembly costs of reduced order models Hale, Jack ; ; Baroli, Davide et al E-print/Working paper (2019) 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: 104 (18 UL)A volume-averaged nodal projection method for the Reissner-Mindlin plate model ; ; Hale, Jack et al in Computer Methods in Applied Mechanics & Engineering (2018), 341 We introduce a novel meshfree Galerkin method for the solution of Reissner-Mindlin plate problems that is written in terms of the primitive variables only (i.e., rotations and transverse displacement) and ... [more ▼] We introduce a novel meshfree Galerkin method for the solution of Reissner-Mindlin plate problems that is written in terms of the primitive variables only (i.e., rotations and transverse displacement) and is devoid of shear-locking. The proposed approach uses linear maximum-entropy approximations and is built variationally on a two-field potential energy functional wherein the shear strain, written in terms of the primitive variables, is computed via a volume-averaged nodal projection operator that is constructed from the Kirchhoff constraint of the three-field mixed weak form. The stability of the method is rendered by adding bubble-like enrichment to the rotation degrees of freedom. Some benchmark problems are presented to demonstrate the accuracy and performance of the proposed method for a wide range of plate thicknesses. [less ▲] Detailed reference viewed: 90 (13 UL)Simple and extensible plate and shell finite element models through automatic code generation tools Hale, Jack ; ; Bordas, Stéphane et al in Computers & Structures (2018), 209 A large number of advanced finite element shell formulations have been developed, but their adoption is hindered by complexities of transforming mathematical formulations into computer code. Furthermore ... [more ▼] A large number of advanced finite element shell formulations have been developed, but their adoption is hindered by complexities of transforming mathematical formulations into computer code. Furthermore, it is often not straightforward to adapt existing implementations to emerging frontier problems in thin structural mechanics including nonlinear material behaviour, complex microstructures, multi-physical couplings, or active materials. We show that by using a high-level mathematical modelling strategy and automatic code generation tools, a wide range of advanced plate and shell finite element models can be generated easily and efficiently, including: the linear and non-linear geometrically exact Naghdi shell models, the Marguerre-von K ́arm ́an shallow shell model, and the Reissner-Mindlin plate model. To solve shear and membrane-locking issues, we use: a novel re-interpretation of the Mixed Interpolation of Tensorial Component (MITC) procedure as a mixed-hybridisable finite element method, and a high polynomial order Partial Selective Reduced Integration (PSRI) method. The effectiveness of these approaches and the ease of writing solvers is illustrated through a large set of verification tests and demo codes, collected in an open-source library, FEniCS-Shells, that extends the FEniCS Project finite element problem solving environment. [less ▲] Detailed reference viewed: 359 (26 UL)The elastic properties of composites reinforced by a transversely isotropic random fibre-network ; ; et al in Composite Structures (2018), 208 This research stems from the idea of introducing a fibre-network structure into composites aiming to enhance the stiffness and strength of the composites. A novel new type of composites reinforced by a ... [more ▼] This research stems from the idea of introducing a fibre-network structure into composites aiming to enhance the stiffness and strength of the composites. A novel new type of composites reinforced by a tranversely isotropic fibre-network in which the fibres are devided into continuous segments and randomly distributed has been proposed and found to have improved elastic properties compared to other conventional fibre or particle composites mainly due to the introduction of cross linkers among the fibres. Combining with the effects of Poisson’s ratio of the constituent materials, the fibre network composite can exhibit extraordinary stiffness. A simplified analytical model has also been proposed for comparison with the numerical results, showing close prediction of the stiffness of the fibre-network composites. Moreover, as a plate structure, the thickness of the fibre network composite is adjustable and can be tailored according to the dimensions and mechanical properties as demanded in industry. [less ▲] Detailed reference viewed: 38 (1 UL)Lack of separation of scales: A view from reduced order modelling and homogenisation
Bordas, Stéphane Speeches/Talks (2018) Detailed reference viewed: 67 (1 UL)Classification of states and model order reduction of large scale Chemical Vapor Deposition processes with solution multiplicity ; ; Beex, Lars et al in Computers and Chemical Engineering (2018), 121 This paper presents an equation-free, data-driven approach for reduced order modeling of a Chemical Vapor Deposition (CVD) process. The proposed approach is based on process information provided by ... [more ▼] This paper presents an equation-free, data-driven approach for reduced order modeling of a Chemical Vapor Deposition (CVD) process. The proposed approach is based on process information provided by detailed, high-fidelity models, but can also use spatio-temporal measurements. The Reduced Order Model (ROM) is built using the method-of-snapshots variant of the Proper Orthogonal Decomposition (POD) method and Artificial Neural Networks (ANN) for the identification of the time-dependent coefficients. The derivation of the model is completely equation-free as it circumvents the projection of the actual equations onto the POD basis. Prior to building the model, the Support Vector Machine (SVM) supervised classification algorithm is used in order to identify clusters of data corresponding to (physically) different states that may develop at the same operating conditions due to the inherent nonlinearity of the process. The different clusters are then used for ANN training and subsequent development of the ROM. The results indicate that the ROM is successful at predicting the dynamic behavior of the system in windows of operating parameters where steady states are not unique. [less ▲] Detailed reference viewed: 21 (3 UL)Improving the conditioning of XFEM/GFEM for fracture mechanics problems through enrichment quasi-orthogonalization ; Bordas, Stéphane ; in Computer Methods in Applied Mechanics and Engineering (2018) Partition of unity enrichment is known to significantly enhance the accuracy of the finite element method by allowing the incorporation of known characteristics of the solution in the approximation space ... [more ▼] Partition of unity enrichment is known to significantly enhance the accuracy of the finite element method by allowing the incorporation of known characteristics of the solution in the approximation space. However, in several cases it can further cause conditioning problems for which a number of remedies have been proposed in the framework of the extended/generalized finite element method (XFEM/GFEM). Those solutions often involve significant modifications to the initial method and result in increased implementation complexity. In the present work, a simple procedure for the local near-orthogonalization of enrichment functions is introduced, which significantly improves the conditioning of the resulting system matrices, while requiring only minor modifications to the initial method. Although application to different types of enrichment functions is possible, the resulting scheme is specialized for the singular enrichment functions used in linear elastic fracture mechanics and tested through benchmark problems. [less ▲] Detailed reference viewed: 55 (2 UL)Identifying fibre material parameter distributions with little experimental efforts Rappel, Hussein ; Beex, Lars ; Bordas, Stéphane Scientific Conference (2018, July 23) Detailed reference viewed: 56 (13 UL)Uncertainty Quantification in Finite Element Models:Application to SoftTissue Biomechanics Hauseux, Paul ; Hale, Jack ; Bulle, Raphaël et al Scientific Conference (2018, July 23) We present probabilistic approaches aiming at the selection of the best constitutive model and to identify their parameters from experimental data. These parameters are always associated with some degree ... [more ▼] We present probabilistic approaches aiming at the selection of the best constitutive model and to identify their parameters from experimental data. These parameters are always associated with some degree of uncertainty. It is therefore important to study how this statistical uncertainty in parameters propagates to a safety-critical quantity of interest in the output of a model. Efficient Monte Carlo methods based on variance reduction techniques (Sensitivity Derivatives Monte Carlo methods [Hauseux et al. 2017] and MultiLevel Monte Carlo [Giles 2015] methods) are employed to propagate this uncertainty for both random variables and random fields. Inverse and forward problems are strongly connected. In a bayesian setting [Matthies et al. 2017], developing methods that reduce the number of evaluations of the forward model to an absolute minimum to achieve convergence is crucial for tractable computations. Numerical results in the context of soft tissue biomechanics are presented and discussed. [less ▲] Detailed reference viewed: 99 (1 UL)Identifying material parameter distributions of fibers with extremely limited experimental efforts Rappel, Hussein ; Beex, Lars ; Bordas, Stéphane Scientific Conference (2018, July 22) Detailed reference viewed: 32 (4 UL)Free boundary problems: numerical methods and data-driven simulations Bordas, Stéphane Presentation (2018, March 05) Detailed reference viewed: 97 (7 UL) |
||