Curved Panels Under Combined

Similar documents
MINISTRY OF SUPPLY AERONAUTICAL RESEARCH COUNCIL REPORTS AND MEMORANDA. from Lag Standpoint. Sheet. Mo FtNE~ BoAo.

Note on the Oharactensc~c Ourve for

Strut Intended for Compressive.eads

Note on the Convergence of Numerical Solutions of the Navier-Stokes Equations

On Certain Types of Boundary-layer Flow with Continuous Surface Suction

On the Design of a Row of Windows in a Pressurized Cylindrical Fuselage

8; Theo U of Sandwich

A with Leadmg-e_~, _e uct~on

P a g e 5 1 of R e p o r t P B 4 / 0 9

/ 4 -DE{ Stress Concentrations at a. Cut-out. R. & M. No (14,~9) (A.R.C. Technical Repot0 MINISTRY OF SUPPLY AERONAUTICAL RESEARCH COUNCIL

The Stress Distribution in Panels Bounded ~ by Constant-Stress

The Flow in an Axially-Symmetric Supersonic Jet from a NearlySonic Orifice into a Vacuum

The Shear Stiffness of a Corrugated Web

A L A BA M A L A W R E V IE W

176 5 t h Fl oo r. 337 P o ly me r Ma te ri al s

Theory of Parachutes with Cords over the Canopy

Flow in a P ipe o f Rectangular Cross-Section.

COMPOSITE PLATE THEORIES

Unit 15 Shearing and Torsion (and Bending) of Shell Beams

UNIT- I Thin plate theory, Structural Instability:

the Position df an -Equivalent

THE POISSON TRANSFORM^)

C~:~nditions for the Existence of an Inertial Subrange in Turbulent Flow

I9 Lateral deflections of circular plates

The Forces on a Cylinder in Shear Flow

ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD

T h e C S E T I P r o j e c t

APPENDIX 1 MODEL CALCULATION OF VARIOUS CODES

On the Flexure of Thin Built-Up Wings

A Simple Method of Computing CD from Wake Traverses at Highsubsonic

Limitations of Use of Busemann's Secondorder Supersonic Aerofoil Theory

Accordingly, the nominal section strength [resistance] for initiation of yielding is calculated by using Equation C-C3.1.

Upwash Interference onwings of Finite Span in a Rectangular Wind Tunnel with Closed Side Walls and Porous-Slotted Floor and Roof

Discontinuous Distributions in Mechanics of Materials

I M P O R T A N T S A F E T Y I N S T R U C T I O N S W h e n u s i n g t h i s e l e c t r o n i c d e v i c e, b a s i c p r e c a u t i o n s s h o

I N A C O M P L E X W O R L D

A *69>H>N6 #DJGC6A DG C<>C::G>C<,8>:C8:H /DA 'D 2:6G - ( - ) +"' ( + -"( (' (& -+" % '('%"' +"-2 ( -!"',- % )% -.C>K:GH>IN D; AF69>HH>6,-+

) sm wl t. _.!!... e -pt sinh y t. Vo + mx" + cx' + kx = 0 (26) has a unique tions x(o) solution for t ;?; 0 satisfying given initial condi

P a g e 3 6 of R e p o r t P B 4 / 0 9

2. Write your full name and section on the space provided at the top of each odd numbered page.

Analytical Strip Method for Thin Isotropic Cylindrical Shells

M5 Simple Beam Theory (continued)

OH BOY! Story. N a r r a t iv e a n d o bj e c t s th ea t e r Fo r a l l a g e s, fr o m th e a ge of 9

ii!i ,..~ ~.in-~shear~ Flow i."!i ' ' c,o,,,,; c,,migl, t ~S,,,,,a R, & M. No s,235) A.R.C. Tedmical Report :,. MINISTRY OF SUPPLY

FREE VIBRATION OF A RECTANGULAR PLATE WITH DAMPING CONSIDERED*

the use f~s Prattles! C~mb ~echmflue

Comb resonator design (2)

, N U No TORSIOI_AL S _TP_ENG_'w OF N!oKEL,4 STEEL AND DURALU}_iN TUBING. AS AFFECTED BY THE RATIO OF DIA}:_ETER TO CCaGE THICKNESS.

Some Additions~ Notes on the eriw.tion of Airworthiness Performance Climb Standards

Elastic buckling of web plates in I-girders under patch and wheel loading

Flexural-Torsional Buckling of General Cold-Formed Steel Columns with Unequal Unbraced Lengths

Chapter 3 Higher Order Linear ODEs

Optimum Height of Plate Stiffener under Pressure Effect

Longitudinal buckling of slender pressurised tubes

Partial Differential Equations (PDEs)

NARAYANA I I T / P M T A C A D E M Y. C o m m o n P r a c t i c e T e s t 05 XII-IC SPARK Date:

Stability of Simply Supported Square Plate with Concentric Cutout

Cyby TISIA E. Document No. - "/r ".. TECHNICAL LIBRARY /, Copy No. /- ' E. E. Secher. 2 uut15. No. EM JReport

ENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1

CONVECTION CURRENTS IN THE EARTH.

Structures. Shainal Sutaria

UNIVERSITY OF SOUTHAMPTON. A foreign language dictionary (paper version) is permitted provided it contains no notes, additions or annotations.

Chapter 3. Load and Stress Analysis. Lecture Slides

Comb Resonator Design (2)

Mechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection

Unit I Stress and Strain

Buckling Analysis of Axially Compressed Square Elastic Tubes with Weakly Supported Edges

2. (a) Explain different types of wing structures. (b) Explain the advantages and disadvantages of different materials used for aircraft

ON LAMINAR FLOW THROUGH A CHANNEL OR TUBE WITH INJECTION: APPLICATION OF METHOD OF AVERAGES*

1-1 Locate the centroid of the plane area shown. 1-2 Determine the location of centroid of the composite area shown.

LAMINAR FLOW BEHIND A TWO-DIMENSIONAL GRID

Torsional Buckling of Double-Walled Carbon Nanotubes

Shafts: Torsion of Circular Shafts Reading: Crandall, Dahl and Lardner 6.2, 6.3

Diffraction by a Half-Plane LL. G. CHAMBERS. {Received 29th March Read 5th May 1950.)

1 Review of di erential calculus

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

7. STRESS ANALYSIS AND STRESS PATHS

Unit 18 Other Issues In Buckling/Structural Instability

/ =0. (c) Section P.3 Functions and Their Graphs. (c) g(-2) = 5-(-2) 2 = 5-4 = 1. (e) g(x) = 0 for x = -I, 1 and 2. 2.

FIXED BEAMS IN BENDING

ri [11111 IlL DIRECTIONAL PROPERTIES Of GLASS-FABRIC-BASE PLASTIC LAMINATE PANELS Of SIZES THAT DO NOT IBUCICLE (P-Q1lAtVjr) No.

By C. W. Nelson. 1. Introduction. In an earlier paper by C. B. Ling and the present author [1], values of the four integrals, h I f _wk dw 2k Ç' xkdx

EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES

EJhsticity of a Sheet Reinforced by Stringers and Skew Ribs~ w~th Applications to Swept ~v'j ngs

APPENDIX A Thickness of Base Metal

Heat-Transfer Calculations for the Constant-Property Turbulent Boundary Layer and Comparisons with Experiment

Combined Flexure and Torsion of Class of Heated Thin Wings: A Large-Deflection Analysis

MATH 251 Final Examination August 14, 2015 FORM A. Name: Student Number: Section:

The Periodic Table. Periodic Properties. Can you explain this graph? Valence Electrons. Valence Electrons. Paramagnetism

Stress Concentration Factors for Reinforced Rounded-Square Holes in Sheets

Executive Committee and Officers ( )

Parts Manual. EPIC II Critical Care Bed REF 2031

COMPRESSIVE EUCICLING CURVES MR SANDWICH PANELS WITH ISOTROPIC FACINGS AND ISOTROPIC OR ORTI1OTROIPIC CORES. No Revised January 1958

CHAPTER 5. Beam Theory

THE AERODYNAMICS OF A RING AIRFOIL*

Basic Energy Principles in Stiffness Analysis

(308 ) EXAMPLES. 1. FIND the quotient and remainder when. II. 1. Find a root of the equation x* = +J Find a root of the equation x 6 = ^ - 1.

Structural Mechanics Column Behaviour

The problem of isotropic rectangular plate with four clamped edges

Sixth Term Examination Papers 9475 MATHEMATICS 3

Transcription:

]; "!" % ~,, R. & M. No. 197S (7514) A.R.~. Technical ~N MINISTRY OF AIRCRAFT PRODUCTION AERONAUTICAL RESEARCH COMMITTEE REPORTS AND MEMORANDA The Initial Buckling o Slightly Curved Panels Under Combined Shear and Compression By D. M. A. LsGG~:rT, Ph.D. o Cro~n C~pyrlgkt R,ser~ed LONDON: HIS MAJESTY'S STATIONERY OFFICE Price 3s. od. net f-

btational A~2}]C;... I~;:C" ;.. g... ~a.~'-a~v~, g,,~, The Initial Buckling o} t l l Curved Panels Under Combined Shear and Compression By D. M. A. L~.QC~TT, Ph.D. COMMUNICATED BY THE DIRECTOR OF SCIENTIFIC RESEARCH, MINISTRY OF x3klrcraft PRODUCTION Reports and Memoranda No. I9 7 2 Decem[ er, 194-3* Summary.--The initial buckling of flat rectangular panels under combined shear and compression has been investigated theoretically in R. & M. 1965. This report extends the results given there to panels which are long and slightly curved. On aircraft with laminar flow wing sections, it is desirable that the wing cover should remain smooth up to a factor of 141g, and to achieve this a possible type of construction is one in which stringers are dispensed with, and the cover is reinforced with closely spaced ribs and stiffeners. These divide the cover into a large number of long and slightly curved panels, and the results given in this report should be of value in estimating the combined shear and compression which such panels can carry without buckling and so developing waviness. LIST OF CONTENTS I. Introduction.................. 2. Statement of Problem and Method of Solution........ 3. Description of Results................ 4. Conclusions.................. List of References.................. Appendix I---The Inftial Buckling of a Long and Slightly Curved Panel in Shear Appendix II--The Initial Buckling of a Long and Slightly Curved Panel in Compression Page I 2 3 4 4 5 12 1. Introduction.--The initial buckling of flat rectangular panels under combined shear and compression has been investigated theoretically in R. & M. 1965, and the purpose of this report is to extend the results given there to panels which are long and slightly curved. For flat panels, it is shown in R. & M. 19651 that the combination of shear stress q, and compressive stress f, which will just cause buckling are connected by the relation-- and for the slightly curved panels considered in this report it is assumed that the same relation is approximately correct. 2 On this basis the problem amounts therefore to finding % and f, for various curvatures and dimensions of panel. (71550) * R.A.E. Report No. S.M.E. 3274, received 6th March, 1944.

2 The buckling of a short cylinder in torsion, of which a long and slightly curved panel in shear can be regarded as a special case (Fig. 1), has already been investigated theoretically by L. H. Donnell. a His treatment, however, involves a number of scarcely justifiable assumptions about boundary conditions, and the more accurate solution developed in this report shows that Fro. 1 Donnell's approximate values for the initial buckling stress are appreciably too high. A comparison of the theoretical results found here with the experimental values obtained by E. E. Lundquist ~ shows that for the small curvatures that exist over the major portion of a wing the agreement is good. For larger curvatures, however, owing to the increasing importance of initial irregularities, experimental values are less than the theoretical, and the latter are to be regarded as an ideal upper limit for panels which are perfectly formed and accurately loaded. The symmetrical buckling of a cylinder under end load, of which a long and slightly curved panel in compression can be regarded as a special case (Fig. 2), has been investigated theoretically by R. V. Southwell, S. Timoshenko ~ and others, and what is done here is to sol;ce the fundamental IITTTII t FIG. 2 equations derived by them for alternative boundary conditions. For slightly curved panels there are very few experimental results available. But such evidence as is supplied by tests on curved panels whose straight edges are much longer than their curved ones, suggests that even for the small curvatures considered in this report experimental values of the buckling stress are somewhat less than the theoretical. 2. Statement of Problem and Method of Solution.--The problem considered is the initial buckling of a long and slightly curved panel under combined shear and compression. The applied shear is constant round the panel, and the compression is uniformly distributed over the two curved edges (Fig. 3). Owing to the length of the panel the type of support for the two straight sides is unimportant, but for the two curved sides results are worked out on the assumption of clamped or simply supported edge conditions. The method of obtaining qc, and fc~ is explained in the appendices. /\ ///4 \\ // \ // \ \ \ 7 FIG. 3

3. Description E ----- V h ---- a b = r ~- K -~ fc, ~- fc,o = f kcrc kc qcr ~c~o q kcrs k, of Results.--The following notation is used throughout :-- Young's modulus Poisson's ratio (taken as 0.25) semi-thickness of panel width of panel, measured along short straight edge length of panel, measured along long curved edge radius of curvature of panel a non-dimensional constant ----[3 (1 -- v ~)/z~'~] 1 i~ a 2/rh stress at which panel buckles under compression alone value of f0r when panel is flat and the edges are simply supported 3 = compression stress at which panel buckles when the shear stress is q =fo,/fo,o ---= stress at which panel buckles under shear alone -~ value of q0, when panel is flat, and the edges are simply, supported ---- shear stress at which panel buckles when the compressmn stress is f = q /%o = q/qo,o The results are shown graphically in Figs. 4 to 7. Fig. 4 refers to the case of pure shear, and gives the variation of k~r, with K for simply supported or clamped edges. Donnell's results are shown by the broken curves and give values for the buckling stress which are between 10 per cent. and 20 per cent. too high. Fig. 5 also refers to the pure shear case and shows how the wave length of the buckles varies with K. Fig. 6 refers to the case of pure compression, and gives the variation of kcrc with K when the edges are simply supported or clamped. Fig. 7 includes the results given in Figs. 4 and 6, and shows, for given K, i.e. for given curvature and size of panel, the various combinations of shear and compression which will just cause the panel to buckle. In Fig. 7 the full and broken curves refer respectively to simply supported and clamped edge conditions. Table 1 shows the values of K which correspond to typical values of a, r and h. TABLE 1.Values ~K I a 2h I r K a 2h r K 6 0.036 100 3.4 9 0.036 200 1.7 300 1.1 100 200 300 7"6 3"8 2'5 0.048 100 200 300 2.5 1.3 0.8 0.048 100 200 300 5.7 2.9 1.9 0.064 100 200 300 1.9 1"0 0.6 0.064 100 200 300 4.3 2-2 1.4 0.080 100 200 300 1"5 0'8 0"5 0.080 100 200 300 3"4 1"7 1"1 All dimensions are in inches. (71559).A 2

4 4. Conclusio~s.--This report shows under what combinations of shear and compression a long and slightly curved panel may first be expected to buckle. In aircraft with laminar flow wing sections, it is desirable that the wing cover should remain smooth up to a factor of lig, and to achieve this a possible type of construction is one in which stringers are dispensed with, and the cover is reinforced with closely spaced ribs and stiffeners. These divide the cover into a large number of long and slightly curved panels, and the results in this report should be of value in estimating the combined shear and compression which such panels can carry without buckling and so developing waviness. No. l Author H.G. Hopkins and B. V. S. C. Rao L~[ST OF REFERENCES Title.. " The Initial Buckling of Flat Rectangular Panels under Combined Shear and Compression ". R. & M. 1965. 1945. F. J. Bridget, C. C. Jerome, A. B. Vosseller " Some New Experiments on Buckling of Thin Walled Construction ". Trans. Am. Soc. Mech. Eng., Vol. 56, p. 569. 1934. 3 L.H. Donnell........ " The Stability of Thin Walled Tubes Under Torsion ". N.A.C.A. Report No. 479. 1933. 4 E.E. Imndquist.... " Strength Tests on Thin WaUed Duralumin Cylinders in Torsion ". N.A.C.A. Technical Note No. 427. 1932. 5 S. Timoshenko........ " Theory of Elastic Stability ". McGraw-Hill, 1936. Chap. 9, p. 81.

5 APPENDIX I The Initial Buckling of a Long and Slightly Curved Panel in Shear FIG. 8 Additional notation used in Appendix I :-- G = shear modulus = E/2(1 + v). The co-ordinate axes OX, 0 Y, OZ are shown in Fig. 8, and are such that OX and 0Y are the generator and line of curvature through the mid point of one of the curved edges, and 07. is normal to the middle surface. Referred to these axes the edges of the panel are X = 0, a ; Y = ± b/2 ; and the equilibrium displacements %, v 0, w 0, are such that u o:0, Vo-- qo, X G ' Wo=0. If this configuration is one of neutral equilibrium Uo+U, Vo+V, Wo+W, are also possible displacements, where u, v, w are indefinitely small but not all zero. Substituting each of these sets of displacements in the shell equations obtained by W. R. Dean* for problems of this kind and assuming that a/r and h/a are small and of the same order, the three fundamental stability equations reduce to 2 ~-X\~x+~-Y +~x ~ + ~-X h2v~ ~, ~w 1 ~V w + v -~ w + 2 (1 -- v~) ~XOY r O-Y r --0, --0. * Proc. Roy. Soc. A, Vol. 107, 1925, p. 734.

Putting these into non-dimensional form by means of the substitutions ax ay X=--,Y-- ~ 7g we have aw ( -- O av Ou... ~x ~x + ~ ~y ~-r + 2 --- ~y -Fx-+ = O, 6 2 ~x ~-x + +b~ Oy ~r + ~-~ =0, gv w+2(1-- v ~) a 2q~, 02 w av aw + v =0. \ ~ E ~xoy ~ar Oy ~r ~-X Introducing a stress function f for the purpose of simplifying the analysis, there result the following five equations (only four of which, however, are independent, since equation (4) is obtained by eliminating u and v from (1)-(3)), (1- v) @ ~v) ~2f... 2-- + ~x = -- ~ gxa-~........... ~u ~ Ov aw ) ~2f ~--x + ~ Oy 7; = ~... Oy2 '........... ( ~v aw ) ~u ~o.f ~y ~- + ~ ~x = ~ --2 Ox............. (1) (2) (3) V4f+ R ~2w ~x 2 -- O, V a P ~2f ~ 2 w w-- +Q - 0, ~x 2 ~x~y (4) (s) wbere Case I. Edges Simply Supported.*--It is now required to find a solution of the equations (1) to (5) which satisfies the boundary conditions U--~ V =0, for x = O, ~ ; and a solution of w andfis 32w.... w -- ~x~. = O, sought for in the form...... (6) w = w 1 cos my + w 2 sin my,)._ f --fl cos'my + f2 sin my. f....... (7) Here Wl, w2, fl, f2, are functions of x only, and m is real but otherwise unspecified (except in the case of a complete cylinder when zrm/a must be an integer). We express w, (s = 1, 2) in the form oo w, -- ~ A~t sin tx~, 0 <~ x <~ t=l * For a more detailed description of the method used, see writer's paper on a similar problem in Proc. Roy. Soc. A, Vol. 162, 1937, p. 62. t By considering the boundary conditions which w~ must satisfy, it can be shown that it is legitimate to differentiate this serie~ term by term four times.

- - RYm and introduce the complex quantities W, F and A,, defined by the relations W ---- w l + iw2, F = fl + if2, At = AI, + ia2,. By this means w and f in (7) can be written in the form where (R) denotes real part. w = (R) We -''~v, f= (R)Fe-'"%...... Substituting for w and f in (4), F is then given by the differential equation oo d~f,~ 2d2F dx* ~m -d-if2 + m4f = R ~ A, t 2 sin tx, t=l 7 (s) the general solution of which is F = {A cosh mx + B sinh rex} + x {C cosh mx + D sinh rex} O0 + R ~ A, t 2 sin tx.. (9),=1. ( + 2........... The next step is to determine the arbitrary constants A, B, C, D in terms of At from the first two boundary conditions in (6). Since, however, these involve u and v, it is first necessary to express u and v in terms of F and W. This is done by solving the equations (1), (2) and (3), and the boundary conditions can then be expressed in the form d2f dx 2 + ~ m ~.F = O, dafdx 3 (2 + v) m 2dF~ + (1 --~2r,,2) a dwdx -- O,...... (10) for x = 0, ~. Now substituting for F and W in (10), the arbitrary constants A, B, C, D are given by AK = 2RXm (1 + v) {(3 --,,) sinh m~ cosh m~r + m~ (1 + ~,)} where -- 2RYm (1 + v) {(3 -- ~) sinh m.~ + rnx (1 + v) cosh m~}, BK -- RXm (1 q- v) {ml~ 2 (1 + v) 2 -- 2 (3 -- v) sinh 2m~} 2~ (1 + v) 2 (1 --,,) sinh m~, CK-----RXm 2 (1 + ~)~ (3- v)sinh ~m~r --RYm3~ (1 + v)3 sinh m~, DK-~ -- RXm ~ (1 + v) 2 {(3 -- v) sinh m:z cosh m~ +mx (1 + v)} + RYm ~ (1 + v) ~ {(3 -- v) sinh m~ + mr: (1 + v) cosh m~}, K = m 2 (1 + v) 2 {(3 -- v) 2 sinh 2mu -- mx ~ (1 + v)2}, X=~=I A, (t~ + m2)2 =-~~=1 AtKt' cx3 t=j. These expressions for and f given by (8) will (-), A, (t2 + m )2 =- (--)' A, K,. t=l A, B, C, D are now substituted in equation (9), with the result that w satisfy all the fundamental equations and boundary conditions except (5).

8 It accordingly remains to substitute for w andf in (5) and equate to zero the coefficient of e -~'~. Doing this we have oo o0 oe l 4 sin tx A t(t 2+m2) 2sintx-imQ~. ta, costx+pr~ A, t=l t=l t=l ( t2 @ 7~2) ) -- Pm 2 {A cosh mx+ B sinh mx+ x (C cosh mx + D sinh rex)} --2Pro (D coshmx - - C sinh rex) = 0,......... (1 1) which must be satisfied for all x in (0, ~). The next step is to express in a Fourier series of sines all those terms in equation (11) which are functions of x, and then to equate to zero the coefficients of the sines. In this way we deduce, after considerable algebraic reduction, the following system of equations Go P Rt4 (l 2 + me) z A t -- imq ~=l c,,, A,~ + (t 2 _}_ mz)zat Here t =- 1,..,o... where and + (PRL t K,,A,, + (_)t ~ (--)"K,,A. = 0..... (12) oo, and the c's and L's are given by ~ [1 - cos (t + n) L 4mKt T L, = (1 + v)' + 1 - cos (t - n) ~] -- 1~ J TN -- {(3 -- v) cosh r~ sinh m~ + m~ (1 +,.)} -- cos t~ {(3 -- v) sinh m~ + m~ (1 + v) cosh m,~}, N = (3 -- v)~sinh 2m~ -- m2~ ~ (1 + v) 2. Since the equations (12) are linear in the A's, their only solution is in general that in which all the A's vanish. If however the determinant formed by eliminating the A's itself vanishes this is no longer true, so that the vanishing of this determinant provides the required equation for qcr. The order of this equation being infinite, qc, has an infinite number of roots, one corresponding to each possible form of instability. But since we are only concerned with that form of instability which is most likely to occur in practice, it is only the smallest value of qc, which is of interest. Introducing Z, defined by the relation PRt 4 Zt _~ (t~. _}_ m 2) 2 _}_ ( t2 + m2)., + 2 PRK, L,, the determinantal equation for the critical shear stress is Zl -- imqc12 2RPKzL 1 -- imqcl4 - imqc21 z 2 - imqc23 2 RPKaL 2 2RP KaLa - irnqc~2 Za -- imqc34 -- imq Q a 2RP K~L 4 -- imqc 4. ~ Z4 =0. (13) To show that the solution given by equation (13) is not merely a formal one, it is necessary to show that the determinant is convergent. The proof however is on exactly the same lines as that given in the paper already referred to* and so is omitted. * Loc cit. p. 6.

9 As it is not possible to solve equation (13) directly, it is necessary to approximate to the solution for Q, i.e. for qo,, by replacing the infinite determinant in that equation by the finite determinant which contains its first n 2 (n = 2, 3,...) elements. The results of the successive approximations, and those obtained by Donnell, are given in Table 2. TABLE 2 K First Approximation 7/q, 2 kc~s Second Approximation m 2 kcr s Third Approximation m 2 kcrs Donnell's Results kc,, Percentage Increase 0 1 2 3 4 5"5 7 10 15 0"59 O" 75 1 "05 1"4 1"7 2"1 2"6 1.046 1.134 1.307 1.507 1.712 2.020 2.325 0"63 0"78 1"1 1 "45 1"8 2"3 2"7 3"6 4"9 1.000 1.071 1.211 1.367 1.524 1.754 1.981 2.420 3.114 m 0.78 1.8 2.7 3-7 5.05 1.070 1. 521 1.976 2.400 3.069 1.003 1.144 1.331 1.515 1.693 1.947 2.189 2.643 3.331 0.3 6.9 10.0 11.0 11.3 11.3 10.8 10.1 8.6 Case II. Edges Clamped.--In this case we again try for a solution of w and f in the form (7), but the boundary conditions to be satisfied are now U -~- V ------ 0, bw w = ~--~= 0, for x = 0, =, and it is no longer possible to express ws as a sine series which can be differentiated term by term four times. Instead, it can be shown by expanding d~ws/dx ~ as a sine series and integrating four times, that the corresponding expression for w, (s = 1, 2) is oo ws = ~, As, sin tx + H, + E,x + J,x 2 + Gsx a....... (14) t=l In (14) the sine series is such that it can be differentiated term by term four times, and Hs, E. J,, Gs, are given by H,=0, z=l {2+(--)'}tA,,, t=l OO E,=--Z,A,,, G, -- ~=1 1 oo ~ ~. {1 + (--)'} IA,,. t=l Introducing E, J and G, defined by the relations E = E 1 @ ie 2, J = J1 + i J2, G = G1 + ig~, we have, after substituting for w and f in (4), the following differential equation for F d4 F,, ~ d2 F { oo ] dx 4 xm -d~ ~ + m4f = R ~. A~t 2 sin tx -- 2.[ -- 6 Gx, t=l

- - imq which has for its general solution, 10 F -- {A cosh mx + B sinh rex} + x {cosh mx + D sinh rex},=1 (t~+m~) 2 m 4 " The general procedure is now very similar to that in Case I. The arbitrary constants A, B, C, D are determined from the boundary conditions, and in place of (11) there ultimately results the equation OO t=l A, (t 2 + m2)2 sin tx + {m ~ (Ex +Jx 2 + Gx a) -- 2m ~ (2J + 6Gx)} oo o~ sin tx -- imq ~. ta, cos tx -- imq (E + 2Ix + 3Gx ~) + PR ~. A,t 4,=i,=i (t ~ + m 2) ~ -- Pm 2 {A cosh mx + B sinh mx+ x (C cosh mx + D sinh rex)} -- 2Pro {C sinh mx + D cosh rex} = O. Expressing the left-hand side of this equation as a sine series, and equating to zero the coefficient of sin tx, it follows that oo 00 O0 (t2 +m2)~a,+ m47~" " e,,a,--4m2~, j,,a,, -- imq.~ c,,~a,, n =1,~ = i n =i ~,,=1 gt,, A, + (t 2 + m2) 2 + PR L, t,=1 A,, (K, + q,,) + (--)~ ~. (--)" A,(K,, + qj -- (U, + S,) ~ A,,p~, + (--)' Z (--)" A,,p,~ -- O,.. (15) where the c's, K's and L's have the same meaning as in Case I, and the e's, ]"s, g's, p's, q's, U's and S's, are given by 4n (_),~, et,, = - ~t312 + (_)t + + 2 ( ),+,~] Jr,," = ~2~ [2 + (-)'+ (-),* + 2 (-)"-J, g,o _ 2~,.~ {1 - (-),+,,} + ~12~ {1 - (-),} {1 + (-)j, 2vn {2 + (--)'} n 6(2+,,) 2,~t [2~y (1 + ~) 2m~ (1 + ~)~ U~:N~(I+~,)L (t 2+m2~ c st~sinhcn~ + (t z+m~)z costzsinhm~ 1 + (t ~ + m2 ) {(3 -- v) (1 --,,) sinh 2mz~ --,n2:~ ~ (1 + v) ~} --2m ~ (1 + v)(3-,,) ] -- (t 2 + m~.)~ sinh 2m~], 4 tm St =xn (t 2 + m ~) [(3 -- ~) sinh 2rex -- m~ (1 + v) cost~ sinh m~].

11 The equation for q,, is now found from the condition that there exists a non-zero solution of the system of equations (15). Introducing N,, defined by the relations PRi 4 N. = (i 2 + mz) 2 + m~tei~--4m~j. + (i 2 + m 2) Ns~ : (s t) + 2PREL,(K,+ q,)--p,(u, + $3" ], m a e,~- 4 m2 ~ this gives rise to the determinantal equation +PR {1 -I-(-)'+')EL, (K, --F q, -- p, ( u, -4- S,)'] Nix -imq (c12 +g1~) NI3 -imq (q4 +g~4) -- imq (c21-j-g21) N22 -- ilgq (c23 +g23) N24 NsI -iraq (ca~+gs~) Na3 -iraq (c34+gaa) -imq (cal+g~a) N42 -imq (ca3+ga8) N~4 ----0. (16) As in Case I, it is now necessary to consider the convergence of the determi~aant in (16), but the proof is omitted as it is very similar to that given in the paper referred to.* The procedure is now the same as in Case I, except that attention is confined to the second approximation, found by taking the first three rows and columns of the determinant. The results, together with those derived by Donnell, are given in Table 3. TABLE 3 Second Approximation D onnell s Results K m ~ kc,. s kcrs Percentage Increase 0 1 2 3 4 5.5 1.5 1.6 1.7 1.9 2.2 2.65 1.664 1.689 1.753 1.846 1.952 2.130 1.680 1.787 1.949 2. 124 2.302 2.560 0.9 5.8 11.2 15.1 18.0 20.2 *Loc. cit. p. 6.

12 APPENDIX II The Initial Buckling of a Long and Slightly Curved Panel in Compression Additional notation used in Appendix II :-- FIG 9. D = flexural rigidity = 2Eh3/3 (1 -- v2). The co-ordinate axes OX, OY, OZ are shown in Fig. 9, and the equation of neutral equilibrium, applicable to types of distortion in which there are no displacements parallel to the Y axis, is-- dhe~ d2w w O -d-fi + 2hfc" dx 2 + 2Eh ~ = 0......... (17) It remains to find the smallest values of fc~ for which there exists a non-zero solution for w satisfying the required boundary conditions.* Case I. Edges Simply Supported.--For this case the boundary conditions are W -- d2w dx ~. --0, forx= ± ~; a mt~x w is A sin -- mt~x (m even) or A cos -- (m odd); and the corresponding values of f,, are given by Dm2 ~ 2 E a 2 fo-- 2ha~ + m2~r2............. (18) Case II. Edges Clamped.--The boundary conditions are now ZO -- dw dx = O' f r x = ± ~-; and the solution is in this case more complicated. * See S. Timoshenko, " Theory of Elastic Stability ", First edition, p. 81.

13 If w is an odd function of x, i.e. if the form of distortion is anti-symmetrical with respect to 0, w is A, sin ~x + B, sin/~x, where A,, B, are constants whose ratio is determined by the boundary conditions, and the equation for fc, is ~2= If..-F {f J-- 2ED h/r 2 }I/,']/D,,': [s..--{s:- h/v, c~a ~a /~a ~a (19) 2 tan 2 -- 2 tan 2---............. If w is an even function of x, i.e. if the form of distortion is symmetrical with respect to 0, w is A c cos c~x + B c cos ~x, where A c, B~ are constants whose ratio is determined by the equation for fc, is boundary conditions, and the -~ ~a tan -~ ~a = --~--tan Sa ~a 2 "" "......... (20) Owing to the different types of distortion which are theoretically possible, the equations for f,,, i.e. (18), (19) and (20), are multi-valued, and the two curves given in Fig. 6 are the envelope of all possible solutions when the edges are either clamped or simply supported.

3.0 14 -I~ EDGES CLAHPED ~ ~ DONNEL,t'5 5OLUTION 8"0 ACCURATE.,,SOLUTION~ ~ ~ ~ ~ ff-~ 1-- ~ i ~DGF..5 51HPL,.Y,~LJP~)RTED, ~.~ [~ONI~E LLtS 5QLU'[ION..~j. ~ ~ ~ ACCURATE ~ - Ze, t ~ J VARIATION OF 6UCKLIM G 5TR~.SS WITH CUR.VATURF" ~v'he~ I..,OADtNC~ IS PURE. SHEAR,. " 'o~ = ~c-,- / t,=.,~ K =,,o,~/.,p ~,, 0 I'0 8-0 3'0 4-'0 5'0 PRACTICAL RJ~t.,JG E 6"0 7.0 8.0 9.0 ~0.0' K FIG. 4.3,0 Z'O RATIO OF WAVE t,. EJqGTH TO WIDTH OF' PANEL 1'5 ED~,E5 SlMPL.V ~JPI:~OR.TED VARIATION OF" WA_Vr_. LENGTH OF" BUCKLE- WITH CURVAYUR.~" WHEN LOADING is PURE. 5HEAR K ={ 3(~ -VO/'~" ~ ~ " 1.0 0"5 WAVE LF-N~TH I'0 ~.0 3"0,, PRACTIGAL RAI~GE 4.0 5.o K 6.0 7.0 8.o 9.0 FIG. 5 I0"0

15 ]~:,O i I0"0 8'0 //\ ' ktaj_u~ Tmxn-~,fc~ EE)GE8 CLANPED~, / / ftl c*c 4.c J / / // // 4. o VARIATION OF BUCKLING.STI~ES,S WITH CUIKVATUR.E WNEN LOADIN G 18 PURE C.OMPRE8810N 2-o J / 0 I'0 "2.O 3.0 <-'O ~.O... PRACTICAl. R,~ kjgf,..... K FIG. 6 14..0 13'0 I '.0 ~ VARIATION OF" IBUCKLING,STRESS WITH ~ CURVATURE %B/HEN LOADING 15 GOmfBIMED " " ~ SHEAR. AND DONPRE.SSIOM.... "'" ~o q4.3.0 2.0 I'0 ~ I~ ~L ", "- FULk CURVES REFERS TO PANF--L i " Q Q "-,% \~. WHOSE. EDGE~S ARE &imply,, x x SUPPORTED.... I... ~ --~~ I ~ ". x-~-... 6I~OKEM EUR"VE.8 REFER 'TO K = s --" ~. ", ", CLAMP~O I --- "<'.. "-.K=4",,K:s I "--'~~--~...,--- -i~--. ~'K:~X" " " " "" "" " ~ -- ~\"'- "- \x x ", "'~'~--{'-\ -X ", ", ", -..q -_~.\. -\ -..., ",. x x~ x N N x \ 0 0'@ 0.4-0"6 O-B (.o "~s I.~ 1.4 b6 1.8 2.0 FIG. 7 (71559) Wt. 9/7116 3/46 Hw. G377/1

Pubhcatlons of t~e J cl Aeronautical Researclh Commktee TECHNICAL REPORTS OF THE AERONAUTICAL RESEARCH COMMITTEE-- x934-35 Vol. I. Aerodynamics. 4os. (4os. 8d.) Vol. II. Seaplanes, Structures, Engines, Materia[s, etc. 4os. (4os. 8d.) I935-36 Vol. I. Aerodynamics. 3os. (3os. 7d.) Vol. II. Structures, Flutter, Engines, Seaplanes, etc. 3os. (3os. 7d.) x936 Vol. I. Aerodynamics General, Performance, Airscrews, Flutter and Spinning. 4os. (4os. 9d.) Vol. II. Stability and Control, Structures, Seaplanes, I937 Vol. Engines, etc. 5os. (5os. Iod.) Io Aerodynamics General, Performance, Airscrews, Flutter and Spinning. 4os. (4os. 9d.) Vol. II. Stability and Control, Structures, Seaplanes, Engines, etc. 6os. (6Is.) ANNUAL REPORTS OF THE AERONAUTICAL RESEA RCH COMMITTEE-- I933-34 Is. 6d. (xs. 8d.) ~934-35 Is. 6d. (xs. 8d.) April r, I935 to December 3 I, I936. 193 7 2s. (2s. 2d.) I938 Is, 6d. (xs. 8d.) 4 s. (4 s. 4d.) INDEX TO THE TE(~HNICAL REPORTS OF THE ADVISORY COMMITTEE ON AERONAUTICSw I9o9-I9Id ) Reports and Memoranda No. r6oo. 8s. (Ss. 5d.) Prices in brackets include postage. Obtainable from H" ~s Maj es ty's St" atloner y Office London W.C.2 : York ttouse, Kingsway Edinburgh 2: I3a Castle Street Manchester 2: 39-4I King Street Cardiff: I St. Andrew's Crescent Belfast: 80 Chichester Street or through any bookseller., S.O. C~de No. 23-i97z