SlideShare a Scribd company logo
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
_______________________________________________________________________________________
Volume: 02 Issue: 12 | Dec-2013, Available @ http://www.ijret.org 652
THE USE OF POLYNOMIAL SHAPE FUNCTION IN THE BUCKLING
ANALYSIS OF CCFC RECTANGULAR PLATE
Ezeh, J. C.1
, Ibearugbulem, O. M.2
, Nwadike, A.N.3
, Maduh, U. J.4
1
Associate Professor, Civil Engineering Department, Federal Unitech Owerri, Imo State, Nigeria,
jcezeh2003@yahoo.com
2
Civil Engineering Department, Federal Unitech Owerri, Imo State, Nigeria, ibeowums@yahoo.co.uk
3
Civil Engineering Department, Federal Unitech Owerri, Imo State, Nigeria, nwadikeamarachukwu@yahoo.com
4
Civil Engineering Department, Federal Unitech Owerri, Imo State, Nigeria,ujmaduh@gmail.com
Abstract
The use of polynomial series function in the buckling analysis of a CCFC is presented. The polynomial series shape function was
truncated at the fifth orthogonal terms, which satisfied all the boundary conditions of the plate to obtain a peculiar shape
function, which was applied in Ritz method. The peculiar shape function is substituted into the total potential energy functional,
which was minimized, and the critical buckling load of the plate was obtained. The critical buckling load is a function of a
coefficient, “K”. The values of K from earlier and the present studies were compared within the range of aspect ratios from 0.1 to
2.0. A graph of critical buckling load against aspect ratio was plotted. It was discovered that for aspect ratios of 0.4, 0.5 and 1.0,
the critical buckling loads coefficients were 26.94, 17.39 and 4.83. It was also observed from the behavior of the graph that as
aspect ratio increases from 0.1 to 2.0, the critical buckling load decreases.
Keywords: Total Potential Energy Functional, Shape Function, Polynomial Series Shape Function, Critical Buckling
load, Boundary Condition, Ritz Method.
--------------------------------------------------------------------***----------------------------------------------------------------------
NOTATION
W is the shape function (deflection function).
V is the shear force.
is Poisson’s ratio
First partial derivative of deflection with respect to R is:
′ = .
First partial derivative of deflection with respect to Q is:
′ =
Second partial derivative of deflection with respect to R and
Q is: ′ =
Second partial derivative of deflection with respect to R
is: ′′ =
Second partial derivative of deflection with respect to Q
is: ′′ =
Third partial derivative of deflection with respect to R and Q
is: " =
=
( )
is flexural rigidity, t is the plate thickness. E is
Young’s modulus.
1.INTRODUCTION
Since the origin of plate, Fourier and trigonometric series
has been used to treat plate problems as shape function.
There are many approaches many scholars employed in
solving the problem of plates like energy method, numerical
method, and equilibrium method. These approaches give
approximate results in most difficult situations.
C.Erdem and Ismail used numerical method to analyze an
isotropic rectangular plate with four clamped edges under
distributed loads of an exact solution of the governing
equation in terms of trigonometric and hyperbolic
function[1]. They also compared the obtained results with
those earlier reported, which shows reasonable agreement
with other available results, but with a simpler and practical
approach. Robert and Govindjee used the classical double
cosine series expansion and an exploit of the Sherman-
Woodburg Formular to solve the problem of clamped
rectangular plate [2]. Ibearugbulem used Taylor-Mclaurin's
series as shape function to analyze the instability of axially
compressed thin rectangular plate with four edges clamped
[5]. He truncated the Taylor-Mclaurins series at the fifth
term which satisfied all the boundary conditions of the plate.
Upon substituting the resulted shape function into the total
potential energy, he obtained the critical buckling load of
the plate. Some scholars used finite element method in
analyzing the problems of rectangular plate like Han and
Petty[3], Tajdari et al [11], An- chien Wu et al [10], Ye [9],
Rakesh [8], Sandeep Singh et al [7], and Ganapathi et al [4]
Although the buckling analysis of rectangular plates has
received the attention of many researchers for several
years, its treatment has left much to be done. Most of the
available solutions do not satisfy exactly the prescribed
boundary conditions or the governing differential agnation
or both. In the past, there were no reference books to which
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
_______________________________________________________________________________________
Volume: 02 Issue: 12 | Dec-2013, Available @ http://www.ijret.org 653
Engineers and students can turn to for a clear and orderly
exposition of the Isotropic rectangular plate with clamped –
clamped – free – clamped edge subjected to dynamic
boundary conditions, which are also known as non-essential
boundary condition. Ibearugbulemsolved CCFC rectangular
plate using Taylor – Mclaurin’s shape function with the aid
of eigen value computer programme[6].
In order to provide the solutions to this problem, this paper
is presented to address the issue. The plate is subjected to in-
plane load in one axis (x-axis) of the principal plane as
shown in figure 1.
2. TOTAL POTENTIAL ENERGY FUNCTIONAL
FOR BUCKLING OF PLATE.
Ibearugbulem [6] derived the total potential energy
functional for a rectangular isotropic plate subjected to in-
plane load in x – direction and the polynomial shape
function as follows;
π
=
D
2b
b
a
!w
#
$ +
a
b
!w
&
$ +
2b
a
!w
# &
$ '
((
∂R ∂Q
−
bN
2a
!w
#
$
((
∂R ∂Q (1)
/ = 0 0 a1R1
. b2Q2
3
24(
3
14(
(2)
Where “a” and “b” are plate dimensions (lengths) along x
and y directions. Nx is the in-plane load in x-direction, D is
flexural rigidity. W is the shape function.π is Total potential
energy functional with load along x – axis; R =
5
; Q =
7
8
; 0 ≤ R ≤ 1; 0 ≤ Q ≤ 1 (R and Q are dimensionless
quantities).
3. SHAPE FUNCTION RESULTING FROM
POLYNOMIAL SERIES.
The boundary conditions for CCFC plate are;
w(# 4 () = 0 ; w (& 4 () = 0; w(# 4 ) = 0 (3)
w #
(# 4 () = 0 ; w &
(& 4 () = 0; w #
(# 4 ) = 0 (4)
V&(& 4 )
= >w′′′&
+ (2 − μ)w # &@(& 4 ) = 0 (5)
Applying these boundary conditions in equation 2 gave;
w = A(R − 2R + R3)(4Q − 4Q + Q3) (6)
4. APPLICATION OF RITZ METHOD
Partial derivative of equation (6) with respect to either R or
Q or both gave the following equations;
w′#
= A(2R − 6R + 4R )(4Q − 4Q + Q3) (7)
w′′#
= A(2 − 12R + 12R )(4Q − 4Q + Q3) (8)
w′&
= A(R − 2R + R3)(8Q − 12Q + 4Q ) (9)
w′′&
= A(R − 2R + R3)(8 − 24Q + 12Q ) (10)
w # &
= A(2R − 6R + 4R )(8Q − 12Q + 4Q ) (11)
Integrating the square of these equations (7), (8), (10), and
(11) partially with respect to R and Q in a closed domain
respectively gave:
!w
#
$ ∂R ∂Q
((
= A (0.019047619)(0.406349206) = 0.007739984876A
(12)
!w
#
$
((
∂R ∂Q
= A (0.8)(0.40634920) = 0.325079365A (13)
(w′′&) ∂R ∂Q
((
= A (1.587301587x10 )(12.8) = 0.02031746A (14)
!w
# &
$
((
∂R ∂Q
= A (0.019047619)(1.219047619) = 0.023219954A
(15)
Substituting equations (12), (13), (14) and (15) into equation
(1), and minimizing it gave the following result;
π
∂A
=
DA
Pa
(0.325079365 + 0.02031746P3
+ 0.023219954P )
−
N A
2P
(0.007739984876) = 0 (16)
Making N the subject of the equation (16) with aspect ratio
of P = a/b
N =
Dπ
b
I
4.255489716
P
+ 0.265968103P
+ 0.303963542J (17)
5. RESULT AND CONCLUSIONS
A graph of critical buckling load against aspect ratio was
plotted. From the graph, it was observed that as the aspect
ratio increases from 0.1 to 2.0, the critical buckling load
decreases. The graph was divided into two segments ranging
from 0.1 to 0.7 In the first segment, the graph is of 6th
degree polynomial with an equation of Y = 21892X6
–
60455X5
+ 68130X4
– 40190X3
+ 13164X2
– 23034X +
b NX NX
y
x
P = a/b
CCFC
a
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
_______________________________________________________________________________________
Volume: 02 Issue: 12 | Dec-2013, Available @ http://www.ijret.org 654
1752. In the second segment, the graph is of 5th degree
polynomial with an equation of Y = -5.704X5
+ 44.58X4
-
139.9X3
+ 223.4X2
– 184.7X + 67.20.
a: First Segment
b: Second Segment
Figure 2: Graph of critical buckling load against aspect
ratio.
The average percentage difference between the solutions
from Ibearugbulem and the present study shown on table 1
was 32.112%. This value is more than the required
percentage because Ibearugbulem used a rigorous computer
matrix program to analyze the plate which has so many
disadvantages. It was observed that as the aspect ratio
increases from 0.1 to 2.0, the present work and
Ibearugbulem start converging.
Table 1: K values for different aspect ratios for CCFC
rectangular plate
Aspect
Ratio
P =
5
8
Z From
Ibearugbul
em (2012)
Z From
Present
Study
Differe
nces
Percentage
Difference
0.1 159.16 425.86 266.69 167.56
0.2 54.24 106.70 52.42 96.72
0.3 35.41 47.61 12.20 34.46
0.4 27.29 26.94 0.35 1.27
0.5 17.73 17.39 0.34 1.90
0.6 12.54 12.22 0.32 2.55
0.7 9.42 9.12 0.30 3.20
0.8 7.41 7.12 0.29 3.88
0.9 6.05 5.77 0.28 4.58
1 5. 08 4.83 0.25 5.01
6. CONCLUSION
The polynomial series as a shape function was applied in
this paper for the buckling behavior of an isotropic
rectangular plate.It was also deduced that the polynomial
series satisfies all the non-essential (Dynamics) boundary
conditions of the graph.
The shape function derived in this approach is a good
alternative to the shape function assumed using
trigonometric series. Therefore, the present method can
accurately predict the critical buckling load of isotropic
rectangular plates.
REFERENCE
[1]. C. ErdemImrak and Ismail Gerdemeli (2006), "The
problem of isotropic rectangular plate with four
clamped edge". Journal of Sadhana, vol. 32, Part 3,
Pp 181-186.
[2]. Robert L.Taylor and Sanjay Govindjee
(2002),"Solution of Clamped Rectangular Plate
Problem”. Journal of Structural Engineering
Mechanic and Materials. Vol.9.
[3]. Han, W.,Petyt, M.(1997)."Geometrically Nonlinear
Vibration Analysis of Thin Rectangular Plates
Using the Hierarchical Finite Element Method".
y = 21892x6 -
60455x5 +
68130x4 -
40190x3 +
13164x2 -
23034x + 1752.
0
50
100
150
200
250
300
350
400
450
0 0.5 1
Criticalbucklingloadcoefficient
Aspect Ratio
CCFC PLATE
y = -5.704x5 + 44.58x4 -
139.9x3 + 223.4x2 - 184.7x +
67.20
0
1
2
3
4
5
6
7
8
9
10
0 1 2 3
Criticalbucklingloadcoefficient
Aspect Ratio
CCFC PLATE
Poly. (CCFC
PLATE)
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
_______________________________________________________________________________________
Volume: 02 Issue: 12 | Dec-2013, Available @ http://www.ijret.org 655
Journals of Computers & Structures, vol. 63, issue
2, Pp 295 - 308.
[4]. Ganapathi, M., Varadan, T.K, and Sharma, B.S
(1991)." Nonlinear Flexural Vibration of
Laminated Orthotropic Plate". Journal of
Computers and Structures. Vol. 39, Issue 6, Pp 685
- 688.
[5]. Ibearugbulem, O.M, Osadebe, N.N., Ezeh, J.C.,
Onwuka, D. O (2013). “Instability of Axially
Compressed CCCC Thin Rectangular Plate Using
Taylor-Mclaurins Series Shape Function on Ritz
Method”. Journal of Academic Research
International, vol. 4, No. 1., Pp 346-351
[6]. Ibearugbulem, O.M. (2012). “Application of a
Direct Variational Principal in Elastic Stability of
Rectangular Flat Thin Plates”. Ph.D. Thesis
submitted to postgraduate school, Federal
University of Technology, Owerri, Nigeria.
[7]. Sandeep Singh et al (2012). “Buckling Analysis of
Thin Rectangular Plates with Cutouts Subjected to
Partial edge Compression Using FEM”. Journal of
Engineering, Design and Technology, vol.10,
issue1.
[8]. Rakesh TimappaNaik (2010).“Elastic Buckling
Studies of Thin Plates and Cold-Formed Steel
Members in Shear”. Report Submitted to the
Faculty of Virginia Polytechnic Institute and State
University in Partial Fulfillment of the
Requirements for the Degree of Master of Science
in Civil Engineering.
[9]. Ye Jianqiao (1994). “Large Deflection of Imperfect
Plates by Iterative BE-FE Method”. Journal of
Engineering Mechanics, vol. 120, No. 3,Pp 431-
445.
[10]. An-Chien Wu, Pao-Chun Lin and Keg-Chyuan Tsai
(2013), “High-Mode Buckling-Restrained Brace
Core Plate”,Journals of the International
Association for Earthquake Engineering.
[11]. Tajdari, M., Nezamabadi, A.R., Naeemi, M and
Pirali, P (2011), “The Effects of Plate-Support
Condition on Buckling Strength of Rectangular
Perforated Plates under Linearly Varying in-plane
Normal Load”, Journal of World Academy of
Science, Engineering and Technology, vol. 54, Pp
479-486.

More Related Content

PDF
The use of polynomial shape function in the buckling analysis of ccfc rectang...
PDF
Stability analysis of orthotropic reinforce concrete shear wall
PDF
Stability analysis of orthotropic reinforce concrete shear wall
PDF
An approximate solution for plates resting on winkler foundation
PDF
DESIGN PARAMETERS OF A DYNAMIC VIBRATION ABSORBER WITH TWO SPRINGS IN PARALLEL
PDF
Galerkin’s indirect variational method in elastic stability analysis of all e...
PDF
Characteristic orthogonal polynimial application to galerkin indirect variati...
PPT
Chapter 2
The use of polynomial shape function in the buckling analysis of ccfc rectang...
Stability analysis of orthotropic reinforce concrete shear wall
Stability analysis of orthotropic reinforce concrete shear wall
An approximate solution for plates resting on winkler foundation
DESIGN PARAMETERS OF A DYNAMIC VIBRATION ABSORBER WITH TWO SPRINGS IN PARALLEL
Galerkin’s indirect variational method in elastic stability analysis of all e...
Characteristic orthogonal polynimial application to galerkin indirect variati...
Chapter 2

What's hot (20)

PDF
Numerical modeling to evaluate pile head deflection under the lateral load
PDF
A study on mhd boundary layer flow over a nonlinear
PDF
Buckling analysis of line continuum with new matrices of stiffness and geometry
PPTX
flexibility method
PDF
Module4 plastic theory- rajesh sir
PDF
Module3 direct stiffness- rajesh sir
PDF
Module1 1 introduction-tomatrixms - rajesh sir
PDF
Evaluation of compressive strength of cement using rayleigh’s dimensional ana...
PDF
Composite_AirFrame_Structures_LeeCatherineRamsay_1103072R
PDF
Numerical study of natural convection in an enclosed square cavity using cons...
PDF
Radiation patterns account of a circular microstrip antenna loaded two annular
PDF
Module2 stiffness- rajesh sir
PDF
PREDICTION OF COMPRESSIVE STRENGTH OF CONCRETE FROM EARLY AGE TEST RESULT
PDF
Predicting 28 Days Compressive Strength of Concrete from 7 Days Test Result
PDF
Module3 rajesh sir
PDF
PREDICTION OF COMPRESSIVE STRENGTH OF HIGH PERFORMANCE CONCRETE CONTAINING IN...
PDF
Computational fracture mechanics
PDF
Coupling_of_IGA_plates_and_3D_FEM_domain
PDF
Mechanics of structures - module1
PPT
Transient response of delaminated composite shell subjected to low velocity o...
Numerical modeling to evaluate pile head deflection under the lateral load
A study on mhd boundary layer flow over a nonlinear
Buckling analysis of line continuum with new matrices of stiffness and geometry
flexibility method
Module4 plastic theory- rajesh sir
Module3 direct stiffness- rajesh sir
Module1 1 introduction-tomatrixms - rajesh sir
Evaluation of compressive strength of cement using rayleigh’s dimensional ana...
Composite_AirFrame_Structures_LeeCatherineRamsay_1103072R
Numerical study of natural convection in an enclosed square cavity using cons...
Radiation patterns account of a circular microstrip antenna loaded two annular
Module2 stiffness- rajesh sir
PREDICTION OF COMPRESSIVE STRENGTH OF CONCRETE FROM EARLY AGE TEST RESULT
Predicting 28 Days Compressive Strength of Concrete from 7 Days Test Result
Module3 rajesh sir
PREDICTION OF COMPRESSIVE STRENGTH OF HIGH PERFORMANCE CONCRETE CONTAINING IN...
Computational fracture mechanics
Coupling_of_IGA_plates_and_3D_FEM_domain
Mechanics of structures - module1
Transient response of delaminated composite shell subjected to low velocity o...
Ad

Viewers also liked (20)

PDF
The pattern and realization of zigbee wi-fi
PDF
Monitoring and control of grain storage using plc
PDF
The velocity control of the electro hydraulic servo system
PDF
Review on tls or ssl session sharing based web cluster load balancing
PDF
Digital signal processing techniques for lti fiber
PDF
Studies on seasonal variation of ground water quality
PDF
An overview of methods for monitoring web services based on the quality of se...
PDF
Supervised machine learning based dynamic estimation
PDF
Usability guidelines for usable user interface
PDF
Seismic pounding of multistoreyed buildings
PDF
Natural disasters vulnerability assessment of gondia district, maharashtra, i...
PDF
Correlation measure for intuitionistic fuzzy multi sets
PDF
Cfd analaysis of flow through a conical exhaust
PDF
A novel rrw framework to resist accidental attacks
PDF
Orchestration of web services based on t qo s using user and web services agent
PDF
Statistical process control
PDF
Brisk and secure ad hoc vehicular communication
PDF
Performance evaluation of uasbr for milk processing
PDF
Arsenic removal from water using activated carbon
PDF
Video inpainting using backgroung registration
The pattern and realization of zigbee wi-fi
Monitoring and control of grain storage using plc
The velocity control of the electro hydraulic servo system
Review on tls or ssl session sharing based web cluster load balancing
Digital signal processing techniques for lti fiber
Studies on seasonal variation of ground water quality
An overview of methods for monitoring web services based on the quality of se...
Supervised machine learning based dynamic estimation
Usability guidelines for usable user interface
Seismic pounding of multistoreyed buildings
Natural disasters vulnerability assessment of gondia district, maharashtra, i...
Correlation measure for intuitionistic fuzzy multi sets
Cfd analaysis of flow through a conical exhaust
A novel rrw framework to resist accidental attacks
Orchestration of web services based on t qo s using user and web services agent
Statistical process control
Brisk and secure ad hoc vehicular communication
Performance evaluation of uasbr for milk processing
Arsenic removal from water using activated carbon
Video inpainting using backgroung registration
Ad

Similar to The use of polynomial shape function in the buckling (20)

PDF
Static and Dynamic Reanalysis of Tapered Beam
PDF
Computer Based Free Vibration Analysis of Isotropic Thin Rectangular Flat CCC...
PDF
25 johnarry tonye ngoji 250-263
PDF
Dp36703707
PDF
Comparative Study of an Educational Building by Linear Static Analysis and Re...
PDF
Static analysis of c s short cylindrical shell under internal liquid pressure...
PDF
Static analysis of c s short cylindrical shell under internal liquid pressure...
PDF
A parametric study of x and v bracing industrial steel structure
PDF
Compósitos de Polimeros - Artigo
PPT
FEM 7 Beams and Plates.ppt
PDF
A COMPARATIVE STUDY OF VARIOUS METHODS TO EVALUATE IMPEDANCE FUNCTION FOR SHA...
PDF
IRJET- Probabilistic Risk Analysis of Seismic Irregular RC Structure using Fr...
PDF
Limit States Solution to CSCS Orthotropic Thin Rectangular Plate Carrying Tra...
PDF
A study on mhd boundary layer flow over a nonlinear stretching sheet using im...
PDF
Performance measurement and dynamic analysis of two dof robotic arm manipulator
PDF
40120130406008
PDF
The numerical solution of helmholtz equation via multivariate padé approximation
PDF
Kk2518251830
PDF
Kk2518251830
PDF
techDynamic characteristics and stability of cylindrical textured journal bea...
Static and Dynamic Reanalysis of Tapered Beam
Computer Based Free Vibration Analysis of Isotropic Thin Rectangular Flat CCC...
25 johnarry tonye ngoji 250-263
Dp36703707
Comparative Study of an Educational Building by Linear Static Analysis and Re...
Static analysis of c s short cylindrical shell under internal liquid pressure...
Static analysis of c s short cylindrical shell under internal liquid pressure...
A parametric study of x and v bracing industrial steel structure
Compósitos de Polimeros - Artigo
FEM 7 Beams and Plates.ppt
A COMPARATIVE STUDY OF VARIOUS METHODS TO EVALUATE IMPEDANCE FUNCTION FOR SHA...
IRJET- Probabilistic Risk Analysis of Seismic Irregular RC Structure using Fr...
Limit States Solution to CSCS Orthotropic Thin Rectangular Plate Carrying Tra...
A study on mhd boundary layer flow over a nonlinear stretching sheet using im...
Performance measurement and dynamic analysis of two dof robotic arm manipulator
40120130406008
The numerical solution of helmholtz equation via multivariate padé approximation
Kk2518251830
Kk2518251830
techDynamic characteristics and stability of cylindrical textured journal bea...

More from eSAT Publishing House (20)

PDF
Likely impacts of hudhud on the environment of visakhapatnam
PDF
Impact of flood disaster in a drought prone area – case study of alampur vill...
PDF
Hudhud cyclone – a severe disaster in visakhapatnam
PDF
Groundwater investigation using geophysical methods a case study of pydibhim...
PDF
Flood related disasters concerned to urban flooding in bangalore, india
PDF
Enhancing post disaster recovery by optimal infrastructure capacity building
PDF
Effect of lintel and lintel band on the global performance of reinforced conc...
PDF
Wind damage to trees in the gitam university campus at visakhapatnam by cyclo...
PDF
Wind damage to buildings, infrastrucuture and landscape elements along the be...
PDF
Shear strength of rc deep beam panels – a review
PDF
Role of voluntary teams of professional engineers in dissater management – ex...
PDF
Risk analysis and environmental hazard management
PDF
Review study on performance of seismically tested repaired shear walls
PDF
Monitoring and assessment of air quality with reference to dust particles (pm...
PDF
Low cost wireless sensor networks and smartphone applications for disaster ma...
PDF
Coastal zones – seismic vulnerability an analysis from east coast of india
PDF
Can fracture mechanics predict damage due disaster of structures
PDF
Assessment of seismic susceptibility of rc buildings
PDF
A geophysical insight of earthquake occurred on 21 st may 2014 off paradip, b...
PDF
Effect of hudhud cyclone on the development of visakhapatnam as smart and gre...
Likely impacts of hudhud on the environment of visakhapatnam
Impact of flood disaster in a drought prone area – case study of alampur vill...
Hudhud cyclone – a severe disaster in visakhapatnam
Groundwater investigation using geophysical methods a case study of pydibhim...
Flood related disasters concerned to urban flooding in bangalore, india
Enhancing post disaster recovery by optimal infrastructure capacity building
Effect of lintel and lintel band on the global performance of reinforced conc...
Wind damage to trees in the gitam university campus at visakhapatnam by cyclo...
Wind damage to buildings, infrastrucuture and landscape elements along the be...
Shear strength of rc deep beam panels – a review
Role of voluntary teams of professional engineers in dissater management – ex...
Risk analysis and environmental hazard management
Review study on performance of seismically tested repaired shear walls
Monitoring and assessment of air quality with reference to dust particles (pm...
Low cost wireless sensor networks and smartphone applications for disaster ma...
Coastal zones – seismic vulnerability an analysis from east coast of india
Can fracture mechanics predict damage due disaster of structures
Assessment of seismic susceptibility of rc buildings
A geophysical insight of earthquake occurred on 21 st may 2014 off paradip, b...
Effect of hudhud cyclone on the development of visakhapatnam as smart and gre...

Recently uploaded (20)

PPTX
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
PDF
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
PPTX
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
PDF
Operating System & Kernel Study Guide-1 - converted.pdf
PPTX
KTU 2019 -S7-MCN 401 MODULE 2-VINAY.pptx
PPTX
Infosys Presentation by1.Riyan Bagwan 2.Samadhan Naiknavare 3.Gaurav Shinde 4...
PPTX
FINAL REVIEW FOR COPD DIANOSIS FOR PULMONARY DISEASE.pptx
PDF
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
PPTX
UNIT-1 - COAL BASED THERMAL POWER PLANTS
PPTX
CYBER-CRIMES AND SECURITY A guide to understanding
PDF
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
PPTX
Geodesy 1.pptx...............................................
PPTX
Welding lecture in detail for understanding
PPTX
UNIT 4 Total Quality Management .pptx
PDF
SM_6th-Sem__Cse_Internet-of-Things.pdf IOT
PDF
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
PDF
The CXO Playbook 2025 – Future-Ready Strategies for C-Suite Leaders Cerebrai...
PPTX
Lecture Notes Electrical Wiring System Components
PDF
Structs to JSON How Go Powers REST APIs.pdf
PPTX
Internet of Things (IOT) - A guide to understanding
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
Operating System & Kernel Study Guide-1 - converted.pdf
KTU 2019 -S7-MCN 401 MODULE 2-VINAY.pptx
Infosys Presentation by1.Riyan Bagwan 2.Samadhan Naiknavare 3.Gaurav Shinde 4...
FINAL REVIEW FOR COPD DIANOSIS FOR PULMONARY DISEASE.pptx
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
UNIT-1 - COAL BASED THERMAL POWER PLANTS
CYBER-CRIMES AND SECURITY A guide to understanding
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
Geodesy 1.pptx...............................................
Welding lecture in detail for understanding
UNIT 4 Total Quality Management .pptx
SM_6th-Sem__Cse_Internet-of-Things.pdf IOT
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
The CXO Playbook 2025 – Future-Ready Strategies for C-Suite Leaders Cerebrai...
Lecture Notes Electrical Wiring System Components
Structs to JSON How Go Powers REST APIs.pdf
Internet of Things (IOT) - A guide to understanding

The use of polynomial shape function in the buckling

  • 1. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 _______________________________________________________________________________________ Volume: 02 Issue: 12 | Dec-2013, Available @ http://www.ijret.org 652 THE USE OF POLYNOMIAL SHAPE FUNCTION IN THE BUCKLING ANALYSIS OF CCFC RECTANGULAR PLATE Ezeh, J. C.1 , Ibearugbulem, O. M.2 , Nwadike, A.N.3 , Maduh, U. J.4 1 Associate Professor, Civil Engineering Department, Federal Unitech Owerri, Imo State, Nigeria, jcezeh2003@yahoo.com 2 Civil Engineering Department, Federal Unitech Owerri, Imo State, Nigeria, ibeowums@yahoo.co.uk 3 Civil Engineering Department, Federal Unitech Owerri, Imo State, Nigeria, nwadikeamarachukwu@yahoo.com 4 Civil Engineering Department, Federal Unitech Owerri, Imo State, Nigeria,ujmaduh@gmail.com Abstract The use of polynomial series function in the buckling analysis of a CCFC is presented. The polynomial series shape function was truncated at the fifth orthogonal terms, which satisfied all the boundary conditions of the plate to obtain a peculiar shape function, which was applied in Ritz method. The peculiar shape function is substituted into the total potential energy functional, which was minimized, and the critical buckling load of the plate was obtained. The critical buckling load is a function of a coefficient, “K”. The values of K from earlier and the present studies were compared within the range of aspect ratios from 0.1 to 2.0. A graph of critical buckling load against aspect ratio was plotted. It was discovered that for aspect ratios of 0.4, 0.5 and 1.0, the critical buckling loads coefficients were 26.94, 17.39 and 4.83. It was also observed from the behavior of the graph that as aspect ratio increases from 0.1 to 2.0, the critical buckling load decreases. Keywords: Total Potential Energy Functional, Shape Function, Polynomial Series Shape Function, Critical Buckling load, Boundary Condition, Ritz Method. --------------------------------------------------------------------***---------------------------------------------------------------------- NOTATION W is the shape function (deflection function). V is the shear force. is Poisson’s ratio First partial derivative of deflection with respect to R is: ′ = . First partial derivative of deflection with respect to Q is: ′ = Second partial derivative of deflection with respect to R and Q is: ′ = Second partial derivative of deflection with respect to R is: ′′ = Second partial derivative of deflection with respect to Q is: ′′ = Third partial derivative of deflection with respect to R and Q is: " = = ( ) is flexural rigidity, t is the plate thickness. E is Young’s modulus. 1.INTRODUCTION Since the origin of plate, Fourier and trigonometric series has been used to treat plate problems as shape function. There are many approaches many scholars employed in solving the problem of plates like energy method, numerical method, and equilibrium method. These approaches give approximate results in most difficult situations. C.Erdem and Ismail used numerical method to analyze an isotropic rectangular plate with four clamped edges under distributed loads of an exact solution of the governing equation in terms of trigonometric and hyperbolic function[1]. They also compared the obtained results with those earlier reported, which shows reasonable agreement with other available results, but with a simpler and practical approach. Robert and Govindjee used the classical double cosine series expansion and an exploit of the Sherman- Woodburg Formular to solve the problem of clamped rectangular plate [2]. Ibearugbulem used Taylor-Mclaurin's series as shape function to analyze the instability of axially compressed thin rectangular plate with four edges clamped [5]. He truncated the Taylor-Mclaurins series at the fifth term which satisfied all the boundary conditions of the plate. Upon substituting the resulted shape function into the total potential energy, he obtained the critical buckling load of the plate. Some scholars used finite element method in analyzing the problems of rectangular plate like Han and Petty[3], Tajdari et al [11], An- chien Wu et al [10], Ye [9], Rakesh [8], Sandeep Singh et al [7], and Ganapathi et al [4] Although the buckling analysis of rectangular plates has received the attention of many researchers for several years, its treatment has left much to be done. Most of the available solutions do not satisfy exactly the prescribed boundary conditions or the governing differential agnation or both. In the past, there were no reference books to which
  • 2. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 _______________________________________________________________________________________ Volume: 02 Issue: 12 | Dec-2013, Available @ http://www.ijret.org 653 Engineers and students can turn to for a clear and orderly exposition of the Isotropic rectangular plate with clamped – clamped – free – clamped edge subjected to dynamic boundary conditions, which are also known as non-essential boundary condition. Ibearugbulemsolved CCFC rectangular plate using Taylor – Mclaurin’s shape function with the aid of eigen value computer programme[6]. In order to provide the solutions to this problem, this paper is presented to address the issue. The plate is subjected to in- plane load in one axis (x-axis) of the principal plane as shown in figure 1. 2. TOTAL POTENTIAL ENERGY FUNCTIONAL FOR BUCKLING OF PLATE. Ibearugbulem [6] derived the total potential energy functional for a rectangular isotropic plate subjected to in- plane load in x – direction and the polynomial shape function as follows; π = D 2b b a !w # $ + a b !w & $ + 2b a !w # & $ ' (( ∂R ∂Q − bN 2a !w # $ (( ∂R ∂Q (1) / = 0 0 a1R1 . b2Q2 3 24( 3 14( (2) Where “a” and “b” are plate dimensions (lengths) along x and y directions. Nx is the in-plane load in x-direction, D is flexural rigidity. W is the shape function.π is Total potential energy functional with load along x – axis; R = 5 ; Q = 7 8 ; 0 ≤ R ≤ 1; 0 ≤ Q ≤ 1 (R and Q are dimensionless quantities). 3. SHAPE FUNCTION RESULTING FROM POLYNOMIAL SERIES. The boundary conditions for CCFC plate are; w(# 4 () = 0 ; w (& 4 () = 0; w(# 4 ) = 0 (3) w # (# 4 () = 0 ; w & (& 4 () = 0; w # (# 4 ) = 0 (4) V&(& 4 ) = >w′′′& + (2 − μ)w # &@(& 4 ) = 0 (5) Applying these boundary conditions in equation 2 gave; w = A(R − 2R + R3)(4Q − 4Q + Q3) (6) 4. APPLICATION OF RITZ METHOD Partial derivative of equation (6) with respect to either R or Q or both gave the following equations; w′# = A(2R − 6R + 4R )(4Q − 4Q + Q3) (7) w′′# = A(2 − 12R + 12R )(4Q − 4Q + Q3) (8) w′& = A(R − 2R + R3)(8Q − 12Q + 4Q ) (9) w′′& = A(R − 2R + R3)(8 − 24Q + 12Q ) (10) w # & = A(2R − 6R + 4R )(8Q − 12Q + 4Q ) (11) Integrating the square of these equations (7), (8), (10), and (11) partially with respect to R and Q in a closed domain respectively gave: !w # $ ∂R ∂Q (( = A (0.019047619)(0.406349206) = 0.007739984876A (12) !w # $ (( ∂R ∂Q = A (0.8)(0.40634920) = 0.325079365A (13) (w′′&) ∂R ∂Q (( = A (1.587301587x10 )(12.8) = 0.02031746A (14) !w # & $ (( ∂R ∂Q = A (0.019047619)(1.219047619) = 0.023219954A (15) Substituting equations (12), (13), (14) and (15) into equation (1), and minimizing it gave the following result; π ∂A = DA Pa (0.325079365 + 0.02031746P3 + 0.023219954P ) − N A 2P (0.007739984876) = 0 (16) Making N the subject of the equation (16) with aspect ratio of P = a/b N = Dπ b I 4.255489716 P + 0.265968103P + 0.303963542J (17) 5. RESULT AND CONCLUSIONS A graph of critical buckling load against aspect ratio was plotted. From the graph, it was observed that as the aspect ratio increases from 0.1 to 2.0, the critical buckling load decreases. The graph was divided into two segments ranging from 0.1 to 0.7 In the first segment, the graph is of 6th degree polynomial with an equation of Y = 21892X6 – 60455X5 + 68130X4 – 40190X3 + 13164X2 – 23034X + b NX NX y x P = a/b CCFC a
  • 3. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 _______________________________________________________________________________________ Volume: 02 Issue: 12 | Dec-2013, Available @ http://www.ijret.org 654 1752. In the second segment, the graph is of 5th degree polynomial with an equation of Y = -5.704X5 + 44.58X4 - 139.9X3 + 223.4X2 – 184.7X + 67.20. a: First Segment b: Second Segment Figure 2: Graph of critical buckling load against aspect ratio. The average percentage difference between the solutions from Ibearugbulem and the present study shown on table 1 was 32.112%. This value is more than the required percentage because Ibearugbulem used a rigorous computer matrix program to analyze the plate which has so many disadvantages. It was observed that as the aspect ratio increases from 0.1 to 2.0, the present work and Ibearugbulem start converging. Table 1: K values for different aspect ratios for CCFC rectangular plate Aspect Ratio P = 5 8 Z From Ibearugbul em (2012) Z From Present Study Differe nces Percentage Difference 0.1 159.16 425.86 266.69 167.56 0.2 54.24 106.70 52.42 96.72 0.3 35.41 47.61 12.20 34.46 0.4 27.29 26.94 0.35 1.27 0.5 17.73 17.39 0.34 1.90 0.6 12.54 12.22 0.32 2.55 0.7 9.42 9.12 0.30 3.20 0.8 7.41 7.12 0.29 3.88 0.9 6.05 5.77 0.28 4.58 1 5. 08 4.83 0.25 5.01 6. CONCLUSION The polynomial series as a shape function was applied in this paper for the buckling behavior of an isotropic rectangular plate.It was also deduced that the polynomial series satisfies all the non-essential (Dynamics) boundary conditions of the graph. The shape function derived in this approach is a good alternative to the shape function assumed using trigonometric series. Therefore, the present method can accurately predict the critical buckling load of isotropic rectangular plates. REFERENCE [1]. C. ErdemImrak and Ismail Gerdemeli (2006), "The problem of isotropic rectangular plate with four clamped edge". Journal of Sadhana, vol. 32, Part 3, Pp 181-186. [2]. Robert L.Taylor and Sanjay Govindjee (2002),"Solution of Clamped Rectangular Plate Problem”. Journal of Structural Engineering Mechanic and Materials. Vol.9. [3]. Han, W.,Petyt, M.(1997)."Geometrically Nonlinear Vibration Analysis of Thin Rectangular Plates Using the Hierarchical Finite Element Method". y = 21892x6 - 60455x5 + 68130x4 - 40190x3 + 13164x2 - 23034x + 1752. 0 50 100 150 200 250 300 350 400 450 0 0.5 1 Criticalbucklingloadcoefficient Aspect Ratio CCFC PLATE y = -5.704x5 + 44.58x4 - 139.9x3 + 223.4x2 - 184.7x + 67.20 0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 Criticalbucklingloadcoefficient Aspect Ratio CCFC PLATE Poly. (CCFC PLATE)
  • 4. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 _______________________________________________________________________________________ Volume: 02 Issue: 12 | Dec-2013, Available @ http://www.ijret.org 655 Journals of Computers & Structures, vol. 63, issue 2, Pp 295 - 308. [4]. Ganapathi, M., Varadan, T.K, and Sharma, B.S (1991)." Nonlinear Flexural Vibration of Laminated Orthotropic Plate". Journal of Computers and Structures. Vol. 39, Issue 6, Pp 685 - 688. [5]. Ibearugbulem, O.M, Osadebe, N.N., Ezeh, J.C., Onwuka, D. O (2013). “Instability of Axially Compressed CCCC Thin Rectangular Plate Using Taylor-Mclaurins Series Shape Function on Ritz Method”. Journal of Academic Research International, vol. 4, No. 1., Pp 346-351 [6]. Ibearugbulem, O.M. (2012). “Application of a Direct Variational Principal in Elastic Stability of Rectangular Flat Thin Plates”. Ph.D. Thesis submitted to postgraduate school, Federal University of Technology, Owerri, Nigeria. [7]. Sandeep Singh et al (2012). “Buckling Analysis of Thin Rectangular Plates with Cutouts Subjected to Partial edge Compression Using FEM”. Journal of Engineering, Design and Technology, vol.10, issue1. [8]. Rakesh TimappaNaik (2010).“Elastic Buckling Studies of Thin Plates and Cold-Formed Steel Members in Shear”. Report Submitted to the Faculty of Virginia Polytechnic Institute and State University in Partial Fulfillment of the Requirements for the Degree of Master of Science in Civil Engineering. [9]. Ye Jianqiao (1994). “Large Deflection of Imperfect Plates by Iterative BE-FE Method”. Journal of Engineering Mechanics, vol. 120, No. 3,Pp 431- 445. [10]. An-Chien Wu, Pao-Chun Lin and Keg-Chyuan Tsai (2013), “High-Mode Buckling-Restrained Brace Core Plate”,Journals of the International Association for Earthquake Engineering. [11]. Tajdari, M., Nezamabadi, A.R., Naeemi, M and Pirali, P (2011), “The Effects of Plate-Support Condition on Buckling Strength of Rectangular Perforated Plates under Linearly Varying in-plane Normal Load”, Journal of World Academy of Science, Engineering and Technology, vol. 54, Pp 479-486.