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 linear smoothed higher-order CS-FEM for the analysis of notched laminated composites ; ; et al in Engineering Analysis with Boundary Elements (2017), 85 Higher-order elements with highly accurate solutions are attractive for stress analysis and stress concentration problems. However, the distorted eight-node serendipity quadrilateral element is known to ... [more ▼] Higher-order elements with highly accurate solutions are attractive for stress analysis and stress concentration problems. However, the distorted eight-node serendipity quadrilateral element is known to yield inaccurate re- sults and sub-optimal convergence rate. In this paper, we present a higher order CS-FEM to alleviate the effect of distorted mesh and guarantee the quality of solutions by employing a linear smoothing technique over eight-node quadratic serendipity elements. The modified. strain matrix is computed by the divergence theorem between the nodal shape functions and their derivatives using Taylor’s expansion of the weak form. The proposed method eliminates the need for isoparametric mapping and numerical studies demonstrate that the proposed method is insensitive to mesh distortion. The improved accuracy and superior convergence rates are numerically demon- strated with a few benchmark problems. The analysis of the stress concentration around cutouts also proves that the present method has good performance for the laminated composites. [less ▲] Detailed reference viewed: 129 (1 UL)A linear smoothed quadratic finite element for the analysis of laminated composite Reissner–Mindlin plates ; ; et al in Composite Structures (2017), 180 It is well known that the high-order elements have significantly improved the accuracy of solutions in the traditional finite element analysis, but the performance of high-order elements is restricted by ... [more ▼] It is well known that the high-order elements have significantly improved the accuracy of solutions in the traditional finite element analysis, but the performance of high-order elements is restricted by the shear-locking and distorted meshes for the plate problems. In this paper, a linear smoothed eight-node Reissner-Mindlin plate element (Q8 plate element) based on the first order shear deformation theory is developed for the static and free vibration analysis of laminated composite plates, the computation of the interior derivatives of shape function and isoparametric mapping can be removed. The strain matrices are modified with a linear smoothing technique by using the divergence theorem between the nodal shape functions and their derivatives in Taylor’s expansion. Moreover, the first order Taylor’s expansion is also employed for the construction of stiffness matrix to satisfy the linear strain distribution. Several numerical examples indicate that the novel Q8 plate element has good performance to alleviate the shear-locking phenomenon and improve the quality of the solutions with distorted meshes. [less ▲] Detailed reference viewed: 133 (2 UL)Error-controlled adaptive extended finite element method for 3D linear elastic crack propagation ; ; et al in Computer Methods in Applied Mechanics and 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: 128 (4 UL)Computational Sciences Year 2016 Activity Report Bordas, Stéphane Report (2016) Born from a bottom-up initiative of Mathematics, Computer Science, Physics and Computational Engineering, Computational Sciences (CoSc) have contributed to create at UL a positive and symbiotic research ... [more ▼] Born from a bottom-up initiative of Mathematics, Computer Science, Physics and Computational Engineering, Computational Sciences (CoSc) have contributed to create at UL a positive and symbiotic research environment relying on a strong fundamental scientific research core. CoSc will continue to rationalize research efforts across a range of strategic innovation domains by centralizing research and development tools and building upon the existing strengths of the Luxembourgish research and socio-economic landscape. [less ▲] Detailed reference viewed: 106 (4 UL)Image to analysis pipeline: single and double balloons kyphoplasty Baroli, Davide ; Hauseux, Paul ; Hale, Jack et al Poster (2016, December 12) In this work, we present a semi-automatic pipeline from image to simulation of a patient fractured vertebra after the kyphoplastic augmentation with two balloons. In this procedure, the CT-scan medical ... [more ▼] In this work, we present a semi-automatic pipeline from image to simulation of a patient fractured vertebra after the kyphoplastic augmentation with two balloons. In this procedure, the CT-scan medical image are pre-processed using open-source software Slice3D for segmentation and 3D reconstruction operation. Then, using geometric processing the 3D surface geometry is enhanced to avoid degenerate element and trigging phenomena on vertebra and cement area. We perform a finite element analysis to evaluate the risk of subsequent vertebral fracture. Finally using Monte-Carlo technique, we assess the propagation of uncertainty of material parameter on the evaluation of this risk. Based on the developed semi-automatic pipelines, it is possible to perform a patient-specific simulation that assesses the successful of kyphoplasty operation. [less ▲] Detailed reference viewed: 230 (35 UL)Numerical methods for fracture/cutting of heterogeneous materials Sutula, Danas ; Agathos, Konstantinos ; Ziaei Rad, Vahid et al Presentation (2016, December) Detailed reference viewed: 204 (15 UL)Shape Optimization Directly from CAD: an Isogeometric Boundary Element Approach Using T-splines ; ; Bordas, Stéphane Report (2016) Detailed reference viewed: 326 (5 UL)Multi-scale modelling of fracture Bordas, Stéphane ; ; et al Speeches/Talks (2016) We present recent models on complexity reduction for computational fracture mechanics Detailed reference viewed: 193 (7 UL)Simulating topological changes in real time for surgical assistance Bordas, Stéphane ; ; et al Speeches/Talks (2016) Detailed reference viewed: 457 (38 UL)Weakening the tight coupling between geometry and simulation in isogeometric analysis Tomar, Satyendra ; ; et al Presentation (2016, June 07) In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry ... [more ▼] In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is employed for the numerical solution of the problem over the domain. However, there are situations where this tight integration is a bane rather than a boon. Such situations arise where, e.g., (1) the geometry of the domain is simple enough to be represented by low order NURBS, whereas the unknown (exact) solution of the problem is sufficiently regular, and thus, the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry, (2) the constraint of using the same space for the geometry and the numerical solution is particularly undesirable, such as in the shape and topology optimization, and (3) the solution of the problem has low regularity but for the curved boundary of the domain one can employ higher order NURBS. Therefore, we propose to weaken this constraint. An extensive study of patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution will be discussed. It will be shown, with concrete reasoning, that why patch test fails in certain cases, and that those cases should be avoided in practice. Thereafter, selective numerical examples will be presented to address some of the above-mentioned situations, and it will be shown that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces, and thus, improved accuracy of the numerical solution. [less ▲] Detailed reference viewed: 179 (9 UL)Weakening the tight coupling between geometry and simulation in isogeometric analysis Bordas, Stéphane ; Tomar, Satyendra ; et al Scientific Conference (2016, June 05) In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry ... [more ▼] In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is employed for the numerical solution of the problem over the domain. However, there are situations where this tight integration is a bane rather than a boon. Such situations arise where, e.g., (1) the geometry of the domain is simple enough to be represented by low order NURBS, whereas the unknown (exact) solution of the problem is sufficiently regular, and thus, the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry, (2) the constraint of using the same space for the geometry and the numerical solution is particularly undesirable, such as in the shape and topology optimization, and (3) the solution of the problem has low regularity but for the curved boundary of the domain one can employ higher order NURBS. Therefore, we propose to weaken this constraint. An extensive study of patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution will be discussed. It will be shown, with concrete reasoning, that why patch test fails in certain cases, and that those cases should be avoided in practice. Thereafter, selective numerical examples will be presented to address some of the above-mentioned situations, and it will be shown that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces, and thus, improved accuracy of the numerical solution. Powered by [less ▲] Detailed reference viewed: 153 (5 UL)Generalizing the isogeometric concept: weakening the tight coupling between geometry and simulation in IGA Tomar, Satyendra ; ; et al Presentation (2016, June 02) In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the ... [more ▼] In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is also employed for the numerical solution of the problem over the domain. However, in certain situations, such as, when the geometry of the domain can be represented by low order NURBS but the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry; or in the shape and topology optimization where the constraint of using the same space for the geometry and the numerical solution is not favorable, this tight coupling is disadvantageous. Therefore, we study the effect of decoupling the spaces for the geometry representation and the numerical solution, though still using the prevalent functions in CAD/CAGD. To begin with, we perform the patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution. This shows that certain cases, perhaps intuitive, should be avoided in practice because patch test fails. The above-mentioned situations are further explored with some numerical examples, which shows that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces. [1] J. Cottrell, T.J.R. Hughes, and Y. Bazilevs. Isogeometric Analysis: Toward Integration of CAD and FEA, volume 80. Wiley, Chichester, 2009. [2] T.J.R. Hughes, J. Cottrell, and Y. Bazilevs. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Computer Methods in Applied Mechanics and Engineering, 194:4135–4195, 2005. [less ▲] Detailed reference viewed: 184 (11 UL)Virtual-power-based quasicontinuum methods for discrete dissipative materials Beex, Lars ; Bordas, Stéphane Scientific Conference (2016, June) Detailed reference viewed: 87 (2 UL)Numerical study of magnetic particles concentration in biofluid (blood) under the influence of high gradient magnetic field in microchannel ; Bourantas, Georgios ; et al Scientific Conference (2016, June) A meshless numerical scheme [1] is developed in order to simulate the magnetically mediated separation of biological mixture used in lab-on-chip devices as solid carriers for capturing, transporting and ... [more ▼] A meshless numerical scheme [1] is developed in order to simulate the magnetically mediated separation of biological mixture used in lab-on-chip devices as solid carriers for capturing, transporting and detecting biological magnetic labeled entities [2], as well as for drug delivering, magnetic hyperthermia treatment, magnetic resonance imaging, magnetofection, etc. A modified one-way particle-fluid coupling analysis is considered to model the interaction of the base fluid of the mixture with the distributed particles motion. In details, bulk flow influences particle motion (through a simplified Stokes drag relation), while it is strongly dependent on particle motion through (particle) concentration. Due to the imposed magnetic field stagnation regions are developed, leading to the accumulation of the magnetic labeled species and finally to their collection from the heterogeneous mixture. The role of (i) the intensity of magnetic field and its gradient, (ii) the position of magnetic field, (iii) the magnetic susceptibility of magnetic particles, (iv) the volume concentration of magnetic particles (nanoparticles) and their size, (v) the flow velocity in the magnetic- fluidic interactions and interplay between the magnetophoretic mass transfer and molecular diffusion are thoroughly investigated. Both Newtonian and non-Newtonian blood flow models are considered, along with different expressions for the concentration and numerical results are presented for a wide range of physical parameters (Hartmann number (Ha), Reynolds number (Re)). A comprehensive study investigates their impact on the biomagnetic separation. For verification purposes, the numerical results obtained by the proposed meshless scheme were compared with established numerical results from the literature, being in excellent agreement. [less ▲] Detailed reference viewed: 386 (11 UL)Linear smoothing over arbitrary polytopes for compressible and nearly incompressible linear elasticity ; Tomar, Satyendra ; Bordas, Stéphane et al Scientific Conference (2016, June) We present a displacement based approach over arbitrary polytopes for compressible and nearly incompressible linear elastic solids. In this approach, a volume-averaged nodal projection operator is ... [more ▼] We present a displacement based approach over arbitrary polytopes for compressible and nearly incompressible linear elastic solids. 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 formulation uses the usual Wachspress interpolants over arbitrary polytopes and the stability of the method is ensured by the addition of bubble like functions. The smoothed strains are evaluated based on the linear smoothing procedure. This further softens the bilinear form allowing the procedure to search for a solution satisfying the divergence- free condition. The divergence-free condition of the proposed approach is verified through systematic numerical study. The formulation delivers optimal convergence rates in the energy and L2-norms. Inf-sup tests are presented to demonstrated the stability of the formulation. [less ▲] Detailed reference viewed: 239 (3 UL)Well Conditioned and Optimally Convergent Extended Finite Elements and Vector Level Sets for Three-Dimensional Crack Propagation Agathos, Konstantinos ; ; et al Scientific Conference (2016, June) A three-dimensional (3D) version of the vector level set method [1] is combined to a well conditioned and optimally convergent XFEM variant in order to deal with non-planar three dimensional crack ... [more ▼] A three-dimensional (3D) version of the vector level set method [1] is combined to a well conditioned and optimally convergent XFEM variant in order to deal with non-planar three dimensional crack propagation problems. The proposed computational fracture method achieves optimal convergence rates by using tip enriched elements in a fixed volume around the crack front (geometrical enrichment) while keeping conditioning of the resulting system matrices in acceptable levels. Conditioning is controlled by using a three dimensional extension of the degree of freedom gathering technique [2]. Moreover, blending errors are minimized and conditioning is further improved by employing weight function blending and enrichment function shifting [3,4]. As far as crack representation is concerned, crack surfaces are represented by linear quadrilateral elements and the corresponding crack fronts by ordered series of linear segments. Level set values are obtained by projecting points at the crack surface and front respectively. Different criteria are employed in order to assess the quality of the crack representation. References [1] Ventura G., Budyn E. and Belytschko T. Vector level sets for description of propagating cracks in finite elements. Int. J. Numer. Meth. Engng. 58:1571-1592 (2003). [2] Laborde P., Pommier J., Renard Y. and Salaün M. High-order extended finite element method for cracked domains. Int. J. Numer. Meth. Engng. 64:354-381 (2005). [3] Fries T.P. A corrected XFEM approximation without problems in blending elements. Int. J. Numer. Meth. Engng. 75:503-532 (2008). [4] Ventura G., Gracie R. and Belytschko T. Fast integration and weight function blending in the extended finite element method. Int. J. Numer. Meth. Engng. 77:1-29 (2009). [less ▲] Detailed reference viewed: 182 (15 UL)Generalizing the isogeometric concept: weakening the tight coupling between geometry and simulation in IGA Bordas, Stéphane ; Tomar, Satyendra ; et al Scientific Conference (2016, May 30) In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the ... [more ▼] In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is also employed for the numerical solution of the problem over the domain. However, in certain situations, such as, when the geometry of the domain can be represented by low order NURBS but the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry; or in the shape and topology optimization where the constraint of using the same space for the geometry and the numerical solution is not favorable, this tight coupling is disadvantageous. Therefore, we study the effect of decoupling the spaces for the geometry representation and the numerical solution, though still using the prevalent functions in CAD/CAGD. To begin with, we perform the patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution. This shows that certain cases, perhaps intuitive, should be avoided in practice because patch test fails. The above-mentioned situations are further explored with some numerical examples, which shows that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces. [less ▲] Detailed reference viewed: 150 (3 UL)Computational mechanics of interfaces Bordas, Stéphane Presentation (2016, May 22) The course will present an overview of recent developments, which will enable students to make informed choices in terms of discretization and model selection in solving numerical problems in mechanics ... [more ▼] The course will present an overview of recent developments, which will enable students to make informed choices in terms of discretization and model selection in solving numerical problems in mechanics. We will cover discretization strategies ranging from the standard finite element method, the smoothed finite element method, the extended finite element method, polygonal and virtual element methods, meshfree methods. The applications range between fracture of heterogeneous materials and biomedical simulations. The intended learning outcomes of the course are such that the students will be: - able to critically assess discretization schemes in mechanics - able to implement simple error estimators for finite element methods - familiar with basic multi-scale methods for fracture and their limitations - able to follow basic literature in model error and model selection, in particular based on Bayesian inference Course participants will learn these topics through lectures and hands-on numerical experiments. [less ▲] Detailed reference viewed: 414 (16 UL)3D Crack Detection Using an XFEM Variant and Global Optimization Algorithms Agathos, Konstantinos ; ; Bordas, Stéphane Scientific Conference (2016, May) Detailed reference viewed: 241 (20 UL)Error estimation and space-time adaptivity for the isogeometric analysis of transient structural dynamics ; ; Bordas, Stéphane et al Scientific Conference (2016, April 01) This paper presents a new adaptive scheme for the error-controlled simulation of transient dynamics problem. We rely on spline bases for the higher-order spatial description of our kinematic fields. Local ... [more ▼] This paper presents a new adaptive scheme for the error-controlled simulation of transient dynamics problem. We rely on spline bases for the higher-order spatial description of our kinematic fields. Local adaptivity is performed by employing a hierarchical T-mesh technology, in combination with geometry independent field approximation. The Newmark algorithm is chosen to solve the semidiscrete equation of motion. We will present some simple local error estimates to drive the adaptivity, and show how we can ensure that the mechanical energy of conservative systems is preserved during the refinement process. [less ▲] Detailed reference viewed: 128 (4 UL) |
||