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 smoothed finite element method for plate analysis ; ; Bordas, Stéphane et al in Computer Methods in Applied Mechanics and Engineering (2008), 197(13-16), 1184-1203 A quadrilateral element with smoothed curvatures for Mindlin-Reissner plates is proposed. The curvature at each point is obtained by a non-local approximation via a smoothing function. The bending ... [more ▼] A quadrilateral element with smoothed curvatures for Mindlin-Reissner plates is proposed. The curvature at each point is obtained by a non-local approximation via a smoothing function. The bending stiffness matrix is calculated by a boundary integral along the boundaries of the smoothing elements (smoothing cells). Numerical results show that the proposed element is robust, computational inexpensive and simultaneously very accurate and free of locking, even for very thin plates. The most promising feature of our elements is their insensitivity to mesh distortion. [less ▲] Detailed reference viewed: 113 (2 UL)Three-dimensional non-linear fracture mechanics by enriched meshfree methods without asymptotic enrichment Bordas, Stéphane ; ; in Proceedings of the IUTAM Symposium on Discretization Methods for Evolving Discontinuities (2007) This paper presents a three-dimensional, extrinsically enriched meshfree method for initiation, growth and coalescence of an arbitrary number of cracks in non-linear solids including large deformations ... [more ▼] This paper presents a three-dimensional, extrinsically enriched meshfree method for initiation, growth and coalescence of an arbitrary number of cracks in non-linear solids including large deformations, for statics and dynamics. The novelty of the methodology fashioned in this work is that only an extrinsic discontinuous enrichment and no near-tip/asymptotic enrichment is required. Instead, a Lagrange multiplier field is added along the crack front to close the crack along the front. This decreases the computational cost and removes difficulties involved with a branch enrichment. Numerical examples treated include the pull-out of a reinforcement bar from a concrete block, and a Taylor bar impact with very large deformation and fragmentation. The results are compared to experimental results, and other simulations from the literature, which shows the robustness and accuracy of the method. © 2007 Springer. [less ▲] Detailed reference viewed: 175 (0 UL)Derivative recovery and a posteriori error estimate for extended finite elements Bordas, Stéphane ; in Computer Methods in Applied Mechanics and Engineering (2007), 196(35-36), 3381-3399 This paper is the first attempt at error estimation for extended finite elements. The goal of this work is to devise a simple and effective local a posteriori error estimate for partition of unity ... [more ▼] This paper is the first attempt at error estimation for extended finite elements. The goal of this work is to devise a simple and effective local a posteriori error estimate for partition of unity enriched finite element methods such as the extended finite element method (XFEM). In each element, the local estimator is the L2 norm of the difference between the raw XFEM strain field and an enhanced strain field computed by extended moving least squares (XMLS) derivative recovery obtained from the raw nodal XFEM displacements. The XMLS construction is tailored to the nature of the solution. The technique is applied to linear elastic fracture mechanics, in which near-tip asymptotic functions are added to the MLS basis. The XMLS shape functions are constructed from weight functions following the diffraction criterion to represent the discontinuity. The result is a very smooth enhanced strain solution including the singularity at the crack tip. Results are shown for two- and three-dimensional linear elastic fracture mechanics problems in mode I and mixed mode. The effectivity index of the estimator is close to 1 and improves upon mesh refinement for the studied near-tip problem. It is also shown that for the linear elastic fracture mechanics problems treated, the proposed estimator outperforms one of the superconvergent patch recovery technique of Zienkiewicz and Zhu, which is only C0. Parametric studies of the general performance of the estimator are also carried out. © 2007 Elsevier B.V. All rights reserved. [less ▲] Detailed reference viewed: 189 (5 UL)Architecture tradeoffs of integrating a mesh generator to partition of unity enriched object-oriented finite element software ; ; et al in European Journal of Computational Mechanics (2007), 16(2), 237-258 We explore the tradeoffs of using an internal mesher in a XFEM code. We show that it allows an efficient enrichement detection scheme, while retaining the ability to have welladapted meshes. We provide ... [more ▼] We explore the tradeoffs of using an internal mesher in a XFEM code. We show that it allows an efficient enrichement detection scheme, while retaining the ability to have welladapted meshes. We provide benchmarks highlighting the considerable gains which can be expected from a well designed architecture. The efficiency of the proposed algorithm is shown by solving fracture mechanics problems of densely micro-cracked bodies including adaptive mesh refinement. [less ▲] Detailed reference viewed: 109 (0 UL)A three-dimensional meshfree method for continuous multiple-crack initiation, propagation and junction in statics and dynamics ; Bordas, Stéphane ; in Computational Mechanics (2007), 40(3), 473-495 This paper proposes a three-dimensional meshfree method for arbitrary crack initiation and propagation that ensures crack path continuity for non-linear material models and cohesive laws. The method is ... [more ▼] This paper proposes a three-dimensional meshfree method for arbitrary crack initiation and propagation that ensures crack path continuity for non-linear material models and cohesive laws. The method is based on a local partition of unity. An extrinsic enrichment of the meshfree shape functions is used with discontinuous and near-front branch functions to close the crack front and improve accuracy. The crack is hereby modeled as a jump in the displacement field. The initiation and propagation of a crack is determined by the loss of hyperbolicity or the loss of material stability criterion. The method is applied to several static, quasi-static and dynamic crack problems. The numerical results very precisely replicate available experimental and analytical results. [less ▲] Detailed reference viewed: 98 (0 UL)A simulation-based design paradigm for complex cast components Bordas, Stéphane ; ; et al in Engineering with Computers (2007), 23(1), 25-37 This paper describes and exercises a new design paradigm for cast components. The methodology integrates foundry process simulation, non-destructive evaluation (NDE), stress analysis and damage tolerance ... [more ▼] This paper describes and exercises a new design paradigm for cast components. The methodology integrates foundry process simulation, non-destructive evaluation (NDE), stress analysis and damage tolerance simulations into the design process. Foundry process simulation is used to predict an array of porosity-related anomalies. The probability of detection of these anomalies is investigated with a radiographic inspection simulation tool (XRSIM). The likelihood that the predicted array of anomalies will lead to a failure is determined by a fatigue crack growth simulation based on the extended finite element method and therefore does not require meshing nor remeshing as the cracks grow. With this approach, the casting modeling provides initial anomaly information, the stress analysis provides a value for the critical size of an anomaly and the NDE assessment provides a detectability measure. The combination of these tools allows for accept/reject criteria to be determined at the early design stage and enables damage tolerant design philosophies. The methodology is applied to the design of a cast monolithic door used on the Boeing 757 aircraft. [less ▲] Detailed reference viewed: 100 (0 UL)Mechanical failure in microstructural heterogeneous materials Bordas, Stéphane ; ; in Lecture Notes in Computer Science (2007), 4310 LNCS Various heterogeneous materials with multiple scales and multiple phases in the microstructure have been produced in the recent years. We consider a mechanical failure due to the initiation and ... [more ▼] Various heterogeneous materials with multiple scales and multiple phases in the microstructure have been produced in the recent years. We consider a mechanical failure due to the initiation and propagation of cracks in places of high pore density in the microstructures. A multi-scale method based on the asymptotic homogenization theory together with the mesh superposition method (s-version of FEM) is presented for modeling of cracks. The homogenization approach is used on the global domain excluding the vicinity of the crack where the periodicity of the microstructures is lost and this approach fails. The multiple scale method relies on efficient combination of both macroscopic and microscopic models. The mesh superposition method uses two independent (global and local) finite element meshes and the concept of superposing the local mesh onto the global continuous mesh in such a way that both meshes not necessarily coincide. The homogenized material model is considered on the global mesh while the crack is analyzed in the local domain (patch) which allows to have an arbitrary geometry with respect to the underlying global finite elements. Numerical experiments for biomorphic cellular ceramics with porous microstructures produced from natural wood are presented. [less ▲] Detailed reference viewed: 102 (2 UL)Linear smoothing over arbitrary polytopes ; ; et al Scientific Conference (n.d.) The conventional constant strain smoothing technique yields less accurate solutions that other techniques such as the conventional polygonal finite element method [1, 2]. In this work, we propose a linear ... [more ▼] The conventional constant strain smoothing technique yields less accurate solutions that other techniques such as the conventional polygonal finite element method [1, 2]. In this work, we propose a linear strain smoothing scheme that improves the accuracy of linear and quadratic approximations over convex poly- topes. The method relies on sub-division of the polytope into simplical subcells; however instead of using a constant smoothing function, we employ a linear smoothing function over each subcell. This gives a new definition for the strain to compute the stiffness matrix. The convergence properties and accuracy of the proposed scheme are discussed by solving few benchmark problems. Numerical results show that the proposed linear strain smoothing scheme makes the approximation based on polytopes able to deliver the optimal convergence rate as in traditional quadrilateral and hexahedral finite elements. The accuracy is also improved, and all the methods tested pass the patch test to machine precision. [less ▲] Detailed reference viewed: 139 (0 UL)Real-time error controlled adaptive mesh refinement in surgical simulation: Application to needle insertion simulation ; Tomar, Satyendra ; et al in IEEE Transactions on Biomedical Engineering (n.d.) This paper presents the first real-time discretisation-error-driven adaptive finite element approach for corotational elasticity problems involving strain localisation. We propose a hexahedron-based ... [more ▼] This paper presents the first real-time discretisation-error-driven adaptive finite element approach for corotational elasticity problems involving strain localisation. We propose a hexahedron-based finite element method combined with local oct-tree $h$-refinement, driven by a posteriori error estimation, for simulating soft tissue deformation. This enables to control the local error and global error level in the mechanical fields during the simulation. The local error level is used to refine the mesh only where it is needed, while maintaining a coarser mesh elsewhere. We investigate the convergence of the algorithm on academic examples, and demonstrate its practical usability on a percutaneous procedure involving needle insertion in a liver. For the latter case, we compare the force displacement curves obtained from the proposed adaptive algorithm with that obtained from a uniform refinement approach. [less ▲] Detailed reference viewed: 664 (63 UL)Linear smoothed extended finite element method ; ; et al E-print/Working paper (n.d.) The extended finite element method (XFEM) was introduced in 1999 to treat problems involving discontinuities with no or minimal remeshing through appropriate enrichment functions. This enables elements to ... [more ▼] The extended finite element method (XFEM) was introduced in 1999 to treat problems involving discontinuities with no or minimal remeshing through appropriate enrichment functions. This enables elements to be split by a discontinuity, strong or weak and hence requires the integration of discontinuous functions or functions with discontinuous derivatives over elementary volumes. Moreover, in the case of open surfaces and singularities, special, usually non-polynomial functions must also be integrated.A variety of approaches have been proposed to facilitate these special types of numerical integration, which have been shown to have a large impact on the accuracy and convergence of the numerical solution. The smoothed extended finite element method (SmXFEM) [1], for example, makes numerical integration elegant and simple by transforming volume integrals into surface integrals. However, it was reported in [1, 2] that the strain smoothing is inaccurate when non-polynomial functions are in the basis. This is due to the constant smoothing function used over the smoothing domains which destroys the effect of the singularity. In this paper, we investigate the benefits of a recently developed Linear smoothing procedure [3] which provides better approximation to higher order polynomial fields in the basis. Some benchmark problems in the context of linear elastic fracture mechanics (LEFM) are solved to compare the standard XFEM, the constant-smoothed XFEM (Sm-XFEM) and the linear-smoothed XFEM (LSm-XFEM). We observe that the convergence rates of all three methods are the same. The stress intensity factors (SIFs) computed through the proposed LSm-XFEM are however more accurate than that obtained through Sm-XFEM. To conclude, compared to the conventional XFEM, the same order of accuracy is achieved at a relatively low computational effort. [less ▲] Detailed reference viewed: 134 (3 UL)Minimum energy multiple crack propagation. Part II: Discrete Solution with XFEM. Sutula, Danas ; Bordas, Stéphane in Engineering Fracture Mechanics (n.d.) The three-part paper deals with energy-minimal multiple crack propagation in a linear elastic solid under quasi-static conditions. The principle of minimum total energy, i.e. the sum of the potential and ... [more ▼] The three-part paper deals with energy-minimal multiple crack propagation in a linear elastic solid under quasi-static conditions. The principle of minimum total energy, i.e. the sum of the potential and fracture energies, which stems directly from the Griffith's theory of cracks, is applied to the problem of arbitrary crack growth in 2D. The proposed formulation enables minimisation of the total energy of the mechanical system with respect to the crack extension directions and crack extension lengths to solve for the evolution of the mechanical system over time. The three parts focus, in turn, on (I) the theory of multiple crack growth including competing cracks, (II) the discrete solution by the extended finite element method using the minimum-energy formulation, and (III) the aspects of computer implementation within the Matlab programming language. This Part-II of our three-part paper examines three discrete solution methods for solving fracture mechanics problems based on the principle of minimum total energy. The discrete solution approach is chosen based on the stability property of the fracture configuration at hand. The first method is based on external load-control. It is suitable for stable crack growth and stable fracture configurations. The second method is based on fractured area-control. This method is applicable to stable or unstable fracture growth but it is required that the fracture front be stable. The third solution method is based on a gradient-descent approach. This approach can be applied to arbitrary crack growth problems; however, the gradient-descent formulation cannot be guaranteed to yield the optimal solution in the case of competing crack growth and an unstable fracture front configuration. The main focus is on the gradient-descent solution approach within the framework of the extended finite element discretisation. Although a viable solution method is finally proposed for resolving competing crack growth in the case of an unstable fracture front configuration, the method is not implemented within the present XFEM code but rather exists as a separate proof-of-concept algorithm that is tested against several fabricated benchmark problems. The open-source Matlab code, documentation and example cases are included as supplementary material. [less ▲] Detailed reference viewed: 1699 (135 UL)Minimum energy multiple crack propagation Part I: Theory. Sutula, Danas ; Bordas, Stéphane in Engineering Fracture Mechanics (n.d.) The three-part paper deals with energy-minimal multiple crack propagation in a linear elastic solid under quasi-static conditions. The principle of minimum total energy, i.e. the sum of the potential and ... [more ▼] The three-part paper deals with energy-minimal multiple crack propagation in a linear elastic solid under quasi-static conditions. The principle of minimum total energy, i.e. the sum of the potential and fracture energies, which stems directly from the Griffith's theory of cracks, is applied to the problem of arbitrary crack growth in 2D. The proposed formulation enables minimisation of the total energy of the mechanical system with respect to the crack extension directions and crack extension lengths to solve for the evolution of the mechanical system over time. The three parts focus, in turn, on (I) the theory of multiple crack growth including competing cracks, (II) the discrete solution by the extended finite element method using the minimum-energy formulation, and (III) the aspects of computer implementation within the Matlab programming language. The key contributions of Part-I of this three-part paper are: (1) formulation of the total energy functional governing multiple crack behaviour, (2) three solution methods to the problem of competing crack growth for different fracture front stabilities (e.g. stable, unstable, or a partially stable configuration of crack tips), and (3) the minimum energy criterion for a set of crack tip extensions is posed as the criterion of vanishing rotational dissipation rates with respect to the rotations of the crack extensions. The formulation lends itself to a straightforward application within a discrete framework for determining the crack extension directions of multiple finite-length crack tip increments, which is tackled in Part-II, using the extended finite element method. In Part-III, we discuss various applications and benchmark problems. The open-source Matlab code, documentation, benchmark/example cases are included as supplementary material. [less ▲] Detailed reference viewed: 3544 (199 UL)Minimum energy multiple crack propagation. Part III: XFEM computer implementation and applications. Sutula, Danas ; Bordas, Stéphane in Engineering Fracture Mechanics (n.d.) The three-part paper deals with energy-minimal multiple crack propagation in a linear elastic solid under quasi-static conditions. The principle of minimum total energy, i.e. the sum of the potential and ... [more ▼] The three-part paper deals with energy-minimal multiple crack propagation in a linear elastic solid under quasi-static conditions. The principle of minimum total energy, i.e. the sum of the potential and fracture energies, which stems directly from the Griffith's theory of cracks, is applied to the problem of arbitrary crack growth in 2D. The proposed formulation enables minimisation of the total energy of the mechanical system with respect to the crack extension directions and crack extension lengths to solve for the evolution of the mechanical system over time. The three parts focus, in turn, on (I) the theory of multiple crack growth including competing cracks, (II) the discrete solution by the extended finite element method using the minimum-energy formulation, and (III) the aspects of computer implementation within the Matlab programming language. The key contributions of Part-III of the three-part paper are as follows: (1) implementation of XFEM in Matlab with emphasis on the design of the code to enable fast and efficient computational times of fracture problems involving multiple cracks and arbitrary crack intersections, (2) verification of the minimum energy criterion and comparison with the maximum tension criterion via multiple benchmark studies, and (3) we propose a numerical improvement to the crack growth direction criterion that gives significant improvements in accuracy and convergence rates of the fracture paths, especially on coarse meshes. The comparisons of the fracture paths obtained by the maximum tension (or maximum hoop-stress) criterion and the energy minimisation approach via a multitude of numerical case studies show that both criteria converge to virtually the same fracture solutions albeit from opposite directions. In other words, it is found that the converged fracture path lies in between those obtained by each criterion on coarser meshes. Thus, a modified crack growth direction criterion is proposed that assumes the average direction of the directions obtained by the maximum tension and the minimum energy criteria. The numerical results show significant improvements in accuracy (especially on coarse discretisations) and convergence rates of the fracture paths. Finally, the open-source Matlab code, documentation, benchmarks and example cases are included as supplementary material. [less ▲] Detailed reference viewed: 1602 (122 UL)A mass conservative Kalman filter algorithm for thermo-computational fluid dynamics ; Baroli, Davide ; et al in Materials (n.d.) Computational fluid-dynamics (CFD) is of wide relevance in engineering and science, due to its capability of simulating the three-dimensional flow at various scales. However, the suitability of a given ... [more ▼] Computational fluid-dynamics (CFD) is of wide relevance in engineering and science, due to its capability of simulating the three-dimensional flow at various scales. However, the suitability of a given model depends on the actual scenarios which are encountered in practice. This challenge of model suitability and calibration could be overcome by a dynamic integration of measured data into the simulation. This paradigm is known as data-driven assimilation (DDA). In this paper, the study is devoted to Kalman filtering, a Bayesian approach, applied to Reynolds-Averaged Navier-Stokes (RANS) equations for turbulent flow. The integration of the Kalman estimator into the PISO segregated scheme was recently investigated by (1). In this work, this approach is extended to the PIMPLE segregated method and to the ther- modynamic analysis of turbulent flow, with the addition of a sub-stepping procedure that ensures mass conservation at each time step and the com- patibility among the unknowns involved. The accuracy of the algorithm is verified with respect to the heated lid-driven cavity benchmark, incorporat- ing also temperature observations, comparing the augmented prediction of the Kalman filter with the CFD solution obtained on a very fine grid. [less ▲] Detailed reference viewed: 131 (6 UL)Bayesian inference for the stochastic identification of elastoplastic material parameters: Introduction, misconceptions and insights Rappel, Hussein ; Beex, Lars ; Hale, Jack et al E-print/Working paper (n.d.) We discuss Bayesian inference (BI) for the probabilistic identification of material parameters. This contribution aims to shed light on the use of BI for the identification of elastoplastic material ... [more ▼] We discuss Bayesian inference (BI) for the probabilistic identification of material parameters. This contribution aims to shed light on the use of BI for the identification of elastoplastic material parameters. For this purpose a single spring is considered, for which the stress-strain curves are artificially created. Besides offering a didactic introduction to BI, this paper proposes an approach to incorporate statistical errors both in the measured stresses, and in the measured strains. It is assumed that the uncertainty is only due to measurement errors and the material is homogeneous. Furthermore, a number of possible misconceptions on BI are highlighted based on the purely elastic case. [less ▲] Detailed reference viewed: 369 (104 UL)Linear smoothed polygonal and polyhedral finite elements ; ; Bordas, Stéphane et al E-print/Working paper (n.d.) It was observed in [1, 2] that the strain smoothing technique over higher order elements and arbitrary polytopes yields less accurate solutions than other techniques such as the conventional polygonal ... [more ▼] It was observed in [1, 2] that the strain smoothing technique over higher order elements and arbitrary polytopes yields less accurate solutions than other techniques such as the conventional polygonal finite element method. In this work, we propose a linear strain smoothing scheme that improves the accuracy of linear and quadratic approximations over convex polytopes. The main idea is to subdivide the polytope into simplicial subcells and use a linear smoothing function in each subcell to compute the strain. This new strain is then used in the computation of the stiffness matrix. The convergence properties and accuracy of the proposed scheme are discussed by solving few benchmark problems. Numerical results show that the proposed linear strain smoothing scheme makes the approximation based on polytopes to deliver improved accuracy and pass the patch test to machine precision. [less ▲] Detailed reference viewed: 464 (10 UL)A well-conditioned and optimally convergent XFEM for 3D linear elastic fracture ; ; Bordas, Stéphane et al in International Journal for Numerical Methods in Engineering (n.d.) A variation of the extended finite element method for 3D fracture mechanics is proposed. It utilizes global enrichment and point-wise as well as integral matching of displacements of the standard and ... [more ▼] A variation of the extended finite element method for 3D fracture mechanics is proposed. It utilizes global enrichment and point-wise as well as integral matching of displacements of the standard and enriched elements in order to achieve higher accuracy, optimal convergence rates and improved conditioning for two and three dimensional crack problems. A bespoke benchmark problem is introduced to determine the method's accuracy in the general 3D case where it is demonstrated that the proposed approach improves the accuracy and reduces the number of iterations required for the iterative solution of the resulting system of equations by 40% for moderately refined meshes and topological enrichment. Moreover, when a fixed enrichment volume is used, the number of iterations required grows at a rate which is reduced by a factor of 2 compared to standard XFEM, diminishing the number of iterations by almost one order of magnitude. [less ▲] Detailed reference viewed: 337 (11 UL) |
||