Implementation of meshless LBIE method to the 2D non‐linear SG problem

In this paper the meshless local boundary integral equation (LBIE) method for numerically solving the non‐linear two‐dimensional sine‐Gordon (SG) equation is developed. The method is based on the LBIE with moving least‐squares (MLS) approximation. For the MLS, nodal points spread over the analyzed domain are utilized to approximate the interior and boundary variables. The approximation functions are constructed entirely using a set of scattered nodes, and no element or connectivity of the nodes is needed for either the interpolation or the integration purposes. A time‐stepping method is employed to deal with the time derivative and a simple predictor–corrector scheme is performed to eliminate the non‐linearity. A brief discussion is outlined for numerical integrations in the proposed algorithm. Some examples involving line and ring solitons are demonstrated and the conservation of energy in undamped SG equation is investigated. The final numerical results confirm the ability of method to deal with the unsteady non‐linear problems in large domains. Copyright © 2009 John Wiley & Sons, Ltd.     

Limit analysis of plates using the EFG method and second‐order cone programming

The meshless element‐free Galerkin (EFG) method is extended to allow computation of the limit load of plates. A kinematic formulation that involves approximating the displacement field using the moving least‐squares technique is developed. Only one displacement variable is required for each EFG node, ensuring that the total number of variables in the resulting optimization problem is kept to a minimum, with far fewer variables being required compared with finite element formulations using compatible elements. A stabilized conforming nodal integration scheme is extended to plastic plate bending problems. The evaluation of integrals at nodal points using curvature smoothing stabilization both keeps the size of the optimization problem small and also results in stable and accurate solutions. Difficulties imposing essential boundary conditions are overcome by enforcing displacements at the nodes directly. The formulation can be expressed as the problem of minimizing a sum of Euclidean norms subject to a set of equality constraints. This non‐smooth minimization problem can be transformed into a form suitable for solution using second‐order cone programming. The procedure is applied to several benchmark beam and plate problems and is found in practice to generate good upper‐bound solutions for benchmark problems. Copyright © 2009 John Wiley & Sons, Ltd.     

Boundary element‐free method (BEFM) for two‐dimensional elastodynamic analysis using Laplace transform

In this paper, we present a direct meshless method of boundary integral equation (BIE), known as the boundary element‐free method (BEFM), for two‐dimensional (2D) elastodynamic problems that combines the BIE method for 2D elastodynamics in the Laplace‐transformed domain and the improved moving least‐squares (IMLS) approximation. The formulae for the BEFM for 2D elastodynamic problems are given, and the numerical procedures are also shown. The BEFM is a direct numerical method, in which the basic unknown quantities are the real solutions of the nodal variables, and the boundary conditions can be implemented directly and easily that leads to a greater computational precision. For the purpose of demonstration, some selected numerical examples are solved using the BEFM. Copyright © 2005 John Wiley & Sons, Ltd.     

A meshless method with conforming and nonconforming sub‐domains

An accurate and easy integration technique is desired for the meshless methods of weak form. As is well known, a sub‐domain method is often used in computational mechanics. The conforming sub‐domains, where the sub‐domains are not separated nor overlapped each other, are often used, while the nonconforming sub‐domains could be employed if needed. In the latter cases, the integrations of the sub‐domains may be performed easily by choosing a simple configuration. And then, the meshless method with nonconforming sub‐domains is considered one of the reasonable choices for computational mechanics without the troublesome integration. In this paper, we propose a new sub‐domain meshless method. It is noted that, because the method can employ both the conforming and the nonconforming sub‐domains, the integration for the weak form is necessarily accurate and easy by selecting the nonconforming sub‐domains with simple configuration. The boundary value problems including the Poisson's equation and the Helmholtz's equation are analyzed by using the proposed method. The numerical solutions are compared with the exact solutions and the solutions of the collocation method, showing that the relative errors by using the proposed method are smaller than those by using the collocation method and that the proposed method possesses a good convergence. Copyright © 2016 John Wiley & Sons, Ltd.     

Boundary element‐free method (BEFM) and its application to two‐dimensional elasticity problems

In this study, we first discuss the moving least‐square approximation (MLS) method. In some cases, the MLS may form an ill‐conditioned system of equations so that the solution cannot be correctly obtained. Hence, in this paper, we propose an improved moving least‐square approximation (IMLS) method. In the IMLS method, the orthogonal function system with a weight function is used as the basis function. The IMLS has higher computational efficiency and precision than the MLS, and will not lead to an ill‐conditioned system of equations. Combining the boundary integral equation (BIE) method and the IMLS approximation method, a direct meshless BIE method, the boundary element‐free method (BEFM), for two‐dimensional elasticity is presented. Compared to other meshless BIE methods, BEFM is a direct numerical method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied easily; hence, it has higher computational precision. For demonstration purpose, selected numerical examples are given. Copyright © 2005 John Wiley & Sons, Ltd.     

An enriched meshless method for non‐linear fracture mechanics

This paper presents an enriched meshless method for fracture analysis of cracks in homogeneous, isotropic, non‐linear‐elastic, two‐dimensional solids, subject to mode‐I loading conditions. The method involves an element‐free Galerkin formulation and two new enriched basis functions (Types I and II) to capture the Hutchinson–Rice–Rosengren singularity field in non‐linear fracture mechanics. The Type I enriched basis function can be viewed as a generalized enriched basis function, which degenerates to the linear‐elastic basis function when the material hardening exponent is unity. The Type II enriched basis function entails further improvements of the Type I basis function by adding trigonometric functions. Four numerical examples are presented to illustrate the proposed method. The boundary layer analysis indicates that the crack‐tip field predicted by using the proposed basis functions matches with the theoretical solution very well in the whole region considered, whether for the near‐tip asymptotic field or for the far‐tip elastic field. Numerical analyses of standard fracture specimens by the proposed meshless method also yield accurate estimates of the J‐integral for the applied load intensities and material properties considered. Also, the crack‐mouth opening displacement evaluated by the proposed meshless method is in good agreement with finite element results. Furthermore, the meshless results show excellent agreement with the experimental measurements, indicating that the new basis functions are also capable of capturing elastic–plastic deformations at a stress concentration effectively. Copyright © 2003 John Wiley & Sons, Ltd.     

Time‐dependent fractional advection–diffusion equations by an implicit MLS meshless method

Recently, many new applications in engineering and science are governed by a series of fractional partial differential equations (FPDEs). Unlike the normal partial differential equations (PDEs), the differential order in an FPDE is with a fractional order, which will lead to new challenges for numerical simulation, because most existing numerical simulation techniques are developed for the PDE with an integer differential order. The current dominant numerical method for FPDEs is finite difference method (FDM), which is usually difficult to handle a complex problem domain, and also difficult to use irregular nodal distribution. This paper aims to develop an implicit meshless approach based on the moving least squares (MLS) approximation for numerical simulation of fractional advection–diffusion equations (FADE), which is a typical FPDE The discrete system of equations is obtained by using the MLS meshless shape functions and the meshless strong‐forms. The stability and convergence related to the time discretization of this approach are then discussed and theoretically proven. Several numerical examples with different problem domains and different nodal distributions are used to validate and investigate the accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of the FADE. Copyright © 2011 John Wiley & Sons, Ltd.     

Topology optimization of structures using meshless density variable approximants

This paper proposes a new structural topology optimization method using a dual‐level point‐wise density approximant and the meshless Galerkin weak‐forms, totally based on a set of arbitrarily scattered field nodes to discretize the design domain. The moving least squares (MLS) method is used to construct shape functions with compactly supported weight functions, to achieve meshless approximations of system state equations. The MLS shape function with the zero‐order consistency will degenerate to the well‐known ‘Shepard function’, while the MLS shape function with the first‐order consistency refers to the widely studied ‘MLS shape function’. The Shepard function is then applied to construct a physically meaningful dual‐level density approximant, because of its non‐negative and range‐restricted properties. First, in terms of the original set of nodal density variables, this study develops a nonlocal nodal density approximant with enhanced smoothness by incorporating the Shepard function into the problem formulation. The density at any node can be evaluated according to the density variables located inside the influence domain of the current node. Second, in the numerical implementation, we present a point‐wise density interpolant via the Shepard function method. The density of any computational point is determined by the surrounding nodal densities within the influence domain of the concerned point. According to a set of generic design variables scattered at field nodes, an alternative solid isotropic material with penalization model is thus established through the proposed dual‐level density approximant. The Lagrangian multiplier method is included to enforce the essential boundary conditions because of the lack of the Kronecker delta function property of MLS meshless shape functions. Two benchmark numerical examples are employed to demonstrate the effectiveness of the proposed method, in particular its applicability in eliminating numerical instabilities. Copyright © 2012 John Wiley & Sons, Ltd.     

The natural neighbour Petrov–Galerkin method for elasto‐statics

In this paper, an efficient and accurate meshless natural neighbour Petrov–Galerkin method (NNPG) is proposed to solve elasto‐static problems in two‐dimensional space. This method is derived from the generalized meshless local Petrov–Galerkin method (MLPG) as a special case. In the NNPG, the local supported trial functions are constructed based on the non‐Sibsonian interpolation and test functions are taken as the three‐node triangular FEM shape functions. The local weak forms of the equilibrium equation and the boundary conditions are satisfied in local polygonal sub‐domains. These sub‐domains are constructed with Delaunay tessellations and domain integrals are evaluated over included Delaunay triangles by using Gaussian quadrature scheme. As this method combines the advantages of natural neighbour interpolation with Petrov–Galerkin method together, no stiffness matrix assembly is required and no special treatment is needed to impose the essential boundary conditions. Several numerical examples are presented and the results show the presented method is easy to implement and very accurate for these problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Least‐squares collocation meshless method

A finite point method, least‐squares collocation meshless method, is proposed. Except for the collocation points which are used to construct the trial functions, a number of auxiliary points are also adopted. Unlike the direct collocation method, the equilibrium conditions are satisfied not only at the collocation points but also at the auxiliary points in a least‐squares sense. The moving least‐squares interpolant is used to construct the trial functions. The computational effort required for the present method is in the same order as that required for the direct collocation, while the present method improves the accuracy of solution significantly. The proposed method does not require any mesh so that it is a truly meshless method. Three numerical examples are studied in detail, which show that the proposed method possesses high accuracy with low computational effort. Copyright © 2001 John Wiley & Sons, Ltd.     

Development of a 4‐node hybrid stress tetrahedral element using a node‐based smoothed finite element method

In this paper, a new 4‐node hybrid stress element is proposed using a node‐based smoothing technique of tetrahedral mesh. The conditions for hybrid stress field required are summarized, and the field should be continuous for better performance of a constant‐strain tetrahedral element. Nodal stress is approximated by the node‐based smoothing technique, and the stress field is interpolated with standard shape functions. This stress field is linear within each element and continuous across elements. The stress field is expressed by nodal displacements and no additional variables. The element stiffness matrix is calculated using the Hellinger‐Reissner functional, which guarantees the strain field from displacement field to be equal to that from the stress field in a weak sense. The performance of the proposed element is verified by through several numerical examples.     

Parallel computing of high‐speed compressible flows using a node‐based finite‐element method

An efficient parallel computing method for high‐speed compressible flows is presented. The numerical analysis of flows with shocks requires very fine computational grids and grid generation requires a great deal of time. In the proposed method, all computational procedures, from the mesh generation to the solution of a system of equations, can be performed seamlessly in parallel in terms of nodes. Local finite‐element mesh is generated robustly around each node, even for severe boundary shapes such as cracks. The algorithm and the data structure of finite‐element calculation are based on nodes, and parallel computing is realized by dividing a system of equations by the row of the global coefficient matrix. The inter‐processor communication is minimized by renumbering the nodal identification number using ParMETIS. The numerical scheme for high‐speed compressible flows is based on the two‐step Taylor–Galerkin method. The proposed method is implemented on distributed memory systems, such as an Alpha PC cluster, and a parallel supercomputer, Hitachi SR8000. The performance of the method is illustrated by the computation of supersonic flows over a forward facing step. The numerical examples show that crisp shocks are effectively computed on multiprocessors at high efficiency. Copyright © 2003 John Wiley & Sons, Ltd.     

The natural radial element method

In this work an innovative numerical approach is proposed, which combines the simplicity of low‐order finite elements connectivity with the geometric flexibility of meshless methods. The natural neighbour concept is applied to enforce the nodal connectivity. Resorting to the Delaunay triangulation a background integration mesh is constructed, completely dependent on the nodal mesh. The nodal connectivity is imposed through nodal sets with reduce size, reducing significantly the test function construction cost. The interpolations functions, constructed using Euclidian norms, are easily obtained. To prove the good behaviour of the proposed interpolation function several data‐fitting examples and first‐order partial differential equations are solved. The proposed numerical method is also extended to the elastostatic analysis, where classic solid mechanics benchmark examples are solved. Copyright © 2013 John Wiley & Sons, Ltd.     

Umbrella spherical integration: a stable meshless method for non‐linear solids

A stable meshless method for studying the finite deformation of non‐linear three‐dimensional (3D) solids is presented. The method is based on a variational framework with the necessary integrals evaluated through nodal integration. The method is truly meshless, requiring no 3D meshing or tessellation of any form. A local least‐squares approximation about each node is used to obtain necessary deformation gradients. The use of a local field approximation makes automatic grid refinement and the application of boundary conditions straightforward. Stabilization is achieved through the use of special ‘umbrella’ shape functions that have discontinuous derivatives at the nodes. Novel efficient algorithms for constructing the nodal stars and computing the nodal volumes are presented. The method is applied to four test problems: uniaxial tension, simple shear and bending of a bar, and cylindrical indentation. Convergence studies at infinitesimal strain show that the method is well‐behaved and converges with the number of nodes for both uniform and non‐uniform grids. Typical of meshless methods employing nodal integration, the total energy can be underestimated due to the approximate integration. At finite deformation the method reproduces known exact solutions. The bending example demonstrates an interesting example of torsional buckling resulting from the bending. Copyright © 2006 John Wiley & Sons, Ltd.     

Combined interior‐point method and semismooth Newton method for frictionless contact problems

In the present paper, a solution scheme is proposed for frictionless contact problems of linear elastic bodies, which are discretized using the finite element method with lower order elements. An approach combining the interior‐point method and the semismooth Newton method is proposed. In this method, an initial active set for the semismooth Newton method is obtained from the approximate optimal solution by the interior‐point method. The simplest node‐to‐node contact model is considered in the present paper, that is, pairs of matching nodes exist on the contact surfaces. However, the discussions can be easily extended to a node‐to‐segment or segment‐to‐segment contact model. In order to evaluate the proposed method, a number of illustrative examples of the frictionless contact problem are shown. The proposed combined method is compared with the interior‐point method and the semismooth Newton method. Two numerical examples that are difficult to solve using the semismooth Newton method are solved effectively using the proposed combined method. It is shown that the proposed method converges within far fewer iterations than the semismooth Newton methods or the interior‐point method. Copyright © 2009 John Wiley & Sons, Ltd.     

Non‐linear version of stabilized conforming nodal integration for Galerkin mesh‐free methods

A stabilized conforming (SC) nodal integration, which meets the integration constraint in the Galerkin mesh‐free approximation, is generalized for non‐linear problems. Using a Lagrangian discretization, the integration constraints for SC nodal integration are imposed in the undeformed configuration. This is accomplished by introducing a Lagrangian strain smoothing to the deformation gradient, and by performing a nodal integration in the undeformed configuration. The proposed method is independent to the path dependency of the materials. An assumed strain method is employed to formulate the discrete equilibrium equations, and the smoothed deformation gradient serves as the stabilization mechanism in the nodally integrated variational equation. Eigenvalue analysis demonstrated that the proposed strain smoothing provides a stabilization to the nodally integrated discrete equations. By employing Lagrangian shape functions, the computation of smoothed gradient matrix for deformation gradient is only necessary in the initial stage, and it can be stored and reused in the subsequent load steps. A significant gain in computational efficiency is achieved, as well as enhanced accuracy, in comparison with the mesh‐free solution using Gauss integration. The performance of the proposed method is shown to be quite robust in dealing with non‐uniform discretization. Copyright © 2002 John Wiley & Sons, Ltd.     

An element‐free method for in‐plane notch problems with anisotropic materials

This study developed an element‐free Galerkin method (EFGM) to simulate notched anisotropic plates containing stress singularities at the notch tip. Two‐dimensional theoretical complex displacement functions are first deduced into the moving least‐squares interpolation. The interpolation functions and their derivatives are then determined to calculate the nodal stiffness using the Galerkin method. In the numerical validation, an interface layer of the EFGM is used to combine the mesh between the traditional finite elements and the proposed singular notch EFGM. The H‐integral determined from finite element analyses with a very fine mesh is used to validate the numerical results of the proposed method. The comparisons indicate that the proposed method obtains more accurate results for the displacement, stress, and energy fields than those determined from the standard finite element method. Copyright © 2013 John Wiley & Sons, Ltd.     

Application of piece‐wise linear weight functions for 2D 8‐node quadrilateral element in contact problems

The present study is a continuation of our previous work with the aim to reduce problems caused by standard higher order elements in contact problems. The difficulties can be attributed to the inherent property of the Galerkin method which gives uneven distributions of nodal forces resulting in oscillating contact pressures. The proposed remedy is use of piece‐wise linear weight functions. The methods to establish stiffness and/or mass matrix for 8‐node quadrilateral element in 2D are presented, i.e. the condensing and direct procedures. The energy and nodal displacement error norms are also checked to establish the convergence ratio. Interpretation of calculated contact pressures is discussed. Two new 2D 8‐node quadrilateral elements, QUAD8C and QUAD8D, are derived and tested in many examples, which show their good performance in contact problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Complex variable moving least‐squares method: a meshless approximation technique

Based on the moving least‐squares (MLS) approximation, we propose a new approximation method—the complex variable moving least‐squares (CVMLS) approximation. With the CVMLS approximation, the trial function of a two‐dimensional problem is formed with a one‐dimensional basis function. The number of unknown coefficients in the trial function of the CVMLS approximation is less than in the trial function of the MLS approximation, and we can thus select fewer nodes in the meshless method that is formed from the CVMLS approximation than are required in the meshless method of the MLS approximation with no loss of precision. The meshless method that is derived from the CVMLS approximation also has a greater computational efficiency. From the CVMLS approximation, we propose a new meshless method for two‐dimensional elasticity problems—the complex variable meshless method (CVMM)—and the formulae of the CVMM for two‐dimensional elasticity problems are obtained. Compared with the conventional meshless method, the CVMM has a greater precision and computational efficiency. For the purposes of demonstration, some selected numerical examples are solved using the CVMM. Copyright © 2006 John Wiley & Sons, Ltd.     

A perturbation method for stochastic meshless analysis in elastostatics

A stochastic meshless method is presented for solving boundary‐value problems in linear elasticity that involves random material properties. The material property was modelled as a homogeneous random field. A meshless formulation was developed to predict stochastic structural response. Unlike the finite element method, the meshless method requires no structured mesh, since only a scattered set of nodal points is required in the domain of interest. There is no need for fixed connectivities between nodes. In conjunction with the meshless equations, classical perturbation expansions were derived to predict second‐moment characteristics of response. Numerical examples based on one‐ and two‐dimensional problems are presented to examine the accuracy and convergence of the stochastic meshless method. A good agreement is obtained between the results of the proposed method and Monte Carlo simulation. Since mesh generation of complex structures can be a far more time‐consuming and costly effort than the solution of a discrete set of equations, the meshless method provides an attractive alternative to finite element method for solving stochastic mechanics problems. Copyright © 2001 John Wiley & Sons, Ltd.