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