The interpolating element-free Galerkin method for three-dimensional elastoplasticity problems

Analysis of two methods based on Galerkin weak form for fractional diffusion-wave: Meshless interpolating element free Galerkin (IEFG) and finite element methods

Shakedown theorems are exact theories of classical plasticity for the direct computation of the load-carrying capacity under varying loads. Based on the classical Melan’s theorem, a numerical solution procedure for determining the shakedown load of elasto-perfectly plastic structure is presented firstly making use of the element free Galerkin (EFG) method. The numerical implementation is very simple and convenient because it is only necessary to construct an array of nodes in the domain under consideration. The reduced-basis technique is adopted here to solve the mathematical programming iteratively in a sequence of reduced self-equilibrium stress subspaces with very low dimensions. The self-equilibrium stress field is expressed by linear combination of several self-equilibrium stress basis vectors with parameters to be determined. These self-equilibrium stress basis vectors can be generated by performing equilibrium iteration procedure during elasto-plastic incremental analysis. The Complex method is used to solve nonlinear programming and determine the maximal load amplifier. The numerical results show that it is efficient and accurate to solve shakedown analysis problems by using the EFG method and nonlinear programming.

The element-free Galerkin (EFG) method is one of the widely used meshfree methods for solving partial differential equations. In the EFG method, shape functions are derived from a moving least-squares (MLS) approximation, which involves the inversion of a small matrix for every point of interest. To avoid the calculation of matrix inversion in the formulation of the shape functions, an improved MLS approximation is presented, where an orthogonal function system with a weight function is used. However, it can also lead to ill-conditioned or even singular system of equations. In this paper, aspects of the IMLS approximation are analyzed in detail. The reason why singularity problem occurs is studied. A novel technique based on matrix triangular process is proposed to solve this problem. It is shown that the EFG method with present technique is very effective in constructing shape functions. Numerical examples are illustrated to show the efficiency and accuracy of the proposed method. Although our study relies on monomial basis functions, it is more general than existing methods and can be extended to any basis functions.

A new formulation of the Element Free Galerkin (EFG) method is developed for the modelling of incompressible materials. Beginning with a mixed variational principle, a selective reduced integration procedure is developed by implementing nodal quadrature. Numerical examples are provided which compare the performance of the proposed technique to the standard EFG approximation. These studies illustrate the capability of the new formulation to eliminate volumetric locking. For the standard method, however, the degree of volumetric locking is shown to be a function of the local support sizes (domains of influence) of the EFG shape functions. Copyright © 1999 John Wiley & Sons, Ltd.

In this paper, a new method for deriving the moving least-squares (MLS) approximation is presented, and the interpolating moving least-squares (IMLS) method proposed by Lancaster is improved. Compared with the IMLS method proposed by Lancaster, a simpler formula of the shape function is given in the improved IMLS method in this paper so that the new method has higher computing efficiency. Combining the shape function constructed by the improved IMLS method with Galerkin weak form of the elasticity problems, the interpolating element-free Galerkin (IEFG) method for the two-dimensional elasticity problems is presented, and the corresponding formulae are obtained. In the IEFG method, the boundary conditions can be applied directly which makes the computing efficiency higher than the conventional EFG method. Some numerical examples are presented to demonstrate the validity of the method.

This study aims to interpret the non-linear steady-state heat conduction for temperature-dependent thermal conductivity ($k$ ($T$)) using Element-Free Galerkin (EFG) method. In this present study, a one-dimensional heat conduction problem with uniform heat generation was explicated. Moving Least Squares (MLS) approximants were applied to estimate the unknown function of temperature $T$ ($x$) with $T^h$ ($x$) using linear basis and weight functions. The variational method has been used to develop discrete equations. Essential boundary conditions are enforced by using the Penalty method. The results have been obtained for the one-dimensional model using essential MATLAB codes. The results obtained by the EFG method are compared with the analytical and finite-element method results. The results are also studied by increasing the number of nodes to study the convergence which indicated that EFG has good convergence behavior. The results have also been obtained for different values of the scaling parameter ($α_s$) and any values of αs between 1.8 and 2.0 were found suitable for providing better results in the EFG method.

Lower bound shakedown analysis by using the element free Galerkin method and non-linear programming

Transient heat transfer analysis of orthotropic materials considering phase change process based on element-free Galerkin method

In this paper, the adaptive nodal generation procedure based on the estimated local and global error in the element-free Galerkin (EFG) method is proposed. To investigate the possibility of h-type adaptivity of EFG method, a simple nodal refinement scheme is used. By adding new node along the background cell that is used in numerical integration, both of the local and global errors can be controlled adaptively. These errors are estimated by calculating the difference between the values of the projected stresses and original EFG stresses. The ultimate goal of this study is to develop the reliable nodal generator based on the local and global errors that is estimated posteriori. To evaluate the performance of proposed adaptive procedure, the convergence behavior is investigated for several examples.

Analysis and application of the interpolating element free Galerkin (IEFG) method to simulate the prevention of groundwater contamination with application in fluid flow

Meshfree analyses of cable-reinforced membrane structures by ALE–EFG method

In this paper, we present a procedure to estimate the error in elliptic equations using the element‐free Galerkin (EFG) method, whose evaluation is computationally simple and can be readily implemented in existing EFG codes. The estimation of the error works very well in all numerical examples for 2‐D potential problems that are presented here, for regular and irregular clouds of points. Moreover, it was demonstrated that this method is very simple in terms of economy and gives a good performance.The results show that the error in EFG approximation may be estimated via the error estimator described in this paper. The quality of the estimation of the error is demonstrated by numerical examples. The implemented procedure of error approximation allows the global energy norm error to be estimated and also gives a good evaluation of local errors. It can, thus, be combined with a full adaptive process of refinement or, more simply, provide guidance for redesign of cloud of points. Copyright © 2001 John Wiley & Sons, Ltd.

Stability of piezoelectric FGM rectangular plates subjected to non-uniformly distributed load, heat and voltage

Based on the improved complex variable moving least-squares (ICVMLS) approximation, the improved complex variable element-free Galerkin (ICVEFG) method for the bending problem of Kirchhoff plate is presented. Compared with the moving least-squares (MLS) approximation, in the ICVMLS approximation, the approximation function of two-dimensional problems can be obtained with one-dimensional basis function, then the computational efficiency of the shape functions is higher. Compared with the meshless methods based on the MLS approximation, under the same node distributions, the ones using the ICVMLS approximation can obtain the solutions with higher computational accuracy; and under the similar computational accuracy, the ones using the ICVMLS approximation have higher computational efficiency. The ICVMLS approximation is used to form the approximation function of the deflection of a Kirchhoff plate, the Galerkin weak form of the bending problem of Kirchhoff plates is adopted to obtain the discretized system equations, and the penalty method is employed to enforce the essential boundary conditions, then the corresponding formulae of the ICVEFG method for the bending problem of Kirchhoff plates are presented. By computing and analyzing four typical examples, it is shown that the ICVEFG method of Kirchhoff plates in this paper is efficient. And the computational precision of the numerical solutions is analyzed to select the basis function, weight function, scaling factor, node distribution and penalty factor in the ICVEFG method. Numerical examples show that the method in this paper has better convergence and higher accuracy. When the quadratic polynomial basis function and the cubic or quartic spline weight function are used, and d m a x = 4.2−4.4 , the ICVEFG method of Kirchhoff plates in this paper can obtain the solutions with high computational accuracy. And more nodes are distributed in the problem domain, higher computational accuracy of the solution will obtained, which shows that the method in this paper has better convergence.

The mathematical model of thermomechanical coupling multi-objective topology optimization for anisotropic structures was established using the element-free Galerkin (EFG) method, and the weighted function of compliance and heat dissipation was defined as the objective function. The proposed model and procedure were verified by the topological results based on the finite element method. The multi-objective EFG optimal topological structures have clearer boundary profiles even without using the sensitivity filtering technique. The effects of the weight coefficient, thermal conductivity factor, Poisson's ratio factor, off-angle and volume fraction on the multi-objective EFG optimal topological structure and multi-objective function were evaluated in detail, and reasonable ranges of the above parameters were recommended to improve the heat dissipation and mechanical performance. The multi-objective optimal results of the anisotropic structure were 3D printed and compared with the isotropic material, and their temperature, displacement and stress were improved, which reflects the advantages of orthotropic structures.