Stéphane Bordas

Emilio Martínez-Pañeda, Sundar Natarajan, Stéphane Bordas

Strain gradient plasticity theories are being widely used for fracture assessment, as they provide a richer description of crack tip fields by incorporating the influence of geometrically necessary dislocations. Characterizing the behavior at the small scales involved in crack tip deformation requires, however, the use of a very refined mesh within microns to the crack. In this work a novel and efficient gradient-enhanced numerical framework is developed by means of the extended finite element method (X-FEM). A mechanism-based gradient plasticity model is employed and the approximation of the displacement field is enriched with the stress singularity of the gradient-dominated solution. Results reveal that the proposed numerical methodology largely outperforms the standard finite element approach. The present work could have important implications on the use of microstructurally-motivated models in large scale applications. The non-linear X-FEM code developed in MATLAB can be downloaded from www.empaneda.com/codes.

Huu Phuoc Bui, Satyendra Tomar, Hadrien Courtecuisse et al.

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. 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.

Peng Yu, Cosmin Anitescu, Satyendra Tomar et al.

This paper presents a novel methodology of local adaptivity for the frequency-domain analysis of the vibrations of Reissner-Mindlin plates. The adaptive discretization is based on the recently developed Geometry Independent Field approximaTion (GIFT) framework, which may be seen as a generalisation of the Iso-Geometric Analysis (IGA). Within the GIFT framework, we describe the geometry of the structure exactly with NURBS (Non-Uniform Rational B-Splines), whilst independently employing Polynomial splines over Hierarchical T-meshes (PHT)-splines to represent the solution field. The proposed strategy of local adaptivity, wherein a posteriori error estimators are computed based on inexpensive hierarchical $h-$refinement, aims to control the discretisation error within a frequency band. The approach sweeps from lower to higher frequencies, refining the mesh appropriately so that each of the free vibration mode within the targeted frequency band is sufficiently resolved. Through several numerical examples, we show that the GIFT framework is a powerful and versatile tool to perform local adaptivity in structural dynamics. We also show that the proposed adaptive local $h-$refinement scheme allows us to achieve significantly faster convergence rates than when using a uniform $h-$refinement.