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 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: 158 (19 UL)Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization ; ; et al in Computer Methods in Applied Mechanics & Engineering (2018), 328 In this paper, we propose a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries. Instead of the computational domain bounded by four B ... [more ▼] In this paper, we propose a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries. Instead of the computational domain bounded by four B-spline curves, planar domains with high genus and more complex boundary curves are considered. Firstly, some pre-processing operations including B´ezier extraction and subdivision are performed on each boundary curve in order to generate a high-quality planar parameterization; then a robust planar domain partition framework is proposed to construct high-quality patch-meshing results with few singularities from the discrete boundary formed by connecting the end points of the resulting boundary segments. After the topology information generation of quadrilateral decomposition, the optimal placement of interior B´ezier curves corresponding to the interior edges of the quadrangulation is constructed by a global optimization method to achieve a patch-partition with high quality. Finally, after the imposition of C1/G1-continuity constraints on the interface of neighboring Bezier patches with respect to each quad in the quadrangulation, the high-quality Bezier patch parameterization is obtained by a local optimization method to achieve uniform and orthogonal iso-parametric structures while keeping the continuity conditions between patches. The efficiency and robustness of the proposed method are demonstrated by several examples which are compared to results obtained by the skeleton-based parameterization approach. [less ▲] Detailed reference viewed: 99 (6 UL)Accelerating Monte Carlo estimation with derivatives of high-level finite element models Hauseux, Paul ; Hale, Jack ; Bordas, Stéphane in Computer Methods in Applied Mechanics & Engineering (2017), 318 In this paper we demonstrate the ability of a derivative-driven Monte Carlo estimator to accelerate the propagation of uncertainty through two high-level non-linear finite element models. The use of ... [more ▼] In this paper we demonstrate the ability of a derivative-driven Monte Carlo estimator to accelerate the propagation of uncertainty through two high-level non-linear finite element models. The use of derivative information amounts to a correction to the standard Monte Carlo estimation procedure that reduces the variance under certain conditions. We express the finite element models in variational form using the high-level Unified Form Language (UFL). We derive the tangent linear model automatically from this high-level description and use it to efficiently calculate the required derivative information. To study the effectiveness of the derivative-driven method we consider two stochastic PDEs; a one- dimensional Burgers equation with stochastic viscosity and a three-dimensional geometrically non-linear Mooney-Rivlin hyperelastic equation with stochastic density and volumetric material parameter. Our results show that for these problems the first-order derivative-driven Monte Carlo method is around one order of magnitude faster than the standard Monte Carlo method and at the cost of only one extra tangent linear solution per estimation problem. We find similar trends when comparing with a modern non-intrusive multi-level polynomial chaos expansion method. We parallelise the task of the repeated forward model evaluations across a cluster using the ipyparallel and mpi4py software tools. A complete working example showing the solution of the stochastic viscous Burgers equation is included as supplementary material. [less ▲] Detailed reference viewed: 1408 (226 UL)Error-controlled adaptive extended finite element method for 3D linear elastic crack propagation ; ; et al in Computer Methods in Applied Mechanics & Engineering (2017), 318 We present a simple error estimation and mesh adaptation approach for 3D linear elastic crack propagation simulations using the eXtended Finite Element Method (X-FEM). A global extended recovery technique ... [more ▼] We present a simple error estimation and mesh adaptation approach for 3D linear elastic crack propagation simulations using the eXtended Finite Element Method (X-FEM). A global extended recovery technique (Duflot and Bordas, 2008) is used to quantify the interpolation error. Based on this error distribution, four strategies relying on two different mesh optimality criteria are compared. The first aims at homogenizing the error distribution. The second minimizes the total number of elements given a target global error level. We study the behaviour of these criteria in the context of cracks treated by an X-FE approach. In particular, we investigate the convergence rates at the element-level depending its enrichment type. We conclude on the most suitable refinement criterion and propose and verify a strategy for mesh adaptation on 3D damage tolerance assessment problems. [less ▲] Detailed reference viewed: 89 (3 UL)FE modelling with strong discontinuities for 3D tensile and shear fractures: Application to underground excavation Hauseux, Paul ; ; et al in Computer Methods in Applied Mechanics & Engineering (2016), 309 Detailed reference viewed: 85 (3 UL)Isogeometric boundary element methods for three dimensional static fracture and fatigue crack growth ; ; et al in Computer Methods in Applied Mechanics & Engineering (2016) Detailed reference viewed: 191 (12 UL)Stable 3D extended finite elements with higher order enrichment for accurate non planar fracture Agathos, Konstantinos ; ; Bordas, Stéphane in Computer Methods in Applied Mechanics & Engineering (2016), 306 An extended finite element method (XFEM) for three dimensional (3D) non-planar linear elastic fracture is introduced, which provides optimal convergence through the use of enrichment in a fixed area ... [more ▼] An extended finite element method (XFEM) for three dimensional (3D) non-planar linear elastic fracture is introduced, which provides optimal convergence through the use of enrichment in a fixed area around the crack front, while also improving the conditioning of the resulting system matrices. This is achieved by fusing a novel form of enrichment with existing blending techniques. Further, the adoption of higher order terms of theWilliams expansion is also considered and the effects in the accuracy and conditioning of the method are studied. Moreover, some problems regarding the evaluation of stress intensity factors (SIFs) and element partitioning are dealt with. The accuracy and convergence properties of the method as well as the conditioning of the resulting stiffness matrices are investigated through the use of appropriate benchmark problems. It is shown that the proposed approach provides increased accuracy while requiring, for all cases considered, a reduced number of iterations for the solution of the resulting systems of equations. The positive impact of geometrical enrichment is further demonstrated in the accuracy of the computed SIFs where, for the examined cases, an improvement of up to 40% is achieved. [less ▲] Detailed reference viewed: 79 (2 UL)A fast, certified and "tuning free" two-field reduced basis method for the metamodelling of affinely-parametrised elasticity problems ; ; Bordas, Stéphane in Computer Methods in Applied Mechanics & Engineering (2016), 298 This paper proposes a new reduced basis algorithm for the metamodelling of parametrised elliptic problems. The developments rely on the Constitutive Relation Error (CRE), and the construction of separate ... [more ▼] This paper proposes a new reduced basis algorithm for the metamodelling of parametrised elliptic problems. The developments rely on the Constitutive Relation Error (CRE), and the construction of separate reduced order models for the primal variable (displacement) and flux (stress) fields. A two field greedy sampling strategy is proposed to construct these two fields simultaneously and in an efficient manner: at each iteration, one of the two fields is enriched by increasing the dimension of its reduced space in such a way that the CRE is minimised. This sampling strategy is then used as a basis to construct goal-oriented reduced order modelling. The resulting algorithm is certified and “tuning free”: the only requirement from the engineer is the level of accuracy that is desired for each of the outputs of the surrogate. It is also shown to be significantly more efficient in terms of computational expense than competing methodologies. [less ▲] Detailed reference viewed: 95 (15 UL)Meshfree volume-averaged nodal projection method for nearly-incompressible elasticity using meshfree and bubble basis functions ; Hale, Jack ; in Computer Methods in Applied Mechanics & Engineering (2015), 285 We present a displacement-based Galerkin meshfree method for the analysis of nearly-incompressible linear elastic solids, where low-order simplicial tessellations (i.e., 3- node triangular or 4-node ... [more ▼] We present a displacement-based Galerkin meshfree method for the analysis of nearly-incompressible linear elastic solids, where low-order simplicial tessellations (i.e., 3- node triangular or 4-node tetrahedral meshes) are used as a background structure for numerical integration of the weak form integrals and to get the nodal information for the computation of the meshfree basis functions. In this approach, a volume- averaged nodal projection operator is constructed to project the dilatational strain into an approximation space of equal- or lower-order than the approximation space for the displacement field resulting in a locking-free method. The stability of the method is provided via bubble-like basis functions. Because the notion of an ‘ele- ment’ or ‘cell’ is not present in the computation of the meshfree basis functions such low-order tessellations can be used regardless of the order of the approximation spaces desired. First- and second-order meshfree basis functions are chosen as a particular case in the proposed method. Numerical examples are provided in two and three dimensions to demonstrate the robustness of the method, its ability to avoid volumetric locking in the nearly-incompressible regime, and its improved performance when compared to the MINI finite element scheme on the simplicial mesh. [less ▲] Detailed reference viewed: 225 (51 UL)A fast, certified and "tuning-free" two-field reduced basis method for the metamodelling of parametrised elasticity problems ; ; Bordas, Stéphane in Computer Methods in Applied Mechanics & Engineering (2015) This paper proposes a new reduced basis algorithm for the metamodelling of parametrised elliptic problems. The developments rely on the Constitutive Relation Error (CRE), and the construction of separate ... [more ▼] This paper proposes a new reduced basis algorithm for the metamodelling of parametrised elliptic problems. The developments rely on the Constitutive Relation Error (CRE), and the construction of separate reduced order models for the primal variable (displacement) and flux (stress) fields. A two-field Greedy sampling strategy is proposed to construct these two fields simultaneously and efficient manner: at each iteration, one of the two fields is enriched by increasing the dimension of its reduced space in such a way that the CRE is minimised. This sampling strategy is then used as a basis to construct goal-oriented reduced order modelling. The resulting algorithm is certified and "tuning-free": the only requirement from the engineer is the level of accuracy that is desired for each of the outputs of the surrogate. It is also one order of magnitude more efficient in terms of computational expenses than competing methodologies. [less ▲] Detailed reference viewed: 369 (11 UL)Stable 3D extended finite elements with higher order enrichment for accurate non planar fracture ; ; Bordas, Stéphane in Computer Methods in Applied Mechanics & Engineering (2015) We present an extended finite element method (XFEM) for 3D nonplanar linear elastic fracture. The new approach not only provides optimal convergence using geometrical enrichment but also enables to ... [more ▼] We present an extended finite element method (XFEM) for 3D nonplanar linear elastic fracture. The new approach not only provides optimal convergence using geometrical enrichment but also enables to contain the increase in conditioning number characteristic of enriched finite element formulations: the number of iterations to convergence of the conjugate gradient solver scales similarly to and converges faster than the topologically-enriched version of the standard XFEM. This has two advantages: (1) the residual can be driven to zero to machine precision for at least 50% fewer iterations than the standard version of XFEM; (2) additional enrichment functions can be added without significant deterioration of the conditioning. Numerical examples also show that our new approach is up to 40% more accurate in terms of stress intensity factors, than the standard XFEM. [less ▲] Detailed reference viewed: 201 (5 UL)Geometry-Independent Field approximaTion: CAD-Analysis Integration, geometrical exactness and adaptivity ; ; et al in Computer Methods in Applied Mechanics & Engineering (2014) In isogeometric analysis (IGA), the same spline representation is employed for both the geometry of the domain and approximation of the unknown fields over this domain. This identity of the geometry and ... [more ▼] In isogeometric analysis (IGA), the same spline representation is employed for both the geometry of the domain and approximation of the unknown fields over this domain. This identity of the geometry and field approximation spaces was put forward in the now classic 2005 paper [20] as a key advantage on the way to the integration of Computer Aided Design (CAD) and subsequent analysis in Computer Aided Engineering (CAE). [20] claims indeed that any change to the geometry of the domain is automatically inherited by the approximation of the field variables, without requiring the regeneration of the mesh at each change of the domain geometry. Yet, in Finite Element versions of IGA, a parameterization of the interior of the domain must still be constructed, since CAD only provides information about the boundary. The identity of the boundary and field representation decreases the flexibility in which this parameterization can be generated and somewhat constrains the modeling and simulation process, because an approximation able to represent the domain geometry accurately need not be adequate to also approximate the field variables accurately, in particular when the solution is not smooth. We propose here a new paradigm called Geometry-Independent Field approximaTion (GIFT) where the spline spaces used for the geometry and the field variables can be chosen and adapted independently while preserving geometric exactness and tight CAD integration. GIFT has the following features: (1) It is possible to flexibly choose between different spline spaces with different properties to better represent the solution of the problem, e.g. the continuity of the solution field, boundary layers, singularities, whilst retaining geometrical exactness of the domain boundary. (2) For multi-patch analysis, where the domain is composed of several spline patches, the continuity condition between neighboring patches on the solution field can be automatically guaranteed without additional constraints in the variational form. (3) Refinement operations by knot insertion and degree elevation are performed directly on the spline space of the solution field, independently of the spline space of the geometry of the domain, which makes the method versatile. GIFT with PHT-spline solution spaces and NURBS geometries is used to show the effectiveness of the proposed approach. Keywords : Super-parametric methods, Isogeometric analysis (IGA), Geometry-independent Spline Space, PHT-splines, local refinement, adaptivity [less ▲] Detailed reference viewed: 1026 (29 UL)A multiscale quasicontinuum method for lattice models with bond failure and fiber sliding Beex, Lars ; ; in Computer Methods in Applied Mechanics & Engineering (2014), 269 Structural lattice models incorporating trusses and beams are frequently used to mechanically model fibrous materials, because they can capture (local) mesoscale phenomena. Physically relevant lattice ... [more ▼] Structural lattice models incorporating trusses and beams are frequently used to mechanically model fibrous materials, because they can capture (local) mesoscale phenomena. Physically relevant lattice computations are however computationally expensive. A suitable multiscale approach to reduce the computational cost of large-scale lattice computations is the quasicontinuum (QC) method. This method resolves local mesoscale phenomena in regions of interest and coarse grains elsewhere, using only the lattice model. In previous work, a virtual-power-based QC framework is proposed for lattice models that include local dissipative mechanisms. In this paper, the virtual-power-based QC method is adopted for lattice models in which bond failure and subsequent frictional fiber sliding are incorporated – which are of significant importance for fibrous materials such as paper, cardboard, textile and electronic textile. Bond failure and fiber sliding are nonlocal dissipative mechanisms and to deal with this nonlocality, the virtual-power-based QC method is equipped with a mixed formulation in which the kinematic variables as well as the internal history variables are interpolated. Previously defined summation rules can still be used to sample the governing equations in this QC framework. Illustrative examples are presented. [less ▲] Detailed reference viewed: 126 (8 UL)Quasicontinuum-based multiscale approaches for plate-like beam lattices experiencing in-plane and out-of-plane deformation Beex, Lars ; ; et al in Computer Methods in Applied Mechanics & Engineering (2014), 279 The quasicontinuum (QC) method is a multiscale approach that aims to reduce the computational cost of discrete lattice computations. The method incorporates small-scale local lattice phenomena (e.g. a ... [more ▼] The quasicontinuum (QC) method is a multiscale approach that aims to reduce the computational cost of discrete lattice computations. The method incorporates small-scale local lattice phenomena (e.g. a single lattice defect) in macroscale simulations. Since the method works directly and only on the beam lattice, QC frameworks do not require the construction and calibration of an accompanying continuum model (e.g. a cosserat/micropolar description). Furthermore, no coupling procedures are required between the regions of interest in which the beam lattice is fully resolved and coarse domains in which the lattice is effectively homogenized. Hence, the method is relatively straightforward to implement and calibrate. In this contribution, four variants of the QC method are investigated for their use for planar beam lattices which can also experience out-of-plane deformation. The different frameworks are compared to the direct lattice computations for three truly multiscale test cases in which a single lattice defect is present in an otherwise perfectly regular beam lattice. [less ▲] Detailed reference viewed: 372 (13 UL)Nitsche’s method method for mixed dimensional analysis: conforming and non-conforming continuum-beam and continuum-plate coupling ; ; et al in Computer Methods in Applied Mechanics & Engineering (2014) A Nitche’s method is presented to couple different mechanical models. They include coupling of a solid and a beam and of a solid and a plate. Both conforming and non-conforming formulations are presented ... [more ▼] A Nitche’s method is presented to couple different mechanical models. They include coupling of a solid and a beam and of a solid and a plate. Both conforming and non-conforming formulations are presented. In a non-conforming formulation, the structure domain is overlapped by a refined solid model which is needed to either get more accuracy or to capture highly nonlinear events. Applications can be found in multi-dimensional analyses in which parts of a structure are modeled with solid elements and others are modeled using a coarser model with beam and/or plate elements. Discretisations are performed using both standard Lagrange elements and high order NURBS (Non Uniform Rational Bsplines) based isogeometric elements. We present various examples to demonstrate the performance of the method. [less ▲] Detailed reference viewed: 263 (17 UL)Explicit finite deformation analysis of isogeometric membranes ; ; et al in Computer Methods in Applied Mechanics & Engineering (2014) NURBS-based isogeometric analysis was first extended to thin shell/membrane structures which allows for finite membrane stretching as well as large deflection and bending strain. The assumed non-linear ... [more ▼] NURBS-based isogeometric analysis was first extended to thin shell/membrane structures which allows for finite membrane stretching as well as large deflection and bending strain. The assumed non-linear kinematics employs the Kirchhoff-Love shell theory to describe the mechanical behaviour of thin to ultrathin structures. The displacement fields are interpolated from the displacements of control points only, and no rotational degrees of freedom are used at control points. Due to the high order Ck (k ≥ 1) continuity of NURBS shape functions the Kirchhoff-Love theory can be seamlessly implemented. An explicit time integration scheme is used to compute the transient response of membrane structures to time-domain excitations, and a dynamic relaxation method is employed to obtain steady-state solutions. The versatility and good performance of the present formulation is demonstrated with the aid of a number of test cases, including a square membrane strip under static pressure, the inflation of a spherical shell under internal pressure, the inflation of a square airbag and the inflation of a rubber balloon. The mechanical contribution of the bending stiffness is also evaluated. [less ▲] Detailed reference viewed: 561 (6 UL)An adaptive singular ES-FEM for mechanics problems with singular field of arbitrary order ; ; Bordas, Stéphane et al in Computer Methods in Applied Mechanics & Engineering (2013), 253 This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded ... [more ▼] This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded linear triangular (T3) elements and a special layer of five-noded singular triangular elements (sT5) connected to the singular-point of the stress field. The sT5 element has an additional node on each of the two edges connected to the singular-point. It allows us to represent simple and efficient enrichment with desired terms for the displacement field near the singular-point with the satisfaction of partition-of-unity property. The stiffness matrix of the discretized system is then obtained using the assumed displacement values (not the derivatives) over smoothing domains associated with the edges of elements. An adaptive procedure for the sES-FEM is proposed to enhance the quality of the solution with minimized number of nodes. Several numerical examples are provided to validate the reliability of the present sES-FEM method. © 2012 Elsevier B.V. [less ▲] Detailed reference viewed: 331 (5 UL)A partitioned model order reduction approach to rationalise computational expenses in nonlinear fracture mechanics ; ; et al in Computer Methods in Applied Mechanics & Engineering (2013), 256 We propose in this paper a reduced order modelling technique based on domain partitioning for parametric problems of fracture. We show that coupling domain decomposition and projection-based model order ... [more ▼] We propose in this paper a reduced order modelling technique based on domain partitioning for parametric problems of fracture. We show that coupling domain decomposition and projection-based model order reduction permits to focus the numerical effort where it is most needed: around the zones where damage propagates. No a priori knowledge of the damage pattern is required, the extraction of the corresponding spatial regions being based solely on algebra. The efficiency of the proposed approach is demonstrated numerically with an example relevant to engineering fracture. © 2012 Elsevier B.V. [less ▲] Detailed reference viewed: 306 (12 UL)Isogeometric boundary element analysis using unstructured T-splines ; ; et al in Computer Methods in Applied Mechanics & Engineering (2013), 254 We couple collocated isogeometric boundary element methods and unstructured analysis-suitable T-spline surfaces for linear elastostatic problems. We extend the definition of analysis-suitable T-splines to ... [more ▼] We couple collocated isogeometric boundary element methods and unstructured analysis-suitable T-spline surfaces for linear elastostatic problems. We extend the definition of analysis-suitable T-splines to encompass unstructured control grids (unstructured meshes) and develop basis functions which are smooth (rational) polynomials defined in terms of the Bézier extraction framework and which pass standard patch tests. We then develop a collocation procedure which correctly accounts for sharp edges and corners, extraordinary points, and T-junctions. This approach is applied to several three-dimensional problems, including a real-world T-spline model of a propeller. We believe this work clearly illustrates the power of combining new analysis-suitable computer aided design technologies with established analysis methodologies, in this case, the boundary element method. © 2012 Elsevier B.V. [less ▲] Detailed reference viewed: 365 (13 UL)Extended finite element method with edge-based strain smoothing (ESm-XFEM) for linear elastic crack growth ; ; Bordas, Stéphane et al in Computer Methods in Applied Mechanics & Engineering (2012), 209-212 This paper presents a strain smoothing procedure for the extended finite element method (XFEM). The resulting "edge-based" smoothed extended finite element method (ESm-XFEM) is tailored to linear elastic ... [more ▼] This paper presents a strain smoothing procedure for the extended finite element method (XFEM). The resulting "edge-based" smoothed extended finite element method (ESm-XFEM) is tailored to linear elastic fracture mechanics and, in this context, to outperform the standard XFEM. In the XFEM, the displacement-based approximation is enriched by the Heaviside and asymptotic crack tip functions using the framework of partition of unity. This eliminates the need for the mesh alignment with the crack and re-meshing, as the crack evolves. Edge-based smoothing (ES) relies on a generalized smoothing operation over smoothing domains associated with edges of simplex meshes, and produces a softening effect leading to a close-to-exact stiffness, "super-convergence" and "ultra-accurate" solutions. The present method takes advantage of both the ES-FEM and the XFEM. Thanks to the use of strain smoothing, the subdivision of elements intersected by discontinuities and of integrating the (singular) derivatives of the approximation functions is suppressed via transforming interior integration into boundary integration. Numerical examples show that the proposed method improves significantly the accuracy of stress intensity factors and achieves a near optimal convergence rate in the energy norm even without geometrical enrichment or blending correction. [less ▲] Detailed reference viewed: 90 (1 UL) |
||