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

Using Bayes' theorem to infer the material parameters of human soft tissue Hale, Jack ; ; Bordas, Stéphane Presentation (2015, October 21) Detailed reference viewed: 132 (5 UL)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)Higher-order quasicontinuum methods for elastic and dissipative lattice models: uniaxial deformation and pure bending Beex, Lars ; ; et al in GAMM Mitteilungen (2015), 38(2), 344-368 The quasicontinuum (QC) method is a numerical strategy to reduce the computational cost of direct lattice computations - in this study we achieve a speed up of a factor of 40. It has successfully been ... [more ▼] The quasicontinuum (QC) method is a numerical strategy to reduce the computational cost of direct lattice computations - in this study we achieve a speed up of a factor of 40. It has successfully been applied to (conservative) atomistic lattices in the past, but using a virtual-power-statement it was recently shown that QC approaches can also be used for spring and beam lattice models that include dissipation. Recent results have shown that QC approaches for planar beam lattices experiencing in-plane and out-of-plane deformation require higher-order interpolation. Higher-order QC frameworks are scarce nevertheless. In this contribution, the possibilities of a second-order and third-order QC framework are investigated for an elastoplastic spring lattice. The higher-order QC frameworks are compared to the results of the direct lattice computations and to those of a linear QC scheme. Examples are chosen so that both a macroscale and a microscale quantity influences the results. The two multiscale examples focused on are (i) macroscopically prescribed uniaxial deformation and (ii) macroscopically prescribed pure bending. Furthermore, the examples include an individual inclusion in a large lattice and hence, are concurrent in nature. [less ▲] Detailed reference viewed: 338 (33 UL)Open-Cell Aluminium Foams with Graded Coatings as Passively Controllable Energy Absorbers ; Beex, Lars ; et al in Materials and Design (2015), 87 Compared to most bulk materials, open-cell aluminium (Al) foams (OCAFs) are light-weight and can absorb a significant amount of energy in compression, e.g. during impact.When coatedwith nickel (Ni), OCAFs ... [more ▼] Compared to most bulk materials, open-cell aluminium (Al) foams (OCAFs) are light-weight and can absorb a significant amount of energy in compression, e.g. during impact.When coatedwith nickel (Ni), OCAFs can absorb even more energy, making them more appropriate for impacts at higher velocities than uncoated OCAFs. When Ni-coated OCAFs experience low-velocity impact however, the stopping distance during the impact is small compared to that of uncoated OCAFs and hence, deceleration occurs fast. This exposes devices (and possibly human beings) protected by OCAFs to large internal forces leading to internal damage. An OCAF that combines the properties of uncoated and coated OCAFs can absorb energy during both low-velocity and high-velocity impact scenarios. This contribution introduces two of such OCAFs which are created by partially and gradually coating OCAFs. The general mechanics of the two OCAFs are revealed using experimental and numerical observation methods. [less ▲] Detailed reference viewed: 628 (27 UL)Isogeometric and multi-scale fracture Bordas, Stéphane Presentation (2015, July 31) Detailed reference viewed: 183 (13 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)A Bayesian inversion approach to recovering material parameters in hyperelastic solids using dolfin-adjoint Hale, Jack ; ; Bordas, Stéphane Presentation (2015, July 01) In the first part of the talk I will describe in general terms the link between classical optimisation techniques and the Bayesian approach to statistical inversion as outlined in the seminal book of ... [more ▼] In the first part of the talk I will describe in general terms the link between classical optimisation techniques and the Bayesian approach to statistical inversion as outlined in the seminal book of [Kaipio and Somersalo, 2005]. Under the assumption of an additive Gaussian noise model, a Gaussian prior distribution and a linear parameter-to-observable map, it is possible to uniquely characterise the Bayesian posterior as Gaussian with the maximum aposteriori (MAP) point equal to the minimum of a classic regularised minimisation problem and covariance matrix equal to the inverse of the Hessian of the functional evaluated at the MAP point. I will also discuss techniques that can be used when these assumptions break down. In the second part of the talk I will describe a method implemented within dolfin-adjoint [Funke and Farrell, arXiv 2013] to quantify the uncertainty in the recovered material parameters of a hyperelastic solid from partial and noisy observations of the displacement field in the domain. The finite element discretisation of the adjoint and higher-order adjoint (Hessian) equations are derived automatically from the high-level UFL representation of the problem. The resulting equations are solved using PETSc. I will concentrate on finding the eigenvalue decomposition of the posterior covariance matrix (Hessian). The eigenvectors associated with the lowest eigenvalues of the Hessian correspond with the directions in parameter space least constrained by the observations [Flath et al. 2011]. This eigenvalue problem is tricky to solve efficiently because the Hessian is very large (on the order of the number of parameters) and dense (meaning that only its action on a vector can be calculated, each involving considerable expense). Finally, I will show some illustrative examples including the uncertainty associated with deriving the material properties of a 3D hyperelastic block with a stiff inclusion with knowledge only of the displacements on the boundary of the domain. J. Kaipio and E. Somersalo, Statistical and Computational Inverse Problems, vol. 160. New York: Springer-Verlag, 2005. S. W. Funke and P. E. Farrell, “A framework for automated PDE-constrained optimisation,” arXiv:1302.3894 [cs], Feb. 2013. H. P. Flath, L. C. Wilcox, V. Akçelik, J. Hill, B. van Bloemen Waanders, and O. Ghattas, “Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations,” SIAM J. Sci. Comput., vol. 33, no. 1, pp. 407–432, Feb. 2011. [less ▲] Detailed reference viewed: 785 (25 UL)fenics-shells: a UFL-based library for simulating thin structures ; Hale, Jack ; Bordas, Stéphane et al Scientific Conference (2015, July 01) Shell, plate and beam (thin) structures are widely used in civil, mechanical and aeronautical engineering because they are capable of carrying high loads with a minimal amount of structural mass. Because ... [more ▼] Shell, plate and beam (thin) structures are widely used in civil, mechanical and aeronautical engineering because they are capable of carrying high loads with a minimal amount of structural mass. Because the out-of-plane dimension is usually much smaller than the two in-plane dimensions, it is possible to asymptotically reduce the full 3D equations of elasticity to a whole variety of equivalent 2D models posed on a manifold embedded in 3D space. This reduction results in massively reduced computational expense and remains a necessity for practical large-scale computation of structures of real engineering interest such as the fuselage of an aircraft. The numerical solution of such mathematical models is a challenging task, especially for very thin shells when shear and membrane locking effects require special attention. As originally noted by [Hale and Baiz, 2013], the high-level form language UFL provides an excellent framework for writing extensible, reusable and pedagogical numerical models of thin structures. To our knowledge fenics-shells represents the first unified open-source implementation of a wide range of thin structural models, including Reissner-Mindlin, Kirchhoff-Love, Von Karman and hierarchical (higher-order) plates, and Madare-Naghdi and Madare-Koiter shell models. Because of the broad scope of fenics-shells, in this talk we will focus on how to cure numerical locking by applying the Mixed Interpolation of Tensorial Components (MITC) approach of [Dvorkin and Bathe, 1986] and [Lee and Bathe, 2010] to a shell with an initially flat reference configuration. The MITC approach consists of an element-by-element interpolation of the degrees of freedom of the rotations onto the degrees of freedom of a reduced rotation space, the latter typically constructed using H(curl) conforming finite elements such as the rotated Raviart-Thomas-Nédélec elements. Then, the bilinear form is constructed on the underlying H(curl) space. Because of the interpolation operator, the original problem is expressed in terms of the degrees of freedom for the rotations only. Within DOLFIN we have implemented this projection operation using two UFL forms within a custom assembler compiled just-in-time using Instant. We show numerical convergence studies that match the apriori bounds available in the literature. E. N. Dvorkin and K.-J. Bathe, “A continuum mechanics based four-node shell element for general non-linear analysis,” Engineering Computations, vol. 1, no. 1, pp. 77–88, 1984. P. S. Lee and K. J. Bathe, “The quadratic MITC plate and MITC shell elements in plate bending,” Advances in Engineering Software, vol. 41, no. 5, pp. 712–728, 2010. J. S. Hale and P. M. Baiz, “Towards effective shell modelling with the FEniCS project” presented at the FEniCS Conference 2013, Jesus College, Cambridge, 19-Mar-2013. [less ▲] Detailed reference viewed: 844 (25 UL)Extended Finite Element Method with Global Enrichment ; ; Bordas, Stéphane et al Scientific Conference (2015, July) A variant of the extended finite element method is presented which facilitates the use of enriched elements in a fixed volume around the crack front (geometrical enrichment) in 3D fracture problems. The ... [more ▼] A variant of the extended finite element method is presented which facilitates the use of enriched elements in a fixed volume around the crack front (geometrical enrichment) in 3D fracture problems. The major problem associated with geometrical enrichment is that it significantly deteriorates the conditioning of the resulting system matrices, thus increasing solution times and in some cases making the systems unsolvable. For 2D problems this can be dealt with by employing degree of freedom gathering [1] which essentially inhibits spatial variation of enrichment function weights. However, for the general 3D problem such an approach is not possible since spatial variation of the enrichment function weights in the direction of the crack front is necessary in order to reproduce the variation of solution variables, such as the stress intensity factors, along the crack front. The proposed method solves the above problem by employing a superimposed mesh of special elements which serve as a means to provide variation of the enrichment function weights along the crack front while still not allowing variation in any other direction. The method is combined with special element partitioning algorithms [2] and numerical integration schemes [3] as well as techniques for the elimination of blending errors between the standard and enriched part of the approximation in order to further improve the accuracy of the produced results. Additionally, a novel benchmark problem is introduced which enables the computation of displacement and energy error norms as well as errors in the stress intensity factors for the general 3D case. Through this benchmark problem it is shown that the proposed method provides optimal convergence rates, improved accuracy and reduced computational cost compared to standard XFEM. [less ▲] Detailed reference viewed: 647 (11 UL)Quasicontinuum methods for planar beam lattices (abstract) Beex, Lars ; ; et al Scientific Conference (2015, July) Detailed reference viewed: 312 (1 UL)Hyperelastic Elastography in a Large-Scale Bayesian Inversion Setting Hale, Jack ; ; Bordas, Stéphane Scientific Conference (2015, July) We consider the problem of recovering the material parameters of a hyperelastic material [1] in the Bayesian inversion setting. In the Bayesian setting we can extract the statistics associated with ... [more ▼] We consider the problem of recovering the material parameters of a hyperelastic material [1] in the Bayesian inversion setting. In the Bayesian setting we can extract the statistics associated with various sources of uncertainty, including noise, experimental deficiencies and incomplete observations of the domain. This will allow medical practitioners to make superior diagnosis decisions when presented with a quantitative measure of uncertainty in the recovered parameters. On the assumption of a Gaussian additive noise model, a Gaussian prior and a linear forward model, the posterior distribution of the material parameters given the observations will also be Gaussian. To ensure that the assumption of a linear forward model is valid, and that the posterior is approximated sufficiently well by a Gaussian distribution, we place a limit on the strain regime in which our current methodology applies. We are developing MCMC methods for exploring the non-Gaussian statistics of the posterior distribution. In the linear case, the covariance matrix of the posterior distribution is then characterised by the inverse of the Hessian of the objective functional evaluated at its minimiser. To extract statistical information from the large and dense Hessian we perform a low-rank approximation of the Hessian [2]. The eigenvectors associated with the lowest eigenvalues are the directions in parameter space that are least constrained by the observations. We implement this work within the dolfin-adjoint [3] software package. We derive the MPI-parallel finite element discretisation of the forward, adjoint (1st and 2nd order), and tangent linear models using the high-level differentiation tools available within the FEniCS project. We show results demonstrating the effects of partial observations and poor experimental design on the reliability of the recovered parameters. [1] N. H. Gokhale, P. E. Barbone, and A. A. Oberai, “Solution of the nonlinear elasticity imaging inverse problem: the compressible case,” Inverse Problems, 10.1088/0266-5611/24/4/045010 [2] H. P. Flath, L. C. Wilcox, V. Akçelik, J. Hill, B. van Bloemen Waanders, and O. Ghattas, “Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations,” SIAM J. Sci. Comput., 10.1137/090780717 [3] P. Farrell, D. Ham, S. Funke, and M. Rognes, “Automated Derivation of the Adjoint of High-Level Transient Finite Element Programs,” SIAM J. Sci. Comput., 10.1137/120873558 [less ▲] Detailed reference viewed: 401 (17 UL)Large scale phase field model of fracture and cutting in soft tissues ; Hale, Jack ; et al Scientific Conference (2015, July) The phase field method has proven to be an important tool in computational mechanics in that it is able to deal naturally with crack nucleation and branching [1]. In this contribution, we demonstrate a ... [more ▼] The phase field method has proven to be an important tool in computational mechanics in that it is able to deal naturally with crack nucleation and branching [1]. In this contribution, we demonstrate a large scale phase field model of fracture and cutting of soft tissues undergoing non-linear deformations with a material law defined by a hyperelastic energy density functional. We will also provide some initial thoughts on the how the effect of a porous medium can be incorporated into the phase field model. We implement this work using the FEniCS project and PETSc software packages [2, 3]. [less ▲] Detailed reference viewed: 260 (9 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)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)XFEM with global enrichment for 3D cracks ; ; Bordas, Stéphane et al Scientific Conference (2015, June 17) We present an extended finite element method (XFEM) based on fixed area enrichment which 1) suppresses the difficulties associated with ill-conditioning, even for "large" enrichment radii; 2) requires 50 ... [more ▼] We present an extended finite element method (XFEM) based on fixed area enrichment which 1) suppresses the difficulties associated with ill-conditioning, even for "large" enrichment radii; 2) requires 50 times fewer enriched degrees of freedom (for a typical mesh) as the standard XFEM with geometrical enrichment (for the same or better accuracy level); 3) increases the accuracy level of the stress intensity factors and leads to "smooth" stress intensity variations along the crack front. [less ▲] Detailed reference viewed: 532 (17 UL)An efficient Computational approach for control of nonlinear transient responses of smart piezoelectric composite plates ; ; et al in International Journal of Non-Linear Mechanics (2015) An efficient computational approach based on a generalized unconstrained approach in conjunction with isogeometric analysis (IGA) are proposed for dynamic control of smart piezoelectric composite plates ... [more ▼] An efficient computational approach based on a generalized unconstrained approach in conjunction with isogeometric analysis (IGA) are proposed for dynamic control of smart piezoelectric composite plates. In composite plates, the mechanical displacement field is approximated according to the proposal model using isogeometric elements and the nonlinear transient formulation for plates is formed in the total Lagrange approach based on the von Kármán strains and solved by Newmark time integration. Through the thickness of each piezoelectric layer, the electric potential is assumed linearly. For active control of the piezoelectric composite plates, a close-loop system is used. An optimization procedure using genetic algorithm (GA) is considered to search optimal design for actuator input voltages. Various numerical examples are investigated to show high accuracy and reliability of the proposed method. [less ▲] Detailed reference viewed: 269 (7 UL)Geometry-Independent Field approximaTion (GIFT) for spline based FEM for Linear Elasticity: a Diffpack implementation ; ; Bordas, Stéphane et al Scientific Conference (2015, June 01) Detailed reference viewed: 365 (7 UL)Isogeometric analysis: an overview and computer implementation aspects ; ; Bordas, Stéphane et al in Mathematics and Computers in Simulation (2015) Isogeometric analysis (IGA) represents a recently developed technology in computational mechanics that offers the possibility of integrating methods for analysis and Computer Aided Design (CAD) into a ... [more ▼] Isogeometric analysis (IGA) represents a recently developed technology in computational mechanics that offers the possibility of integrating methods for analysis and Computer Aided Design (CAD) into a single, unified process. The implications to practical engineering design scenarios are profound, since the time taken from design to analysis is greatly reduced, leading to dramatic gains in efficiency. In this manuscript, through a self-contained Matlab⃝R implementation, we present an introduction to IGA applied to simple analysis problems and the related computer implementation aspects. Furthermore, implementation of the extended IGA which incorporates enrichment functions through the partition of unity method (PUM) is also presented, where several examples for both two-dimensional and three-dimensional fracture are illustrated. We also describe the use of IGA in the context of strong-form (collocation) formulations, which has been an area of research interest due to the potential for significant efficiency gains offered by these methods. The code which accompanies the present paper can be applied to one, two and three-dimensional problems for linear elasticity, linear elastic fracture mechanics, structural mechanics (beams/plates/shells including large displacements and rotations) and Poisson problems with or without enrichment. The B ́ezier extraction concept that allows the FE analysis to be performed efficiently on T-spline geometries is also incorporated. The article includes a summary of recent trends and developments within the field of IGA. [less ▲] Detailed reference viewed: 1115 (11 UL)Control of Flame Spray Pyrolysis synthesis of Li4Ti5O12: Experimental and Computational study ; ; Bourantas, Georgios et al Poster (2015, May 14) Lithium titanate (Li4Ti5O12, LTO) is a promising anode material for the next generation of lithium ion batteries. Its physical properties and morphology (which consequently affect its electrochemical ... [more ▼] Lithium titanate (Li4Ti5O12, LTO) is a promising anode material for the next generation of lithium ion batteries. Its physical properties and morphology (which consequently affect its electrochemical performance) highly depend on its synthesis method. Flame spray pyrolysis (FSP) is an attractive process for the controlled one-step synthesis of functional multicomponent oxides from low cost precursors. The main aim of this study is to control the growth process of LTO by FSP in order to maintain the desired particle properties. LTO nanoparticles of different sizes are synthesized by variation of the FSP processing conditions and characterized accordingly. Numerical simulations based on Population Balance Models are also implemented in order to investigate the evolution of primary and agglomerate particle growth. [less ▲] Detailed reference viewed: 351 (1 UL)The stable GFEM. Convergence, accuracy and Diffpack implementation ; ; Bordas, Stéphane et al Presentation (2015, May 12) Detailed reference viewed: 95 (1 UL) |
||