CompositeStructureDesignAndAnalysis PDF
CompositeStructureDesignAndAnalysis PDF
CompositeStructureDesignAndAnalysis PDF
4.1 General
4.1.1 Overview
© Chemical Industry Press, Beijing and Springer Nature Singapore Pte Ltd. 2018 353
X.-S. Yi et al. (eds.), Composite Materials Engineering, Volume 1,
https://doi.org/10.1007/978-981-10-5696-3_4
354 Z. Shen et al.
In the USA and Europe, development of advanced composite materials began in the
1960s and was first applied in the early 1970s. Fighter jets in service in the 1980s
used advanced composite materials in their wing and tail structures, comprising 20–
30% of the total aircraft weight. The wing structures of the stealth B-2 developed in
the 1980s was made up of *60% composite materials.
Applications of advanced composite materials in civil aircraft proceeded more
cautiously because of safety and cost considerations. To promote confidence in
composite technologies for civil aerospace applications, many programs were
implemented in the USA in the 1970s. These include the Aircraft Energy Efficiency
(ACEE) and Advanced Composite Technology (ACT) programs, and the
Composite Affordability Initiative (CAI). In Europe, the Technology Application to
the Near-Term Business Goals and Objectives (TANGO) project was initiated with
the aim of producing composite wing and fuselage components at a competitive
cost. In the recently launched A380, A350, and Boeing 787, contributions of
advanced composite materials to the structural weight were 22, 52, and 50%,
respectively.
In helicopters, the use of composites has reached 50–60% of the total structure
weight in military helicopters. Advanced composite materials make up 41% of the
structural weight of the US RAH-66. The vertical landing and tilt rotor V22 Osprey
contains up to 51% composite materials by weight and is considered to be a
full-composite vehicle. To date, many small full-composite airplanes have been
launched. Among these, the well-known passenger-goods dual-purpose “ship star”
has passed the airworthiness certification. The world-famous “voyager” has set the
world record for a continuous flight around the world without refueling or landing.
These successes are an excellent showcase for the effectiveness of advanced
composite materials.
In China, research and development of advanced composite materials and their
applications in aircraft structures were initiated in the late 1960s. In the mid-1970s,
the first composite structural components were successfully used in the air duct wall
of aircraft fighter. In 1985, a jet fighter with a composite vertical tail made its first
flight. In 1995, a composite wing with an integrated fuel tank was successfully
developed. This advance marked a new milestone in composite material applica-
tions in aircraft structures in China. Almost all aircraft currently in service now use
composite parts to some extent. In 2000, the commercial airplane-Y7, with a
4 Composite Structure Design and Analysis 355
composite vertical tail, passed the airworthiness certification ushering in a new era
of composite material applications in civil aircraft. Many large-scale commercial
aircrafts currently being planned will feature large proportions of advanced com-
posite materials. Applications of composite materials in helicopters have also
considerably advanced. Imported manufacturing technologies have been replaced
by native design approaches to next generation helicopters in China.
Metal airplane structures usually consist of parts such as the skin, beams, stringers,
rib, and frame. These parts are mechanically assembled with a large quantity of
fastening components. Metal parts can be processed by machining, rolling, forging,
casting, and welding.
For composites, the material and structure can be manufactured simultaneously.
Structural elements can be connected simultaneously with material processing
through co-curing, knitting, braiding, and Z-pin technologies. Thus a large and
integrated structure can be designed and manufactured in one time. The number of
elements and fastening parts as well as the machining and assembling work can be
greatly reduced. Hence, the structure weight and production costs can be consid-
erably reduced. Owing to the above features, integration of design and processing is
a serious concern for composite structures.
composite molding (LCM), such as resin transfer molding (RTM), resin film
infusion (RFI), and vacuum assisted RTM (VARTM). More affordable ACM
structures have received considerable attention for integrated structure design and
manufacture.
The static strength design requirements of metal structures can also be applied to
composite structures. Additional special requirements for composite structures
include: the combined effects of operation temperature and moisture absorption to
362 Z. Shen et al.
determine the allowables for the composite structures. The uncertainty coefficient
used in this case (the original safety coefficient) will remain as 15. The strength of
composite structures is related to the laminating code, geometrical shape, and
applied load. Hence, these factors should be also considered in the determination of
B-allowables. Large numbers of complex structural tests are used to determine
which of the B-allowables of composite structures are not applicable. It is
acceptable to use a specimen that can replicate the structure laminating code,
geometrical shape, and load to determine the B-allowables.
4.2.2.2 Durability
Table 4.3 Low-energy impact damage and durability testing requirements (tool impact)
Zone Damage Damage level Requirement
source
Zone 1 Impactor Impact energy less than 8.1 No functional problems, no need
Easily diameter, in J, or visible damage for structure repair, and no water
impacted 12.7 mm low (energy not below 5.5 J) leakage over two designed life
speed cycles
Vertical to No visible damage caused by a
surface single impact of 5.5 J
Zone 2 Same as zone Impact energy less than No functional problem over two
Not 1 8.1 J, or visible damage designed life cycles, no visible
easily damage, no water leakage after site
impacted repair
Table 4.4 Low-energy impact damage and its durability testing requirement (hail stones and
runway chippings)
Zone Damage source Damage level Requirement
All vertical Hail: Uniformly No functional problem over two
surfaces and Diameter: 20.3 mm distributed designed life cycles, no need for
upward Specific Central distance structure repair
horizontal weight = 0.9 between impact No visible damage
surfaces 27.4 m/s points 20 mm
Vertical to
horizontal surface
45° angle with
vertical surface
Possibly Runway: None No functional problem over two
impacted Diameter: 12.7 mm designed life cycles, no visible
structures Specific weight = 3 damage, no water leakage after
Equivalent to site repair
aircraft speed
Easy to reach from top, i.e., two feet can be placed on the structure.
Note: It is necessary to define the contact zone, position, and weight
related to each of the conditions above.
(3) Honeycomb sandwich structure control
This type of structure should be used in design elements requiring very light
weight, and in parts that can be easily replaced or be accessible for repair. These
structures are not permitted for use in environments exposed to water or for parts
that require immediate replacement. The impact energy level acceptable for this
kind of structure could be appropriately decreased, for example, to the level where
no visible damage 2.5 mm in depth or penetration could be produced by an impact
of 0.5–0.7 J.
(4) Damage sensitive zone and details
Special care should be paid to some damage sensitive zones of the airplane,
including the lower fuselage and the radome, inside the flap lower surface, and
cabin doors. These zones need to be strengthened with thicker structures, possibly
using glass-fiber to replace carbon fiber. Furthermore, the potential for a tire burst
necessitates special attention to be paid to the damage sensitivity of the tire zone.
Similar structures include components near the jet thrust reverser, which are sen-
sitive to damage from ice and debris on the runway.
① Minimum-weight structure: Components such as the radome structure
may not operate if the weight is designed to be too small. Another
example is sandwich structures with low core density. The surface plates
should have the minimum required thickness built into the design to
prevent water invading the core. Surface coatings should not be con-
sidered as the only water resistance measure. Coatings may become
corroded or abraded exposing the component to water.
② Joints of thin skin honeycomb sandwich structures are easily damaged
during assembly and disassembly. Thus, it is necessary to use solid
laminated structures in appropriate joining areas.
③ The rear edge of control surfaces is very sensitive to damage, the area
located 102 mm from the back may be easily damaged by ground
impact, loading and unloading impact, or by lightning strikes. These
components are difficult to repair because of the need to strengthen the
skin and the rear edge. It is acceptable in the design to add a load-bearing
element to resist the load ahead of the rear edge. The rear edge itself or
its surface should be reinforced with a material that is easily repaired,
which will not endanger the functional parts when damaged.
Considering the possibility of cracking leading to sealing problems,
sealing agents should not be used at the ends of components.
④ The edge of laminates should not be directly exposed to an air stream to
avoid delamination.
Possible measures include:
4 Composite Structure Design and Analysis 365
Damage tolerance design principals of metal structures can generally be adopted for
composite aircraft structures. The major source of damage to composite structures is
impact damage. The feature of impact damage is that damage may not be visible on
the component’s surface, but underlying delamination can considerably decrease
the compression strength of the component. Composite structures are defined as
undetectable structures in design. In this section, some particular requirements of
composite structures are outlined [1–3].
(1) Assumption of defect dimensions
The damage tolerance of a component containing a defect should be sufficient that
the component has enough residual strength over the specified operation period.
Here defect refers to both initial and service defects. Defects caused by foreign
low-speed impacts can occur at any time during service, including directly after
being put into service. It is very difficult to detect such defects visibly over the
service time, and thus, this type of defects is defined as an initial defect for the
purposes of damage tolerance analysis and certification.
① Initial defect size assumption: Initial defects may be divided into three
types: impact damage, delamination, and scratches.
a. Impact damage: Considering the features of composite impact
damage, composite structures are defined as undetectable structures
to ensure structure safety. The initial defect size can be determined by
the basic concept defined as barely visible identification (BVID). The
following standards can be used in general cases:
(a) The BVID dent caused by a half-sphere impactor with
25.4-mm diameter (depth less than 2.5 mm).
(b) The damage caused by a half-sphere impactor, with
25.4 mm diameter, at the maximum energy possibly
encountered in service (typically less than 136 J).
The required energy for a BVID dent is related to the structure thick-
ness. The former of the conditions mentioned above is used for medium
thickness structures (less than 6 mm), while the latter is used for thicker
structures (greater than 6 mm). For general cases, the smaller value is
selected from the above two conditions.
366 Z. Shen et al.
The requirements for military aircraft structure design are generally suitable for
civil aircraft with the following differences.
Some items involved in AC 20-107A “Composite Structure” are given below. (The
following has been simplified, the reader can refer to the original.)
5.d Impact damage is generally accommodated by limiting the design strain
level.
6.g It should be shown that impact damage that can be realistically expected
from manufacturing and service (but not more than the established threshold of
detectability) for the selected inspection procedure will not reduce the structural
strength below ultimate load capability.
7. Verification of structure fatigue/damage tolerance: The evaluation of composite
structure should be based on the applicable requirements of FAR 23.571, 23.572,
25.571, 27.571, and 29.571. The following considerations are unique to the use of
composite material systems and should be observed for the method of substantiation
selected by the applicant. When selecting the damage tolerance or safe life approach,
attention should be given to geometry, inspectability, good design practice, and the
type of damage/degradation of the structure under consideration.
(1) Damage tolerance (fail-safe) evaluation.
① Structural details, elements, and subcomponents of critical structural
areas should be tested under repeated loads to define the sensitivity
of the structure to damage growth.
② The extent of initially detectable damage should be established and
be consistent with the inspection techniques employed during
manufacture and in service. Flaw/damage growth data should be
obtained by repeated load cycling of intrinsic flaws or mechanically
introduced damage.
③ The extent of damage for residual strength assessments should be
established. Residual strength evaluation by component or sub-
component testing or by analysis supported by test evidence should
be performed considering that damage.
④ An inspection program should be developed consisting of frequency,
extent, and methods of inspection for inclusion in the maintenance
plan.
⑤ The structure should be able to withstand static loads (considered as
ultimate loads) which are reasonably expected during a completion
of the flight on which damage resulting from obvious discrete
sources occurs (i.e., uncontained engine failures).
368 Z. Shen et al.
verify the adherence of any joint regarded as a key safety flight element and any
joints that have a load-bearing ability not lower than their design limit load:
① Use analysis, testing or both to determine the maximum acceptable
debonding area of each joint under the load in the most serious case. The
detailed design should avoid the possibility of debonding over a greater
area.
② Proof testing should be performed for each end product, i.e., the largest
potential load should be applied to each key adhered joint.
③ The load spectrum and load removal method as well as all other aspects
related to damage tolerance evaluation should be written into testing
programs, and submitted for FAA certification.
④ The B-allowable repeated load verification on flight safety redundant
structures as well as the single-load transferring route structures should
have Failure-Safety capacity. Here, the Failure-Safety capacity refers to
the capacity to withstanding the designed limit load after a major part of
a single-load transferring route structure is damaged.
⑤ Ice hail impact should be included in damage analysis items.
Data sources for mechanical properties of composite materials system and sandwich
composite should be authorized and approved for design purposes. The adopted
data should meet the following requirements:
① Performance characterization based on statistics should be used to
establish design data.
② All factors affecting the required strength, stiffness, and reliability should
be adequately considered in performance characterization, especially the
combined action of moisture and high temperature, as well as the effect
of defects/damage possibly occurring during manufacture or in service
(mainly open hole and low-speed impact).
(b)
Table 4.6 Comparison of commonly used resin matrices and their performance
Performance Resin
Epoxy Phenolic BMI Polyimide High
performance
thermoplastics
Processing ability Excellent Good Good Bad Good
Mechanical Excellent Medium Excellent Good Excellent
property
Heat resistance (°C) Below 130 Below 177 Below 230 Above 288 Above 120
Toughness Good Bad Good Bad Excellent
Dimension stability Excellent Excellent Excellent Excellent Excellent
Cost Low Low Medium High High
dependent. Toughened resin can improve the damage resistance of a composite, and
absorb energy simultaneously.
Composite performances such as heat resistance, aging resistance, fire retardant,
hot/wet behavior, anti-corrosion, and electromagnet performance are mainly
determined by the resin matrix.
Composite processing abilities such as rheological properties, tack and drape, gel
time, prepreg shelf stability, processing temperature, pressure, and time are also
directly dominated by the resin matrix. Common requirements of advanced com-
posite resin matrices are: high strength, good toughness, medium-proofing, low
processing temperature and pressure, long prepreg shelf life, a wide applied pres-
sure window, low cured shrinkage, low toxicity, and suitability for prepreg solution
or melting preparations. The performances of commonly available resin matrices
are given in Table 4.6.
Reinforcing fibers offer composites with high strength and moduli and can also
greatly increase damage resistance. Fiber selection should be optimized through a
combined evaluation of performance and cost. Table 4.7 lists some aspects to
consider for fiber selection. Table 4.8 includes data on environmental effects on the
materials.
Structural processing ability includes the ability to use cure processing and
assembling processing methods. Different processing methods have different
requirements and different structural processing capabilities. Critical issues in
considering structural processing capability will vary depending on the structure, so
it is important to select a proper curing processing method during design. The
curing process of a composite structure can be simultaneously completed with the
4 Composite Structure Design and Analysis 375
material processing. Hence, the design of the structure and its processing are linked.
Structure design should consider the feasibility of processing methods, particularly
for integrated part processing.
Allowables are mainly used to characterize material performance and are defined as
characteristic values of material performance, with a certain confidence and relia-
bility, determined by statistical analysis of coupon testing data under certain loads
and environmental conditions. Design allowables are mainly used in structural
design. These are defined as the design limits determined on the basis of material
allowables and the testing results of typical coupons, elements, and structures.
These allowable may be determined for project requirements and structure integrity,
and also based on previous design and operation experiences. Designers should
follow requirements for structural integrity, such as static strength, stiffness,
durability, and damage tolerance to specify design allowables. The design allow-
ables of a composite system are commonly expressed as strain values. The com-
ponents are subject to the testing verification performed on subcomponents and
full-size scale parts, to guarantee that structures designed by following the design
allowables can meet these requirements [1, 2].
Material allowables are characteristic mechanical properties of material systems,
and mainly used for material selection, acceptance, and equivalence evaluation.
These values are usually not directly used in design except for the modulus, which
needs detailed analysis of open hole, delamination, and impact damage in analytical
models. When performing structural strength verification, design allowables should
be used rather than material allowables.
4 Composite Structure Design and Analysis 377
On the basis of documents available from aircraft companies, for current carbon
fiber resin matrix composites, the allowable strain values under the designed ulti-
mate loads are as follows:
Compression [ec] = 4000 le; tensile [et] = 5500 le; Shear [c] = 7600 le。
Table 4.9 indicates some design allowables of current aircraft structures for
reference.
Currently used design allowables have been lowered, to further reduce weight and
fully utilize the high specific strength and specific stiffness of composite materials [7].
In aircraft primary structures, it is necessary to increase the structural design
allowables, especially the compression design allowables. This approach can be
summarized as:
378 Z. Shen et al.
Owing to the lack of effective analysis and extensive design and use the experience
of composite structures the structural integrity needs to be ensured in some way
[1–3]. The BBA consists of multiple level testing and certification, at the level of
coupons, structural elements, subcomponents, components, and finally a complete
full-scale product. This approach is used to solve (and certify) challenging issues
such as environmental effects and damage. This approach can also reduce the
complexity of full-size tests, qualify the safety of cost reduction measures, and
ensure a pass on full-size verification.
For widely used composite/metal structures, the BBA verification approach can
also be applied for testing below the substructure component level, and used to
complement verification of composite durability and damage tolerance.
For the static strength verification tests of full-size parts, it is difficult to simulate
a full range of combined ambient effects, thus the verification may not be fully
comprehensive. However, the overall integration can be guaranteed for composite
structure static verification, if a BBA verification test program is properly
organized.
On the basis of the maturity of materials and design, and the accumulation of
experience, the testing content and coupon numbers may be reduced to decrease
development costs.
The application of the BBA has not yet been standardized. Relationships
between numbers of specimens and material bases values are well defined for
specimen tests at the lowest level (see part 15), the numbers of specimens required
at higher levels of complexity are somewhat arbitrary and largely based on his-
torical experience, structural criticality, engineering judgment, and economic con-
siderations. Thus, there is currently no standardized methodology for statistically
validating each level of the process, although some attempts have been made to
develop models that relate specimen quantities to overall reliability.
The sensitivity to out-of-plane loads, failure mode multiplicity, sensitivity to
ambient operation conditions, and the absence of mature and reliable analysis for
comparing composite and metal structures, motivate the use BBA for structural
certification of composite structures [1–3].
The multiplicity of potential failure modes is perhaps the main reason that the
BBA is essential for the development of composite structural substantiation. The
380 Z. Shen et al.
many failure modes in composite structures are mainly caused by defects, envi-
ronmental effects, and out-of-plane sensitivities of the materials.
Figure 4.4 shows the building block integration. BBA analysis/testing verification
is generally divided into five steps described as follows [3]:
① Coupon: Small dimension test specimens used for the evaluation of
laminar and laminate performances, as well as general structural char-
acteristics. Test specimens including commonly used laminate strips,
and adhered or mechanically joined strip joints, may be used to gain
material allowables. This approach is also used for evaluation of material
notch sensitivity, ambient effects, and specimen failure modes. At this
level, the number of specimens is usually large, sometimes on the order
of thousands.
② Element: Includes elements and typical structures, i.e., typical
load-bearing units in complex configuration structures, such as skin,
spars, shear panels, laminate and varied joints. Elements also include
weaker components in complex configuration structures such as spe-
cially designed complex configuration joints, mechanical fasteners,
stranger ends, and large inspection ports. The purpose is to verify the
load-bearing ability, failure modes, environmental effects, and analysis
methods for a variety of elements. Design allowables will be generated
at this testing level on the basis of the material allowables, with a large
number of specimens on the order of hundreds.
③ Subcomponent: Large 3-dimensional structures that feature all the
properties of an integral structure section, such as a box segment, frame
segment, wing panel, fuselage panel, wing rib, cabin sector, or frame.
These features all represent a typical section of an integral structure and
their load-bearing ability, environmental effects, damage tolerance,
durability, and analytical methods can be verified. The determination of
design allowables is also generated at this testing level. Fewer specimens
are used at this level, usually of the order of tens, but occasional
hundreds.
④ Component: Refers to the main structural section of a wing, fuselage,
vertical tail, and horizontal stabilizer, which can be used as an integral
aircraft body structure in verification tests. A small number of specimens
will be used at this level, usually one or more. Considering current
application status, verification of structural integrity can be performed at
this level.
⑤ Full-size scale: At this level, full-size scale aircraft structures are used for
analysis and verification and full-size scale static and fatigue tests should
be performed to verify structural integrity and the internal load distri-
bution, deflection, and entire structure failure modes predicted by finite
element analysis (FEA). The requirements for full-size tests should be
proposed by the customer, verification agency, or by the airworthiness
certification body, based on safety and durability requirements. The
numbers of specimens at this level is usually one or two. For
metal/composite combined structures, the main purpose of verification at
this level is to verify the metal structure. Full-size tests can also include
testing content at the component level, which may overlap with lower
levels.
The BBA has been accepted by all composite designers and manufacturers
globally; however, different individuals may use different procedures. On the basis
of the different tasks, verification levels, and test items, the numbers of specimens
may change. Furthermore, the testing sequence may change or overlap. For
example, testing and analysis at the element or subcomponent levels should be
performed in advance, to gain knowledge of potential risks as soon as possible. In
this sense, the current BBA/verification approaches are not yet fully standardized or
documented and certain flexibility is allowed in practical use.
382 Z. Shen et al.
In Boeing 777, 9.9 t composites have been used to make up 11% of the total
structural weight [3]. The horizontal and vertical stabilizers of Boeing 777 are all
made of carbon fiber-reinforced composites. The horizontal tail is a main torque
box main torque box structure consisting of two spar and multi-rib construction.
The BBA for Boeing airplane primary structures (Fig. 4.5) include the following
tests.
(1) Specimens and elements;
(2) Subcomponents(Table 4.10);
(3) Components;
(4) 777 Pre-production horizontal stabilizer test (Table 4.11);
(5) Stabilizer root attachment test;
(6) 777 Horizontal stabilizer tests;
(7) 777 Vertical stabilizer test.
Fig. 4.5 Building block approach for commercial aircraft primary structure
4 Composite Structure Design and Analysis 383
After the determination of laminate total thickness and its local variation based on
the application requirements, ply laminate design should be performed [2, 3, 11,
13]. This step will mainly involve:
① Selection of proper ply angles;
② Determination of the ply ratio at different angles;
③ Provision of proper laminating stacking sequence (LSS).
In addition, ply design also includes local ply alternation, such as the ply design
in connecting zones and at the edges of open hole, as well as the ply transitions for
sudden changes in the thickness of the structure. It may be said that ply design is
one of main factors influencing the characteristics of composite materials. Rational
design will directly affect the composite structural strength, stiffness, stability, and
other important properties such as delamination, damage, failure, and dimensional
stability, as well as processing performance. These factors have direct implications
for load-carrying capacity and operation functions. The critical issues of laminate
ply design are as follows:
(1) A LSS is said to be homogeneous if the ply angles are evenly distributed
throughout the laminate thickness. Homogeneous LSS is recommended for
strength controlled designs.
(2) A laminate is considered to be symmetric if the plies and property parameters
are symmetrical to the mid-plane. A balanced laminate is defined as having
equal numbers of +h and −h plies, where h is measured from the primary
4 Composite Structure Design and Analysis 385
mxy ðlaminate 1Þ mxy ðlaminate 2Þ\0:1
(11) Usually, ply decreasing is used to realize laminate thickness alternation. The
step length of each ply transition zone should not be less than 2.5 mm, or be
followed by l/t > 10 to realize a ply decrease, as illustrated in Fig. 4.6. In
general, there should be no more than two ply numbers simultaneously
decreasing and no ply decreasing is allowed on the width direction of a raised
386 Z. Shen et al.
Laminate stiffness is not only dependent on its thickness, but also related to its
constituent lamina and their stacking sequence. The stiffness of laminate can be
calculated by classical laminate theory [8–10, 12].
(1) Lamina on-axis and off-axis stiffness
Lamina on-axis stiffness (also known as on-axis modulus) refers to the lamina
stiffness along the fiber direction (1). The matrix for laminar on-axis stress–strain is
given as:
8 9 2 1 38 9
< e1 = E1 Ev22 0 < r1 =
6 7
e2 ¼ 4 Ev11 1
0 5 r2 ð4:1Þ
: ; E2
1 : ;
c12 0 0 G
s12
12
where e1, e2, c12 are on-axis strains; r1, r2, s12 are on-axis stresses; E1, E2, m1, G12
are the four independent stiffness coefficients;
v1 v2
¼
E1 E2
where Q11, Q22, Q12, Q21, and Q66 are the lamina on-axis stiffness coefficients:
Q11 ¼ mE1
Q22 ¼ mE2
M ¼ 1=ð1 v1 v2 Þ
Q66 ¼ G12
Lamina off-axis stiffness (also known as the off-axis modulus) refers to the
laminar stiffness along the direction forming an angle h with the fiber direction. The
off-axis stiffness of a laminar at angle h can be expressed by the matrix:
2 3
11
Q 12
Q 16
Q
4
½Q ¼ Q21 22
Q 26 5
Q ð4:3Þ
61
Q 62
Q 66
Q
11 ; Q
where Q 12 ; Q
22 ; Q
16 ; Q
26 and Q
66 are lamina off-axis stiffness coefficients.
is a symmetrical matrix, so that the off-axis stiffness coefficient Q
½Q 12 ¼ Q
21 ;
16 ¼ Q
Q 61 ; Q
26 ¼ Q
62 .
The conversion between the off-axis stiffness coefficient and the on-axis stiffness
coefficient can be expressed by trigonometric functions of the off-axis angle as
given in the table:
11 , Q
② Q 22 , Q
12 , and Q 66 are even functions and thus independent of the
polarity of the angle h, Q 16 and Q
26 are odd functions, and thus
dependent on the polarity of the angle h.
11
③ @@h
Q
¼ 4Q 16 ; @ Q 22 ¼ 4Q
26 :
@h
Hence, only four of the six off-axis stiffness coefficients are independent.
Off-axis stiffness coefficients can also be expressed as a combination of the
on-axis stiffness coefficients as listed in the table:
1 U2 U3
11
Q U1 cos2h cos4h
22
Q U1 −cos2h cos4h
12
Q U4 −cos4h
Q66 U5 −cos4h
16
Q 1 sin4h
2 sin 2h
26
Q 1 −sin4h
2 sin 2h
1
U1 ¼ ð3Q11 þ 3Q22 þ 2Q12 þ 4Q66 Þ
8
1
U2 ¼ ðQ11 Q22 Þ
2
1
U3 ¼ ðQ11 þ Q22 2Q12 4Q66 Þ
8
1
U4 ¼ ðQ11 þ Q22 þ 6Q12 4Q66 Þ
8
1 1
U5 ¼ ðQ11 þ Q22 2Q12 þ 4Q66 Þ ¼ ðU1 U4 Þ
8 2
From the lamina stiffness expression, it can be seen that the off-axis stiffness
coefficient Q ij is a function of the off-axis angle h. Thus, the off-axis angle h can
determine the lamina off-axis stiffness. In general, i.e., when h 6¼ 0° or 90°, the
lamina off-axis coupling stiffness coefficients Q 16 and Q
26 will not be zero. This
means that tensile–shear coupling or compression–shear coupling will exist, as
shown in Fig. 4.7.
The off-axis angle also refers to the ply laminating angle. The design the lamina
off-axis stiffness based on the laminating angle is the main approach to tailoring
laminate properties.
4 Composite Structure Design and Analysis 389
where
{N} laminate mid-plane resultant force array;
{M} laminate mid-plane resultant moment array;
{e0} laminate mid-plane strain array;
{k} laminate mid-plane curvature array;
[A] laminate tensile stiffness matrix (33 matrix), where each tensile stiffness
coefficient can be calculated from the equation:
X
N X
N
Aij ¼ ij Þ ðZk Zk1 Þ ¼
ðQ ij Þ hk l
ðQ ð4:5Þ
k k
1 1
where ðQ ij Þ is the off-axis stiffness coefficient of layer k and hkl is the area of layer k.
k
[B]—laminate coupling stiffness matrix (33 matrix), wherein each coupling
stiffness coefficient can be calculated by the equation:
1XN
ij Z 2 Z 2
XN
ij Zk hk l
Bij ¼ Q k k k1 ¼ Q k
ð4:6Þ
2 1 1
1XN
3 XN
h3k
Dij ¼ Qij k Zk Zk1 ¼
3
Qij k þ Zk h k l ð4:7Þ
2 1 1
12
þ Zk hk l is the inertia moment from the layer k to mid-plane.
h3k
where 12
The strain on laminate cross section {e}Z is expressed as:
t t
fegZ ¼ e0 þ Z fkg Z ð4:8Þ
2 2
The stress of each lamina depends on the lamina strain and stiffness. The stress
in layer k frgk can be expressed as:
4 Composite Structure Design and Analysis 391
fe g ¼ ½ Q
frgk ¼ ½Q fe0 g þ ½Q
Zfkg ðZk1 Z Zk Þ ð4:10Þ
k k k k
It can be seen that the stress in cross section is not constant and may be sum-
marized based on the analysis of stiffness matrices [A], [B], and [D]:
① The expected laminate stiffness and strength can be derived from the
proper selection of each lamina laying-up angles, ply percentages and
stacking sequences, similar to cloth tailoring, and so-called tailoring
design.
② Coupling is a unique feature possessed by laminates, and it is the basis of
aeroelastic tailoring design of aircraft wing surfaces.
Laminate stiffness design analysis generally begins from the matrix [A] and goes
through parameter adjustment (ply sequence, laying-up angle and axis moving) to
generate matrices [D] and [B].
With the use of symmetrical laminating sequences, where each lamina is of the
same material and has the same thickness, then [B] = 0; however, in-plane tensile–
shear coupling (A16 6¼ 0, A26 6¼ 0) and bending–twisting coupling (D16 6¼ 0,
D26 6¼ 0) will exist.
If an asymmetrical laminating sequence is used, when each lamina has the same
material and thickness, then [B] 6¼ 0. Thus, in-plane tensile–shear and out-of-plane
bending–twisting coupling will take place in the laminates.
If a symmetrical and balanced stacking sequence is used, such as [+a/−a], or
[±a]ns, [0/90/±a]s, where each lamina is composed of the same material and has
the same thickness, symmetrical laminating will make [B] = 0. Simultaneously,
balanced laminating will make A16 = A26 = 0. The values of D16 and D26 will
depend on the number of plies (i.e., the more the plies, the smaller the coupling
effect). When the ply number is equal to or greater than 12, D16 and D26 0. In
this type of laminate, there will be no in-plane tensile–shear coupling, and no
out-of-plane bending–twisting coupling. It is also possible to obtain larger in-plane
shear stiffness and out-of-plane twisting stiffness. Hence in structural design,
symmetrical and balanced laminating is usually preferred.
Symmetrical laminating or symmetrical unbalanced laminating may be used to
give unique composite coupling deformation features, and gives great freedom to
composite structural designs, for example, to meet the aerodynamic tailoring
requirements.
(3) Typical laminate stiffness
Typical laminate stiffness is an example of how to describe the designed ability of
laminate stiffness. Common laminates are listed below to illustrate the concept of
laminate stiffness.
① Stiffness of quasi-isotropic laminates: For example, laminates with a ply
sequence of [0/±45/90], [0/±60] are quasi-isotropic laminates. The
stiffness coefficients A11, A22, A12, and A66 in the in-plane stiffness matrix
will correlate as:
392 Z. Shen et al.
1 1
A66 ¼ U5 t ¼ ðU1 U2 Þt ¼ ðA11 A12 Þ
2 2
2 3
A11 A12 0 0 0 B16
6 A12 A22 0 0 0 B26 7
6 7
60 0 A66 B16 B26 0 7
6 7
60 0 B16 D11 D12 0 7
6 7
40 0 B26 D12 D22 0 5
B16 B26 0 0 0 D66
r1 Xt jr1 j Xc
r2 Yt jr2 j Yc ð4:11Þ
js12 j S
The Tsai−Hill criterion can give a smooth and continuous failure envelope,
with small differences between the theoretical values and testing results.
The main problem is that the interaction terms r1 and r2 are only related to
4 Composite Structure Design and Analysis 397
Fij ri rj þ Fi ri ¼ 1; ð4:14Þ
In the equation, there are four quadratic strength parameters and two linear
strength parameters. Of these six material strength parameters, five can be
determined by simple tensile, compression, and shear tests.
1 1 1
F11 ¼ F1 ¼
XX 0 X X0
1 1 1
F22 ¼ F2 ¼
YY 0 Y Y0
1
F66 ¼
S2
Compared with the Tsai–Hill criterion, the Tsai–Wu criterion can give
more information about strength standards, such as the inequality between
tensile and compression strength and how the interaction terms r1 and r2
depend on X, X′, Y, and Y′ values.
Tsai–Wu criterion can also be expressed in terms of strain parameters (refer
to Tsai’s publication for details).
(4) Strength ratio: Strength criterion gives the material failure criterion under
operation stress. To describe the safety margin of materials under propor-
tional loading conditions, the concept of the strength ratio is introduced.
The strength ratio is defined as the ratio between allowable and operation
stress:
398 Z. Shen et al.
ri ðaÞ
R¼
ri
When the strength ratio is used to describe the safety margin under pro-
portional loads, it is assumed that the deformation of a composite is linear
until its final failure.
If R = 1, ri = ri(a), failure occurs.
When R > 1, ri < ri(a), R denotes an operation stress ri smaller than the
allowable stress ri(a). Hence the factor of proportional load increases,
within the safety margin. When R = 2, the operation stress doubles and
failure occurs. Hence, if a load is proportionally applied and doubled
material failure will occur.
When R < 1, rI > ri(a) material failure occurs, or the situation is not
physically possible.
The strength ratio can be also used in strain space. The Tsai–Wu criterion
in strain space can be expressed as:
where
G12 ¼ F11 Q11 Q12 þ F12 ðQ11 Q12 þ Q212 Þ þ F22 Q12 G22
G1 ¼ F1 Q11 þ F2 Q12
G2 ¼ F1 Q12 þ F2 Q22
(4) Determination of laminate ultimate load: Fig. 4.14 shows a frame chart of
the process from ply successive failure to total damage. This is an iterative
calculation process.
Note:
① Ply failure criterion is considered suitable for other laminate
plies.
② Assume that Kirchhoff hypothesis can be always applicable in
laminate successive ply failure.
③ Decide and select stiffness modification models based on the ply
failure modes.
Fig. 4.15 Schematic of rudder canard wing with symmetrical and unbalanced full-height
honeycomb
Fig. 4.16 Conversion of coordinate system during symmetrical and balance skin design
402 Z. Shen et al.
deformation control. Unbalanced ply design can be realized with unequal numbers
of plies with ply angles +h and −h, or by changing the ply angle in balanced
laminates. These choices will depend on the requirements of the design program.
In the above discussion on forward-swept wing design, the focus is on a tech-
nical approach to increase the divergence speed. In fact, aeroelastic tailoring is the
height of a combined design approach, giving many benefits in aeroelastic per-
formances, besides increasing divergence speed. For forward-swept wings, this
approach can also increase buffet speed, improve operational safety, reduce motor
load, and improve the lift-to-drag ratio in aerodynamic quality curves. These per-
formances are closely related to aeroelastic deformation control, reflecting wash-in,
washout, and chord wise deflection, simultaneously. There is often conflict between
these technical approaches, for example, the contradiction between divergence and
buffet speed. In general, wing washout is useful for divergence prevention, while
wash-in is useful to increase buffet speed. Hence aero elastic tailoring should try to
reach a combined optimized design.
404 Z. Shen et al.
may be classed as longitudinal (L) and transverse (W), and have very low
stiffness in the LW plane, i.e., GLW, EL, EW = 0, and give a definite value to ET
and GLT and GWT. The load-bearing conditions are shown in Fig. 4.22, surface
panels will withstand tensile, compression, and shear loads in the xy plane,
where the xy axes are in the same plane as LW. Core materials provide support
to the surface panels and can only withstand transverse shear load and loads
vertical to the xy plane [2].
The failure modes of sandwich structures include: total buckling, surface panel
wrinkling and buckling, surface bending failure, transverse shear failure, local
crash, and impact damage failure. Several failures can occur at the same time
during practical applications, and hence strength corrections should be per-
formed for several failure modes.
406 Z. Shen et al.
Ei1 t1 Ei2 t2 h2 1
Di ¼ þ ðEi1 t13 þ Ei2 t23 Þ ð4:18Þ
ðEi1 t1 þ Ei2 t2 Þk 12k
where k = (1 − mxymyx),
mxy, myx are the panel Poisson’s ratio, and subscripts 1 and 2
denote the top and bottom surface panels. When the panel is
very thin, the second term in the equation can be ignored.
(b) The shear stiffness of sandwich panels: The transverse shear
stiffness of a unit width of panel (N/mm) can be calculated by
the following equation:
408 Z. Shen et al.
U ¼ h2 Gc =tc ð4:19Þ
where Gc is the core shear modulus. In some special cases (Gc = GTL
or GTW)
U ¼ h2 GTW =tc
8
q ¼ pffiffiffi qm ðtm =bÞ ð4:21aÞ
3 3
2tm
qc ¼ q ð4:21bÞ
b m
where
qm —core cell wall material density, kg/m3;
tm —core cell wall material thickness, mm;
b —cell side length, mm.
½rc =rc ¼ 3
sc ¼ V=ðhbÞ ð4:23Þ
(c) Deflection is given by:
2kb PL3 k kb PL
D¼ þ ð4:24aÞ
E f t f h2 b hGc b
or
kb PL3 kb PL
D¼ þ ð4:24bÞ
D hGc b
2Ef tf 2
rcr ¼ ð4:25Þ
S S
(e) The buckling stress on a surface panel is given by:
12
Ec tf
rcr ¼ 0:82Ef ð4:26Þ
Ef tc
Two-point
Fixed-support
Uniform force
0.5P 0.25PL 0.02083 0.25
Simple support
Concentrated force
0.5P 0.125PL 0.00521 0.25
Two-point
Fixed-support
Concentrated force
P 0.5PL 0.125 0.5
(P = qL)
Suspend beam
Uniform force
P PL 0.3333 1
Suspend beam
Concentrated force
P 0.3333PL 0.06666 0.3333
(P = 0.5qL)
Suspend beam
Triangle distribution
0.625P 0.125PL 0.005403 0.07042
(P = qL)
where
{Ai} is the thermal expansion coefficient,
ΔT is the temperature difference.
The subscripts 1–6 of the stiffness coefficient Gij are used for the
axes L, W, T, LW, WT, and TL of the core material, respectively. In
the case of a core special body unit, then: G33 = Ec, G55 = GWT,
G66 = GTL, Ec is the compression modulus, GWT and GTL are the
412 Z. Shen et al.
rz fc ½rc
where
[rc]—core material allowable compression strength,
fc—modification factor, depending on the core partition layer
number as well as loading conditions usually fc = 1.0–2.0,
In the case of a core “single-layer” partition and a normally dis-
tributed load applied to a sandwich panel, then fc 2.0.
Core shear strength can be corrected by the following equation:
Composite body structures should be able to provide at least the same level of crash
safety as metal structures. The design of energy absorption components provides a
basis for materials selection and selection of configuration parameters in crash
absorption structure designs. Thus, composite crash/absorption component design
is an important part of the design process. Anti-crash and energy absorption
components can be classified as follows:
(1) Metal crash absorption components
Structural metal materials such as Al alloy are the tough materials, which can
produce large deformations to absorb crash impact energies. The toughness of
metals can be used to design crash/absorption body structures.
(2) Composite crash/absorption components
Fiber-reinforced composites show linear elastic behaviors at 0° tensile and com-
pression loads with a small failure strain. However, nonlinear ±45° off-axis tensile
and compression loads have a large failure strain and show tough material behavior.
Thus, in composite crash/absorption component design, a design scheme with tube
and wave beams with ±45 plies as the base, and 0° plies as supplements are
preferred, as shown in Fig. 4.24. Tube components may be easily manufactured at
low cost, and the test result analysis is straightforward.
Composite sine-wave web beams are a high stiffness and stability structural
component with both load-bearing and energy absorption abilities. The shapes of
sine waves and ply stacking can be designed to give good processing ability. An
energy absorption component made of a sandwiched web beam with a ladder core
is shown in Fig. 4.25. This structure is used in the helicopter NH90. It absorbing
ability is equivalent to that of a sine-wave web beam, but it can be processed more
easily.
Owing to the dimensional increase in the thickness direction, the stress component
in this direction cannot be ignored. Hence, 3D stress analysis should be performed
on thick laminates. Even under a single in-plane load, thick laminates will also
show a 3D stress distribution. Any stress components reaching a critical state can
result in thick laminate failure. Therefore, the above-mentioned thin laminate 2D
stress analysis and the corresponding failure criteria are not appropriate for thick
laminates.
In addition, the 3D effects in thick laminate composite are more significant than
those in uniform isotropic materials. The strength along the thickness direction is
very low, and has a high sensitivity to matrix cracking and delamination. Thus, it is
necessary to perform 3D stress analysis to establish the failure criterion for thick
laminates. The failure modes governed by the fiber, matrix, and interface should all
be considered.
Many new problems will be encountered in thick composite processing, such as
decreasing residual stress, reducing void content, and ensuring full curing. To
minimize these effects, it is necessary to use special resin matrices, processing
techniques, modes, and curing conditions. Special attention should be given to two
main issues in thick composite processing: A low-level residual stress should be
achieved; the production efficiency should be high, i.e., the time required for full
curing should be as short as possible. Rapid heating and cooling can reduce the
curing time, but can also induce higher residual stress. Slow curing cycles will
result in a low production rate and high cost; however, a fully cured part can be
expected. Cure modeling is very important for thick composite manufacture and can
provide a good understanding of the cure kinetics and instant cure degree in the
cure cycle. This knowledge is useful for predicting the processing stress and is an
important approach to guaranteeing processing quality.
In thick composite laminate analysis and design, it is necessary to understand the
multi-axial strength and stiffness to fully take advantage of thick composites.
Currently, there is a lack of studies on thick composite design, analysis, and
materials testing.
In Fig. 4.30, a flowchart of thick composite analysis is given.
where
m12 m21 m13 m31 m23 m32
¼ ; ¼ ; ¼ ð4:30Þ
E1 E2 E1 E3 E2 E3
420 Z. Shen et al.
Nine independent elastic parameters are involved: E1, E2, E3, G12, G13, G23, m12,
m13, and m23.
If a difference exists between the tensile modulus and compression modulus, the
mean value of the two should be used for small differences; for large differences,
the applied external load should be used for the tensile and compression moduli.
Symmetrical and balanced laminates can be treated as orthotropic laminates, by
changing the subscripts 1, 2, 3 in the above equation to x, y, z coordinate axes to
derive the laminate stress–strain relationship. Hence, the laminate takes the nine
independent elastic parameters: Ex, Ey, Ez, Gxy, Gxz, Gyz, mxy, mxz, myz.
Table 4.13 Typical 3D test data of medium modulus C-fiber/epoxy unidirectional laminates
Item xt E1t e1t xc E1c e1c Yt E2t e2t
Property 1720 114 15,200 1170 114 10,300 55.2 9.65 5700
Item Yc E2c e2c S G12 r12 Zt E3t e3t
Property 207 9.65 21,500 103 6.0 17,000 55.2 9.65 5700
Item zc E3c e3c S13 G13 r13 S23 G23 c23
Property 207 9.65 21,500 82.7 6.0 4000 82.7 3.8 22,000
Note 1. Units: Strength in MPa, modulus in GPa, strain in le; 2. Laminate thickness: >6.35 mm; 3.
Assume transverse isotropy in 2–3 plane
Table 4.14 Typical 3D test data of medium modulus C-fiber/epoxy cross-ply laminates [03/90]n
Item rxt Ext ext rxc Exc exc ryt Eyt eyt
Property 965 103 9330 765 88.9 8600 241 39.0 12,900
Item ryc Eyc eyc Sxy Gxy cxy rzt Ezt ezx
Property 503 39.0 12,900 105 4.8 22,000 23.4 7.72 3040
Item rzc Ezc ezc Sxz Gxz cxz Syz Gyz cyz
Property 414 11.3 3600 28.0 3.7 7700 42.4 4.6 9300
Note 1. Units: same as above; 2. Laminate thickness: 15 mm; 3. Vf = 61.4%, void content 0.04%
422 Z. Shen et al.
E2
E3 ¼ E2 ; G13 ¼ G12 ; n13 ¼ n12 ; G23 ¼ :
2ð1 þ m23 Þ
The failure mode of thin panel structures of composite materials under compressive
or shear loads is an instability known as buckling. Therefore, stability analysis is
required to design these structures [2, 14].
For analysis the structure may be simplified as three components: ① rectangle
laminates; ② stiffened stringers; ③ stiffened laminates.
Rectangular flat plates are widely used in numerous aerospace structures in the
form of unstiffened panels and panels between the stiffened stringer of a stiffened
panel, elements of a stiffened stringer and the skin of the air foil. The bending of the
air foil skin is usually ignored in the analysis. The results of simulations with this
assumption are relatively safe, but not conservative. The bending of fuselage skin
cannot be ignored; however, in this section, only the air foil structure is discussed.
The stiffened stringer is an important component for enforcing the stability of the
air foil skin. Commonly used section configurations include angle-, T-, Z-, I-,
channel-, and hat-shaped stiffened stringers. It can be assumed that there is no shear
load on the stiffened stringer; hence, only the compressive stability needs to be
considered.
The skins of airfoils and the empennage are usually made of stiffened laminates.
Hence these are the most widely studied components in stability analysis. Although
global analysis is highly complex, programs based on the FEM are frequently used
for calculations performed by computer. The performance of a preliminary design
can be estimated by considering the rectangular plate and stiffened stringer
separately.
Next, the stability analysis of three typical components/elements will be
introduced.
buckling load has a closed analytical solution. Therefore, the four edges of
orthotropic laminated plate in supported conditions can be properly simplified: The
ideal conditions are simply supported, fixed supported, and free boundary condi-
tions. Engineers can apply existing closed formulae to calculate the buckling loads.
The calculation of buckling load in unbalanced and asymmetrical laminates is
difficult. Calculations based on numerical methods are commonly used.
The next eight sections introduce calculations used to determine the buckling
load of orthogonal anisotropic laminates under different loads and boundary
conditions.
(1) Uniaxial load, rectangular flat plate with all sides simply supported
In the case of a rectangular flat plate with all sides simply supported and a com-
pressive pressure applied equally to the two edges of a rectangular flat plate
(Fig. 4.34), the formula of the buckling load is:
"
#
p2 b 2 2 a 2 D
22
Nxcr ¼ 2 D11 m þ 2ðD12 þ 2D66 Þ þ ð4:31Þ
b a b m2
In this formula:
Nxcr —axial compressive buckling load per unit length;
m —buckling half-wave number along the x-axis of plate;
a, b —length and width of the plate;
Dij (i, j = 1, 2, 6) —bending stiffness factor of plate.
The parameter m can take the values 1, 2, 3, … in the calculation, to determine a
corresponding set. The minimum value of the set is the buckling load of the
laminate, Nxcr.
(2) Uniaxial load, laminates with loaded edges simply supported and unloaded
edges fixed
The case of a uniaxially loaded plate with the loaded sides simply supported and
unloaded sides fixed is considered in Fig. 4.35. In this case, the calculation formula
of the buckling load is:
8 9
> D b 2 m2 þ 2:67D12 þ >
p < 11 a
2 =
Nxcr ¼ 2 h i ð4:32Þ
b >: 5:33 D22 a 2 þ D66 12 >
;
b m
2p2 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
Nxcr ¼ 2
D11 D22 þ D12 þ 2D66 ð4:33Þ
b
This formula can also be applied under the conditions when the two compressive
pressured edges are fixed.
Supporting experiments have demonstrated that the error of this calculation is
within 10% for a long plate with a width-to-thickness ratio of b/t > 35; however, in
the case of a narrow plate with b/t < 35, the transverse shear effect must be con-
sidered and the calculation results should be revised.
(4) Uniaxial load, long plate with all sides fixed
In the case of a long plate with a length-to-width ratio a/b > 4 and all edges fixed,
when an even compressive pressure is applied to two edges (Fig. 4.37), the
buckling load can be calculated from the formula:
4 Composite Structure Design and Analysis 427
p2 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
Nxcr ¼ 2
4:6 D11 D22 þ 2:67D12 þ 5:33D66 ð4:34Þ
b
This formula can also be applied in the situation of two simply supported loading
edges.
In the case of a narrow flat plate with a width-to-thickness ratio of b/t < 35, it is
also necessary to consider the transverse shear effect and correct the calculation
result.
(5) Uniaxial load, long plate with three edges simply supported and one
unloaded edge free
When an equal compressive pressure is applied to two edges of a long plate with a
length-to-width ratio of a/b > 4 having three edges simply supported and one free
unloaded edge (Fig. 4.38), the calculation to determine the buckling load is:
12D66 p2 D11
Nxcr ¼ þ ð4:35Þ
b2 a2
For a narrow flat plate, with a width-to-thickness ratio b/t < 20, it is also nec-
essary to correct the calculation for the transverse shear effect.
(6) Biaxial load, rectangular flat plate with all edges simply supported
In the case of a rectangular flat plate with all edges simply supported and the short
edges under an equal longitudinal compressive pressure Nx and the long edges
under and equal transverse compressive pressure Ny (Fig. 4.39), the buckling load
can be calculated from the formula:
428 Z. Shen et al.
Nxcr ¼ pb2
2
ð4:36Þ
D11 ðbaÞ m4 þ 2ðD12 þ 2D66 ÞðbaÞ m2 n2 þ D22 n4
4 4
ðbaÞ m2 þ un2
2
Nycr ¼ uNxcr
u ¼ Ny =Nx
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi b pffiffiffi
a¼ D11 D22 =D3 ; b ¼ 4D11 =D22 ; D3 ¼ D12 þ 2D66 :
a
Rx þ R2xy ¼ 1 ð4:38Þ
Rx ¼ Nx =Nxcr
0
; Rxy ¼ Nxy =Nxycr
0
0
Nxcr 0
, Nxycr —buckling load of a laminate with uniaxial and pure shearing buck-
ling load, respectively.
The safety margin of shearing and compression of laminates MS can be calcu-
lated from the method illustrated in Fig. 4.44.
ON
MS ¼ 1
OM
A flat laminate with all edges supported, which has already undergone buckling
may still have some load-bearing ability, which is known as post-buckling strength.
This parameter must be considered for supporting components/elements, and more
details on post-buckling strength will be discussed in Sects. 4.7.1.2 and 4.7.2.3 of
this chapter.
2p2 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
Nxcr ¼ D11 D12 þ D12 þ 2D66 ð4:39Þ
b2
The buckling load of the flange of long lath with one free edge is given by:
12D66 p2 D11
Nxcr ¼ þ ð4:40Þ
b2 L2
In these formula, the length of the stiffened stringer L is used instead of the
length of the lath a, which appeared in the initial formula.
Equations (4.39) and (4.40) are only suitably for use with orthogonal anisotropic
laminates and do not consider the effects of bending and torsion rigidity of lami-
nates, i.e., D16 and D26. However, these two formulae are correct for most
Nxcrk
rstcr ¼ st
Exc ð4:42Þ
tk Exck
where:
Nxcri buckling load of ith plate element;
Nxcrk buckling load of kth plate element;
rcri buckling stress of ith plate element;
rcrk buckling stress of kth plate element;
Exck equivalent compressive modulus along x-axial of kth plate element.
1 A212k
Exck ¼ A11k
tk A22k
st
Exc —equivalent compressive modulus along x-axial of stiffened stringer;
n
P
A2
A11i A12i
22i
bi
st
Exc ¼ i¼1 P
n
bi t i
i¼1
Fig. 4.50 Crippling curve of plate element with one free edge
The crippling curves in Figs. 4.51 and 4.52 are corrected crippling curves of
two kinds of laminates with one free edge and no free edge, respectively.
The y-coordinate is the dimensionless crippling stress rrcucc EExc ; the x-coordinate
qffiffiffiffiffiffiffiffiffiffiffiffiffi
is the dimensionless ratio of width to thickness t Exc pffiffiffiffiffiffiffiffiffi , and
b E r cu
Exc Eyc
¼ 12D11 1 mxy myx
E ð4:43Þ
t 3
436 Z. Shen et al.
where
D11 —bending stiffness coefficient of laminate;
mxy, myx —Poisson’s ratio of the laminate;
Exc, Eyc —equivalent longitudinal and transverse compressive moduli of
laminate;
b —laminate width;
t —laminate thickness;
rcu —compressive strength (ultimate);
rcc —crippling strength (stress).
The values of Ex ðExc Þ, Ey ðEyc Þ, mxy and myx of a symmetric laminated plate can
be calculated from the following formulae.
1 A2
Ex ¼ A11 12
t A22
4 Composite Structure Design and Analysis 437
1 A212
Ey ¼ A22 ð4:44Þ
t A11
A12
mxy ¼
A22
A12
myx ¼
A11
rcu ¼ Exc ec
where
Exc —equivalent longitudinal modulus of elasticity of laminate;
ec —compressive strain design allowable value of composite laminate.
④ Apply the following formula to determine the weighted contribution of
crippling stress of each plate element of the composite laminated plate.
The crippling stress of a stiffened stringer can be determined from:
P
N
rcci bi ti
i¼1
rstcc ¼ ð4:45Þ
PN
bi t i
i¼1
where
rcci —crippling stress of ith stringer;
bi —width of ith stringer;
ti —thickness of ith stringer;
N —number of stringer.
438 Z. Shen et al.
Note: if the value rcc of one plate element is higher than its rcu , the
whole calculation should use rcu instead.
(3) Some issues should be considered when calculating the crippling stress of a
stiffened stringer.
① The data in Figs. 4.51 and 4.52 are taken from experiments based on
a stiffened stringer with uniform thickness. If the thickness of a
single plate element in the stiffened stringer differs greatly, the
thicker plate elements will give greater resistance than the thinner
plate elements and enhance the buckling and crippling stress of the
thinner plate elements. However, the buckling and crippling stress of
the thicker plate elements will also be reduced. Therefore, the
crippling stress of the affected plate element should be modified. The
crippling stress of a stiffened stringer depends on the plate elements
that will undergo buckling or crippling first.
② Consideration of fillets: As shown in Fig. 4.53, a 0°-material is used
to fill the corners of stiffened stringers, and these materials can boost
the crippling rigidity of the stiffened stringer. The area of the fillets
under pressure is directly proportional to the square of the corner
radius. Thus, the larger the corner radius, the greater the enhance-
ment on the crippling rigidity. The following formula can be used to
estimate the enhanced crippling stress:
0 1
1 þ PfE fb t
EA
stcc ¼ @ Arst
i i i
r ð4:46Þ
1 þ P fb t
A cc
i i
where
rstcc crippling stress of stiffeners (without fillets);
stcc
r crippling stress of stiffeners (with fillets);
Af cross-sectional area of fillets;
Ef equivalent longitudinal modulus of elasticity of fillets.
where (EAÞst and ðEIÞst are the tensile or compression rigidity and the
bending rigidity, respectively. These two values may be calculated
from Eqs. (4.53) and (4.54), given in Sect. 4.7.2.3 of this chapter.
The critical stress of a stiffened stringer is:
" 0
2 #
rcc L
rcr ¼ rcc 1 2 ð4:48Þ
4p Exc q
If the value of L0 =q is greater than 12, this formula may require some
modifications.
④ In Fig. 4.54, the broken line represents the calculated initial buckling
stress compared with experimental results from plate elements with
one free edge and with no free edges. The solid lines show crippling
stress data from corresponding experiments. The calculated initial
buckling stress is smaller than the experimental values when b/t is
great, i.e., in the case of a thin plate. The value of the calculated
initial buckling stress is larger than the experimental value when b/
Stiffened stringers are a typical component used in airfoil structures. The stability of
a stringer is enhanced when it is reinforced with a cover. Part of the stiffened
stringer between two wing ribs and two wing spars should be analyzed to consider
the stability of the whole design. For convenience, the structure can be simplified as
a set of parallel stringers and the dimensions and materials (i.e., layering) of a
section plane of the stringers are same, with equal spacing. The slight lateral cur-
vature of a stiffened stringer may be ignored, as shown in Fig. 4.55. This approach
is widely accepted by engineers.
The loading situations of stiffened stringers can be divided into three categories:
axial compression loading (along the length direction of the stiffened stringer),
shear loading, and combinations of shear and compression loadings.
Failure modes of instability can be divided into four categories: ① covers
between stringers or parts of the stringers buckling; ② general instabilities of the
stiffened stringer. A long stiffened stringer can be considered as a wide column with
the use of Euler instability under axial compression; ③ crippling damage that may
occur in a short stiffened stringer under axial pressure; ④ a combination of the
previously mentioned Modes ① and ②.
Structural stability depends on the rigidity of the structure and the rigidity of the
support conditions (i.e., the boundary supporting conditions). The structural sta-
bility of laminated plate is closely related to the layering order. Hence, it is nec-
essary to consider the influence of layering order on stability in the design phase.
(1) Influence of layering order on buckling of laminated plate
The buckling load of laminated plate is related to layering order, loading
environment, geometric dimensions, and boundary supporting conditions.
Thus, there are no general rules for setting the best layering order of a laminated
plate. Specific analyses are needed for specific loading situations, geometric
dimensions, and boundary conditions.
To enhance the buckling load of a laminated plate, the following observations
may help guide layering design:
(1) Symmetrical and balanced laminated layering are adopted in most cases,
except for situations with special requirements, such as requirements for
aeroelastic tailoring. To avoid plate deflection caused by coupling of
flexural tension and bending, let Bij ¼ 0, D16 0 and D26 0: This
deflection is equal to the amount of initial deflection of a laminated plate
and it will decrease the buckling load.
(2) For a rectangular laminated plate that is under pressure along its length, a
higher buckling load may be achieved when ±45° plies are layered on the
surface of the laminated plate.
(3) For a rectangular laminated plate that is under pressure along its width, a
higher buckling load may be achieved when 0° plies are layered on the
surface of the laminated plate.
(4) The maximum buckling load of a laminated plate under a given shear stress
is achieved when ±45° plies are layered on the surface of the laminated
plate. The buckling load value of a plate under positive shear stress is
lower, than that of a plate under negative shear stress, as shown in
Fig. 4.56. This effect is attributed to the D16 and D26 values of the plate.
The buckling loads of an orthotropic plate are the same no matter if the
plate is under positive or negative shear stress.
Fig. 4.56 Rules for shear plates under positive and negative shear stress
4 Composite Structure Design and Analysis 443
Fig. 4.57 Illustration of the 45°-surface fibers along the compression direction under combined
shear stress
For symmetrical laminated plates, where D16 6¼ 0 and D26 6¼ 0, the fibers
on the outer surfaces may allow for a higher buckling load when the plate is
under a combination of shear stress in the pressure direction, as illustrated
in Fig. 4.57.
(5) The behavior of a laminated plate under combined stress from pressure and
shear loading is an unusual situation because of the effects of D16 and D26.
In Sect. 4.1.1 of this chapter, a pressure and shear stress formula is pre-
sented for orthotropic laminated plates (D16 = 0, D26 = 0):
Rx þ R2xy ¼ 1
Thus, a parabola may be defined in this the coordinate plane, to describe the
pressure load-to-shear load ratio with Rx and Rxy as coordinates.
A related parabola for shear buckling of a symmetrical laminated plate under
pressure and shear stress loading can also be defined. When D16 6¼ 0 and
D26 6¼ 0, the parabola may be distorted becoming more prominent or concave
according to the different direction of the shear stress (i.e., positive or negative
shear stress). Figure 4.58 shows buckling curves of a symmetrical laminated
plate, with D16 > 0 and D26 > 0, undergoing combined pressure and shear
stress loading in the positive and negative directions. As shown in the figure:
① In the case of D16 > 0 and D26 > 0, a negative shear stress
makes the parabola more prominent, indicating an enhancement
of the buckling load under axial pressure. A positive shear stress
makes the parabola concave and reduces the buckling load under
axial pressure. Furthermore, higher values of D16 and D26 have a
more prominent effect on the concaving of the parabola.
② In the case of D16 < 0 and D26 < 0, the influences of negative or
positive shear stress on the buckling curves have the opposite
effect. Thus, negative shear stress makes the parabola concave
and decreases the buckling load under axial pressure; however, a
444 Z. Shen et al.
The classic theory of linear buckling has been used to analyze the stability of
structures in engineering. According to this theory, when a structure has achieved
the critical state of initial buckling, its normal deformation (deflection) suddenly
increases arbitrarily. This means that the structure loses its load-bearing capacity. In
practice, when the skin of a thin-walled stiffened structure of a plane features local
buckling, the structure generally maintains the ability to bear load, which is known
as post-buckling strength. For structures designed according to their initial local
buckling stress as the limiting allowable stress, the post-buckling strength of the
structure is not used. Thus, the potential load-bearing capability of a structure is not
fully accounted for [3, 16].
To explain the differences between the practical stabilities of structures and the
stability calculated based on the theory of linear buckling, nonlinear large deflection
buckling theory has been proposed. This theory is based on in-depth theoretical and
experimental studies of post-buckling behavior of structures.
446 Z. Shen et al.
In this section, nonlinear large deflection buckling theory and linear buckling theory
are compared in terms of analysis, processing, and the solutions derived. The basic
concepts and features of post-buckling issues are introduced.
(1) The post-buckling problem involves analysis of a structure from initial buckling
to damage and failure.
Linear buckling theory analysis indicates that when a structure has achieved the
critical state of initial buckling, its deformation (deflection) increases arbitrarily,
and the load-bearing capacity is suddenly lost. It necessary to determine the
load and buckling mode of the initial buckling of a structure.
Nonlinear post-buckling theory can be used to solve the deformation and forces
acting on a structure from the initial buckling to damage to the failure. This
approach involves both stability analysis and requires judgement of the failure
related to the intensity of the damage. Thus, analysis of post-buckling unifies
the analysis of the stability and the issue of strength. In the analysis, many
factors that affect the stability and strength of the structure should be consid-
ered, including: the impact of damage, initial defects, temperature and
humidity, and guidelines of material damage.
(2) In the analysis of post-buckling of a structure, the impact of a large deformation
needs to be considered to accurately describe the state and strength
characteristics.
4 Composite Structure Design and Analysis 447
Linear buckling theory analysis establishes the equilibrium equations for the
initial position and shape of the structure and, therefore, does not reflect the
impact of structural deformation on the equilibrium state.
In practice, a structure under a load undergoes some deformation. After the
initial post-buckling deformation, the structure will enter a buckled state.
Analysis by theory of nonlinear buckling considers the structure and processes
that might change the position and shape of the structure from their equilibrium
values. This analysis allows for a more accurate description of the structure and
the forces acting on it and can more truly reflect the characteristics of the
system.
(3) Post-buckling analysis of a structure is calculated from progressive data sets in
the moment after a load is applied to the structure. The structure’s stiffness after
deformation as well as changes in its position and shape are recalculated in
iterations. This analysis can determine whether a process will undermine the
strength of a structure but requires an understanding of the structure and the
acting forces.
(4) The use of FEMs for linear buckling analysis of a structure can be reduced to
solving a set of linear algebraic equations equal to zero for the determinant of a
coefficient matrix of the eigenvalue problem.
The FEM and nonlinear buckling analysis require the solution of the nonlinear
algebraic equations in repeated iterations. Accurate calculations and conver-
gence are not always achieved. Thus, nonlinear analysis calculations are a
specialized research field.
Fig. 4.61 Damage path in clamped laminated square plates composed of three layers with axis
15°, 45°, 30°
4 Composite Structure Design and Analysis 449
while the dashed line shows a calculation considering the path of an internal
damage. The laminate may undergo gradual failure indicated by the dashed line
showing a gradual downward trend or a more sudden failure indicated by a sharp
downward trend. Analysis of the former case suggests that failure occurred owing
to tension caused by destruction of fibers. The latter case reflects tension (pressure)
caused by the destruction of the matrix.
Figure 4.62 shows the load–deflection curve (Nx−wc) of a vertical reinforced
composite laminate material and skin layer, clamped at both ends under axial
compression with two simply supported edges (taken from a NASA report). The
focal points for deflection of the skin map, respectively, are given for a thin mesh
(solid line), a dense grid (dashed lines), and theoretical calculations and experi-
mental measurement points (triangles).
There is clearly a large difference between the theoretical values and the test
results. In the theoretical analysis and experimental measurements of the reinforced
laminates, the presence of geometric defects and internal damage, or improper
handling of boundary conditions will cause errors in the results of theoretical
calculations and experimental measurements.
be ¼ ub ð4:49Þ
where be is the width and u is the effective width coefficient, determined from
experimental data.
A relationship for estimating the post-buckling load-bearing capacity of a stiff-
ened plate of a given width is presented in subsection (3).
(2) Estimation of reinforced post-buckling load-bearing capacity:
A reinforced Be in the buckling and pressure loss analysis may be divided into
two types of stiffened plates. FEMs and experimental studies of the two types of
plate elements have been used to study the buckling pressure loss after destruction
in pressure loss curves. The Be of components in a reinforced plate element under
pressure loss can be considered to be a stress-weighted sum of estimates of the
post-buckling load-bearing capacity.
(3) Estimation of stiffened panel post-buckling load-bearing capacity:
Here, two pilot projects based on this estimation method are introduced.
① The subsection approach used for metal plates can be applied to com-
posites subject to axial compression. A long board is divided into shorter
board panels based on the slenderness ratio (L′/q). Figure 4.64 illustrates
three regimes for division of boards.
pffiffiffiffi
In a stiffened panel L0 ¼ L= C for an effective column length C, where the end
of the stiffened plate support profile or q factor can take C = 1–4, although it is
generally assumed that C = 2.0. The value of q for a stiffened plate radius of
gyration can be determined by the following equation:
sffiffiffiffiffiffiffiffiffiffi
ðEIÞ
q¼
ðEAÞ
where (EA) and (EI) are the stiffened plate tensile (compression) stiffness and
bending stiffness, respectively, according to Eqs. 4.53 and 4.54.
In the pilot study and mechanical analysis:
(a) In the D-E section the short-board features pressure damage, where 0 < L’/
q 20.
(b) In the B-A section a long board features damage leading to overall instability.
(c) In the D-B section of medium and long boards, before damage to the rein-
forcement between the skins local buckling occurs first. Thus, it is necessary to
account for the post-buckling load-bearing capacity.
The D-E between each section, where L’/q = 20, can be defined separately for
the D-B sections and the B-A cutoff points between sections. For B the skin between
the reinforced parts determines the initial buckling stress.
The actual structures of a stiffened panel include medium and long boards and
stiffened plates. Thus, these are the focus of post-buckling load-bearing capacity
analysis. Test results show that in the D-B section of a stiffened plate, the
post-buckling load-bearing capacity and average failure stress can be fitted by a
parabola. The vertex of the parabola is D, the other point is B. This allows esti-
mation of the post-buckling load-bearing capacity of reinforced pressed plates from
the equation:
rcr rcr
co ¼ 1 1
r cc
r ð4:50Þ
cc rr
r
where
co stiffened panel average failure stress;
r
cc Type of short stiffened plate (0 < L’/q 20) average pressure loss of the
r
failure stress;
rcr reinforcement between the skin of the initial local buckling stress;
rr A factor to discount the skin or be reinforcement after the effects of local
buckling decreases the stiffness. In the calculation of the overall instability of
stiffened plate stress, for be reinforcement of more than 4, the system can be
considered a side support with the width determined by the Euler column
formula.
In preliminary design, the following simplified formula are used for preliminary
estimates
co ¼ r
r cc ð x Þ
rcc Þ2 ðL0 =qÞ2 =ð4p2 E ð4:51Þ
where
x stiffened plate x direction equivalent elastic modulus;
E
q stiffened plate section radius of gyration;
X
m
A¼ bi t i ð4:52Þ
i¼1
Xm
A2
ðEAÞ ¼ A11i 12i bi ð4:53Þ
i¼1
A222i
Xm
A2 D2
ðEIÞ ¼ A11i 12i bi ðzi zc Þ2 þ D11i 12i bi ð4:54Þ
i¼1
A22i D22i
where
bi first plate element of width i;
ti thickness of the ith plate element;
A11i, A12i, A22i first plate element of the ith plane stiffness coefficient;
D11i, D12i, D22i first plate element i of the bending stiffness coefficient;
(zi−zc) first section i of a plate element on the neutral axis from the
center;
Zc stiffened plate section on the neutral axis position (from
calculation of the distance between the reference axis).
P
m
Exi bi ti zi
zc ¼ i¼1
Pm
Exi bi ti
i¼1
where
Exi first plate element i in x direction of the equivalent modulus of elasticity;
zi part i of a plate element calculation of the reference section of the center
distance from the axis;
In general, Eq. (4.54) is used, when the second part is negligible compared with
the first; see “Stability Analysis of Composite Structures Guide” in Appendix A.
② The effective width method for skin damage occurring prior to local
buckling of a stiffened panel, can be used to estimate the post-buckling
load-bearing capacity of the skin. For be reinforcement, the buckling
454 Z. Shen et al.
where
P stiffened plate load damage;
Exs direction of the skin equivalent elastic modulus;
Exst reinforced
be equivalent
x direction modulus of elasticity,
A2
Exs ¼ 1t A11 A12
22
;
t thickness of skin;
F reinforced area profiles;
eb reinforcement of the buckling pressure loss or strain;
be effective width of skin.
A11, A12, A22 skin stiffness coefficient of the plane.
When the computation can be divided into m articles reinforcing a symmetric
laminated plate element, the following equations may be used:
m
1X A212i
Exst ¼ A11i bi
F i¼1 A22i
X
m
F¼ bi t i
i¼1
where
bi reinforcement be of the first plate element of width i;
ti reinforcement be of the first i of the thickness of a plate element;
A11i, A12i, A22i reinforced articles in the first i-plane of the plate element stiffness
coefficient.
The be reinforcement between the skin of the effective width can be determined
by the following equations:
be ¼ ub
u ¼ n þ ð1 nÞescr =eb
4 Composite Structure Design and Analysis 455
a4
n ¼ 1 2= 3 þ g ð4:56Þ
b
where
a, b be reinforcement between the length and width of skin;
u be reinforcement of a skin with the effective width coefficient;
escr local buckling of the strain skin;
eb reinforcement of the buckling strain;
η anisotropy degree of the skin, η = A22/A11.
In addition, the articles reinforcing the effective width between the skin can also
be determined from the following equation:
b rscr
be ¼ 1þ f ð4:57Þ
2 rcc
where
be reinforcement between the effective width of the skin;
b reinforcement between the width of the skin;
rscr local buckling stress of the skin;
rfcc pressure loss stress of skin attached to the end of the reinforced section.
If Eq. (4.57) is used, to estimate the damage to a stiffened plate, be can be
determined from eb of the load P, with Eq. (4.55) where the response of the
pressure loss is given by:
uA ¼ vA ¼ 0; vB ¼ 0 ðor uD ¼ 0Þ:
Other supported boundary conditions can be used with this method. In the case
of a support for an elastic boundary, the corresponding w and hxi, hyi values are
given by the stiffness of the elastic support. For the sides of a fixed supported plate,
hxi or hyi is 0; for the sides of simply supported plate wi = 0 or replaced by the
stiffness of the elastic support.
③ Loading, as shown in Fig. 4.66. Put in-plane load Nx, Ny, Nxy to the
nodes of each side, with upper and lower points corresponding to the
same node. If the load changes along the edge, the load of each node can
be not same.
④ Critical buckling load: Nicr = kminNi
where kmin—minimum eigenvalue;
Ni—stress of analysis with design load.
(2) Local buckling analysis of laminate
In the local buckling analysis of laminates, the panel can be considered to be a
beam support for the core, or flexibility base. The flexibility base has the bending
stiffness and shear stiffness of the core.
Local buckling failure modes of laminates can be divided into three types:
damage to laminated panels, damage to the core, and damage to the interface.
Failure modes of laminated panels include: single-layer instability, folding of the
intergrid, laminated panel buckling.
Failure modes of the core: sandwich core crush, shear failure of sandwich core;
interface damage of the sandwich core and panel debonding from the core.
The calculation methods of various failure modes are introduced as follows:
4 Composite Structure Design and Analysis 457
① Single-layer instability
rjcr nj (2 nj ) nj \1
¼ j ¼ x; y; xy ð4:59Þ
Gz 1 nj [ 1
where
rjcr buckling stress of single-layer, j = x, y, xy, compressive buckling stress and
shear buckling stress along the x- and y-axis. respectively;
Gz interlaminar shear modulus;
pffiffiffiffiffiffi
nj DBj
stiffness parameter, nj ¼ Sj , j = x, y, xy;
Bj base stiffness, Bj ¼ LbEt0 , j = x, y, xy;
z
E22C
Lx ¼ SPL0 þ SPL9 þ 0:5 SPL45
E11C
E22C
1 þ 2:5
E11C
E22C
Ly ¼ SPL9 þ SPL0 þ 0:5 SPL45
E11C
E22C
1 þ 2:5
E11C
E22C
Lxy ¼ SPL45 1 þ þ 0:25
E11C
E22C
ðSPL0 þ SPL9Þ 1 þ 2:5
E11C
where
458 Z. Shen et al.
2Ec
B¼ ð4:60Þ
tc
Dx ¼ D11
Dy ¼ D22
Shear stiffness Sj , j = x, y, xy
D211
Sx ¼
DEN1
D222
Sy ¼ ð4:62Þ
DEN2
D2xy
Sxy ¼ pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
DEN1 DEN2
where
X
N
DEN1 ¼ CAðiÞ tðiÞ=G13
1
4 Composite Structure Design and Analysis 459
X
N
DEN2 ¼ CBðiÞ tðiÞ=G13
1
X
i
CA(i) ¼ Q11 (k) AB2(k) t(k)
1
imax ¼ N; k ¼ 1; 2; . . .i
X
i
CB(i) ¼ Q22 (k) AB2(k) t(k)
1
imax ¼ N; k ¼ 1; 2; . . .i
X
k1
AB2(k) ¼ 0:5 h þ ðtf t(i) 0:5 t(k))
1
k ¼ 1; 2; . . .N
Ef0 tf 2
rcr ¼ 2 ð4:63Þ
k Sc
pffiffiffiffiffiffiffiffiffiffiffiffiffi
where Ef0 ¼ E1f E2f
k ¼ 1 t12 t21
Vx
rxcr ¼ pffiffi ð4:64Þ
d0 B7 4B=D11
1þ hs13b
Vy
rycr ¼ p ð4:65Þ
4 ffiffiffiffiffiffiffiffiffi
d0 B8 B/ D22
1þ hs23b
Vxy
rxycr ¼ pffiffiffiffiffiffiffiffiffiffi p ffiffiffiffiffiffiffiffiffi ð4:66Þ
d0 B7 B8 4 B/Dxy
1þ pffiffiffiffiffiffiffiffiffiffiffiffiffiffi
h s13b s23b
B7 ¼ DD211 =DDEN1
B8 ¼ DD222 =DDEN2
where DD11 and DD22 are equivalent conversion bending stiffness of the composite
panel, the equivalent conversion formula is:
Z 00:5tf
ðKÞ
Aij ¼ Qij dz i,j ¼ 1, 2, 6
00:5tf
where
h
DDEN1 ¼ 2 DEN1 þ CA(N)2
GC13
h
DDEN2 ¼ 2 DEN2 þ CB(N)2
GC23
Vxy — shear; (i.e., the lowest value of the following: overall critical buckling
stress of panel, critical buckling stress of layer, critical core-wise buckling
stress of panel);
④ Sandwich core crush
Vj
rjcr ¼ ; j ¼ x; y; xy ð4:67Þ
1 þ B rdcc0
where
rcc compressive strength of sandwich core;
d0, B, Vj, j = x, y, xy see definition above.
Vj
rjcr ¼ ; j ¼ x; y; xy ð4:68Þ
1 þ B rdbt0
where
rbt bonding strength of sandwich structure panel;
d0, B, Vj, j = x, y, xy definition see above.
The local buckling analysis described above is compiled in the calculation
software BUCKLSCP.
almost always produce local highly stressed areas. Stress concentration in com-
posites is not only more severe but also more complex than that in metals. Stress
concentration in metals depends only on geometry; however, composites are
affected by the layering pattern as well as geometric parameters. Well-established
joining technologies for metallic structures are not directly applicable to
composites.
Stress concentration in mechanically fastened joints is particularly severe
because the load transfer between the elements of the joint has to take place over a
fraction of the available area.
Composite joint strength is closely related to the layering pattern, load direction,
and environment. There are more failure modes of composite joints, and moreover,
strength prediction is more difficult. These complicating factors require careful
consideration.
This section deals with the joining of advanced fiber composites, mainly
focusing on mechanically fastened and adhesively bonded joints.
There are two methods of advanced composite joining: adhesive bonded and
mechanical fastening [2, 13, 17, 18].
(3) The precision of the fastener fit with the hole is important. A poor fit will
increase the shear deformation of the joint, resulting in shear failure of the
bond-line, and induce shear failure of the fasters and bearing failure of holes.
Hence, there may be no net benefit to the use of fasteners and bonding.
The selection of the joining methods should seek to take advantage of the respective
features of joint types. In general, some basic principles should be followed:
(1) Bonded joints are generally suitable for thin structures with low running loads
(load per unit width, i.e., stress element thickness) or structures carrying
shear load. The main advantages of bonded joints are their lightweight nature
and high joint efficiency. Thus, bonded construction tends to be more prevalent
in small light aircraft and secondary aircraft structures. Well-designed, bonded
joints can also transmit large loads;
(2) Mechanically fastened joints are mainly used in structures where concentrated
loads occur or an emphasize on high reliability is required. Bolted joints can
transfer greater loads than riveted joints. Thus, bolted joints are mainly used in
primary aircraft structural components. The main disadvantage of mechanically
fastened joints is the decrease in the strength of the basic laminate owing to the
fastener holes;
(3) Combined joints are generally suitable for jointing places requiring greater
margins and for medium thickness laminates.
Bonded joints have advantages in terms of their lightweight and high joint effi-
ciency; thus, their use in aircraft structural components has grown. For example, the
spar of the B-737 horizontal stabilizers; the root-stepped joins of the F-14
all-movable horizontal stabilizers; the joints of wing panel-to-root rib of the F-15
aircraft; Joints of fuselage panels to frame and joints of skin-to-skin for the Lear Fan
2100 all composite plane; joints of the skin of the pelvic fin of the clapboard for the
Y7-200B; the skin-stringer joint of the Y7-FC vertical stabilizer; joints of the
p-stringer to panels, and the p-stringer to web for the DC-10 vertical stabilizer wall.
Bonded step lap joints are used in the attachments for the F-14 and F-15 horizontal
stabilizers as well as the F-18 wing root fitting, and the majority of the airframe
components in the Lear Fan and the Beech Starship [2, 17–24].
4 Composite Structure Design and Analysis 465
The main factors affecting adhesive joints strength include: material of the adher-
ends, stiffness ratio and thermal expansion coefficients of the adherends, joint
configuration and geometry, fiber orientation of the bond-line, temperature and
moisture, adhesive, and manufacturing procedure.
(1) Effects of unbalanced adherend stiffness: All types of joint geometry are
adversely affected by unequal adherend stiffness, where the stiffness is defined
as the axial or in-plane shear modulus multiplied by the adherend thickness.
As an example, for single-lap joints, if the stiffness of the adherends is bal-
anced, the bending moments at two ends of the joint will be the same and the
deformation of the adherends will be equal. If the stiffness of the adherends is
unequal, the bending moment at two ends of the joint will be different and a
higher deformation will generally occur at the loaded end of the more flexible
adherend.
Where possible, the stiffness of adherends should be kept approximately equal.
For example, for step lap and scarf joints between quasi-isotropic carbon/epoxy
and titanium (Young’s moduli: 55 and 110 GPa, respectively) ideally, the ratio
of the maximum thickness (the thickness just beyond the end of the joint) of the
composite adherend to that of the titanium should be 110/55 = 2.0.
466 Z. Shen et al.
Any bond defects will result in load redistribution along the entire bond-line
and stress from discontinuity of the bond-line will increase. When the defect
size of debonding and cracks is small compared with the length of the
bond-line, any increase in stress will not be obvious. The stress will increase
markedly as the defect size increases. Thus, it is necessary to establish stan-
dards for bond quality.
4.8.2.3 Adhesives
(1) General requirements of adhesives
Adhesives should have the following features:
(1) Compatibility with the adherends and high bonded strength, such that
bond-interface failures will not occur;
(2) The curing temperature should be as low as possible;
(3) The thermal expansion coefficient of the adhesive should be nearly identical to
that of the adherends;
(4) Temperature effects should be minimal;
(5) Good mechanical properties;
(6) Simple processing;
(7) The durability of the bond should be greater than the anticipated life of the
structure.
(2) Types of adhesive and their selection
Adhesives can be broadly classified into two major groups on the basis of their
stress–strain curve, i.e., ductile and brittle adhesives (Fig. 4.67). The limit of shear
strain of a ductile adhesive is greater than 0.05, whereas that of brittle adhesive is
typically far less than 0.05.
As shown in Fig. 4.67, the shear strength of a brittle adhesive is higher than that
of a ductile adhesive. However, peel stresses can be eliminated from consideration
by approaches such as adherend tapering. The static shear strength of the bonded
joint does not depend only on a single parameter and is determined by the strain
energy to failure of the adhesive under a shear load (i.e., the area under the curve).
Therefore, joints based on ductile adhesives have greater strength. From the
viewpoint of fatigue performance, a brittle adhesive will rupture near the inflexion
and its fatigue life is lower. The ultimate strain of a ductile adhesive is also greater.
Ductility in aerospace adhesives is beneficial in reducing stress peaks in the
adhesive, i.e., lowering the stress concentration. If higher fatigue stresses can be
withstood, the fatigue life will be longer. When the environment temperature does
not exceed 70 °C, ductile adhesives should be used as far as possible.
Near the engine or in ultrasonic airplanes, high operating temperatures neces-
sitate that the brittle adhesives are used despite the loss of strength.
It is necessary to consider the effects of temperature. If the temperature remains
below the glass transition temperature of the adhesive, the bond strength will not be
sensitive to temperature effects. However, the strength will be reduced at low
temperature.
Materials commonly used in structural adhesive bonding of composite structures
are thermosetting resins, which can be subdivided into four basic chemical classes:
epoxy, polyimide, phenolic, and silicone.
(1) Epoxy: The advantages of epoxy resins include its high strength and modulus,
low levels of volatiles, excellent adhesion, low shrinkage, low moisture
absorption, good adhesion, good chemical resistance, and ease of processing.
Therefore, epoxy resins are the most widely used structural adhesives. Forms
are packed with resin and curing agents, which are mixed and cured with heat.
The major disadvantages of epoxy resins include brittleness, generic hardness,
low thermal strength, and poor wear characteristics. The curing is usually
accomplished by the application of heat under pressure. For example, a cure
will typically be performed at 145 °C and 0.7 MPa and be complete within
20 min. Some cures will also be completed at room temperature.
(2) Epoxy–Phenolic: This class of adhesives are a modified epoxy, which can be
completed within 60 min at 250–350 °C. Its advantages include high strength
and good performance at low temperatures; its major disadvantages are the
need to heat during curing, porosity of the bond, and poor electrical
performance.
(3) Polyimides: This class of adhesives requires high temperature curing, usually
between 250−400 °C. A post-cure is also required to attain maximum strength.
The highest operating temperatures of these adhesives are in the range of 250–
400 °C. Advantages of this class of adhesives include their resistance to tem-
perature, moisture, fire, and corrosion as well as their low coefficient of thermal
expansion. Disadvantages of polyimides include their high cost, porosity, and
corrosiveness.
4 Composite Structure Design and Analysis 469
(4) Phenolic: Mixed resin adhesives are usually composed of a phenolic resin
mixed with another resin. The advantages of such mixtures include high
thermal strength, acid resistance, low cost, and good electric performance. Their
major disadvantages are the need for high curing temperatures, high shrinkage,
and corrosiveness. Common used resins include:
① Phenolic polyamide: Shear strength can be as high as 36 MPa and
maintains excellent strength at high temperature.
② Phenolic ethylene: Shear strength can be as high as 30 MPa and can
operate at very low temperatures. Performance is rapidly degraded
above 100 °C.
(5) Silicone has good resistance to heat, cold, radiation, and good isolation;
however, its strength is low. Therefore, joints requiring high stability and the
high mechanical strength may be achieved with the use of this resin in com-
bination with others. Epoxy–silicone can be used continuously at temperatures
as high as 340 °C and discontinuously at temperatures up to 510 °C.
(3) Adhesives suitable for bonding different materials
Adhesives for bonding different materials may be selected as outlined in
Table 4.16. Blank entries for material/adhesive combinations in the tables indicate
that it may be difficult to achieve bonding.
Adhesives suitable for aeronautic structures are listed in Table 4.17.
(4) Measurements of the mechanical properties of the bond-line
Stress–strain characterization of adhesive films and their mechanical performances
form the basis of static strength design for adhesive bonded joints. Because the
bond-line is very thin, the interface will have some influence and the specimens
used for testing must have the same configuration as that of the actual part.
Measurement results show that the actual stress–strain curve (Fig. 4.68) is
complicated and may be difficult to apply directly for joint analysis. Equivalent
elastic–plastic and bilinear stress–strain curves are commonly used simplified
models. The elastic–plastic curve is particularly useful and the simplification allows
closed form analytical solutions to be obtained. The principles of this simplification
are that any adhesive is defined by two straight lines having the same strain energy
and failure stress and strain. The peak allowable shear stress should be multiplied
by a factor of 0.8 to account for both bonding defects and the differences between
laboratory and production fabrication. The peel strength and other data needed for
design can be measured from related test standards.
470
MPa 25 °C
33
20
30
33
30
28 24.5
80 °C
21(70 °C)
12
21 (70 °C)
100 °C
18
15
150 °C
10 (130 °C)
20
15 13.3
175 °C
18
90° peel 55 °C
4.5②
4.5②
strength/(kN/m) 25 °C
6.0②
5.0
6.0②
5.9
7.0 7.5
150 ° C
3.9 4.0
Manufacturers BIAM BIAM HIP HIP BIAM HIP HIP
Overseas similar Metland Metland Redux 319A
trademarks 1113.06 1113.06 Redux 119
Metlbond Metlbond 6726
6726
Overseas corresponding DHS172−292 DHS172−292 DHS174−292
materials standard DHS186−211 DHS186−211 DHS186−231
471
Note ①. All data of J−116B adhesives are B−basis; ②. Bell peel strength
472 Z. Shen et al.
The overall purpose of these principles is to ensure that the strength of the
bonded layer is higher than or close to that of the adherends. Therefore, it is
necessary to adopt measures to ensure that the configuration and geometric
parameter satisfy these requirements.
(2) Failure modes of bonded joints
Basic failure modes for adhesively bonded composite joints are as follows (shown
in Fig. 4.69):
(1) Tension (or tension-bending) failure of adherend;
(2) Shear failure of glue-line;
(3) Peel failure of bond-line and adherends.
Alongside these three basic failure modes, combined modes may also occur. The
failure modes of bonded joints will depend on the joint configuration, geometric
parameters, fiber direction near the glue-line, and loading properties. The adherend
thickness is the most important geometric parameter, as outlined for the following
cases:
(1) Tension (or tension-bending) failure of adherend will occur when adherends are
very thin, and joint strength is sufficient;
474 Z. Shen et al.
(2) Shear failure of the glue-line will occur when the adherends are thick and the
eccentric moment is small;
(3) Peel failure will occur under eccentric moments when the adherend thickness
reaches a certain value and the bond length is not long. The interlaminar tension
strength of CFRP is very low; thus, composites are prone to interlaminar ten-
sion failure. Peel failure will reduce the load capability greatly and should be
avoided.
(3) Selection of basic joint configuration for bonded joints
Figure 4.70 shows some basic joint configurations for panel components of aircraft.
The selection of a joint configuration is key for bonded joint design. Joints must
be designed to transfer their maximum load in the shear direction with smaller loads
in other directions. This will avoid the occurrence of large peel stress. Figure 4.71
illustrates the strengths of basic joint classifications as a function of the adherent
thickness. Each curve shown represents the best strength that can possibly be
obtained for each joint type.
(1) Single-lap joints may be used when adherends are thin ( 1.8 mm). Note that
additional bending moments, caused by eccentricity of the load path, will
result in very high peel stress at both ends of the bonded joint, which will
reduce the joint strength. Therefore, it is necessary to increase the
overlap-to-thickness ratio. Bending moments may be alleviated through the
use of a high ratio L/t = 50–100. When adherends feature an imbalance of
stiffness, eccentricity effects will be greater. The use of single-lap joints should
be avoided. However, in a single-lap joint supported against bending,
eccentricity effects may be alleviated and deformations restricted. Such joints
4 Composite Structure Design and Analysis 475
Fig. 4.73 Stress and strain distribution of balanced stiffness double-lap joint
between adherends is required for thin adhesive layers such that a very thin
glue-line may not be possible.
The design of simple, bonded splices of uniform thickness for near
quasi-isotropic carbon/epoxy is simple. Use a 30t-overlap for double shear, 80t-
overlap for single-lap joints, and a 1-in-50 slope for scarf joints. Overlapping ends
should taper to 0.51 mm with a slope of 1/10.
Analysis and test results indicate that the shear stress distribution is not uniform
throughout the bonded area under an applied load. Most of the load is transferred
through two end zones which form a low stress elastic trough (Fig. 4.73).
Because of the presence of these elastic troughs, the load carrying capacity of the
bonded joint increases gradually in the beginning. However, the width and depth of
the elastic trough only increases continuously when the length of the overlap attains
a certain value. Increases in overlap length above this value do not add to the joint’s
load carrying capacity. From the viewpoint of static strength, there is no need to
increase the overlap length; however, service life and durability should also be
considered and longer overlaps are often used. For very short overlap and transfer
of large loads, the minimum shear stress and strain in the middle of the overlap area
are nearly equal to that at both ends. Thus, the entire bond is in a plastic state. When
the load is removed, the adhesive in the middle cannot recover, and the joint will
fail soon. Analysis results demonstrate that at a minimum stress equal to 10% of the
maximum stress the glue-line can recover its original state. For double-lap joints, an
elastic trough width of 6/k is sufficient to ensure a minimum adhesive shear stress
distribution, which is no greater than 10% of the maximum stress.
(5) Fiber orientation of the bond surface
The surface fiber direction of the laminate should be in the primary load direction or
45° to the load direction, but not perpendicular to the load direction, to prevent
adherend premature interlaminar tension (peel) failure.
(6) Surface preparation of adherends
The bonding of adhesive is a complicated activation process between the adherends
and adhesives. It is important to prepare a quality adherend surface for good quality
bonds in terms of static strength and durability. The bond should meet prescriptive
technical specifications. Strict quality control and inspection should be performed in
the bond processing. Nondestructive inspection should be performed for all
important parts. Surface preparation deficiencies are particularly troublesome
because there are currently no nondestructive evaluation techniques for detecting
low interfacial strength between the bond and adherends.
For bonds between carbon–epoxy composites, solvents may be used to clean the
surface together with mechanical abrasion of the surface. For bonds between
composites and metal, in addition to the surface preparation the metal will require a
surface treatment. Corrosion barriers (such as fiberglass and sealants) are placed at
the interfaces between the composites and aluminum or steel to prevent galvanic
corrosion.
4 Composite Structure Design and Analysis 477
Thin section joints can only transfer small loads; however, it is possible to transfer
larger loads through thick section joints. Failure will occur preferentially in the
adhesive for thick adherends in a simple joint configuration. Thick adherends
cannot perform effectively and joint efficiency will be low. To ensure the glue-line
is not a weak link in the joint and to make full use of the load-bearing capacity of
the adherends and avoid premature failure, the bond surfaces should be increased
and peel stress reduced. Complex stepped and tapped joints are typically used.
(1) Selection of stepped and scarf joints: Thick adherends under a large load are
suitable for stepped or tapped joints. The use of stepped or scarf joints is
effective for reducing peel stress. The advantages of stepped joints over scarf
joints are their ease of fit and high strength achievable by adjusting structural
parameters. Therefore, higher joint efficiency may be attained.
Composite-to-titanium stepped joints are used extensively throughout the
aerospace industry for high load transfer.
Stepped and scarf joints are appropriate for highly loaded thick plate bonded
joints. The use of scarf and stepped joints is effective for reducing peel stress.
Unlike scarf lap joints, stepped joints have simple processing and can achieve
high strength by adjusting structural parameters. Stepped-lap joints are com-
monly used for joining cover panels and titanium structures, see Fig. 4.74.
(2) Strain-level requirements: In the design of strain levels for thick adherend
structures, values should be properly lowered considering the need for future
repair. It is impractical to repair thick structures by bonding because of the taper
ratio requirement, i.e., 1:50. When there is no need for repairs such as one-shot
Detail design of composite bonded structures should not only consider the static
strength of the bonded structure but also the durability, bonding technology, and
cost of the bonded structure. In addition to the aforementioned basic principles, the
following issues should be noted for detailed design of composite bonded
structures.
(1) Selection of bonded joint configuration: The configuration of a bonded joint is
a critical design aspect. The load-bearing capabilities of bonded joints work
best in the shear direction and have poor resistance to peeling. The maximum
load should be transferred in the shear direction and minimum loading should
be induced in other directions. The use of stepped-lap or scarf-lap joints is
effective for reducing peel stress and ensuring the joint strength is not lower
than that beyond the joint.
4 Composite Structure Design and Analysis 479
(2) Procedures for reducing stress concentration and peel stress for adhesive
joints: Whenever possible, the peel stress in the structure should be reduced by
induction of eccentricities in the load path and asymmetry. For example, the
use of a symmetric double-lap joint increases the bending stiffness of the outer
adherend, and tapering of the edges of the overlap in single joints. Three
procedures for decreasing stress concentration are illustrated in Fig. 4.75.
A high stress concentration at both ends of the joint occurs owing to eccen-
tricity in the load path for single joints. Thus, the peel stress can result in the
premature failure of the adherends. The load-bearing capacity will be
improved by supporting a joint with a transverse fully stiffened restraint
(Fig. 4.76).
(3) Skin (web) strengthening procedure: The use of T-type components rather
than angled components is recommended. Peeling at the corner of the angle
reinforcement can occur easily if there is a tension force. Premature peeling
may be prevented with the use of filler at the corner (Fig. 4.77). When a T-type
element is used, peeling can be prevented if there are tension forces. Peeling
will be improved by edge tapering of the profile element and balancing the
stiffness between the profile element and web. The stress distribution of T
strengthened elements is shown in Fig. 4.78.
(4) Thermal stress of bonded structures: In the bonding of carbon/epoxy,
boron/epoxy composites to metals, such as titanium and steel, Thermal
residual stress arises due to differences in thermal expansion coefficient of the
materials. In particular, the thermal stress is proportional to the difference
between the operating and cured temperature. Thermal stress can be reduced
through the use of laminate layering design.
(5) Avoidance of galvanic corrosion: In the bonding of composites to metal,
galvanic corrosion may occur owing to differences between the electrode
potentials of materials. Surface treatment of metal elements should be per-
formed. Whenever possible, direct bonding of carbon composites to aluminum
should be avoided and an isolating layer should be placed between the
materials.
Carbon fibers must be isolated from aluminum or steel through the use of an
adhesive layer and/or a thin glass-fiber ply at such interfaces. The galvanic
interaction between carbon and aluminum or steel will cause corrosion of the
metal.
(6) Prevention of moisture entering adhesive layer: Unlike metal adherends,
composite adherends are subject to the effects of moisture diffusion. As a
result, moisture is more likely to affect the whole component rather than be
confined near the exposed edges of the joint in the case of metal adherends.
The response of adhesives to moisture is an important issue for composite
joints.
(7) Tooling design and manufacture: The quality of bonded joints is influenced
greatly by tooling. Therefore, careful attention is required for the design and
processing of bonded joints to minimize thermal deformation and residual
stress. Tools should be applied under uniform pressure to the adherends.
(8) Quality control: Adhesive quality should be controlled based on allowable
values for defects prescribed for different positions.
4 Composite Structure Design and Analysis 481
Stress analyses of adhesive joints range from very simplistic ‘P over A’ formula-
tions in which only the average shear stress in the bond layer are considered, to
extremely elegant elasticity approaches that consider fine details—for example,
calculation of stress singularities by applications of fracture mechanics concepts.
A compromise between these two extremes is desirable, because the adequacy of
structural joints does not usually depend on knowledge of their details at the
micromechanics level, but rather at the scale of the bond thickness. Practical
considerations require bonded joints to incorporate adherends, which are thin rel-
ative to their dimensions in the load direction; hence, the stress variation through
the thickness of the adherend and the adhesive layer tend to be moderate. Such
variations do tend to have a great effect on polymer matrix composite adherends
because of their relative softness with respect to transverse shear and thickness
482 Z. Shen et al.
The following points should be noted for mechanical joints [2, 13, 17]:
(1) Owing to the brittle nature of composite materials, multiple fastener joint load
distributions are non-uniform. The stress and strain of basic laminates will be
lower when joints fail;
(2) The bolted joint strength of laminates with a certain content of 0°-plies is less
than the unnotched laminate strength;
(3) The load-carrying capability of joints does not show a directly proportional
increase with the end distance;
(4) Bolted joints should be designed to carry a load such that the bolt is under a
shear force rather than tension. Bolt bending in composites is more common
than that in metals.
4 Composite Structure Design and Analysis 483
There are many more factors that affect the mechanical joint strength of composites
than those affecting metals. It is important to understand and consider all factors in
design.
These factors can be classified into the following five types:
(1) Material type and form: unidirectional tape or woven fabric fibers, resin type,
fiber orientation of, fiber volume fraction, and laminate pattern;
(2) Processing methods: prepreg, RFI, RTM, curing, and consolidation processes
(vacuum bag molding and oven and autoclave curing);
(3) Configuration: joint types (single or double lap), geometry (pitch, space, edge
distance, side-end distance, thickness, hole diameter and tolerance, hole pat-
terns, and washer size);
(4) Fastener types (hexagonal head bolt, big foot bolt, blind fastener, protruding
and countersunk head fastener), clamp-up force;
(5) Load: static, dynamic, fatigue load, load direction, loading rate;
(6) Environment: temperature and humidity.
(1) Laminate pattern
Laminates used in aerospace structures are generally composed of layers in the
0°, ±45° and 90° directions with respect to the axes of the laminate. The percentage
of ±45° plies has an important effect on laminate bearing strength. Shear-out or
cleavage failure can occur more readily when the ±45°-ply content is less than that
of the 0° plies. Unlike metals, shear-out failure can only be prevented by increasing
the end distance of holes. It is more important that a proper percentage of ±45°
plies is maintained. Bearing strength increases with the percentage of ±45° plies.
The recommended layering ranges to achieve maximum strength in joint areas
are ±45 plies
40%, 0° plies
30%, 90° plies in the range 10–25%, with
variations of 5% allowed. Bearing strength will decrease as the percentage of
the ±45° layers is increased further.
Characteristics of ±45° layer content
50% are as follows:
(1) Joint strength is less sensitive to load direction;
(2) Initial failure strength may occur earlier;
(3) Shear load-bearing ability is stronger and tension load-bearing ability is lower.
Particular care should be given to the tension in multi-row fastener joint design.
(2) Ply stacking sequence
The stacking sequence is a special parameter effecting the mechanical nature of
composites. Laminates of the same ply numbers and proportions can have various
stacking sequences, which can change interlaminar stresses, and the mechanical
natures of laminates may be affected. Whenever possible, maintain a well-dispersed
stacking sequence and avoid grouping similar plies.
484 Z. Shen et al.
te ¼ t for t d;
te ¼ d for t [ d:
(a) (b)
(1) The percentage of ±45°, 0°, and 90° plies should not be less than 40%, 30%,
and 10%, respectively. This is particularly important for mechanical joint
design;
(2) Extremely thin laminates should be reinforced locally at the attachment area to
provide greater thickness. This reinforcement will avoid the reduced bearing
allowables that result from a D/t ratio greater than four. The general rules
D/t
1 should be followed to avoid failure of the fastener;
(3) In areas of load induction there should be equal numbers of +45° and −45°
plies on each side of the mid-plane;
(4) Butt-splined fibers should be avoided in join areas.
(5) Geometry requirements
To prevent low strength failure and ensure high strength of mechanical joints,
geometric parameters of jointed plates should be selected according to Table 4.18.
Definitions of the geometric parameters are shown in Fig. 4.81.
In addition, the geometric size of joints should consider future repair demands.
The next largest size fastener should be useable after the repair.
(6) Fastener requirements
To prevent galvanic corrosion, fasteners made from titanium, titanium alloy,
stainless steel, and Monel should be used because the electrode potentials of these
alloys are close to those of the composites.
(1) Principles for selecting fastener diameter: General guidance for selecting fas-
tener diameters are as follows:
① Sufficient bearing strength of jointed component should be ensured.
where
D —fastener diameter, mm;
t —laminate thickness, mm;
rbr —allowable laminate bearing strength, MPa;
[sb] —allowable fastener shear strength, MPa.
Proper tightening torque can increase bolt-joint strength. The torque moment
should be selected based on the relational standard for various material, diameter,
and the bolt type. If there are no special requirements, the tightened torque may be
selected according to Table 4.19.
(7) Requirements of galvanic corrosion resistance for mechanical joints
The three conditions which lead to galvanic corrosion should be excluded in the
design: potential differences between materials, presence of an electrolyte, and the
electric connections. The following measures of corrosion prevention should be
used:
① Material matching can prevent galvanic corrosion. Metals that have
electrode potentials that match those of carbon/epoxy composites
include: titanium alloy and stainless steel.
② Prevention of electrolyte accumulation should be considered in the
design; sealing of joints should be performed to prevent infiltration of
electrolyte and avoiding corrosion battery formation.
③ For materials unsuitable for direct contact, an insulating layer of
glass/epoxy or aramid/epoxy should be used. At important joint sites
which may be predisposed to corrosion, full sealing of the joints should
be used to prevent corrosion.
④ Joints can be installed wet with sealant, in addition to insulation. In
riveted joints, it is important to wet set with sealant, to prevent galvanic
corrosion but also compensate for any manufacturing damage.
⑤ Fiber laminates in direct contact with aluminum and aluminum and
cadmium plated steel components should be avoided to prevent galvanic
corrosion. Otherwise, an insulating layer should be added. Carbon fiber
4 Composite Structure Design and Analysis 491
point of reducing weight and cost. However, these must be wet set and strict
measures taken to prevent galvanic corrosion prevention.
(5) Avoid buck rivets in composite structures. Squeeze rivets can be used if
washer is installed on the tail side.
(3) Measures to improve pull-out strength
For outer surface of structures, such as the rudder, a hole cap and countersink can
be used to strengthen the structure with titanium alloy or stainless steel to improve
the pull-out strength.
(4) Reliable measures for galvanic corrosion prevention
(5) Riveting processing requirements
(1) Riveting should follow technology specifications. Strict quality control and
inspection should be conducted during hole drilling, countersinking, and
riveting. Nondestructive evaluation should be conducted for important parts;
(2) Damage to the exit site of the drill should be prevented by coating the com-
posite with a layer of film adhesive, glass-cloth, or a pad plate;
(3) When composites come into direct contact with metallic components, under
structure permissive conditions, snap the head of the rivet at the metallic
surface whenever possible. If the snapped head of the rivet is on a composite
surface, a pure titanium, titanium alloy, or stainless steel washers must be
placed on the snap head;
(4) Whenever possible, squeeze rivets should be used for parts requiring common
solid rivets. Bull rivets may be considered where squeeze riveting cannot
conducted. Strong power rivets should be avoided.
Mechanically fastened joints are the main joint type used in primary composite
structures. To meet structural integrality requirements, in addition to meeting
strength and stiffness requirements, fatigue, damage tolerance, and functional
requirements must also be satisfied. Stress concentration in mechanical joints can
create fatigue weak points in the primary composite structure. Fatigue strength is
determined mainly by testing now, because methods for pre-estimating the life time
of composite joints are not mature, and are complicated by environmental
conditions.
Three fatigue failure criteria should be considered in the rational design of
mechanical joints under wet-heat conditions and different load spectrums.
Tension, shear-out, and bearing failure of fasteners loaded hole; permanent
elongation deformation of fastener holes exceeding allowables; residual strength of
joints is lower than the design requirements. The joint life will fail when any one of
aforementioned items occurs. Generally, permanent elongate deformation of load-
ing holes is the first limiting value.
4 Composite Structure Design and Analysis 493
One major difference in the mechanical behavior of composite materials and metals
is that composites are brittle and anisotropic; while metals are plastic. Metal has the
capability to redistribute load, thus allowing each of the fastener holes of a multirow
joint to uniformly carry the load distribution. However, for brittle composite
materials this is not the case [26–31].
Composite (fiber-dominated) laminates generally show linear behavior up until
failure. The material will not yield locally and redistribute stress. Effective joint
design should adopt measures to reduce the bolt bearing stress in the most critically
loaded locations. Even if at ultimate load non-uniformity of the fastener load dis-
tribution shows little improvement in comparison with the initial load for steel or
titanium fasteners. Effective joint design requires that the greatest load-bearing
fastener row should be reduced.
The strength of multirow bolted joints in composite structures is governed by
associated bearing–bypass load interactions under tensile or compressive loads. The
key to obtaining high operating strain in bolted joints in fibrous composite lami-
nates is to restrict the bolt bearing stress in the most critically loaded locations. By
tailoring the joint geometry, a bolt load distribution can be generated which
maintains low bearing high bypass conditions in the first or outermost row of
fasteners. With efficient joint design, cross-section strain in basic skin laminates can
reach 0.005 in room-temperature tests.
The laminate fiber pattern is a design variable and optimizing the joint for
maximum strain does not guarantee the highest strength or the most weight-efficient
design. The principle design parameter governing the design of composite joints is
the amount of load that must be transferred rather than the operating strain level of
the adjacent structure.
494 Z. Shen et al.
fiber pattern for optimal layer compositions. For carbon–epoxy laminates, the
optimum w/d is likely to be in the range 4–5 for multirow joints.
(7) Adequate consideration of the bolt diameter-to-laminate thickness ratio (or
more appropriately, bolt bending stiffness-to-laminate thickness ratio) is
warranted in joint design to assure that fasteners are the weak link.
Fastener bending elastic deformation may decrease the clamp force and
allowable bearing stress and should thus be avoided. Therefore, selection of
fastener sizes should not be based only on the rated shear strength of the
fasteners but should also consider the fastener stiffness.
(8) Interference fit systems with a sleeve of fasteners having the same outer
diameter as the sleeve, generally do not feature increased strength (strength
may actually decrease slightly). This is because any potential benefits are
negated by recurrent bolt bending failures.
(9) Materials should be selected to take advantage of their strengths while
avoiding their weakness. Metals should be used in parts for which composite
materials are unsuitable.
Metal materials are selected for splice plate members for several reasons.
If protruding head fasteners are used in the subcomponent tension joints of
tapered composite splice plates, tapered members require either spot-facing
of the splice plate surface or the use of tapered washers under the fastener
heads and nuts. These features may cause premature failure owing to the high
peel stress and interlaminar forces. The use of tapered washers also increases
the cost and complexity of the assembly procedure. Thus, metallic splice
plates with spot-facing on tapered surfaces are used to accommodate the
fastener seating. The use of metallic splice plates is the simplest and most
cost-effective way of avoiding these potential failure modes.
Composite materials are not well-suited to applications where high out of
plane forces are present. The T-splice members are likely to encounter such
forces, and the magnitude of the forces is very difficult to predict analytically
or measure experimentally. The fabrication of the corner fittings based on
composite materials would be impractical for similar reasons and
cost-prohibitive compared with the use of aluminum parts.
The splice plates may be slightly heavier, owing to the use of metals;
however, any small extra weight in the splices (or fasteners) is compensated
by maximizing the efficiency of the large heavy skins. For a large airplane,
the weight of the splicing elements as a percentage of the total wing weight is
small, and splice efficiencies should be evaluated solely on the basis of the
minimum splice and fastener weight.
(10) Joint strength is typically greater under compression than under tension
loading. An example of the application of these principles is presented in
Fig. 4.82. An optimum splice structure is represented, including a cover of
uniform thickness, tapered splice plates and varying diameter fasteners. The
bolt diameter of the inner most row near the cover butt is largest, S/D = 3.
There are no bypass loads on the skin. The combination of maximum bearing
496 Z. Shen et al.
and bypass loads act on the splice plate of the innermost skin such that the
splice plate thickness must be properly increased. Research results indicated
that the optimum splice plate thickness is 1.5 times the basic plate thickness.
In the example, the basic plate thickness is 12.7 mm, and the total plate
thickness of the taped splice is 19.1 mm, including the thickness of both the
top and down taped splice (9.5 mm each). The diameter of the middle two
row bolts has an intermediate S/D value of 4. The diameter of the outermost
row of bolts is smallest with S/D = 5. A low thickness of the splice plate
outer end may result in shear failure of laminates under a large load, which
should be avoided.
Static analyses of mechanical joints generally include the following three aspects
[2, 17, 25–31]:
(1) Exterior forces acting on the mechanical joint are determined from overall
structural analysis of the whole joint.
(2) These forces are then used to determine individual fastener loads and bypass
forces acting at each fastener hole of the joint.
(3) Joint strength can be assessed by applying two methods: one is the semiem-
pirical failure envelope method; another is to use material failure criteria and
characteristic curves.
4 Composite Structure Design and Analysis 497
design purposes, the following relation of bolt flexibly can be used, which is
simple and satisfies engineering precision requirements.
where
Ks —shear stiffness of the fastener;
As —shear area of the fastener;
G —shear modulus of the fastener;
L —effective length of the fastener.
For single shear, the effective length l can be assumed to be one-fourth the
combined thicknesses of the attached sheets. The effective length in double shear
can be approximated as half the single shear value. Equation 4.8.2 is used for
fasteners where only shear is accounted for; fastener bending and rigid body
rotation (in a single-lap joint) are not considered. The fastener load distribution
derived from these relations will be slight conservative.
(2) Nonlinear analysis: Load-deflection (P–d) curves from single fastener joint tests
can be modeled as bilinear curves, as shown in Fig. 4.83. The nonlinear
strength analysis should permit some bolts to fail while the structure should still
be able to carry loads. Nonlinear analysis can provide more exact load-sharing
analysis and ultimate strength predictions.
Stress analysis methods of single fastener joints are described in detail in Ref. [1].
After the load distribution is determined (i.e., the bearing and bypass loads of
fastener holes) analytical methods for single fastener joints can be used to calculate
the detailed stress and strain near the fastener hole. Finally, joint strength and failure
modes can be assessed by applying material failure criteria or characteristic curves.
Theoretical analytical methods of single fastener joints mainly depend on ana-
lytical and FEMs. In finite element analysis a fine mesh must be used in regions of
high stress gradients, such as around the cutouts and at ply and stiffener drop-offs.
Joint analysis should include the effects of shimming to the limits permitted by
drawings. The effects of shimming may reduce joint strength. The effects of per-
missible manufacturing parameters should be considered, for example, hole per-
pendicularity (±10°), shimming, and loose holes.
Semiempirical Methods
where
rb —unnotched laminate tensile strength;
rbr —loaded hole bearing stress;
rnet —laminate net-tension stress caused by bypass loads;
Kbc —composite bearing stress concentration factor, with respect to bearing
stress;
500 Z. Shen et al.
Otherwise, joint-bearing failure will occur when the bearing stress achieves the
bearing strength:
The information presented here is not only applicable to single fastener joints, but
also useful for determination of multirow joint strength. All joint strength data are
developed from tensile test results, and the results will be conservative for use in
compression loads [2, 17].
(1) Allowable bearing stress
4 Composite Structure Design and Analysis 503
To ensure structure integrity the loading should generally not be greater than the
initial bearing failure stress. Therefore, selection of an allowable bearing stress
strongly depends on the failure criteria. The initial bearing stress is very different to
the definition of failure. Failure criteria can be classified in different ways: One
approach is to base failure on stress, which guarantees that structures have sufficient
strength; another is based on deflection of the loaded hole, which guarantees that
structure have sufficient stiffness. One frequently used approach is based on the
degree of hole deformation. However, the failure deflection limits of loaded holes
selected by various countries and departments are very different, ranging from 0.5
to 6%.
The following are recommended criteria for determining the initial bearing
failure stress of a loaded hole: the lowest value between the first slope inflexion
point and bearing deformation of 4% in the load deformation curve. Experience
indicates that the minimum initial bearing failure stress can be considered to be half
of the ultimate bearing strength rbru.
The ratio of the initial bearing failure stress to the ultimate bearing strength
depends on the material system and laminate pattern. Generally, the ratio decreases
as the proportion of ±45° plies is increased. For laminate frequently used in joint
areas, the ratio is in the range 0.55–0.66. Selection of allowable bearing stress
should also consider joint importance, structural characteristics, load type, dura-
bility and service life and environmental effects.
Allowable bearing stress can be determined from the following:
where
Cw correlation factor for width;
Ce correlation factor for end distance;
Cp correlation factor for load direction;
Cd correlation factor for hole diameter;
Cs correlation factor for single shear;
Cen correlation factor for environment;
K factor considering initial failure, durability, aging, and technological quality.
The value of K is typically in the range 0.50–0.66;
rbru bearing strength, MPa.
The bearing strengths of several composite laminates are illustrated in Fig. 4.86.
Various correlation factors are shown in Fig. 4.86 for W/D
6, e/D
4,
D/t = 1.0–2.0, D = 5 mm, double shear, torque 4 Nm, at room temperature, in dry
conditions.
For laminates typically used in joint areas (i.e., 0°-plies = 25–
50%, ±45°-plies
40%, 90°-plies = 10–25%), the allowable bearing stress for
HT3/QY891 and HT3/4211 can be taken as 600 and 500 MPa, respectively.
Hence, the formula (4.78) is a concise, convenience, and effective model. The
effects of many parameters have been considered in various correlation factors, and
therefore, numerous procedures can be avoided. Traditionally, both the bearing
strength and tension strength as well as shear strength would require checking. This
method has been successfully used in joint design for many aircraft structures.
(2) Bearing strength
To fully develop the bearing capability, joint geometry selection requires that
bearing failure or combined failure modes depending on bearing failure are con-
sidered. Full load-bearing failure strengths are the foundation of joint design.
Failure modes are dependent not only on geometric parameters but also the fiber
pattern. Full bearing failures typically occur when W/D = 6 and e/D = 4 in the
laminate pattern range of joint areas. The bearing strengths given in this paragraph
are equal to the ultimate load divided by the bearing area Dt.
4 Composite Structure Design and Analysis 505
where
Pbr —fastener load, N;
D —hole diameter, mm;
te —plate effective thickness, defined as:
te = t, when t D,
te = D, when t > D;
[rbr] —allowable bearing value, MPa.
Note that tension and shear strength will be satisfied automatically without
checking because the effects of width and edge distance have been considered
in the allowable bearing strength value.
(2) Checking of fastener strength: The shear strengths of single fastener joints can
be checked as follows:
(2) Tensile strength checking: For strength checking, the tensile loading of a
multirow joint, rb, of the right side of Eq. (4.71) can be replaced by the
allowable tension stress of the laminate, i.e.,
where
[r] —allowable tension stress of laminate, MPa, [r] = Ext[e];
Ext —longitudinal tensile elasticity modulus of laminate, MPa;
[e] —allowable tensile strain of laminate.
Design allowable strains can be classified on A-basis and B-basis. The use of
either basis depends on the structure design criteria of the practical engineering
project. Generally, for components without a structure test or single path
transfer component, A-basis is used; B-basis is used for multi-path transfer or
fail-safe components. For carbon fiber resin matrices composites, allowable
tension strains are [eA] = 0.0082 for A-basis and [eB] = 0.0090 for B-basis.
Shear failure will not occur within the ply range recommend for joint areas
when the pitch is not less than 4D and the edge distance is not less than 3D.
(2) Compressive load
For strength checking of multirow joints under a compressive load, rbru on the
right side of Eq. (4.71) may be replaced by the allowable bearing stress of the
laminate, [rbr], i.e.,
4.9.1 Overview
Damage can be divided into two types according to its source: manufacturing
defects, which cover structural abnormalities caused by production, and operational
damage, which covers structural abnormalities caused in service [1, 2].
Manufacturing defects can usually be divided into two categories: First, lamination
and part curing processes may create defects such as voids, delamination,
debonding, inclusions, resin-rich or resin-poor areas, improperly cured resin,
deviation of fiber orientation (fiber bending), layering sequence errors, and gaps
between fibers. Second, defects may be produced in machining, packing and
delivery such as scratches, abrasion, improper hole drilling, and torque and impact
damage.
external bodies such as bullets, non-inclusive engine fragments and bird-strikes are
classed as high-energy impacts. These events are also known as high-speed impacts
and can produce penetrating damage with a certain amount of delamination.
Lightning can also break through the structure of the skin and produce deep
delamination and burning. This type of damage is visibly inspectable and can be
detected, allowing the part to be replaced or repaired. During production and
maintenance, low-energy impacts include events such as: tool dropping; impact
with maintenance facilities such as forklifts, trucks, and work platforms; damage by
personnel standing on structures; impacts caused by stones, screws, and tire frag-
ments during taking off or landing; impact of hail stones.
In fact, impact damage modes depend not only on external impact energy, but
also the laminate thickness. For thin skins or thin surface panels, impact damage
mainly results in fiber fracture, or penetration, resulting in decreased compressive
and shear strength. Furthermore, after such damage water may diffuse into the
sandwich core and causing durability issues. For medium thickness laminated
structures (less than 6 mm), impact damage may not be visible from the surface.
However, damage may be induced inside the laminate in the form of delamination
or matrix cracking. Such damage will greatly reduce the compression strength of
the component and presents damage tolerance safety issues.
Great attention has been paid to the effects of defects/damage on the strength of
composites. Since their initial use in aircraft primary structures in the 1970s, many
tests and investigations have been performed on the effects of damage on com-
posites. On the basis of test data derived from various composite material systems
(mainly carbon/epoxy, and carbon/BMI systems), and studies on the effects of
defects/damages on the static strength and fatigue strength of specimens under
different ambient conditions (room temperature/dry, hot/wet, cold/dry), the effects
of defects/damage on composite strength have been established as follows:
(1) Tensile loading: Many dangerous defects, such as cuts and slots, are inspectable
to some extent. The residual strength of laminates containing cuts will mainly
depend on the width, and is basically independent of the cut shape. Test results
for an open hole (typically 6.35 mm in diameter) can be used to consider the
strength reduction associated with an edge cut of similar size, when the
structure design allowable values are to be determined.
(2) Compression loading
4 Composite Structure Design and Analysis 513
d pffiffiffiffiffiffiffiffiffiffiffiffi
ry ðx; 0Þ 1 þ a 2x=W ¼0 ð4:84Þ
dx
Where
where
where
r0y ðx; 0Þ —the normal stress distribution of 0° plies on the notch cross
section without considering damage zone influence;
where
ry(x, 0) —stress distribution in Y direction of the minimum cross section
with a hole;
R —hole radius, and half length of the central crack;
a0 —characteristic length determined by testing.
2rb ð1 n2 Þ
r1
c ¼ ; ð4:88Þ
2 n22 þ ðKT1 3Þðn62 n82 Þ
n42
where n2 ¼ R þR a0 :
② Point stress criterion (PS): This criterion assumes that failure of
the laminate will occur if the stress ry at a point d0, a charac-
teristic distance, reaches the ultimate strength rb of a notch-free
laminate (Fig. 4.92), that is:
2rb
r1
c ¼ ð4:90Þ
2 þ n24 þ 3n44 ðKT1 3Þð5n64 7n84 Þ
where n4 ¼ R þR d0 :
③ Characteristic length a0 and d0: The characteristic length a0 and
d0 in average stress criterion and point stress criterion are
determined by testing. A number of specimens with different
hole sizes and crack lengths are used for tensile failure testing to
obtain a set of residual strength data ðr1 c ÞT These data are
substituted into the residual calculations by Eqs. (4.88)
and (4.90) based on the average stress criterion and point stress
criterion. The finite width correction and notch-free specimen
tensile strength r0, and a set of a0 and d0 values corresponding
to hole diameter, and crack length can be derived. Their average
values will be the characteristic lengths a0 and d0.
④ Finite width correction: The above-mentioned open hole lami-
nate or cracked laminate residual strength r1 c is the stress of a
laminate with an infinite width. Thus, corrections should be
performed for finite width laminates. Let rc be the residual
strength of a finite width laminate, such that:
r1
c ¼ grc ð4:91Þ
2 þ ð1 2R=WÞ3
g1 ¼ ð4:92Þ
3ð1 2R=WÞ
For laminates with an ellipse hole (long axis is 2a, short axis is
2b):
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
2
k2 1 2k 2a
g2 ¼ þ 1 þ ðk 2
1Þ M
ð1 kÞ2 ð1 kÞ2 W
"
2 #1=2 ð4:93Þ
2
k2 2a 2a
M 1 þ ðk 1Þ
2
M
1k W W
where
0vu
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
" #ffi 1
u 3ð1 2a=wÞ
M 2 ¼ @ t1 8 1 1A=2ð2a=wÞ2 ; k
2 þ ð1 2a=wÞ3
b
¼ :
a
the impact damage, such as shape, size, and distribution along the thick-
ness. This information can be derived from testing and inspection (such as
nondestructive CT-scans and X-ray methods), or from quantitative ana-
lytical estimations. The analytical estimation of composite laminate impact
damage includes two parts: a) analysis of the impact transient response of
laminates, b) the use of appropriate failure criteria to calculate the impact
zone, which will be mainly discussed in this section.
Impact damage of composite laminates includes matrix cracking, fiber
rupture, and delamination. In the following section, methods for calculating
impact damage size will be discussed based on delamination failure crite-
rion, which can be used in composite structural design.
① Delamination failure criterion: In terms of bending strain energy
density delamination failure criterion, If impact delamination of
composite laminates is dominated by matrix strength and
interlaminar strength, initial delamination can be derived from
the criterion:
2
2
2
rL2 sU sL31
eD ¼ f 1 þ f2 23
þ f3
1 ð4:97Þ
YL SU
23 Si
t Q11 QnU
f1 ¼ 11
ð4:98Þ
tU QL11 QU
22
1 GU GU
f2 ¼ 1 þ 23 cos 2
ðDhÞ þ 23
sin 2
ðDhÞ ð4:99Þ
2 GL23 GL31
15 GU GU
f3 ¼ 1 þ 31 cos 2
ðDhÞ þ 31
sin 2
ðDhÞ ð4:100Þ
16 GL31 GL23
tdmax ð1 h Þ
tdmax
a¼ 2 1 ð4:101Þ
h
Composite laminated structures can offer excellent fatigue performance. For common
fiber-dominatedmulti-directionlaminates(includingspecimenswithholes),thetensile–
tensilefatiguelifeis106 cyclesunderamaximumstressequalto80%oftheultimatetensile
strength. In the case of tensile–compression fatigue, the fatigue strength will be slightly
4 Composite Structure Design and Analysis 523
lower,withthefatiguestrengthequalto50%ofthecorrespondingstaticstrengthafter106
cycles. In particular, for specimens with impact damage, the fatigue strength will not be
lower than 60% of the corresponding static strength after 106 cycles. In thermoplastic
composites, values may reach up to 65% under the same conditions. Currently used
allowablesincompositestructuredesignmainlydependondamagetoleranceallowables.
Under such strain levels, composite structures can have infinite service life, which is
so-called static cover fatigue. Special attention should be paid to adhesive structures
becausefatiguefailuremayoccurifthedesignisperformedincorrectly.Currently,fatigue
failure is not a critical problem in design; however, no mature analytical methods for
durabilityarecurrentlyavailable.
low modulus ±45° plies (at a percentage of 70–80%), and contains a certain
amount of glass-fibers in some local zones, which can bear shear load and the
internal pressure of an oil tank for example. Laminates with a small ratio of 0°,
90° plies, such as a (10/80/10) ply ratio, can also be used to ensure local
strength and structural stability. Stiffeners mainly use 0° plies that are orientated
along the wing span direction and can be used to withstand tensile and com-
pression loading in wing panels. The skin and stiffeners are mechanically joined
or co-cured to form the wing panel. In some design programs (such as for body
panels), a certain proportion of 0° plies are embedded into the soft skin at
certain intervals as additional reinforcing elements (crack-blocking zones). This
approach is mainly used in shear-bearing transportation aircraft wings.
(3) Film enclosure
A layer of adhesive film may be introduced in between laminate plies to
increase the interlaminar damage resistance or to reduce the interlaminar stress
concentration for easily impacted structures. Epoxy films (such as FM series
films) or thermoplastic films (such as HXT series films and PEEK film) can be
inserted between the carbon fiber plies to increase damage resistance. A new
generation of interlaminar enclosed films can be made by spraying toughened
particles on prepreg tapes, which can largely increase the interlaminar tough-
ness and compression strength after impact without increasing the thickness
between plies.
(4) 3D reinforcing (Z-axis reinforcing)
3D reinforcing is mainly used to inhibit the delamination growth caused by
impacts and to increase the composite structural damage tolerance. Approaches
include reinforcement braiding in the thickness direction (such as 3D braiding
and Z-axis knitting performed in combined RFI and RTM processing), as well
as fasteners and Z-axial pin joining. Among these methods, dry/knitting and 3D
braiding/RTM show great potential for applications in improving damage
tolerance.
Z-pin joining is another mechanical joint for Z-axial reinforcement, other than
the use of metal fasteners. A foam preform (made of FM) containing small
strong carbon/epoxy pins is placed on a laminate structure. These preforms will
be pressed into laminates during hot pressing. This approach can be used for
reinforcement, locally or over the entire component. This approach can also be
used to replace metal fasteners used to fix frame construction. Test results
indicate that it may be possible to reduce delamination size and increase
damage tolerance with this method.
(5) Other approaches for durability/damage tolerance improvement
Three approaches can be used for durability/damage tolerance improvements:
use of special designs methods to inhibit damage growth and increase residual
strength. On the basis of the analysis of failure mechanisms of laminated
structures containing damage (including impact damage), the composite per-
formances can be improved to increase their damage tolerance. Namely, the
residual strength can be improved when laminates contain damage of the same
4 Composite Structure Design and Analysis 525
4.9.5.1 Review
In fact, the most effective approach to increase composite structure damage toler-
ance and damage resistance is to develop new material systems with high damage
tolerance and high damage resistance. The traditional method to evaluate the
composite damage tolerance is use of compression strength after impact (CAI) as
detained in NASA RP1142 and SACMA SRM 2R−94; in the recent studies, it has
been indicated that CAI obtained in such approaches can only evaluate damage
resistance, rather than damage tolerance, which needs a proper evaluation approach
[38–45].
In ASTM D 3878−07 Standard Terminology for Composite Materials, defini-
tions on damage resistance and damage tolerance are given, and their differences
are also discussed.
have a smaller effect on the damage area and damage width. The residual compression
strength is directly related to the damage width (damage area); thus, a knee point on the
impact energy (dent depth) versus compression failure strength (strain) curve will
result. Before the knee point, the compression failure strength (strain) will rapidly
decrease as the impact energy (dent depth) increases. After the knee point, the com-
pression failure strength (strain) will not change any further or only show a small
change.
On the basis of the physical consequences of a knee point, typical values taken
from areas adjacent to the knee point can be used to characterize the damage
resistance and damage tolerance of composite laminates.
For a long time, CAI values obtained from NASA RP 1142 or SACMA SRM 2R
−94 have been considered to be the main specifications for characterizing damage
tolerance. In the NASA standard, an impactor 12.7 mm in diameter with an impact
energy of 27 J (about 4.45 J/mm) is used. In the SACMA standard, an impactor
16 mm in diameter with an impact energy of 6.67 J/mm is used.
Here, the obtained damage tolerance values represent the corresponding com-
pression failure strength obtained under testing conditions of 27 J (NASA standard)
or 6.67 J/mm (SACMA standard). In fact, the impact energy cannot reflect damage
parameters, such that values derived from these methods cannot be used to evaluate
the damage tolerance behavior of composite systems perfectly. In Fig. 4.101, the
relationship between the damage width versus compression failure strain of two
different toughened composite systems is shown. In terms of damage tolerance, the
composite systems IM6/3501−6 (brittle epoxy) and IM7/8552−1 (toughened
epoxy) have similar damage tolerance behavior, but their CAI values are quite
different [15]. According to the above analysis, the composite systems with the
higher CAI value at knee point may produce a larger maximum damage area (or
diameter) than composite systems with a lower CAI values. Composite systems
with higher CAI values may show lower CAIT as given in Fig. 4.102. In
Table 4.23, some test results are listed.
4 Composite Structure Design and Analysis 531
4.10.1 Introduction
(3) The service moisture absorption content, final moisture content, steady condi-
tions, and hygrothermal allowables of chosen materials and configurations
should be determined.
Aircraft composite structures are susceptible to impact damage from tool dropping,
runway detritus, hail stones, and ground service vehicles. Other considerations
include lighting strikes, bird-strikes, and bullet damage. In general, lightning strikes
will result in visible damage and local ablation of the composite structures,
Here, the aging environments used to study composite structures are discussed,
including corrosive liquids (fuel, hydraulic fluid, and antifreeze), ultraviolet radi-
ation, weathering, and sand and rain erosion.
Table 4.25 Maximum and minimum air-temperature of every month in typical areas °C
Area 1 2 3 4 5
Month Dry-cold Basic warm Hot–wet Warm coast Wet-hot
inland coast
Tmax Tmin Tmax Tmin Tmax Tmin Tmax Tmin Tmax Tmin
1 −7.4 −21.8 1.4 −9.9 9.6 3.2 3.5 −2.6 20.0 12.0
2 −4.0 −18.4 2.5 −8.5 9.8 3.9 4.4 −0.2 20.0 12.0
3 4.0 −10.4 10.5 −0.5 15.7 6.8 8.1 1.4 24.0 16.0
4 12.6 −0.8 18.6 6.0 20.2 12.8 11.6 5.7 30.0 22.0
5 18.9 5.5 27.0 15.0 26.7 18.8 18.2 11.2 31.0 23.0
6 23.5 10.1 30.0 19.0 28.3 20.5 22.1 16.7 32.0 24.0
7 25.5 12.6 30.8 21.5 31.9 23.5 25.5 21.3 33.0 24.0
8 23.8 11.6 29.5 20.5 31.1 22.5 28.2 23.7 34.0 24.0
9 18.2 5.4 23.5 4.5 29.3 19.5 25.1 19.7 32.0 22.0
10 13.2 −2.4 18.0 6.0 21.4 4.1 20.2 15.4 28.0 22.0
11 1.6 −13.2 9.30 −3.0 16.4 9.3 13.0 6.2 24.0 18.0
12 0.2 −20.2 2.0 −10.0 11.5 3.5 6.9 0.5 21.0 13.0
Table 4.26 Statistical results of relative humidity in basic warm and dry-cold areas
Area Relative humidity/%
Average monthly Average monthly Monthly Monthly
maximum minimum maximum minimum
Basic warm 85 30 100 0
Dry–cold 78 25 100 0
534 Z. Shen et al.
One important key point for composite structure design is to consider the influence
of the hygrothermal environment on structural performance. The resin matrix has
the ability to absorb moisture, and moisture diffusion can result in a distribution of
moisture content in the structure. Thus, both the anti-corrosion resistance of fibers
and the glass transition temperature Tg might decrease. The structural stiffness and
strength of the composites might also be reduced through these effects. At all stages
of material and configuration selection, detailed design and testing of the composite
structure should account for environmental response of the system [1, 2, 13].
Theoretic Predictions
@C @2C
¼ Dz 2 t [ 0; z 2 ½h=2; h=2 ð4:102Þ
@t @Z
Cðz; 0Þ ¼ C0
where
C moisture concentration;
C0 , C1 initial concentration and equilibrium concentration;
Dz moisture diffusivity through the thickness direction;
t time;
h laminate thickness;
z coordinate in the thickness direction.
The total moisture content is:
where
M(t) moisture of laminate;
M0, M∞ initial moisture content and equilibrium moisture content
kn ¼ 2p þ n; n ¼ 0; 1; 2. . .
h2 Mðt1 Þ Mðt2 Þ 2
D¼p pffiffiffiffi pffiffiffiffi ð4:104Þ
4ðM1 M0 Þ t1 t2
The disadvantage of this model is the use of fixed boundary conditions, and the
shape of the moisture curve at the initial stage are not accurately described.
Therefore, the initial moisture absorption is greatly over estimated and the diffu-
sivity is under estimated.
(2) Vapor boundary model
When a solid absorbs or desorbs water vapor from the atmospheric environment,
the Fickian model produces large deviations. A proportionality constant F may be
introduced. Hence, F is defined as the moisture absorption gradient and is pro-
portion to the difference between the actual surface concentration and equilibrium
concentration. The diffusion equation is unchanged from that given in (4.102);
however, the boundary conditions are modified as:
C ðz; 0Þ ¼ C0 ð4:105Þ
D @C
Cð h=2; tÞ ¼ C1 þ ð h=2; tÞ
F @z
Thus, in the limit F ! ∞, the vapor boundary model degrades to the Fickian
model.
The moisture content is:
X
1
MðtÞ ¼ M1 ðM1 M0 Þ
n¼0
( ) ð4:106Þ
2
2 sin bn
expð4b2n Dt=h2 Þ
bn ðsinbn cosbn þ bn Þ2
bn : b tanb ¼ hF=2D n ¼ 0; 1; 2. . .
n p o
bn 2 np; þ np
2
as F
D ! 0,bn ! np;
as DF ! 1, bn ! p2 þ np;
An iteration of the following form can be adopted.
Table 4.27 Equilibrium moisture content and diffusivity of materials in specified environment
Material D D me
30 °C, 95%RH 50 °C, 95%RH 50 °C, 95%RH
T300/5405 1.15 10−7 mm2/s 3.788 10−7 mm2/s 0.85%
T300/QY8911 3.50 10−7 mm2/s 7.043 10−7 mm2/s 1.35%
F ¼ 2ðC11C0 Þ
ð4:107Þ
M0 ðt2 t1 Þ Mðt1 Þ Mðt1 Þ 2 Þt1
2M
t2
0
t1 t2 þ t1 þ t2 t1 t2Mðt
ðt2 t1 Þ
This model corresponds to a situation in which moisture enters the material from
the ambient environment. The initial rate parameters can be obtained easily and
have clear physical meanings. Results predicted by this model for high relative
humidity are consistent with experimental findings, and at low relative humidity the
predictions can satisfy engineering requirements.
In conclusion, for environments with low relative humidity, the moisture
absorption diffusion process in composite laminate can be described by the Fickian
model; however, under environments of high relative humidity, the moisture
absorption diffusion process is better described by the vapor boundary model owing
to swelling of the composite laminate.
t1 eC=T2 /2
k¼ ¼ ð4:108Þ
t2 eC=T1 /1
where
K accelerated time coefficient;
t1 actual exposed time;
t2 time after acceleration;
T1 /1 temperature (°C) and relative humidity of actual exposure environment;
T2 /2 temperature (°C) and relative humidity of accelerated environment.
(2) The accelerated time coefficient K can be estimated from the ratio of dif-
fusivity in the different environments:
t1 D2
K¼ ¼
t2 D1
D ¼ D0 expðC=T Þ ð4:109Þ
where
4 Composite Structure Design and Analysis 539
Table 4.28 Diffusion constants of different materials under different environments at room
temperature
Material T300/1034 AS/3501-5 T300/5208 934 3501-5 5208
Constant D0 C D0 C D0 C D0 C D0 C D0 C
Distilled 16.3 6211 768 7218 132 6750
water
Saturated 5.85 6020 5.38 6472 6.23 5912
brine
Moist air 2.28 5554 6.5 5722 0.57 4993 4.85 5113 16.1 5690 4.19 5448
where
Tg glass transition temperature of matrix under certain moisture content;
bm wet swelling coefficient of matrix under certain moisture content;
bf wet swelling coefficient of fiber under certain moisture content, usually equal
to zero;
Vf fiber volume content under certain moisture content;
Tgm glass transition temperature of matrix under certain moisture content;
Tgf glass transition temperature of fiber under certain moisture content.
Tgf —glass transition temperature of fiber under certain moisture content.
The changes of glass transition temperature with moisture content for some
composite material systems are given in Table 4.29. The experimentally determined
variation of the glass transition temperatures of three material systems is presented
in Fig. 4.103. This figure shows that for polymer matrix systems the Tg declines by
approximately 25 °C for at a moisture content of 0.5%. For further increases in the
moisture content over 1.2% there is only a slight decrease of Tg. For cyanate esters
matrix composites, Tg declines by approximately 20 °C when the moisture content
is greater than 0.3%. For a BMI matrix composite, moisture content has hardly any
effect on Tg.
(2) Wet swelling coefficient and thermal expansion coefficient
The wet swelling and thermal expansion coefficients for some materials are
shown in Table 4.30. The change of the thermal expansion coefficient with mois-
ture content can be predicted by the following equation:
1
Tgw T 2
aðTÞ ¼ aðRTÞ ð4:111Þ
Tgd TRT
where
a(T) thermal expansion coefficient at temperature T under a certain moisture
content;
a(RT) thermal expansion coefficient at room temperature;
Tgw glass transition temperature at certain moisture content;
Table 4.29 The change of glass transition temperature with moisture content for some materials
Material 914 T300/914C T300/5222 T300/4211 4211
4 Composite Structure Design and Analysis
Moisture 0 0.9 0.6 2.9 4.7 0 0.9 1.6 0 1.0 0 1.0 0 1.02
content/%
Tg/°C 208 188 178 160 137 219 169 153 245 195 99 84 156 132
541
542 Z. Shen et al.
Table 4.30 Wet swelling and thermal expansion coefficients of some materials
T300/5208 B/5505 AS/3501 Scotch/1002 Kevlar/epoxy
aL/10−6K−1 0.02 6.1 −0.3 8.6 −0.4
aT/10−6K−1 22.5 30.3 28.1 22.1 79.0
bL 0.0 0.0 0.0 0.0 0.0
bT 0.6 0.6 0.44 0.6 0.6
where
T operation temperature;
Tg glass transition temperature;
Tr reference temperature or room temperature.
When experimental results on the influence of hygrothermal environment on
composite mechanical performance are lacking, the following equations may be
used to make predictions:
Xt =Xt0 ¼ Vf =Vf0 ðT Þ0:04
E1 =E10 ¼ Vf =Vf0 ðT Þ0:04
XC =XC0 ¼ Vf =Vf0 ðT Þ0:04 ðT Þ0:05
h i0:5
E2 =E20 ¼ Tg T = Tg0 Tr ð4:113Þ
h i0:5
G12 =G012 ¼ Tg T = Tg0 Tr
S=S0 ¼ ðT Þ0:2
h i0:5
m=m0 ¼ a Tg T = Tg0 Tr
where
the superscript 0—represents the dry state;
a—parameter related to moisture content, for carbon fiber-reinforced composite
(when moisture is not more than 1.0%, a 1:0);
Vf—fiber volume content.
Tension and compression tests have been performed on different multi-laminate
T700S/5405 material systems of different thicknesses immerged in 70 °C distilled
water for 3 weeks. The ply ratio of these laminates was 0°-plies 40%, 45°-plies
50%, and 90°-plies 10%. The test results are shown in Table 4.32. The following
conclusions could be drawn:
① Moisture absorption alone has little influence on the tension strength and
modulus.
② Moisture absorption alone has little influence on the compression
modulus; however, the compression strength of 1.5-, 2.5-, and
3.0-mm-thick laminates dropped 2.2%, 8.6%, and 5.0%, respectively.
4 Composite Structure Design and Analysis 545
Table 4.32 Tension and compression properties of T700S/5405 system laminate at different
thickness
Properties Nominal thickness of specimen
1.5 mm Cv / 2.5 mm Cv / 3.0 mm Cv /
% % %
Tension strength/MPa Dry 1410.0 3.1 1151.8 7.5 1291.1 6.0
Wet 1326.0 6.6 124.8 3.9 1308.5 6.3
Tension modulus/GPa Dry 81.77 3.5 64.71 3.3 71.92 3.3
Wet 81.72 4.6 66.31 1.5 75.56 6.2
Tension Poisson’s ratio Dry 0.34 5.6 0.52 3.5 0.52 5.9
Wet 0.36 6.7 0.53 2.8 0.53 3.3
Tension extensibility/ Dry 2.03 5.0 2.09 6.6 2.11 5.1
% Wet 1.90 8.8 2.17 3.6 2.01 5.9
Compression Dry 560.28 7.9 716.81 7.1 752.66 6.9
strength/MPa Wet 548.07 6.8 654.99 4.1 715.10 6.6
Compression Dry 79.19 8.7 61.32 2.7 67.71 5.1
modulus/GPa Wet 80.97 10.3 64.00 9.6 73.79 5.2
Compression Poisson’s Dry 0.31 8.5 0.48 8.3 0.49 4.1
ratio Wet 0.32 9.4 0.48 13.4 0.50 10.6
Table 4.33 Mechanical properties of T300/QY8911-III system stitched laminate under different
thickness
Properties Layering and nominal thickness of specimen
Type A Type B
3.0 mm 4.0 mm 4.5 mm 3.0 mm 4.0 mm 4.5 mm
Tension strength/MPa Dry 631.03 565.16 601.80 498.55 671.29 707.15
Wet 672.38 573.61 609.30 – – –
Tension modulus/GPa Dry 63.21 58.59 56.88 54.97 68.11 69.85
Wet 64.03 64.92 62.86 – – –
Tension Poisson ratio Dry 0.41 0.49 0.42 0.42 0.54 0.38
Wet 0.42 0.52 0.45 – – –
Tension extensibility/% Dry 1.22 1.15 1.25 1.08 1.16 1.21
Wet 1.20 1.05 1.16
Compression Dry 613.41 659.89 610.26 435.86 648.97 693.03
strength/MPa Wet 519.59 637.94 557.81 – – –
Compression Dry 51.81 52.99 51.42 48.20 65.31 66.57
modulus/GPa Wet 54.38 54.07 52.48 – – –
In-plane shear Dry 247.92 276.32 289.41
strength/MPa Wet 219.92 271.09 294.88
In-plane shear Dry 13.13 16.77 17.47
modulus/GPa Wet 15.37 18.11 18.25
Flexural strength/MPa Dry 713.33 699.77 74.61 783.13 801.58 765.60
Wet – 700.75 – 793.79 – 771.36
Flexural modulus/GPa Dry 47.93 40.66 48.30 52.60 44.18 52.03
Wet – 40.58 – 53.41 – 51.99
Flexural failure Dry 10.20 10.70 7.38 10.55 10.58 7.86
deformation/mm Wet – 10.73 – 10.45 – 8.00
Interlaminar shear Dry 53.61 64.28 60.69 57.03 66.23 67.37
strength/MPa Wet 53.19 62.74 58.48 – – –
Table 4.34 Moisture contents of nine materials systems after immersion in 70 °C distilled water
for three weeks
Material Moisture content/% Material Moisture content/%
T700S/5428 0.58 T800/QY9511 0.95
T700S/5429 0.46 T700S/QY9511 0.58
T700S/5405 0.67 T300/QY8911 1.11
T700S/5228 1.00 Stitched T300/QY9512 1.45
T300/5405 0.66
Table 4.35 Test results of tension and compression and tension and compression with open hole
specimens for nine material systems
Material Test Test type
environment Tension with Tension Compression Compression
an open hole with an open hole
rkt rt Et mt rkc/MPa rc Ec
/MPa /MPa /GPa /MPa /GPa
T700S/5428 −55 °C 628.3 991.5 63.36 0.46 347.2 574.5 53.50
23 °C, 50% 675.5 1071.4 59.58 0.51 332.2 549.4 59.39
RH
80 °C, 50% 694.8 1045.6 3.48 0.54 283.8 466.3 54.04
RH
125 °C 709.4 908.4 72.12 0.55 248.0 419.2 54.54
150 °C 717.5 934.0 68.29 0.48 204.1 358.4 54.98
170 °C 728.5 870.1 69.32 0.53 146.0 260.1 55.31
T700S/5429 −55 °C 622.0 1037.2 59.78 0.48 348.4 543.1 50.31
23 °C, 50% 636.9 1129.1 58.51 0.49 304.1 491.3 50.06
RH
80 °C, 50% 626.9 1047.8 57.03 0.53 276.6 495.8 50.73
RH
125 °C 654.3 990.5 68.54 0.60 263.3 433.4 47.63
150 °C 637.1 900.0 68.66 0.59 165.3 360.0 53.55
170 °C 636.4 863.0 68.27 0.53 135.4 234.7 44.57
T700S/5405 −55 °C 706.8 1261.2 65.60 0.50 411.8 525.9 54.84
23 °C, 50% 732.5 1268.7 68.48 0.53 355.4 493.3 51.08
RH
80 °C, 50% 874.0 1151.4 63.46 0.56 306.8 485.9 52.78
RH
T700S/5405 125 °C 722.1 1045.8 69.81 0.60 239.0 411.5 50.69
150 °C 683.2 1004.3 64.55 0.54 120.8 409.6 52.42
170 °C 668.9 942.3 50.69 0.45 92.0 250.7 51.04
T700S/5228 −55 °C 522.9 905.7 56.36 0.53 363.0 528.4 45.44
23 °C, 50% 573.5 832.6 49.99 0.51 319.3 497.6 44.15
RH
80 °C, 50% 577.7 899.3 52.65 0.54 258.1 500.5 44.58
RH
125 °C 532.5 843.6 69.08 0.52 180.2 436.3 44.32
150 °C 507.7 754.2 59.76 0.50 108.2 315.1 47.22
170 °C 533.9 685.5 59.70 0.52 81.3 215.4 48.60
T300/5405 −55 °C 355.2 768.2 64.28 0.49 449.9 644.4 61.37
23 °C, 50% 340.3 702.4 69.01 0.50 356.1 570.9 56.12
RH
80 °C, 50% 352.3 688.8 61.60 0.52 299.3 487.7 54.14
RH
(continued)
4 Composite Structure Design and Analysis 549
Fig. 4.107 Influence of temperature on tension strength with an open hole for nine material
systems
Fig. 4.108 Influence of temperature on tension strength for nine material systems
Fig. 4.109 Influence of temperature on compression strength with an open hole for nine material
systems
Fig. 4.110 Influence of temperature on compression strength for nine material systems
EC /% 88.4 100 85.3 78.8 74.8 82.0 88.2 100 95.4 99.1 101.2 99.7
Residual Stitched T300/QY9512
performance −55 °C 23 ° 80 ° 125 °C 150 °C 170 °C −55 °C 23 °C, 80 °C, 125 °C 150 °C 170 °C
C, C, 50% 50%
50% 50% RH RH
RH RH
rkc /% 110.3 100 73.0 37.2 22.4 16.6
rc /% 104.8 100 78.6 41.4 25.3 4.6
EC /% 106.9 100 95.1 93.5 101.3 90.1
553
554 Z. Shen et al.
The hygrothermal environment not only affects the physical and mechanical
properties of composite laminate, but also affects failure modes. The failure modes
at low temperatures and in the dry state are related to basic failure of the matrix
itself. The failure modes of wet composites at room temperature involve hybrid
failure of a matrix/interphase. The failure modes of wet composites at elevated
temperature involve failure of the fiber/matrix interface.
eT ¼ aDT
eC ¼ 0 ð4:114Þ
ZZ
fRT ge ¼ Dt ½BT ½DfeT gdxdy
0 1 0 1
e1T a1 DT
@ e2T A ¼ @ a2 DT A
0 0
0 1 0 1
e1C b1 DC
@ e2C A ¼ @ b2 DC A
0 0
ZZ
fRT ge ¼ Dt ½BT ½T T ½QfeT gdxdy ð4:115Þ
e
ZZ
e
fRT g ¼ Dt ½BT ½T T ½QfeT gdxdy
e
ZZ
fRC ge ¼ Dt ½BT ½T T ½QfeC gdxdy
e
where
a —thermal expansion coefficient of material, 1/°C;
DT —increment of temperature, °C;
DC —increment of moisture content;
a1 —longitudinal thermal expansion coefficient of laminate, 1/°C;
a2 —transverse thermal expansion coefficient of laminate, 1/°C;
b1 —longitudinal wet swelling coefficient of laminate;
b2 —transverse wet swelling coefficient of laminate;
eT —initial strain caused by temperature;
eC —initial strain caused by moisture absorption;
{RT}e —equivalent thermal load at element node;
{RC}e —equivalent wet swelling load at element node.
The strength and stiffness performance of resin matrix composites will vary con-
siderably with extended usage-time, especially in certain hygrothermal environ-
ments. Hygrothermal aging of fiber-reinforced composites is a gradual degradation
process caused by the combined action of moisture uptake, temperature, and stress.
Fibers and the fiber/matrix interface are degraded by physical/chemical reactions.
During the moisture absorption process, a swelling stress will be introduced to the
interior of composites. A greater swelling stress might be introduced owing to rapid
desorption of the surface layer of wet structures under thermal spiking. Under this
repeated interior stress, at a certain threshold stress, cracking will occur followed by
crazing. The moisture re-absorption and re-desorption rates will be affected by the
556 Z. Shen et al.
crazing, and finally, macroscale cracks will form. Therefore, the hygrothermal
aging response of any selected composite material system should be investigated at
the design stage. However, this theoretical analysis is difficult owing to uncer-
tainties of the environmental and the coupling effect between hygrothermal stress
and exterior loading. In general, experimental methods are used to study these
factors based on ground environmental aging, accelerated laboratory aging and
aging in actual flights. The experimental data are globally analyzed to obtain design
criterion for hygrothermal aging [1, 2, 13].
Tgd Tg0
a1 ¼
Tgf Tg0
m0 m
a2 ¼
m0 mf
where
am ′ matrix thermal expansion coefficient after accounting for
aging effects;
am matrix thermal expansion coefficient at room temperature
in dry state;
a1 degradation parameter of the crosslinking mechanism;
a2 parameter of matrix mass change;
Tgd measured glass transition temperature at room temperature
in dry state;
T0g measured glass transition temperature at initial aged state;
Tgf measured glass transition temperature at final aged state;
m0 mass of a small neat matrix specimen at initial aged state;
m mass of a small neat matrix specimen at room temperature
and dry state;
mf mass of a small neat matrix specimen at final aged state.
Daa1 , Daa2 , naa1 , naa2 fitting parameters based on the change of a1 and a2 data.
The matrix wet expansion coefficient can also be modified by an analogous
methodology.
responses of BMI matrix composite are shown in Figs. 4.114 and 4.115. When the
moisture content was less than 0.6% (corresponding to 70 h of aging), there was
little change in the interlaminar shear strength; however, for a moisture content
greater than 0.6%, the interlaminar shear strength showed a marked decrease. The
interlaminar shear strength of the moisture saturation state dropped by approxi-
mately 50% compared with that of the dry state.
4 Composite Structure Design and Analysis 559
This may be explained by the reaction between the water and matrix requiring a
certain time. Thus, controlling moisture content is a design criterion for application
of MBI matrix composites in structures.
Physical Aging
When a polymeric matrix material is used below its glass transition temperature for a
long time, the mechanical properties will change markedly. This change is termed
physical aging. During the physical aging process, the material becomes stiffer, with
decreased compliance and an increased modulus. The physical aging responses of
resin matrix composites have received considerable research attention. The results of
various investigations have shown that the matrix-dominated properties of continuous
fiber-reinforced composite (e.g., the shear and transverse responses) are most seri-
ously affected by physical aging in a similar manner to that of a pure polymer.
d lg ate
l¼
d lg te
where
S(t) compliance at time t;
S0 initial compliance;
te aging time;
s relaxation time;
T time;
teref reference aging time;
ate aging time shift factor at aging time;
b shape parameter;
l shift rate.
4 Composite Structure Design and Analysis 561
Note from Eq. (4.117) that if the initial compliance S0, shape parameter b, shift
rate l and relaxation time s at reference aging time teref of the polymer are known,
then the compliance at any time may be determined. If te > teref, then ate < 1;
otherwise, ate > 1. In the case of te > teref, the relaxation time at te is also greater
than that at the reference aging time [s (te) > s(teref)]. This relation shows that at a
given time the modulus of the material is higher and the compliance is lower [S(t;
te) < S(t;teref)].
For most polymers, if the material is being used at temperatures close to its Tg,
the material changes into an equilibrium state in a relative short time. The time
required to achieve the equilibrium state is known as the equilibrium aging time.
The shift factor, l characterizes the influence of aging on material properties. For
the same aging time a larger value of l indicates a smaller compliance change.
Thus, the shift factor l can be used as a screening parameter for selection of
materials. Materials with larger l values should be chosen. In general, experimental
results have shown that before the aging equilibrium l 1 and after the aging
equilibrium l 0.1. These results indicate the dramatic change of materials in the
aging equilibrium state. Thus, in the design stage of polymer matrix composite
structures, materials with larger l value should be selected, while avoiding aging to
an equilibrium state during the full life period, particularly in structures for use in
high temperature applications.
(2) Influence of aging temperature
Although the shift factor l over a large temperature range is constant, in fact,
both l and ate are functions of temperature. In general, the relationship of the time
temperature-aging time shift factor can be expressed as:
where aT is a time temperature shift factor, i.e., a function of temperature and aging
time at temperatures below Tg.
The relationship between aT and l (T) can be expressed as:
lðT2 ÞlðT1 Þ
atTe21 =T2 te2
¼ ð4:119Þ
atTe11 =T2 te1
where atTe11 =T2 —time temperature shift factor between temperature T1 and T2 at aging
time te1 .
Thus, if the l(T) value or its expression and the time temperature shift factor at a
single aging time are given, then at at any aging time can be calculated. Therefore,
the shift factor, l, and the time temperature shift factor, at , have an effect on aging
time.
562 Z. Shen et al.
where
½S—effective compliance matrix;
[S]—compliance matrix with respect to the fiber coordinate system;
[T]—transformation matrix.
The elastic stress–strain relation under in-plane loading is given by:
2 3 2 3
exx rxx
4 eyy 5 ¼ S 4 ryy 5 ð4:121Þ
exy rxy
Experimental studies of polymer matrix composites have shown that the trans-
verse compliance S22 and the shear compliance S66 are related to time temperature
and subject to physical aging. Their values can be determined by Eq. (4.122), such
that in a functional form:
Note from Eq. (4.122), the transverse and shear compliance of composite uni-
directional laminates are independently described by four viscoelastic parameters,
namely the initial compliance, shape parameter, relaxation time at a given reference
aging time, and shift factor. For any given material the four independent parameters
may be determined by short-term aging tests in the laboratory.
Equation (4.121) is rewritten, accounting for time relativity as:
2 3 2 3
exx ðtÞ rxx
6 7 6 7
6 eyy ðtÞ 7 ¼ ½sðtÞ6 ryy 7;
4 5 4 5 ð4:123Þ
exy ðtÞ rxy
Sij ðtÞ ¼ f ðh; Sij Þ
Table 4.39 Viscoelastic parameters, transverse and shear compliance of IM7/8320 material
system
Viscoelastic parameter S22 S66
l 0.77 0.93
b 0.416 0.456
s 1.19 106 s 4.31 105 s
S0 750 10−9 Pa−1 1364 109 Pa−1
teref 3.24 104s 3.24 104 s
Elastic parameters: S11 = 5.75 10−91/psi t12 = 0.348
required test methods, were exposed for approximately 13 years. The exposure tests
included ground exposure, flight travel\exposure, and accelerated laboratory aging,
which simulated the change of temperature, humidity, and pressure during aircraft
flight. The ground-based exposure tests were performed at Dallas, NASA Dryden,
Honolulu, and Wellington. Aloha Airlines, Air New Zealand Ltd., and Southwest
Airlines were selected for flight exposure studies. The laboratory accelerated aging
results are shown in Fig. 4.117.
On the basis of global analysis of the test results, the following recommendations
for composite structural design were proposed:
① The tension and flexure strength at room temperature after aging for the
three materials showed a slight overall increase. At elevated tempera-
tures, the results were mixed. For the T300/5209 and T300/934 systems
the flexure and tension strength decreased slightly. For the T300/5208
system, both these properties were greater their baseline strength. The
T300/934 tension strength also increased. However, in all cases, the
differences were relatively small.
② Room-temperature compression strength dropped in general. At the end
of 10 years’ exposure, all three materials showed decreases of approx-
imately 30%. The elevated temperature residual strength was likely
seriously decreased; however, the exact test data could not be deter-
mined owing to the grab-tab failure.
③ The short-beam shear displayed a peculiar pattern for residual strength in
both room and elevated temperature tests on all three material systems.
The drop of the shear strength was largest after 1, 2, and 3 years of
exposure; lesser degradation was found after 5 years exposure; however,
the room-temperature residual strength increased slightly, and strength at
elevated temperatures remained at or near their baseline levels after
10 years of exposure.
4 Composite Structure Design and Analysis 565
In general, the designed life of fighter aircraft is 5000 flight hours, or 20–30 years.
Complete simulation of both the mechanical and environmental loading history is
the most credible evaluation; however, this would be impractical. Acceleration of
the actual temperature/humidity time history can be used to obtain the accelerated
hygrothermal aging results. Therefore, a large amount of comparable data can be
accumulated and the development period for structures can be shortened and the
test costs reduced. Furthermore, individual test results can be very easily interpreted
and estimated.
In this section, based on the flight environment and service mission of aircraft in
China, accelerated hygrothermal aging and test results for composite components of
certain fighter aircraft are introduced.
(1). Basis for Establishing Scheme
On the basis of a typical mission profile, involving 5000 flight hours over
20 years, an accelerated aging program is developed.
(2). Developing Requirements and Basic Rules
① Accelerated tests should yield the same results for composite degrada-
tion and residual strength compared with that resulting from the real-time
history, or give more conservative results than those from experiments.
② The actual aircraft usage environment should be reflected reasonably.
③ For accelerated aging, the response of thermal spiking caused by aero-
dynamic heating should be considered because an elevated temperature
environment will have a considerable influence on composite properties
over the long term.
④ The greatest test acceleration may be achieved by compressing the
simulated ground standing time as much as possible. The ambient
environmental exposure over 20 years may be simulated by accelerated
tests over one year.
⑤ The selected accelerated conditions should not have any additional
effects on composites. For BMI matrix systems, an accelerated temper-
ature of 70 °C is appropriate.
566 Z. Shen et al.
Table 4.41 Accelerated hot–wet spectrum aging results of laminates and sandwich construction
Strength properties Baseline value Accelerated hot–wet spectrum aging
Compression/MPa 525.3 500.4
Tension/MPa 606.7 660.1
SBS/MPa 68.2 46.7
Interlaminar shear/MPa 4.3 13.3
Bearing/MPa 104.1 851.1
Flatwise tension/MPa 2.1 2.1
Core shear/MPa 1.1 1.0
Core shear modulus/MPa 39.6 38.7
(5) Properties of the sandwich construction after aging were unchanged with the
exception of the facing modulus.
immersion). On the basis of the acceleration principle, 1 year of the natural envi-
ronment could be simulated by 28 days of in the accelerated environment. The
effects of midday solar radiation in June, July, August, and September were con-
sidered. If the irradiation time each day is 2 h, and the irradiation temperature is 50
°C, and then the total irradiation time of each month is 60 h, or 240 h in 4 months.
Therefore, 1 year of ambient environmental aging could be simulated using our
accelerated testing method for 38 days. This scheme involved the following steps:
4 þ 1 þ 4 þ 1 þ 4 þ 2 þ 4 þ 2 þ 4 þ 2 þ 4 þ 1 þ 4 þ 1 ð4:124Þ
where
4—Immersion in 70 °C/100% RH (70 °C distilled water immersion) for 4 days;
• —50 °C heating for 1 day;
②—50 °C heating for 2 days.
For an aging duration of 3 years, the above spectrum may be repeated three
times.
Ambient environmental aging for 1 and 3 years was performed at an environ-
mental experiment field in Hainan Province in China. The corresponding acceler-
ated aging was performed in a laboratory of the Aircraft Strength Research Institute.
The test results are shown in Tables 4.42 and 4.43. The test results indicated:
① Ambient aging has no measurable influence on tensile and compressive
strength.
② Ambient aging has a large influence on short-beam shear (SBS) strength.
The SBS strength after 1 and 3 years of ambient aging decreased by
26.5% and 37.0%, respectively. The shear strength of single- and
double-lag bonds after 1 and 3 years of ambient aging decreased by
45.1% and 12.7%, and 51% and 30.2%, respectively. The tension–shear
strength after 1 and 3 years of ambient aging decreased by 9.8% and
30.2%, respectively.
Table 4.42 Ambient aging and accelerated laboratory aging results of laminates
Strength Baseline 1 year 3 years 3 years 1 year 3 years
properties value ambient ambient ambient accelerated accelerated
aging aging aging aging aging
(unpainted) (painted)
Compression/MPa 525.3 533.3 659.7 695.6 482.8
Tension/MPa 606.7 662.4 682.6
SBS/MPa 68.2 50.1 43.0 65.9 55.1
Bonded I/MPa 37.7 20.7 18.5
Bonded II/MPa 32.9 26.3
Tension– 4.3 12.9 11.8
shear/MPa
570 Z. Shen et al.
Table 4.43 Ambient aging and accelerated laboratory aging results of sandwich construction
Strength properties Baseline 1 year 3 years 3 years 1 year 3 years
value ambient ambient ambient accelerated accelerated
aging aging aging aging aging
(unpainted) (painted)
Flatwise 3.62 3.49 2.96
compression
(H = 19)/MPa
Flatwise 2.88 3.03 2.97
compression
(H = 44)/MPa
Core shear 1.21 1.23 0.77 0.76
(H = 19)/MPa
Core shear 1.11 1.04 1.03 1.06
(H = 44)/MPa
Core shear modulus 39.2 33.7 45.0 44.4
(H = 19)/MPa
such as its size and velocity. When the angle of incidence of a raindrop is
90° and the raindrop velocity is more than 200 m/s an anti-rain erosion
protective coat should be applied to the surface of the composite
structure. At the design stage, the layout of components should attempt
to minimize the area of impacted surfaces and lower the angle of rain
drop incidence in addition to protective coatings. Anti-erosion paint and
metal or ceramic protective coatings may be adopted as anti-rain water
erosion measures.
(3) Corrosive Control
There are two main principles for controlling composite corrosion caused by
environmental media.
(1) Enhance innate material corrosion resistance: In some cases it is possible to
improve the crystallinity, tropism grade, or crosslinking density of com-
posites. The matrix compactness can also be enhanced to reduce the dif-
fusion coefficient and penetrative coefficient of the medium. A surface
cleanup solvent may be used to enforce the adherence strength between the
reinforced fiber and matrix, reduce the interface clearance, and enhance
impermeability.
2) Use of protective coatings: A protecting coat is sprayed on the composite
surface to avoid direct corrosion of the composite by environmental media.
(4) Biological Corrosion and its Control
Biological corrosion occurs mainly at the fuel box position of aircraft structures.
The combination of moisture and other impurities in fuel can provide appropriate
conditions for biological organisms to grow. The main microorganisms that might
affect composites are germs, epiphyte, and mildew. Such microorganisms may
reproduce and excrete acidic substances, such as lactic acid and grass fungus. These
acidic substances might react with composites. Biological corrosion of composite
structures used in a sea environment can pose a threat. Composite destruction
caused by oceanic organisms can occur. The composite might be bitten away by
hexapods and chisel-ship worms, and the above-mentioned microbe encroachment
can become more serious. The following steps should be taken to control biological
corrosion:
① Fuel quality should be controlled. The content of moisture and impurities
in fuel and possible pollution during the fuel transport process should be
minimized to remove the conditions necessary for microbial growth.
② An effective drainage system should be designed in the fuel tank.
A fluent-fuel mouth should be installed at the lowermost position and
water should be drained at intervals.
③ Anti-bioerosion protective coatings may be applied. A coating in com-
mon use is SF-9 epoxy, which can effectively prevent biological growth.
574 Z. Shen et al.
When two types of materials with different electrode potentials are directly con-
nected or in contacted through an electrolyte, accelerated corrosion might be caused
in the lower potential material. This is known as galvanic or electrical dipolar
erosion. Carbon has good electrical conductivity and a relatively high electrode
potential. Carbon fiber-reinforced composites under general environmental condi-
tions show inert behavior similar to that of noble metals with high electrode
potentials. Thus, when carbon fiber-reinforced composites are joined with metal the
cathode-like behavior could accelerate corrosion of the metal. The electrode
potential difference between carbon fiber-reinforced composite and most metals is
0.5–1.0 V, and in some cases, may reach as high as 1–2 V. Therefore,
anti-electrical dipolar erosion steps must be adopted in areas of connected metals
and composites. The generation of electrical dipolar erosion requires three condi-
tions: an electrode potential difference, electrolyte, and an electrical conductive
connection. Protective measures against galvanic erosion should consider these
three aspects.
(1) Structural Design
The accumulation of electrolyte can be prevented to a large extent by careful
structural design to avoid formation of corrosion batteries.
① Attention should be paid to structural seals, to avoid infiltration of
rainwater, fog, and seawater. Holes and places where contamination may
accumulate should be reduced. Countermeasures should be mounted at
positions subject to seepage.
② Small metallic elements surrounded by a large area of composite should
be avoided. Strict protecting steps should be adopted for mechanical
fastener joints.
③ The lumen and blind holes should be designed with perforation to pre-
vent condensation water cohesion.
(2) Selection of Materials
Consistent materials should be selected to prevent galvanic erosion.
① Anticorrosive materials and materials with a low potential difference
compared with that of the composite should be chosen. Pay special
attention to small parts such as fasteners.
② Insulated and closed down materials should be non-hygroscopic and
should not contain any corrosive components. When a single layer of
4 Composite Structure Design and Analysis 575
Protective coatings for composite components can not only improve the appearance
of faces but are also important for slowing moisture absorption and aging processes
of the material. The application of an anti-friction dope on the front structural
features can improve resistance to sand and rain erosion. Application of anti-friction
dope on interior surfaces can prevent direct contact of composites with metal and
avoid electrical dipolar erosion.
(1) Cover Coatings and their Effects
Different cover coat materials are used at different positions and fall mainly into
the following types:
① For interior surfaces and end faces a protective coat formed by priming
should be applied.
② For general exterior surfaces with ornamental protective coatings,
priming, sealants, transition priming, and surface paint may be applied.
③ For front structural features, an anti-friction and antiscouring protective
coating system formed by priming, an elastic anti-friction dope, and
surface paint should be applied.
④ For upright surfaces or other surfaces requiring anti-static protection,
anti-static protective coatings system formed by priming, and anti-static
or elastic anti-static dopes should be applied.
⑤ For the interior surfaces of the fuel tank an anti-static protective coat
system formed by priming, and anti-static and oil proof dopes should be
applied.
⑥ For exterior surfaces of the fuel tank, an electrical protective coating
system formed by priming, an electrical layer, and painting of the sur-
faces should be applied.
(2) Surface Dopes and Coating Systems
Currently, the most widely used dopes are epoxy- and polyurethane-based.
Epoxy dopes are strongly adherent, show low contraction, and good toughness.
4 Composite Structure Design and Analysis 577
Epoxy dopes are versatile and can be applied with many types of surface paint.
However, because this dope has poor gloss retention after film formation and it is
prone to pulverization, it is mainly applied in priming coats.
A large variety of polyurethane dopes are available, which shows excellent
performance, strong adherence, and high rigidity. Polyurethane dopes form bright
films with excellent oil and moisture proofing, heat endurance, and wear and
chemical resistance. Therefore, these dopes are widely applied in the aviation
industry.
Dopes and coating systems in common use for aircraft composite structures are
shown in Tables 4.45 and 4.46.
The relationships between three different types of aging methods are investigated to
determine an optimal method for studying aging response. It is important for the
method to not only reflect the real history of an aircraft but also to be convenient
and simple. On the basis of real investigations, the following recommendations are
given.
① For the composite structures of military aircraft, the accelerated testing
scheme shown in Fig. 4.119 may be considered to be a standard
accelerated hygrothermal scheme. An actual history of 5000 flight hours
and 20-year service life can be simulated in approximately 1 year. The
highest temperature of thermal spiking can be determined based on the
type of fighter plane (M = 2.0, T = 110 °C; M = 2.2, T = 125 °C).
② In general, the compression, interlaminar shear and compression after
impact strength are sensitive to hygrothermal aging. Thus, the afore-
mentioned properties must be tested for all material selected during the
structural design phase.
In particular, the residual strength should be tested at the operating
temperatures.
③ Hygrothermal aging has an influence on the facing properties of sand-
wich constructions, but it has no obvious influence on other properties
providing that the facing is undamaged and sealed.
④ Ultraviolet radiation and rain erosion have no obvious effects on com-
posite properties providing that a protective coating is maintained.
⑤ Cryogenic temperatures and changes of air pressure have little effect on
the properties of composites and may be neglected to reduce the test
costs.
Table 4.45 Paint and coatings commonly used for composites
578
Number Title Sign Main Elasticity Impact Grade of Properties Dry criterion
Standard constituent /mm /Ncm adherent T/°C t/h
force
1 Epoxy H01-101H E-20 epoxy and 1 490 1 Paint film tenacity, 18–25 24 – 26
polyamide Q/6s72-80 polyamide resin good adhesiveness 50–60 6–8
varnish with composites, 110 – 1–2
liquid proof, 120
waterproof, low air
permeability
2 Epoxy H01-102H E-20 epoxy and 1 490 1 Liquid proof, 18–25 24–26
polyamide Q/6s455-85 polyamide matches with many 50–60 6–8
priming resin, strontium other surface 110– 1–2
yellow paints, high 120
interlaminar
adhesion force,
good adaptability
3 Polyurethane 13-2 Acrylic resin, 1 490 1 Good 18–25 24–26
surface paint HOI, paint, adhesiveness, 50–60 6–8
padding waterproof, 110– 1–2
hygrothermal 120
resistance,
excellent resistance
to radiation and
weatherability
4 Elastic HTY/B-80-15(1) TDI mixed 1 490 180° peel off Paint film tenacity 18–25 24–26
polyurethane polyamine
39.2 N/cm and good elastic 50–60 8–12
magnetism Prepolymer, behavior, good
paint MOCA curing wear, rain, and
agent, paint scouring
padding resistance, good
Z. Shen et al.
adhesiveness to
priming
(continued)
Table 4.45 (continued)
Number Title Sign Main Elasticity Impact Grade of Properties Dry criterion
Standard constituent /mm /Ncm adherent T/°C t/h
force
5 Epoxy Q/6sz358-83 Epoxy Good adhesiveness 18–25 36–48
polyamide polyamide to fundus varnish 50–60 6–8
sealant varnish, talcum and transition
powder priming, easily
scraped, easily
polished, low
contraction
6 Thinner X-7 Xylene, Use to dilute epoxy
X-10 positive butyl type paints
alcohol, Use to dilute
4 Composite Structure Design and Analysis
4.11.1.1 General
Composite materials are widely used in modern structures for their high perfor-
mance and reliability. However, because these structures usually operate in hostile
and variable service environments, it is difficult to predict their structural perfor-
mance. In addition, experiments show that composite structural behavior exhibits a
wide scatter as a result of the inherent uncertainties in design variables.
Design variables, known as primitive variables, include: the fiber and matrix
material properties at the constituent level; fiber and void volume ratios; ply
misalignment and ply thickness; the fabrication process; size of random structures;
boundary conditions; loadings; and the operating environment.
The full range of structural behavior cannot be computationally simulated by
traditional deterministic methods, which use a safety factor to account for uncertain
structural behavior. Thus, the true structural reliability cannot be discerned.
A probabilistic design methodology is needed to accurately determine the structural
reliability of composite structures.
For the purposes of structural reliability analysis, it is necessary to distinguish
between at least three types of uncertainty: physical uncertainty, statistical uncer-
tainty, and model uncertainty.
dF ðtÞ
f ðt Þ ¼ ð4:126Þ
dt
The reliability function RðtÞ, which is the probability that the system will still be
operational at time t is given by
Variables that need to be considered in the stiffness reduction model can be clas-
sified into three categories:
4 Composite Structure Design and Analysis 583
distribution by the usual scale ðbÞ and shape ðaÞ parameters, the threat is charac-
terized by two impact energy levels. These are the model energy level associated
with a high possibility of occurrence ðXm Þ, and the high energy level associated
with a low probability of occurrence ðXP Þ: The relationships between the energy
parameters and the Weibull scale and shape parameters can be expressed by the
following two equations.
a 1 1=a
Xm ¼ b ð4:128Þ
a
and
Fig. 4.125 Cumulative damage tolerance strength reliability R
Fig. 4.126 Relation between Stress and Cumulative Damage Tolerance Strength Reliability
(In the Chapter, 4.1–4.6 and 4.9 were translated by Jianmao Tang; 4.7 was
translated by Jiahui Xie.)
(4.1–4.6 and 4.9 were translated by Jianmao Tang, 4.7 was translated by Jiahui
Xie.)
4 Composite Structure Design and Analysis 587
References
29. Bunin BL (1985) Critical joints in large composite primary aircraft structures vol II —
Technology demonstration test report. NASA-CR-172587
30. Bunin BL, Sagui RL (1985) Critical joints in large composite primary aircraft structures vol
III—Ancillary data test results. NASA-CR-172588
31. Nelson WD, Bunin BL, Hart-Smith LJ (1983) Critical joints in large composite primary
aircraft structures. NASA-CR-3710
32. Hortin RE, McCarty JE (1987) Damage tolerance of composites. In: Engineered materials
handbook Vol.1 Composites. ASM International, Ohio, USA, pp 259–267
33. Tang XD, Shen Z, Chen PH, Gaedke M (1997). A methodology for residual strength of
damaged laminated composites. AIAA Paper 97–1220
34. Chen PH (1999). Damage tolerance analysis for composite laminates and stiffened panel. PhD
thesis of Nanjing Aeronautical and Astronautical University
35. Chen PH, Shen Z, Wang JY (2001) Damage tolerance analysis of cracked stiffened composite
panels. J Compos Mater 35(20)
36. Chen PH, Shen Z, Wang JY (2001). Impact damage tolerance analysis of stiffened composite
panels. J Compos Mater 35(20)
37. Chen PH, Shen Z, Wang JY (2001) Prediction of the strength of notched composite laminates.
Compos Sci Technol 61(10)
38. Chen PH, Shen Z, Wang JY (2002) A new method for compression after impact strength
prediction of composite laminates. J Compos Mater 36(5)
39. ASTM D 3878-07 (2004) Standard terminology for composite materials
40. ACEE Composites Project Office (1985) NASA/Aircraft industry standard specification for
graphite fiber/toughened thermoset resin composite material. NASA RP 1142
41. SACMA 2R-94 (1994) SACMA recommended test method for compression after impact of
oriented fiber-resin composites
42. Dost EF, Avery WB, Finn LB, Ilcewicz LB, Scholz DB, Wishart RE (1993) Impact damage
resistance of composite fuselage structure, Part 3. In: Paper presented at fourth ACT
conference
43. Shen Z, Zhang ZL, Wang J, Yang SC, Ye L (2004) Characterization on damage resistance
and damage tolerance behavior of composite laminates. Acta Materiae Composite Sinica 21
(5):140–145 (in Chinese)
44. Chen PH, Shen Z et al (2006) Failure mechanisms of laminated composites subjected to static
indentation. J Compos Struct 75(1–4):486–495
45. Shen Z, Yang SC, Chen PH (2008) Experimental study on the behavior and characterization
methods of composite laminates to withstand impact. Acta Materiae Composite Sinica 25
(5):125–133 (in Chinese)
46. Li Y (1998). Effect of hygrothermal spectrum aging on composite laminates. NCCM-10.
Hunan Science and Technology Press, Changsha, pp 143–146
47. Kan HP, Whitehead RS, Kauts E (1992). Damage tolerance certification methodology for
composite structures, N92-32579
48. Bai G (1997). Reliability evaluation of composite laminates under impact threat. MSC thesis
of Northwestern Polytechnical University
49. Ma ZK, Yang L (1998) Analysis method of cumulative damage tolerance reliability for
composite laminated structures. Design and Research of Commercial Aircraft
50. Chen PH, Shen Z (2004). Statistical analysis and reliability evaluation of post- impact
compressive residual strength. Acta Aeronautica Et Astronautica Sinica 25(6)
51. Tong MB, Chen P H, Shen Z (2004). Reliability analysis of impact damage tolerance for
composite materials. Acta Materiae Composite Sinica, 21(6)
52. JSSG-2006 (1998) Joint service specification guide—Aircraft structures