A.- .. ._. ‘_____. _‘r -
LISRARY
Michigan State
University
This is to certify that the
thesis entitled
PRE- AND POST-BUCKLING BEHAVIOR
OF PLATES OF VARIABLE STIFFNESS
USING FINITE DIFFERENCES
presented by
Mohammad Ali Barkhordari
has been accepted towards fulfillment
of the requirements for
__£h.._D_.__degree in Weering
( . ’D
a . _ ) , . ,
“21/ Aft (we a (34;..[4'1
Major professor <
Date December 24L 1980
0-7 639
_4_ __1.7. '
OVERDU FINES:
..J ' 25¢ pen-dumm-
‘ if" '\\\\ f , RETURNXQ ugnm M'rsagus:
Place in book return to move
chum from circulation records
‘KJ;\ -“‘“"” »-
PRE-AND POST-BUCKLING BEHAVIOR
OF PLATES OF VARIABLE STIFFNESS
USING FINITE DIFFERENCES
By
Mohammad Ali Barkhordari
A DISSERTATION
Submitted to
Michigan State University
in partial fulfillment of the requirements
for the degree of
DOCTOR OF PHILOSOPHY
Department of Civil and Sanitary Engineering
1980
ABSTRACT
PRE-AND POST-BUCKLING BEHAVIOR
OF PLATES OF VARIABLE STIFFNESS
USING FINITE DIFFERENCES
By
M. Ali Barkhordari
The Von Karman large deflection equation is applied to plates
of variable stiffness. Equilibrium equations and the in-plane
compatibility equation are derived. The ordinary finite difference
technique is employed to solve the nonlinear coupled partial differental
e(Nations. Iva different methods of formulation are considered:
a) In terms of lateral displacement w and a stress function,
5) In terms of displacement components u, v and w.
Stiffness variation can be implemented in two different ways, either
by varYiug the thickness of the plate with constant E, or by taking
a uniform thickness plate of variable E. Both types of stiffness
variation are considered. The nature of in-plane displacements on
the boundary is a significant factor in postbuckling. This effect
is examined by considering plates with different in-plane displacement
boundary conditions.
Several problems with different stiffness variation and
bOundary conditions are solved. The applicable computer program is
utilized to carry out the numerical solutions. In each case the
problem is investigated for different stages of loading as follows:
o
«I . ..- p,
n
I- u.-. o -
'
.. ,
..: ‘.e..':.
“-0.1.-
0...... 'v"
. ' -
Ks “
“" s, o
a
“..
“OJ 1
, .
F “V‘ a
Nu
m A. -
'
a
on-
-
.~‘
in.
. ‘
N . :‘ u',
H‘
a.
.‘4
‘o
v “.
. < _,
5“
...:~-._
‘I
a). Membrane solution analyzes the behavior of in—plane
forces and displacements for undeflected plates,
b) Stability analysis investigates the buckling and effect
of stiffness variation on critical loads and buckling
modes,
c) Postbuckling discusses the behavior of various aspects
of the problem due to edge loads or displacements higher
than critical values.
For clarity, the results are always accompanied by graphical illustrations
of membrane and bending stress as well as displacement components.
The accuracy of the solution is evaluated by comparison of the results
obtained with results from past studies and exact results, where
these results are available. The influence of the grid-spacing on
the accuracy of the results is investigated by taking successively
finer grid-spacings. The numerical results are analyzed and the effect
of stiffness variation on different aspects of the problem discussed.
One objective is to design a plate with stiffness variation
such that it be optimum in some respect. Some possible cases of
optimization are discussed and,as examples,some problems related to
buckling are solved. The results indicate that a considerable weight
and/or material savings can be achieved by using an efficient stiffness
variation pattern.
ACKNOWLEDGMENTS
The author wishes to express his most sincere gratitude to
Ids major Professor, Dr. W. A. Bradley, Professor of civil engineering
for his encouragement and constant help and guidance in the author's
academic development and preparation of this dissertation.
Thanks are also expressed to the other members of the guidance
committee - Dr. N. Altiero, Dr. J. L. Lubkin and Dr. R. K. Wen -
for their guidance and encouragement. The author also owes his
appreciation to Mrs. Clara Hanna for typing of this dissertation.
Special appreciation is also due his wife and son for their patience
and understanding.
ii
TABLES OF CONTENTS
Page
LIST OF TABLES-o.oooooooooo0.0000000000000000. v
LIST OF FIGURESOOOCOOOOOOOIOOOOOOOOOOOOOOO0.00 Vi
Chapter
I INTRODUCTION.................................. 1
1.1 General Remarks.......................... 1
1.2 Previous Developements................... 3
1.3 Present Investigation.............. ..... . 8
1.4 Notations................................ 10
II THEORETICAL DERIVATIONS 13
2.1 General.................................. 13
2.1.1 Nonlinear Equilibrium Equations... 14
2.1.2 Relation Between Stress Resultants
and Displacements................. 18
2.2 Formulation in Terms of Stress Function
and w.................................... 19
2.3 Formulation in Terms of Displacements
u, v and w............................... 21
2.4 Finite Difference Approximations......... 22
2.4.1 Principles of Finite Differences.. 22
2.4.2 Finite Difference Approximation
to Method Discussed in Section
(2.2)............................. 26
2.4.3 Finite Difference Approximation of
Method Discussed‘fin Section (2.3). 37
2.5 Boundary Conditions...................... 45
2.5.1 Some Examples of Practical B.C.... 48
2.6 Summary.................................. 48
2.6.1 Membrane Solution................. 49
2.6.2 Lateral Loading................... 49
2.6.3 Stability Analysis................ 50
2.6.4 Postbuckling...................... 51
III APPLICATION...................................' 52
3.1 Force Boundary Conditions................ 56
3.1.1 membrane Solution................. 58
3.1.1.1 Analysis of Results From
Membrane Solution......... 70
3.1.2 Buckling.......................... 77
iii
Chapter Page
3 2.1 General Buckling.............. 81
3.2.2 Analysis of the Results. ...... 88
3.1.3 Postbuckling. ...................... 93
3. .3.1 General Procedure. ............ 93
3. 3. 2 Numerical Solution............ 96
a) Square Plate with R a 1/10. 96
b) Square Plate with Uniform
Stiffness, R = 1........... 104
c ,d) Square Plate with R = 1/2
and R = 10. ......... ..... 110
3.1.3.3 Analysis of Postbuckling
Results....................... 119
3.1.3.3.1 Simply-Supported Edges...... 119
a) Uniform Stiffness Plate.... 119
b) Variable Stiffness, R = 1/10 124
c) Variable Stiffness, R 8 1/2 126
d) Variable Stiffness, R = 10.. 126
3.1.3.3.2 Clamped Edges................ 129
3.1.4 Optimization........................ 147
3.1.5 Summary............................. 154
3.2 Displacement Boundary Condition........... 156
3.2.1 Membrane Solution.................. 159
3. 2. 2 Buckling Solution.................. 164
3.2.2.1 Convergence Check and
Comparison.................... 164
.1.
l.
Idld
3.2.2.2 Optimization Analysis......... 169
3.2.2.3 Analysis of Buckling Mbdes.... 172
3. 2.3 Postbuckling....................... 175
3.2.3.1 Uniform Thickness Plate,
Simply Supported.............. 175
3.2.3.2 Simply Supported Square Plate
of Variable Thickness......... 181
3.2.3.3 Uniform Thickness Plate,
Clamped Boundary.............. 189
3.2.3.4 Clamped Plate with Variable
Thickness..................... 194
3.3 Comparison of Two Methods................. 202
IV CONCLUSIONooooooooooooocoo-000......ooooooooso. 204
4.1 The Problem Summary....................... 204
4. 2 Conclusions............................... 206
4.3 Recommendations........................... 208
BIBLIOGRAPE‘IY.oooooooooo'oo0.000000000000000ooooo 210
APPENDICES.°O............°.°.'...Ooooooooooooooo 214
iv
Table
3.1
3.2
3.3
3.4
3.5
3.6
3.7
3.8
3.9
3.10
3.11
3.12
3.13
3.14
LIST OF TABLES
Stress Function and Stress Resultant Ratios
Square Plate with R = 1/10.. ....................
Stress Function and Membrane Forces, R = 1/2....
Stress Function and In-plane Resultant Ratios,
R = 1 ...........................................
Stress Function and Membrane Forces, R = 10 .....
Convergence of the Solution, R = 1/10... ........
Eigenvalues of Simply-Supported Square Plate,
R 8 1/10, with Different Grid Spacings.... ......
Comparison of First Eigenvalues of Different
solutions (R=I)OOOOOOOOOOOCOOOOOOOOOOOOOOOO .....
Critical Loads of Simply-Supported Square Plate
R81 ...... 0.0... ........ .0 OOOOOOOOOOOOOOOOOOOOO
Critical Load of Clamped Square Plate Under
Bi-axial Uniform Load, R = 1............ ........
Boundary Displacement for Square Plate of
Figure 3.62................ .....................
Critical Displacements for a Simply-Supported
Plate using Different Mesh Sizes, Uniform
Stiffness Plate ..... ........ ....................
Critical Displacements for a Clamped Square
Plate using Different Mesh Sizes, Uniform
Stiffness .......... . ........... . ................
Convergence of Postbuckling Solution with Mesh
Size. Simply-Supported Plate...................
Convergence of Postbuckling Solution with Mesh
Size. Clamped PlateOOOIOOO.COOOOOOOOIOOOOOOOOOO
Page
65
67
68
68
69
80
80
82
91
158
165
166
176
189
-170
no"
.\'
q...
.o..
a...
“II.
J"‘
g.‘
C“--
t
7
LIST OF FIGURES
Figure Page
2.1 Rectangular Flat Plate .......................... 15
2.2 Plate Element dxdy in Undeformed Configuration. . 15
2.3 Schematic Illustration of Internal Forces and
Moments on the Element of Middle Surface in
Deformed Configuration ........ . ................. 17
2-4 Function f(x) . . ................................. 23
2-5 Nodal Pattern..l........ ......................... 24
2.5 Two Dimensional Operator for fxy ................ 25
2.7 Nodal Arrangement. . . . . . . . . ...................... 27
2.8 Difference Operator for Left Hand Side of
Equilibrium Equation (2.14) ..................... 28
2.9 Difference Operator for Right Hand Side of
Equation (2.14) . . ............................... 30
2.10 Difference Operator for Left Hand Side of
Compatibility Equation (2.15) ................... 33
2.11 Operators for x—Equilibrium Equation (2.16) ..... 39
2.12“ Operators for y-Equilibrium Equation (2.17) ..... 41
3.1 Geometry and Stiffness Variation of Square
Plate-coo ooooooooooooooo Oooooooo ooooooooooooooo o 54
3-2 Stiffness Variation....... ...................... 57
3.3 Geometrical Plan and (9 Function ..... . ........... 61
3.4 Stiffness Ratio at Nodes and Intermediate
Nodes, Square Plate, R = 1/10 ......... . ......... 64
3.5 Convergence of Membrane Solution.. .............. 70
Force Distribution for Undeflected Square
3.6
Plate, R - 1...........
Figure Page
3.7 Contours of In-plane Displacement,
Undeflected Square Plate, R = 1 ................ 72
3.8 Distribution of In-plane Force and Displace-
ment, Square Plate, R = 1/10 ................... 73
3.9 Distribution of In-plane Force and Displace-
ment, Square Plate, R = 1/2 .................... 75
3.10 Distribution of In-plane Force and Displace-
ment, Square Plate, R = 10 ..................... 76
3.11 Convergence of Eigenvalue, R = 1/10 ............ 79
3-12 Modes of Buckling of Square Plate, R = 1/10,
Simply-Supported ....................... . ....... 84
3-13 Modes of Buckling of Square Plate, R = 1/10,
Clamped........‘. ............................... 84
3.14 Modes of Buckling of Square Plate, R = 1/2,
Simply-Supported ........ . ...................... 85
3.15 Modes of Buckling of Square Plate, R = 1/2,
Clamped. . . . . ................................... 85
3.16 Modesof Buckling of Square Plate, R = 1,
Simply-Supported..............................
3.17 Modes of Buckling of Square Plate, R = 1,
Clamped..0.O.C00.0.0...OOOIOOOOOOUOOOCOOO0.0..
3.18 Modes of Buckling of Square Plate, R = 10,
Simply-Supported” . ...................... 87
3.19 Modes of Buckling of Square Plate, R = 10,
Clamped........_ ........................ . ....... 87
3.20 Deflected Shape for s—s and Clamped Plate ...... 90
3.21 Flow Chart of Iterative Procedure .............. 94
3.22 Plots of w, U, Nx and M x’ Square Plate, 8-8,
R8 llloooooooooooo o oooooooo oo ooooooo 000000000. 100
3.23 Plots of Stress Components, Square Plate, s-s,
Rs 1/10....... .......... ..... ......... o ....... 101
3 24 Contours of In-plane Force and Profiles of
N ’ 5-3, R = 1/10.. .......................... 102
x
3.25
3.26
3.27
L28
L29
L30
3.31
L32
L33
L34
L35
L36
3.37
IL38
L39
IL40
3.41
3.42
3.43
In-Plane Displacement, Square Plate, s-s,
R=1/10 ooooooooooooooooooooooo oooooo oooooooooooo
Plots of w, U, Nx and M x’ Square Plate, s-s,
R = 1 ......................................... ...
Plots of Stress Components, Square Plate, s-s,
R=1000 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
Contours of In-Plane Force and Profiles of N x’
s—s, R = 1 .......................................
In-Plane Displacement, Square Plate, s—s, R = 1..
Square Plate Under Uniform Load, s-s, R = 1 ......
Plots of w, U, Nx and M x’ Square Plate, s-s,
R = 1/2 ................................ . .........
Plots of Stress Components, Square Plate, s-s,
Ra 1/20 oooooooo o oooooooooooooooooooooooooooooooo
Contours of In-Plane Force and Profiles of N x’
S-S’R=1/20 ..... 00...... ..... O... ....... O ......
In-Plane Displacement, Square Plate, s-s, R = 1/2.
Plots of w, U, Nx and M x’ Square Plate, s-s,
R=100 0...... 000000 O ..... O OOOOOOOOOOOOOOOOOOOO
Plots of Stress Components, Square Plate, s-s,
R=100000000000 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
Contours of In—Plane Force and Profiles of N x’
8-3 , =10 oooooooooooo o ooooooooooooooooooooooooo
In-Plane Displacement, Square Plate, s-s, R = 10.
Convergence of the Solution with Mesh Size .......
Deflected Shape for Various R Values... ..........
Plots of w, U, Nx and Mx’ Square Plate, Clamped,
R=1/10ooooooo ooooo ooooo oooooooooooooooooooooooo
Plots of Stress Components, Square Plate,
Clamped, R = 1/10. .................... . ..........
Contours of In-Plane Force and Profiles of N x’
Clamped, R = 1/10 ooooo o ooooooooooooooooooo o ......
105
107
108
109
111
112
114
117
118
120
132
133
‘ A
3.44
3.45
3.46
3.47
3.48
3.49
3.50
3.51
3.52
3.53
3.54
3.55
3.56
3.57
3.58
3.59
3.60
3.61
In-Plane Displacement, Square Plate, Clamped,
R - 1/10 ...................... . .................. 134
Plots of w, U, Nx and Mx , Square Plate, Clamped,
RF 1.......... ................................... 135
Plots of Stress Components, Square Plate,
Clamped, R= l ........ . ............ .. ........... 136
Contours of In-Plane Force and Profiles of N x’
Clamped, R . 1 ...................... . .......... 137
In—Plane Displacement, Square Plate, Clamped,
R. 1... .......... ...... ..... . ........... . ..... 138
R-l/z ooooo oooo oooooooooooooooooooo oooo ooooooo 139
Plots of Stress Components, Square Plate, Clamped
R'l/Zooo ..... ........ ............ . ........... 140
Contours of In-Plane Force and Profiles of N ,
Clamped, up 1/2.............................¥.. 141
In-Plane Displacement, Square Plate, Clamped,
R.1/2....OO.......0...O....0.000....0.0...... 142
Plots of w, U, Nx and M x’ Square Plate, Clamped,
R a 10... ............. . ............. . .......... 143
Plots of Stress Components, Square Plate,
Clamped, R a 10... ............. . ...... . ........ 144
Contours of In—Plane Force and Profiles of
Nx’ Clamped, R = 10......... ................... 145
In-Plane Displacement, Square Plate, Clamped,
R810. ....... .0.............I. ....... . ........ 146
Thickness Variation... ..... ..... ..... .... ...... 148
Critical Load vs. RT, Square Plate, Simply—
SupportedOOOOOOO0..............I..O0.00.00 ..... 152
Critical Load vs. RT, Square Plate, Clamped.... 153
Central Deflection for Different R Values,
8-8 and Clampedooooooooo ooooooo ooo ooooooooooooo 155
Plan.......0... ....... OOOOOOOOOOOOOI0.0000 ...... 156
3.62 Node Arrangement, Square Plate, h = a/8 ........ 157
3.63 Contours of Membrane Force, Principal Stress
and U-displacement, Undeflected Plate, RT = 1.. 160
3.64 Contours of Membrane Force, Principal Stress
and U-displacement, Undeflected Plate, RT=1/4... 162
3.65 Contours of Membrane Force, Principal Stress
and U-displacement, Undeflected Plate, RT = 2... 163
3.66 Convergence of Buckling Solution. ............... 165
3.67 Convergence of Buckling Solution ....... . ........ 166
3.68 Plan ................................... . ........ 168
3.69 Variation of Critical Displacement Versus RT.... 170
3.70 Variation of Critical Displacement Versus RT.... 171
3.71 Buckling Modes of Simply-Supported Square Plate,
RT!l..................... ....... . 0000000000000 174
3.72 Buckling Modes of Clamped, Square Plate, RT = 1.. 174
3.73 Convergence of Center Deflection, s-s, Square
Plate... ............................. .. ........ . 176
3.74 Plots of Stress Components, Square Plate,
S-S, RT=loooo ooooooo ooooooooooo ooooooo o ooooooo 178
3.75 In-Plane Displacement, Square Plate, s—s,
RT:1..0............O........... ............... 179
3.76 Square Plate Under Uniform Lateral Load......... 180
3.77 Central Deflection Versus Edge Displacement for
Different RT Values, Simply-Supported Square
Plate........ ...... ........O.............. ..... O 182
3.78 Plots of Stress Components, Square Plate,
Simply-Supported, RT = 1/4................. ..... 183
3.79 Plots of Stress Components, Square Plate, Simply-
Supported, RT = 2. ...... ................ ........ 185
IL80 Principal Stress Versus Edge Displacement for
Different RT Values, Simply-Supported Square
Plate...... .............................. .. ..... 187
In
II)
3.81
3.82
3.91
A3
A4
BI
32
InéPlane Displacement, Square Plate, Simply-
Supported, RT = 1/4 ............................. 188
In-Plane Displacement, Square Plate, Simply-
Supported, RT = 2 ............................... 188
Convergence of Solution ......................... 190
Central Deflection Versus Edge Displacement for
Different RT Values, Clamped Edges .............. 190
Plots of Stress Components, Square Plate, Clamped,
RT = 1 .......................................... 192
In-Plane Displacement, Square Plates, Clamped,
RT = 1.... ..... . ................................ 193
Plots of Stress Components, Square Plate, Clamped,
RT = 1/4 ........................................ 196
Plots of Stress Components, Square Plate, Clamped,
RT = 2 ........................................... 197
Max. Principal Stress Versus Edge Displacement
for Different RT Values, Square Plate, Clamped... 199
In—Plane Displacement, Square Plate, Clamped,
RT = 1/4 ........................................ 201
In-Plane Displacement, Square Plate, Clamped,
RT 8 2 ..... . .................................... 201
Difference Operator for Left Hand Side of Equi-
librium equation (2.24), a = 1.................. 215
Difference Operator for Left Hand Side of Equation
(2.14), a - 1, Uniform Stiffness P1ate.......... 216
VFinite Difference Operator for Left Hand Side of
Compatibility Equation (2.15), a = l............ 217
Finite Difference Approximation of Equation (2.29),
a.l.........l...0...............0.0.......0.0. 218
Nbde Arrangement for Force Boundary Condition, Square
Plate,h-8/8.................................. 220
Nede Arrangement for Force Boundary Condition, Square
Plate,h.all-6.00.00...000......000.0.0.0.0.000 221
xi
0
B3 Node Arrangement for Displacement Boundary
condition,ha8/12....0.0.....0.0.0.......O.I.0. 222
B4 Nbde Arrangement for Displacement Boundary
Conditon, h = a/l6............................... 223
C1 Nede Numbers..................................... 227
C2 Arrangement of Vector K......................... 228
C3 Arrangement of Nodes in Vector KK............... 231
xii
CHAPTER I
INTRODUCTION
1.1 GENERAL REMARKS
The widespread use of plate elements in many engineering
structures such as buildings, bridges, pavements, missiles, containers,
ship structures and space structures has made plate analysis the
subject of scientific investigation for more than 200 years. Because
of their two dimensional action, the mechanical behavior of plates
under thrust loads is completely different from beam elements. In
contrast to beam elements, in which buckling is usually associated
with collapse of the structure, the buckling of a plate is not an
end point in the serviceability of the structure.
The capability of a plate to carry load after buckling
is an interesting subject which has motivated many investigators
to study posbuckling behavior of plates, especially in connection
with weight-sensitive space applications. Most of the plate analyses
involve-uniform stiffness plates. However, elastic plates of variable
stiffness are used in many engineering structures such as aircraft
wings, turbine disks, etc. The need to conserve material and/or
minimize weight motivates the designers to make optimum.use of
the material.
‘..
,.:
we
From the structural point of view, knowledge of critical
buckling loads is of great importance. To make an optimum design
with respect to some variables, an extensive analysis of the variable-
stiffness plate is necessary. The failure strength of a thin plate
can exceed the buckling strength appreciably. In many cases, the
structure is not sensitive to large deflection. Thus, it is of
technical importance to consider the postbuckling behavior of plates
(especially the variable stiffness plate) in order to optimize
the design.
Although a considerable amount of work has been done in
the area of variable stiffness plates, most studies have achieved
solutions by analytical methods which are restricted to some specific
geometry and boundary conditions. (See Section 1.2)
The purpose herein is to investigate the behavior of a'variable
stiffness plate so that a full history of the stress and strain
components of plates with different stiffness variation can be
presented. Such a history will help give the designer a better
understanding of the behavior of the plate and the effect of stiff-
ness variation on various aspects of the problem so that a more-
nearly optimum design may be achieved. Two different types of
variation in stiffness are possible: one with uniform thickness and
varying E, such as reinforced concrete or fiber-reinforced plastic.
The other has variable thickness and constant E. Both cases are
considered and analyzed.
4
u~§l
t-uos
5‘...
.5, ‘
n“...
L:
.-
1.2 PREVIOUS DEVELOPMENTS
The study of plate theory began in the l7601s. Euler (17)
presented the first mathematical approach to plate studies in
1766.
In 1815, Sophie Germain (22) presented a fairly satisfactory
fundamental equation for the flexural vibrations to the French
Institut as the result of her investigation during the 1809 to 1815
period. Within the same period (in 1811), Lagrange arrived at his
equation, which is known, therefore, as Lagrange's equation for the
flexure and the vibration of plates. Kirchhoff (1824-1887) is
considered the founder of the extended plate theory which takes into
account combined bending and stretching. In 1910, Von Karman
introduced a set of differential equations valid for plates subject
to large deflection. These equations are referred to in the liter-
ature as the large deflection equations.-
The development of the modern aircraft industry directed
the attention of many scientists and researchers toward the study
of plate vibration, plates subject to in-plane loads and postbuckling
behavior of plates. The earliest solution of a flat plate stability
problem apparently was given by Bryan (lO)in 1891.
The ability of a plate to carry additional load after
buckling was apparently discovered in the late 1920's through
experimental studies made in connection with the design of air-
planes. In 1929, wagner (49) studied a shear web and based on his
findings, established a criterion for postbuckling strength of
the web. In 1942, Levy (29) presented solutions to the plates with
l.2 PREVIOUS DEVELOPMENTS
The study of plate theory began in the 176013. Euler (17)
presented the first mathematical approach to plate studies in
1766.
In 1815, Sophie Germain (22) presented a fairly satisfactory
fundamental equation for the flexural vibrations to the French
Institut as the result of her investigation during the 1809 to 1815
period. Within the same period (in 1811), Lagrange arrived at his
equation, which is known, therefore, as Lagrange's equation for the
flexure and the vibration of plates. Kirchhoff (1824-1887) is
considered the founder of the extended plate theory which takes into
account combined bending and stretching. In 1910, Von Karman
introduced a set of differential equations valid for plates subject
to large deflection. These equations are referred to in the liter-
ature as the large deflection equations.-
The development of the modern aircraft industry directed
the attention of many scientists and researchers toward the study
of plate vibration, plates subject to in-plane loads and postbuckling
behavior of plates. The earliest solution of a flat plate stability
problem apparently was given by Bryan (lO)in 1891.
The ability of a plate to carry additional load after
buckling was apparently discovered in the late 1920's through
experimental studies made in connection with the design of air-
planes. In 1929, wagner (49) studied a shear web and based on his
findings, established a criterion for postbuckling strength of
the web. In 1942, Levy (29) Presented solutions to the plates with
_———_-—‘-
large deflections under combined edge compression and lateral
loading. His investigation was based on analytical solutions using
Fourier series. He also considered postbuckling analysis of plates.
In 1970, Supple (42) analyzed a rectangular plate with constant
in-plane compressive loads on opposite edges using the out-of—plane
deflection, w, and the Airy stress function as variables.
A considerable amount of work has been done on plate
analysis by different methods; among these are solutions of the
equilibrium equations by series expansion, energy methods , and Vlasov's
method (46). Until relatively recent times, however, the investigations
have centered on analytical solutions, which are in most cases limited
to relatively simple geometry, load, and boundary conditions.
In many particular cases where these conditions are more complex,
the analysis via the classical route becomes increasingly difficult
and is often impossible. In such cases, the use of an approximate
approach becomes more practical due to the flexibility and quick
results.
Near the end of World War II, the invention of digital
computers,with their capability of processing large numerical
problems, caused rapid development of various numerical techniques.
Of these, the finite element, finite difference and boundary integral
methods are of most general use.
Although previous analysis of variable stiffness plates
has been limited, there has been a considerable amount of work
done on uniform stiffness plates using finite element techniques,
and dealing with stability and postbuckling of plates. The finite
element method waSintroduced by Turner, Clough, Martin, and Topp
(45) in 1956. Argyris (A) and Zienkiewicz (53) have made numerous
contributions in this field. Gallagher (20) and Hartz (24) also
have made great contributions in improving the method and including
nonlinear terms. A series of studies considering postbuckling
behavior of plates was made in the 1970's (13, 19, 52) using the
finite element technique. Murray and Wilson (31) have conducted
research on postbuckling of plates,considering various aspect
ratios and applying the finite element method. Other significant
contributions include papers by Conner (l4) and Yang (52). They
applied the finite elementnethods to solve postbuckling plate problems.
The finite difference method is also one of the general
numerical solution techniques which has been frequently used. The
finite difference method was first used by N. J. Neilsen (33) for
analysis of plates in 1920.
The first finite difference solution of the large deflection
of plates is due to Kaiser (26). More recently, Basu and Chapman
(5) contributed to this study. Kaiser also carried out some ex-
perimental tests which verified the theoretical results. Both
aforementioned investigations were formulated in terms of lateral
deflection, w, and a stress function. The finite difference ap-
proach to large deflection of plates was also used by Brown and
Harvey(9) who have studied large deflection of plates subject to
lateral pressure combined with different ranges of edge loadings.
More recently a new method of solving finite difference equations—
namely,the dynamic relaxation method was described by Otter (35). The
————_—l
basis of the method is to add dynamic terms to the equations. The
addition of dynamic terms such as acceleration and viscous damping
makes the problem analogous to a vibration problem. The damping
coefficients are taken corresponding to critical damping resulting
in a motion which dies out quickly. Thus, the solution to the
static problem is obtained. Rushton (38,39,40) has published papers
applying the dynamic relaxation method to large deflection of plates
subject to lateral load and to postbuckling of plates under in-plane
loads. Since the method is an alternative technique to solution
of the finite difference equations, it has the advantage that
variations in stiffness can be included. Rushton has stated that,
with appropriate time increment and damping coefficient, a solution
can be obtained with no difficulties.
The Boundary integral equation (BIE) method has also proved
to be successful in solving plate bending problems. Jaswon (25)
and Haiti (30) introduced the direct method of solution and recently
.Altiero and Sikarskie (3) presented the indirect method of solution
fifliich proved to be more efficient. In 1980, Wu (50) modified the
“method by moving the integration contours outside the real boundaries.
Tile plate of interest is embedded in a fictitious plate for which
Ifle Green's function is readily known. Fictitious forces and moments
ire then applied outside the real boundary and the solution can be
>btained by finding the magnitude of these fictitious loads such
111st the original boundary conditions are satisfied. The method
r213 proved to be very efficient in general plate bending problems.
Particularly pertinent to this study is a paper by Prabhakara
.). In his paper, he considered postbuckling of orthotropic
xtes. Recently (in 1980), Kennedy and Prabhakara (27) have
1died the postbuckling behavior of orthotropic skew plates and
:ained solutions to some problems using a series expansion method.
L.3 PRESENT INVESTIGATION
In the study of thin plates subject to lateral and edge
loading, especially in the postbuckling range where the deflections
are not small, the Kirchhoff theory (which neglects stretching and
shearing in the middle surface) can not yield satisfactory results.
In this case the Von Karman large deflection equation can be employed
to obtain more accurate results.
In Chapter II, a brief review of the theoretical background
is given and the derivation of the compatibility and equilibrium
equations is first presented. Next, by applying the ordinary finite
difference method, the required operators are derived and the pro-
cedures for solution of different problems are briefly discussed.
Two different alternative methods of formulation are considered:
a) in terms of lateral displacement, w, and a stress function,
b) in terms of the displacement components, u, v, and w.
Yhe solution procedures for both methods are also discussed. A few
imamples of practical boundary conditions are listed and theoretical
ialations for each boundary condition are mentioned.
Chapter III includes numerical solutions and analysis
f the results. A computer program and the required subroutines
iseee computer program in Appendix C) have been developed to facilitate
P1>lication of procedures discussed in Chapter II.
In order to provide a more complete view of the variable
tiliffness plate and it's behavior relative to the uniform stiff-
3538 plate, several different types of variation in stiffness are
)rlsidered. Uniform stiffness plate results are given for comparison.
_-——!QM
For clarity, in the procedure presented, the results are always
accompanied by graphical illustrations of membrane and bending stress
as well as displacement components. The behavior of those graphs
and their relations with applied load is discussed.
The accuracy of the solutions is evaluated by comparison
of the results obtained with results from past studies and exact
results, where these results are available.
Convergence of the solutions is examined by using different
umsh sizes with extrapolation.
Results obtained for the effect of stiffness variation on
in-plane forces, bending moments, in-plane displacements and lateral
deflection, provide a good source of information for optimization
in each case. Although the optimization procedure is straight-forward,
the stability optimization with respect to amount of material used
is presented as an example. Two computer programs are provided,
One for force boundary conditions and the other for displacement
‘50undary conditions. Both programs are listed in appendix (C).
It was found that convergence was easily obtained for the range of
lxaading less than the second critical load because the assumed
Single-wave buckled shape is the only possible pattern of stable
equilibrium other than the flat plate. For loading beyond the
E3€=cond critical load, due to different possible equilibrium states,
time problem does not converge easily. For solution beyond that
range a proper deflection shape must be enforced, as appropriate
:EKDr the physical conditions of the problem.
10
NOTATIONS
The symbols are properly identified when first introduced;
for the reader's convenience, symbols are tabulated here.
Side length of square plate
a
c1,c2,.. Constants
Et3
D =-———-—7f- Flexural rigidity of the plate
12(1-v )
D Stiffness of a uniform stiffness, unit thickness .}
0 plate
Dr Reference stiffness = stiffness at center of the plate
E Young's modulus
F Force function
h,k' Mesh intervals in x and y
K = Et Membrane rigidity
Kr Reference membrane stiffness = Et at center of the
plate
K3 Membrane rigidity of a uniform stiffness, unit
thickness plate
*1 ,M.,M Bending and twisting moments per unit width
1‘ y xy
of plate
“ 2
(D1 - M ; M ) = (Dat )(Mx; M.; M ) Dimensionless moments per
1 y unit width
be Applied edge force per unit width of plate
In-plane stress resultants per unit width of plate
2
5e
(1W - N*; N* ) = (%—)(N ; N ; N ) Dimensionless membrane forces
0 x y xy per unit width of plate
(IV - ". " = . .
N , N ) (Nx’ Ny’ ny)/N Membrane force ratios
11
(N;; N'; N; ) = (Nx; Ny; ny)/Ncr Membrane force ratios
Y
q Lateral distributed load per unit area of plate
6': qaa/Doti Dimensionless lateral load per unit width of plate
Q Transverse shear per unit width of plate
_ edge stiffness .
- central stiffness Stiffness ratio
_ edge thickness . .
- central thickness Thickness ratio
t Plate thickness
ti Unit thickness
T Temperature
15v,w Displacement components in x,y, and 2 directions
3 Volume
uo Edge displacement
(U; V) = (u; v) a/t2 Dimensionless displacement
i .
U a
U/Ucr
9 KO
U a ”if Dimensionless displacement
U0 Dimensionless edge displacement
[1* = U
cr
W a 55—— Dimensionless lateral deflection
i
":37,z Cartesian coordinates
(I(; Y) = (x; y)/a Dimensionless coordinates
0'- ==£~ Grid size ratio
81 = 51- Membrane stiffness ratio
Kr
oflr Coefficient of thermal expansion
6 Di
1 § -D—— Flexural stiffness ratio
12
.1,A2,... Eigenvalues
Ll,A2,... Eigenvectors
,s ,8 Strain components
{ Y KY
,0 Components of normal stress
i Y
2
'. ' = a . o 0
3x, 0y) (D t.)(ox’ 0y) Dimen51onless stress components
,T ,... Shear stress components
Ky xz
= t@' Stress function
Airy's stress function
= cp/Na2 Dimensionless stress function
i] Coefficient matrix for m
1w], [bw] Coefficient matrices for w
:Aw]; [Bw]) - ([aw]; [bw])xh4 Coefficient matrices for w
.u], [Bu] Coefficient matrices for u
v], [Bv] Coefficient matrices for v
ul], [AuZ] Coefficient matrices for u
V1], [Av2] Coefficient matrices for v
CHAPTER II
THEORETICAL DERIVATIONS
2.1 General
In this chapter, the equilibrium and compatibility equations
of the plate based on the theory of elasticity are first derived. Then,
the finite difference approximationsto these equations are developed.
These will be used to facilitate numerical solutions of those equa-
tions,for which,in most of the cases, closed form solutions,if not
impossible,are very tedious.
Thin plate theory is applied and homogeneous, isotropic
material is assumed.
Depending on the boundary conditions, two different approaches
are possible. Here, both approaches will be diScussed.
Geometrical and material nonlinearity can arise in plate
Prxiblems. In this study only geometrical nonlinearity will be con-
Sid ered.
Figure (2.1) shows the geometry and orientation of a plate
it: the cartesian coordinate system. The x-y plane lies in the middle
Filaine of the plate and z is normal to the middle plane.
Internal forces and moments acting on the edges of a dx by dy
pl«ate element, as shown in Figure (2.2), are related to the internal
Stresses by the equations:
13
14
t/Z t/z
Nx = f oxdz N = f 0 dz
-t/2 y ’t/z y
t/ t/
2 2
N 8 f T dz N = f I dz (2.1)
t/2 t/z
Q =f r dz Q =[ r dz
x -t/2 XZ y -t/2 yz
’ t/z t/z
Mx = f oxzdz My = f o zdz
't/Z -t/2
t/2 t/
2
M =f ‘I zdz M =f r zdz
xy -t/2 xy yx -t/2 yx
where Nx’ N , N , Nyx = in-plane normal and shearing stress resultants.
Qx’ Qy = transverse shearing stress resultants.
M , M = bending moments.
X Y
M , M = twisting moments.
KY YX
2.1.1 NONLINEAR EQUILIBRIUM EQUATIONS
In the literature, nonlinear behavior is commonly classified
as either
1) Material nonlinearity
2) Geometric nonlinearity
Material nonlinearity may arise in case of time-dependent
lfilteria]. or materials with nonlinear stress-strain relations (plastic,
tlastoplastic, viscoe'lastic, etc.) .
Geometric nonlinearity is usually associated with large
iJBplacements. It may also occur for small displacement if the
15
Figure 2.1 Rectangular flat plate
Figure 2.2 Plate element dxdy in undeformed configuration
16
behavior is such that variation in the applied load alters the
distributions of displacement.
In this report only geometric nonlinearity is considered
and the material is assumed linear elastic, isotropic,
To determine equilibrium equations applicable to moderately
large deformations, they must be derived using slightly deformed
configurations. Figure (2.3) represents stress resultants and
internal moments for an element dx by dy in the deformed con-
figuration. B and B are rotations in the xz
x 3’ 3N
x
-—-dx etc.
3x
and yz planes
respectively, and N: denotes Nx +
Summation of forces in the x-direction gives:
8N
__§
-Nxdy + (Nx + 3x
3N
_ +.__Z§_ a
dx)dy Nyxdx + (Nyx 3y dy)dx 0
which simplifies to
aux 3N x
—3—}—{ +—X-ay =0 (2.2)
Similarly, summation of forces in the y-direction leads to:
EN EN
is); +—]—3y =0. (2.3)
Fromsummation of forces in the z-direction, we obtain
3Q 3Q 2 2 2
- x-—1=q+N3—1"—-+N a"’+2N3—"—- (2.4)
EX 3y x ax2 x ay2 Ky 3x3y
SLnnmation of moments about
x and y axes will result in:
and homogeneous.
17
3M 3M
Qy 3 3y + 3x
BMX 6M x
Q = + —X— (205)
x 3x 3y
Figure 2.3. Schematic illustration of internal forces and moments
on the element of middle surface in deformed configuration.
18
2.1.2 Relation between stress resultants and displacements.
From Hooke's
where:
For moderately large
law, we have
N = (e + v e )
x l-v2 x y
Ny — Et (5 + v ex)
1-v
xv YX2(l+v) xy
_ Bu 22.2
8x - 3x + l/2(8x)
3v 3w 2
=-——+12_
6y 8y / (3y )
_ Bu 3v 3w 3w
Y .___ .__. .___
xy ‘ ay 3x 3x 3y
displacements, the relations between moments
and lateral displacement are:
32w 32
M = -D(—— +v—3)
x 2 2
9x 8y
2
M = 4,49% + 1%)
y 3y 3x
M = M =-1)(1-\»)32w
xy yx 3x3y
Et3
where D = —-—-—2- is the flexural rigidity of the plate.
12(l-v )
(2.6)
(2.7)
(2.8)
19
2.2 Formulationimxterms of stress function and w
By substitution of Equation (2.8) into (2.5) and (2.4),
we obtain the equilibrium equations in the z-direction in terms of
membrane resultants and lateral displacement w:
2 2 32D 32w 32D 32w 32D 32w
V (W w) ‘ (1“) ‘2—3 ‘2 3x3 3x3 2 2
8x 3y y y 3y 3x
2 2 2
_ 8 w 3 w 3 w
3x By
The compatibility equation for mid—plane strains is:
2 2 ' 2
3 Ex 3 e a Y 32w 2 82W 32W
8y2 + 2 - Bxay = (axay) - 2 2 (2'10)
3x 3x By
and from (2.6), theastrainsin terms of membrane forces are
0)
ll
1
Et (Nx-vNy)
x
e = JL—-(N -vN ) ' (2 11)
y Et . y x '
- 5.131.
ny Et ny
Now, we define a stress function, Q, similar to Airy's
stress function, so that:
2
N .152
x 2
8y
32
N -——‘§ (2.12)
y 3x
2
N -12..
\
l
I
20
The Airy's stress functions is defined as
2.52.
O a
x 2
3y
32 '
0y = ——‘§L (2.13)
8x
- - 239:
Txy 3x3y
are N = t ox = t-BJEE etc.
X By
is is suitable for a uniform thickness plate. However, in the case of
:iable thickness, if we define m', as (2.13), it will complicate
a formulation. For example, substitution in equilibrium equations
.2), would result in
‘33—:32 2 J" 2 3:3 ‘ta 2=0
3y axay y axay axay
[ch in the case of uniform thickness, leads to 0 = 0.
For our purpose the definitions of (2.12) will be used.
)stitution of (2.12) into (2.9), will result in the equilibrium
Jation,in terms of m and w, as
2 2 32D 32w 32D 82w 82D 32w
V (nv ‘0 - (I'V) 7 7 "'2 axay axay + —-2' "-2
3x 3y 3y 8x
_ 329 32w+ azgz 3w 322 32w
- Q+ .__2+ 2 3y2 2 Bxay axay ° (2'14)
3y 3x 8x
substitution of (2.12) into (2.11) and then into (2.10), we
tain the compatibility equation in terms of m and w:
21
2 2
3 a 2(l+v)
ayz [t yy xx aXZE 3X y fipxy
Z 2 2
3x3y2 3 2 3y2
32
where cp .- __52 , etc.
xx 3x2
2;; Formulation in terms of 3 displacements u,v,gand w
Substituting equation (2.7) into (2.6) and then into (2.2)
and (2.3),along with substitution of (2.8) into (2.5) and then into
(2.4),results in 3 equilibrium equations in terms of u,v, and w.
The equation of equilibrium in the x. direction (2.2)
becomes
2 2
Et 1:... 32V wit-k v(32v + 32w 3):) l-v Et 3 u 8 v
2
l-v2 3x2 3x2 3:: My axay ay 2 1w 3y2 3an
2 2
+§ifl+flu +_}___9EE)_ _3_1_1_+_];(_3_w >24. \)g_( E)... %(_:_W_)2
axay 3y ax 3y2 1_\)2 ax ax 2 3x
+ 1"“2 30%) fl+fl+flflJ .. 0. (2.16)
2(1_\, ) 3y 3y 3)! 3}! 3y
Similarly, the equation of equilibrium in the y-direction (2.3),
becomes
Et 32v 32w aw azu 32w aw l-v Et 32v azu
2 2 + 2 a—y ”(any + axay E) 2 2 2 + axay
1"v 33' 3y l-v ax
82w aw aw 32w 1 3(Et) 3v 1 3w 2
+axa E+T—2’+ 2 7+2?) +
y 3,, 1_\, 3y 5* Y
Bu v 3w 2 l-v 3(Et) 8V 311 3W W
+—-— -— —— -— —— I
v5 2(3x) :] +2(1-v2) 3x [3x + 3y + 3:: 3y] 0 (2.17)
22
Interchanging u and x with v and y respectively in
equation (2.17) results in equation (2.16).
The equilibrium equation in the z-direction (out of plane)
can be expressed in terms of 3 displacements by substituting (2.7)
into (2.6) and then the results into (2.9). We obtain:
2 2 2 2 2
2
2 2 3 D 8 w 3 D 8 w 8 D 8 w
V(DVw)-(l—v)—--—-- +——-
3x2 ay2 axay 3x3y 3y2 3x2
2
Et Bu 13w2 av v3w2 3w Et 3v 13w2
=q+—— —-+—(—-—) +v—+—(-— ]—+—[—+-(—)
1_\’2 [3x 2 8x 3y 2 3y 3x2 1-v2 8y 2 3y
+ Bu +3(__3_3)2 32w + Et(l-v)fl1_+_ay_+_3l§3 32w (2 18)
3x 2 3x 3y2 l-v2 3y 3x 3x 3y axay °
NOte: Since in this approach we are working with displacements,
compatibility need not be checked.
.L-é FINITE DIFFERENCE APPROXIMATION
So far we have derived the necessary equations for analysis
of tflle plate, but solving these coupled nonlinear partial differential
equations analytically may be difficult.
Here we employ finite difference techniques to transform
the differential equations into ordinary algebraic equations in terms
of values of the functions themselves at certain specified points.
'ggizl;‘flgrinciple of finite differences
The derivation of finite difference expressions is based on a
Tayl‘DIT series expansion. We expand the function at some
successive grid points, truncate higher order terms, and
. 801.
"ee for desired derivatives, we can obtain approximate expressions
'23
for first, second, or higher order derivatives in terms of values
of function at the discrete points.
Y
f’,. f(x)
m;2 m-l m m+1 m+2
AL»
Figure 2.4. Function f(x).
19 one dimensional cases we obtain the following approximation of
the derivatives of the function:
. ._1_ - - l 2 m
f'(x)m 2(Ax) (f1n+1 fm—l) 6(Ax) fm +....
n _1__ _ _ l— 2 .V
f (x)m a 2 (Em+1 2£m + fm_l) 12(Ax) f; + (2.19)
(AX)
' ._l;___ - _ _ l_ 2 v
f"(x)m 3 (fm+2 2 fIn+l + me-l fm-Z) 4(Ax) fm +...
2(Ax)
iv
1
f(x)In - 837; (fm+2 -4 fm+l + 6 fm - 4 fm_1 + fm-Z)
l 2 vi
-g(Ax) fm + ....
24
In practice, we truncate the terms following the parentheses.
These represent the error in the approximation. We will refer to these
error terms later in the discussion of accuracy.
x
NN
O
I
l
'N
Y NW:|3— ——O———ONE
I, I
I l l
| I
wwo———‘"O———Oo——-(l)§-— 05:
l I I
I l
I ., !
Sw3——_¢_)_..._ SE
l
|
555 b
Figure 2.5
To determine the finite difference approximation for two dimensional
problems, we consider Figure (2.5), and the fact that similar re-
lations for the approximations to the derivatives in the X-direction
hold also in y-direction.
Thus,we will be able to derive expressions for any order
of derivatives in x and y and combinations of x and y derivatives.
For example, if we consider grid points of Figure (2.5) ,
the derivatives with respect to x and y at point 0 are
l
fx- 2h (fE-f
)
W
fy g'2k (f8 - fN)
f = i— (f -2f + f )
xx 2 E 0 W
h
fyy = 7 (fS-ZfO-l- fN)
h
f _ 1
xy ' 4hk (fSE fsw' fNE+ wa)
or if £ = a then
fxy = F (fSE-fSW-fNE + wa) etc.
Often, these formulas for derivatives are represented
geometrically by stencil patterns, such as in Figure 2.6.
--—®
igure 2.6. Two dimensional operator for fx
(9-- -®- -j®
772- CD" (ID-- (:9
C9-- {'9
(2.20)
26
2.4.2 FINITE DIFFERENCE APPROXIMATION OF METHOD DISCUSSED IN SECTION(2.2)
In section (2.2), we derived equations of equilibrium and
compatibility as well as the components of internal forces and dis-
placements, in terms of lateral displacement, w, and stress function,
Q. In this section, we will discuss numerical solution of those
equations using finite difference techniques.
To find a solution to a plate problem, we must satisfy
both equilibrium in the z direction, and compatibility.
a) Equilibrium Ermation I
For this purpose, we will apply relations (2.20) to equilibrium
equation (2.14), and the results will be represented in two dimensional
Operator form. Introducing
Dr *3 flexural stiffness of the plate at some reference point
(center of the plate in this case),
U
- __i_
61 ‘ D
r
A finite difference operator for the left hand side of the equilibrium
6.and6
c
e‘l'aation (2.14) is given in Figure 2.8 where Ga, 6,), (1
refer to midpoints a, b, c and d of Figure (2.7) as explained in
reference ( 8).
27
Figure 2.7
Note: All terms in the operator of Figure (2.8) are coefficients of w.
If the grid spacing is the same in the x and y directions (a = l),
oPet'ator (2.8) will be simplified to the operator given in appendix
(A. 1) .
In case of a uniform stiffness plate, where 6, =- D— = 1
for all points, the operator reduces to the usual finite difference
OI’erator for V4w, as given in appendix (A.2).
‘
28
.A¢H.~v coauaauu assessaaauo ac cuss can; some now “compose ouaououuan .o.~ masses
oauuu Paoonqom I > . MI 5
u
6.33 coo—6:: I lab I «a
a m
o 6
u e
unwound: mucououou I a
d
:aisai 2;? 32-5.. 2353....
d
a are N was $3.3 $33 $8.. necréfiua
u 6
2.? in? 3 €ud~a~+
v a I u A I a
2 2. 3? dues. z m 3 u 2 0+ 3? Gag. NI
u . S :7 . 2+ 3.3+ 3.? . A: .a: 3:: .. an
A a+ avau +H Auap~+sa+~v as 9+ N u
Seize: Iain??? 222333
6 I z o I 0 I
a? 33% A 0+ 21533 $3 2., same A
20 a
,oc
.o-o
29
For the right hand side, if cp values are known, we can
derive the finite difference operator as coefficients of w, as
shown in Figure (2 .9 .a) , where
32 .. (FE-2'90 + W
<9 = —%
xx 3x h2
2
fie .. (phi-2'90 + CPS _ 0‘ (IN-2'90 + c"3) (2.21)
‘Pyy ' ‘ 2 ‘ 2
3y k h
a 3259 .. CPSE + (PNw'cpsw'cPNE g c‘(CPSE + ‘PNw'q’sw'q'Nfl
.2ny axay 4hk 2
4h
For a = l, we get the expression shown in Figure 2.9(b).
flLUTION OF THE EQUILIBRIUM EQUATION
To solve the equilibrium equation for w, we must have either
t _
he cp values, or the in plane forces (cpxx, cpy and cpxy).
Y
Substitution of these values into the operator of Figure
2.9 (a) or (b) will result in a known Operator at each node.
By applying the operator at each node, we will be able to
5011:: a matrix of coefficents of w, which along with the qi vector
will form the right hand side of equation (2.14) as
{q} + wa] {w}
where {q} is the lateral load vector and [bw] is the coefficient matrix
containing constants. Similarly, application of the operator
of Figure (2.8) will result in the formation of matrix DranHw} in
the left hand side of equation (2.14). [aw] is a constant coefficient
matrix.
Therefore, we can represent the equilibrium equation in
numerical form as :
30
_ _ ._ 2 _ _
' I ‘ny °‘ ‘9... I + 9‘2" ‘ny
I I I
I I I
I I -I
l ("W L' " 'I -mpyy + 0‘ 'pxx) I' — " (PYY
2
h I I I
I I l
4 I I
2
+%'ny I-“ 0' cpxx ___ -%CPXY
(a)
1 ' _ __ _ _1 -
’FW’SEWNw'q’NE‘IIIF '7 ‘PE'Z‘I’oJ'q’w EWSEWNW (ENE-(98W)
I I I
I I I
I I I
-2(
2 ML; -v a—%>1 (2.27)
3x2 3y2 3y (Et) y 3y 3x
+ zuwflagig 3(Et) 2 _ _1_ 3 203133329 }.
3y (Et)2 Et Bxay ax3y
Now by applying relations (2.20) and adding all contributions
at each node, we will obtain the finite difference operator of
Figure (2-10):
‘where
Kr: Et = in—plane rigidity of the plate at center.
K1 = in-plane rigidity at point i.
If a --£-- 1, this operator will be simplified to the one
given in Appendix (A.03)
An alternate approximation to the compatibility equation
(Z-GIS) in finite difference form can be developed.
Denoting
1
Fl - K ((Pyy “V (9“)
F=-1-(
+AVIA2 I A2IA2 2m.V.2I 12 2+. 2AI 22V ? A22I22I22¢22V 35% +
. A . .2
A....|.2.IM..AII.A....2.. . ..A 21:2,..- 2.... III. .- .2... .223? -
2I2 I2 NII. 2I2 .I 22II2INI.
A 2 2+22 2VN2 2A2A A 2 2VAM$AVA2+AA£AVA2 2.2 A 2+ 2 22 2VA2 + 2A..A
2 2 . z . 2 ,
2 2 . I2 . . I2VI2INIA2+2AI2V3+
AAI 28W» IA 2+ 2A .2VA3+ 2+ 2AI22V AaA? AaVA+ AA22I22VAIAIIA2I I AA A2,. A 2
22A 22+. 2 2 V A22IA2V .MIA +AA22+. 2AI22VI Iml + 2....
AA 2I 22V IA 1 2AI 2VI 22- IleflalIlVI IA22+.2AI2 2VA ASIAVA+ A 22I22 A
A22I22VAA£AV+AA£AV.22I A 2$A _A.2+2V.2 A22I22VAA2+AVIAA2+AV.2...I . .2
'M I.
AA. I22I22+A2VA>+AYAI+ AA22I22 2.2V .mIA 2+A2 2+. 2AI 222V 21 A22I22Iu2+a2VA>+AVAa
.2A A
I AiAVAI I I .45... +
HAzuI22VA3MI22VH_ NA=I 22VII IA: 2+ 22 222V 3+ HA22I2 22VAH2 I2 2VHA «a
A 22+ 22.22-22VI_A. I 2.A_.A A2 I22VA Aa+AV A_.IA AAo+AV _. .2..- A2 2-2 2+2 2- 28h. +.2A_.A
A N
NO
34
Equation (2.15) can be written as
2 2 2
a F 3 F a F 2 2 2
1 2 3 3 w 2 3 w 3 w
8y?" 3x2 8x8y axay 8x2 3y2
Now we can approximate derivatives of F1, F2, F3 as
a2Fl F1 -2F'l + F1N
2 S 0 etc. , where Fl , according to (2.27),can be
ayz k2 S
approximated as
cp-ZCP +4) (I? -2¢ +2)
[53 S 41-» SE 3 SW ,etc.
1
F = _—
K k2 h2
l
S S
For a = £ = 1, this approximation results in the operator
given in Appendix (AJI)
W TO THE COMPATIBILITY EQUATION
To be able to solve the compatibility equation, we must have
Values of w at nodal points. Then, we are able to compute the
right hand side at each node using expressions (2.26).
To determine the left hand side, we apply the operator of
Figure (2.10) or the one in Appendix (A.3) or (A.4) at each node.
By adding all contributions, a coefficient matrix will be
formed, Thus, we have:
[AJICPI = {w} (2.30)
where the w vector is known, and solutions of this system of
e
quations results in the cp values at prescribed nodes.
35
c) Approximation to other equations
After solving equilibrium and compatibility equations, we
may be interested in calculating in-plane forces and displacements
as well as bending stresses. Finite difference approximation to
some of these equations will be discussed below.
Insplane forces
By definition (2.12) we have:
- L22 2 28.20 + ‘PN
N
x ayZ k2
322 cI’E'ZCI’O + qu
N = 2 = 2 (2.31)
y 3x h
__ 32g) __. cPSI: + (PNW-‘%I3 ‘PNE
xy axay 4hk
infiplane Displacements
Comparing equations (2.7) and (2.11)
-22122
l
E:x - 3x 2(3x) 2 Et1+<—2¥>(——§-fi—>a< 2h>(2h )+
K -K 2 w 2w w -w
S N a S N 2 v E W 2
a(——§H—)D§‘G‘jfir§ +‘5I 2h ) J
Finally, the equilibrium equation in the y-direction can
{be represented in the following operator scheme:
(u-operator) u + (v-operator) v + yw function = 0 (2.38)
‘22~§Eugilibrium in z-direction
The equation of equilibrium in the out-of-plane direction
1
s 1Jltroduced in Equation (2.18).
41
_ 1+v
L‘1—§”-> K. — - -<-l-I-"-> Kc
I I I
I I I
I I I
_g_ l;v-(KS-KN) '- — O " ‘ V - — 202.222) — $3422.22
I I I
I I I
, I I I
K -
b) V-OPERRTOR
FIG 2.12 OPERFITORS FUR Y-EQUILIBRIUI’I EQ.(2.17)
P
42
Approximation to the left hand side of this equation has
already been explained in Section (2.3.2-a).
As for the right hand side, we approximate the derivatives
of displacements to get
-v 2 w -w
(RHS) = q + -—K—21II(—-u-§—huw) +-12— (32-21% + 2 (vfih N) + V—g‘ (3h “>21
<-5::9-+:-) + [III— S ——N-h> +£33- (%)2+ 23%;“) +3 <———- 23h ”)23
22(——$-::-w§9-:) + (I-II) [AIS—u: :1”) + (3%?) + 2(WE;:W)(WE:N)J
a(WSE+WN:;:sw’wNE)} (2.39)
Finally, equilibrium in the z-direction can be represented
[Aw]{w} = RHS (2.40)
where matrix [Aw] consists of contants.
<1) Solution_procedure using this method
In order to perform numerical analysis in computer programs,
we arrange the equations in suchaway that the equations can be
represented in matrix form.
To find a solution by this method, we must satisfy all three
equilibrium equations. These contain three unknowns, u , v, and w.
There is no simple technique providing a direct solution to these
coupled nonlinear equations; thus, we employ an iterative technique
to solve them.
If w is known (or assumed),the first two equations of
e
quilibrium will become two uncoupled equations in u and v
¥ __
43
and can be solved as follows.
Applying the u—operator corresponding to equilibrimn in x
at each node and adding the contributions, and doing the same to the
v-operator, we get a matrix representation. of equation (2.16) as
[AulJNxN{u}le + [Avl]NxN {VIle = -{xw Function}le (2.41)
Repeating the same procedure for y-equilibrium equation (2.17) we
obtain
I:Au2]NxN{u}le + [Av2]NxN{v}le = -{yw Function}le (2.42)
Where [Aul], [Avl], [Au2] and [Av2] are constant coefficient matrices.
Both equations are coupled in u and v, and each contains
N equations in 2N unknowns.
One way of approaching this problem is to try to solve the
equations by iteration until reaching a solution that satisfies
both equations.
An easier approach can be employed if we realize that,
although the equations are coupled in u and v, there are no
In"Lited terms containing both u and v.
Therefore, we can combine the two, to get
Aul Avl u xw Function
-------- ' ' ---------- (2.43)
Au2 Av2 v yw Function
J
which is not only more efficient in computer programing but leads
t
o a. unique solution for the u and v displacements at specified
11 .
Ode points,based on an assumed (or known) w.
44
To solve the z-equilibrium equation, we substitute values
of u, v and w, in the right hand side (2.39), and solve equation
(2.40) for new values of w. The iteration will continue until the
new values of u, v and w, are equal to or very close to old ones.
A practical use of this method, including the details, will be dem-
onstrated in chapter 3.
POST SOLUTION DETAILS
After a solution is found for the three displacements
11, v and w, any components of stress and strain can be computed.
Average strain
Average strain at each node can be found by the finite
difference approximation of equations (2.7).
“E'uw 1 WE'ww 2
5 3 2h +‘2'(2h )
X
vs"'N 22 wS'WN 2
6y “—23-— +7 (‘73—) (“4)
uS-uN V 'VW W -‘Ww Ink-“W
E + II E M 3 N)
ny=°‘(2h )+ 2h (2h 2h
EEEEEQEane Forces
Substitution of equations (2.7) into (2.6) and approximating
the derivatives by finite differences, will result in
u- w-w v-v zw-w
N 2 K Eu” 2.2112 L! m 8 N2
X l_v2[2h +2(2h)+w(2h)+2(2h)]
v-v 2w-w - w-w
N a_BK SN g_ 5 N2 “E“w 3 Ew2 (2.45)
y 122M 2h>+2 (22 > +II(———2h )+2<2h >1
III-u V-‘V W'W W‘W
N aRSI-v) S N EW Ew S N
20 2) [a(-——2h )+(2h )+a( 2h )( 2h )1
-\)
45
BENDING MOMENTS
The bending moment approximation is given in equations (2.35).
2. 5 BOUNDARY CONDITIONS
In the small deflection theory of.plates, we consider only
out of plane (or flexural) boundary conditions because the effect
of in-plane displacements on the boundary is negligible. However,
-they become the chief factor in large deflections behavior and in
the postbuckling range. Thus, we discuss in-plane boundary conditions
as well as out of plane conditions.
The flexural boundary conditions, as commonly discussed in
elementary plate theory, are:
a) Simply-supported boundary
w=0
32
M =—D (.__w_ +
x 2
8x
2
v 2%) = 0 (on boundaries parallel to y)
3y
b ) Fixed boundary :
w=0
3w (n, normal to the boundary)
_20
an
c) Free boundary :
d) Others, such as elastic support, or partially fixed support, etc.
46
For the in-plane boundary conditions, there is a variety of
possible combinations that may occur in the postbuckling range.
On each edge, either 11 or v displacements can either be unrestrained
or have some specified values; also, there could be restrictions on
derivatives of either u or v, or both.
In terms of in-plane boundary conditions, we can classify
the problem in three major types.
1) - Force boundarL conditions (i.e. in—plane forces are specified
on the boundary).
If applied forces Nx’ Ny and N are known, we can use
to choose values of the Q functions on the
A practical
relations (2.12)
boundary points so that they satisfy boundary conditions.
example of this nature is discussed in chaper 3.
2; Displacement boundary condition Some possible cases are:
u and/or v are specified on the boundary,in which case the
a)
vallies of displacements would be '
x
assrlgned to boundary points.
b) Edge remains straight and parallel to y. (u-displacement is
constant all along the x = 0 edge).
e) Edge remains straight with no shear force along the edge; in this
Case, from equations (2.7) and (2.6), we have along edge x 2 0,
47
.__E_t__ 22 a_v 3222 .
ny 2(l+v) [3y + 3x + 3x 3y] 0 (2°46)
On supported edges (fixed or hinged), %§;- 0, thus
.. A! 2 22 . (1)
3x 3y
If the x-edge (x = constant) is straight with u constant, then
Bu 8v
-53; = 0; therefore, equation (1) results in 3;-= 0 along the edge-
d) Other conditions include possible restrictions on u, v or
their derivatives which lead to particular relations between
displacements or their difference approximations. For
example,the edge can be subjected to thermal expansion
(see section 3.2) such that %§" ex - constant along edge y = 0.
1). - Mixed boundary conditions. i.e.,case 1 applies to part of the
boundary, while the rest of the boundary is defined by case 2.
The computer program developed can solve either case 1 or 2.
Therefore, to handle a problem with mixed boundary conditions, we
can solve the problem by trial and error, as follows.
i) Assign some fictitious displacement values to the points
at which forces are specified, and solve the problem as
one with displacement boundary conditions.
ii) Compute forces at the boundary points.
iii) Compare with actual forces at the points.
iv) Correct previous fictitious displacements in such a way
that the solution is improved.
v) Repeat steps (ii) to (iv) until the computed forces are
equal to or close enough to the actual ones.
48
2.5.1 SOME EXAMPLES OF PRACTICAL B.C.
subject
1.
2.
Following is a list of some practical examples of plates
to various loading and boundary conditions.
Window glass can undergo large-deflection under lateral
wind pressure; the out-of-plane boundary condition is in
most cases simply-supported or sometimes built-in. Either
case may be accompanied by:
a) in—plane displacement possible.
b) in-plane displacement restricted.
Plates on stringers forcing the plate edges to remain straight,
as in many ship and aircraft sections surrounded by stringers.
Mechanical and instrumental plate elements subject to tens
perature change will be subjected to tension or compression
on some or all edges due to temperature change in surrounding
elements. Various combinations of boundary conditions are
possible.
The webs of structural steel profiles used in construction can
be categorized as plates subject to in-plane shear and normal
forces along the edges.
6 SUMMARY
\—
We will summarize the theoretical formulations discussed
11‘ ‘211apter 2, and mention procedures of solving some problems.
.Among several types of problems which can be solved numerically
based on the finite difference approximations shown, and using the
computer program which has been written, are:
49
2.6.1 MEMBRANE SOLUTION
For flat plate with w = 0 everywhere, the equilibrium in the
z-direction (normal to the plane of the plate) is trivial; to
find a solution for in-plane resultants and displacements,
a) In the case of force boundary condition, we solve compat-
ibility equation (2.30) with the w vector equal to zero.
[AJICPI = 0 (2-47)
The solution results in Cp values at discrete nodes, which can be
used in equations (2.31), (2.32), and (2.33) to find in—plane forces
and displacements.
b) If the displacements are specified on the boundaries, we
solve equation (2.43). Considering w - 0 everywhere,
we have
...... <..- =2.» (2.48)
for which the solution results in displacement values at the nodes.
Application of equation (2.45) then leads to the membrane resultants.
2-\6.2 LATERAL LOADING
a) Force boundary conditions.
1) Small displacement.In this case, the (9 values and
in—plane resultants are known from (2.5.1), so we can
solve the equilibrium equation (2.25) to obtain w.
Then (2.33), (2.34) and (2.35), can be applied to
‘_—;_g____ A
50
compute in-plane resultants, displacements and bending
moment 3 .
ii) Large deflection.
Compatibility equation (2.30)-and equilibrium equation
(2.24) could be solved iteratively, and the force resul-
tants and displacements can be calculated as discussed in a).
b) Displacement boundary conditions
i) For small deflections, we can solve equation (2.48),
neglecting the effect of w on in—plane solutions; then,
solve the z-equilibrium equation (2.40) for w, by ignoring
wbterms in the RHS.
ii) Large displacement problem - this requires an iterative
solution of equations (2.40) and (2.43) as discussed
in 2.3.3 (d).
L623 STABILITY ANALYSIS
a) Force boundary conditions.
The in-plane forces and m values are known from part
2'5-1 -a; then,we can use the equilibrium equation (2.25). If
q " 0, this will result in a characteristic matrix, the eigenvalues
of Which lead to the critical forces and the eigenvectors represent
the buckling modes.
1)) Displacement boundary condition.
Using in-plane resultants and displacements obtained from
(2""5-1.b) and forming the R.H.S. as a coefficient matrix for w,
with _ 0
q ~ , results in the characteristic matrix equation
51
{[Aw] - [Bw]} {w} = 0, (2-49)
for which the eigenvalues and eigenvectors lead to critical
boundary displacements and the mode shapes, respectively.
2.6.4 POSTBUCKLING
Since after buckling the plate takes on a state of stable
equilibrium, we can analyze the plate as a regular large deflection
C388 .
a) Force boundary condition.
We solve the equilibrium equation (2.24) and the
compatibility equation (2.30) iteratively.
b) Displacement boundary condtions.
In this case we employ the iteration technique to solve
the z-equilibrium equation (2.40) and the in-plane
equilibrium equation (2.48).
CHAPTER III
APPLICATION AND RESULTS
In this chapter, the theory and the methods developed in the
preceding chapters are applied to a variety of problems. A computer
program has been developed which is applicable to rectangular plates
with different boundary conditions and variation in stiffness. The
objective is to illustrate the application of the method to plates with
several types of variations in stiffness, as well as to the uniform stiff-
ness plates. Since solutions to the uniform stiffness plate are
known, it provides a good measure for “verifying the accuracy of the
SOlution procedure. For the plates considered, the solution is
Obtained for a few problems for all successive steps of loading from
zero load up to secondary buckling and the results are analyzed.
Convergence of the solution is checked and accuracy of the results
13 examined via comparison with known results when possible. Some
oPtilnzization problems are presented at the end of each section.
Chapter III is divided into two sections. Section 3.1
deals with force boundary conditions . In section 3.2, plates with
displacement boundary conditions are considered.
A square plate with a symmetrical variation in stiffness is
considered. Stiffness is symmetric with respect to both centerlines
a
nd diagonals as shown in Figure 3.1; thus only a quadrant of the
Pl
ate will be considered. The variation in stiffness is such that
52
53
in quadrant (I) of the plate (see Figure 3.1) the stiffness is a
function of x only. For example,in section 3.1 a parabolic variation
in stiffness is considered; the flexural stiffness can be
represented by a parabolic equation:
D(x) - DrER + 4(1-R)(§)21 (3.1)
where D is stiffness at point x
Dr is stiffness at center of the plate
R is ratio of edge stiffness to center stiffness
a is length of each edge of square plate
Ekate:
: Et3
Since bending stiffness, D =-—————-—- , and membrane stiff-
12(l-v )
ness, K - Et, are both present in the plate equations, it will make a
difference whether the stiffness variation is due to a variation in
13 (or in t.
a)
For the variation in E, with t constant, the D
variation and K variation will have the same pattern.
K1
Let Bi - KE" the ratio of membrane stiffness at point
r
i to membrane stiffenss at reference point, and let
D
6 = -2- be the bending stiffness ratio at corresponding
1 D
r
points. Then,
81 Bit/Err
a 2 3 3 2 1 °r 81 ' 51
i E t /E t
i r
54
‘X
\_ /
\ /
\ /
\ /
\ /
\ /
/
3 ,
A - - - A o
\
/ \
/ \
I \
/ \
/ \
/ \
1’.
A a A]
Y
I D
edge
Dcenter
Stiffness variation at section A-A
F
igure 3.1. Geometry and stiffness variation of square plate .
b) For the variation in
s =3
1 K
1'
D
51:3—
r
i-
55
t, with E constant,
Et
_123
Et t
r I'
Eti
12(1-02)
Et3
r
12(1-v2)
Both cases can be considered without any major difficulties.
In section 3.1 (excluding 3.1.4) case (a) is assumed,
and in 3.1.4 and the entire section 3.2, the variable thick-
ness case (b) is considered.
In terms of boundary conditions, two separate classes of
PrOblems are considered :
l - Force boundary condition
2 - Displacement boundary conditon
In order to avoid computational difficulties,the following
not1‘-d:l.mensional variables are introduced and frequently used in the
anal3'sis .
t'-
t/ti
== 3
D a
DI
1{/a
‘Y " 37/a
w 22 W/ti
where t
= unit thickness.
flexural stiffness of a uniform stiffness, unit thickness plate.
Inembrane rigidity of a uniform stiffness, unit thickness plate.
U = ua/ti 56
K0
'-
U “N"
2 2
V va/ti
N'; N '; N' = - -
( x y xy) (Nx, Ny, ny)/Ncr
* 2 k 2
N. N' N )=(a_.)(N. N' N )
( x’ y’ xy Do x’ y’ xy
N- N; N = N; N; N N
(x. y xy> ()2( y WM
—0 -o - = a o o
(Mx, My, MU) (Dot1)(Mx’ My, Mxy)
a2
0' = ( )o
Doti
‘6'= qa4/D t
o i
<9 cp/Nz
3.1 Force Boundary Condition
The previously mentioned plate case (a) is considered sub-
ject to a compressive normal in-plane force 'N' on all edges
(at):restrainton.in-plane displacements).The solution is obtained
for: three successive ranges of loading (membrane, buckling, and
POStbuckling), for different variations in stiffness, and for both
sinlPly'supported and fixed edges. The effect of a transverse load
is also examined.
The behavior of the plate within each range is observed.
For each range the in-plane stress resultants, the in-plane dis-
plaCements and the out-of-plane displacement are calculated and
Plots are given. The results of solution for variable stiffness
plates are compared with those for the uniform stiffness plate.
S
ection 3.1 deals with different phases of the behavior, as follows:
3.1.1 Membrane solution
3.1.2 Stability analysis
3.1.3 Postbuckling behavior
_w
.o.|
‘7
~,
,
‘L:~
o.
\.
. I
i.
‘
v‘.‘ ,
b
V“:
57
In order to have a common base in different cases of stiff-
ness variation for comparison of the results, the reference stiff-
ness, Dr’ is taken such that the volume under the stiffness curve
be constant for all cases. (i.e.,mean stiffness is constant)
The stiffness variation in this section is
x 2
D(x) - DrER + 4(1-R) (E) J a/2
+——-—1
as introduced in equation X
(3.1). The volume under
a/ 2
the curve over a quadrant II.
of the plate (see Figure
3/2
3.2) is vol.- 2] 9. D(x)dx
Substituting for 2 and
{r
90‘) we obtain De].- [D
D(x r
a/
vol . zpf 2(-a~x)[R+4(1-RX£)2]dx.
r o 2 a
Integrating produces Figure 3.2
2
V°1 =- ————a Dr (1 + SR
24 )'
This volume is to be constant for all cases; thus, the variation of
Dr W1th R must have the form,
D a constant
1: 1+5R
D
'I'
he base stiffness is chosen to be 63
o
= l for a uniform stiffness
plate. and the Dr for different stiffness ratios are tabulated below:
1/10
1/2
10
3.1 .1 Membrane Solution
58
1.7142857
.117647058
In this part, a flat plate with no initial deflection is
considered and the solution obtained.
Since w is zero, the solution
can be obtained by solving the compatibility equation (2.47) only.
Boundary Conditions
To determine the value of the stress function, <9, along
the boundary, recall Equations (2.12); along the x-edges, we have:
32
NX III J2 s-N (1)
3y
N 3.2};
KY “Exay a 0 (11)
From (1) we have
a
3 (~33) = -N
int eg: at ing
a
$5 =' ~Ny+cl +f1(x)
2
(P ‘3 lNZ + c1
b
“t from (11)
y + y fl(x) + f2(x) + c2
59
3.39:. 29.
By (3x) 0 + 3x constant along the edge.
df1(x) + df2(x)
dx dx
.-. _nggy
= constant
8x
From the above we can conclude:
2
=-N — _—
¢ -%-- + cly + c2
similarly, along y-edges, we will get:
-N 2
¢ = g + c3x + c4
The constants must be chosen such that the given boundary
conditions are satisfied. Since the second derivatives of Q
determine the resultants, and, in this case, the plate is symmetrical
with respect to its centerlines, ¢ can be chosen such that it will
be symmetric about both centerlines. Thus, arbitrarily choosing
w - 0 at the corners leads to:
along edge yi= 0
since at corner x = 0, w = 0, and at corner x = a, m — 0
or
¢ =«g-(ax - x2) (3.2)
60
figfi -
3x 2 (a 2x)
A‘E=N_a =
8x 2 at x O
fla-flé. =
8x 2 at x a
similarly, along edge x = 0 we get:
N 2
and (x = a,%¥ = - 5%);
similarly from equation (3.3) (y = 0, g; = if) and
(y =
Now, we are able to compute (p values at each boundary
node. Figure 3.3 (a) shows the geometry of the plate and location
of node points.
The
=' 5.9879807 Na
A
V
N
2.5852510 -58.708963 122.8929425 Q3
L a ‘
£11.4951922)
64
10 9 8 7
9 5 5 4 s 5
II II II II
a 3 ,2 5 5
II I I II
7 4 4 1 4
II I I II
5 3 2 3 5
II II . II 11
5 5 4 5
K
Node No %—--'EE
1 i
1 l.
2 3.07692
3 3.07692
l\/j 4 10'
5 10.
6- 10.
I 1.649484
II 6.4
Figure 3.4. Stiffness ratio at nodes and intermediate nodes
square plate R = 1/10.
65
solving the equations, we get:
The membrane resultants can be computed using Equation (2.12)
For example:
N
x(3)
”.28204873
< 62 ? < .23348206
.199145001
\
3y2 (3) h2
1.136927044N
Na2
= .233482-2(.199l45)+3/32
a 2
(29
=
similarly, stress resultants at each node are calculated;the results
for nodes shown, are given in Table 3.1.
Table 3.1 Stress function and stress resultant ratios at each.node.
Square plate with R = 1/10, v = .25, h = a/4.
Stress function
Node Cp/Na2 Nx/N Ny/N ny/N
1 .28204873 1.55413 1.55413 0.0
2 .23348206 1.09879 .95865 0.0
3 .100145001 1.13693 1.13693 .03205
4 .12500 1.0 .52857 0.0
5 .0937500 1.0 .67636 0.0
6 0 1.0 1.00 0.0
6
5 3
4 1
66
Equilibriumis satisfied along any section of the plate.
For example, along x = .25a, using the block approximation:
P = (Nx)5(%)(2) + (Nx)3(%)(2) + (Nx)2(-:‘-) = -.9999999 Na 3 -Na
and along x = .5a:
p = (Nx)4(%)(2) + (Nx)2(%)(2) + (Nx)l(%) = -.9999999 Na
which are both very close to in-plane resultant at the edge of the
block along x = 0:
P = -N x a = -Na
3.1.1.b The Same Problem With R - .5
The same square plate is considered except the ratio of
g edge stiffness a 5
stiffness 13’ R center stiffness
Resulting values of the stress function m, and the calculated
in-plane stress resultants at the nodes shown below, are listed in
Table 3.2.
67
Table 3.2 Stress function and membrane forces
(R = .5, h = a/4, v = .25)
Stress function
P°int gnga2 Nx/N Ny/N ny/N
1 .26057639 1.19625 1.19625 0.0
2 .22319361 1.04536 .97297 0.0
3 .19052623 1.02574 1.02574 9.01058
4 .125000 1.000 .8578 0.0
5 _.0937500 1.000 .90316 0.0
6 0.00000 1.000 1.000 0.00
As in the preceding problem, equilibrium is satisfied along all
sections of the plate.
3.1.l.c Uniform Stiffness Plate (R = l)*
It is anticipated that as the stiffness of the plate approaches
uniformity, the solution will converge to the known results for uni-
form stiffness plate. Thus, this case is considered as a measure
of verification. The results shown in Table 3.3, are as expected.
The membrane resultant ratios are equal to 1.0 everywhere and
shear stress is zero; these are exact values.
*Both E and thickness are uniform all over the plate.
68
Table 3.3 Stress function and in-plane resultant ratios.
(R - 1, h =-%)
Point gglNaZ Nx/ N Ny/ N ny/ N
1 .2500000 1.00 1.00 0.0
2 .21875000 1.00 1.00 0.0
3 .1975000 1.00 1.00 0.0
4 .125000 1.00 1.00 0.0
5 .0937500 1.00 1.00 0.0
6 .0000 1.00 1.00 0.0
L1.1.d The Same Problem as in 3.1.1.a With R = 10
In contrast to the previous problems, in this case the
stiffness is increasing from center to edges and at the edges it
Results are shown in Table 3.4.
18 ten times stiffer than at the center.
Equilibrium is satisfied along all sections. 2
4 2 1
Table 3.4 Stress function and membrane force ratios, square plate
R " 1.0, h a % 9 V a '25
\
._Efleggpc g/Na2 Nx/N Ny/N ny/N
.1 .21667651 .32022 .32022 0.0
2 .20666961 .78443 1.14660 0.0
3 .18215623 1.02229 1.02229 -.03332
4 .125000 1.00 1.38657 0.0
5 .09375000 1.00 1.171 0.0
L11-ji .000 1.00 1.00 0.0
69
Improvement of the solutions
Since this is an approximate method, it is desirable to study
the convergence of the solution with increasing numbers of node points
(decreasing grid spacing). In this section the same
problems are solved using finer grid spacings (h = g at h = i— .
Solutions to all four problems are obtained. The node arrangements
are shown in appendices 3.1 and 3.2.
The convergence of the solution for the case R = 1/10
is shown in Table 3.5 and illustrated in Figure 3.8.
Table 3.5 Convergence of the solution, square plate R . 1/10, v = .25
Grid ‘65 -values by extra) N , at extrapolation 3 pt.
spacing at node 1 olation noée 1 results * extrap-
h/a olation
% .28204873 1.55413
.27472129 1.59060
% .27655315 1.58149 1.59026
.27465677 1.59028
T15 . 27513087 1 . 58809
L\ _L
* Iiiczhardson's extrapolation.
70
2 .
cp/Na Nx/N
.28 «-
«1.65
T 1
.275 ~ 1.6
NX
41.55
.27 5
4 8 l6 a/h
Figure 3.5 Convergence of membrane solution
§;J;Ll:l. Analysis of Results from Membrane Solution.
i) Uniform stiffness plate (R = 1)
Compared to theoretical values, in this case the finite
difference solution gives accurate results. The
difference operators agree exactly with the usual
difference operator for uniform stiffness plate.
Since the m function is parabolic in this case, the
difference approximationsto the second and higher
derivatives of' m are exact, and the solution is exact,
a
even for h =-Z. The in-plane stress resultants are uniform
71
over all the plate and there is no shear stress, as
expected. Distribution of Nx/N is shown in Figure
(3.6). The displacements vary linearly from the symmetric
centerlines (see Figure 3.7), and the strain is constant,
which agrees with elasticity theory. The solution
obtained with the grid-spacing, h =-% is exact, as are
solutions with finer grid-spacings.
(L.
niece
N/N=1
X
(a) Plan ' (b) Contours of Nx/N
F — r“ w r-
9-9 8-3 C-C 0-0 E-E
N
(c) Profiles of-i?
2 1 -
SCHLE
FIEUINe 3.6. Force distribution for undeflected square plate R = 1,
Figure 3.7. Contours of in-plane displacement (U' =
ii)
72
.375 .28 .1875.094 0.
K
N) for
undeflected square plate, R = 1, v = .25.
Square plate with R = 1/10.
In this case, because of variation in stiffness, the
stresses vary over the plate. Thus,the m function
is not a smooth, parabolic one as it was in the uniform
stiffness case and the solution would be an approximate
one. Solutions with finer grid spacing were compared
(see Table 3.5) and they show fairly good convergence.
Study of in-plane resultants in Figure (3.8) shows the
expected behavior,with a shifting of the load toward
the stiffer parts of the plate.
Equilibrium is satisfied along any arbitrary section
of the plate. Shear stress is zero at points of symmetry
but at nonsymmetric points, because of the load shifting
process, there is a small shear force created, as expected..
73
D...—
m
a) Plan b)
—
L.
8-9 8-8 C-C
c) Profiles of Nx/N
Contours of Nx/N
1
[3-0
.5.4 .3 .2 .1
d) U'-displacement
.75
'1.0
\ ‘2
\ ~——
\ 1.
. P
1 (7115.0.
25
F
igul‘e 3.8. Distribution of in-plane force (Nx/N) and displacement
K
(U' = U 3;) square plate, R = 1/10.
74
In-plane displacement patterns are illustrated in
Figure(3-8 d)-It is observed that the displacement
normal to the edge is increasing as we approach the corners,
at which stiffness is less than at the center. This behavior
seems reasonable since there is no displacement restraint
along the edges.
iii) Square plate with R = 1/2
iv)
This problem can provide a good check on solutions,
because it lies between two previous cases. As
we go from the R = 1/10 case to the uniform stiffness
case (R = 1), we would expect the solutions to approach
the uniform stiffness results. Investigation of the
results in Figure(3.9) along with the results found by different
grid-spacings and comparing with cases (i) and (ii), indicates:
a) The convergence as the grid spacing becomes finer
b) Results for the stresses and displacements trend
toward those for the uniform thickness case as R is
changed from 1/10 to 1/2.
Square plate with R = 10.
Solution for this case shows convergence as the grid
spacing is decreased, and it also agreeswith previous
results in that:
a) Load in the plate shifts toward the edges which
are stiffer as illustrated in Figure (3.10).
.b) Displacement normal to the edge is a little larger
at center of the edge and decreases toward the
corner; (see Figure 3.10 d).
75
' J
1.0
1.10
Q l 1 1 F
b) Contours of Ng/N
a) Plan
7"
_J
R-H B-B C-C 0-0 E-E
c) Profiles of Nx/N
. .__J
.35.3 .2 .1
d) U'-disp1acement
Figure 3.9. Distribution of in-plane force (Nk
K
(U' = U 7;) square plate, R = 1/2.
/N) and displacement
fr
II
‘1.
‘ V
5116 3
I
76
I I 1
\\
\ 1.25
‘~
1.0
..75
l l ///’/"-—~‘;:::EEEE .50
9 édd 4'—
R
b) Contours of Nx/N
a) Plan
_J _J
8-9 8-8 8-8 . 00-0 E—E
c) Profiles of Nx/N
2 13-1-8
SCHL
d) U'displacement
Figure 3.10. Distribution of in-plane force (Ng/N) and displacement
K
(U' = U'jf) square plate R = 10.
77
3.1.2 Buckling
To study buckling and determine the critical value of the
applied compressive force, NC , the equilibrium Equation (2.14)
r
is employed. The left hand side of this equation is approximated
by the operator Figure (2.8) and the right hand side is represented
by the difference operator Figure (2.9). The stress function,
9, values obtained in the membrane solution for the plate with
no lateral displacement are used in the equilibrium equation. This
is acceptable,because in the stability analysis we are seeking
bifurcation of an initially undeflected plate.
Boundary conditions
In the case of buckling and postbuckling, the out-of-plane
displacement, w, has a considerable effect on the solution and
out-of-plane boundary conditions must be considered. In the simply
supported case, along edge x = 0, we have:
i) w = 0
82w 32w
11) M --D(—-+ v —) = 0
x 2 2
3x 8y
But --—5 = 0, so that (ii) becomes ——§-= 0. Using a difference
3y 3x
approximation,
32w _ W1+1 ' 2"1 + w1-1
8x2 h2
we will get at point B (Figure 3.3),
since wB
In the case of fixed support, the conditions are:
i) w B 0
11) gg=o
using a central difference approximation for slope we have:
8x 2h b
3.1.2.a Buckling of Plate with R = 1/10.
Let us consider problem a) again (R
supported edges and h 31%.
= 1/10) with simply
In the right hand side of equation (2.14) the m values have
already been obtained for the initially flat plate.
Applying operator Figure (2.8) at each point, and substituting
Q values in the right hand side, the following eigenvalue problem
is obtained.
P
1.071875 -5.5375 6.5725 w3
_J 1 J b
2
or: calling A = 32-
D
r
r
14.9375-.388533ZA ’ -20.52500+ .38853321 4.2875 1
L 1.071875+3004006l -5.5375+ .1421158A
-5.13125+.06867411 10.86875- .25717881 -5.5375+ .11983.GK
F' )
14.9375 -2o.52500 4.28751rw 3.885332 -.3885332
-5.l3125 10.86875 -5.5375 =r -.0686741 .2571788
-.0040060 -.l421158
6 . 57 25- . 284231?
0 ‘1 w
.1198306
.2842316
W
W
l
2
3
79
Solving leads to these eigenvalues,
f
5.20748
20.42753
>4
II
A
46.71190
\
01'
5.20749
N =1—‘2' = 20.42753
Nln
46.71190
The corresponding eigenvectors are:
11_ .12. .12.
1.0 1.0 -.69810
.84441 .36007 1.0
.63179 -.57565 .03122
To check convergence, the same problem was solved with finer grid
a a
spacings, (h -'5) and (h 16).
The eigenvalues shown in Table 3.6 were obtained.
The con-
vergence is fairly good and extrapolation flnproves the results
further.
A1 6.
5. ~—————___11 ._1
4._
'r. . r
4 8 16
Figure 3.11. Convergence of eigenvalue, R = 1/10.
a/h
80
Table 3.6 Eigenvalues of s-s square plate, R = 1/10,with different
grid spacings.
h/a A1 2 pt. extrapolation 3 pt. extrapolation
1
1/4 5.20748
4.727158
1/8 4.84724 4.855814
4.847773
1/16 4,84764
The first mode shapesin the three solutions are very close to each other.
3.1.2.b Buckling of Uniform Stiffness Plate
As before, a plate with uniform stiffness is considered;
the solution can be used for an accuracy test since the exact solution
is known. The first eigenvalue obtained using different grid
spacings is tabulated below and compared with the exact solution.
Table 3.7 Comparison of first eigenvalues of different solution (R - 1)
h/a Al 21 exact difference 2 Pt. extrqr- 3 Pt. extrap-
olation Aolation
1/4 18.74517 19.7392088 5%
19.734057
1/8 19.48684 1.2
N
19.739197
19.738873
1/16 19.67587 .3%
81
It can be seen that the results are very close to the exact values
given in Reference (44);
Na 2
N = 2V or, A =-——— = 2N = 19.7392088
The solution is satisfactorily converging to the exact value.
Three point extrapolation results in an error of less than 10—6.
Table (3.8) shows the critical loads obtained for two dif-
ferent loading cases and gives comparison with previous work as
well as exact values. It can be seen that even with 8 x 8 nodes
the results are satisfactorily within engineering accuracy. Mbre
accurate results can be obtained by extrapolation. In Table 3.8,
the results given by Clough and Felippa are obtained by the finite
element method and Dawe used the "discrete element displacement
method", which in principle is the same as finite element method.
Eigenvalue problems for other cases (R = 10, and R = 1/2)
were solved, and the convergence was examined with increasingly
finer grids. Details will be discussed later.
3.1.2.1 General Buckling
So far, the problem was considered symmetric with respect
to both centerlines and the diagonals. Therefore, the solutions
are limited to symmetrical modes of buckling only and nonsymmetrical
82
cowumaoamuuxo usaoalm mo oasmom
um
«smoooo.~
ommm.a oH\H
mmmmmmmm.a
.N wmm.a omqqnm.a w\H
mumma.~ scammoumaoo
mwm.H mmmmm.a «\H Hoaxmlam enemas:
oH\H
.c Hmo.e mmm.m om.m m\a cofimmouaaoo
oNH.¢ m~m.m «\H Hmfixmasa auouwsa
coquaaom Ado mmmwaom soHumHo nuanmom m\:
Hmowmmmao can nwsoao AnHV mama Iamuuxm hm uawmmum ouam 5mm: ommu mafivooq
. nu: Ho
mmm. u 9 H u an .mml I «z .mumaa mumavm wouuonmsm manawm mo mumoa Hmowufiuo w.m NHANH
83
modes are absent. This means that 12 in this solution is not
the eigenvalue corresponding to the second mode, but it represnts
the eigenvalue corresponding to the second symmetric mode of buckling.
In this particular case the 2nd symmetric mode corresponds to the
5th general buckling mode.
To obtain more accurate results, the plate was considered
without imposing any symmetry. Thus, all possible degrees of
freedom were allowed, within the restrictions imposed by the choice
of grid spacing. The solution for each case was Obtained and the
buckling modes and corresponding critical loads are studied in the
next section.
3.1.2.143 General Buckling_of Square Plate with R = 1/10.
1) Simply-supported boundaries.
Problem a is solved for general buckling (no
symmetry imposed) with h = a/8 and assuming simple
support along all edges. The first few modes of buckling
and the corresponding eigenvalues are shown in Figure
(3.12)
ii) Clamped edge.
The same problem as in (i), but with all edges clamped,
was solved. The first six modes of buckling and the
corresponding eigenvalues are shown in Figure (3.13).
84
-— +
+
+ .—
* * 3 44 * 44
N1 = 19.39 N2 9. N3 39.
+ _ + _
+ +
- +
* 61,61 N* 75 15 N* 76 52
N4 . 5 . 6 .
Figure 3.12. Modes of buckling of square plate, R = 1/10,simp1y supported.
+
+
- +
* 42 24 N* - 66 60 N* — 66 60
N1 . 2 . 3 .
+
+ - + _ +
_ + - + _
* 91 60 N* 103 32 * 3
N4 . 5 . N6 . 118.28
Figure 3.13. Modes of buckling of square plate, R = 1/10, clamped,
(N* - Na2/Do).
cur
85
+
+ + -
* 20 04 N* 46 73 N* 46 73
N1 . 2 . 3 .
+
+ - +
— + +
N* 74 79 N* - 88 78 N* - 89 28
4 O 5 _ O 6 - 0
Figure 3.14. Modes of buckling of square plate, R = l/2,simply supported.
- +
+ +
N* 49 56 N* 82 0 *
1 . . 2 . 3 N3 - 82.03
+ +
+ —
"' " +
* 114 20 N* 12 o *
N4 . 5 7.3 N6 = 137.84
Figure 3.15- Moges of buckling of square plate, R = 1/2, clamped,
(N a NaZ/Do).
86
' *
N: = 19.48684 N: = 47.23375 N3 = 47.23375
+ — + +
- +
+
* * *
N4 = 74.98066 N5 = 88.75994 N6 . 88.75994
Figure 3.16. Modes of buckling of square plate, R = 1,simply supported.
+
+ + -
N* 49 567 N* 82 131 N* - 82 131
l O 2 . 3 .- O
+ - + +
- + +
* * *
N4 = 112.693 N5 = 124.778 N6 = 137.127
Figure 3.17. Moges of buckling of square plate, R = 1, clamped,
(N a NaZ/Do).
*
N
l
= 18.42
* 70 04
N4 .
Figure 3.18. Modes of buckling of square plate R = 10,simp1y supported.
* 52 09
Nl .
N* 101 78
4 .
Figure 3.19. Modes of buckling of square plate,R =
(N* = Na2/Do)
87
+
N* - 47 00
2 " o
+
+
N* - 83 01
5 - .
iéz‘é'.
+
N* - 81 37
2 - o
+
+
N* 11 02
5 3.
*
N = 47.00
N* 90 24
6 .
+
* 8
N3 - 1.37
+
N* 1 3 16
6 3 .
10, clamped,
A‘s _
.5 . «a
o .
a.
d— 44“
88
3.1.2.1(b,c,d)
The plates with R = 1/2, R = l and R = 10 were solved
for both simply supported and clamped boundaries; the buckling
modes and critical loads are illustrated in Figure (3.14)
through (3.19).
3.1.2.2 Analysis of the Results
A. Simply-supported edges:
1) Uniform stiffness plate.
As discussed earlier, the critical load obtained for
this problem agrees very well with the exact solution.
Mbdes of buckling based on theoretical solutions are
combinations
°. 10
‘7‘ 1
1/2
1/10
Figure 3.20.
Deflected shape for s-s and clamped plate.
(N = 2.40 Ncr)
91
B. Clamped Boundaries.
1) Uniform stiffness plate.
The first eigenvalue is found to be 49.56763 which
is very close to the exact value, obtained by Levy (28),
of 5.0378
2 = 49.71319.
Table 3.9 shows convergence and accuracy of the results,
and gives comparison with some previous works and the exact value.
Extrapolation of the results shows an accuracy of about .22.
Table 3.9. Critical load ‘N* = N
2
a
or NZD
bi-axial uniform load, R = l, v = .316.
of clamped square plate under
Grid-spacing» A Extrapolation Levy (28) Clough & Classical
(hlg) 1 Felippa Solution
1
‘2 5.625
.% 5.02225 5.037 5.399 5.31
5.29921
1
IE’ 5.22997
The first buckling mode agrees almost exactly with the theoretical
mode shape (1- cos 22E5, for m = 1.
ii) Case R = 1/10.
111)
See dashed curve in Figure (3.20).
Similar to the s-s case, in the central region the plate
remains flat and sharper curvature occurs near the edges.
Case R = 1/2.
As anticipated, results obtained for this problem lie
between cases (i) and (ii) supporting the validity of
the solutions.
92
iv) Case R = 10.'
As before, sharp curvature is observed in the central
region of the plate due to smaller stiffness of that
region.
afzer
be in
Theref
be an;
sewn:
and t
the e
93
3.1.3 POSTBUCKLING
3.1.3.1 General Procedure
As discussed in the preceding chapter, a plate will stabilize
after first buckling. Immediately after buckling, the plate would 4
be in a state 0f Stable equilibrium with moderately large deflection.
Therefore, the large deflection theory discussed in chapter 2 can
be employed to study postbuckling behavior of the plate up to the
secondary buckling point.
In this section, the procedure followed will be discussed
and the results will be analyzed.
For a solution to the large deflection behavior of a plate,
the equilibrium and compatibility conditions must be satisfied;
both equations are coupled in w and m. In this case, in addition
to the variation in stiffness, geometrical nonlinearity caused by
large deflection will also be involved.
To solve these coupled, nonlinear equations, an iterative
technique is employed. A schematic flow chart of the procedure
is given in Figure (3.21). The steps of the procedure in Figure (3.21)
are as follows:
Step 1 - Solve the equilibrium equation (2.14)
The left hand side is approximated by the operator of Figure
(2.8),which will be applied at each node to form the matrix
[Aw]. To calculate the right hand side of this equation,
initial values for m and w are needed. w is assumed
consistent with the first buckling mode shape, and the
94
Apply load N, (N > Ncr)
Input @ and w based on previous
’ solutions
Solve equilibrium equation
(2.14)
A A Get new w
Solve compatibility equation
(2.15)
Get new m
‘ no ////:::{:\\\\\g Yes
+1
convergenc
Increase Load
N = N + AN
Yes Nx
axial compression téI Itc
. t(x)
be maximum for a
Figure 3.57.
constant amount of material.
For a square plate of thickness, t, and sides, a, the total
volume of material used is
?'= azt-
For a variable thickness plate, if the variation in thick-
ness is linear, the thickness t(x) at point x is
t-t
c e
t(x) te + a/2 x
where
tC = thickness at center
149
te = thickness along the edges
I:
Let '32 = RT= Ratio of edge thickness to thickness at
c
center. Then:
_ £15
t(x) - tCERT + (l-RT) a
and the volume is
a
- 5. 2x 2x
v 4 f0 (a-Zx)tc(RI+-:r--RT:;)dx
azt
_._ 2 l + ZRT _ c
v - 4a tCG-jf§- )- 3 (l + ZRT). (i)
If the volume is limited to the original value, then
t —7I;2L__— ; for '3 = azti (ti = unit thickness)
c a (1+2RT)
3t1
t:c a (1+2RT) (11)
To maintain a constant volume, the center thickness must
vary with ratio RT, according to equation (ii).
For example:
RI tc/t1 v
2
.l 2.5 a T1
_1_ 2 11
2 2
1. 1. "
2 .6 II
150
It is evident that as the ratio, RT, and the center thickness,
tc, change, the in-plane and flexural stiffness of the
plate at each node will change.
A subprogram is PrOVidEd to compute center thickness and
stiffness at each node in each trial for a given RT.
The buckling problems were solved for a family of given
ratios, and the variation of critical load with ratio RT,
is shown in Figure (3.58). In this problem the plate is
simply supported along all edges and v = .316.
Graph (3.58 a) indicates that the maximum critical edge
load for a given volume of material occurs with RT : .2._
To find a more accurate value, the trial is continued with
finer intervals between RT = .15 and RT - .25. The larger
scale graph (3.58 b) shows that the'maximum
critical load is (N* )
cr imun = 25.66, corresponding
nl n
m
to RT= = .22. Thus, consider a simply supported square
C
plate under bi-axial compression.
From the stability point of view,the minimum material can
be used if RT ‘ .22 for the linear thickness variation
25.66-19.48
19.48
over the buckling load for a uniform thickness plate of the
introduced in Figure (3.57). An increase of - 31.6%
same volume results.
151
11935
i) The result obtained above is not the absolute optimum.
For the given amount of material, a still greater critical
edge force could no doubt be obtained with a thickness
variation other than linear. This form was used for
practical simplicity.
ii) For other types of loading or boundary conditions,
appropriate forms could be proposed and analyzed by
trial.
EXAMPLE 2
The plate in this example is the same as in example 1,
but with clamped edges; results are shown in Figure (3.59) .
In this case, the maximum critical load is found for an RT
of approximately 0.8.
Figure (3.59 b) is obtained by taking finer intervals
(AR)?8 .01) between .6 and .8. This graph indicates the maximum
critical load is (Nzr)max. = 473.962, corresponding to RT =.71.
51.20-49.56
49.56
over the plate of uniform thickness. In the case of fixed support,
The results show an increase of = 3.3%
from the stability point of view, it is not worthwhile constructing
-a plate of variable thickness. The effect of variation in thickness
on other aspects of the problem such as internal forces, bending
moments and lateral deflection will be analyzed later in this chapter.
152
28.00 32.00
4.00
25.1
20.00
c115.00
.00 0.40 0180 11120 1160 2100 2140
a)
26.80
2§.SO
.12 0116 0120 0124 0128
RT
c26.40
b)
Figure 3.58. Critical load vs RT,for a simply supported square
plate under Bi-axial compression, N.
153
9-00 47.00 55.00
31.00
c23.00
.20 1100 1180 2160 3140 4120 43100
RT
53.80
53.40
1
53.20
.50 0165 0170 0175 0.00
b) RT
c153.00
t
Figure 3.59. Critical load VS RT= 7?-, for square plate of constant
C
volume under Bi-axial compression;clamped edges.
154
3.1.5 SUMMARY
In Section 3.1 a method was developed for force boundary
condition problems, and the related computer program was applied to
a variable stiffness plate and the results were discussed extensively.
It seems necessary to emphasize that the plates discussed
in Section 3.1.1 through 3.1.3 have uniform thickness with varying
E, so that patterns of membrane and bending stresses follow exactly
the pattern of in-plane forces and bending moments respectively.
Solutions to similar problems are not available in the literature;
however, comparison with the uniform stiffness plate,as a special
case, supports the accuracy of the solutions. Convergence
of the solutions with an increasing number of nodes strengthens con-
fidence in the method.
In Section 3.1.4 the weight-saving advantage of a variable
thickness plate,from the stability point of view, was demonstrated as
an example. However, one can apply the optimization to any possible
aspects of stress or strain as desired.
Figure (3.60) shows variation of central deflection with load
for all cases. The graphs show that the plates with less stiff
edges undergo larger lateral deflection because,in the postbuckling
range,the main portion of the in-plane load is carried by portions
of the plate near the edges.
155
0.00
400
I
‘b.00 0100 110 .40 0120
D 2
CENTRHL w t 1
a) s-s
100.00
A
00.00
094m
1
40.00
20.00
cp.00
0100 4100
i
.00 1100 2100
CENTRRL w / t
b) Clamped
F
iiiure 3.60. Central deflection for different R values, s-s and
clamped edges.
156
3.2 DISPLACEMENT BOUNDARY CONDITION
In this section, problems with specified in-plane displace—
ments along the boundaries will be considered.
A.computer program hasbeen developed to solve this type of
problem based on the procedure discussed in section 2.3.
Plates with uniform thickness are analyzed, as well as variable
thickness plates. The results are compared with previous works or
exact solutions when they are available. As an example, optimization
of the thickness variation, from the stability point of view, is
also considered.
a - Geometry and Loading;Conditions
For an example, let us take a plate under uniform edge
displacement due to thermal load, and examine the membrane, buckling,
and postbuckling behavior of plates of different thickness variation,
‘with both simply-supported and clamped 'boundaries. Figure (3.61)
shows a plate, surrounded by a
g 0 .4
rigid frame that undergoes a r‘ ’T
////////j//Z I "
temperature change of either / /
/ /
expansion or contraction. g a
/
/
The strain in the frame will a 5
M /
impose a state of displace- 5 j
/ /
ment on the plate edges. //////[///[]/ I
The strain in the frame is
s: = 011(A'1‘)
Figure 3.61. Plan
157
where aT= coefficient of thermal expansion of the frame and
AT a temperature change.
Egg
1) It is assumed that the effects of the reaction forces
of the plate edges are negligible, so that the stresses
created within the frame have negligible effect on the
strain and displacements of the frame.
ii) Because of symmetry, we analyse only a quadrant of
the plate.
31 __32___33___34 35 36
r' I” T-_'T—_‘
: 1 I I l i
301__2 22| 2l 24 25
I
1
|
29L__29 13 14 15 lb
1
l
28i..__12 17- 7 8 9
|
|
27 I 18
r—-- 11 6 3 4
l
I
264 17. 10 5 2 1 _
Figure 3.62. Node arrangement, square plate, h = a/8
Since, because of symmetry, the centerlines of the
plate will not move, displacements enforced along the
edges will be as follows:
1 - Along edge x - 0, u displacement would be constant
0.
and equal to - ETaAT.
Table 3.10 .
2 — Because of constant strain (in frame),
158
the v-displace-
ment along x-edge will 'be linear with respectto
y, so that v =0TAT(y - g).
The boundary displacements are tabulated in Table
(3.10).
Boundary displacementfbr square plate of Figure 3.62
Point u/-JamTAT v/- aaTATfi
__7r—_ “3F“
17 1. 0.
18 1. .25
19 1. .5
20 1. .75
21 1. l.
22 .75 1.
23 .5 1.
24 .25 1.
25 0. 1.
Based on the given data,
analysed step-by step in the following pages.
the problem is solved and the results
159
3.2.1 MEMBRANE SOLUTION
In the solution for in-plane forces and displacements,
assuming zero w—displacement, the in—plane equilibrium equations
can be approximated by finite differences to obtain a set of linear
equations in u and v. (See equation (2.43)).
The equations are solved and in-plane forces and displace-
ments lead us to following conclusions.
a) In case of a uniform thickness plate, the equilibrium
equations (2.16) and (2.17) are linear differential
equations, and in this symmetric case the solution leads
to exact values of membrane stress resultants and dis-
placements, even with very coarse mesh sizes.
The theoretical solution predicts constant strain in both
" directions.
-uo -2uo a
.y = 8.. = .72 = T W (81> : 8)
and the in-plane stress resultants are
-2uo _ 110
N = N a Etz [e + ve ] = Et [(1+v)(————)] = ZEt --
y X l-V x y 1_v2 a (l-v) a
* 32 -2Eti no 82 (1-V2)12 “ca
3 N —— = __ = _ _—
Nx x Do (l-v) a. ( 3 ) (1+v)(24) t2
Eti i
uoa
or, calling -—§-= UO
‘1
*
Nx
U0 = -24(l+v) *
N
and for v = .316, x = 31.58
C.‘
O
160
Similarly,
N:
I = -31.58
N = 0
X?
and
a 2uO a
u = u(x) = 8(3 -x) = T (-2- -x)
a 2v0 a
v = v = e<§ -y) = —a— (5 -y)
These values agree exactly with the solution obtained with
the computer program listed in Appendix:C. Contours of the membrane
force, N* , the principal stress and U-displacement are shown in Figure
(3.63). N& and v-displacement can be obtained considering symmetry.
N; = 31.58 o'= 31.53
U0 U0
1. 25 50 25 o
0) 11151101111111: 18011c1:.11’;‘{/Uo b1 911111019111. STRESS o'/U° c) U-OISPLRCEHENT
FIGURE 3.63.CONTOURS OF HEHSRRNE FORCE,PRINCIPRL STRESS,RNO U-DISPLRCEHENT
UNDEFLECTED SDURRE PLRTE. RT=I.
161
SOLUTION OF PLATES WITH VARIABLE THICKNESS
To observe the behavior of a variable thickness plate,two
different problems were solved. One with RT = .25, and the other had
RT - 2, and the results studied. Parameter RT =
edge thickness
center thickness
is the thickness ratio.
b)
Square plate with RT 8 1/4
Figure (3.64) illustrates membrane forces, principal
stresses and in-plane displacements;from the figures, we can
conclude the following.
i)
ii)
iii)
Figure (3.64 a) shows shifting of the load carrying
toward the center where the plate is thicker. The
shifting of the load becomes greater as we move
toward the center, and it is nearly uniform along
the edge.
Although there are larger in-plane stress resultants
in the central region, because the thickness is
large,the stress is smaller. Figure (3.64 b) shows
the principal stress within the plate. It is
observed that the minimum stress exists at the
center and it becomes larger toward the thin edges.
In Figure (3.64 c), the wide-spacing of contours of in-
plane displacement U in the central region, indicates
small strain corresponding to smaller stress in that
region.
20
25
,30
[1..
01 1151101111112 FORCEM; /Uo
162
50 40
b1 PRINCIPRL STRESS 0'71]o
Q—SO
4O
30
.25
r—jZO
l. .75 .50 .25
c1 U-D ISPLRCEHENT
FIGURE 364. CONTOURS OF HEHBRRNE FORCE.PRINCIPRL STRESS.RNO U-DISPLRCEI‘IENT
UNDEFLECTED SDURRE PLRTE. RT=I/4
c) Square plate with RT - 2.
The same problem is considered except with
RT . edge thickness
center thickness
are shown in Figure (3.65).
- 2. and the results obtained
following conclusions.
This leads to the
1) Supporting discussion ofthe preceding section, less
load is carried by the thin central region, and as
we move toward the edges, more load is transmitted.
11) Figure (3.65 b)
region.
shows larger stresses in the central
iii) Contours of in-plane displacements are consistent
with part (iix i.e.,larger strain occurs in the central
region,due to larger stresses.
Clearly, the out—of-plane support condition has
no effect on the membrane solution of the plate.
163
It should be noted that Figures (3.63), (3.64) and (3.65)
correspond to plates with RT = 1, 1/4 and 2, respectively,
but with thickness such that total material volume is the
same in each case.
//’-_______.——~*—135 I‘I““‘30
30 35
A 25 , f 40 _
1 . .75 .5 .25
01 1151101111111: F°R°E"‘;/Uo b1 PRINCIPHL STRESS o"/U° c1 u-msrmceneur
FIGURE 3.65.CONTOURS OF HEHBRRNE FORCE.PRINCIPRL STRESS.RND U-DISPLRCEHENT
UNDEFLECTED SDURRE PLRTE. RT=2.
164
3.2.2 BUCKLING SOLUTION
As is discussed in Section (3.1), the buckling solution
is based on the eigensolution of the out-of-plane equilibrium
equation (2.40) (h1the right hand side of this equation,the effect
of w on in-plane forces is considered to be zero. Thus, in the matrix
[Bw] of Equation (2.49), the in—plane forces obtained from the
membrane solution of the unbent plate will be used. The eigenproblem
will be solved, giving the critical edge displacements and the
corresponding buckled mode shapes.
Following-are solutions to some stability problems.
3.2.2.1 Convergence Check and Comparison
In order to check the convergence and consequently the
accuracy of the buckling solution, a uniform thickness plate
(problem 3.2 a) was solved using different mesh sizes. The critical
loads obtained in each solution are shown in Table (3.11) for
s-s boundaries. Also,the values of critical displacements are
compared to exact values.
Convergence of the solution is graphically illustrated in
Figure (3.66) for s-s boundaries.
Table (3.12) and graph (3.67) show the convergence pattern
of critical displacement for the case of boundaries fixed against
out-of-plane displacement.
The convergence tables show that the result obtained using
h 8:- in the simply supported case and h - {12- in the clamped
case are sufficiently accurate for engineering purposes.
165
More accurate results can be predicted by extrapolation.
The critical edge displacement obtained from extrapolation of the
last 3 lines has an accuracy of £H)0016 in the simply-supported
case and .0013 in case of clamped edges, compared to exact results.
Table 3.11 .
using different mesh sizes.
Critical displacements for a simply-supported plate
Uniform stiffness plate.
v = .316.
Mesh size Present solution Exact(1) Difference Two point 3 point
(h/a) Uo % extrap- extrap-
ér olation olation
1- 5935 5
4 .
.6248
§ .61698 1.2 .624972
.62497 .62495
1
"1‘2— “62141 5 .62498
.624975
1
16' .62297 .3
ID
0
61
. . . . Exact
O
(D
o.‘ //_——‘
‘2
3
O"
8
“11.00 4100 0100 12.00 {0.00 a/ h
Figure 3.66.
Convergence of buckling solution.
166
Table 3. 12 . Critical displacements for a clamped square plate
using different mesh sizes. Uniform stiffness. v = .316.
Mesh size Present Solution Exact (1) Difference 2 point 3 point
(h/a) Ufi Z extrap- extrap-
r 'olation olation
% 1.34062 19.2
1.6456
_1_
8 1.56939 5.4 1.65620
1.6593 1.6550
12 1.61697 2.5 1.65709
1.6566
1. 1.6343 1.5
16
2
£1
Exact
ID
"3.
2
3
ID
‘3 .
“0.00 4100 0100 12.00 16.00 a/11
Figure 3.67. Convergence of buckling solution.
(1)The exact value of critical displacement is derived on page 167.
buc‘:
pla
all
for
bu
V
t:
167
Since it is assumed that the plate is perfectly flat before
buckling, the in-plane displacements are linear within the entire
plate (constant strain), and the membrane forces will be the same at
all points. Thus, critical displacements can be related to edge
forces as
along x-edges
1
ex Et [Nx - vNy]
but in this case
-u
cr _ g
(Ex)cr a/2 ’ and Nx - Ny Ncr
therefore
-a
ucr 2Et [l-VJNcr
anDr
For s-s square plate under bi-axial compression Ncr = --§—-.
a
[see Table 3.8]
Thus
u :3 l2 D— (1-\,) = i112..— .
cr a Et 12a(l+v) ’
u §__= .822467 .
or t2 1+V
.822467
Uo = -I:;- or for v = .316, UO = .62497
cr cr
nzDr
similarly, considering Ncr = 5.31 2 for clamped boundaries
a
[see Table 3.9 ],the critical edge displacement is
Uo = 1.6593 for v = .316
cr
"|\I
.
l
:2:—
I"
U'fi
168
EXAMPLE 1
Tofurther check the accuracy of the method, the solution
to problem (9.13) of Reference (44 ) is examined. In this example,
along the edges x = 0, and x = a, the
u-displacement is prevented. In the
y-direction, constant strain, a,
is assumed. Therefore, v = 6(3 - y) and
assuming that centerlines of the plate
coincide with axes of symmetry
during deformation, (6 a 353-. /
--_-——————
No u-displacement is
Y
allowed alon ed es = 0,
g g y Figure 3.68. Plan
and y = a, and with regard
to out of plane displacements, all edges are taken as simply supported.
The boundary data.are tabulated below. For the node arrange-
ment, see Figure ( 3.62).
Node u v w
17 0 0 0
l8 0 .25vo 0
19 0 .5 v0 0
20 0 .375vo 0
21 0 v0 0
22 0 v0 0
23 0 v0 0
24 0 v 0
25 0 v 0
169
Using the data for input, the buckling solution was obtained.
The critical uniaxial edge displacement was found to be
2
v0 = 1.23033 Ea— .
cr
This is to be compared to Timoshenko's results (44), obtained by
au1 energy method. Remember that he was taking the total edge strain to be
2v
ey, while in our case the compressive strain e = ——9
9
cr a
2
e . 2v = 2.46066 L.
cr __.°cr 2
a a 2
. _ h
Timoshenko found ecr - .632 —a—2 but he used a plate
of sides 2a=§'and thickness h. By substituting t. for h, and
-' 2 2
éfora in Timoshenko's result it becomes e = .632 h— - .632 t—-—-— =
2 2 cr 2 - 2
1: «'8 (a/2)
2'528122 . The results differ by 2.6%, fairly close
a
for such a coarse mesh size (h = 3).
3.2.2.2 Optimization Analysis
Using the same data as in problem 3.2.a, and solving for
the eigenvalues, two series of solutions are obtained for simply-
supported and for fixed out-of-plane boundary conditions—considering
different thickness variations. The variation in thickness is
taken to be linear, as in Section 3.1.4. The variation in thickness
is to be optimized in order to give the greatest edge displacement
at buckling for a given volume of material.
In the case of the simply-supported edge, Figure (3.69 a)
shows that the optimum variation corresponds to
a edge thickness
center thickness
RT = .15, with an increase in buckling displacement
170
0.80 0.90
0.70
O
9
a.
C3
“3
c13.00 0140 0160 12120 1160 2100 572140
a) Simply supported
2.00
1-75
1
.50
fl
ID
‘1‘
0
°.
'b.00 0140 0160 1120 1160 2100 2140
RT
b) Clamped
Figure 3.69. Variation of critical displacement versus RT.
171
0.65 0.70
0.60
D
0.55
.00 4100 0100 .00 16.00 20.00 21.00
cp.60
a) Simply supported
1-40 1.60
.20
D...‘
1.00
.00 4100 6100 .00 16.00 20.00 21.00
cp.60
12
RT
b) Clamped
Figure 3.70. Variation of critical displacement versus RT.
172
of '8119igiiél6 = 33% with respect to uniform thickness plate
of the same volume.
Figure (3.69 b) illustrates the variation of critical displace-
ment with respect to thickness variation for a clamped square plate.
The maximum critical displacement can be achieved if the edge thickness
1.589-1.569
1.569
over the critical edge displacement for the plate of uniform thickness.
is .85 times the center thickness with an increase of a 1.2%
To check the buckling behavior of the plate, within a wider
range of variation in thickness, the critical displacement for
different RT values up to 20 is obtained and plotted in Figure (3J0 );
the plot shows a decrease in critical displacement all the way up
‘to RT = 20,
Note: Examination of mode shapes shows that as the center of the
plate gets thinner, beyond some point, the buckling mode associated
with the lowest critical displacement is not a single concave shape. The
Ibuckling mode corresponding to the lowest eigenvalue for a clamped plate
with RT > 2.2 consists of more than one buckled wave. These results are
not shown in this thesis.
3.2.2.3 Analysis of Buckling_Modes
In this section, the shape of buckling modes will be reviewed
and compared with the exact shape in those cases where the exact
solutionsare available.
Up to this point, all solutions were based on the assumption
of symmetry about both axes of the plate. Thus, only a quarter of the
plate was considered.
173
Therefore, the nonsymmetric modes are missing. To obtain
all modes of buckling, a solution is obtained by considering every
node as an independent degree of freedom, from the out-of-plane displace-
ment point of view, while symmetry and anti-symmetry in u and v
are assumed, as before. Solution for the simply-supported plate
of uniform thickness shows the first mode shape to consist of half-sine
waves in both directions, and the second mode to consist of two
half-sine waves in one direction and one in the perpendicular direction.
The first few mode shapes for the simply-supported plate
are shown in Figure (3.71) and for clamped edges in Figure (3.72) .
174
+
+ + -
0* - 619 U* - 1 51 0* - 1 51
1 ' ' 2 ‘ ' 3 ‘ °
+
+ —
+ — +
- +
0* - 2 458 0* - 2 75 * - 2 95
5 - . 6 - . U6 - .
Figure 3.71. Buckling modes of simply supported square plate, RT= 1.
0* - 1 56 * 2 58 * 2 58
1 - . 02 — . U3 - .
.1.
+
+ —
- + +
0* - 3 55 0* - 3 94 * - 4 3o
4 — . 5 - . 06 - .
Figure 3.72. Buckling modes of clamped square plate, RT = 1.
175
3.2.3 POSTBUCKLING
In this section, a study is made of the postbuckling behavior
of uniform and variable thickness plates under uniform bi—axial edge
displacement compressed beyond the critical displacement. Convergence
of the solution is examined by solving the problem with different
mesh sizes. Variation of in-plane forces as well as in-plane displace—
ments is also studied.
3.2.3.1 Uniform Thickness Plate, s-s Boundaries
A simply-supported uniform thickness plate is compressed
beyond the critical displacement. and the solution is obtained.
Following are some results from the solution.
8)
anvergence of the solution
The problem is solved,successively taking
h/a = 1/4, 1/8, 1/12 and 1/16, and the results
are compared. Table (3;E3) shows the central deflection
of the plate due to edge displacements of 1.26 E;—, which
is almost two times the critical displacement.
Table 3.13 .
176
Convergence of postbuckling solution with mesh size
Mesh size IN-center Difference % Extrapolation N:-center Mgécenter
(h/a)
%- .96795 -118.99 100.80
9.7 .853016
%- .88175 -13o.17 98.55
1.5 .8643366
%5- .86869 -132.98 98.27
.1 .8701166
%g" .86976 -136.74 98.88
Review of Table (3.13 ) indicates that the solution obtained
by only 8 x 8 nodes is satisfactorily close to the converged solution,
keeping in mind that the accuracy of the iterative solution is set
to be one percent.
is small.
Figure 3.73.
O
O
.11
0.96
v
l
l
dues
I
4.00 0100
12.00 10.00
The results can be improved by extrapolation.
a/h
The difference for mesh sizes finer than 8 x 8
Convergence of center-deflection, s-s square plate.
(U = 2 U
c
r)
177
.’—--~ —_ Nx
_God \\\ ----NY
“ \
g; \\
'2' ‘ \x ..
E ‘ '5
.1 \ a!
ELJN. \ 12
u: \ g;
a: \ c1
8 \ =
In . \ <3
:5. \ E
0) In-.. ------ \ g
g—zo. “~- 1. A
a." \”“-~--h .4 -
\
if x
. x
\
\
“~11;
‘—12 r T’ 1 I
0.25 0.50
(X/OI
0) NX RNO NY RLONG RXIS Y=0/2
90.
70 ..
b' '3
a) 'I .6
a: a:
E g
“' SO 3 r-
00 j 40
.1 .1
a: 4:
0. a.
0-0 _. H
c0 c0
12 12
H H
a: 4:
°- 00. ‘-
2A
10
I f I o :25 I I I 0 .ko
(X/O)
c1 Pmucmu. STRESS 610110 111110 1:012
.___.fix
12. ""’"Y
03
J
I
0.50
I I I I I r
0 .26
(X/al
5) SENDING HOHENTS RLONG RXIS Y=°/2
so.
10.
00.
-I
00..
d
10
0.0 '110 1 2'0 ' 030 ' 470
0
019911401991. STRESS v0. U
NOTRTION FOR CURVE IDENTIFICRTION. R) UEI.O IUNDEFLECTED) S) 031.20 C) 022.40 .
FIGURE 3.74. PLOTS OF STRESS COHPONENTS.SDURRE PLRTE.SIHPLY SUPPORTEO.RT=1
b)
e)
178
Stress Analysis
The distribution of the membrane forces will change
as the lateral deflection becomes larger.
Also,the bending moments will vary with variation of the
deflected shape of the plate.
Graphs of IFigures (3.74 a) and (3.74 b) show the variation of
center-line membrane forces and bending moments as the
edge compression varies. Study of these graphs leads
to the following conclusions.
1) As the center of the plate deflects transversely,
ii)
iii)
iv)
the in-plane load shifts with more of the load
being carried by the portion of the plate near the
edges.
The bending moment, which is maximum at the center,
is increasing with the increase in deflection w.
Along the centerline parallel to x, the membrane
force Ny is increasing toward the edge, while Nx
is almost constant.
Bending moments Mx and My both are maximum
at the center and dimimish to zero along the edges,
as expected.
Infplane displacements
The distribution of in-plane displacements in the plate
is shown in Figure.(3. 75) at different stages of loading.
179
Analysis of these leads to the following conclusions:
1) Before buckling, the plate is perfectly flat.
Since there is no w effect, in-plane displacements
are linear, as expected.
ii) After buckling, the plate will undergo more contraction
near the edges and the central region carries less
compressive load; less inrplane displacement
occurs there.
iii) As the load increases this phenomenon becomes more
visible, so that at large edge displacements, the
curves show very little u and v displacement
in the central regions.
.625 .468 .312 .156 0 J5 50 .25 .1. 0 1.5 1. .5 .25 0
a) Undeflected (U/Ucr=1'0) b) Buckled (u/Ucr=1'20) c) Buckled (U/Ucr=2.40)
Figure 3.75. In-plane displacement (U = ua/ti), square plate,
simply-supported, RT = 1.
180
d) simply—supportedJ square_plate, loaded by uniform
normal pressure
To check the reliability of the method by comparing with
previous results, a problem similar to one presented by
Levy is considered. In Reference (29), Levy has solved
a square plate under a uniform normal pressure, with zero
in-plane displacement along the edges. iLevy applies the large
deflection equations and uses the series expansion method.
In this section, a solution is obtained for the same
problem to compare with Levy's. The plot of (3.76 )
shows the variation of w with lateral load as deter-
mined by Levy and by the difference method.
c:
C)
‘¢_
Present result
c:
c3
c6“ A Levy's result
:3
c:
2.3“,-
c:
c?
c:
c?
I
°0.00 10.00 20.020 30.00 70.00
x10
4
_£Ei__
D t
o 1
Figure 3.76. Square plate under uniform lateral load;no in-plane
displacements on boundary;v = .316.
181
3.2.3.2 Simply Supported Square Plate of Variable Thickness
To study the variable thickness plate, a solution was obtained
edge thickness
center thickness
for two different cases. In the first, RT = = %3
for which the membrane stiffness at the edge is l-that at the center,
4
while the flexural stiffness at the edge is %z-of the center stiffness.
In the second, RT - 2, with the edge flexural stiffness being-g-of
the central stiffness. These two opposite variations are chosen
so that the results, along with those for the uniform thickness
plate, would give an idea about the effects of variation in thick-
ness.‘ The discussion follows.
Figure (3.77 ) shows the central deflection of the
plate with variable thickness and also the result for a uniform
thickness plate. Comparison of the three plots leads to the conclusion
that, corresponding to the same edge compression, less deflection
occurs in the plate with RT 8 %-and the plate with RT - 2 under-
goes larger deflection. It should be noted that all these plates
contain the same volume.
This behavior is expected because, in a simply-supported
plate, more bending is occurring in the central region. Thus, plates
with thicker central region will experience less deflection and
plates with thin central regions are more likely to have greater
curvature and deflect large amounts.
182
2.00
1020 1060
Ran/4
0.80
0.40
—l
.00 0'.50 14.00 {.50 2.00
cp.00
Figure 3.77. Central deflection versus edge displacement for different
a)
b)
RT values, simply supported square plate.
Solution Procedure
In (3.2.3.1 a) it was shown that a grid spacing of
h =‘%g-will be accurate enough for engineering design
use. Both problems were solved with a grid spacing
of h =-%g- in determining the results plotted.
Stress Analysis
Distribution of in-plane forces, bending moments and
principal stress,along axes of the plate, are shown
1
in Figure (3.78 ) for RT = 2"
The Graph of Nx and Ny shows a decrease in membrane
forces with an increase in edge displacement in the central
183
I" ....... ux fix
42. x “x “““ "Y :0.
*A
'2 ‘ .. ‘
C’ h:
.— no
.1 {'3
033-24- 2 20d
01 01
K s
33 1:
U .1 o -
a. ‘z
.— H
°° s
u: u1
E"‘~ ‘9 10..
.4
“r
‘2
H A
-8 I 1 I t 1 r 1 i 0 1 1 l l I I I 4
0.28 0.50 0.28 0.80
(X/O) (Kid)
0) NX HNO NY HLGNG 9X18 Y=alé bl SENDING flOflENTS HLONO HXIS Y=G/2
ISO- 100.
120.. 120..
a: ' a: 7
a: in
01 In
a. 1m
'5 00.4 s; .0.-
4 .4
c: is
9.: 2:
c1 ' La ‘
.2. E.
m: a:
a. ‘0. 0. ‘0.-
0 0 a
' ' ' 0125 ' ' ' 0.h.0 .0 l 1'.0 ' ¢.'0 ' sfa ' 411
(11101 F
c) PRINCIPSL STRESS RUONO flXlS Y=°/2 leRINCIFHL STRESS VS. U
NOTRTION FOR CURVE IOENTIFICRTION. B) 021.0 (UNOEFLECTEO) B) UEI.ZS C) 052.75 .
FIGURE 3.78. PLOTS OF STRESS COHPONENTS.SOURRE PLHTE.SIHFLY SUPPORTEO.RT=1/4.
184
region, while in a region close to the edge, Ny is large.
As discussed before, the central region will carry less
load because of the transverse deflection and the load
will be shifted toward edges.
The moments are maximum at the center and vanish at
the edge, but fairly large moments occur at almost
halfway from center to edges. This is due to the relatively
larger thickness at center. Because most of the bending will
occur in the outer region, creating large curvature and
resulting in large moments.
Since we are concerned with the state of stress within
the plate and not necessarily membrane force or bending
moments individually, in Figure (3.78c) , the variation
of principal stress along the axis of plate is plotted
The curve shows maximum stress occuring approximately at a
point x =-%- on this axis. It should be noted that this
stress is also the absolute maximum for the entire plate.
Figure (3.79 ) show the same variables for a plate
with RT = 2.
In this case. because of the thicker edge region, beyond
buckling, the membrane load is sharply shifted toward
the edges.
185
‘----‘\ ——NX
-SO.¢ \ -----NY
\\
A \
*2 J \\
- \
.— \
Ei’i \ '5
_1 \ 0‘)
s x g
0.1 . \ 1.1
K M.-.‘ \ g
0) ‘ ~‘s \ a
33'2“ “x \ a
0: ‘~‘ \ 2
B “’5- “‘ \ 3
1.1 “ ---~.‘"--~. “Q E
z x --~-.-'=‘§ m
a: m
.1 x‘
*‘52‘ ‘
E.
-‘ I r r I I I I I f r I r j
0.25 0.50 0.25 0.50
(X/O) (X/O)
o) NX RNO NY RLONG 8X13 Y:°/2 b) SENDING HOHENTS RLONG flXIS Y=Ol2
90.4 c 90.
J ..
S 7°~ g 10..
a: 00
a: o:
:2 + s -
I- I—
a: 00
.1 50.1 .1 00..
a: a:
4‘: * e:
u u
s. . s -
a: a:
a. 0..
30... 30.
A
V __// d
10 10
' I ' 0128 f ' ' 0.k0 0.0 ' 1'.0 ' 210 ' 3.11 ' 4.0
(X10) U
c) PRINCIPRL STRESS SLONG RXIS Y=°/2 leRINCIFflL STRESS VS. U
NOTFITION FOR CURVE IDENTIFICRTION. fl) 17:1.0 (UNOEFLECTEO) S) 0:1.25 C) 17:2.20 .
FIGURE 3.79. PLOTS OF STRESS COHPONENTsoSOURRE PLRTE.STHPLY SUPPORTEO.RT=2 .
186
The bending moments have a completely different pattern
than for the case RT =-%, because the more flexible
central region results in maximum curvature and maximum
moments in this central region. Also, because of smaller
thickness and smaller flexural rigidity, the maximum
principal stress is always at the center point. In
Figure (3.80 )shows the variation of the maximum principal
stress with postbuckling edge compression for the
different variation in thickness. It can
be seen that the least stress occurs for the uniform
thickness plate,and the plate with RT a-l- is subject to
4
largest stress at points away from center.
c)Inep1ane Displacements
Contours of the in-plane displacement u, in Figure (3.81 )
for RT 8 %-, and Figure (3.82 ) for RT - 2., show
a decreasing displacement in the central region due to
w deflection and corresponding decrease in membrane
forces. In the case of RT- %, we observe very small
displacements in the central region, while the more
closely-spaced contours in the vicinity of the edge
indicate very large strain in this region.
Comparison with the uniform thickness plate shows it
to be between these two variations, as expected.
187
100.00 125.00
78.00
PRINC IPHL STRESS 0'
Q
Do .00 0'.so 1100 1'.50 21.00
Figure 3.80. Principal stress versus edge displacement for different
RT values, simply-supported square plate.
188
1/5/25/
a) Undeflected (U/Ucr-l.) b) Buckled (U/Ucr-1.SO) c) Buckled (U/Ucr-2.50)
W
O 2. .50 .25 O
.8 .6 .4 .2 0
Figure 3.81. In-plane displacement (U=ua/t:), square plate, simply-supported,
1
RT 2'.
.54 .4 .27 .13 0 .7 .6 A1 2. 0 l. .75 .50 .25 0
a) Undeflected (U/Ucr= 1.) b) Buckled (U/Ucr-1.30) c) Buckled (UfUcr-1.85)
Figure 3.82. In-plane displacement (U - ua/ti), square plate, simply-
supported, RT = 2.
189
3.2.3.3 Uniform Thickness Plate, Clamped Boundaries
The problem studied in 3.2.3.1 was solved with clamped
boundaries, and the results of the solution are as follows.
a) Convergence of the solution.
Table<314) presents results for different mesh sizes corresponding
to a boundary displacement of (i°:59 = 1.93 times the
critical displacement.)
Table 3.14. Convergence of postbuckling solution with grid spacing:
3153.21.13 w-center Difference Extrapolation 1*? -center H -center
hla % x x
%. 1,9657 , 314.73 295.45
34.4 1.295366
-% 1.4632 _ 277.99 263.23
4.6 1.376613
-%5 1.39826 289.56 262.96
2.2 1.357113
'%5 1.3674 285.00 260.576
Considering that the accuracy test of the iterative solution
was set at 12 in successive trials, the convergence as illustrated
in graph (3. 83) is good.
Plot of central deflection with edge compression is shown
in Figure (3.84).
190
8
“"1
g
:4
3
S
:1
B
“0.00 4100 0'.00 1'2.00 16.00 a/h
Figure 3.83. Convergence of solution.
‘
g. .7) 7'
m 4: 4.?
a .2
:3 «‘5
m-
[D
c?
N‘
2'
D
'9
1.0
0'
a-
D
c?
I
c0.00 0.50 Lino 1'.50 2.00
Figure 3.84. Central deflection versus edge displacement for
different RT values, square plate, clamped edges.
191
b) Stress Analysis
Distribution of membrane forces, bending moments, and
principal stress along the axis of the plate, are plotted
in Figure (3. 85). These figures show:
i)
ii)
iii)
Larger membrane force occurs at the edge and the
central membrane forces decrease with an increase
in postbuckling edge compression.
Bending moment, Mx’ is positive at the center and
along the edge it is negative with an absolute
value larger than the central moment for the larger
deflection.
Principal stress is almost equal at center and edge
in the early postbuckling stage, increasing on the
edge with larger deflection.
c) Ineplane Displacements
Contours of in-plane displacement, u, shown in Figure
(3.86 ) indicate the following.
i)
ii)
111)
Before buckling, in the flat plate, in-plane dis-
placements are linear everywhere as expected.
Immediately after buckling, as the deflection starts to
increase, the u—displacement tends to decrease in the
central region.
The pattern of displacement is qualitatively similar
to that for the simply supported case.
192
’1’, ~‘\\ .___—_Nx
-100‘II \\ ----- NY
A \
\
1: J \
" \
SE \ -
5— w
.J 0‘)
D I-
O) 2
1.1 1.1
K I:
O
0) t
0)
01 c:
0: 2
h- 0.0
0) O
2
01 01
2 111
S
.J
0..
I
2
H
-20 T’ f 1’ I I I
0.25 0.00
(X/a)
0) NX 0N0 NY RLONG RXIS Y:°/2
fi
200.
b b
so no
8 C 8
0: 1:
h- F-
0) a:
A .J
C: e a:
L 0.
H 0..
L.) U
2 2
H H
a 120.. a:
Q. 0.
.1
A
‘0 ' 4T ' 0125 ' ' ' 0.00
(X/O)
C) PRINCIPRL STRESS RLONG RXIS Y=0V2
30.1
-80 .1
b1 SENDING noncurs RLONO 0x13 7:012
ZOO-J
200 .4
150 .1
40
I I I I
0.0 1.0 210
u
leRINCIPRL STRESS vs. U
310 ' 4.0
"(ITRTION FOR CURVE IDENTIFICRTPON. B] 021.0 (UNOEFLECTED) B) 051.20 C) 052.25 .
FIGURE 3.85. PLOTS OF STRESS COflPONENTS.SOUflRE PLHTE.CLRHPED. RT=! .
193
1.56 1.25 .33 .41 0 2. 1. f o 3. 2. 1. .5 o
a) Undeflected (U/chfllj b) Buckled (U/Ucr=l'20) c) Buckled (U/Ucr32.2)
Figure 3.86. In-plane displacement (U = ua/ti), square plate, clamped
edges, RT - l.
194
3.2.3.4 Clamped Plate with Variable Thickness
Again,the plates discussed in Section 3.2.3.2 (RT = 2%, and
RT = 2)are solved with clamped boundaries along all edges.
For comparison,the variation of w with edge compression
is plotted in.Figure 3.84 for plates with RT =-l- and RT = 2,
4
along with uniform thickness plate.
It can be seen that in the case of the clamped plate, due to
edge displacement appreciably greater than critical displacement,
the plate with RT =-% undergoes less deflection than either of the
plates with RT = 1 and RT = 2. This is similar to the case of
simple support, but in the early stages of postbuckling, the plate
with RT = 1 has the smaller lateral displacement.
195
a) Stress Analysis
Figure €3.87) illustrates the distribution of membrane
forces, bending moments, and principal stresses for' the
1
plate with RT = 2:. Analysis leads to the following
conclusions:
1)
ii)
iii)
iv)
While the membrane forces are decreasing with in-
crease in edge displacement at the center, these forces
increase sharply on the edges.
Bending moments are maximum.at the center, and the
negative moments along the edge are small because
of small flexural rigidity.
Principal stress is maximum on the edge and in-
creases with an increase in edge displacement.
Maximum principal stress occurs at the center of
the edge.
Study of membrane forces, bending moments, and principal
stresses for' the plate with RT = 2, in Figure (3.88), results
in the following observations.
1)
Ny decreases in the central region as the edge
displacement increases,while it increases sharply
near the edge.
ii) Nx also decreases at the center and decreases along
the edge.
196
1 ‘1
L
O
i
.1.
O
P
-0500100 000601810)
IN-PLRNE STRESS RESULTRNT
~10
I I I I I r I I
0.28 0.80
(X/Ol
0) NX RNO NY HLONG RXIS Y=0/2 b) SENDINC HOHENTS RLONG RXIS Y=°/2
90°? 8007
4801 180..
‘6 b
a: 7 a: 7
33 33
15 55
a, 300. a, 3003
.1 _J
C C
2: 2:
u " (J I
35 55
z.‘ 1:
“- 180- “- 180..
a
1\
0 A 0
r ' ' 01.28 ' ' ' 0.150 0.0 ' 1r0 ‘ 210 ' 330 r 4.11
(XI 1 U
:1 901001901 STRESS 01000 0x18 7:012 01901001901 STRESS vs. 0
00101100 900 cuave 10001191001100. 01 051.0 1000091501001 01‘Ue1.30 01 Ue2.30 .
FIGURE 3.87. PLOTS OF STRESS COHPONENTS.SOURRE PLRTE.CLRHPED. RT=1/4 .
197'
’ ‘-’---- -
------- x
‘1' s‘ -——N
”
I
O
l
SENDING flOflENTSIfiI
I
a
O
1
'2‘ .1
10-91009 810988 098011001 (0‘1
a W I I I I I I I -‘OJ
0.28 0.80
1x101
01 0x 000 01 01000 01118 1:0/2 111 8900100 11009018 01000 01118 1:0/2
228- 226.
I75... ‘75-1
128 .. 128.1
901001901 810988 10”)
901001901. 810988 (0'1
75.1 75.1
28 I I I I as j I I I I I I
0.28 E80 0.0 1.0 2.0 3.0 470
(X101 - U
:1 901001901 810988 01000 01118 Y=°/2 01901001901 810988 V8. 0
NOTRTION FOR CURVE-IOENTIFICRTION. R) UEI.O (UNOEFLECTEO) B) 051.30 C) U22.40 .
FIGURE 3.88. PLOTS OF STRESS COHPONENTS.SOURRE PLRTE.CLRHPED. RT=2 .
198
iii) Bending moments increase in magnitude at the
center and edges as the edge displacement increases,
as expected. The negative moments at the edges are
much larger, because the flexural rigidity is larger there.
iv) The graph of the principal stress shows an increase everywhere
as the edge displacement increases,but it is alwaysaa
maximum at the center of the plate.
v) Figure (3. 89) shows the variation of maximum
principal stress with edge compression, 0 , for the
three different RT values and for plates of constant
volume.
It can be seen that the plate with RT = %—
always undergoes larger stress. Although the
uniform thickness plate is less highly stressed
than the case RT = 2 for smaller edge displacements,
in the higher range of edge compression it is more highly
stressed than the plate with RT = 2.
NOte: It should be noted thatfor RT = l, and RT =~£,
the maximum.principal stress is located at the center of the
edge while for case of RT = 2, the location is .at the center
of the plate.
199
{x
62.50
J
'810‘
37.60 69.00
J
25.00
J
RT=2
CIPHL STRESS cr
PRIN
250
j
I I I fi
000 1000 ‘ 2000 3000 4000
ALDO
Figure 3.89. Max. Principal stress versus edge displacement for
different RT values, square plate, clamped boundaries.
200
b) Ineplane Displacements
Contours of in-plane displacement, U, are plotted in
Figure (3.90 ) for plate with RT = %-, and in Figure
(3.91 ) for plate of RT = 2. These contours show the
following:
1) Before buckling, the contours are exactly the same as
those for simple support, because in the undeflected
position the out-of-plane boundary condition has no
effect on the solution.
ii) After buckling, as usual, the in-plane displacement
in the central region is smaller compared to undeflected
case (i.e.,the contours are expanding at the center
and compacted contours are located away from center
toward the edge depending on the RT values.)
iii) In the case RT = %- closer contours are located
in the vicinity of the edge, while for the plate
RT = 2, because of stiff edges these concentrated
contours are seen to be close to the center.
iv) The behavior of the uniform thickness plate falls between
the cases RT =-%- and RT = 2.
201
/ W
1.25 1. .5 .25 o 1.5 1. .50 o 3
0
a) Undeflected (U/Ucral.) b) Buckled (U/Ucr=l.35) c) Buckled (U/Ucr=2.30)
Figure 3.90. In-plane displacement (U = ua/ti), square plate, clamped,
RT = 1/4.
1.40 1.05 .70 .35 0 1.50 1. .50 0 3. 2. 1. 0
a) Undeflected (UflUCEI1.) b) Buckled (UflJcr-l.25) c) Buckled (UYUcr-2.15)
Figure 3.91. In-plane displacement (U - ua/ti), square plate, clamped,
RT = 2.
202
3.3 COMPARISON OF TWO METHODS
In order to compare results from the two methods discussed-
formulation in terms of the stress function or in terms of displacements-
two problems are solved by both methods and compared. Their results
can also be used as a measure of the accuracy of either of the methods.
a) Simply Supported
i) A simply-supported uniform stiffness plate, with
ii)
iii)
no restraint on in-plane boundary displacement and
loaded beyond the critical load was solved using the method
discussed in Section (2.2) (in terms of m and w).
The solution includes in-plane displacements, u
and v, on the boundary.
The problem is solved using the method of Section (2.3)
(in terms of three displacements, u, v and w),
by applying ‘the boundary displacements obtained in (i)
as boundary conditions.
If both methods are correct, we expect (1) and
(ii) to result in the same solutions, and.they did
turn out to be very close. The central deflection (w/t)
is .9613 in (i) and .9736 in (ii), witha difference
of 1.3% . This is very good agreement. Other
components of stress and displacement are also
very close.
b)
203
Clamped Edges
A uniform stiffness plate with clamped boundaries was
also examined in the same way as discussed in part (a).
The central deflection (w/t) was found to be 1.67121,
using the method of Section 2.2 (i.e.,formulation in
terms of stress function). It was 1.68578 using the formulation
in terms of displacements, with a difference of only .82.
The results are considered to be very good.
CHAPTER IV
CONCLUSION
4.1 THE PROBLEM SUMMARY
In the preceding chapters, the behavior of variable stiffness
plates was studied in the prebuckling and postbuckling range. In-
stability criteria were also examined. The work used the ordinary
finite difference technique. No results for similar variable stiffness
plates are available in the literature to confirm the validity of the
solutions. However, a uniform stiffness plate was included in each
case to serve as a control problem, and the results obtained by the
difference method were compared with those of analytical solutions
and other published results.-
Since the main purpose of varying stiffness is to optimize
the plate with respect to some design variables, some optimization
examples were presented in the buckling analyses.
In order to provide a better perspective on the change in
behavior of plates due to stiffness variation, the different problems
considered were assumed to contain a constant mean stiffness
or a constant amount of material (see Section 3.1.4). Two different
approaches were discussed:
1 Formulating in terms of stress function, Q, and the lateral
deflection, w.
204
205
2. Formulating in terms of three displacement components,
u, v and w.
First the applicable difference operators were worked out; then
the solution procedures were presented in detail. Finally, an
increasing number of nodes were used to check the convergence of
the solutions and also to provide guidance in selecting an
appropriate mesh size to obtain the desired accuracy.
The behavior of different stress and displacement components
was illustrated in suitable graphs and the results were analyzed.
In Section (3.1.1), uniform thickness plates with different
variations in E were considered,and the behavior of in-plane forces
and displacements was analyzed in Section (3.1.1.1). The buckling
of those plates was considered in Section (3.1.2.1) and the effect of
variation in stiffness on stability criteria was discussed in Section
(3.1.2.2),In Section (3.1.3), postbuckling behavior of those plates
‘was examined and the results were analyzed in Section (3.1.3.3).
Sections (3.1.4) and (3.2.2.2) deal with optimization of
variation in thickness from the stability point of view. The remainder
of Section 3.2. shows the effect of variation in thickness (with
constant B) on displacement, forces and moments as well as buckling
behavior.
206
4.2 CONCLUSION
The comparison of problems solved in Section 3.1 (with free
in-plane movement on the boundaries) with those of Section 3.2
(restricted in-plane displacement on the boundaries), shows the
great effect of in-plane boundary restraints in the postbuckling state.
Study of figures (3.60), (3.77) and (3.84), leads to
the conclusion that the behavior of central deflection due to vari-
ation in stiffness is not only dependent on the out of plane boundary
conditions but also greatly affected by in-plane boundary conditions.
Figures (3.60), (3.77) and (3.84) show that the trend of
central deflection of plates with different variation in stiffness
is not the same for all postbuckling ranges. For example, in Figure
(3.77), we observe that corresponding to the same load, a simply
supported plate with RT = l undergoes larger deflection than a plate
with RT =-% . Study of Figure (3.84) shows that the same plates
with clamped boundaries exhibit different behavior. (Although, for highly
compressed edges, less deflection is observed for the plate with
RT = %3 when the edge compression is only slightly above critical
displacement, the larger deflection corresponds to plate with RT = %).
The method and corresponding computer program utilized for
force boundary condition (Section 2.2) has been shown to give more
accurate results than that developed for the displacement boundary
condition (Section 2.3). This is due to the difference in the order of
derivatives involved in the formulation. In the force boundary
condition formulation, only second and higher order derivatives
of the two functions, 0 and w, are involved in the equilibrium and
207
compatibility equations (equation (2.14) and (2.15). The displace-
ment formulation, however, involves first and higher order derivatives
of the three functions, u, v, and w. Hence accumulative errors
could be expectedfiln difference approximation equations (2.19), the first
error term in the first order derivative includes the third derivative
of the function, and the error term in the second derivative includes
the fourth derivative of the function, etc).
Investigation of convergence indicates that for engineering
purposes, grid spacings for h =-% in the force formulation (see
Figure 3.39) and h =-%§ in the displacement formulation (see Table
3.13) are reasonably accurate; more accurate results can be obtained
by using finer grids and applying Richardson's extrapolation to the
results. Comparison of the results with known values supports the
reliability of the solutions.
208
4.3 RECOMMENDATIONS
The objective of this work was to examine application of the
two formulations in general, and to investigate the behavior of plates
with different boundary conditions and stiffness variations. Re-
finements and extensions of the method which are possible and de-
sirable include:
1.
Application of improved difference methods involving
more accurate approximations to the derivatives.
Inclusion of more nonlinear terms in the strain components
and a study of their effect on the results.
The difference operators can be revised to make them
applicable to orthotropic plates, and the computer
program improved so that it can be applicable to orthotropic
and nonhomogeneous materials.
Due to the absence of experimental sources to guarantee
the accuracy and practicability of the results,and re-
cognizing the advantages in the use of variable stiffness
plates, an experimental study of such plates from the
stability point of view and in the postbuckling range
should be very useful.
The computer programs developed here were mainly aimed
to solve particular problems. Although they are
more general than needed for the problems solved here,for
applications to loading and geometry different from the
ones presented here, the programs should be used with
caution and appropriate changes made. Also, the efficiency
of the programs can be improved.
209
Application of other methods such as the finite element
method using:
a) Elements with variable stiffness within the element.
b) Constant stiffness within an element but variation
of stiffness from element to element,
The boundary integral method might also be considered.
Consideration of the same cases and comparing numerical
results, convergence, computer cost, etc., would be of
interest.
B IBLOGRAPHY
10.
11.
12.
BIBLIOGRAPHY
Aalami, B., and Chapman, J.C., "Large Deflection Behavior of
Rectangular Orthotropic Plates and Transverse and In-plane
Loads",Institution of Civil Eng. Proceedings 42, pp 347—382
(Mar. 1969).
Allen, D.N. De G., and Windle, D.N., "The Finite Difference
Approach",Stress Analysis,edited by O.C. Zienkiewicz and
G. S. Holister, John Wiley and Sons,(1965).
Altiero, N.J., and Sikarskie, D.L.,"A Boundary Integral Method
Applied to Plates of Arbitrary Plan Form",Computer and
Structures (1979). '
Argyris, J.H.,"Energy Theorems and Structural Analysis",Aircraft
Eng. 26, 1954, 27,(1955).
Basu, A.K., and Chapman, J.C., "Large Deflection Behavior of
Transversely Loaded Rectangular Plates",Proc. Inst. of
Civil Eng. 35, p 79 (1966).
Bathe, K.J., Oden, J.T. and Wunderlich, W., "Formulation and
Computational Algorithms in Finite Element Analysis",
MIT Press (1977).
Billington,D.P., "Thin Shell Concrete Structures", McGraw-Hill (1965).
Bradley, W.A., "An Alternate Finite Difference Operator for
Plates of Variable Stiffness",Int'l Journal for Numerical
Methods in Eng. 11, 6, 1044-1048 (1977).
Brown, J.C. , and Harvey, J.M., "Large Deflection of Rectangular
Plates Subjected to Uniform Lateral Pressure and Compressive
Edge Loading", Journal of Mech. Eng. Science, vol. 11,
No. 3, (1969).
Bryan, G.H. "On the Stability of Plane Plates Under Thrust
in it's Own Plane with Application to the Buckling of the
Side of the Ship",Proc. London Math. Soc., vol. 22, pp 54-67,
(1891).
Clough, R.W., and Felippa, C.A., "A Uniform Quadrilateral Element
for Analysis of Plate Bending", Proceeding of the Second
Conference on Matrix Methods in Structural Mechanics,
AFFDL-TR-68-150,(Dec. 1969) pp 399-440.
Coan, J.M., "Large Deflection Theory for Plates with Small
Initial Curvature Loaded in Edge Compression", Journal
of App. Mech. Tran. ASME, vol. 18 series E, pp 143-150,
(1951).
210
211
13. Colville, J., Beker, E.B., and Eurlang, R.W., "Large Displacement
Analysis of Thin Plates", J. Struct. Div., ASCE, vol 59,
No.ST3,pp 349-364,(Mar. 1973).
14. Conner, J., Jr., Logcher, R.D., and Chan, S.C., "Nonlinear
Analysis of Elastic Framed Structures", J. Struct. Div.,
ASCE, vol. 94, No. ST6, Proc. Paper 6011, pp 1525-1547
(June 1968).
15. Dawe, D.J.,"Application of the Discrete Element Method to the
Buckling Analysis of Rectangular Plates under Arbitrary
Membrance Loading", Aeronautics Quarterly, vol. 20, pp
114-128 (1969).
16. Dombourian, E.M., Smith Jr., C.V., and Carlson, R.L., "A
Perturbation Solution to a Plate Postbuckling Problem",
Int. J. Nonlinear Mech.,vol. 11, pp 49-58 (1976).
17. Euler, L.,"De Motu Vibratorio Tympandrum", Novi Commentari
Acad. Petropolit., 10,pp 243-260 (1766).
18. Galerkin, B.C., "Series Solution of some cases of Equilibrium
of Elastic Beams and P1ates§(in Russian), Vestnik, Inzhenernov.
1, pp 879—903 (1915).
19. Gallagher, R.H.,"Geometrically Nonlinear Finite Element Analysis",
Proc. Specialty Conference on Finite Element Method in
Civil Eng.,McGill University, Montreal, Canada; pp l-33(1972).
20. Gallagher,R.H., Lien, 8., and Mau, S.T.,"Finite Element Plate
and Shell Pre-and Postbuckling Analysis",Proceedings,
3rd Air Force Conference on Matrix Methods of Structural
Mechanics, Dayton, Ohio, (Oct. 1971) pp 857-879.
21. Gass, N., and Taborrok, B.,"Large Deformation Analysis of
Plates and Cylinderical Shells by a Mixed Finite Element
Method", Int. J. for Num. Methods in Eng., vol. 10,pp 731-746
(1976).
22. Germain, Sophie, "Ibcherches Sur La Theorie des Surfaces
elastiques",Paris (1821).
23. Green, A.E., "On Reissner's Theory of Bending of Elastic Plates”,
Appl. Math., 7, pp 223-228 (1949).
24. Hartz, B.J., "Matrix Formulations of Structural Stability
Problems", J. Struct. Div., Proc. ASCE,vol. 91, No ST6,
pp 141—157 (1965).
25. Jaswon, M.A. and Maiti, M., "An Integral Equation Formulation
of Plate Bending Problems",J. Eng. Mech., II (1968).
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
212
Kaiser, R., "Theoretical and Experimental Determinations of
Deflections and Stresses in Square Plate Freely Supported
at the Edges, with Uniform Load Distributions and Large
Deflection",Z. Angew. Math. Mech. 16, 73 (1936).
Kennedy, J.B., and Prabhakara, M.K., "Postbuckling of Orthotropic
Skew Plate Structures", J. Struct.Div., ASCE,pp 1497-1511
(July 1980).
Lagrange, J. Ann. Chim. vol. 39 pp 149 (1828).
Levy, S.,"Bending of Rectangular Plates with Large Deflections".
NACA, TR737,pp 139-157 (1942).
Maiti, M. and Chakrabarty, S.K.,"Integral Equations for Simply
Supported Polygonal Plates", Int. J. Mech. Sci., 12 (1974).
Murray, D.W. and Wilson, E.,"Finite Element Postbuckling
Analysis of Thin Elastic Plates", AIAA J.,vol. 7, No. 10,
pp 1915-1920 (1969).
Murthy, S. D.N., and Sherbourne, A.N., "Nonlinear Bending of
Elastic plates of Variable Profile", J. Eng. Mech. Div.,
Proc. ASCE, 100, EM2, 251-165 (Apr. 1974).
Neilsen, N.J., Bestemmelse af Spaendinger i plader ved Anvendelse
af Differensligninger, Copenhagen (1920).
Oden, J.T., "Numerical Formulation of Nonlinear Elasticity
Problems", J. Struct. Div., ASCE, vol 93, No. ST3 pp 235-255
(June 1967).
Otter, J.R.H., "Computations for Prestressed Concrete Reactor
Pressure Vessels using Dynamic Relaxation", Nuclear Structural
Eng., 1, 61 (1965).
Prabhakara, M.K.,and Chia,C.Y., "Postbuckling Behavior of
Rectangular Orthotropic Plates", J. Mech. Eng. Sci.,vol.
15, No.1, pp 25-33 (1973).
Ramberg, W., McPherson, A.E., and Levy, 8., "Normal Pressure
Tests of Rectangular Plates", NACA TN748,pp 349-371 (1942).
Rushton, K.R., "Dynamic Relexation Solutions of Elastic Plate
Problems", J. Strain Analysis,pp 3-23 (1968).
Rushton, K.R., "Large Deflection of Variable Thickness Plates",
Int. J. Mech. Sci.,vol 10, pp 723-735 (1968).
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
213
Rushton, K.R., "The Dynamic Relaxation Method used for Stress
Analysis", Conf. Recent Adv. Stress Analysis, Royal Aeronaut.
Soc. 3, (1968).
Shye, K.Y., and Colville, J., "Postbuckling Finite Element
Analysis of Flat Plates", J. Struc. Div., ASCE, vol. 105,
No. 8T2, pp 297-311 (Feb. 1979).
Supple, W.J., "Change of Wave Form of Plates in the Postbuckling
Range", Int. J. Solids Structures 6, 1243 (1970).
Szilard, R., "Theory and Analysis of Plates", Prentice-Hall
(1974).
Timoshenko and Gere, "Theory of Elastic Stability", 2nd edition,
McGraw-Hill (1961).
Turner, M.J., Clough, R.W., Martin, G.C. and Topp, L.J.,
"Stiffness and Deflections Analysis of Complex Structures",
J. Aero. Sci., 23 pp 805-823 (Sept. 1956).
Vlasov, V.Z., "Some New Problems on Shells and Thin Structures",
NACA, Tech. Memo. No. 1204 (Mar. 1949).
Von Karman, T.,"Festigkeitsprobleme hm Maschinenbau", vol.
IV,Pt. 4, of Encyk. der Math. Wiss., 1910, Art. 27, pp
311-388, ASME Trans. Apm 54-5, vol. 54, No.2 pp 53-57
(Jan. 1932).
V08, G.R., "Finite Element Analysis of Plate Buckling and
Postbuckling", Ph.D. Thesis, Rice University (1971).
Wagner, H.,"Ebene Blechwandtrager mit sehr Dfinnem Stegblech",
Z. Flugtech. Motorluftschiffahrt vol. 20, pp 200-306 (1929).
Wu, B.C., "A New Method for the Solution of Thin Plate Problems",
Ph.D. Thesis,Michigan State University (1980).
Yamaki, N.,"Postbuckling of Rectangular Plates with Small
Initial Curvature Loaded in Edge Compression", Trans.
ASME J. App. Mech., vol. 26,pp 407-414 (1959).
Yang, T.Y., "A Finite Element Procedure for Large Deflection
Analysis of Plates with Initial Deflections", AIAA Journal,
vol. 9, No. 8, pp 1468-1473 (Aug. 1971).
Zienkiewicz, O.C., "Finite Element Procedures in the Solution
of Plate and Shell Problems", Stress Analysis,edited by
O.C. Zienkiewicz, and G.S. Holister, John Wiley and Sons (1965).
APPENDICES
APPENDIX A
APPENDIX A
Finite Difference Operators
214
215
.a . a .Aea.~. 860888.. aaauaaaaauo so .80. ease aeoa
you hauouoeo mucouomwan
Ame+m8va+ Haso+oeciwuac+ Aze+meva+
.8.»-..N .m.+....5+~..~1 s..>-a.~
e .1 infant—3.3.2-3? o a
2 9. 38.? z m z u 2 9. 3.10..
m . 8+ 8+ 0+ 9+ 3 .
A 8+ s...+~..~1 .a+~..ee A 8+ o..:+a..~1
.we+zs.>+ H.88+u...>-av+ Azo+38.>+
.23: .22....308- 6. .3:
z
“I
216
one: a . .
uuwu Because H I s .AeH NV cowunsuo we 01«J.vann awed you nousuono mucououuua .N < shaman
.oucda
cu
”I
217
.H I s .AmH.~V cowussco >ufiaanfiuoeaoo we «can can: uuos now unusuoeo menopauuav waded» .n< ouswuh
AAzrmeAaruswam +
Au To a... 9... $I~I +
9:
Amrurzessm ta
«A2 Tm 31'.“ +
A3809 mm .. Azanaumeu
A
As... a-.. 8 3A2 a... e ~+. 8..
z m s u .31
2 9 a: 9 s 95
84.4.2648.
Az em 1: a... 8% +. a
«Azrmsnfi .. $8.808...
.a
A281... 3 9A
N
NAsses Mm, -Az¥.a-me$
. .1 6 .181
Annulm+mv>+u+
A21... 9.353% tam
«AssumaTmN. 1A..$.a~..m 8..
A3... a-.. 8 1303?. 8-
«Asrmshw i=8. 0...... A2 a... 8-...eAaéaJA31881M - "AaruaTalm 138.808.
Anaheim..- A=$.s-ueA10N+.8~ Anaheim..-
z m s u .3 Azeums.3+ m a u .
:1 a: a. seam- a a :21 8:... 84.1%..
0 a
A a. fiefiaeha +
A2121... fuem +. a
.H I a .Am~.~v nouuoaco no sowusaauouees ousououuav ouweum
AL.”
218
.4. 60.6.0
A
iul-lllll.
a + AA. 8.. we
Alia-cm Tm + .4195? a; N +
u a m . . 3. m
m u . . a.
4.0 A-m+1m:10- . Im-fl.
U 0 n— 6 o
u o m m a a m
a a .. Al-+-l+11+-I.A9..3~+.A-IJ1- u a
u z m . a . a .
LL- LLLL .900..-
.AHI. AVA>S~ .. a+ A... H+H AH+ AVA>d~
. 0
v a
Q v 0 +.
in u + A-wTIMSAidT 5: N
m z z . 3 2
ii..- 1..-«:38- 1-1.-
HI..."
APPENDIX B
APPENDIX B
The plan and node arrangement on the portion of plate considered
in each case is shown. It should be noted that the second rows of
exterior nodes are auxiliary nodes for defining the K vector at edge
nodes only, and those nodes do not participate in any calculations. Thus,
'tlie node number for them could be any number or repetition of previous
Ones.
219
220
22
D
\
I \
22I \
[— "2'F\. .
| l \ ’
1 I \
221-__2