TABLE OF CONTENTS Two main purposes of developing meshfree methods are
to relieve the users or the analysts from the trivial and time-consuming task of meshing, and
to solve problems with extremely large deformation and even breakage.
The computer implementations of meshfree methods are, in general, more difficult than that in a grid/mesh-based method, simply because there is no predefined grid/mesh to use in establishing the discrete system equations. Among the meshfree methods, SPH method is relatively simple to code. A detailed discussion on issues and techniques related to computer implementation of a meshfree method has been given by Liu (2002) in his recent monograph.
In this chapter, issues related to the computer implementation of meshfree particle methods are discussed. Computer implementations of meshfree particle method over serial and parallel computers are also briefly addressed. The source code of a standard serial 3D SPH code is provided. The main features of the SPH code, detailed descriptions and the source code of the related subroutines are also provided.
The programs demonstrate most of the concepts and techniques related to the SPH method. The code can be easily extended to other corrective or modified versions of SPH with proper treatment either on the kernel approximation or the particle approximation. The 3D SPH code can solve compressible flows with real viscosity, and can be readily modified for applications to hydrodynamics with material strength with a proper constitutive model and an equation of state. As a structured code, it can be readily modified or imported to...
