RBF-Based Meshless Method for Large Deflection of Elastic Thin Rectangular Plates with Boundary Conditions Involving Free Edges

Hussein Al-Tholaia, Mohammed M.; Al-Gahtani, Husain Jubran

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

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusions Acknowledgments References Copyright Related Articles

Research Article | Open Access

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

RBF-Based Meshless Method for Large Deflection of Elastic Thin Rectangular Plates with Boundary Conditions Involving Free Edges

Mohammed M. Hussein Al-Tholaia¹and Husain Jubran Al-Gahtani¹

Academic Editor: Evangelos J. Sapountzakis

Received17 Feb 2016

Revised22 Apr 2016

Accepted12 May 2016

Published09 Jun 2016

Abstract

An RBF-based meshless method is presented for the analysis of thin plates undergoing large deflection. The method is based on collocation with the multiquadric radial basis function (MQ-RBF). In the proposed method, the resulting coupled nonlinear equations are solved using an incremental-iterative procedure. The accuracy and efficiency of the method are verified through several numerical examples. The inclusion of the free edge boundary condition proves that this method is accurate and efficient in handling such complex boundary value problems.

1. Introduction

The nonlinear deflection problems arise when using the light weight structures. In several practical light weight structures applications such as fluid tank, aircraft, and hydrospace, flexible plates are commonly used and undergo large deflection. Consequently, the analysis of the plates becomes more difficult. The governing equations become coupled and highly nonlinear [1]. The available analytical methods are based on simplified assumptions and are limited to simple loading and boundary conditions [1–5]. For such complicated problems, numerical methods offer convenient and reliable solutions. The ideal numerical method for the solution of nonlinear partial differential equations (PDEs) such as the one considered here should be high-order accurate, flexible with respect to the geometry, computationally efficient, and easy to implement. The conventional numerical methods that are commonly used usually fulfill one or two of the above criteria, but not all. Finite difference methods (FDM), finite element methods (FEM), and boundary element methods (BEM) have been the dominating methods for the numerical solution of PDEs [6–13]. Referring to the most dominant approach, that is, FEM, it is highly flexible, but it is hard to achieve high-order accuracy and both coding and mesh generation become increasingly difficult as the problem dimension increases. The use of a mesh implies that specific procedures have to be devised just to define the mesh. Also, and to keep the order of the local approximation within reasonable limits, the element size has to be reduced, whenever better approximations are pursued. The extraordinary amount of work, which has been put into FEM research since its early years, has, one way or another, circumvented these and other problems associated with the existence of a mesh and made FEM the dominant approach for most problems in computational mechanics. Accordingly, many sophisticated powerful codes (e.g., ANSYS, ABACUS, and COMSOL) have been established and have proven to be reliable in solving almost any computational mechanics problem. FDM can be made high-order accurate in resolving PDEs but require a structured grid (or a collection of structured grids), which makes it difficult to model features of irregular domain. Furthermore, solutions of PDEs using FDM can be derived from the assumptions of the local interpolation schemes and require a mesh to support the localized approximations; however, the construction of a mesh in two or more dimensions is a nontrivial problem. In recent years, BEM has become a powerful alternative to FEM and FDM, especially for problems involving high gradients and stress concentrations. It has been successfully applied to solve the problems of large deflection of thin elastic plates. However, this was possible by devising some techniques to overcome the inherent deficiency of BEM as a self-standing numerical method in handling nonlinearities.

Nevertheless, the possibility of obtaining numerical solutions for PDFs without resorting to element frame has been the goal of many researchers throughout the computational mechanics community for the past three decades or so. One of the earliest meshless methods is the smoothed particle hydrodynamics (SPH) which was developed in 1977 by Lucy et al. [13, 14]. Since then, many versions and types of the meshless methods have been proposed. Examples are the radial basis function-based collocation method (RBF) by Kansa [15, 16], the diffuse element method (DEM) by Nayroles et al. [17], the element-free Galerkin (EFG) method by Belytschko el al. [18], the reproduced kernel particle method (RKPM) by Liu et al. [19], the finite point method by Onate et al. [20], the meshless local Petrov-Galerkin method (MLPG) by Atluri and Zhu [21], and the mesh-free weak-strong form (MWS) by Liu and Gu [22].

RBF method, as one of the efficient meshless methods, has attracted attention in recent years especially in the area of computational mechanics. In general, RBF method expands the solution of a problem in terms of RBFs and chooses expansion coefficients such that the governing equations and boundary conditions are satisfied at some selected domain and boundary points. This method does not require mesh generation which makes it advantageous for nonlinear problems that require frequent remeshing such as the one considered in this study. The roots of RBF go back to the early 1970s when it was first used for fitting scattered data [23]. In the early 1980s, it was coupled with BEM in a technique called dual reciprocity-boundary element method where the RBF was employed to transform the domain integrals into boundary integrals [24]. Thereafter, many researchers have used RBF in conjunction with BEM to solve various problems in computational mechanics. The method, however, has not been applied directly to partial differential equations until 1990 by Kansa [15, 16]. Since then, many researchers have suggested several variations to the original method [25–31].

Most of previous applications of meshless methods to large deflection of plates are limited to simple and clamped boundary conditions. The inclusion of free edge boundary condition poses a challenge for both analytical and numerical methods. The difficulty of the free edge boundary condition is due to the enforcement of zero shear on the free boundary of the plate which, for large deflection, involves lengthy and highly nonlinear differential operators as explained in Section 2. In this paper, a multiquadric radial basis function- (MQ-RBF-) based meshless model is developed for the solution of large deflection of thin plates with different boundary conditions including free edges. The proposed method has the advantages of simplicity and ease of code implementation. The accuracy of the model is validated through several numerical examples.

2. Governing Equations

The governing equilibrium equations of a plate undergoing large deflection can be expressed in terms of the transverse displacement and a stress function or in terms of and the in-plane displacements and . The summary of both formulations and their corresponding boundary conditions are given below.

2.1. Formulation

In this formulation, the governing equations for large elastic deflection of thin plates are represented by two coupled nonlinear differential equations, known as Von Kárámn equations, which are expressed in terms of the displacement and a stress function [1]:where is the distributed load, is the plate thickness, and is the flexural rigidity of the plate having elastic constants and , and is a nonlinear differential operator given byThe stress function is related to the membrane forces , , and by the following differential operators:The details for deriving (1) and (2) are given in the classical book by Timoshenko and Woinowsky-Kreiger [1]. The first equation describes the transverse or bending action whereas the second equation represents the in-plane action of the plate. In general, there are two types of in-plane boundary conditions, namely, movable and immovable. The formulation can be utilized directly in the former type where the movable boundary condition is represented by . However, the immovable boundary condition () cannot be accurately expressed in terms of the stress function . To overcome this difficulty, (2) can be replaced by two partial differential equations in terms of the in-plane displacements, and , as described in the following section.

2.2. Formulation

Excluding the body forces, the equilibrium equations along and are given bythe membrane strains are related to the membrane through the following equations:The membrane strains are also related to the displacements by the following equations:Using (7) into (6) and then in (5), we getTo complete the formulation, (1) is rewritten after replacing the stress function by the proper terms involving the displacements , , and . Equation (1) becomesEquations (8) and (9) can be expressed in compact forms as follows:whereIn the above equations, subscript variables denote partial differentiation with respect to the variables.

2.3. Boundary Conditions

The general boundary conditions for large deflection of plates can be classified into two types.

(1) The first is transverse boundary conditions which are encountered in both small and large deflection formulations. For this type, we will assume that, at each boundary point, there are two prescribed boundary conditions:(a) (b) where is the normal bending moment which is given byThe shear force is composed of a linear and a nonlinear part; that is, , where is given byThe nonlinear part of the shear, , is given byThe above nonlinear formula for the nonlinear part of the shear force is suitable for the formulation. The nonlinear part of the shear corresponding to the formulation can be obtained by replacing the stress function with the proper derivatives of , , and . The result iswhere and are the and components of the unit vector normal to the boundary, respectively.

(2) The second is in-plane boundary conditions which have to be addressed in the case of large deflection formulation. This type of boundary condition is further classified into two types of edge conditions:(a)Movable edge (free to move in the in-plane directions) which is represented by the following equation: (b)Immovable edge which is represented byIn summary, the plate problem with movable edges can be handled by the formulation represented by the governing equations (1) and (2), the lateral boundary conditions (14) and (15), and the in-plane boundary conditions (20), while the plate problem with immovable edges can be handled by the formulation represented by the governing equations (8) and (9), the lateral boundary conditions (14) and (15), and the in-plane boundary conditions (21).

3. RBF Formulation

Consider the 2D computational domain (Figure 1) that represents the plate geometry. For collocation, we use node points distributed both along the boundary and over the interior . Let , so that the total number of points called poles is .

The transverse deflections , the stress function , and the in-plane deflections and are interpolated linearly by suitable radial basis functions, respectively:where is the magnitude of the concentrated load acting at the th point with a multiquadratic radial basis function (MQ-RBF), is the shape factor, and , , , and are unknown coefficients to be determined by applying the governing equations at the domain points and satisfying the boundary conditions at the boundary points.

For the formulation, the collocation of (1) and (2) yields the following coupled nonlinear algebraic equations:Similarly, the collocation of (10), (11), and (12), for the formulation, yields the following coupled nonlinear algebraic equations:where the collocations of (10) and (11) have been combined to get (26). For the solution of (25), (26), and (27), the loads are applied incrementally and then solved by iterations until the convergence is satisfied for each incremental load. The nonlinear terms and − in (25) and , , and , in (26) and (27), are put equal to zero in the first iteration of the first incremental load and then (25), (26), and (27) are solved accordingly to get , , , and functions after the nonlinear terms are evaluated to be used in the next iteration. This procedure is iterated until convergence is satisfied; otherwise, the load increment is reduced. The iterative procedure is repeated by adding load increments and performing iterations until the total load is applied.

4. Numerical Examples

In order to examine the effectiveness of the proposed RBF method, the following three examples are considered. In all examples, the loads are assumed to be either uniformly distributed = or concentrated = , Poisson ratio is assumed 0.3, and the analysis was performed for several combinations of boundary conditions and shapes. For generality of the solutions, all results are made dimensionless, so that the coordinates, the load, the deflection, and the stress are represented by , , , ; , , respectively. The shape factor, , of the RBF is changed according to the boundary conditions and the type of load. Its optimum value ranged between 0.3 and 0.8.

Example 1. Consider a square plate with three simply supported edges and the fourth edge free (denoted by SSSF). All edges are allowed to move in the in-plane directions and hence, the formulation is used. The plate is analyzed for two different load types: a uniformly distributed load ranging from 1.5 to 12 and a central concentrated load ranging from 1 to 8. The plate is discretized using a uniform node distribution consisting of 36 boundary nodes and 81 domain nodes as shown in Figure 2. The results of the deflections as obtained by RBF and FEM at the center and at the mid-point of the free edge are given in Tables 1 and 2 for the uniform and central concentrated load, respectively. The maximum differences in the deflection at the center of the plate are 1.25% and 3.95% for the uniform and concentrated load, respectively, whereas the differences at the free edge of the plate are 1.26% and 4.73% for the uniform and concentrated loads, respectively. The maximum differences in the bending stresses at the center of the plate are 2.93% and 4.49% for the uniform load and concentrated load, respectively. For the stresses at the mid-point of the free edge, the maximum differences are 3.75% and 4.95% for the uniform and concentrated load, respectively.

Example 2. Consider the plate of Example 1 with its simply supported edges replaced by clamped edges and let us denote it by CCCF. The uniform load is incrementally increased from 5 to 40 while the central concentrated load is incrementally increased from 2 to 16. The plate is discretized using a uniform node distribution similar to the one employed in Example 1. The results of the analysis of this example are given in Tables 3 and 4 for the uniform load and Tables 5 and 6 for the central concentrated load. The results show a good agreement between RBF and FEM solutions for the uniform load with a maximum difference of 1.61% in the deflection results and a maximum difference of 4.99% in the stress results. For the concentrated load, the results show a reasonable agreement between the two solutions with a maximum difference of 4.37% in the deflection results and a maximum difference of 4.88% in the stress results.

Example 3. Consider the plate of Example 2 but with immovable edges (CCCF-IM). The solution is obtained using formulation. The plate is analyzed under two different load types: uniformly distributed load ranging from 15 to 75 and a central and concentrated load ranging from 4 to 20. The results (Tables 7–10), once again, show a good agreement between RBF and FEM solutions. For the uniform load, the maximum difference is 4.66% in the deflection results and 4.92% in the stress results while, for the concentrated load, the maximum difference is 4.25% in the deflection and 4.69% in the stress results.

5. Conclusions

An RBF-based meshless method has been presented for the solution of large deflection of thin plates with different boundary conditions including free edges. The method is based on collocation with MQ-RBF. The load types considered in the study are uniform and concentrated. The boundary conditions covered are clamped, simply supported, free, movable, and immovable conditions. The FEM-based software COMSOL has been used to verify the RBF results. The proposed RBF meshless method has the advantage of being simple, easy to implement in a computer code, and reasonably accurate and therefore can qualify as an alternative numerical technique for the solution of highly nonlinear boundary value problems.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors would like to express their appreciation to King Fahd University of Petroleum & Minerals for supporting this study.

References

S. P. Timoshenko and S. Woinowsky-Kreiger, Theory of Plates and Shells, McGraw-Hill, New York, NY, USA, 1959.
C. Augural, Stresses in Plates and Shells, McGraw-Hill, New York, NY, USA, 1999.
G. H. Little, “Efficient large deflection analysis of rectangular plates with general transverse form of displacement,” Computers and Structures, vol. 71, no. 3, pp. 333–352, 1999.
View at: Publisher Site | Google Scholar
G. H. Little, “Large deflections of rectangular plates with transverse edges remaining straight,” Computers and Structures, vol. 71, no. 3, pp. 353–357, 1999.
View at: Publisher Site | Google Scholar
L. S. Ramachandra and D. Roy, “A novel technique in the solution of axisymmetric large deflection analysis of a circular plate,” Journal of Applied Mechanics, vol. 68, no. 5, pp. 814–816, 2001.
View at: Publisher Site | Google Scholar
G. Papakaliatakis and T. E. Simos, “A finite difference method for the numerical solution of fourth-order differential equations with engineering applications,” Computers & Structures, vol. 65, no. 4, pp. 491–495, 1997.
View at: Publisher Site | Google Scholar | MathSciNet
M. Duan and M. Mahendran, “Large deflection analyses of skew plates using hybrid/mixed finite element method,” Computers and Structures, vol. 81, no. 13, pp. 1415–1424, 2003.
View at: Publisher Site | Google Scholar
R. Szilard, Theories and Applications of Plate Analysis, John Wiley & Sons, New Jersey, NJ, USA, 2004.
View at: Publisher Site
J. T. Katsikadelis and M. S. Nerantzaki, “Non-linear analysis of plates by the analog equation method,” Computational Mechanics, vol. 14, no. 2, pp. 154–164, 1994.
View at: Publisher Site | Google Scholar
A. El-Zafrany, S. Fadhil, and M. Debbih, “An efficient approach for boundary element bending analysis of thin and thick plates,” Computers & Structures, vol. 56, no. 4, pp. 565–576, 1995.
View at: Publisher Site | Google Scholar
A. C. J. Luo and J. Mote, “Asymmetric responses of rotating, thin disks experiencing large deflections,” Computers & Mathematics with Applications, vol. 45, no. 1–3, pp. 217–228, 2003.
View at: Publisher Site | Google Scholar | MathSciNet
E. Byklum and J. Amdahl, “A simplified method for elastic large deflection analysis of plates and stiffened panels due to local buckling,” Thin-Walled Structures, vol. 40, no. 11, pp. 925–953, 2002.
View at: Publisher Site | 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: Publisher Site | Google Scholar
E. J. Kansa, “Multiquadrics—a scattered data approximation scheme with applications to computational fluid-dynamics. I. Surface approximations and partial derivative estimates,” Computers & Mathematics with Applications, vol. 19, no. 8-9, pp. 127–145, 1990.
View at: Publisher Site | Google Scholar | MathSciNet
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 & Mathematics with Applications, vol. 19, no. 8-9, pp. 147–161, 1990.
View at: Publisher Site | Google Scholar | MathSciNet
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 | MathSciNet
E. Onate, S. Idelsohn, O. C. Zienkiewicz, and R. L. Taylor, “A finite point method in computational mechanics. Applications to convective transport and fluid flow,” International Journal for Numerical Methods in Engineering, vol. 39, no. 22, pp. 3839–3866, 1996.
View at: Publisher Site | Google Scholar | MathSciNet
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
G. R. Liu and Y. T. Gu, “A meshfree method: meshfree weak-strong (MWS) form method, for 2-D solids,” Computational Mechanics, vol. 33, no. 1, pp. 2–14, 2003.
View at: Publisher Site | Google Scholar
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
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 | Zentralblatt MATH
C. M. C. Roque, A. J. M. Ferreira, A. M. A. Neves, C. M. M. Soares, J. N. Reddy, and R. M. N. Jorge, “Transient analysis of composite plates by radial basis functions in a pseudospectral framework,” Computers & Structures, vol. 89, no. 1-2, pp. 161–169, 2011.
View at: Publisher Site | Google Scholar
C. M. C. Roque, D. Cunha, C. Shu, and A. J. M. Ferreira, “A local radial basis functions—finite differences technique for the analysis of composite plates,” Engineering Analysis with Boundary Elements, vol. 35, no. 3, pp. 363–374, 2011.
View at: Publisher Site | Google Scholar | MathSciNet
S. Xiang, H. Shi, K.-M. Wang, Y.-T. Ai, and Y.-D. Sha, “Thin plate spline radial basis functions for vibration analysis of clamped laminated composite plates,” European Journal of Mechanics—A/Solids, vol. 29, no. 5, pp. 844–850, 2010.
View at: Publisher Site | Google Scholar
A. J. M. Ferreira, C. M. C. Roque, A. M. A. Neves, R. M. N. Jorge, C. M. M. Soares, and J. N. Reddy, “Buckling analysis of isotropic and laminated plates by radial basis functions according to a higher-order shear deformation theory,” Thin-Walled Structures, vol. 49, no. 7, pp. 804–811, 2011.
View at: Publisher Site | Google Scholar
H. J. Al-Gahtani and M. Naffa'a, “RBF meshless method for large deflection of thin plates with immovable edges,” Engineering Analysis with Boundary Elements, vol. 33, no. 2, pp. 176–183, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
M. Naffa and H. J. Al-Gahtani, “RBF-based meshless method for large deflection of thin plates,” Engineering Analysis with Boundary Elements, vol. 31, no. 4, pp. 311–317, 2007.
View at: Publisher Site | Google Scholar
M. M. Hussein Al-Tholaia and H. Jubran Al-Gahtani, “RBF-based meshless method for large deflection of elastic thin plates on nonlinear foundations,” Engineering Analysis with Boundary Elements, vol. 51, pp. 146–155, 2015.
View at: Publisher Site | Google Scholar | MathSciNet

Copyright

Copyright © 2016 Mohammed M. Hussein Al-Tholaia and Husain Jubran Al-Gahtani. 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

1261

Downloads

876

Citations