Research Article | Open Access

Volume 2015 |Article ID 246159 | https://doi.org/10.1155/2015/246159

Chia-Cheng 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

Accepted29 Dec 2014
Published01 Mar 2015

#### Abstract

Analytical particular solutions of the polyharmonic multiquadrics are derived for both the Reissner and Mindlin thick-plate models in a unified formulation. In the derivation, the three coupled second-order 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

Boundary-type 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) , 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  and Karur and Ramachandran  applied the theory of radial basis functions (RBFs) to the DRM. Among these RBFs, Hardy’s multiquadrics (MQ)  and Duchon’s augmented polyharmonic spline (APS)  are the most popular ones. For example, Golberg [12, 13], Chen , and Karur and Ramachandran  demonstrated the superiority of the APS over the ad hoc radial function. Then, Golberg et al.  further improved the accuracy of the approximated particular solution by utilizing the MQ. Recently, Tsai  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 thick-plate 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.  derived the analytical particular solutions for the Laplace equation. Samaan and Rashed [17, 18] and Tsai and Hsu  found the analytical particular solutions, respectively, for the two- and three-dimensional elasticity problems, which can be converted to a biharmonic equation using the Galerkin-Papkovich vector . Basically, the applicability of MQ is limited to the harmonic and biharmonic operators  before Tsai  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 thick-plate model , Wen et al.  and Tsai and Wu  found the analytical particular solutions for the ad hoc RBF and the APS, respectively. In addition, the analytical particular solutions for the Chebyshev polynomials  and the APS  were derived for the Reissner thick-plate model [27, 28]. Basically, a product operator of biharmonic and Helmholtz operators resulted when applying the Hörmander operator decomposition technique  to the coupled PDEs of the Reissner and Mindlin thick-plate models as demonstrated in the previous studies . This seems to infer that the accuracy improvement of the PMQ over the APS is not possible for the thick-plate models.

In this study, a method is introduced to eliminate the Helmholtz operator, which makes the application of the PMQ for the thick-plate 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 MFS-DRM. 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  thick-plate 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 second-order PDEs for the three unknown functions , and, therefore, three boundary conditions are required as where are given boundary data.

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 . 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 .

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 even-coefficient 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 floating-point 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 e-mails.

#### 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 , which begins with the definition of the adjoint operator as where is the determinant of and is the three-by-three 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. MFS-DRM Formulation

Now, we are in a position to review the application of the MFS-DRM procedure [24, 26] for solving the well-posed thick-plate 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 MFS-DRM 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  and the Reissner  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 , Chebyshev polynomials , 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 two-stage method  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 high-order PMQ over the traditional MQ can also be observed. Overall, the best accuracy of the high-order PMQ is close to that of the high-order APS while the former RBF does not need augmented polynomials. In practice, it can be cumbersome to implement the particular solutions of the high-order augmented polynomials and thus the superiority of the PMQ over the APS becomes more relevant.

 = 6 × 6 = 8 × 8 = 10 × 10 = 12 × 12 Chebyshev 2.10E − 06 8.62E − 09 2.18E − 09 3.81E − 10 MQ 6.10E − 06 4.30E − 07 4.30E − 07 2.74E − 08 PMQ ( = 2) 6.04E − 06 2.17E − 07 2.05E − 08 5.94E − 09 PMQ ( = 4) 4.51E − 06 1.20E − 07 5.12E − 09 1.66E − 09 APS (order = 1) 2.22E − 04 9.38E − 05 4.96E − 05 2.85E − 05 APS (order = 3) 3.71E − 05 9.08E − 06 3.02E − 06 1.49E − 06 APS (order = 5) 4.88E − 06 7.35E − 07 1.90E − 07 5.90E − 08 APS (order = 7) Div 8.14E − 08 1.27E − 08 3.20E − 09 APS (order = 9) Div Div 2.44E − 09 4.20E − 10 APS (order = 11) Div Div Div 3.81E − 10

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 ill-conditioned . Then, a typical error contour is given in Figure 5 for the solution obtained by the MFS-DRM based on the MQ. In the figure, the larger errors are on the boundary as expected.

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.

 = 6 × 6 = 8 × 8 = 10 × 10 = 12 × 12 Chebyshev 5.36E − 06 7.98E − 08 4.43E − 08 1.07E − 08 MQ 5.51E − 05 1.03E − 05 1.42E − 06 6.88E − 07 PMQ ( = 2) 4.97E − 05 5.07E − 06 5.51E − 07 1.15E − 07 PMQ ( = 4) 3.12E − 05 3.30E − 06 1.64E − 07 4.21E − 08 APS (order = 3) 1.85E − 03 4.77E − 04 1.64E − 04 6.88E − 05 APS (order = 5) 1.11E − 04 1.13E − 05 2.27E − 06 6.49E − 07 APS (order = 7) Div 4.20E − 06 7.06E − 07 1.63E − 07 APS (order = 9) Div Div 5.66E − 08 1.20E − 08 APS (order = 11) Div Div Div 1.08E − 08

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.

#### 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 second-order governing equations is transformed into a single sixth-order 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 103-2221-E-022-017-MY2.

1. 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
2. G. Fairweather and A. Karageorghis, “The method of fundamental solutions for elliptic boundary value problems,” Advances in Computational Mathematics, vol. 9, no. 1-2, pp. 69–95, 1998. View at: Publisher Site | Google Scholar | MathSciNet
3. 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.
4. C.-S. Liu, “A highly accurate solver for the mixed-boundary 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
5. 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
6. P. W. Partridge, C. A. Brebbia, and L. C. Wrobel, The Dual Reciprocity Boundary Element Method, Computational Mechanics, Southampton, UK, 1992. View at: MathSciNet
7. 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
8. 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
9. 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
10. 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
11. J. Duchon, “Splines minimizing rotation-invariant semi-norms 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
12. 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
13. 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
14. C. S. Chen, “The method of fundamental solutions for non-linear thermal explosions,” Communications in Numerical Methods in Engineering, vol. 11, no. 8, pp. 675–681, 1995. View at: Publisher Site | Google Scholar | MathSciNet
15. 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
16. 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
17. M. F. Samaan and Y. F. Rashed, “BEM for transient 2D elastodynamics using multiquadric functions,” International Journal of Solids and Structures, vol. 44, no. 25-26, pp. 8517–8531, 2007. View at: Publisher Site | Google Scholar
18. 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
19. C.-C. Tsai and T.-W. Hsu, “Multiquadric and Chebyshev approximation to three-dimensional 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
20. 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. 4-5, pp. 377–387, 2001. View at: Publisher Site | Google Scholar
21. 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.
22. 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
23. 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
24. 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
25. 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
26. 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
27. 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
28. E. Reissner, “On the theory of bending of elastic plates,” Journal of Applied Mechanics, vol. 23, pp. 184–191, 1945. View at: Google Scholar
29. L. Hörmander, Linear Partial Differential Operators, Springer, Berlin, Germany, 1963. View at: Publisher Site
30. H. Wendland, “On the smoothness of positive definite and radial functions,” Journal of Computational and Applied Mathematics, vol. 101, no. 1-2, pp. 177–188, 1999. View at: Publisher Site | Google Scholar | MathSciNet
31. 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
32. 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

#### More related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.