Abstract
To analyze thin and thick plates, the paper presents two rectangular finite elements with high accuracy. In these elements, the proposed formulations of the displacement field utilize the Bergan-Wang approach, which depends only on one variable: the plate lateral deflection. This approach ensures that shear-locking problem will not happen as thickness decreases. The degrees of freedom of the proposed elements are twenty-four for the first element and it is named BWRE24, while the second one has thirty-six degrees of freedom and is named BWRE36. To evidence the efficiency of the two elements, a series of numerical examples for an isotropic plate subjected to various loadings and with different boundary conditions have been analyzed. Very good results are obtained suffering no numerical difficulties in case of very thin plates.
Keywords:
Bergan-Wang approach; Thin and Thick plate; Displacement finite element formulation
1. Introduction
Plates have wide applications in different engineering fields. There are many theories for the analysis of plate bending problems. The classical plate theory (CPT) based on the assumption that the perpendicular line to the middle plane of the un-deformed plate remains straight and perpendicular to the deformed mid-surface. The only variable in CPT is the lateral deflection of the plate middle plane. Because of neglecting the effect of shear deformation, the (CPT) applies mainly to thin plates [1[1] Gou, Y., Y. Cai, and H. Zhu, A Simple High-Order Shear Deformation Triangular Plate Element with Incompatible Polynomial Approximation. Applied Sciences, 2018. 8(6): p. 975., 2[2] Noor, A.K.; Burton, W.S. Assessment of shear deformation theories for multilayered composite plates.Appl. Mech. Rev. 1989, 42, 1-13.].
For moderately thick and thick plates with thickness to length ratio over 0.05, it is necessary to take into account the effect of transverse shear deformations [3[3] Fuh-Gwo, Y. and R.E. Miller, 1988, A rectangular finite element for moderately thick flat plates. Computers & structures. 30(6): p. 1375-1387.]. The mostly used theory in engineering applications is that of Reissner - Mindlin plate theory. They modified the Kirchhoff's normality assumption by allowing uniform rotation (constant shear- Strain and stress) for the normal line about the middle surface after deformation [4[4] Lee, K., G. Lim, and C. Wang, Thick Lévy plates re-visited. International Journal of Solids and Structures, 2002. 39(1): p. 127-144.].
According to Reissner-Mindlin plate theory, There are three independent variables: the lateral deflection of the plate middle plane and the two rotations of the normal line to the plate and. These three variables are governed by a coupled system of partial differential equations of six order. Therefore, it is necessary to satisfy three boundary conditions along the edges of a plate while only two boundary conditions are required in the classical plate theory [5[5] Soh, A.-K., et al., 2001, A new twelve DOF quadrilateral element for analysis of thick and thin plates. European Journal of Mechanics-A/Solids. 20(2): p. 299-326.
[6] Reissner, E. The effect of transverse shear deformation on the bending of elastic plates. J. Appl. Mech. 1945, 12, 69-77.-7[7] Mindlin, R.D. Influence of rotatory inertia and shear on flexural motions of isotropic elastic plates. J. Appl. Mech. 1951, 18, 31-38.].
Despite the simple formulation of The Reissner-Mindlin theory, the discretization by means of finite elements based on standard low-order schemas show an intense lack of convergence when the thickness of the plate is too small compared with the other dimensions of the plate. This unwanted phenomenon, known as shear locking, can be overcome by reducing the effect of the shear energy [8[8] Da Veiga, L.B., et al., An isogeometric method for the Reissner-Mindlin plate bending problem. Computer Methods in Applied Mechanics and Engineering, 2012. 209: p. 45-53., 9[9] F. Auricchio, C. Lovadina, Analysis of kinematic linked interpolation methods for Reissner-Mindlin plate problems, Comput. Methods Appl. Mech. Engrg. 190 (2001) 2465-2482.].
P.G. Bergan and X. Wang [10[10] Bergan, P.G. and X. Wang, Quadrilateral plate bending elements with shear deformations. Computers & Structures, 1984. 19(1-2): p. 25-34.] H. Abdalla and K. Hassan [11[11] Abdalla, H. and K. Hassan, On the Bergan‐Wang approach for moderately thick plates. Communications in applied numerical methods, 1988. 4(1): p. 51-58.] have employed an approximation that have led to general expressions of rotations of the normal as functions of displacement. Therefore, the strain energy is function of only one variable Second, third, and fourth derivatives of this variable appear in the strain energy and consequently, finite elements of displacement type require polynomials of order larger than or equal to four. We notice that P.G. Bergan and X. Wang [10[10] Bergan, P.G. and X. Wang, Quadrilateral plate bending elements with shear deformations. Computers & Structures, 1984. 19(1-2): p. 25-34.] used quadrilateral finite element, which is only complete for a polynomial of degree three.
Liu et al. [12[12] Ine-Wei Liu, Tienfuan Kerh, Chien-Chang Lin, A conforming quadrilateral plate bending element with shear deformations, Computers & structures, 1995, 56(1), 93-100.] have also followed the same approximation. They present the formulation and assessment of a quadrilateral finite element for thin and moderately thick plates (NCQS). The NCQS element appears to be a simple and consistent 16-degree-of-freedom (DOF) conforming quadrilateral plate-bending element, but also complete for a third-degree polynomial only.
In this paper, new rectangular finite elements of a complete, forth degree polynomial and fifth degree polynomial based on Bergan - Wang principal approximation. The elements have been named as BWRE24 and BWRE36 respectively. We have tested the performance of these elements on different boundary condition of thin and thick plate.
For the reader's convenience, we start by briefly summarizing the theoretical formulation of the strain energy expression introduced by Bergan-Wang approximation for the case of isotropic material. Followed by the introduced displacement finite elements with the formulation of corresponding stiffness matrices. Accuracy illustrated on several numerical examples. Work is in progress to extend the formulation to anisotropic plates.
2. The governing equations of the first order shear plate theory
In Mindlin-Reissner theory, the normal on the middle un-deformed plane of the plate may rotate without remaining normal as shown in Fig. 1-a. Therefore, the displacement fields [13[13] Ozkul, T.A. and U. Ture, The transition from thin plates to moderately thick plates by using finite element analysis and the shear locking problem. Thin-Walled Structures, 2004, 42(10): p. 1405-1430.] are.
Where, and are the rotations of the transverse line about x and y-axes, respectively. Moreover, and are the shear strain of the mid-plane of the plate.
The strain deformation vector of the plate due to bending and shear can be written as.
The stress vector due to bending and shear effects depending on stress-strain relationship in the isotropic homogenous elastic plate can be written as.
Where, the flexure bending and shear rigidity of isotropic elastic plate is defined as.
In Eq. (5), is the shear correction factor, which is taken as 5/6 [6[6] Reissner, E. The effect of transverse shear deformation on the bending of elastic plates. J. Appl. Mech. 1945, 12, 69-77.], is the Poisson’s ratio, is the Young’s modulus and is the modulus of rigidity.
By integrating the planar stresses in Eq. (3) over the thickness of the plate, bending moments and twisting moment per unit length in Fig.1-b can obtained as
Moreover, the foundational equations of the shear force per unit length in Fig. 1-b are
Moreover, according to the principle of virtual work [14[14] Reddy, J.N., Energy principles and variational methods in applied mechanics, John Wiley & Sons, 2001.], we get:
(a) Un-deformed and deformed geometries of an edge of the plate according to Reissner and Mindlin theories. (b) Bending moments, shear forces and edge forces acting on the plate.
The governing differential equations of plate can written as
3. Bergan-Wang energy expression.
Bergan and Wang obtained a mathematical expression for the shear strains in terms of the transverse displacement only by substituting the shear forces and from Eq. (9.b and 9.c) into Eq. (7) the shear strain can written as [12[12] Ine-Wei Liu, Tienfuan Kerh, Chien-Chang Lin, A conforming quadrilateral plate bending element with shear deformations, Computers & structures, 1995, 56(1), 93-100.]:
A similar expression can be obtained for
Bergan and Wang assumed that, the shear deformations of the plate and can be approximated by piecewise linear polynomials. This assumption is suitable for finite element technique and leads to the following expressions [10[10] Bergan, P.G. and X. Wang, Quadrilateral plate bending elements with shear deformations. Computers & Structures, 1984. 19(1-2): p. 25-34.]:
Where,
Substituting Eq. (11.a and 11.b), into the expressions for the bending and twisting moments (Eq. (6) and Eq. (7)) leads to:
Finally, the strain energy expression per unit area, which is the sum of bending energy and shear energy according to Bergan and Wang approach [10[10] Bergan, P.G. and X. Wang, Quadrilateral plate bending elements with shear deformations. Computers & Structures, 1984. 19(1-2): p. 25-34., 11[11] Abdalla, H. and K. Hassan, On the Bergan‐Wang approach for moderately thick plates. Communications in applied numerical methods, 1988. 4(1): p. 51-58.], can be written as:
Where,
4. Formulation of finite element method.
The basic thought of the present approach in Eq. (16) and Eq. (17) can achieved by using a higher order polynomial approximation (higher than the third degree) of the displacement field. Therefore, two new rectangular finite elements of a complete fourth and fifth degree polynomial [15[15] P.G. Ciarlet, The finite element method for elliptic problems, North Holland, 1975., 16[16] O.C. Zienkiewicz, R. Taylor, J.Z. Zhu, The Finite element method: its basis and Fundamentals, 7 edition, Butterworth-Heinemann, 2013.] have used in the following plate bending problems and named as BWRE24 and BWRE36 respectively.
4.1 Derivation of element transverse displacement function
In both above plate elements, the transverse displacement function can expressed as generalized displacement modes
The number of terms for polynomial must be equal to the total number of degrees of freedom that each element has. The transverse displacement coefficients vector (where, from 1 to number of degrees of freedom of one element) can be evaluated by substituting with the local coordinates of each node of the element into the degrees of freedom mathematical formula for each plate element BWRE24 and BWRE36. This leads to a square matrix of size , consisting of numbers only [17[17] AboElsooud M.T., Vibration control of plates using periodically distributed shunted piezoelectric patches. (Doctoral dissertation) University of Maryland, College Park. 2003.]. Finally, the unknown coefficients vector can obtained from
Therefore, the transverse displacement function can written as
4.2 Derivation of element stiffness matrices and load vector
According to the strain energy expression (Eq. (16), (17)) and the generalized strains in Eq. (2), the generalized stiffness matrix of a thick isotropic plate element given by a bending part plus a shear part [18[18] Oñate, E., Structural analysis with the finite element method. Linear statics: volume 2: beams, plates and shells, Springer Science & Business Media, 2013.].
Where,
The generalized bending strain vector and shear strain vector can written according to the approach presented in Eq. (13, 14) and according to the finite element formulation of the transverse displacement function in Eq. (20) as
Finally, the finite element formulation of the element stiffness matrix is expressed by:
Similarly, according to the principle of virtual work (Eq. (8)), and the finite element formulation of the transverse displacement function (Eq. (20)), the finite element formulation of the element load vector [19[19] Ugural, A.C., Stresses in plates and shells, 2nd revised edition,McGraw-London, 1998.] can written as
After assembling the elements stiffness matrix and the element force vector to obtain the formation of the overall plate stiffness matrix and the resultant systems force vector, the final equilibrium equation can written as
5. The selection of the displacement field.
Consider a rectangular plate of dimensions and thickness (h) have been meshed by two models of rectangular finite element which the origin of its local coordinate system is located at the lower-left node of the element.
5.1 Finite element formulation of BWRE24
The first finite element model, BWRE24 has four nodes and the degrees of freedom at each node are the lateral deflection, two rotations and, the two degrees of freedom and which are proportional to the internal bending moments per unit length and , and finally which is proportional to the internal twisting moment per unit length. So, the generalized displacement field at each node of the element can expressed as:
Because of the six degrees of freedom at each node, only twenty-four terms of Pascal's triangle are selected [20[20] Popplewell, N. and D. McDonald, Conforming rectangular and triangular plate-bending elements. Journal of Sound and Vibration, 1971. 19(3): p. 333-347.], Such that as tends to zero BWRE24 will be as 24-dof displacement type finite element model [20[20] Popplewell, N. and D. McDonald, Conforming rectangular and triangular plate-bending elements. Journal of Sound and Vibration, 1971. 19(3): p. 333-347.]. Therefore, the deflection function as a complete polynomial of degree four can expressed as:
5.2 Finite element formulation of BWRE36
For the second finite element model, BWRE36 has nine nodes as shown in Fig. 2 and the degrees of freedom at each node are the lateral deflection, two rotations and, and which is proportional to the internal twisting moment per unit length. So, the generalized displacement field at each node of the element can be expressed as:
Because of the four degrees of freedom at each node, only thirty-six terms of Pascal's triangle have selected, Such that as tends to zero, BWRE36 will be as the finite element model in [16[16] O.C. Zienkiewicz, R. Taylor, J.Z. Zhu, The Finite element method: its basis and Fundamentals, 7 edition, Butterworth-Heinemann, 2013.]. Therefore, the deflection function as a complete polynomial of degree five can expressed as:
6. Numerical tests
6.1 Bending of a square plate under uniformly distributed load with various boundary conditions
In order to verify the validity of the two finite element models BWRE24 and BWRE36, four cases of a thin and moderately thick square plate of side length ,having 0.3, subjected to uniform load of intensity have been studied.
6.1.1 The results shown in Table 1 for case of two opposite edges simply supported and the other two edges clamped SCSC. With an illustration of the geometric shape of deflection and slopes for the plate at a certain case when h/a=0.1 as shown in Fig.3. And, comparing the plate’s central deflection and slopes results from BWRE24 and BWRE36 with these results from the Reissner-Mindlin plate theory as shown in Fig.4 and Fig.5
The non-dimensional central deflections and bending moments for a uniformly distributed load with two opposite edges simply supported and the others are clamped SCSC.
Comparison between R-M, BWRE24 and BWRE36 in Deflection w of SCSC plate at midpoint with different h/a ratio.
Comparison between R-M, BWRE24 and BWRE36 in rotations and of SCSC plate at centroidal axes.
6.1.2 The results shown in Table 2 for case of two opposite edges simply supported and the other two edges free SFSF. With an illustration of the geometric shape of deflection and slopes for the plate at a certain case when h/a=0.1 as shown in Fig.6. And, comparing the plate’s central deflection and slopes results from BWRE24 and BWRE36 with these results from the Reissner-Mindlin plate theory as shown in Fig.7 and Fig.8
The non-dimensional central deflections and bending moments for a uniformly distributed load with two opposite edges simply supported and the others are free SFSF.
Comparison between R-M, BWRE24 and BWRE36 in Deflection w of SFSF plate at midpoint with different h/a ratio.
Comparison between R-M, BWRE24 and BWRE36 in rotations and of SFSF plate at centroidal axes.
6.1.3 The results shown in Table 3 for case of three edges simply supported and the other edge clamped SCSS. With an illustration of the geometric shape of deflection and slopes for the plate at a certain case when h/a=0.1 as shown in Fig.9. And, comparing the plate’s central deflection and slopes results from BWRE24 and BWRE36 with these results from the Reissner-Mindlin plate theory as shown in Fig.10 and Fig.11
The non-dimensional central deflections and bending moments for a uniformly distributed load with two opposite edges simply supported and the others are clamped SCSS.
Comparison between R-M, BWRE24 and BWRE36 in Deflection w of SCSS plate at midpoint with different h/a ratio.
Comparison between R-M, BWRE24 and BWRE36 in rotations and of SCSS plate at centroidal axes.
6.1.4 The results shown in Table 4 for case of three edges simply supported and the other edge free SFSS. With an illustration of the geometric shape of deflection and slopes for the plate at a certain case when h/a=0.1 as shown in Fig.12. And, comparing the plate’s central deflection and slopes results from BWRE24 and BWRE36 with these results from the Reissner-Mindlin plate theory as shown in Fig.13 and Fig.14
The non-dimensional central deflections and bending moments for a uniformly distributed load with two opposite edges simply supported and the others are clamped SFSS.
Comparison between R-M, BWRE24 and BWRE36 in Deflection w of SFSS plate at midpoint with different h/a ratio.
Comparison between R-M, BWRE24 and BWRE36 in rotations and of SFSS plate about x and y axes that passes through point ().
6.2 C. Chinosi and C. Lovadina Problem
Both finite element models BWRE24 and BWRE36 have performed on a model problem for which the exact solution explicitly known [24[24] C. Chinosi, C. Lovadina, Numerical analysis of some mixed finite element methods for Reissner-Mindlin plates, Comput. Mech. 16 (1995), 36-44.]. This test consists of a unitary square block [0, 1]2 with clamped boundary conditions on all four sides and a distributed load given by the following function:
The exact solution for the displacement, rotations and bending moments have given in [24[24] C. Chinosi, C. Lovadina, Numerical analysis of some mixed finite element methods for Reissner-Mindlin plates, Comput. Mech. 16 (1995), 36-44., 25[25] J. Kiendla, F. Auricchio, L. Beir˜ao da Veigac, C. Lovadinad, A. Realia, Isogeometric collocation methods for the Reissner-Mindlin plate problem, Computer Methods in Applied Mechanics and Engineering, 284 (2015), 489-507.] as:
In the following table 5, 6 this plate problem is solved according to Reissner-Mindlin theory and due to Bergan-Wang approach by the two elements BWRE24 and BWRE36, With an illustration of the geometric shape of deflection and slopes for the plate at a certain case when h/a=0.1 as shown in Fig. 15. For each study a thin plate with h = 0.001 and a thick plate with h = 0.1 are considered, and comparing the plate’s central deflection and slopes results from BWRE24 and BWRE36 with these results from the Reissner-Mindlin plate theory as shown in Fig. 16 and Fig. 17 In order to test the performance of each element and investigate the presence of shear locking.
The central deflections, rotation and bending moments for a functional distributed load with clamped boundary conditions on all four sides (h/a = 0.001), 10.92 and 0.3
The central deflections, rotation and bending moments for a functional distributed load with clamped boundary conditions on all four sides (h/a = 0.1), 10.92 and 0.3
Comparison between exact solution, R-M, BWRE24 and BWRE36 in Deflection w of CCCC plate at midpoint with h/a =0.001 and h/a = 0.1.
Comparison between R-M, BWRE24 and BWRE36 in rotations and of CCCC plate at centroidal axes.
7. Conclusion
In this paper, the isotropic thick plate-bending problem has been analyzed by using two new elements. The formulations of these elements is based on a general expression for the strain energy of the plate as a function of plate deflection only. Therefore, a complete polynomial of fourth degree for the first element (BWRE24) and a complete polynomial of fifth degree for the second element (BWRE36) have taken in the constructing shape functions for both elements. It emphasized that the shear-locking problem in the thin plates has been avoided.
According to this hypothesis, the internal displacement fields for both thin and thick plate have calculated under different loads by using the same three boundary conditions of Mindlin Reissner plate theory. Numerical results show that the two elements offer good accuracy and convergence rate for both thin and thick plate bending problem. And of course, using complete polynomial of degree five (BWRE36 element) has a better performance than using complete polynomial of only degree four (BWRE24).
Applying this approach to orthotropic material and composite material is a relatively easy task. For general domain, triangular elements such as Bell and Argyris elements [15[15] P.G. Ciarlet, The finite element method for elliptic problems, North Holland, 1975.] are good candidates, one needs to modify them using the same degrees of freedom as BWRE24.
Acknowledgment
The authors thank professors: H. Asks of University of Sheffield, F. F. Mahmoud of Zagazig University, H. Abd Allah of Military Technical College and M. El Shazly of British University in Egypt for useful discussions and comments.
References
-
[1]Gou, Y., Y. Cai, and H. Zhu, A Simple High-Order Shear Deformation Triangular Plate Element with Incompatible Polynomial Approximation. Applied Sciences, 2018. 8(6): p. 975.
-
[2]Noor, A.K.; Burton, W.S. Assessment of shear deformation theories for multilayered composite plates.Appl. Mech. Rev. 1989, 42, 1-13.
-
[3]Fuh-Gwo, Y. and R.E. Miller, 1988, A rectangular finite element for moderately thick flat plates. Computers & structures. 30(6): p. 1375-1387.
-
[4]Lee, K., G. Lim, and C. Wang, Thick Lévy plates re-visited. International Journal of Solids and Structures, 2002. 39(1): p. 127-144.
-
[5]Soh, A.-K., et al., 2001, A new twelve DOF quadrilateral element for analysis of thick and thin plates. European Journal of Mechanics-A/Solids. 20(2): p. 299-326.
-
[6]Reissner, E. The effect of transverse shear deformation on the bending of elastic plates. J. Appl. Mech. 1945, 12, 69-77.
-
[7]Mindlin, R.D. Influence of rotatory inertia and shear on flexural motions of isotropic elastic plates. J. Appl. Mech. 1951, 18, 31-38.
-
[8]Da Veiga, L.B., et al., An isogeometric method for the Reissner-Mindlin plate bending problem. Computer Methods in Applied Mechanics and Engineering, 2012. 209: p. 45-53.
-
[9]F. Auricchio, C. Lovadina, Analysis of kinematic linked interpolation methods for Reissner-Mindlin plate problems, Comput. Methods Appl. Mech. Engrg. 190 (2001) 2465-2482.
-
[10]Bergan, P.G. and X. Wang, Quadrilateral plate bending elements with shear deformations. Computers & Structures, 1984. 19(1-2): p. 25-34.
-
[11]Abdalla, H. and K. Hassan, On the Bergan‐Wang approach for moderately thick plates. Communications in applied numerical methods, 1988. 4(1): p. 51-58.
-
[12]Ine-Wei Liu, Tienfuan Kerh, Chien-Chang Lin, A conforming quadrilateral plate bending element with shear deformations, Computers & structures, 1995, 56(1), 93-100.
-
[13]Ozkul, T.A. and U. Ture, The transition from thin plates to moderately thick plates by using finite element analysis and the shear locking problem. Thin-Walled Structures, 2004, 42(10): p. 1405-1430.
-
[14]Reddy, J.N., Energy principles and variational methods in applied mechanics, John Wiley & Sons, 2001.
-
[15]P.G. Ciarlet, The finite element method for elliptic problems, North Holland, 1975.
-
[16]O.C. Zienkiewicz, R. Taylor, J.Z. Zhu, The Finite element method: its basis and Fundamentals, 7 edition, Butterworth-Heinemann, 2013.
-
[17]AboElsooud M.T., Vibration control of plates using periodically distributed shunted piezoelectric patches. (Doctoral dissertation) University of Maryland, College Park. 2003.
-
[18]Oñate, E., Structural analysis with the finite element method. Linear statics: volume 2: beams, plates and shells, Springer Science & Business Media, 2013.
-
[19]Ugural, A.C., Stresses in plates and shells, 2nd revised edition,McGraw-London, 1998.
-
[20]Popplewell, N. and D. McDonald, Conforming rectangular and triangular plate-bending elements. Journal of Sound and Vibration, 1971. 19(3): p. 333-347.
-
[21]Wang, C.M., et al., Relationships between bending solutions of Reissner and Mindlin plate theories. Engineering structures, 2001. 23(7): p. 838-849.
-
[22]Timoshenko, S.P. and S. Woinowsky-Krieger, Theory of plates and shells. McGraw-hill, 1959.
-
[23]Siva Srinivas Kolukula, isoparametric Q4 elements have used to dicretize the plate according to Mindlin-Reissner plate theory, Indira Gandhi Center for Atomic Research, http://www.mathworks.de/matlabcentral/fileexchange/32029-plate-bending, 2011.
» http://www.mathworks.de/matlabcentral/fileexchange/32029-plate-bending -
[24]C. Chinosi, C. Lovadina, Numerical analysis of some mixed finite element methods for Reissner-Mindlin plates, Comput. Mech. 16 (1995), 36-44.
-
[25]J. Kiendla, F. Auricchio, L. Beir˜ao da Veigac, C. Lovadinad, A. Realia, Isogeometric collocation methods for the Reissner-Mindlin plate problem, Computer Methods in Applied Mechanics and Engineering, 284 (2015), 489-507.
Edited by
Editor:
Publication Dates
-
Publication in this collection
18 Sept 2020 -
Date of issue
2020
History
-
Received
26 June 2020 -
Reviewed
13 July 2020 -
Accepted
04 Aug 2020 -
Published
08 Aug 2020