Chapter 1
Introduction
âMeshlessâ methods are alternative techniques to the finite element method in solving partial differential equations. While the finite element method derives an approximation based on the elements, using shape functions, the meshless methods allow us to derive an approximation at any point; thanks to the information provided by the surrounding nodes. In these approaches the concept of element is thus not used any more. Connectivity between the nodes is not defined any more by the mesh but only by the concepts of âvicinityâ or âfield of influence.â These methods were developed with the aim of avoiding the numerical problems involved in mesh construction. These problems have been discussed in many studies; it is, forexample, a question of simulation of manufacturing processes such as extrusion, injection, or setting forms by removal of matter where it is necessary to face extremely large distortions of the mesh. In other processes such as foundry, drilling, or laser welding, precisely knowing the position of the interface between the solid phase and the liquid phase is essential. In the simulation of processes such as cutting by adiabatic shearing which involves a localized deformation, possibly accompanied by the propagation of a fissure, it is necessary to carry out the simulation without the mesh being conceived influencing the direction of propagation of the shear band or the fissure. The appearance of a localized deformation requires a finer representation of the solution in certain areas of the domains, and it is thus necessary to be able to refine the mesh easily without the geometrical constraints known within the framework of finite elements (mainly in 3D) and the problems related to precise projection of the fields between the two meshes. The objective of the meshless methods is to eliminate the structure of the mesh and to build the approximation starting only from the nodes. Although structures with a geometrical character are necessary (to build node connectivity for the integration of the weak form associated with the equation to be solved and so on), these do not interfere, in general, with the quality of the solution and thus can be built independently. Even after being proposed at the end of 1970s, the âmeshlessâ methods had to wait approximately 15 years before having a real development and an interest within the scientific community.
In the interval, little passion had been shown for them because of the numerous difficulties presented by the first techniques. The first âmeshlessâ method seems to be the so-called smooth particle hydrodynamics (SPH) method (Lucy 1977), which was initially used to model astronomical phenomena in unbounded domains. This method, based on an approximation using the properties of the convolution product, has two disadvantages: low consistency and difficulty associated with the imposition of boundary conditions. In 1992, Nayrolles, Touzot and Villon proposed using a local approximation of least squares in a new method called the âdiffuse elements methodâ (DEM). In 1994, Belytschko et al. proposed the âelement-free Galerkinâ (EFG) method based on the same principles as the preceding one but using âexactâ derivatives of the shape functions. The method known as the âreproducing kernel particle methodâ (RKPM) introduced by Liu et al. in 1995 is an extension of the SPH method but with the reproduction of linear fields or polynomials of higher order being introduced, thanks to the correction function affecting the kernel function used in SPH method. Finally, the so-called partition of unity method introduced by Babuska in about 1996 is a general principle allowing us to enrich any function associated with a problem involving known physics, within the framework of finite elements and of meshless methods, by adding additional unknowns in the global system of equations. Thus, particular functions such as discontinuous functions and singular functions can be reproduced.
Lastly, more recently, the natural element method (NEM) rests on principles completely different from the previous ones. This method is halfway between meshless methods and the finite element method. The NEM proposes an interpolation based on the concepts of the Voronoi diagram and its natural neighbors. The Voronoi diagram associated with a cloud of nodes distributed over the domain to be studied is the Delaunay dual mesh. Thus, a mesh is being used for the construction of the interpolation. However, as the examples presented in this chapter show, the quality of interpolation produced does not depend on the form of the triangles (2D problems) or tetrahedrons (3D problems) present in the Delaunay mesh. The latter is built in a systematic way without requiring repositioning of nodes. With NEM the choice of support of shape functions is automatic and optimal in the sense that node vicinity is taken into account as much as possible to define the interpolation. Wit...