Face Sheet Wrinkling NASA
Face Sheet Wrinkling NASA
Face Sheet Wrinkling NASA
Facesheet Structures
Wrinkling
in Sandwich
Robert Northrop
P.
Ley,
Weichuan
Lin,
and
Uy El
Grumman
Corporation,
January
1999
The
NASA
STI Program
Office
... in Profile
Since its founding, NASA has been dedicated to the advancement of aeronautics and space science. The NASA Scientific and Technical Information (STI) Program Office plays part in helping NASA maintain this important role. a key
CONFERENCE
PUBLICATION.
Collected papers from scientific and technical conferences, symposia, seminars, or other meetings sponsored co-sponsored by NASA. SPECIAL technical,
or
The NASA STI Program Office is operated by Langley Research Center, the lead center for NASA's scientific and technical information. The NASA STI Program Office provides access to the NASA STI Database, the largest collection of aeronautical and space science STI in the world. The Program Office is also NASA's institutional mechanism for disseminating the results of its research and development activities. These results are published by NASA in the NASA STI Report Series, which includes the following report types: TECHNICAL PUBLICATION. Reports of completed research or a major significant phase of research that present the results of NASA programs and include extensive data or theoretical analysis. Includes compilations of significant scientific and technical data and information deemed to be of continuing reference value. NASA counterpart of peer-reviewed formal professional papers, but having less stringent limitations on manuscript length and extent of graphic presentations. TECHNICAL MEMORANDUM.
NASA programs, projects, and missions, often concerned with subjects having substantial public interest. TECHNICAL TRANSLATION. language translations of foreign and technical material pertinent NASA's mission. Englishscientific to
Specialized services that complement the STI Program Office's diverse offerings include creating custom thesauri, building customized databases, organizing and publishing research results ... even providing videos. For more Program information Office, about the NASA STI
Access the NASA STI Program Page at http'//www.sti.nasa.gov E-mail your question [email protected]
to
Fax your question to the NASA Desk at (301) 621-0134 Phone the 621-0390 Write to: NASA STI Help Desk
STI Help
at (301)
Scientific and technical findings that are preliminary or of specialized interest, e.g., quick release reports, working papers, and bibliographies that contain minimal annotation. Does not contain extensive analysis. REPORT. Scientific and
NASA STI Help Desk NASA Center for AeroSpace 7121 Standard Drive Hanover, MD 21076-1320
Information
CONTRACTOR
NASA/CR-1999-208994
Facesheet Structures
Wrinkling
in Sandwich
Robert Northrop
P.
Ley,
Weichuan
Lin,
and
Uy El
Grumman
Corporation,
National
Aeronautics
and
Space Administration Langley Research Center Hampton, Virginia 23681-2199 Prepared for Langley Research under Contract NAS1-19347 Center
January
1999
Available
NASA Center for AeroSpace 7121 Standard Drive Hanover, MD 21076-1320 (301) 621-0390
ABSTRACT
The analysis
purpose
is to provide mode
a concise
summary
of failure
in sandwich
an exhaustive number
related to provide
of key papers
analysts
to guide
of a wide problem.
of available
wrinkling
is applicable
to a specific
design
CONTENTS
Section 1 2 INTRODUCTION ASSESSMENT FACESHEET 2.1 FACESHEET 2.1.1 ........................................ OF THE STATE-OF-THE-ART WRINKLING WRINKLING Structures FAILURE FOR PREDICTING ......................... .................... Facesheets and 7 Facesheets and 10 Composite Facesheets... 12 13 15 Correlation with 15 with 19 28 33 34 36 3 3 Page 1
Sandwich Solid
Cores
2.1.2
Sandwich Cellular
Sandwich Summary
............................ ...................
2.3 2.4 3 4
EFFECTS EFFECTS
IMPERFECTIONS LOADS
........................
CONCLUDING REFERENCES
................................
..........................................
ii
ILLUSTRATIONS
Figure Page
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
Global
Modes
in Sandwich
Structures
................
4 4 5 7 9 11 14
Uniaxial Strut
Load
........................................
of Gough,
of Hoff and Mautner 7 ..................................... Mode in Sandwich Mathematical (1) Failure (1) Failure (2) Failure (2) Failure Panels Models ...................... .................... Core Model Core Model Core Model Core Model ........ ...... ........ ......
Non-Harmonic Summary Reference Reference Reference Reference Reference Reference Reference Reference Reference Reference
Wrinkling
of Facsheet
Wrinkling
11 Test Data Mode 11 Test Data Mode 11 Test Data Mode 11 Test Data Mode
18 18 20 20 21 21 23 23
12 Test Data - Isotropic 12 Test Data - Anti-Plane 20 Test Data - Isotropic 20 Test Data - Anti-Plane 21 Test Data - Loading 21 Test Data - Loading
Core Model Core Model Core Model Core Model Parallel Normal
25 25
Reference 12 Test Data Plot of k z Versus Ratio of Wrinkling Wavelength to Core Cell Size ..................................................
26 in 31
19
iii
TABLES
Table Page
Summary of Theoretical Wrinkling Stresses for Sandwich Cores ............................................................................. 2 Important Assumptions Underlying Facesheet Wrinkling Predictions ....................................................................... Published Wrinkling Test Data Correlation to Theoretical Expressions Stress .................................................................
Struts With
Thick 14
Theoretical 15 for 29
iv
NOMENCLATURE
2)
Ao b
(see Figure
strut or panel width facesheet D12, facesheet laminate laminate bending bending stiffness stiffnesses in the direction of the applied load
DI
Dll,
EI
Young's flatwise
(isotropic)
G_
kl, k2
K0 L m P
Pcr
of wrinkling
half-waves
RE
Ps
S
2)
tc tl
W o
(see Figure
2) 2) 21)
thickness
(see Figure
initial imperfection
NOMENCLATURE
(CONT'D)
x,y,z 60
sandwich
strut or panel
coordinate
system
2)
amplitude
(see Equation
half-wavelength half-wavelength major Poisson's ratio in the direction of the applied load
_'cr
critical facesheet
'V
Vc
vl
(_Tdimp
compressive stress
stresses allowables
respectively respectively
compressive smallest
of Crwr _ and Crwry stress compressive compressive compressive compressive shear stress stress allowable stress stress stress allowable
"_core
through-the-thickness
vi
SECTION
INTRODUCTION
The the cost development of aerospace for primary of fiber structures structures. facesheets from composite has materials resulted and the drive in renewed in this paper, to reduce both the weight use and
interest
in the
construction
a sandwich
structure thick,
facesheets mainly
buckling
individually.
structure made
or cellular
may
be used. to is
exhibit
high
be manufactured In addition,
more
than equivalent
metallic
laminates
can be cured
in a single
for complex
Cocuring
for complex
furthermore, to laminated
these
stiffeners
composite
structures. and lightweight wrinkling. Since failure cores are prone to a type of local may exhibit by facesheet wrinkling Throughout is this
structures asfacesheet
sandwich of these
structures structures
carrying
Hence, and
prediction sandwich
of facesheet structures.
of reliable to modes
having sense,
wavelengths referring
in a more
general
wavelength. purpose of this paper wrinkling summary the purpose is to provide mode a concise summary of the state-of-the-art The objective wrinkling a smaller of the here is
of failure
in sandwich related
published number of
information understanding
to provide should
designers
a working to guide
of the when
their judgment
deciding
of a wide variety of available facesheetwrinkling designformulas is applicableto a specific designproblem. This work was performedunder Task 13of NASA ContractNAS1-19347,in supportof NASA's Environmental ResearchAircraft and SensorTechnology(ERAST) Program. The technicalmonitorwasMr. JuanR. Cruz.
SECTION ASSESSMENT
This assessment of treating describes to calculate wrinkling facesheet the results necessary as a strength
as a short
performed Section
"knockdown problem,
by considering a brief
manufacturing effects
irregularities loads.
(or imperfections).
2.4 contains
discussion
of the
of combined
2.1
FACESHEET Facesheet
WRINKLING has
ANALYSES been treated as a local, modes short wavelength by buckling sandwich mode in of
traditionally
phenomenon structures results all but wrinkling Results procedure loads stress starts is then (see
of possible small
buckling
exhibited in the
Figure
buckling
wrinkling
insensitive useful
conditions used
Many
to predict flat
model
loaded complex
strut. usual
extrapolated
to more
used with
to predict
the onset
of wrinkling
structures facesheet
the calculation
compared structures
derived
facesheets,
evidence
indicates wrinkling
stress
model
in reasonably
wrinkling
load estimates. a sandwich a length strut loaded a, and in uniaxial compression as shown in Figure 2. The strut The
to have
supported
at boundary
conditions
at its ends.
infinite
varies
half-wavelength,
A.
GLOBAL
BUCKLING
MODE
B.
SHEAR
CRIMPING
MODE
C.
FACESHEET
DIMPLING
SYMMETRIC D.
B.Ley-97.07/BTI -03
Figure
1. Global
and Local
Buckling
Modes
in Sandwich
Structures
tF
3_
z Y
.2_P
/4
B.Ley-97-07/BTI -04
_/
Figure
2. Sandwich
Strut Under
Uniaxial
Load
PCR
Curve 1 Curve 2 _
Curve 3 --.-------------------Sandwich Model, Symmetric Wrinkling Curve 4 Curve 5 PS _ _ Sandwich Model, GC, EC Finite, Facesheets Extensible Sandwich Model, GC, EC Finite, Facesheets Inextensible
B. Ley-97-07/BTI -05
_,
Strut Euler relative strut Pcr, can column. In most sandwich and
in the (Ref.
sandwich load,
2 as a from the
Timoshenko following
2, pp.
the buckling
be determined
1
Per
1
PE
+
tbGc
(1)
where
PE is the Euler
load, modulus
b is the width
wavelength
2 in Figure 1 vanishes
zero,
This shear
load
is usually mode
referred
to as the lb.
shear Since
crimping shear
load
in sandwich is actually
structures. a short
The
crimping
is shown
in Figure
crimping
wavelength
form of antisymmetric wrinkling, it can be calculated either using Equation 1 with short wavelengthsor Equation 2. In general,Equation 1 is generally usedfor the calculation of buckling loadsassociated with longerwavelengthmodesonly, while Equation2 is usedfor the calculation of the crimping load. Up to this point, somewhatsimple expressionshave been derived for buckling of the sandwich strut owing to the fact that the through-the-thickness Young's modulus of the core has been assumedto be infinite. The problem becomes substantiallymorecomplicatedwhenthis assumption is eliminated. If the through-the-thickness stiffnessof the core is consideredin the buckling analysis, short wavelengthbuckling modes,usually known asfacesheet wrinkling modes, arise. Two
facesheet symmetric simple wrinkling case, surface wrinkling with cases must be considered. to the resting varies middle surface In the first case, of the sandwich foundation. appears that the wrinkling and of the facesheets using is the respect can be predicted how 3.
model load
wrinkle If the
is antisymmetric
stiffness
antisymmetric antisymmetric
not occur
wavelength ld.
wrinkling
modes
in Figure wavelength
Antisymmetric accounting
wrinkling
buckling
mode shear
while
for the core through-the-thickness wrinkling wavelength load appears appears buckling
and transverse
how the antisymmetric of buckling facesheets inextensible), for long respect load during versus
wavelengths
antisymmetric
is neglected cannot
facesheets overall
the resulting
approximation
be used to predict
buckling wrinkling
wavelengths. to wavelength
to the antisymmetric
described
upon
the assumption
honeycomb are
or other sufficiently
cellular buckle
sandwich
structures,
facesheets
facesheets
type of local
instability
is known
asfacesheet
dimpling
(see Figure
2.1.1
With studies
and Solid
Cores were such performed study was were (3) the of the
struts first
facesheets and
solid, 4
Gough,
de Bruyne. directly
assumed surface
(1) the
facesheets and
attached
to the middle
of the core compressive could the be neglected. core having solving was the the
stresses Gough,
of the applied
4.
The core
in a state Elam,
boundary
of Gough,
and de Bruyne 3.
4 as well
results
are summarized
164 of Reference
(1)
B .Le y.-97.0 7/BTI
(2)
ANTIS YMMEYRIC WRINKLING MODE
-06
Figure
4. Wrinkling
Models
of Gough,
Elam,
and de Bruyne 4
When
thick,
that is when
(3)
where
of the facesheet,
tc is the thickness
modulus
Young's of infinite
as if it rests would 4 as
theoretically
Crwr,assuming
the core
Poisson's
v c, is zero is
O'wr _ 0.794(E:EcG
c)1/3
(4)
G c is the
core
transverse
shear
modulus.
Note
the
lack
of dependence thickness.
of the core
infinite
with thinner
cores,
(5)
there
is some
interaction in this
between case
the faces
sides and
While
load
properties of t: and
3 and 5, a theoretical
v c is zero is given
in Reference
O'wr
--0.630(E:EcGc)
1/3
(6)
of Gough, would
Elam, appear
facesheets
to be inextensible.
Legget, buckling
and
Hopkins
5 were
the
first
to solve of a strut
the
more
general core
problems
of Their as
wrinkling
of finite
thickness.
flexibilities using
of the core,
well as the stretching be similar attached core axial general short to curves to the middle compressive allowed
this analysis
would
be neglected.
model
Williams,
5 concluded
always
symmetric of Gough,
in sandwich and
isotropic
analysis
de Bruyne 4 produced
estimates
of the
facesheet using
wrinkling
an approach
Hopkins
cores,
the facesheet
O'wr _ 0.760(EIE
cG c)1/3
(7)
Hoff
models and
wrinkling They
for
sandwich
solid
depicted zone
assumed to be
to zero within
w that is chosen
or a value
calculated
using as well
energy
formulation. of the
extensional neglected
as the axial
strain
energy
formulation;
the theory
is related
in Figure
by Hoff based
on the
symmetric reasonable
wrinkling estimate
strut
w < tc/2)
This value
of O'eris given
O'wr =0.910(EIE_Gc)
v3
(8)
S
LIMIT OF CORE SHEAR AND EXTENSION
i
'%
LIMIT OF CORE EXTENSION
7
(2)
(1)
SYMMETRIC WRINKLING ANTISYMMETRIC
WRINKLING
B.I_ey-97-O7/BTI -07
Figure
5. Wrinkling
Models
2.1.2
Sandwich Based
Structures
with
Isotropic Williams
Facesheets 8 reasoned
and Cellular
Cores wrinkling to
that in order in Figure the assumption that the axial were based
sandwich enough
strut
depicted
core-to-facesheet core of
to render
on the assumption
that this
of zero core
stress
is known
anti-plane
assumption
as opposed
with cellular
analysis occur:
following
expression
for the
at which
wrinkling
0.825
](
1/3
(9)
where
ratio of the facesheet. that the form of the wrinkling deformations consistent assumed with by Williams the He anti-plane considered 8 result in
core
stresses
that
violate
the
equilibrium was
equations first
stress both
presented
of a sandwich equations
of finite
thickness. stress
solution
the
anti-plane
to expressions
polynomials
coordinate,
wrinkling
in the symmetric
by Hemp 9 is given
_wr
(10)
Note shown
that
a simple
beam stress
on an elastic
foundation
with
foundation
modulus
2Ec/tcCan
be
to have a buckling
O'wr=O.82Ef
IEctz E/tc
(11)
10
theoretical
value cores.
is suggested After
in Reference
considering stress
wrinkling stress
in sandwich with
cores given
times
associated
wrinkling
in the symmetric
mode
in Equation Norris
the analysis
of Gough,
4 to include section
allowance or a
for an orthotropic, state of plane flatwise wrinkling Goodier facesheets to results neglected. anticipated consideration Goodier could and occur strain.
a state of plane
of the core
Norris,
modulus developed
to the axial
Neou _3 evaluated
stress
by comparing that
assumption
tests indicating
at loads
well below
wrinkling
of wrinkling
Hsu
TM
were
showed
of the sandwich
in a mode
localized based
6 below)
the load
predicted
using
an analysis
harmonic
mode
(1)
ENDS RESTRAINED AGAINST ROTATION-HARMONIC WRINKLING
B .Le y.-97.0 Z,BTI -08
(2)
ENDS FREE TO ROTAqENON HARMONIC WRINKLING MOD E LOCALIZED ATEND S
MODE
Wrinkling
Mode
in Sandwich
by Yusuff
15 are functions
5), over
width
W is calculated
shear
and
strainenergies storedwithin thewidth W to the energystoredin anequivalentextensional spring of modulusK, thenminimizing K with respectto W. For thick cores,W < 0.5tc, and
Cr_r=0.961(EIEcGc) 1/3 (12)
(13)
where
W = 0.5t c (14)
2.1.3
With
Laminated
Composite
focused
on sandwich
structures
facesheets. facesheets
earliest
theoretical by Pearce
wrinkling Hemp's
sandwich and
orthotropic of the
They so that
anti-plane
accounted as well
wavelength symmetric to be
wrinkling
analysis. was
of a sandwich
orthotropic
facesheets
shown
Dllm2+ 2(/)12 +
2066
]{a_
Jib]
2 .t_ O22(l!{a_4"
\m
]_b)
2Ec a2 mZ YgZ t ft c
(15)
where
the D_,
D12 , D22 ,
and
D66
are
laminate
bending
It can be shown
15 reduces facesheets
facesheet
isotropic
bending-extension
TM
extended beams
the subject
analysis
of Pearce They
and
Webber
_6 to study
the stress
facesheet
wrinkling
of sandwich
to bending.
compared
the facesheet
12
side
of the
beam loaded
to cause (both
wrinkling facesheets
to the
wrinkling
on a uniaxially in Reference
stress
calculated stress
using calculated
approximately
16% higher
strut model.
of Gutierrez
facesheets
wrinkling stress by
on the wrinkling
panels
with very
may be very conservative. of an isotropic plane with was strain those the first facesheet elasticity calculated one resting model using on a semi-
infinite,
solid,
isotropic stresses
on an elastic deformation
foundation.
of shear results
in addition beam
stress
generally stress
indicated
the simple
estimates
wrinkling
facesheets
resting
on solid cores.
2.1.4
Summary Theoretical studies performed since wrinkling in Figure 7. 1940 have mode These were yielded equations and used to design on one sandwich of the three
structures mathematical
against
are based
models describing
the Reference
how these
equations
in Table by
1. Generally,
axial compressive
stress
in the facesheet
(16)
for sandwich
isotropic
cores,
or by
(17)
13
THIN THIN PLATE REACTING CORE CORE TRACTIONS FROM CORE ELASTICITY EQUATIONS REACTIN
PLATE G CORE
uc=uo
Z Assumed Deformation Core -_X
_X
TRACTIONS
SOLUTION
._ .................._I w
(1)
GENERAL THIN REACTING PLAFE COR E LIMITED CORE (2) ZONE OF
DEFORMATION
i
x, Oy, Txy = NEGLECT IN PLANE OF CORE STRESSES SOLUTION ELASTICITY B.Ley-9707/BTI -09
(3)
ANTI PLANE CORE
Figure
7. Summary
of Facesheet
Wrinkling
Mathematical
Models
Table
1. Summary
Swr
of Theoretical
Wrinkling
Stresses
For Sandwich
Struts
With
Thick
Cores
0.79(EIEcGc) v3
0.7 6(Ef 0.91(ERE 0.96(El E, G_ y/3 G_ y/3 E c G c )V3
(1)
0.82(EI<G_y/3 0.82E I I_
Eftc
W = t]2
(2)
Yusuff,
15
>
t J2
00.82EI
(2) + (3)
Williams s
(3)
Hemp 9
14
for sandwichstruts with coresfor which the anti-planeassumptionis valid (e.g., honeycomb cores). Equation 17was derived assumingthe core to be in an anti-planestateof stresswhile Equation16wasderivedwithout the anti-planestressassumption.For sandwichstructures with orthotropic laminatedcompositefacesheets, an equivalentmembraneYoung's modulusshould not beused. Rather,E I should
be replaced by
where design
laminate
bending
stiffness
load.
For
constants
have 2.2.
determined
experimentally.
Experimental
2.2
RESULTS to evaluate a highly localized instability stresses 17. failure such as facesheet by tests will empirically to with with cores). that
wrinkling
is a difficult
task.
facesheet
wrinkling
determined In particular,
in Equations
16 and
indicated
in Equations Riddell
16 and 6 (thick
sandwich
and de Bruyne 4 (thin sandwich sandwich predictions, are valid in Table design for structures
suggested
by Yusuff
_5 (thick
with cellular
In order
to correlate
with theoretical
to ensure test.
made more
theoretical
development appears
A list of is the
important
assumptions it is necessary
to either
or refine possible,
to allow in this 2.
of the in
section
be evaluated
of the relevance
of the assumptions
2.2.1
Exhibiting
Reasonably
Good
CorrelationWith
Theoretical
7 tested were
sandwich
panels
and
10.5-in (papreg)
long.
These
made
plastic
presence
of large
15
Table2. ImportantAssumptionsUnderlyingFacesheet Wrinkling TheoreticalPredictions ASSUMPTIONNO. 1 2 ASSUMPTION Thewrinkling load is independent of boundaryconditionsandhasan associated harmonicmodeshape. Symmetricwrinkling failure in sandwichstructures alwaysoccursat loadslower thanthosenecessary to causeantisymmetricwrinkling failure. Thecoreprovidescontinuoussupportto the facesheets. Neitherthe sandwichnor the individual facesheets exhibit shearextension, bending-twisting,or membrane-bending materialcoupling behavior. Any effect of the facesheet-to-core adhesivelayermaybe neglected. Thecoremay betreatedasif it wasattached to the middle surfaceof the facesheets.
3 4
5 6
difficulties encountered in creatinga uniform stateof stressin the specimens.Furthermore,test datausedto determinethe propertiesof the materialsusedto fabricatethe specimens exhibited extremelyhigh scatter. Using the resultsof the 12 valid testsandacknowledgingthe variability in the materialpropertiesof their specimens, Hoff andMautner 7suggested using a valueof k 1 =
0.50 in Equation 16. This value of kl is 34% lower than the theoretical value, k_ = 0.76, suggested bound by Cox and Riddell k_ = 0.63, derived 6 for thick from core sandwich structures, Elam, and 20% lower and than the lower
value,
of Gough,
de Bruyne 4 for thin cores. of the core and a conservative The facesheet
Note that the formula facesheet approximation wrinkling of sandwich Norris and core stress thicknesses; of the given
given
16 is independent 16 represents
in Equation
in 1945, is widely
in the design
structures
with solid cores. hundreds facesheets of sandwich were made struts made of various steel, combinations and glass cloth of facesheet laminates. hard failure
of aluminum, cellulose
acetate,
modes
of failure
during
were:
Elastic facesheet
wrinkling material
of the facesheets
at stresses
below
the
proportional
limit
of the
Core failure
due to "initial
irregularities
in the facesheets"
16
3. 4.
at stresses strength
above failure
facesheets
wrinkling. focus is on wrinkling here. failure, the results in Reference (3) was of tests on specimens 11, the authors poor, of these observed that failed that the
the present
As described in mode
generally several
yielding specimens
unacceptable failed in an
strength. from
In addition,
that failed
in mode
(3) will
not be considered
All failed
specimens
made
and
steel
with
solid
granulated
cork
cores and
in mode
(1)--the determined
for which
between
the theoretically
experimentally wrinkling
of wrinkling
stress extracted
is expected from
stresses
of the mode
Tables
versus
in Figure value
8. The
of a straight to Equation
indicates
of kl to be wrinkling
16, thus
facesheet
Furthermore,
made
of solid
cork lower
core, than
of the
material
to four the
Hence,
it is anticipated stress
anti-plane
are valid
so that
in Equation
of the specimens
indicates
yielding
a second
semi-empirical value
vI = 0) provides
conservative
of the wrinkling
17
0.0025
0.0020
OWR
K I(E
FE CGC )1/3
0.0015
aWR EF 0.0010
0.00050
EF
Figure
0.0025
8. Reference
(1) Failure
- Isotropic
Core Model
Owr
O0 0
0.0020
_ o
-
o o
_q_ "
0.0015 aWR EF
0.0010
0 0
0.00050
I
0 0.00050 0.0010 A_ctF 0.0015 0.0020
I
0.0025
Figure
9. Reference
(1) Failure
-Anti-Plane
Core Model
18
specimens
that failed
in mode
by failure
of the core due to what A plot of the wrinkling 11, normalized by the
(EIE c GJ/3/E_
modulus stresses
theoretical
for a single
used.
With
it can be seen
k_ = 0.50 estimate
k2 = 0.60
by Hoff
and Mautner,
lower
of the wrinkling
to
Similarly, another
Figure estimate
the factor
Equation
17
stress. performed performed by Norris et al. H to struts cores from (EIE _ by from to core
Boller,
VOSS
12
struts
facesheets
3 of Reference in Figure
modulus
appears
Young's
modulus
in Figure that
12 and instead
13 in cases
where
the specimens
of facesheet
wrinkling.
line passing
through 16.
in Figure straight
to Equation
value ofk_
ofk 2 to
to Equation
Figure
= 0.76
unconservative of k_ = 0.63
however,
k 2 = 0.82.
2.2.2
Exhibiting
Generally
Poor
Correlation
With
Theoretical
studies
wrinkling
sandwich
panels
were
and Kuenzi
2 tested
facesheets
on aluminum of 168
They
of 28 different
configurations
for a total
19
0.0040
0.0030
0 0
_'_K 0
I = .63 --
(JWR EF 0.0020
8
0
KI
= .50
0 0.0042
I
0.0044
I
0.0046 (E v EC GC ) 1/3
I
0.0048
!
0.0050
B. Ley-97-08/BTI
-o3
EF
Figure
0.0040
10. Reference
m
(2) Failure
- Isotropic
Core Model
E C tF
// "_
OWR 0.0030
K2EF
_E--_
09/
o'-'/
aWR EF 0.0020
0.0010
0 0
I 0.0010
I 0.0020
I 0.0030
I 0.0040
I 0.0050
B. Ley-97-08/BTI -04
Figure
11. Reference
(2) Failure
- Anti-Plane
Core Model
20
0.0060
o
aWR = K I(E FE CGC)I/3
0.0050
0.0040
aWR EF 0.0030
0.0020
0.0010
EF
Figure
0.0060 --
12. Reference
Core Model
O'WR = K2 EF 0.0050 --
_/E
c tv EF tc
O .,,_q,
._
//
0.0040
--
O'WR EF 0.0030
0.00 20
O
0.0010
I 0.0070
-06aBLey'97"OS/BTI
Figure
13. Reference
Core Model
21
tests.
configurations wrinkling
tested, away
Jenkinson
and
Kuenzi
2 reported
that
only
10
by facesheet
A plot 20,
stresses
of these Young's
from
normalized
of the wrinkling
(E c tf El
of the specimens
normalized slopes
Young's
modulus through
these
indicate Equation
to Equation 2 suggested
14. A differ 2
approximately
as shown and
These
theoretical
values
of k_ = 0.76
k 2 =
Jenkinson
Kuenzi
these
of the facesheets. Norris et al. H and Norris, Boller, flatwise strength, and
Voss
12
that
was a function
data
wrinkling
strength
data
in column flatwise
12 of Table
in the measured
facesheetfor
strengths. one
an average deviation
flatwise
of the test
Flatwise
specimens
attributable available
adhesive
bonding
of modern
processes technology
bondline
strengths.
more
by Harris and
panels
having
0.020cores.
fiberglass were to
on 0.40a more
to 1.00-in-thick semi-empirical
aluminum analysis
reliable
facesheet
waviness
and to
possible
differences direction
in facesheet
wrinkling core
on panels
measured
are presented;
of scatter
in Reference 1 of Reference
Plots
loading
parallel
to the
appear
22
0.0040
0 0
0 0.0010
OWR
K 1 (E FE CGC )1/3
I
0 0.020 0.040 (Er, E c Gc)
B.Ley-97-O8/BTI -07
I
0.060 1/3
!
0.080
EF
Figure
0.0040
14. Reference
Core Model
0.0010
J
J
7
......
K2EF
_ EFtf C
I
0
B. Ley-97-O8/BTI -08
I
0.020 Ec tF EF tc
!
0.030
0.010
Figure
15. Reference
Core Model
23
Figures16and 17,respectively.Straightlines boundingthesedataalsoappearin Figures16and 17. The slopeof a straightline passingthroughthe dataplotted in Figures16 and 17 indicates the appropriate valueof k 2 to be used in Equation 17.
Two bound 0.82. conclusions can be drawn in these stresses from figures Figures 16 and 17. First, both the upper the theoretical and lower ofk 2 = are,
values Second,
value
direction
Since
in core
transverse
G c, Harris by Norris,
21 chose
the equations
derived
semi-empirical stress
the theoretical
mode
the theory
Boller,
_5 equation
of G c. Harris
of the panels
in an antisymmetric
than
of the and
wrinkling
stress
with
honeycomb
better
by Jenkinson
and Crisman2_?
attributed
to manufacturing
be found
by investigating information
Using
the specimen
provided
and Kuenzi,
and Crisman,
2_ the critical
wrinkling
half-wavelengths
can be calculated
(19)
Equation using
19 is derived expression
load
of the specimens 22
the following
4 of Reference
(20)
where
s is the
cell since
size
of the
honeycomb
core
material.
Equation in cell
it does
not take
into account
the walls
the facesheet
to the core.
24
0.020
/ (JWR ---- K2 EF i/EC tF
J
_f
0.015
Owr
EF 0.010
0.0050
I 0 0.020 0.040
Ec tF EF tc
0.060
0.080
21 Test Data-
0.010
(JWR
EF
0.0050
,F
B.Ley-97-08/BTI Y "_F
tF
_"C
21 Test Data-
Equation
19, it was
determined
reported was
and Crisman,
21 the critical
of a single
furthermore,
calculated
using
_,/s < 1.0, use of a furthermore, is very the strong. 17 to with that for the the
Young's between
of an interaction
modes given
in Equation
by Norris,
_,/s appears
predictions and
Hence, and
performed
Kuenzi
2 and
Harris
Crisman
assumption
(3) indicated
1.5
0 0
1.0
0
VALUE THEORETICAL O O []
K 2
[] []
SET
1 psi 2
EC = 68,600 SET
E C = 16,700psi
O'WR
= K2E F
I/
-_-_ EFt C
I 0.5
I 1.0
I 1.5
I 2.0
CRITICAL
WRINKLING
HALF-WAVELENGTH/CORE
CELL
SIZE
Figure
18.
Plot of
k2
Versus
Ratio
of Wrinkling
Half-Wavelength
to
26
possible reported
explanation
of of
by Jenkinson
buckling specimens
in fact failed
This could
in an antisymmetric
of initial why
could load of
explain loaded
in failure load
normal
of specimens of and
parallel shear
modulus.
in Table
in to 0.02-in-thick cores
aluminum
honeycomb
to be used
to validate
the theory
in previous were
studies observed
7, 11, 12, 20, and 21) with immediately of Pearce stable upon 3
only
on the
loaded
objective supported
of the experimental to see if panels tested but below was never occurred
all four
of all four
(symmetric)
directly; occurred
however, close
of local
final
of the The
indeed
stable.
evidence were
Similar
difficulties
encountered
during
24 tested
supported
panels
with
0.024-in-thick
quasi-
glass-polyimide The
along
predictions Note
of wrinkling
Equation
11, appear
facesheet reflect
load same
(7.538 laminate
material
properties stiffness
modulus
of the facesheet
27
than the membranestiffness is 11.405 Msi resulting in a 23% increasein the theoretical wrinkling loads listed in Table 7 of Reference24. Basedon the new theoreticalwrinkling load estimates, it canbeseenfrom Table7 thatthe measured wrinkling loadsareprecisely50%lower thanthe theoreticalpredictions. This is in striking contrastto the resultsreportedby Pearceand Webber. 23 Camarda 24also statesthat all of the specimenswrinkled very close to a supported edge.The authorwentto greatlengthsin designinga testfixture thatwould imposeno rotational restrainton the paneledges.As pointedout in Section2.1.2,Goodierand Hsu TM showed that a
nonsinusoidal formulas mention sinusoidal violated based is made mode wrinkling mode can occur near a supported mode edge under at one-half such support the load predicted conditions. Since by no on the assumption of Goodier was missed. of a sinusoidal
14 work
in Reference likely
24, it is possible
that a critical,
non_4
it is highly in Table 2.
performed
by Camarda
assumption
(1) indicated
and Pfeifer 25 conducted thin laminated communications a loading concluded about condition
a theoretical-experimental composite satellites. absent facesheets Their from value work every
of honeycomb of those
typical
experimental
17 is between
were
to have that
be neglected.
No specific were
specimens laminates
the facesheets
themselves
bending-stretching assumption
coupling. in
the tests
performed
by Bansemir
and Pfeifer
(4) indicated
A summary expressions
described
to the theoretical
for wrinkling
2.3 As
EFFECTS was
OF INITIAL
IMPERFECTIONS 2.2.2, poor correlation loads These has between generally theoretical been are estimates and
mentioned measurement
in Section
of facesheet in the
wrinkling facesheets.
attributed random
to initial in nature;
imperfections
as a linear
wrinkling
mode
shapes shape
are orthogonal.
assumption
made
28
Table
3. Published
Test Data
Correlation
Expressions
Stress
SOURCE Theory
SUGGESTED kl* 0.76 (upper bound) 0.63 (lower bound) 0.50 0.63-0.76 0.50 0.63 0.044
Uncertain
material
properties.
facesheets, cores.
and
Honeycomb
Tests probably violated assumptions (2) and (3) in Table 2. Tests probably violated assumptions (2) and (3) in Table 2. Tests probably violated assumption (1) in Table 2. Tests violated assumption (4) in Table 2.
--
---
wr
wrinkling
load
is the dominant
term
combination. of a sandwich
In
can be reasonably
approximated
w = 6sin(
_x ]
\_,cr]
(21)
(waviness),
wavelength. failure
observations by causing
either failure
(symmetrical
imperfections) strut
(antisymmetrical in the
Yusuff 26)
of a wrinkling expression
mode
is
it has been
shown
that (see, is
the resulting
for the
displacement
of the facesheet
29
6 sin / W-Dcr 1
P where P is the applied _,. load If the and Pcr is the wrinkling facesheet contains load associated with a wrinkling in the _, is given form by mode
(22)
of
of a symmetric
the maximum
facesheet-to-core 2EcW
tc
oz -
2E c 6
k P (Pcr -1) tc
(23)
where
Young's
modulus
and tc is the core thickness. in the form 1, page of an antisymmetric wrinkling mode,
the maximum
core shear
stress
163)
"lJ ....
kT/)
---_--k
tc
)l_-pc-__i_
(24)
kkP
where G c is the core flatwise shear modulus and tI is the facesheet the allowable core thickness. shear
]
If either o-z exceeds panel the will
allowable fail.
stress
or "Vcore exceeds
stress,
the sandwich
of the tested
criticality
in the
27 The
shape author
of the tested
with notches
built into the facesheets. during the fabrication Given be calculated this approach these
depth, then
fabrication. 6 o could
two
parameters ratio.
wavelength,
by a simple given
correlation flatwise
using
strength
be located
attempts
to estimate
6 in Equation
21 based
on direct
of surface
profiles
and by comparisons
predictions.
30
A o
B.Ley-97-O8/BTI -12
.,9------XCR-----._
Figure
19. Typical
Facesheet an
Manufactured for
longer
imperfections. 6 was
sandwich to the
honeycomb i.e.,
Williams
8 suggested
proportional
6 -
KZcr
(25)
where
that
6 was to the
proportional facesheet
proportional
(26)
Note
loads
associated 25 and
variety
of different
modes of failure
Equations
that
wrinkling wrinkling
load may be invalid wavelength form may have for the initial theorized
associated
imperfection
a different The
imperfection
authors
during
cores
be more
the pressure
than more
Consistent
Norris
et al. H found
it convenient
31
(27)
whereF_ is
the flatwise
strength
of the core. on the determination that all three facesheet of the constant forms K 0 by a large number of
25 through Test
27 rely
tests. of theory
results
indicate By assuming
good the
to test.
that
allowable formula
flatwise
stress
wrinkling imperfection:
structures
that includes
the effect
_O'wr
0.82
I_ O/ [ IT-'c rT" \ I
/ Ect:
f_Eftc (28)
1+
kF_to)
Note mode
that
precipitates
fracture
failure
rather
In order
to eliminate to directly
measure
at the wrinkling
to be used number
the result
of tests. from
the largest
available between
reasonable that
theory
however,
it is likely theory
effects
between
wrinkling with
and
and cell
Crisman size
2_
assumed core.
wrinkled
of the 6 o,
They
an empirical
parameter,
28 using Harris
the results
of a large
tests.
allowed
2_ to obtain 2.2.2,
between phenomena,
as discussed
it is possible been
of Equation
responsible made
with theoretical
wrinkling
predictions
17 with
of caution assumption
is warranted of criticality
here.
Expressions
such
in Equation wrinkling
in the
However,
imperfections
in an antisymmetric
not be ignored.
In Equation
that the
core
shear
stress
generated
with
an antisymmetric values
is
proportional
stresses
crimping.
crimping
In reality,
is actually on the
of this shear
imperfection,
below
crimping
load given
in Equation
2.4
EFFECTS A rigorous
OF COMBINED treatment
LOADS of sandwich As discussed the largest panels in Section principal stresses When subjected to combined 1 showed stress loads is
of the wrinkling
suggest using
allowable
wrinkling equations.
them
in interaction
principal equation
compressive,
it is suggested
in Reference
interaction
(29)
where
principal
stress
and the
corresponding
allowable
compressive working
stress with
wrinkling system
with axes
directions,
then applying
the interaction
equation
of the form,
facesheet
( 2
0
q- Oyy
3)
+ _ = 1 (30)
Kow_
where direction
stresses
such
is aligned
with the
wrinkling K=I.0
along is parallel
the core
applied
compressive
to the core
direction,
else K=0.95
33
SECTION CONCLUDING
There predicting can occur have been has been extensive analytical in sandwich and structures
3 REMARKS
experimental subjected mode shape treatment to compressive as shown of the loads. problem Wrinkling 1. Analyses The first of
facesheet either
wrinkling
in Figure
developed involves
treatment applicable
solution made
is most
to sandwich
structures stress
wrinkling
a function on assumed
Young's
functions a small
exponentially
adjacent
stress
problem
anti-plane panel
assumptions.
It is assumed loads
in an anti-plane shown
sandwich that
It has been
sandwich mode
anti-plane
wrinkle
at a lower
symmetric stress
mode. loaded
expression
wrinkling
(symmetric
mode)
sandwich
core is of the
k 2 is a constant. of tests of uniaxially cores loaded sandwich to determine of far fewer and supported to control possible failure and struts with isotropic wrinkling
honeycomb
conducted
facesheet
wrinkling
to combined extremely
since structures
catastrophic impair
failure
in any
of several of the
accurate
mode
in which
experimental shown
been
shown cases,
good
cases, cores
to be surprisingly is generally
honeycomb
are used.
correlation
to small However, of
wrinkling
sandwich
structures
34
honeycombcoresis dueto lack of validity during test of one or more of the basicassumptions, listed in Table 2, that facesheet wrinkling theoryis generallybasedupon. Test resultsreported for sandwichstrutshaving(1) predictedwrinkling wavelengths greaterthanthecorecell size,(2) predictedfacesheetdimpling loads much higher than the predicted wrinkling loads, and (3) predictedantisymmetricbuckling loads much higher than the wrinkling loads generally show correlationwith theoreticalpredictionsto within 20%or better. Initial facesheetimperfections can seriously lower facesheetfailure loads. Analyses incorporating the effects of initial facesheetwaviness generally rely on a parameter that characterizesthe amplitude of this waviness. Attempts have been made to measurethis parameterdirectly; however, the best results have been obtained when the parameter is determinedempirically from a largenumberof destructivetests. Thesetestsshouldincludethe effects of combined loads, especially in the case of sandwich structures with anisotropic facesheets.Rigoroustheoreticalandexperimental treatmentof the wrinkling of sandwichpanels subjectto combinedloadsis conspicuously absentfrom the literature. There aremany waysin which the state-of-the-art of predictingfacesheet wrinkling loads in sandwichstructures may beadvanced.Amongthe effectsthat shouldbe investigatedfurther are:
.
of initial
imperfections loads.
in the symmetric
mode.
2. 3. 4.
coupling. flexibility.
35
SECTION REFERENCES
.
F. J., Sandwich
Construction,
New York,
2.
S. P. and Gere,
J. M., Theory
McGraw-Hill,
Allen, 1969.
H. G., Analysis
and Design
of Structural
Sandwich
Panels,
Pergamon
Press,
Oxford,
Gough,
N. D., "The
Stabilization Aeronautical
of a Thin Society,
Sheet Vol.
by a 44,
Continuous
of the Royal
Williams, Compressive
D.,
D. M. A., ARC
and
Hopkins, Report
H. G., R&M
"Flat
Sandwich
Panels
Under
Technical
1987,
Construction Report
Technical "Buckling
Hoff,
N.J.
of Sandwich-Type
Journal
of the
Aeronautical
.
12, 1945, pp. 285-297. Construction: 2, 1947. of Sandwich Construction," ARC Technical Report R&M A Practical Approach for the Use of Designers,"
"The
Basics
of Bonded
Sandwich
TSB
124, 1986.
W. S., March,
Constructions
12.
Norris,
and
Voss,
A. W.,
"Wrinkling
of the
Facings
of Sandwich 1953. in
Compression," Evaluation
J. N. and Plates,"
Neou, Journal
Compression 664.
of the Aeronautical
Sciences,
36
14.
Goodier, Journal
J. N. and
Hsu,
Buckling
Modes
of Sandwich
Plates,"
Vol. 21, 1954, pp. 525-532. in Sandwich Construction," Journal of the Royal
15.
Yusuff,
of Wrinkling
Vol. 59, 1955, pp. 30-36. J. P. H., "Buckling of Sandwich Panels with Laminated Face
Vol. 23, 1972, pp. 148-160. I. B., "A Theoretical Plates Using Solution for the Buckling System," of Sandwich No. 407,
a Computer
Algebra
Report
Wrinkling of Solids
Sandwich Vol.
Beams
International
Buckling of Applied
of an Elastic Mechanics,
Layer Vol.
Bonded 61,
to
an Elastic pp. 231-235. 20. Jenkinson, Steel 21. Harris, ASCE 22.
1994,
E. W., "Wrinkling
of the Facings
Sandwich
Compression," Wrinkling
of Sandwich
Mechanics
Division, Sandwich
Department Washington,
"Structural
Composites,"
23.
Pearce, With
Buckling
Loads
of Sandwich
Panels
24.
Camarda, Edgewise
25.
Stability
Structures of the
With SAMPE
Thin
Fiber
Proceedings England,
European
26.
Yusuff, Royal
in Sandwich
37
27. Rogers, C.W., "Face Wrinkling As a Function of Surface Waviness," Report ERR-FW-196, GeneralDynamicsCorp., 1964. 28. Wan, C. C., "FaceBuckling andCoreStrengthRequirements in SandwichConstruction,"
Journal of the Aeronautical R. T., Smith, Plates Sciences, Vol. 14, 1947, pp. 531-539. for Structural Stability Analysis of
29.
Sullins, Sandwich
G. W., and
Spier,
and Shells,"
NASA
CR-1457, of Flight
30.
Bruhn,
E. F., Analysis
and Design
Structures,
Tri-State
Offset
Company,
38
REPORT
DOCUMENTATION
PAGE
Form
Approved
Public reporting burden for t_is collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing t_e collection of information. Send comments regardingtNs burden estimate or any other aspect of t_is collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE
Janua_
4. TITLE AND SUBTITLE
1999
3. REPORT
DATES
COVERED NUMBERS
5, FUNDING
Facesheet
Wrinkling
in Sandwich
6. AUTHOR(S)
Robert
P. Ley, Weichuan
7. PERFORMING
ORGANIZATION
NAME(S)
AND
ADDRESS(ES)
Northrop Grumman Corporation Military Aircraft Systems Division One Hornet Way El Segundo, CA 90245-2804
9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES)
National Aeronautics and Space Administration Langley Research Center Hampton, VA 23681-2199
NASA/CR-
1999-208994
12a.
DISTRIBUTION/AVAILABILITY
STATEMENT
12b.
DISTRIBUTION
CODE
Nonstandard
The purpose of this paper is to provide a concise summary facesheet wrinkling mode of failure in sandwich structures. This published research related to facesheet wrinkling. Instead, a smaller provide designers and analysts with a working understanding of the use this survey to guide their judgment when deciding which one of design formulas is applicable to a specific design problem.
of the state-of-the-art for the analysis of the document is not an exhaustive review of the number of key papers are reviewed in order to state-of-the-art. Designers and analysts should a wide variety of available facesheet wrinkling
14. SUBJECT TERMS Sandwich Structures, Facesheet Wrinkling, Honeycomb Core, Foam Core
15. NUMBER
OF
PAGES
49
18. PRICE CODE
A03
17. SECURITY CLASSIFICATICN OF REPORT 18. SECURITY OF THIS CLASSIFICATICN PAGE 19. SECURITY CLASSIFICATION OF ABSTRACT 20. LIMITATION OF ABSTRACT
Unclassified
NSN 7540-01-280-5500
Unclassified
Unclassified
Standard Prescribed 298-102
UL
Form 298 (Rev. 2-89) by ANSI Std. Z-39-18