Split-Deflection Method of Classical Rectangular Plate Analysis

International Journal of Scientific and Research Publications, Volume 6, Issue 5, May 2016
ISSN 2250-3153
147
Split-Deflection Method of Classical Rectangular Plate

Analysis
1
Ibearugbulem, Owus M., 2Ibearugbulem, C. N, 3Habib Momoh and 3Asomugha, U.C

1,2,
Civil Engineering Department, Federal University of Technology, Owerri, Nigeria

3
Civil Engineering Department, Imo State University, Owerri, Nigeria
Abstract- This paper presents an split-deflection method of

classical rectangular plate analysis. In this method, the deflection
was split into x and y components of deflection. It was assumed
that the deflection of the rectangular plate is the product of these
two components. With this assumption, total potential energy
functional was derived from principles of theory of elasticity. By
direct variation of the potential energy functional, direct
governing equation was obtained. Polynomial deflection and
trigonometric deflection were respectively used for the x and y
components of the deflection for ssss plate. This deflection
components were substituted into this direct governing equation
and the coefficient of the deflection was obtained to be
0.0132629qa4/D. With this coefficient, the maximum central
deflection of the plate was obtained to be 0.0041446qa4/D. This
was compared with the maximum central deflection of Naviers
(0.00416qa4/D) and
Levis(0.00406qa4/D) plates. It was
observed that the value from the present study makes lower
bound and upper bound differences of 0.37% and 2.08% with the
values from Navier and Levi respectively. This shows that this
present method is reliable.
Index Terms- Direct variation, split-deflection, total potential
energy functional, polynomial function, trigonometric function,
deflection
The main reason for this modification is to help the analysis

who may have difficulty in obtaining orthogonal function for a
plate of a particular boundary condition. In this case, the analyst
who may have easy access to deflection equations for beams of
any boundary condition can find the proposed method quite
useful and handy.
SPLIT-DEFLECTION
The assumption here is that the general deflection, w is split
into wx and wy . That is:
Where the x and y components of the deflection are defined as:
Substituting equations (2) and (3) into equation (1) gives:
IN-PLANE DISPLACEMENTS
From the assumption that vertical shear strains are zero for
classical plate, we obtain:
I. INTRODUCTION
he common energy methods in classical plate theory analysis

include Galerkin's, Raleigh, Raleigh-Ritz, Ritz methods etc.
(Njoku et al., 2013, Ibearugbulem et al., 2014, Ibeabuchi, 2014).
These methods use orthogonal functions in form double Furrier
series (that product of two mutually perpendicular trigonometric
functions) or orthogonal polynomial. The energy equations of the
methods used earlier than now are based on the orthogonal
deflection functions (Hutchinson, 1992, Jianqiao, 1994, Ugural,
1999, Ventsel and Krauthammer, 2001, Wang et al., 2002,
Taylor and Govindjee, 2004, Szilard, 2004, Jiu et al., 2007,
Erdem et al., 2007, Ezeh et al., 2013, Ibearugbulem, 2014).
However, none of the previous energy methods have used spited
deflection functions - that is a deflections function that is
typically separated into two independent distinct functions (w =
w x * w y ). In this case wx may be polynomial and wy may be
trigonometry. In similar way stress and strain terms are also
based on individual components of w (deflection) for example
using:
STRAIN DEFLECTION RELATIONSHIP

The three in-plane strains of classical plate are obtained by
using equations (5) and (6):
STRESS STRAIN RELATIONSHIP

The in-plane constitutive equations are:
www.ijsrp.org
ISSN 2250-3153
148
STRESS DEFLECTION RELATIONSHIP

Substituting equations (7), (8) and (9) into equations (10), (11)
and (12) where appropriate gives:
Now, equation (20) can be written in non dimensional axes

R and Q.
TOTAL POTENTIAL ENERGY

The strain energy is defined as:
For pure bending analysis , the external work is given as:
Where a, b and P are the plate lengths in x and y axes and

long span- short span aspect ratio respectively.
Substituting equations (21), (22) and (23) into equation (20)
gives:
That is
Substituting equations (10) to (15) into equation (16) gives

strain energy deflection relationship as:
DIRECT VARIATION OF TOTAL POTENTIAL ENERGY

Equation (24) shall be differentiated with respect to the
deflection coefficient, A gives:
That is
Subtracting equation (17) from Equation (18) gives the total

potential energy functional as:
Substituting equations (1) and (2) into equation (19) gives:
That is
This equation (25) is the direct governing equation of

rectangular plate under pure bending from this present method.
Rearranging equation (25) and making the coefficient of
deflection the subject of the formula gives:
www.ijsrp.org
ISSN 2250-3153
149
II. RESULTS AND CONCLUSIONS
Where
NUMERICAL EXAMPLE
Analyze a classical rectangular thin isotropic plate with all
the four edge simply supported subjected to uniform distributed
load, q. Let the x and y components of deflection be:
From equations (30) (31), h 1 and h 2 are:
With these we obtain:
The result of this present study is presented on table 1. A

close critical observation of the table reveals that the maximum
recorded percentage difference is 0.276. The result from past
work was obtain from the work of Ibearugbulem et al. (2014).
They used complete orthogonal polynomial deflection equation
where as the present study used a part as polynomial function
and the other as trigonometric function. This difference in the
deflection equations between the previous study and the present
one may be responsible for the maximum recorded percentage
difference. Again, for aspect ratio of 1.0, the result from the
present study is 0.00414qa4/D. On comparing this with values
from
works
of
Naviers
(0.00416qa4/D)
and
Levis(0.00406qa4/D), it was observed that the value from the
present study makes lower bound and upper bound differences of
0.37% and 2.08% with the values from Navier and Levi
respectively.
Thus, an inference that the results obtained herein is close to
those of previous study is drawn. Hence, the present method as
present herein is reliable and is recommended for use in classical
plate analysis. For future studies, it is recommend that other
aspects of thin plate analysis like buckling and free vibration
should be carried out using the present method. It is also
recommended that this present method is extended to refined
plate theory analysis (RPT).
Table 1: Center deflection of ssss isotropic thin plate under
uniform distributed load
Substituting equations (34) to ( 29) into equation (37) gives
The maximum deflection shall occur at the center of the

rectangular plate. At the center of the circle, R = Q = 0.5.
Substituting o.5 into equations (32) and (33) gives:
Aspect
ratio, P
Deflection at center (qa4/D)

Past
(Ibearugbulem et
Present
al., 2014)
Percentage
difference
0.00414
0.000
0.00414
1.1
0.00497
0.00496
0.202
1.2
0.00577
0.00576
0.174
1.3
0.00654
0.00653
0.153
1.4
0.00727
0.00725
0.276
1.5
0.00794
0.00793
0.126
1.6
0.00857
0.00856
0.117
1.7
0.00915
0.00913
0.219
1.8
0.00968
0.00966
0.207
1.9
0.01016
0.01015
0.099
0.01061
0.01059
0.189
REFERENCES
Substituting equations (38), (39) and (40) into equation (4)
gives the equation for maximum deflection as:
[1]
[2]
Erdem, C. Imrak and Ismail Gerdemeli (2007).The problem of isotropic

rectangular plate with four clamped edges.Sadhana Vol. 32, Part 3, pp.
181186.
Ezeh, J. C., Ibearugbulem, O. M., Njoku, K. O., and Ettu, L. O. (2013).
Dynamic Analysis of Isotropic SSSS Plate Using Taylor Series Shape
Function in Galerkins Functional. International Journal of Emerging
Technology and Advanced Engineering, 3 (5): 372-375.
www.ijsrp.org
ISSN 2250-3153
[3]
[4]
[5]
[6]
[7]
[8]
Hutchinson, J. R. (1992). On the bending of rectangular plates with two

opposite edges simply supported. J. Appl. Mech. Trans. ASME 59: 679
681.
Ibeabuchi, V. T. (2014). Analysis of Elastic Buckling of Stiffened
Rectangular Isotropic Plates Using Virtual Work Principles. A Masters
Thesis submitted to the Post graduate School, Federal University of
Technology, Owerri, Nigeria.
Ibearugbulem, O. M. (2014), Using the product of two mutually
perpendicular truncated polynomial series as shape function for rectangular
plate analysis, International Journal of Emerging Technologies and
Engineering (IJETE) ISSN: 23488050, ICRTIET-2014 Conference
Proceeding, 30th -31st August 2014, 1-4
Jiu, Hui Wu, A. Q. Liu, and H. L. Chen (2007).Exact Solutions for FreeVibration Analysis of Rectangular Plates. Journal of Applied Mechanics
Vol. 74 pp. 1247-1251.
Njoku, K. O., Ezeh, J. C., Ibearugbulem, O. M., Ettu, L. O., and Anyaogu,
L. (2013). Free Vibration of Thin Rectangular Isotropic CCCC Plate Using
Taylor Series Formulated Shape Function in Galerkins Method. Academic
Research International, 4 (4): 126-132.
Szilard, R. (2004). Theories and Applications of Plate Analysis. New
Jersey: John Wiley & Sons Inc. Taylor, R. L.and S. Govindjee
(2004).Solution of clamped rectangular plate problems. Communi. Numer.
Meth. Eng. 20: 757765.
150
[9]
Ugural, A. C. (1999). Stresses in plates and shells, 2nd ed. Singapore:

McGraw-hill.
[10] Ventsel, E. and T. Krauthammer (2001).Thin Plates and Shells: Theory,
Analysis and Applications. New York: Marcel Dekker.
[11] Wang, C. M., Y. C. Wang, and J. N. Reddy (2002). Problems and remedy
for the Ritz method in determining stress resultant of corner supported
rectangular plates. Comput.Struct. 80: 145154
[12] Ye, Jianqiao (1994).Large deflection of imperfect plates by iterative BE-FE
method.Journal of Engineering Mechanics, Vol. 120, No. 3 (March).
AUTHORS
First Author Ibearugbulem, Owus M, Civil Engineering
Department, Federal University of Technology, Owerri, Nigeria
Second Author Ibearugbulem, C. N, Civil Engineering
Department, Federal University of Technology, Owerri, Nigeria
Third Author Habib Momoh, Civil Engineering Department,
Imo State University, Owerri, Nigeria
Fourth Author Asomugha, U.C, Civil Engineering
Department, Imo State University, Owerri, Nigeria
www.ijsrp.org

Split-Deflection Method of Classical Rectangular Plate Analysis

Uploaded by

Split-Deflection Method of Classical Rectangular Plate Analysis

Uploaded by

International Journal of Scientific and Research Publications, Volume 6, Issue 5, May 2016

Split-Deflection Method of Classical Rectangular Plate

Ibearugbulem, Owus M., 2Ibearugbulem, C. N, 3Habib Momoh and 3Asomugha, U.C

Civil Engineering Department, Federal University of Technology, Owerri, Nigeria

Abstract- This paper presents an split-deflection method of

The main reason for this modification is to help the analysis

Substituting equations (2) and (3) into equation (1) gives:

he common energy methods in classical plate theory analysis

STRAIN DEFLECTION RELATIONSHIP

STRESS STRAIN RELATIONSHIP

STRESS DEFLECTION RELATIONSHIP

Now, equation (20) can be written in non dimensional axes

TOTAL POTENTIAL ENERGY

For pure bending analysis , the external work is given as:

Where a, b and P are the plate lengths in x and y axes and

Substituting equations (10) to (15) into equation (16) gives

DIRECT VARIATION OF TOTAL POTENTIAL ENERGY

Subtracting equation (17) from Equation (18) gives the total

Substituting equations (1) and (2) into equation (19) gives:

This equation (25) is the direct governing equation of

II. RESULTS AND CONCLUSIONS

From equations (30) (31), h 1 and h 2 are:

With these we obtain:

The result of this present study is presented on table 1. A

Substituting equations (34) to ( 29) into equation (37) gives

The maximum deflection shall occur at the center of the

Deflection at center (qa4/D)

Erdem, C. Imrak and Ismail Gerdemeli (2007).The problem of isotropic

Hutchinson, J. R. (1992). On the bending of rectangular plates with two

Ugural, A. C. (1999). Stresses in plates and shells, 2nd ed. Singapore:

You might also like