A Meshless Radial Basis Function Based on Partition of Unity Method for Piezoelectric Structures

Li, Sen; Yao, Linquan; Yi, Shichao; Wang, Wei

doi:https://doi.org/10.1155/2016/7632176

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusions Acknowledgments References Copyright Related Articles

Special Issue

Advanced Mesh-Based and Particle-Based Numerical Methods for Engineering and Applied Mathematics Problems

View this Special Issue

Research Article | Open Access

Volume 2016 | Article ID 7632176 | https://doi.org/10.1155/2016/7632176

A Meshless Radial Basis Function Based on Partition of Unity Method for Piezoelectric Structures

Sen Li,¹Linquan Yao,²Shichao Yi,³and Wei Wang²

Academic Editor: Lihua Wang

Received27 Aug 2016

Accepted20 Oct 2016

Published12 Dec 2016

Abstract

A meshless radial basis function based on partition of unity method (RBF-PUM) is proposed to analyse static problems of piezoelectric structures. The methods of radial basis functions (RBFs) possess some merits: the shape functions have high order continuity; -adaptivity is simpler than mesh-based methods; the shape functions are easily implemented in high dimensional space. The partition of unity method (PUM) easily constructs local approximation. The character of local approximate space can be varied and regarded as -adaptivity. Considering the good properties of the two methods, the RBFs are used for local approximation and the local supported weight functions are used in the partition of unity method. The system equations of the RBF-PUM are derived using the variational principle. The field variables are approximated using the RBF-PUM shape functions which inherit all the advantages of the RBF shape functions such as the delta function property. The boundary conditions can be implemented easily. Numerical examples of piezoelectric structures are investigated to illustrate the efficiency of the proposed method and the obtained results are compared with analytical solutions and available numerical solutions. The behaviors of some parameters that probably influenced numerical results are also studied in detail.

1. Introduction

Piezoelectric material is a new advanced material which can be used for developing smart structures with self-controlling and self-monitoring capabilities. The essential feature of piezoelectric material is the capability of energy transformation between electrical energy and mechanical energy. Hence, piezoelectric materials have been widely used as intelligent control devices such as sensors, actuators, ultrasonic transducers, and active damping devices.

Due to the complicated properties of piezoelectric materials, various approaches including analytical and numerical methods have been proposed to deal with the problems involving piezoelectric materials. However, the analytical solutions are rather restricted in general and difficult to find unless some geometry and boundary are relatively simple. Accurate, efficient, and robust numerical methods are required to simulate piezoelectric problems. Some numerical methods have been proposed to analyse piezoelectric problems such as the finite element method (FEM) [1–4], the boundary element method (BEM) [5–7], and the smoothed FEM [8]. Among them, the FEM is considered to be an effective method for mechanical and electromechanical coupling problems. Since the FEM is a mesh-based method, the FEM may lead to large errors when mesh distortion appears or low quality meshes are adopted. Mesh distortion may occur in the problems of large deformation, crack propagation, and so forth. An alternative method is called meshless method which can eliminate mesh-based problems. In the meshless method, only a set of scattered nodes is required to represent problem domain and boundary and therefore the adaptivity is simpler than mesh-based methods.

The initial ideal of meshless method dates back to the smooth particle hydrodynamics (SPH) method proposed by Lucy [9] and Gingold and Monaghan [10] for modeling astrophysical phenomena. Then several meshless methods including the diffuse element method (DEM) [11], the element-free Galerkin (EFG) method [12], the reproducing kernel particle method (RKPM) [13], the - clouds method [14], the point interpolation method (PIM) [15–17], the radial basis function point interpolation method (RPIM) [18, 19], the moving Kriging (MK) meshless method [20], and the meshless local Petrov-Galerkin (MLPG) method [21–25] have been proposed to solve partial differential equations and mechanical problems. Most of these methods are carried out based on weak form equations and a set of background cells or local subdomains is required to implement numerical integration. In the MLPG method, the moving least square (MLS) approach [26] is used to construct approximate functions and the major drawback is that the MLS shape functions lack interpolation property. The boundary conditions cannot be imposed directly and some methods should be proposed to overcome this problem, such as the Lagrange multiplier method [12] and the penalty function method [18]. In the RPIM, radial basis functions (RBFs) are used to construct the shape functions. The variable is only related to the distance between interpolation point and node and therefore the method is easily implemented in high dimension space. The behaviors of radial basis functions [27–32] have been widely investigated. Radial basis functions are also used in the meshless collocation method [33–35] which is based on the strong form equations. These methods based on the RBFs have obtained high accurate solutions and have been shown excellent interpolation property.

Some methods coupling with the partition of unity method (PUM) [36, 37] have also been investigated in mechanical problems in recent years. The main advantages of the PUM are the flexibility in choosing the local approximation which can be used the approaches of the MLS, the PIM, and the RPIM. Rajendran and Zhang [38] developed a new “FE-Meshfree” method using the concepts of partition of unity method. In this method, the FEM shape functions are used for weight functions and the meshless approach is used to construct local approximation. Some of important works including Q4-LSPIM [38, 39], Q4-RPIM [40], and T3-RPIM [41] have been presented for static and vibration analysis of mechanical problems and obtained very accurate results.

In this paper, a meshless radial basis function based on partition of unity method (RBF-PUM) is proposed to analyse 2D piezoelectric structures. The RBFs are used for local approximation and Shepard’s method [44] is used for the construction of weight functions in the partition of unity method. The RBF-PUM method has been studied by Wendland [45] and Fasshauer [29] and successfully applied in the convection diffusion equations [46] using the collocation method. The field variables are approximated using the RBF-PUM shape functions which inherits all good advantages of the PUM and RBFs. The performance of the RBF-PUM shape functions will be discussed in detail.

The rest of this paper is organized as follows. The formulation of the RBF-PUM method is presented in Section 2. The governing equations and variational equations for 2D piezoelectric problems are given in Section 3. Numerical examples are considered to verify the effectiveness of the proposed method in Section 4. Finally, we end this paper with some conclusions in Section 5.

2. The Formulation of the RBF-PUM Method

2.1. The Partition of Unity Method

This section recalls briefly the partition of unity method for 2D problem. Let be a problem domain and be an open cover of satisfying some mild overlapping conditions such thatwhere is an index set of the patches in which the interpolation point locates and is the total number of all patches.

For every patch , we construct a local approximation . According to the concept of the partition of unity method, the global approximation function can be expressed aswhere are weight functions. The weight functions are compactly supported, nonnegative, and continuous and satisfy the partition of unity:

Equation (3) shows that for . Therefore, (2) can be rewritten as

2.2. Construction of Weight Functions with Partition of Unity Property

For the partition of unity method, the patch can be of any shape, such as square, circular, spherical, or hyperspherical. The basic requirement for all patch domains is that they cover problem domain and boundary. In this paper, the circular patch is used for the analysis in the PUM.

Let be a set of central points of patches and be the radiuses as shown in Figure 1. Nonnegative and compactly supported weight functions can be constructed using Shepard’s method [44] as follows:where function is given byIn the present study, Wendland function [47] is used for function :in which

2.3. Radial Basis Function Method

In the present method, radial basis functions with polynomial terms are used to construct the local approximation of the partition of unity method. Using radial basis functions with polynomial functions, the local approximation is given bywhere and represent radial basis and polynomial functions, respectively. is the total number of nodes in patch and represents the number of polynomial terms. and are unknown vectors to be determined.

In the radial basis function , the variable is only related to the space distance between interpolation point and node and therefore the RBF method is easily implemented in high dimensional problems. Some radial basis functions such as multiquadric radial basis function (MQ-RBF), thin plate spline radial basis function (TPS-RBF), and Gaussian radial basis function (EXP-RBF) have been widely investigated in [18, 27–32]. For the present method, the MQ-RBF is adopted and given bywhere and and are two shape parameters.

Enforcing (9) to pass through all node values in patch , the following equations are obtained:where

Adding polynomial terms is an extra requirement. In order to guarantee unique approximation [48], the following conditions are usually imposed:which can be rewritten in matrix form

Vectors and can be obtained by (11) and (14), and then substitute into (9), the local approximation function is eventually expressed aswhere

The existence of inverse matrix for arbitrary scattered nodes has been given by Wendland [30]. Finally, the shape functions associated to all nodes in the patch can be expressed as

2.4. RBF-PUM Shape Functions and Mathematical Properties

The radial basis functions are selected as local approximation and the global approximation in (4) can be expressed aswhere is nodal displacement vector composed of all the nodes in the union of the patches and is the vector of the RBF-PUM shape functions given bywhere is the total number of nodes in domain . It is noted that the dimension of the RBF shape functions may be less than , but they can be expanded to vector with dimension .

The derivative functions of the RBF-PUM shape functions can be obtained using Leibniz’s rule:where are multi-indices and is a spatial differential operator.

The RBF-PUM method inherits the advantages of the RBF meshless method and the following are the mathematical properties of the RBF-PUM shape functions:(1)The RBF-PUM shape functions possess the delta function property:(2)The RBF-PUM shape functions possess the property of partition of unity:(3)If the RBFs include at least linear polynomial terms, the RBF-PUM shape functions can ensure an exact reproduction of linear polynomials; that is,(4)Local compactly supported property.(5)Especially for , the RBF-PUM shape functions are the RBF shape functions as well as their derivative functions.

The RBF-PUM shape functions and their derivative functions are shown in Figures 2 and 3 for 1D and 2D space, respectively. In Figure 2, six regularly distributed nodes in interval are used for the construction of the shape functions and plotted in different colours. In Figure 3, the shape function and its derivative of point are plotted in interval . It can be observed that the RBF-PUM shape functions have interpolation property and local compactly supported property.

3. Variational Equations for Piezoelectric Structures

The variational equations for linear piezoelectricity are presented in this section. Considering a two-dimensional piezoelectric problem in domain with boundary , the governing equations are given bywhere , are stress and electric displacement vector, respectively. is body force density. The operators and in plane are given by

The boundary conditions are given bywhere and denote prescribed displacements and electric potential, respectively. is the unit outward normal vector. and are the surface traction and charge, respectively.

The strain vector related to displacement isand the electric field vector iswhere

The constitutive relations for piezoelectric structures are given bywhere , , and are the elastic stiffness matrix, piezoelectric matrix, and dielectric matrix, respectively.

The constitutive relations can also be expressed explicitly in matrix form:

The mechanical and electrical constitutive relations can also be given bywhere is the elastic compliance matrix, is the piezoelectric strain matrix, and is the dielectric matrix. The relationships between these constant matrixes can be transforming each other.

The generalized energy functional is determined by a summation of strain energy, electric energy, and energy of external work as follows:

Using the generalized variational principle, the variational form of piezoelectric structures can be derived:Substituting (27)–(30) into (34), the linear piezoelectric equations can be obtained:where

4. Numerical Examples

Numerical examples for piezoelectric structures are used to study the RBF-PUM. The choice of the shape parameters and is an open issue and has been discussed in the RBF-based methods [18, 31]. When shape parameter is chosen as an integer, the moment matrix is ill-conditioned. Therefore, shape parameter cannot be an integer and, generally, the optimal parameters may change with the practical problems. The effects of the shape parameters on numerical results are also discussed in this paper. In some problems, unless otherwise stated, the two shape parameters and investigated the performance of the present method. Gauss integration scheme is used for all numerical examples: the boundary and each background cell adopt 4 and Gauss integration, respectively.

In the RBF-PUM, the patch fill distance is used as a measure of the density of central points and defined aswhich indicates how well the central points fill in problem domain. For regularly distributed central points, the patch fill distance is proportional to the distance of central points.

The support radius or the patch size is defined aswhere is a nondimensional coefficient. The patch size usually controls the number of nodes in each patch and therefore parameter can be adjusted according to node density and the patch fill distance.

For the purpose of error estimation and convergence studies, the relative error is defined aswhere “” and “” represent exact and numerical solutions, respectively.

4.1. Shear Deformation of Piezoelectric Strip

The shear deformation of a square piezoelectric strip polarized in -direction is taken into consideration. The piezoelectric strip is made of material PZT-5 and its material properties are given as follows: (mm)²/N (mm)²/N (mm)²/N (mm)²/N N/(mm)² N/(mm)² mm mm/V mm/V mm/V N/V² N/V² V mm

The strip is subjected a uniform load on the top and bottom boundaries and an applied voltage on the left and right boundaries as shown in Figure 4. Since the applied electric field is perpendicular to polarization direction, it induces a shear strain. The overall deformation of the piezoelectric strip is a superposition of the deformation due to the shear strain and the applied load.

The following boundary conditions of the piezoelectric strip are given by

The analytical solutions for this shear problem are given by [49]

The problem domain is represented by both regular nodes and 121 irregular nodes in Figure 5. The PU cover consists of regular circular patches and the coefficient is taken here. square background cells are used for numerical integration.

(a)

(b)

The obtained results of the present method are plotted in Figure 6, in which symbols and denote regular and irregular nodes, respectively. Figure 6(a) shows the horizontal displacement along and the deflection along is plotted in Figure 6(b). The variation of the electric potential is depicted in Figure 6(c). The results of the RBF method are also presented in these figures and the size of influence domain of the RBF method is taken as 4.0 times average nodal distance. The results of the RBF-PUM method and the RBF method are very accurate comparing with the analytical solutions. Moreover, the present results are slightly better than those of the RBF method.

(a)

(b)

(c)

The computed mechanical deformation results of the present method for irregular nodes are plotted in Figure 7, which seems quite reasonable for the shear deformation. The convergence curves and convergence rates of displacements and and electric potential using the linear regression are shown in Figure 8 for regular nodes. It can be seen that the present results are very accurate even for coarse node distribution. The relative errors will decrease with the increase of nodes. It is due to the fact that the number of nodes using in the RBF-PUM shape functions increases with the node density and the numerical results become better.

In the RBF-PUM, the common requirement for all patches is that they cover the problem domain and boundary. Numerical results are usually influenced by the patch size, since the number of nodes using in the construction of the RBF-PUM shape function are influenced by the patch size. It is great of interest to study the effects of the patch size on numerical results. The results of displacements and electric potential at point (, ) are listed in Table 1. The present method gives stable and accurate results for different patch sizes and nodal patterns. The solutions will converge to the analytical solutions for a larger patch size. It is due to the fact that increasing the patch size will lead to an increase of the number of nodes in each patch and the RBF-PUM shape function for a given interpolation point contains more data information for a larger patch size. However, it should be noted that the optimal parameter may change with the patch fill distance and the coefficient can be adjusted according to the patch fill distance.

4.2. Bending Deformation of Piezoelectric Strip

The same material and geometry of the piezoelectric strip in bending are studied. The top and bottom surfaces are subjected to an applied voltage and the right side is applied to a linear stress as depicted in Figure 9. The analytical solutions of displacements and electric potential are given by [50]and the boundary conditions for this problem are given by

To demonstrate the efficiency of the proposed method for irregular central points, 25 circular patches with irregular central points are used for the bending analysis of the piezoelectric strip. The central points and nodes are plotted in Figure 10, in which the patch fill distance is about and the coefficient is also taken here.

(a)

(b)

The solutions of the RBF-PUM, the RBF method, and the analytical method are shown in Figure 11. The displacement and electric potential along are plotted in Figures 11(a) and 11(c), respectively. The variation of the deflection along the bottom of the piezoelectric strip is depicted in Figure 11(b). It can be observed that the present results match the analytical solutions very well even for irregular central points and the present results are more accurate than the RBF results by comparing with the analytical solutions.

(a)

(b)

(c)

4.3. A Parallel Piezoelectric Bimorph Beam

Consider a two-layer parallel piezoelectric bimorph beam as shown in Figure 12. The beam is made of material PVDF with the same thickness and polarization direction and the material properties are summarized as follows: Pa C/m² C/m² F/m F/m

For this problem, a unity voltage (1 V) is applied across the thickness between the surfaces and the intermediate electrode. In this study, a plane stress problem is assumed. For an applied voltage , the tip deflection can be approximated as [51]where is the length of the beam, is the total thickness, and is the piezoelectric strain constant. The size of the beam is taken as mm and mm and the tip deflection computed by (44) is m.

Considering the thickness-length ratio of the beam, a total of regular patches and regular nodes are used to represent the problem domain. Due to its geometric symmetry, only one-half of the patches of the top beam are shown in Figure 13, where the central points are at the upper and down quarter through the thickness direction and the coefficient is taken.

The tip deflections for different node densities are listed in Table 2 and the present results are very close to the results given by the ABAQUS [19]. The effect of the patch size on the tip deflection is investigated in Table 3. As the patch size increases, the number of nodes in each patch will increase as well as the nodes using in the construction of the RBF-PUM shape functions. The obtained solutions will converge for a larger patch size. The shape parameters and are also studied in Tables 4 and 5, respectively. It can be observed that increasing the parameter leads to an increase of the tip deflection, but it changes very slowly and the tip deflection keeps a constant for different values of parameter . All the results computed by present method are very stable and slightly larger than the analytical solution. It may be due to using the equations derived by the Euler-Bernoulli model and the numerical method results in a less stiff model [19].

4.4. Piezoelectric Cook’s Membrane

The piezoelectric Cook’s membrane problem is solved using the present method. The geometry and applied unity tip load are shown in Figure 14. The lower surface of the membrane is prescribed by zero voltage. Cook’s membrane is made of PZT-4 material provided as follows: N/mm² N/mm² N/mm² N/mm² pC/mm² pC/mm² pC/mm² pC/GVmm pC/GVmm

The analytical solutions of this problem are unavailable and the best known numerical results are given by Long et al. [42] using the finite element method with a fine mesh. The results of the displacement, the electric potential, the first principal stress, and the electric flux density at nodes , , are

The problem is represented by circular patches and regular nodes as shown in Figure 15. The patch fill distance is about and the coefficient is . The results and relative errors percentage of the present method for different nodal densities are listed in Table 6 together with the results of the smoothed finite element method [8]. In the smoothed FEM, a fine mesh with elements is used for this problem. It can be seen that the present results of the displacement and the first principal stress are in very good accord with the results given by the FEM. Although the errors of the electric potential and the electric flux density are slightly lager, the present method still achieves reasonable prediction.

4.5. A Piezoelectric Bimorph Optical Microscanner

The purpose of this problem is to simulate the linear title angle of the reflected light through a mirror of an optical microscanner. The device is composed of two parallel bimorphs connected by a mirror at the tip center as depicted in Figure 16. Each bimorph has length μm and height μm and the length of the mirror is μm. The two-layer piezoelectric bimorph is made of PVDF material and its material properties are summarized as follows: N/μm² N/μm² N/μm² N/μm² N/Vμm N/Vμm N/V² N/V²

When a voltage is applied on the surfaces, the bimorphs deflect in opposite direction so that the mirror moves up and down, resulting in the change of the tilt angle . The direction of the reflected light can be changed with applied different voltages and the tip deflection can be calculated by the present method. For sufficient small rotation, the angle can be approximated from the tip deflection

The boundary conditions of the bottom layer areAnd the boundary conditions of the top layer are

Considering the geometry of the optical microscanner device, node arrangements and patches are very similar to the parallel piezoelectric bimorph beam. The central points are at the upper and down quarter through the thickness direction. patches and are used for this problem.

The tip displacements calculated for several applied voltages with different nodal densities are listed in Table 7, in which the first two rows represent background cells and nodal densities, respectively. The numerical results obtained by the meshless point collocation method (PCM) [43] are also given in the table. In the PCM, the reproducing kernel particle approximation is used for the displacement field and electric potential. The governing equations and boundary conditions are satisfied only at specific nodes. In the present method, the system equations are derived by the variational principle and numerical integration is required. Although the present results are very close to the PCM solutions, there are still some differences between these two different numerical methods.

From the tip deflection, the title angle of the mirror can be calculated by (46). The title angles varying with different applied voltages for regular nodes are shown in Figure 17. The results obtained by the node-based element (NSPE-T3, NSPE-Q4), the cell-based element (SPQ4), and the FEM-T3 [52] are also shown in Figure 17. The linear variation of the angle is in good agreement with the PCM.

5. Conclusions

In this paper, a meshless radial basis function based on partition of unity method is presented for studying 2D piezoelectric structures. The multiquadric radial basis functions are used for local approximation and Shepard’s method is used for the construction of weight functions in the partition of unity method. The weak form equations are derived by the generalized variational principle. The mechanical displacements and electric potential are approximated using the RBF-PUM shape functions which inherit all advantages of the RBF shape functions. The RBF-PUM shape functions have been shown excellent interpolation property and therefore the boundary conditions can be imposed easily without special treatments.

Numerical examples are calculated to illustrate the efficiency of the proposed method for static problems. The effects of some parameters including the shape parameters, the patch size, and nodal density have also been studied in detail. The obtained results are compared with the solutions of the analytical method and other numerical methods such as the RBF method, the finite element method, and the smoothed finite element method. The present results are very accurate and are in excellent agreement with the analytical solutions.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (nos. 11172192 and 11572210) and the Graduate Innovation Project of Jiangsu Province (KYLX16_0095).

References

C. D. Mercer, B. D. Reddy, and R. A. Eve, “Finite element method for piezoelectric media,” Applied Mechanics Research Unit Technical Report 92, University of Cape Town/CSIR, 1987.
View at: Google Scholar
H. Allik and T. J. R. Hughes, “Finite element method for piezo-electric vibration,” International Journal for Numerical Methods in Engineering, vol. 2, no. 2, pp. 151–157, 1970.
View at: Publisher Site | Google Scholar
K. Y. Sze, X.-M. Yang, and L.-Q. Yao, “Stabilized plane and axisymmetric piezoelectric finite element models,” Finite Elements in Analysis and Design, vol. 40, no. 9-10, pp. 1105–1122, 2004.
View at: Publisher Site | Google Scholar
J. Mackerle, “Smart materials and structures—a finite element approach—an addendum: a bibliography (1997–2002),” Modelling and Simulation in Materials Science and Engineering, vol. 11, no. 5, pp. 707–744, 2003.
View at: Publisher Site | Google Scholar
F. Garcła-Sánchez and C. Zhang, “2D transient dynamic analysis of cracked piezoelectric solid by a timedomain BEM,” Computer Methods in Applied Mechanics and Engineering, vol. 197, pp. 3108–3121, 2008.
View at: Google Scholar
A. Sáez, F. García-Sánchez, and J. Dominguez, “Hypersingular BEM for dynamic fracture in 2D piezoelectric solids,” Computer Methods in Applied Mechanics and Engineering, vol. 196, pp. 235–246, 2006.
View at: Google Scholar
Ch. Zhang, F. García-Sánchez, and A. Sáez, “Time-domain BEM analysis of cracked piezoelectric solids under impact loading,” in Computational Mechanics, pp. 206–218, Springer, 2007.
View at: Google Scholar
H. Nguyen-Van, N. Mai-Duy, and T. Tran-Cong, “A smoothed four-node piezoelectric element for analysis of two-dimensional smart structures,” Computer Modeling in Engineering and Sciences, vol. 23, no. 3, pp. 209–222, 2008.
View at: Google Scholar
L. B. Lucy, “A numerical approach to the testing of the fission hypothesis,” The Astronomical Journal, vol. 82, pp. 1013–1024, 1977.
View at: Publisher Site | Google Scholar
R. A. Gingold and J. J. Monaghan, “Smoothed particle hydrodynamics: theory and application to non-spherical stars,” Monthly Notices of the Royal Astronomical Society, vol. 181, no. 3, pp. 375–389, 1977.
View at: Google Scholar
B. Nayroles, G. Touzot, and P. Villon, “Generalizing the finite element method: diffuse approximation and diffuse elements,” Computational Mechanics, vol. 10, no. 5, pp. 307–318, 1992.
View at: Publisher Site | Google Scholar
T. Belytschko, Y. Y. Lu, and L. Gu, “Element-free Galerkin methods,” International Journal for Numerical Methods in Engineering, vol. 37, no. 2, pp. 229–256, 1994.
View at: Publisher Site | Google Scholar | MathSciNet
W. K. Liu, S. Jun, and Y. F. Zhang, “Reproducing kernel particle methods,” International Journal for Numerical Methods in Fluids, vol. 20, no. 8-9, pp. 1081–1106, 1995.
View at: Publisher Site | Google Scholar
C. A. Duarte and J. T. Oden, “An $h$ - $p$ adaptive method using clouds,” Computer Methods in Applied Mechanics and Engineering, vol. 139, no. 1–4, pp. 237–262, 1996.
View at: Publisher Site | Google Scholar | MathSciNet
G. R. Liu and Y. T. Gu, “A point interpolation method for two-dimensional solids,” International Journal for Numerical Methods in Engineering, vol. 50, no. 4, pp. 937–951, 2001.
View at: Publisher Site | Google Scholar
G. R. Liu and Y. T. Gu, “Local point interpolation method for stress analysis of two-dimensional solids,” Structural Engineering and Mechanics, vol. 11, no. 2, pp. 221–236, 2001.
View at: Publisher Site | Google Scholar
G. R. Liu, K. Y. Dai, K. M. Lim, and Y. T. Gu, “A point interpolation mesh free method for static and frequency analysis of two-dimensional piezoelectric structures,” Computational Mechanics, vol. 29, no. 6, pp. 510–519, 2002.
View at: Publisher Site | Google Scholar
G. R. Liu and Y. T. Gu, An Introduction to Meshfree Methods and Their Programming, Springer, Amsterdam, The Netherlands, 2005.
G. R. Liu, K. Y. Dai, K. M. Lim, and Y. T. Gu, “A radial point interpolation method for simulation of two-dimensional piezoelectric structures,” Smart Materials and Structures, vol. 12, no. 2, pp. 171–180, 2003.
View at: Publisher Site | Google Scholar
T. Q. Bui, M. N. Nguyen, C. Zhang, and D. A. K. Pham, “An efficient meshfree method for analysis of two-dimensional piezoelectric structures,” Smart Materials and Structures, vol. 20, no. 6, Article ID 065016, pp. 1–16, 2011.
View at: Publisher Site | Google Scholar
S. N. Atluri and T. Zhu, “A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics,” Computational Mechanics, vol. 22, no. 2, pp. 117–127, 1998.
View at: Publisher Site | Google Scholar | MathSciNet
J. Sladek and V. Sladek, “Meshless local Petrov-Galerkin method for plane piezoelectricity,” CMC-Computers Materials & Continua, vol. 4, pp. 109–118, 2006.
View at: Google Scholar
J. Sladek, V. Sladek, and C. Zhang, “Dynamic crack analysis in functionally graded piezoelectric solids by meshless local petrov-galerkln method,” Key Engineering Materials, vol. 348-349, pp. 149–152, 2007.
View at: Publisher Site | Google Scholar
J. Sladek, V. Sladek, C. Zhang, P. Solek, and L. Starek, “Fracture analyses in continuously nonhomogeneous piezoelectric solids by the MLPG,” Computer Modeling in Engineering and Sciences, vol. 19, no. 3, pp. 247–262, 2007.
View at: Google Scholar
J. Sladek, V. Sladek, C. Zhang, and P. Solek, “Application of the MLPG to thermo-piezoelectricity,” CMES-Computer Modeling in Engineering and Sciences, vol. 22, no. 3, pp. 217–233, 2007.
View at: Google Scholar
P. Lancaster and K. Salkauskas, “Surfaces generated by moving least squares methods,” Mathematics of Computation, vol. 37, no. 155, pp. 141–158, 1981.
View at: Publisher Site | Google Scholar | MathSciNet
E. J. Kansa, “Multiquadrics-A scattered data approximation scheme with applications to computational fluid-dynamics-I surface approximations and partial derivative estimates,” Computers and Mathematics with Applications, vol. 19, no. 8-9, pp. 127–145, 1990.
View at: Publisher Site | Google Scholar
E. J. Kansa, “Multiquadrics—a scattered data approximation scheme with applications to computational fluid-dynamics—II solutions to parabolic, hyperbolic and elliptic partial differential equations,” Computers and Mathematics with Applications, vol. 19, no. 8-9, pp. 147–161, 1990.
View at: Publisher Site | Google Scholar
G. E. Fasshauer, Meshfree Approximation Methods with MATLAB, World Scientific Publishing, 2007.
H. Wendland, “Error estimates for interpolation by compactly supported radial basis functions of minimal degree,” Journal of Approximation Theory, vol. 93, no. 2, pp. 258–272, 1998.
View at: Publisher Site | Google Scholar | MathSciNet
J. G. Wang and G. R. Liu, “On the optimal shape parameters of radial basis functions used for 2-D meshless methods,” Computer Methods in Applied Mechanics and Engineering, vol. 191, no. 23-24, pp. 2611–2630, 2002.
View at: Publisher Site | Google Scholar | MathSciNet
R. Schaback and H. Wendland, Characterzation and Construction of Radial Basis Functions, Cambridge University Press, Cambridge, UK, 2001.
X. Zhang, K. Z. Song, M. W. Lu, and X. Liu, “Meshless methods based on collocation with radial basis functions,” Computational Mechanics, vol. 26, no. 4, pp. 333–343, 2000.
View at: Publisher Site | Google Scholar
H. Y. Hu, J. S. Chen, and W. Hu, “Weighted radial basis collocation method for boundary value problems,” International Journal for Numerical Methods in Engineering, vol. 69, no. 13, pp. 2736–2757, 2007.
View at: Publisher Site | Google Scholar | MathSciNet
F. Chu, L. Wang, and Z. Zhong, “Finite subdomain radial basis collocation method,” Computational Mechanics, vol. 54, no. 2, pp. 235–254, 2014.
View at: Publisher Site | Google Scholar | MathSciNet
J. M. Melenk and I. Babuska, “The partition of unity finite element method: basic theory and applications,” Computer Methods in Applied Mechanics and Engineering, vol. 139, no. 1–4, pp. 289–314, 1996.
View at: Publisher Site | Google Scholar | MathSciNet
I. Babuska and J. M. Melenk, “The partition of unity method,” International Journal for Numerical Methods in Engineering, vol. 40, no. 4, pp. 727–758, 1997.
View at: Publisher Site | Google Scholar | MathSciNet
S. Rajendran and B. R. Zhang, “A ‘FE-meshfree’ QUAD4 element based on partition of unity,” Computer Methods in Applied Mechanics and Engineering, vol. 197, no. 1–4, pp. 128–147, 2007.
View at: Publisher Site | Google Scholar
B. R. Zhang and S. Rajendran, “'FE-Meshfree' QUAD4 element for free-vibration analysis,” Computer Methods in Applied Mechanics and Engineering, vol. 197, no. 45–48, pp. 3595–3604, 2008.
View at: Publisher Site | Google Scholar
J. P. Xu and S. Rajendran, “A partition-of-unity based ‘FE-Meshfree’ QUAD4 element with radial-polynomial basis functions for static analyses,” Computer Methods in Applied Mechanics and Engineering, vol. 200, no. 47-48, pp. 3309–3323, 2011.
View at: Publisher Site | Google Scholar
Y. Yang, D. Xu, and H. Zheng, “A partition-of-unity based ‘FE-meshfree’ triangular element with radial-polynomial basis functions for static and free vibration analysis,” Engineering Analysis with Boundary Elements, vol. 65, pp. 18–38, 2016.
View at: Publisher Site | Google Scholar | MathSciNet
C. S. Long, P. W. Loveday, and A. A. Groenwold, “Planar four node piezoelectric elements with drilling degrees of freedom,” International Journal for Numerical Methods in Engineering, vol. 65, no. 11, pp. 1802–1830, 2006.
View at: Publisher Site | Google Scholar
R. R. Ohs and N. R. Aluru, “Meshless analysis of piezoelectric devices,” Computational Mechanics, vol. 27, no. 1, pp. 23–36, 2001.
View at: Publisher Site | Google Scholar
D. Shepard, “A two-dimensional interpolation function for irregularly spaced points,” in Proceedings of the 23rd ACM National Conference (ACM '68), pp. 517–24, 1968.
View at: Publisher Site | Google Scholar
H. Wendland, “Fast evaluation of radial basis functions: methods based on partition of unity,” in Approximation Theory X Wavelets Splines & Applications, pp. 473–483, 2002.
View at: Google Scholar
A. S. Vaighani, A. Heryudono, and E. Larsson, “A radial basis function partition of unity collocation method for convection-diffusion equations arising in financial applications,” Journal of Scientific Computing, vol. 64, no. 2, pp. 341–367, 2015.
View at: Publisher Site | Google Scholar
H. Wendland, “Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree,” Advances in Computational Mathematics, vol. 4, no. 4, pp. 389–396, 1995.
View at: Publisher Site | Google Scholar | MathSciNet
M. A. Golberg, C. S. Chen, and H. Bowman, “Some recent results and proposals for the use of radial basis functions in the BEM,” Engineering Analysis with Boundary Elements, vol. 23, no. 4, pp. 285–296, 1999.
View at: Publisher Site | Google Scholar
V. Z. Paarton and G. K. Mikhailov, Applied Mechanics: Soviet Reviews, vol. 2 of Electromagneoelasticity, Hemisphere, 1989.
P. Gaudnzi and K.-J. Bathe, “An iterative finite element procedure for the analysis of piezoelectric continua,” Journal of Intelligent Material Systems and Structures, vol. 6, no. 2, pp. 266–273, 1995.
View at: Publisher Site | Google Scholar
Product Catalogue, Piezo Systems Inc., Cambridge, Mass, USA, 1995.
H. Nguyen-Van, N. Mai-Duy, and T. Tran-Cong, “A node-based element for analysis of planar piezoelectric structures,” Computer Modeling in Engineering & Sciences, vol. 36, no. 1, pp. 65–96, 2008.
View at: Google Scholar | MathSciNet

Copyright

Copyright © 2016 Sen Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1886

Downloads

1038

Citations