14-2-2011.pdf 156 The primary objective is to determine whether the structural efficiency of plates can be improved with variable thickness. The large displacement analysis of steel plate with variable thickness at - direction is obtained numerically, using finite differences. The effects of boundary condition, tapering ratio, type of tapering equation and plate aspect ratio on large deflection behavior of rectangular plates are investigated. Numerical results for rectangular steel plate are presented for the different effects. This study showed that the large deflection behavior is very sensitive for thickness variation (tapering ratio) where the maximum deflection will increase about 5% for slenderness ratio ( / =100) and tapering ratio (1-2) of simply supported plate. : ـل .فعالیة لتحلی . ل . ، نسبة، . ) ()5% ()( / =100 )٢-١ ( . 157 Plate dimension in and - directions respectively. = clamped edge. =Slope coefficient of the Tapered plate. 23 112 Modulus of Rigidity. Modulus of Elasticity. , , In-plane forces (per unit width). Transverse load (per unit area). = Simple supported edge. Plate thickness. Thickness at the side = . = Average thickness (( + )/2). Thickness at the side =0. Poisson’s ratio. Out-of-plane displacement (or deflection). With the increasing use of rolled, machined, or them-milled skins in aircraft and missile designs, analyses of plates tapered in thickness are becoming the rule rather than the exception. Design data for the large displacement of such plates are limited(2). Nonlinear structural problems usually fall into one of the following main categories: (a) large deflections; (b) finite strains; (c) nonlinear material properties; (d) deformation-dependent interactions between structural parts; (e) combination of nonlinear material behavior with one of the other categories. The problem of geometrical nonlinearity is of considerable practical interest for aerospace engineers and naval architects. In civil engineering, hanging roofs, suspension bridges, etc. constitute the most important class of structures which display pronounced geometrical nonlinearities(2). The function of a thin plate element is generally to withstand a distributed lateral pressure, or to act with the adjoining structure in sustaining in-plane forces, or both. Although the equations of the large deflection behavior of plates were first derived by Von Karman, it is only through recent advances in the development of numerical methods that the general problem of plates has been treated satisfactorily. The early investigators used infinite double Fourier series (Levy[12]) and finite differences (Basu and Chapman[3]). They studied the large deflection behavior of thin plate under uniformly distributed loading. Aalami and Chapman[1] used the finite difference method to obtain solutions for a number of isolated plates under uniformly distributed loading and with simple boundary conditions. Ueda, et al.[4] studied the large deflection behavior of a rectangular plate by an efficient semi-analytic method. An incremental form of the governing differential equations of plates and stiffened plates with initial deflection had been derived. For each load increment, these equations were solved by Galerkin method with special consideration of simple supported boundaries. Recently, Jayachandran, et al.[5] derived incremental matrices for thin initially imperfect plates with a small out-of-flatness by using minimum potential energy principles. Explicit coefficients of the displacement gradient tensor had been evaluated. These matrices were used in combination with any thin plate element. The formulations were incorporated in software plot-cold. The basic concepts of three methods (or steps) for the solution of the large deflection behavior of plates together with a description of the boundary conditions of the cases examined are summarized in the following. 158 Starting from the equilibrium of compatibility of a thin plate element and expressing the strains and curvatures as functions of the stress resultants, the following equations are presented(2): - The stress-strain relationships become as follows: - where ; in which ; and denote the thickness at the sides 0 and , denote to type of equation. After derivation of strains and substitute the derivations into Equation (1) and express it as a stress function resultant where the compatibility equation becomes: . (3) By similar algebraic steps it is possible to write the equilibrium equation in terms of and Φ, thus: - . (4) Where )( )( (5) 159 And the variable is Airy’s stress function and is defined such that: - ;; (6) All the other quantities related to the bending and membrane actions of the plate can be expressed in terms of two variables ( ) and ( )(2). The boundary conditions considered for = , , as shown in Figure (1), are as follows: - 1-Boundary on rigid supports (7) 2- For rotationally free cases (on edge parallel to -axis) (8) For rotationally fixed cases(on edge parallel to -axis) (9) For membrane action cases of fully free is considered. 3-For fully free condition of zero membrane direct stress (on edge parallel to -axis) (10) Zero membrane shear stress (on edges parallel to or -axis) (11) Equations (10) and (11) refer to the boundary condition of an isolated free plate(2). For solution, the plate is subdivided into a graded mesh. The grading is chosen to be fine under distributed loading and for the fixed boundary conditions for increased accuracy. Equations (3) and (4), together with the boundary conditions are expressed at the nodes in terms of central finite difference expressions, as shown in Figure(2). The resulting equations may then be arranged into the following coupled matrix form: - (12) (13) in which [ ] = a square matrix with constant coefficients depending on ; [ ] = a square matrix depending on ; [ ] = a square matrix for bending effect obtained from the left hand side of Equation (13); [ ( )] = a square matrix depending on ; { } and { } = column matrices of the unknown variables and , respectively; { } = the applied transverse loading. For numerical evaluation, it is convenient to rewrite Equation (10) as follows: - 1 (14) 160 Now for any specified conditions, Equations (12) and (13) can be solved for and using an iterative procedure. There are various schemes used for the solution of the present type of coupled equations. In the present study, the successive iteration procedure is considered. The solution to be obtained for a given applied loading { } is achieved by assigning a value to (which may be assumed equal to zero for first loading, and equal to the previous values for subsequent loading). Thus Equation (12) can be solved for { }, from which [ ] is evaluated, and subsequently used in Equation (13) to get new values of . The procedure is repeated until a desired degree of accuracy is reached. There are certain refinements in the iteration, which are employed to ensure a rapid convergence. A computer program was written by Amash [2] and developed to take the out of plane loading and multi equations of variable thickness. The program is written in FORTRAN 90 language. All the solutions presented were obtained with a specified degree of accuracy ( ) of 0.1% in the iterative procedure. In most cases the plate was divided into 16 divisions in any direction. This was considered adequate for obtaining deflections which are accurate enough for practical purposes. Ideally it would be desirable to compare the theoretical predictions of the program with the results of carefully controlled experiments and other theoretical results. The accuracy of the results of the present study in the analysis of real panels is compared with the theoretical results obtained by Levy (11) [1942] on simply supported panels. The numerical analysis of Levy based on the infinite double Fourier series for the non-linear analysis of general steel-plate. In the present study, this plate is analyzed based on the prescribed procedure and it is divided into (16 16) divisions. Figure (3) shows a comparison between the theoretical results of Levy’s study and the present study for the out-of-plane displacements. The curves shown relate the applied load ( 44 ) on the vertical side to the non-dimensional maximum deflection ( ) on the horizontal side. A similar format is adopted for all other figures. It is clear from this figure that good agreement between the results by the present method and the experimental results and the theoretical results is obtained for a simply supported thin plate. Figure (4) presents the load-deflection curve of a square simply supported plate under transverse uniform load. The effect of a variation in thickness is considered. The values of tapering ratio are taken to be ( / =1.0, 1.5, and 2.0). So, the modulus of elasticity for all plate in this study was taken equal to 200 GPa. This figure shows that the effect of tapering ratio on large deflection behavior is considerable when the applied load is increasing. The increase in the effect of tapering ratio will be more appearing when the effect of membrane action becomes more appearing. Figure (5) presents the load-deflection curve of a square thin plate with all edges clamped. This figure shows a comparison between the constant thickness and the variable thickness in perfect plate. The effect of variation in thickness is very clear on large deflection behavior of thin plates. This effect is appearing in the beginning of loading of plate by transverse uniform load. The effect 161 of material nonlinearity is neglected in the present study. These results are compared for a plate with constant thickness (ordinary plate) and with average thickness which is equivalent thickness of ordinary plate. Figure (6) presents the load deflection curve of a square plate with two edges simply supported and other edges clamped and under transverse uniform load. This figure shows the effect of boundary condition and tapering ratio on the large deflection behavior of a square plate under uniform lateral load. Figures (7) and (8) presents the load-deflection curve of a rectangular simply supported thin plate under uniform lateral load with aspect ratios ( / =2.0 and 0.5), respectively. These figures show the effect of tapering ratio on the large deflection behavior of a rectangular thin plate with increasing and decreasing in the length in -direction (parallel to the variation in thickness). From Figures (7) and (8), the following is noticed: 1- The values of the deflection of the thin plate with aspect ratio ( / =2.0) and under uniform load are less than the values of the deflection of the thin plate with aspect ratio ( / =0.5). This difference is due to the effect of tapering ratio on the plate where the deflection becomes larger for the plate with aspect ratio ( / =0.5). 2- The effect of tapering ratio on large deflection behavior of a thin plate with aspect ratio ( / =0.5) is more than the effect of tapering ratio on large deflection behavior of a thin plate with aspect ratio ( / =2.0). Figures (9) presents the load-deflection curve of a rectangular simply supported thin plate under uniform lateral load with aspect ratios ( / =1.0), slenderness ratio ( / =100), tapering ratio ( / =1.5) and varying values of tapering equation ( from 1-4). This figure shows the effect of tapering equation on the large displacement behavior of a rectangular thin plate. from this figure, can be noticed that the effect of tapering decrease with increasing the order of tapering equation. A simplified computational procedure is used to study the large deflection analysis of a rectangular thin tapered plate under lateral uniform load. Approximate values can be obtained with a good accuracy when compared with other works. The effects of boundary condition, tapering ratio and plate aspect ratio on the large deflection behavior of rectangular thin tapered plates are studied and presented in graphs. The tapering ratios are taken to be ( / =1.0,1.5,2.0). It is shown that the large deflection behavior is dependent on the tapering ratio. So the effect of order of tapering ratio shows that the deflection of steel plate increases with increase the order of tapering ratio. The effect of tapering ratio is dependent on plate aspect ratio. In the present study, the material nonlinearity, initial imperfection, free and rotationally boundary condition are not considered. Aalami,B., and Chapman, J.C., “Large Deflection Behavior of Orthotropic Plates under Transverse and In-plane Loads”, Proc. ICE, London, England, Vol.42, Mar., 1969, pp.347-382. Amash,h.k., “Post-buckling and Post-yielding analysis of Imperfect thin plate by finite difference method.” M.Sc. thesis, University of Babylon, February 2003. 162 Basu, A.C., and Chapman J.C., “Large Deflection Behavior of Transversely Loaded Rectangular Orthotropic Plates”, Proc. ICE, London, England, Vol.35, Sept., 1966,pp.79-110. Bjelajac, N., “Evaluation of Post-buckling Equilibrium Branches for perfect and Imperfect Plates by Finite Difference Method”, Numerical Methods in continuum Mechanics, Liptovskk َy Jaَn, Slovac Republic, 2000, pp.1-18. Chajes, A., “Principles of Structural Stability Theory”, Prentice- Hall, Inc., Englewood Cliffs, New Jersey, 1974. Chia, C.Y., “Nonlinear Analysis of Plates”, McGRAW- Hill International Book Company, 1980. Coan, J.M., “Large Deflection Theory for Plates with Small Initial Curvature Loaded in Edge Compression”, J. Appl. Mech., Vol.18, June, 1951, pp.143-151. Colville, J., Becker, E.B., and Furlong, R.W., “Large Deflection Analysis of Thin Plates”, ASCE, J. Struc. Div., Vol.99, No.ST3, Mar., 1973, pp.349-364. Husain, H.M., Al-Daami, H., and Amash, H.K., “Buckling Behavior of Rectangular Plate with Variable Thickness.” Engineering & Technology (Baghdad) Vol.21, No.10, 2002, pp.736-745. Klein, B. "Shear Buckling of Simply Supported Rectangular Plates Tapered in Thickness", J.F.I.,June 1957, pp537-541. Jayachardalan, S.A., Gopalakrishnan, S., and Narayana, R., “Explicit Incremental Matrices for the Post-buckling Analysis of Thin Plates with Small Initial Curvature”, An International Struct. Eng. And Mech., Vol.12, No.30, 2001, pp.283-295. Levy, S., “Bending of Rectangular Plates with Large Deflections” NACA Report 737, (U.S.A.) 1942. Murray, N.W., “Introduction to Theory of Thin Walled Structures”, Oxford Univ., Press, New York, 1986. Rushton, K.R., “Large Deflection of Plates with Initial Curvature”, Int. J. Mech. Sci., Vol.12, 1970, pp.1037-1051. Salvadori, M.G., “Numerical Computation of Buckling Loads by Finite Differences”, ASCE, Trans., Dec., 1949, pp.590-636. Saxena, H.C., “The Calculus of Finite Differences”, S. Chand & Company, LTD, 1980. Szilard, R., “Theory and Analysis of Plates: Classical and Numerical Methods”, Prentice- Hall, Inc., Englewood Cliffs, New York, 1974. 163 Timoshenko, S.P., and Woinowshy, S., “Theory of Plates and Shells”, 2nd Ed, McGraw-Hill Book Co., Inc., New York, 1959. Timoshenko, S.P., and Gere, J.M., “Theory of Elastic Stability”, 2nd Ed, Mc Graw-Hill Book Co., Inc., New York, 1961. Ueda, Y., Rashed, S.M.H., and Paik, J.K., “An Incremental Galerkin Method for Plates and Stiffened Plates”, Comp. & Struct., Vol.27, No.1, 1987, pp.147-156. Ugral, A.C., “Stresses in Plates and Shells”, McGraw-Hill, Inc, 1981. Wang, P.C., “Numerical Matrix Method in Structural Mechanics”, John Wiley & Sons, Inc., New York, London, Sydney, 1966. Williams, D.G., and Aalami, B., “Thin Plate Design for In-Plane Loading”, Halsted Press, A Division of John Wiley & Son, Inc., New York, 1979. 4 2423 23 22 222 42 234 )(2 42 234 )(2 44 2423 23 4 Plate equation in finite difference molecule form(2) 164 Load-deflection curve of a square simply supported thin plate under transverse uniform load s s s s Central deflection of simply supported square plate versus uniform lateral load 165 Load-deflection curve of a square thin plate with all edges clamped and under transverse uniform load Load-deflection curve of a square thin plate with two edges simply supported and other edges clamped under transverse uniform load 166 Load-deflection curve of a rectangular simply supported plate with aspect ratio ( / =2.0) and under transverse uniform load Load-deflection curve of a rectangular simply supported plate with aspect ratio ( / =0.5) and under transverse uniform load 167 Load-deflection curve of a rectangular simply supported plate and under transverse uniform load with different values of order tapering equation