Research Article  Open Access
ChiaCheng Tsai, "Polyharmonic Multiquadric Particular Solutions for Reissner/Mindlin Plate", Mathematical Problems in Engineering, vol. 2015, Article ID 246159, 12 pages, 2015. https://doi.org/10.1155/2015/246159
Polyharmonic Multiquadric Particular Solutions for Reissner/Mindlin Plate
Abstract
Analytical particular solutions of the polyharmonic multiquadrics are derived for both the Reissner and Mindlin thickplate models in a unified formulation. In the derivation, the three coupled secondorder partial differential equations are converted into a product operator of biharmonic and Helmholtz operators using the Hörmander operator decomposition technique. Then a method is introduced to eliminate the Helmholtz operator, which enables the utilization of the polyharmonic multiquadrics. Then, the analytical particular solutions of displacements, shear forces, and bending or twisting moments corresponding to the polyharmonic multiquadrics are all explicitly derived. Numerical examples are carried out to validate these particular solutions. The results obtained by the present method are more accurate than those by the traditional multiquadrics and splines.
1. Introduction
Boundarytype numerical methods have been emerged as a popular research field because only boundary discretizations are required when they are applied to solve homogeneous problems. These methods include the boundary element method (BEM) [1], method of fundamental solutions (MFS) [2, 3], and Trefftz methods (TM) [4, 5]. When they are applied for solving an inhomogeneous partial differential equation, the dual reciprocity method (DRM) can be used [6, 7].
In the early development of DRM, the ad hoc radial function, , was exclusively used. In order to improve the accuracy of the computation, Golberg and Chen [8] and Karur and Ramachandran [9] applied the theory of radial basis functions (RBFs) to the DRM. Among these RBFs, Hardy’s multiquadrics (MQ) [10] and Duchon’s augmented polyharmonic spline (APS) [11] are the most popular ones. For example, Golberg [12, 13], Chen [14], and Karur and Ramachandran [9] demonstrated the superiority of the APS over the ad hoc radial function. Then, Golberg et al. [15] further improved the accuracy of the approximated particular solution by utilizing the MQ. Recently, Tsai [16] generalized the MQ to the polyharmonic multiquadrics (PMQ) and showed further accuracy improvement of the PMQ over the traditional MQ. In this study, the accuracy improvement of the PMQ over the MQ and APS is demonstrated for the Reissner and Mindlin thickplate models.
When the MQ is adopted in a DRM procedure, the applicability depends on the availability of the analytical particular solution of the basis function associated with the partial differential operator of a given problem. Golberg et al. [15] derived the analytical particular solutions for the Laplace equation. Samaan and Rashed [17, 18] and Tsai and Hsu [19] found the analytical particular solutions, respectively, for the two and threedimensional elasticity problems, which can be converted to a biharmonic equation using the GalerkinPapkovich vector [20]. Basically, the applicability of MQ is limited to the harmonic and biharmonic operators [21] before Tsai [16] derived the analytical particular solutions of MQ associated with the polyharmonic operators, which are denoted as the PMQ in the following.
When applying the DRM for solving the Mindlin thickplate model [22], Wen et al. [23] and Tsai and Wu [24] found the analytical particular solutions for the ad hoc RBF and the APS, respectively. In addition, the analytical particular solutions for the Chebyshev polynomials [25] and the APS [26] were derived for the Reissner thickplate model [27, 28]. Basically, a product operator of biharmonic and Helmholtz operators resulted when applying the Hörmander operator decomposition technique [29] to the coupled PDEs of the Reissner and Mindlin thickplate models as demonstrated in the previous studies [24–26]. This seems to infer that the accuracy improvement of the PMQ over the APS is not possible for the thickplate models.
In this study, a method is introduced to eliminate the Helmholtz operator, which makes the application of the PMQ for the thickplate models possible. Numerical experiments are carried out to demonstrate the superior accuracy of the PMQ over the traditional MQ and the first few orders of the APS.
This paper is organized as follows: the unified Reissner/Mindlin plate model is introduced in Section 2. Then, a review of the PMQ is given in Section 3. And the particular solutions of the PMQ associated with the Reissner/Mindlin plate model are derived in Section 4. Section 5 gives the formulation of MFSDRM. Finally, some numerical experiments are carried out to validate the particular solutions in Section 6 and the conclusions are drawn in Section 7.
2. Reissner/Mindlin Plate Model
In the following, indices are in the range and indices are in the range . As described in Figure 1, we consider a plate of uniform thickness with its middle plane being a domain with boundary in the plane and thickness coordinate . The plate is subjected to a normal loading with intensity . Then, the equations of equilibrium are given by where is the transverse shear force and is the bending or twisting moment.
In the Reissner [27, 28] and Mindlin [22] thickplate theory, represents an average slope angle and is the lateral deflection of the plate in the middle surface. According to the theories, the constitutive equations are given by with for the Mindlin plate model and for the Reissner plate model. Combining (1)~(2) results in the following governing equation: withwhere and are Young’s modulus and Poisson ratio, respectively.
Equation (5) is a coupled system of three secondorder PDEs for the three unknown functions , and, therefore, three boundary conditions are required as where are given boundary data.
3. Polyharmonic Multiquadrics
Before deriving the particular solutions of PMQ associated with the Reissner/Mindlin plate model, we give a short review on the PMQ previously derived by the method of undetermined coefficients together with the Laurent series in [16]. Considering the th order PMQ , it is governed by the following hierarchical relation: Here, (10) is introduced so that the PMQ are hierarchically unique and is a shape parameter to tune the numerical accuracy [10].
Then, by using the method of undetermined coefficients, we can find the solution of (8)~(10) as with for any positive integer . And, the unknown in (11) can be solved by using the uniqueness condition (10) as with In (14), the double factorial is defined as The PMQ defined in (11) can be expanded into the Maclaurin series as Observing (16), we can conclude that the proposed PMQ are infinitely differentiable since they only consist of evencoefficient polynomials [16, 30]. Then, substituting (16) into (8) and (9) and performing some mathematical manipulations can result in which can be used for solving the unknown coefficients and in (11) and (16). In practice, we simply enforce equal to in (17) which results in linear equations for solving the unknown coefficients and of the PMQ. After and are obtained, (13) can be used to have . Then, the PMQ can be computed simply by (11). For the cases when is very small or equal to zero, one can alternatively use the Maclaurin series (16).
For the purpose of a later derivation of the particular solutions associated with the Reissner/Mindlin plate model, we need to find the formulas when the radial differential operator is applied onto (11) and (16). This can be done for arbitrary nonnegative integer since the proposed PMQ are infinitely differentiable. Then, the required formulas for (11) and (16) can be derived, respectively, as with In deriving (18), we have used And, in deriving (20) and (21), we have used the product rule of differentiation as follows: with and being two arbitrary radial functions. Finally, additional formulas are required to complete (20) and (21), respectively, as follows: with
The proposed formulas can be implemented into a floatingpoint subroutine for evaluating with arbitrary nonnegative integers and . Therefore, the subroutine can be invoked when implementing numerical methods. In addition, the source codes can be obtained from the author via emails.
4. Analytical Particular Solution of Plate Model
Analytical particular solutions are required when the DRM is applied for solving a Reissner/Mindlin plate problem under an arbitrary loading. In this section, the analytical particular solutions of the PMQ associated with the Reissner/Mindlin plate model will be derived, which are governed by where
Formally, the particular solutions , governed by (27), should be obtainable by using the Hörmander operator decomposition technology [29], which begins with the definition of the adjoint operator as where is the determinant of and is the threebythree identity matrix. Here, both and can be obtained by linear algebraic manipulations, respectively, defined asThen, the particular solutions can be found by assuming where is an unknown function to be determined. In (32), the required operator can be obtained as
After substituting (32) into (27) and using (29) and (30), the resulted equation indicates that the following particular solution is required: However, the analytical particular solution of the above equation is not known due to the existence of the Helmholtz operator.
In this study, we find that the Helmholtz operator can be eliminated by rewriting (33) as where the operator is defined by Then, (29) and (35) can be combined to have The above equations have suggested us to modify the Hörmander operator decomposition technology by assuming where is an unknown function to be determined. Substituting (38) into (27) and using (37) result in By using the PMQ definitions (8) and (9), we can have the particular solution of (39) as Substituting (40) into (38) and using (36) can result in the desired analytical particular solutions of displacements as where we have used In addition, the required formulas of and have been provided in the previous section. Then, the corresponding shear forces and bending/twisting moments can be obtained by using the definitions Equations (41), (42), (44), and (45) are sufficient to obtain the analytical particular solutions of shear forces and bending/twisting moments, respectively, as follows: where the additional terms and have been given in the previous section.
This completes the derivation of the analytical particular solutions of the PMQ associated with the Reissner/Mindlin plate model.
5. MFSDRM Formulation
Now, we are in a position to review the application of the MFSDRM procedure [24, 26] for solving the wellposed thickplate problem governed by (5) and (7). First of all, we need the principle of superposition as where the particular solution satisfies without specifying any boundary condition. In addition, the homogeneous solution satisfies with modified boundary conditions where and .
In order to apply the DRM for solving the particular solution , the loading intensity needs to be first approximated by the PMQ as follows: where is the distance between and the th DRM field point, , as depicted in Figure 2. In order to solve the unknown coefficient , (52) should be collocated on the DRM field points as for and is the distance between the th and th DRM field points. Then, the corresponding particular solution can be approximated by with the particular solutions given by (41) and (42). Also, the particular solutions of shear force and bending or twisting moment can be approximated, respectively, as In (55) and (56), the particular solutions and are given in (46) and (47), respectively.
After the particular solutions are approximated by the DRM, the modified boundary conditions (51) become well defined. Therefore, it is ready to use the MFS for solving the homogeneous solution. Formally, the homogeneous solution can be approximated by where are the fundamental solutions defined byand are source points outside the plate domain as described in Figure 3. In (57), the unknowns can be solved by collocating the modified boundary conditions (51) on boundary points. More details of the MFS and the explicit fundamental solutions can be found in [24, 26]. After both the particular and homogeneous solutions are solved, the principle of superposition (48) can be applied for obtaining the desired solution.
6. Numerical Results
In order to validate the derived particular solutions, three numerical cases are considered. Typically, and are selected in these numerical studies. Furthermore, the MFS is always arranged to be sufficient and accurate such that the numerical errors are majorly controlled by the DRM. The root mean square errors (RMSEs) of the numerical solutions are defined as where is the numerical solution obtained by the MFSDRM at the th point of the considered positions and is the corresponding analytical solution of the problem. Here, will be chosen to be sufficiently large to ensure the utilization of the RMSEs.
In our numerical experiments, we consider the following Reissner/Mindlin plate problem: This problem has been considered for both the Mindlin [24] and the Reissner [26] cases and their analytical particular solutions can be unified as
Example 1. First, we consider a Reissner plate defined in subjected to clamped boundary conditions, which are set according to the analytical solutions in (61).
The particular solutions of the present problem are solved by the DRM based on the APS [26], Chebyshev polynomials [25], and PMQ. Table 1 gives the RMSEs for the solutions obtained by the DRMs with different bases and . Basically, the solutions obtained by the Chebyshev polynomials are the most accurate and the application is straightforward if the arbitrary loadings are given by functions. However, when the loadings are measured on some scattered points, a twostage method [31] is required which is relatively more difficult compared to the direct application of either APS or PMQ. In the table, the superior accuracy of the MQ over the first few orders of the APS can be observed. In addition, the accuracy improvement for the numerical solutions obtained by the highorder PMQ over the traditional MQ can also be observed. Overall, the best accuracy of the highorder PMQ is close to that of the highorder APS while the former RBF does not need augmented polynomials. In practice, it can be cumbersome to implement the particular solutions of the highorder augmented polynomials and thus the superiority of the PMQ over the APS becomes more relevant.

For both MQ and PMQ, the optimal shape parameters are searched as demonstrated in Figure 4. Basically, the numerical accuracy improves with increasing value. However, as is continuously made larger, the interpolation matrix becomes illconditioned [32]. Then, a typical error contour is given in Figure 5 for the solution obtained by the MFSDRM based on the MQ. In the figure, the larger errors are on the boundary as expected.
(a)
(b)
Finally, the RMSE comparison for different numbers of nodes is given in Figure 6, in which higher resolution gives better accuracy. Overall, these results should have validated the polyharmonic multiquadric particular solutions of .
Example 2. Then, we consider a similar Mindlin plate problem by changing the boundary conditions to the free boundary condition on one side. Table 2 addresses the RMSEs for the numerical results in this example. The results are very similar except that the errors are slightly larger due to the numerical approximations on the free boundary condition.

In addition, the optimal shape parameter searches and the error comparison for different numbers of nodes are given in Figures 7 and 8, respectively. The results are similar to the previous example. Therefore, the polyharmonic multiquadric particular solutions of and are also validated.
(a)
(b)
7. Discussions
In this paper, analytical particular solutions of the polyharmonic multiquadrics associated with Reissner and Mindlin thick plates in bending are derived. In the derivations, the coupled system of the three secondorder governing equations is transformed into a single sixthorder product operator of biharmonic and Helmholtz operators. Then, a method is introduced to eliminate the Helmholtz operator, which enables the use of the polyharmonic multiquadric. Numerical results are performed to validate the derived particular solutions and to demonstrate the superiority of the polyharmonic multiquadrics over the traditional multiquadrics and the augmented polynomial splines.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The National Science Council of Taiwan is gratefully acknowledged for providing financial support to carry out the present work under Grant no. MOST 1032221E022017MY2.
References
 A. H.D. Cheng and D. T. Cheng, “Heritage and early history of the boundary element method,” Engineering Analysis with Boundary Elements, vol. 29, no. 3, pp. 268–302, 2005. View at: Publisher Site  Google Scholar
 G. Fairweather and A. Karageorghis, “The method of fundamental solutions for elliptic boundary value problems,” Advances in Computational Mathematics, vol. 9, no. 12, pp. 69–95, 1998. View at: Publisher Site  Google Scholar  MathSciNet
 M. A. Golberg and C. S. Chen, Eds., The Method of Fundamental Solutions for Potential, Helmholtz and Diffusion Problems, Computational Mechanics Publications, Southampton, UK, 1999.
 C.S. Liu, “A highly accurate solver for the mixedboundary potential problem and singular problem in arbitrary plane domain,” CMES: Computer Modeling in Engineering and Sciences, vol. 20, no. 2, pp. 111–122, 2007. View at: Google Scholar
 C.S. Liu, “A highly accurate MCTM for direct and inverse problems of biharmonic equation in arbitrary plane domains,” CMES: Computer Modeling in Engineering and Sciences, vol. 30, no. 2, pp. 65–75, 2008. View at: Google Scholar
 P. W. Partridge, C. A. Brebbia, and L. C. Wrobel, The Dual Reciprocity Boundary Element Method, Computational Mechanics, Southampton, UK, 1992. View at: MathSciNet
 D. Nardini and C. A. Brebbia, “A new approach to free vibration analysis using boundary elements,” Applied Mathematical Modelling, vol. 7, no. 3, pp. 157–162, 1983. View at: Publisher Site  Google Scholar
 M. A. Golberg and C. S. Chen, “The theory of radial basis functions applied to the bem for inhomogeneous partial differential equations,” Boundary Elements Communications, vol. 5, pp. 57–61, 1994. View at: Google Scholar
 S. R. Karur and P. A. Ramachandran, “Radial basis function approximation in the dual reciprocity method,” Mathematical and Computer Modelling, vol. 20, no. 7, pp. 59–70, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 R. L. Hardy, “Multiquadric equations of topography and other irregular surfaces,” Journal of Geophysical Research, vol. 76, no. 8, pp. 1905–1915, 1971. View at: Publisher Site  Google Scholar
 J. Duchon, “Splines minimizing rotationinvariant seminorms in sobolev spaces,” in Constructive Theory of Functions of Several Variables, vol. 571 of Lecture Notes in Mathematics, pp. 85–100, Springer, Berlin, Germany, 1977. View at: Publisher Site  Google Scholar
 M. A. Golberg, “Numerical evaluation of particular solutions in the BEM—a review,” Boundary Elements Communications, vol. 6, no. 3, pp. 99–106, 1995. View at: Google Scholar
 M. A. Golberg, “The method of fundamental solutions for Poisson's equation,” Engineering Analysis with Boundary Elements, vol. 16, no. 3, pp. 205–213, 1995. View at: Publisher Site  Google Scholar
 C. S. Chen, “The method of fundamental solutions for nonlinear thermal explosions,” Communications in Numerical Methods in Engineering, vol. 11, no. 8, pp. 675–681, 1995. View at: Publisher Site  Google Scholar  MathSciNet
 M. A. Golberg, C. S. Chen, and S. R. Karur, “Improved multiquadric approximation for partial differential equations,” Engineering Analysis with Boundary Elements, vol. 18, no. 1, pp. 9–17, 1996. View at: Publisher Site  Google Scholar
 C.C. Tsai, “Analytical particular solutions of multiquadrics associated with polyharmonic operators,” Mathematical Problems in Engineering, vol. 2013, Article ID 613082, 11 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 M. F. Samaan and Y. F. Rashed, “BEM for transient 2D elastodynamics using multiquadric functions,” International Journal of Solids and Structures, vol. 44, no. 2526, pp. 8517–8531, 2007. View at: Publisher Site  Google Scholar
 M. F. Samaan and Y. F. Rashed, “Free vibration multiquadric boundary elements applied to plane elasticity,” Applied Mathematical Modelling, vol. 33, no. 5, pp. 2421–2432, 2009. View at: Publisher Site  Google Scholar
 C.C. Tsai and T.W. Hsu, “Multiquadric and Chebyshev approximation to threedimensional thermoelasticity with arbitrary body forces,” Engineering Analysis with Boundary Elements, vol. 37, no. 10, pp. 1259–1266, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A. H.D. Cheng, C. S. Chen, M. A. Golberg, and Y. F. Rashed, “BEM for theomoelasticity and elasticity with body force—a revisit,” Engineering Analysis with Boundary Elements, vol. 25, no. 45, pp. 377–387, 2001. View at: Publisher Site  Google Scholar
 G. Yao, The method of approximate particular solutions for solving partial differential equations [Doctor of Philosophy], Department of Mathematics, The University of Southern Mississippi, 2010.
 R. D. Mindlin, “Influence of rotary inertia and shear on flexural motions of isotropic, elastic plates,” Journal of Applied Mechanics, vol. 18, pp. 336–343, 1951. View at: Google Scholar
 P. H. Wen, M. H. Aliabadi, and A. Young, “Transformation of domain integrals to boundary integrals in BEM analysis of shear deformable plate bending problems,” Computational Mechanics, vol. 24, no. 4, pp. 304–309, 1999. View at: Publisher Site  Google Scholar
 C.C. Tsai and E. M. Wu, “Analytical particular solutions of augmented polyharmonic spline associated with mindlin plate model,” Numerical Methods for Partial Differential Equations, vol. 28, no. 6, pp. 1778–1793, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 C.C. Tsai, “The particular solutions of chebyshev polynomials for reissner plates under arbitrary loadings,” CMES: Computer Modeling in Engineering and Sciences, vol. 45, no. 3, pp. 249–271, 2009. View at: Google Scholar
 C. C. Tsai, M. E. Quadir, H. H. Hwung, and T. W. Hsu, “Particular solution of polyharmonic spline associated with reissner plate problems,” Journal of Mechanics, vol. 27, no. 4, pp. 493–501, 2011. View at: Publisher Site  Google Scholar
 E. Reissner, “The effect of transverse shear deformation on the bending of elastic plates,” Journal of Applied Mechanics, vol. 12, pp. 69–76, 1945. View at: Google Scholar  MathSciNet
 E. Reissner, “On the theory of bending of elastic plates,” Journal of Applied Mechanics, vol. 23, pp. 184–191, 1945. View at: Google Scholar
 L. Hörmander, Linear Partial Differential Operators, Springer, Berlin, Germany, 1963. View at: Publisher Site
 H. Wendland, “On the smoothness of positive definite and radial functions,” Journal of Computational and Applied Mathematics, vol. 101, no. 12, pp. 177–188, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 S. Y. Reutskiy and C. S. Chen, “Approximation of multivariate functions and evaluation of particular solutions using Chebyshev polynomial and trigonometric basis functions,” International Journal for Numerical Methods in Engineering, vol. 67, no. 13, pp. 1811–1829, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 R. Schaback, “Error estimates and condition numbers for radial basis function interpolation,” Advances in Computational Mathematics, vol. 3, no. 3, pp. 251–264, 1995. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2015 ChiaCheng Tsai. 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.