Stéphane Bordas

Review: Meshless methods: A review and computer implementation aspects

The aim of this manuscript is to give a practical overview of meshless methods (for solid mechanics) based on global weak forms through a simple and well-structured MATLAB code, to illustrate our discourse. The source code is available for download on our website and should help students and researchers get started with some of the basic meshless methods; it includes intrinsic and extrinsic enrichment, point collocation methods, several boundary condition enforcement schemes and corresponding test cases. Several one and two-dimensional examples in elastostatics are given including weak and strong discontinuities and testing different ways of enforcing essential boundary conditions.

Meshless methods: review and key computer implementation aspects

The aim of this manuscript is to give a practical overview of meshless methods (for solid mechanics) based on global weak forms through a simple and well-structured MATLAB code, to illustrate our discourse. The source code is available for download on our website and should help students and researchers get started with some of the basic meshless methods; it includes intrinsic and extrinsic enrichment, point collocation methods, several boundary condition enforcement schemes and corresponding test cases. Several one and two-dimensional examples in elastostatics are given including weak and strong discontinuities and testing different ways of enforcing essential boundary conditions.

NURBS-based finite element analysis of functionally graded plates: static bending, vibration, buckling and flutter

In this paper, a non-uniform rational B-spline based iso-geometric finite element method is used to study the static and dynamic characteristics of functionally graded material (FGM) plates. The material properties are assumed to be graded only in the thickness direction and the effective properties are computed either using the rule of mixtures or by Mori–Tanaka homogenization scheme. The plate kinematics is based on the first order shear deformation plate theory (FSDT). The shear correction factors are evaluated employing the energy equivalence principle and a simple modification to the shear correction factor is presented to alleviate shear locking. Static bending, mechanical and thermal buckling, linear free flexural vibration and supersonic flutter analysis of FGM plates are numerically studied. The accuracy of the present formulation is validated against available three-dimensional solutions. A detailed numerical study is carried out to examine the influence of the gradient index, the plate aspect ratio and the plate thickness on the global response of functionally graded material plates.

Isogeometric analysis: An overview and computer implementation aspects

Isogeometric analysis (IGA) represents a recently developed technology in computational mechanics that offers the possibility of integrating methods for analysis and Computer Aided Design (CAD) into a single, unified process. The implications to practical engineering design scenarios are profound, since the time taken from design to analysis is greatly reduced, leading to dramatic gains in efficiency. In this manuscript, through a self-contained Matlab® implementation, we present an introduction to IGA applied to simple analysis problems and the related computer implementation aspects. Furthermore, implementation of the extended IGA which incorporates enrichment functions through the partition of unity method (PUM) is also presented, where several examples for both two-dimensional and three-dimensional fracture are illustrated. We also describe the use of IGA in the context of strong-form (collocation) formulations, which has been an area of research interest due to the potential for significant efficiency gains offered by these methods. The code which accompanies the present paper can be applied to one, two and three-dimensional problems for linear elasticity, linear elastic fracture mechanics, structural mechanics (beams/plates/shells including large displacements and rotations) and Poisson problems with or without enrichment. The Bezier extraction concept that allows the FE analysis to be performed efficiently on T-spline geometries is also incorporated. The article includes a summary of recent trends and developments within the field of IGA.

Nitsche's method for two and three dimensional NURBS patch coupling

We present a Nitche’s method to couple non-conforming two and three-dimensional non uniform rational b-splines (NURBS) patches in the context of isogeometric analysis. We present results for linear elastostatics in two and and three-dimensions. The method can deal with surface-surface or volume-volume coupling, and we show how it can be used to handle heterogeneities such as inclusions. We also present preliminary results on modal analysis. This simple coupling method has the potential to increase the applicability of NURBS-based isogeometric analysis for practical applications.

Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework†

Partition of unity methods, such as the extended finite element method, allows discontinuities to be simulated independently of the mesh (Int. J. Numer. Meth. Engng. 1999; 45:601-620). This eliminates the need for the mesh to be aligned with the discontinuity or cumbersome re-meshing, as the discontinuity evolves. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity is commonly adopted. In this paper, we use a simple integration technique, proposed for polygonal domains (Int. J. Nuttier Meth. Engng 2009; 80(1):103-134. DOI: 10.1002/nme.2589) to suppress the need for element subdivision. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics and a multi-material problem show that the proposed method yields accurate results. Owing to its simplicity, the proposed integration technique can be easily integrated in any existing code. Copyright (C) 2010 John Wiley & Sons, Ltd.

Natural frequencies of cracked functionally graded material plates by the extended finite element method

In this paper, the linear free flexural vibration of cracked functionally graded material plates is studied using the extended finite element method. A 4-noded quadrilateral plate bending element based on field and edge consistency requirement with 20 degrees of freedom per element is used for this study. The natural frequencies and mode shapes of simply supported and clamped square and rectangular plates are computed as a function of gradient index, crack length, crack orientation and crack location. The effect of thickness and influence of multiple cracks is also studied.

Real-time simulation of contact and cutting of heterogeneous soft-tissues

This paper presents a numerical method for interactive (real-time) simulations, which considerably improves the accuracy of the response of heterogeneous soft-tissue models undergoing contact, cutting and other topological changes. We provide an integrated methodology able to deal both with the ill-conditioning issues associated with material heterogeneities, contact boundary conditions which are one of the main sources of inaccuracies, and cutting which is one of the most challenging issues in interactive simulations. Our approach is based on an implicit time integration of a non-linear finite element model. To enable real-time computations, we propose a new preconditioning technique, based on an asynchronous update at low frequency. The preconditioner is not only used to improve the computation of the deformation of the tissues, but also to simulate the contact response of homogeneous and heterogeneous bodies with the same accuracy. We also address the problem of cutting the heterogeneous structures and propose a method to update the preconditioner according to the topological modifications. Finally, we apply our approach to three challenging demonstrators: (i) a simulation of cataract surgery (ii) a simulation of laparoscopic hepatectomy (iii) a brain tumor surgery.

Isogeometric Analysis of Laminated Composite Plates Using the Higher-Order Shear Deformation Theory

Isogeometric analysis (IGA) aims at simplifying the computer aided design (CAD) and computer aided engineering (CAE) pipeline by using the same functions to describe the geometry (CAD) and the unknown fields (Analysis). IGA can be based on a variety of CAD descriptions, the most widely used today being non-uniform rational B-splines (NURBS). In this article, the suitability of NURBS-based isogeometric analysis within a third-order shear deformation theory for the simulation of the static, dynamic, and buckling response of laminated composite plates is investigated. The method employs NURBS basis functions to both represent the geometry (exactly) and the unknown field variables. One of the main advantages of the present method is directly inherited from IGA, that is to easily increase the approximation order. To avoid using a shear correction factor, a third-order shear deformation theory (TSDT) is introduced. It requires C1-continuity of generalized displacements and the NURBS basis functions are well suited for this requirement. Several numerical examples are used to demonstrate the performance of the present method compared with other published ones.

On the approximation in the smoothed finite element method (SFEM)

This letter aims at resolving the issues raised in the recent short communication [1] and answered by [2] by proposing a systematic approximation scheme based on non-mapped shape functions, which both allows to fully exploit the unique advantages of the smoothed finite element method (SFEM) [3, 4, 5, 6, 7, 8, 9] and resolve the existence, linearity and positivity deficiencies pointed out in [1]. We show that Wachspress interpolants [10] computed in the physical coordinate system are very well suited to the SFEM, especially when elements are heavily distorted (obtuse interior angles). The proposed approximation leads to results which are almost identical to those of the SFEM initially proposed in [3]. These results that the proposed approximation scheme forms a strong and rigorous basis for construction of smoothed finite element methods.