Radial Point Interpolation Technique Research Articles

3D meshless modeling of piezoelectric structure based on the radial point interpolation method

Piezoelectric materials find widespread applications in high-precision actuators and sensors. However, the traditional finite element method falls short in meeting the simulation needs of piezoelectric structures due to complexities in mesh generation and precision requirements for accurate simulations. This paper focuses on adapting and generalizing the meshless modeling technique based on the radial point interpolation method to address three-dimensional simulations of piezoelectric structures and accurately model positive and inverse piezoelectric behaviors. The piezoelectric body is discretized into point cloud data, and displacement and potential are approximated by interpolating values through the radial point interpolation method. Subsequently, a 3D meshless model for general piezoelectric structures is formulated based on the principle of minimum potential energy. The influences of key factors on interpolation accuracy and computational cost were discussed. Numerical examples, featuring a 3D piezoelectric bimorph beam as actuators and sensors, validate the efficacy of this extended approach. Results affirm that the three-dimensional meshless method, rooted in the radial point interpolation technique, achieves superior accuracy in predicting piezoelectric behavior compared to the traditional finite element method. This adaptation not only enhances the applicability of meshless modeling but also represents a significant advancement in accurately simulating complex three-dimensional piezoelectric structures.

A local meshless method for solving multi-dimensional Galilei invariant fractional advection–diffusion equation

In this article, a local meshless technique is applied for numerical simulation of multi-dimensional Galilei invariant fractional advection–diffusion model on regular and irregular computational domains. In the suggested method, a second-order Crank–Nicolson scheme along with the second-order weighted and shifted Grünwald difference (WSGD) formula, is used to discretize the time derivatives of the model. This time-discretization scheme is unconditionally stable and convergent with order O(τ2). To approximate the spatial derivatives of this model, a local radial point interpolation technique is employed. Finally, to prove and demonstrate the validity of the proposed algorithm, various one, two and three-dimensional problems are investigated on regular and irregular computational domains.

Mesh Free Radial Point Interpolation Based Displacement Recovery Techniques for Elastic Finite Element Analysis

The study develops the displacement error recovery method in a mesh free environment for the finite element solution employing the radial point interpolation (RPI) technique. The RPI technique uses the radial basis functions (RBF), along with polynomials basis functions to interpolate the displacement fields in a node patch and recovers the error in displacement field. The global and local errors are quantified in both energy and L2 norms from the post-processed displacement field. The RPI technique considers multi-quadrics/gaussian/thin plate splint RBF in combination with linear basis function for displacement error recovery analysis. The elastic plate examples are analyzed to demonstrate the error convergence and effectivity of the RPI displacement recovery procedures employing mesh free and mesh dependent patches. The performance of a RPI-based error estimators is also compared with the mesh dependent least square based error estimator. The triangular and quadrilateral elements are used for the discretization of plates domains. It is verified that RBF with their shape parameters, choice of elements, and errors norms influence considerably on the RPI-based displacement error recovery of finite element solution. The numerical results show that the mesh free RPI-based displacement recovery technique is more effective and achieve target accuracy in adaptive analysis with the smaller number of elements as compared to mesh dependent RPI and mesh dependent least square. It is also concluded that proposed mesh free recovery technique may prove to be most suitable for error recovery and adaptive analysis of problems dealing with large domain changes and domain discontinuities.

Analysis of antisymmetric cross-ply laminates using high-order shear deformation theories: a meshless approach

For many years finite element method (FEM) was the chosen numerical method for the analysis of composite structures. However, in the last 20 years, the scientific community has witnessed the birth and development of several meshless methods, which are more flexible and equally accurate numerical methods. The meshless method used in this work is the natural neighbour radial point interpolation method (NNRPIM). In order to discretize the problem domain, the NNRPIM only requires an unstructured nodal distribution. Then, using the Voronoi mathematical concept, it enforces the nodal connectivity and constructs the background integration mesh. The NNRPIM shape functions are constructed using the radial point interpolation technique. In this work, the displacement field of composite laminated plates is defined by high-order shear deformation theories. In the end, several antisymmetric cross-ply laminates were analysed and the NNRPIM solutions were compared with the literature. The obtained results show the efficiency and accuracy of the NNRPIM formulation.

Geometrically nonlinear dynamic analysis of FG graphene platelets-reinforced nanocomposite cylinder: MLPG method based on a modified nonlinear micromechanical model

The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.

Numerical simulation of electromagnetic wave scattering from perfectly conducting cylinders using the local radial point interpolation technique

ABSTRACTThis paper is devoted to the investigation of the electromagnetic scattering problems from infinite perfectly conducting cylinders with arbitrary cross-sections. The problems can be mathematically modeled as surface integral equations. Both the electric field integral equation (EFIE), magnetic field integral equation (MFIE) and their combined form are studied. An efficient computational technique based on the local radial point interpolation method is performed to investigate the models. Some test problems with various cross-sections have been analyzed to validate the capability of the proposed technique in calculating the current distribution induced on the scatterer and measuring the Radar Cross Section.

Variational multiscale element-free Galerkin method combined with the moving Kriging interpolation for solving some partial differential equations with discontinuous solutions

Solving partial differential equations with discontinuous solutions is an important challenging problem in numerical analysis. To this end, there are some methods such as finite volume method, discontinuous Galerkin approach and particle technique that are able to solve these problems. In the current paper, the moving Kriging element-free Galerkin method has been combined with the variational multiscale algorithm to obtain acceptable and high-resolution solutions. For testing this technique, we select some PDEs with discontinuous solution such as Burgers’, Sod’s shock tube, advection–reaction–diffusion, Kuramoto–Sivashinsky, Boussinesq and shallow water equations. First, we obtain a time-discrete scheme by approximating time derivative via finite difference technique. Then we introduce the moving Kriging interpolation and also obtain their shape functions. We use the element-free Galerkin method for approximating the spatial derivatives. This method uses a weak form of the considered equation that is similar to the finite element method with the difference that in the classical element-free Galerkin method test and trial functions are moving least squares approximation (MLS) shape functions. Since the shape functions of moving least squares (MLS) approximation do not have Kronecker delta property, we cannot implement the essential boundary condition, directly. Thus, we employ the shape functions of moving Kriging interpolation and radial point interpolation technique which have the mentioned property. Also, in the element-free Galekin method, we do not use any triangular, quadrangular or other type of meshes. The element-free Galerkin method is a global method while finite element method is a local one. This technique employs a background mesh for integration which makes it different from the truly mesh procedures. Several test problems are solved and numerical simulations are reported which confirm the efficiency of the proposed schemes.

Error analysis of a meshless weak form method based on radial point interpolation technique for Sivashinsky equation arising in the alloy solidification problem

In this paper, meshless weak form techniques are applied to find the numerical solution of nonlinear biharmonic Sivashinsky equation arising in the alloy solidification problem. Stability and convergence analysis of time-discrete scheme are proved. An error analysis of meshless global weak form method based on radial point interpolation technique is proposed for this nonlinear biharmonic equation. In addition, a comparison between meshless global and local weak form methods is done from the perspective of accuracy and efficiency. The main purpose of this paper is to show that the meshless weak form techniques can be used for solving the nonlinear biharmonic partial differential equations especially Sivashinsky equation. The numerical results confirm the good efficiency of the proposed methods for solving this nonlinear biharmonic model.

The use of element free Galerkin method based on moving Kriging and radial point interpolation techniques for solving some types of Turing models

In this paper two numerical procedures are presented for solving a class of Turing system. Firstly, we obtain a time discrete scheme by approximating time derivative via finite difference technique. Then we introduce the moving Kriging interpolation and radial point interpolation and also obtain their shape functions. We use the element free Galerkin method for approximating the spatial derivatives. This method uses a weak form of the considered equation that is similar to the finite element method with the difference that in the classical element free Galerkin method test and trial functions are moving least squares (MLS) approximation shape functions. Since the shape functions of moving least squares (MLS) approximation do not have Kronecker delta property, we cannot implement the essential boundary condition, directly. Thus we employ the shape functions of moving Kriging interpolation and radial point interpolation technique which have the mentioned property. Also, in the element free Galerkin method, we do not use any triangular, quadrangular or other type of meshes. The element free Galerkin method is a global method while finite elements method is a local one. This technique employs a background mesh for integration which makes it different from the truly mesh procedures. The coefficient matrix of the element free Galerkin is symmetric. Also, using numerical algorithms, we can conclude that the eigenvalues of the coefficient matrix are positive. Thus, for solving the obtained linear system of equations from the discretization, we use the conjugant gradient method. To keep away from solving a nonlinear algebraic system of equations and obtaining the acceptable numerical results, we use a predictor–corrector algorithm. Several test problems are solved and numerical simulations are reported which confirm the efficiency of the proposed schemes.