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

Real-time error control for surgical simulation Bui, Huu Phuoc ; Tomar, Satyendra ; et al Poster (2016, December 12) Objective: To present the first real-time a posteriori error-driven adaptive finite element approach for real-time simulation and to demonstrate the method on a needle insertion problem. Methods: We use ... [more ▼] Objective: To present the first real-time a posteriori error-driven adaptive finite element approach for real-time simulation and to demonstrate the method on a needle insertion problem. Methods: We use corotational elasticity and a frictional needle/tissue interaction model based on friction. The problem is solved using finite elements within SOFA. The refinement strategy relies upon a hexahedron-based finite element method, combined with a posteriori error estimation driven local $h$-refinement, for simulating soft tissue deformation. Results: We control the local and global error level in the mechanical fields (e.g. displacement or stresses) during the simulation. We show 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. Conclusions: Error control guarantees that a tolerable error level is not exceeded during the simulations. Local mesh refinement accelerates simulations. Significance: Our work provides a first step to discriminate between discretization error and modeling error by providing a robust quantification of discretization error during simulations. [less ▲] Detailed reference viewed: 246 (20 UL)Bayesian inference for parameter identification in computational mechanics Rappel, Hussein ; Beex, Lars ; Hale, Jack et al Poster (2016, December 12) Detailed reference viewed: 181 (10 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: 226 (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: 202 (15 UL)Uncertainty quantification for soft tissue biomechanics Hauseux, Paul ; Hale, Jack ; Bordas, Stéphane Poster (2016, December) Detailed reference viewed: 232 (20 UL)Shape Optimization Directly from CAD: an Isogeometric Boundary Element Approach Using T-splines ; ; Bordas, Stéphane Report (2016) Detailed reference viewed: 323 (5 UL)Isogeometric finite element analysis of time-harmonic exterior acoustic scattering problems ; ; Bordas, Stéphane E-print/Working paper (2016) We present an isogeometric analysis of time-harmonic exterior acoustic problems. The infinite space is truncated by a fictitious boundary and (simple) absorbing boundary conditions are applied. The ... [more ▼] We present an isogeometric analysis of time-harmonic exterior acoustic problems. The infinite space is truncated by a fictitious boundary and (simple) absorbing boundary conditions are applied. The truncation error is included in the exact solution so that the reported error is an indicator of the performance of the isogeometric analysis, in particular of the related pollution error. Numerical results performed with high-order basis functions (third or fourth orders) showed no visible pollution error even for very high frequencies. This property combined with exact geometrical representation makes isogeometric analysis a very promising platform to solve high-frequency acoustic problems. [less ▲] Detailed reference viewed: 195 (16 UL)Bayesian inference for material parameter identification in elastoplasticity Rappel, Hussein ; Beex, Lars ; Hale, Jack et al Scientific Conference (2016, September 07) Detailed reference viewed: 247 (29 UL)Stochastic FE analysis of brain deformation with different hyper-elastic models Hauseux, Paul ; Hale, Jack ; Bordas, Stéphane Scientific Conference (2016, September) Detailed reference viewed: 266 (29 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: 191 (7 UL)Numerical studies of magnetic particles concentration in biofluid (blood) under the influence of high gradient magnetic field in microchannel ; Bourantas, Georgios ; et al Scientific Conference (2016, July 15) Detailed reference viewed: 161 (9 UL)Linear elastic fracture simulation directly from CAD: 2D NURBS-based implementation and role of tip enrichment ; ; et al in International Journal of Fracture (2016) Detailed reference viewed: 286 (21 UL)Simulating topological changes in real time for surgical assistance Bordas, Stéphane ; ; et al Speeches/Talks (2016) Detailed reference viewed: 445 (38 UL)A Bayesian approach for parameter identification in elastoplasticity Rappel, Hussein ; Beex, Lars ; Hale, Jack et al Scientific Conference (2016, June 09) Detailed reference viewed: 149 (15 UL)Phase field approach to fracture: Towards the simulation of cutting soft tissues Ziaei Rad, Vahid ; Hale, Jack ; et al Scientific Conference (2016, June 08) Detailed reference viewed: 333 (11 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: 176 (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: 151 (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)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: 179 (15 UL) |
||