Aero Struct
Aero Struct
Aero Struct
Aerospace Structures-
an Introduction to
Fundamental Problems
Purdue University
This is an introduction to aerospace structures. At the end of one semester, we
will understand what we mean by the “structures job” and know the basic
principles and technologies that are at the heart of aerospace structures design
and analysis. This knowledge includes basic structural theories, how to choose
materials and how to make fundamental design trades, make weight estimates
and provide information for decisions involved in successful aerospace structural
design and development.
While the wisdom of doing this has been hotly Similar experiences interacting with other
debated, I believe that it is a wise academic engineers, both in the working world and the
choice. On the other hand to provide knowledge research world has convinced me that there is no
about structures, over and above that found in a substitute for a blackboard/whiteboard, good
standard Strength of Materials course, with a sketches of configurations and Free Body
one-semester structures course is challenging. Diagrams. There is also no substitute for basic
Creating this course is a design problem; each mechanics; you will find a review of mechanics
instructor will have a different set of prioritized integrated into the truss analysis review in
requirements. The results will differ depending Chapter 2.
on the background of the instructor.
External and internal load computations are
The basic question that has to be answered in the extraordinarily important. “Loads” is a broad
beginning is this: ” if I want my students to have area that could easily occupy an entire semester
a skill set, a vocabulary and a level of literacy or two. One loads subject often overlooked in
that will enable them to be competitive, the educational system is inertia loading. Inertia
productive product team members for aerospace loads have destroyed more than one vehicle.
employers, what must I absolutely include in This subject is considered in Chapter 3.
this class?” My answers are contained in this set
of notes. Knowledge of the problems and issues related to
materials and material processing is important.
My experience during the past four decades has An aerospace engineer must have knowledge of
biased my vision. I have taught this course or materials, their properties, how they are used,
similar courses at The University of Maryland, and how they are manufactured. Few aerospace
Virginia Tech and at Purdue over a span of 38 curricula include such courses. Chapter 5 is a
years. During this time I also spent time on summary of considerations and choices that
i
must be made when selecting aerospace majoring in structures is free to take any or all of
materials. Chapter 5 draws very heavily on the these classes. But for the student who does not
educational approach used by Michael Ashby feel a calling to the structures and structural
and his colleagues in Great Britain. At Purdue mechanics world, this class will enable them to
we use his CES Edupack computer aided design be conversant and aware of the problems we
software. face in developing new structures and sustaining
old structures. If this goal is reached then all of
Stress analysis is important. The most important the work that it takes to write up this set of notes
part of stress analysis is the development of (and wrestle with Microsoft Word) is worth it.
simplified models that can be used independent
of elaborate computer models. Beam analysis is
briefly reviewed in Chapter 6. Thin-wall beam Terry A. Weisshaar
approximations are excellent ways to describe
these models and are included in Chapter 7. West Lafayette, Indiana
ii
Table of Contents
CHAPTER 1 – AN OVERVIEW WITH A PLAN ............................................................................................................ 1
1.1 INTRODUCTION........................................................................................................................................................1
1.2 WHAT IS AN “AEROSPACE STRUCTURE?” SOME ESSENTIAL FEATURES OF AIRCRAFT AND SPACECRAFT SYSTEMS. ..........................2
1.2.1 Spacecraft structures ...................................................................................................................................3
Spacecraft structures summary ............................................................................................................................6
1.2.2 Aircraft structures - Background .................................................................................................................7
The legacy – the first two decades of flight ........................................................................................................13
A quick trip to the present day and a look to the future .....................................................................................15
The past is prologue – new materials for the future ...........................................................................................19
1.3 THE STRUCTURES GROUP – WHAT DO THEY DO AND HOW DO THEY DO IT?........................................................................21
1.4 SUMMARY AND COURSE GOALS ...............................................................................................................................24
1.5 BIBLIOGRAPHY – CHAPTER 1 ....................................................................................................................................24
1.6 - PROBLEMS TO CONSIDER ........................................................................................................................................26
CHAPTER 2 – A REVIEW OF STATICS AND MECHANICS OF SIMPLE STRUCTURES .................................................. 27
2.1 INTRODUCTION – MECHANICS, EQUILIBRIUM AND FREE BODY DIAGRAMS ........................................................................27
2.2 APPLIED, INTERNAL AND REACTIVE FORCES .................................................................................................................29
2.3 EQUIVALENT FORCES AND MOMENTS .........................................................................................................................30
2.4 DRAWING FREE BODY DIAGRAMS .............................................................................................................................31
2.5 TRUSS ANALYSIS REVIEW - THE METHOD OF JOINTS.......................................................................................................32
2.6 STATICALLY DETERMINATE AND INDETERMINATE STRUCTURES..........................................................................................33
2.7 THE METHOD OF SECTIONS ......................................................................................................................................34
2.8 A WORKED TRUSS ANALYSIS EXAMPLE ........................................................................................................................35
2.9 A SECOND EXAMPLE-BAR ELEMENTS CONNECTED TO A BEAM...........................................................................................36
2.10 DESIGN TERMINOLOGY AND WEIGHT CALCULATION .....................................................................................................38
2.11 ESTIMATING STRUCTURAL WEIGHT ..........................................................................................................................40
2.12 ENERGY STORED IN THE TRUSS STRUCTURE – STRAIN ENERGY AND DEFLECTION .................................................................43
2.13 SUMMARY ..........................................................................................................................................................44
2.14 – HISTORICAL NOTES .............................................................................................................................................44
2.15 – REFERENCES .....................................................................................................................................................45
2.16 HOMEWORK PROBLEMS FOR CHAPTER TWO .............................................................................................................46
CHAPTER 3 – INERTIA LOADS ............................................................................................................................... 49
3.1 INTRODUCTION ......................................................................................................................................................49
3.2 CRITICAL LOADS......................................................................................................................................................49
3.3 D’ALEMBERT’S PRINCIPLE – INERTIA LOADS INSIDE THE VEHICLE ......................................................................................50
3.4 AN EXAMPLE – CARRIER LAUNCH ...............................................................................................................................51
3.5 LOAD FACTORS AND “G-LOADING” .............................................................................................................................52
3.6 SUMMARY – THE PROCEDURE FOR COMPUTING APPARENT WEIGHTS USING LOAD FACTORS ..................................................53
3.7 EXAMPLE – AIRCRAFT TOUCHDOWN ..........................................................................................................................54
3.8 CALCULATING INTERNAL REACTION FORCES AND MOMENTS ............................................................................................55
3.9 ADDING IN THE ANGULAR ACCELERATION “FORCES” ......................................................................................................56
3.10 A SECOND EXAMPLE - MISSILE JOINT BENDING MOMENTS IN A MISSILE WITH ROTATIONAL ACCELERATION ..............................59
3.11 EXAMPLE – ANOTHER IDEALIZED MISSILE IN FLIGHT .....................................................................................................60
3.12 DISTRIBUTED INERTIA FORCES–A TOPPLING SMOKESTACK .............................................................................................62
3.13 SUMMARY ..........................................................................................................................................................65
3.14 CHAPTER 3 - HOMEWORK PROBLEMS ......................................................................................................................66
CHAPTER 4 - MATRIX METHODS .......................................................................................................................... 70
4.1 REVIEW - STATICALLY INDETERMINATE STRUCTURES ......................................................................................................70
iii
4.2 THERMAL LOADS....................................................................................................................................................71
4.3 INDETERMINATE STRUCTURES-A SUMMARY ..................................................................................................................72
4.4 FINITE ELEMENT METHODS.......................................................................................................................................73
4.5 MATRIX METHODS AND THE FINITE ELEMENT MODEL METHOD ........................................................................................73
Deriving the element stiffness matrix .................................................................................................................74
4.6 ASSEMBLY OF MATRIX EQUATIONS.............................................................................................................................77
4.7 BOUNDARY AND SUPPORT CONDITIONS ......................................................................................................................78
4.8 ASSEMBLING TRUSSES .............................................................................................................................................79
4.9 MATRIX OPERATIONS AND MAGIC TRICKS ....................................................................................................................82
4.10 THERMAL LOADS RE-VISITED ...................................................................................................................................83
4.11 TRICKS OF THE TRADE-MATRIX PARTITIONING ............................................................................................................89
4.12 MATRIX METHODS IN GENERAL ...............................................................................................................................92
4.13 CHAPTER 4 - HOMEWORK PROBLEMS .....................................................................................................................94
CHAPTER 5 -–MATERIALS AND PROCESSES .......................................................................................................... 98
5.1 INTRODUCTION .....................................................................................................................................................98
5.2 MATERIALS FAMILIES AND CLASSIFICATION ................................................................................................................102
5.3 MATERIAL PROPERTIES..........................................................................................................................................104
5.4 MATERIALS METRICS AND BUBBLE CHARTS ...............................................................................................................105
5.5 MATERIALS INDICES FOR STIFFNESS DESIGN ...............................................................................................................106
5.6 SHAPE SELECTION AND OPTIMIZATION ......................................................................................................................107
5.7 THE RELATIONSHIP BETWEEN STRESS AND STRAIN .......................................................................................................109
5.8 FATIGUE.............................................................................................................................................................113
5.9 SPECIAL PROBLEMS IN AEROSPACE STRUCTURES ..........................................................................................................116
5.10 HOMEWORK – CHAPTER 5 ..................................................................................................................................116
CHAPTER 6 - A REVIEW OF BEAMOLOGY-BENDING AND SHEAR STRESSES IN CLASSICAL BEAMS ....................... 120
6.1 INTRODUCTION – THE FLEXURE FORMULA ..........................................................................................................120
6.6.1 – Moment of inertia example...................................................................................................................122
6.1.2 – Moment of inertia computation – Example 2 .......................................................................................123
6.2 BEAM SHEAR STRESS ............................................................................................................................................126
6.3 – SUMMARY - A FORMULA FOR SIMPLE SECTIONS ........................................................................................................130
6.3 PROBLEMS .........................................................................................................................................................132
6.4 REVIEW: AREA MOMENTS OF INERTIA AND THE PARALLEL AXIS THEOREM.........................................................................134
6.4.1 Example 6-1-moments of inertia .............................................................................................................135
6.4.2 Example 6-2, calculation of 2nd area moments by subsection analysis .................................................137
6.4.3 Example 6-3 ...........................................................................................................................................138
6.4.4 PROBLEMS ..............................................................................................................................................139
CHAPTER 7 REINFORCED, THIN-WALLED WING & FUSELAGE STRUCTURE IDEALIZATIONS ................................. 140
7.1 - INTRODUCTION ..................................................................................................................................................140
7.1.1 Semi-monocoque aircraft structural design ............................................................................................140
7.1.2 Wings, tails and other aerodynamic surfaces: .....................................................................................142
7.1.3 Fuselage ...............................................................................................................................................142
7.1.4 Skins .....................................................................................................................................................143
7.1.5 Spar webs .............................................................................................................................................143
7.1.6 Longitudinal elements-stringers/stiffeners/longerons ........................................................................143
7.1.7 Spar caps ..............................................................................................................................................143
7.1.8 Transverse members (ribs, frames, bulkheads) ...................................................................................143
7.2 THE SKIN-STRINGER IDEALIZATION - FUNDAMENTAL ASSUMPTIONS ................................................................................144
7.3 THE ONE-HALF AND ONE-SIXTH RULES FOR SKIN-STRINGER IDEALIZATION ........................................................................145
7.4 EXAMPLE 7-1 .....................................................................................................................................................147
7.5 EXAMPLE 7-2 .....................................................................................................................................................148
iv
7.6 HOMEWORK PROBLEMS .......................................................................................................................................151
7.7 SHEAR FLOWS IN THIN WEBS ................................................................................................................................152
7.7.1 –Homework problems ..............................................................................................................................158
7.8 COMPUTATION OF THE SHEAR CENTER FOR THIN-WALL OPEN SECTION BEAMS...............................................................159
7.9 THIN-WALL CLOSED BOX BEAM ANALYSIS ..................................................................................................................161
7.10 THIN-WALL BEAM DEFLECTION AND TWIST ...............................................................................................................163
7.11 – THE SHEAR CENTER...........................................................................................................................................166
7.12 ESTIMATING THE EI AND GJ FOR A WING .................................................................................................................168
7.13 HOMEWORK PROBLEMS .......................................................................................................................................168
CHAPTER 8 – STABILITY – UNDER CONSTRUCTION ............................................................................................. 172
APPENDIX – WORKED HOMEWORK PROBLEMS ................................................................................................. 180
Chapter 2 - Trusses ............................................................................................................................................180
Chapter 3 – Inertia forces and load factors .......................................................................................................181
Chapter 4 – Matrix methods .............................................................................................................................186
Chapter 5 -Materials .........................................................................................................................................210
Chapter 6 –Beam Bending Review ....................................................................................................................220
Chapter 7 – Thin-wall shear flow ......................................................................................................................223
More Chapter 2 solutions ..................................................................................................................................244
v
Introduction
Chapter 1 Page 1
environmental loads to support the performance surfaces that may be difficult to design if the
and mission of the aerospace vehicle of which travel is extreme.
they are part; 2) are reliable during all mission
phases; and, 3) can be maintained and repaired. In addition, since aircraft cost is a function of
We also add that aerospace structural weight and complexity, airplane cost almost
components should be as light-weight as always goes up when weight increases. Not
possible and cost effective to construct, only is the total weight important, but the weight
assemble and maintain. Aerospace structural distribution is important since it controls sizing
product development requires analysis and test, for control surfaces as the result of its effects on
together with in-service monitoring. c.g. positions and range.
Let‟s consider what aerospace structures do and All structures, from bicycles to buildings, are in
how they perform their missions. some way cost effective, but nowhere else is
such a premium placed upon weight as it is in
1.2 What is an “aerospace structure?” Some aerospace design. Aerospace structures
essential features of aircraft and spacecraft development not only emphasizes design of light
systems. weight components, but it has to consider a wide
If someone were to ask you “what makes variety of special loads, such as inertia and
aerospace structures different from other vibratory loading during launch, take-off or in-
structures?” or “what is unique about an flight maneuvering. If you drop your computer
aerospace structure?” what would your answers it is ruined. If you drop an airplane on to the
be? runway it taxis to the terminal.
Chapter 1 Page 2
Introduction
Our most common aircraft construction material Each of these component elements are building
is aluminum. This is hardly an exotic material blocks that have evolved over time to create the
since we use far more aluminum in soft-drink wonderful lightweight efficient structures we
cans each year that we do in airplanes. have today. Some of these elements trace their
However, I have never seen titanium, beryllium origin to antiquity, others have been used by
or magnesium in a beer can, but they are there in nature for millions of years. Still others are of
aircraft and spacecraft. recent origin, a response by innovative designers
to the demands of flight in the atmosphere and
We also see more “multi-functional” structures Space.
in aerospace design. The skin that covers
airplane wings furnishes wing strength and 1.2.1 Spacecraft structures
stiffness but also forms the aerodynamic shape The Space Age began October 4, 1957 with the
required to generate lift forces. This dual launch of Sputnik I. This satellite operated until
functionality is not seen in brick buildings, but it January 3, 1958. It consisted of an aluminum
is present in automobile construction which uses sphere with a 58 centimeter diameter, two radio
unitized body construction to shape low drag transmitters, four antenna and total weight of
cars. Combining structural design with shaped 83.6 kg. Contrast this with the Chandra system,
components is an essential task for an aerospace launched in 1999.
structural sub-system designer.
Today, space systems come in all sizes and
Structural failure can be catastrophic, shapes; they operate in orbits near the Earth, in
spectacular and costly, whether it is a building geosynchronous orbit and around other planets
collapse or an airplane crash. What makes the in our solar system. Some systems probe the
reliability and safety of aerospace products more outer reaches of the solar system itself. Some
challenging is the variety of operational loading large systems, like the International Space
conditions. The loads on a skyscraper, even if Station, are products of multiple launches.
we throw in resisting loads from high winds and
earthquakes, are nowhere near those experienced
by a spacecraft whose journey to space is
marked by extreme vibrational and acoustic
loading and then followed by environmental
conditions due to radiation and extreme cold.
Chapter 1 Page 3
vacuum; radiation, solar UV, and atomic used to impart energy to cause movement of a
oxygen; temperature extremes and thermal device or subsystem for deployment as well as
cycles. It also includes unique and diverse providing continuous motion and positioning.
loading such as: pre-launch handling, Mechanisms also provide means to lock a device
transportation, test, and deployment; launch with or subsystem.
axial and lateral accelerations, acoustic, staging
separation, shock, and vibration; and, on-orbit The main or “primary” bus structure transmits
vibration, thermal management and loads to and from the base of the satellite during
maneuvering. launch through structural attachments and
components such as tubes, honeycomb plates
and panels and trusses.
Chapter 1 Page 4
Introduction
Chapter 1 Page 5
“bubble chart” to help the material selection Titanium replaces aluminum in higher-
process. temperature environments; it has the ability to
remain strong and stiff at temperatures up to
1,200º F. Unfortunately, titanium becomes
brittle at low temperatures and does not have
good fatigue properties. It is also more difficult
to weld.
Aluminum, magnesium, titanium, and beryllium Galileo Spacecraft, launched to Jupiter from a
are the elements that make up the major Space Shuttle in October 1989. This spacecraft had
lightweight alloys used in space vehicles. They a probe that separated from the main craft and
are all much lighter than steel and are penetrated the Jovian atmosphere.
nonmagnetic. Aluminum alloys are the most
widely used structural materials. Spacecraft structures summary
From this brief introduction we can conclude
that there are several special features of
spacecraft structures. First of all, there are few
requirements for the structural form of the
spacecraft structure. Its efficiency is determined
by the ease of attaching other elements of the
bus to it. On the other hand, there is a need to
design folding or latched structures that deploy
on station. This structure/mechanism design can
be challenging.
Chapter 1 Page 6
Introduction
required by the payload and by the Early airplane builders saw structural design and
communication systems. This leads to a large construction as a minor problem, easily solved
number of beam-like or shell-like deployable using commonly available materials and trial-
structures. Deployable structures experience and-error procedures. Of the conventional
high shock loads; for example during latch-up of materials available for construction of bridges
hinges in their final positions. In addition, and other stationary structures, only steel and a
lightweight structures with relatively large few varieties of wood were considered suitable
surface areas (such as solar panels and for flight structures. Familiarity with
spacecraft antennas) are sensitive to acoustic woodworking and wood construction made
loads at launch. Thermal deformation and wood the natural choice of aviation pioneers.
stresses due to temperature gradients in the Wood was readily available, inexpensive,
structure must be calculated to check alignment required few specialized tools and could be
requirements. easily worked by those with limited construction
skills.
An important output of the structural design and
bus integration is sub-system mass. The bus or Early aviation construction could be called the
“platform” mass is the total, on-orbit dry mass era of wood and piano wire since joining
(mass without propellant) of the spacecraft sub- wooden trusses and other components together
systems, structure, thermal control, power, to form a truss-like structural framework was the
attitude control, telecommunications and data easiest solution to early aircraft structural
sub-systems. problems. After all, if this type of construction
worked for bridges, why not airplanes?
1.2.2 Aircraft structures - Background
Aircraft design is driven by the goal to transport Octave Chanute (1832-1910) is the father of
cargo and passengers quickly from point to airplane structures. Chanute was a successful
point. Since aircraft are highly integrated civil engineer who, after his retirement from a
systems, developments in propulsion, successful career as a railway engineer, became
aerodynamics and structures feed on each other interested in flight. An additional item of
to produce aircraft that operate from relatively interest for Indiana residents is the fact that his
low speeds to hypersonic speeds. glider experiments were conducted at the
Indiana Dunes in Gary, Indiana.
This section reviews some of the historical
developments in aircraft structural design. Its
purpose is to acquaint you with the rich history
of aeronautical vehicles. This section will set
the stage for the material in the chapters to
follow.
Chapter 1 Page 7
winning truss railway bridges – to create the first the Wrights built a wind tunnel to create their
biplane surfaces. This type of aerodynamic and own data set.
structural design was adopted quickly by the
Wright Brothers. The Wright Brothers were relatively secretive
about their efforts because they feared ridicule.
Chanute‟s glider experiments were not without Typical of this ridicule was a comment on an
controversy. On September 11, 1896, an unsuccessful Langley effort, published in a New
unnamed Chicago newspaper wrote: “Flying- York Times editorial published on October 9,
machines are being tested in Indiana, it being an 1903, only a few months before the Wright‟s
axiom among confidence men that a scheme that first flight.
will go anywhere will go in Indiana.” On the
other hand, the Chicago Times Herald wrote on The source of the Times‟ wrath was the
September 12, 1897: “Those who have by unsuccessful flight of Langley‟s airplane, an
casual reading formed the opinion that a few airplane project with substantial U.S.
cranks had hidden themselves among the sand Government backing. When Langley‟s airplane
hills at the foot of Lake Michigan for fear failed at take-off from a houseboat on the
somebody might steal the pattern of their airship Potomac River near Washington, the
have made a wide mistake...” condemnation from members of Congress and
the New York Times was loud and furious.
Truss construction, fabricated from simple “two-
force” members that could be easily replaced
when damaged, creates a stiff, lightweight
structure. Reliable methods for analysis and
design of beams and trusses had been available
for decades.2 Truss structures were ideal for bi-
plane and tri-plane aerodynamic designs.
The Wright Brothers were the first true Don’t believe everything you read in the newspapers
aeronautical engineers. They began their quest or see on the Internet. There is a big difference
for flight by mailing a letter to another aviation between people who know how to do things and
pioneer, Dr. Samuel Langley, on Memorial Day, people who don’t.
1899. Their first flight, only a little over three
and one-half years later in 1903, launched the Congressman James Robinson (from Indiana) is
age of aviation. quoted on January 24, 1904 (a month after the
first flight) “Here is $100,000 of the people’s
The Wright Brothers built on the success of money wasted on this scientific aerial navigation
others such as Chanute. Using other people‟s experiment because some man, perchance a
data was not without its own difficulties. Some professor wandering in his dreams was able to
of Chanute‟s aerodynamic data was incorrect so impress the [military] officers that his aerial
scheme had some utility.”3 (Robinson (1861-
2 3
S.P. Timoshenko, History of Strength of Materials, Roger D. Launius, “Innovation and the
McGraw-Hill Book Company, New York, 1953. Development of Flight,” 1999.
Chapter 1 Page 8
Introduction
1942) was from Ft. Wayne, Indiana. He was an The Wright‟s structural engineering tasks
unsuccessful candidate for reelection in 1904, included flight safety. Two entries in Wilbur
but probably not because of his failure to Wright‟s diary are typical: "I am constructing
embrace aeronautics.) my machine to sustain about five times my
weight and am testing every piece." (1900) and
“We also hung it by the tips and ran the engine
screws with the man also on board. The strength
of the machine seems 'OK'." (December 2,
1903)
5
A. W. Judge, The Design of Aeroplanes, Isaac
4
Richard Hallion, Taking Flight: inventing the Pitman and Sons, Ltd., 1917, p. 156. This 242
aerial age from antiquity through the First World page textbook by Judge is one of the early
War, Oxford University Press, 2003, page 156. textbooks on airplane design.
Chapter 1 Page 9
Because of this haphazard approach to design, and the landing gear below. Bleriot‟s designs
early aircraft had a high incidence of accidents added an additional post, half-way along the
caused by structural failures. As early as 1911 wing span. This gave his designs increased
structural failure was identified as the principal strength and permitted of a greater wing span
cause of flight accidents. and thus reduced drag.
Chapter 1 Page 10
Introduction
Loftin8 notes “Aircraft design during World The wood layers were glued together with their
War I was more inventive, intuitive and daring grains running in different directions to
than anything else… The principles of strengthen the skin. Ruchonnet called the new
aerodynamics that form so important a part of technique monocoque construction, or single-
aircraft design today were relatively little shell construction.
understood by aircraft designers during the
War… In the area of engineering in which
structural strength, lightweight and
aerodynamic efficiency are so important, it is
indeed surprising that a number of relatively
good aircraft were produced.”
Early airplane structures used trusses and frames Bechereau used the monocoque fuselage,
to resist bending and shear forces. This was together with a monoplane design the
particularly effective in creating lightweight Deperdussin streamlined, externally braced,
fuselage structures. mid-wing monoplane. This airplane won the
Gordon Bennett Cup in 1913 with a speed of
In the beginning aircraft structures were little 125 mile per hour – only a decade after the
more than cloth covered trusses and frames. Wright‟s first flight!
However, in 1911, Louis Bechereau built on an
idea from Eugene Ruchonnet, a Swiss enginner
who had worked at as a shop foreman for the
Antoinette Airplane Company in France. Rather
than simply mounting fabric on to a truss,
Ruchonnet‟s idea was to use the airplane‟s skin
carry the structural load and eliminate the heavy
truss. He formed the fuselage out of multiple
layers of wood to create a streamlined shape.
Deperdussin Racer design with monocoque shell
structure fuselage, 1912. Compare this design to
8
L.K. Loftin, Quest for Performance: The the Wright Brothers airplane only nine years
Evolution of Modern Aircraft, NASA SP-468, earlier.
Washington, D.C., 1985, pp. 8-9. Later, when airplane size grew the monocoque
9
J. D. Anderson, A History of Aerodynamics, design was not suitable because it buckled under
Cambridge University Press, 1997. high loading. A variation of the concept, a
Chapter 1 Page 11
reinforced or semi-monocoque structure was streamlined flying wing, and then developed an
created to create a light-weight durable structure. internally braced cantilever monoplane design.
This type structure is still used today, both for Junkers‟ monoplane used combinations of steel
fuselage design and wing design. beams interconnected by load-bearing skin to
create a semi-monocoque wing structure.
A box spar was created by constructing the main
supporting member of the wing as a long,
narrow box beam that provided the strength and
also torsional stiffness. While cantilevered
wings first covered with fabric in the same way
as biplane wing structures, they began to use
skins made of thin wooden veneers. These
wooden-skinned wings became known as
stressed-skin wings because the veneers added
strength while they shaped the wing surface Interior of the Zeppelin Rigid Airship, the
itself. Bodensee, showing how ring frames and
longitudinal girders form hull framework
Three men contributed to the development of
revolutionary aircraft structures and airplanes Junkers primary interest was in using metal
with these structures before and after World War rather than wood for construction of low-drag
I. All three were German. monoplanes with high durability. Steel
structures were too heavy so he used aluminum
The first two pioneers, Claude Dornier and alloys developed for Zeppelin airships. His first
Adolf Rohrbach, began their careers as structural all-metal, aluminum, cantilever wing airplane
engineers for the Graf Zeppelin Company. Both was the Junkers J.1.
had extensive experience designing the internal
metal framework for the large airships
developed before World War I. Rohrbach‟s
academic training was in shipbuilding. Both
eventually went on to found their own aircraft
companies.
Chapter 1 Page 12
Introduction
Dornier achieved early fame by designing The legacy – the first two decades of
seaplanes. His early efforts were done outside flight
of Germany because the Versailles Treaty The past few pages have discussed aircraft and
ending World War I prohibited production of aircraft structures from the beginning with the
aircraft within Germany. Dornier‟s Wal Wright Brothers to after World War I. We have
(Whale) aircraft was used by Norwegian polar seen that the first airplane builders used tried
explorer Roald Amundsen in his unsuccessful and true methods and materials for aircraft
attempt to reach the North Pole in 1925. structures. The most common structural form
was a truss composed of wooden members and
wires. Very little analysis was used, but
structural testing was the rule of the day. From
the beginning, emphasis was placed on light
weight construction.
Chapter 1 Page 13
come in the first 15 years of flight. Of the three, In England, government sponsorship led to the
the German Fokker D VII is most interesting founding of the Royal Aeronautical
from a structures standpoint. This aircraft used Establishment at Farnborough, formerly the
fully cantilevered wings, but disguised their Royal Balloon Factory. This effort produced
structure, not to fool the enemy, but to appease superior British fighter airplanes for World War
the German Army buyers. All of these aircraft I, but also produced technical data that replaced
were biplanes and all were wood and wire intuition with scientific fact.
construction.
In Germany, academic leadership was furnished
by Ludwig Prantdl, who built large laboratories
in Göttingen. This not only allowed
development of new theories but also served as a
clearinghouse for information exchange.
Chapter 1 Page 14
Introduction
In 1916 the first textbook dedicated to aircraft A list of useful historical references is provided
structural design appeared (see footnote 6). This in Section 1.5. These articles and books will
text focused on beam and truss analysis as well give an interested reader an historical
as engineering mechanics. perspective to the social and technical events
that influenced aircraft design in general and
Today, texts still stress beam and truss analysis - structural design in particular from 1920 to the
the size and speed of the airplane has changed present day.
but the fundamental construction layout has not.
On the other hand, aircraft today make use of On the other hand, a few points need to be made
reinforced shell, semi-monocoque construction to understand how we got to this point in the
so the so-called skin-stringer analysis (Chapter structural design and development business.
7) has been emphasized during the past 70 years. The 1920‟s were tumultuous in that biplanes and
In addition, the use of matrix methods (Chapter monoplanes continued to be developed. Both
4) for solution of stresses and deflections has wood and metal competed equally as the
largely displaced classical methods used 60 material of choice. The recognition of the
years ago. importance of streamlining began to take hold,
as did the development of new, more powerful
During the first twenty years of aviation aircraft piston engines.
engines improved and became more powerful,
enabling higher speeds. Monoplane design, In America, the belated interest by the U.S.
aided by the semi-monocoque structural government led to subsidies for airlines. The
concept, contributed to the success of low drag greater travel distances and the promise of
monoplane. The concepts largely replaced shorter travel times led to a birth of passenger
biplanes so that by 1935 new designs were aircraft constructed by American designers such
solely monoplane designs. In addition, metallic as Glenn Curtiss, Donald Douglas and Jack
construction replaced wood and cloth. Northrop.
However, beams and trusses remained the basic
internal structural components, augmented by The tension between biplanes and monoplanes is
semi-monocoque structures. exemplified by the Curtiss Condor, first flown in
1933.
These components required analysis. Some of
the analysis was already well-known and present
in civil engineering texts and papers. As
configurations became more complex analysis
had to respond to new structural forms and new
regulations mandated by government agencies.
In general these new requirements were in
response to increased airspeed and airplane
accidents.
A quick trip to the present day and a look Curtiss AT-32 Condor – 1933 -this airplane had
modern features such as a streamlined, semi-
to the future
monocoque fuselage and a retractable landing gear,
By 1920 all of the essential structural concepts but only held about 12 passengers. It was the last
were in place so that structural design could biplane passenger airplane produced in the U.S.
enable better, faster aircraft. In particular, semi-
monocoque structural design had enabled The Condor‟s streamlined fuselage looks out of
streamlined fuselages and low-drag cantilever place with the strut-braced biplane wings. Only
wings. Metallic designs created strong, robust, about 40 were built but they were in service with
durable structures. Eastern Airlines and American Airlines.
Chapter 1 Page 15
The commercial transport biplane‟s demise was was banning wooden aircraft from U.S.
hastened by two factors. The first was passenger service.
innovative monoplane designs such as the
Boeing 247 which also appeared in 1933. This
airplane used semi-monocoque metal
construction, together with features such as
retractable landing gear.
Chapter 1 Page 16
Introduction
The advent of the economical airliner in the engine with swept wings to produce lethal
United States, led by the DC-3, transformed fighter/interceptors, fortunately too late to be a
aviation and made the United States a leader. factor in the European air war.
This could not have happened without enablers
like metallic semi-monocoque structural After the war, a team of American engineers was
concepts. As we have seen these concepts were sent to Germany to collect aeronautical data.
not invented in the United States but were Some of this data related to German swept wing
substantially advanced by American designers. research and development. This data was
instrumental in creating the Boeing B-47 jet
The NACA played a significant role in airplane bomber.
structures development in the 1920‟s and
1930‟s. They developed new structural analysis
methods and provided valuable test results to
American industry. As a result, American
propeller driven commercial transports were the
best in the world.
Chapter 1 Page 17
revolutionized air travel for the last half of the
20th century. Boeing again became the premier
air transport company.
Chapter 1 Page 18
Introduction
with the variable sweep Bell X-5, the variable consisted of six electric actuators on each wing
sweep aircraft enjoyed a two decades long portion and weighed several tons.
popularity run as a strike fighter that
incorporated the best features of a low speed The past is prologue – new materials for
unswept wing with a high speed sweptback the future
wing. Aircraft such as the Navy‟s F-14 Tomcat, In Chapter 5 we consider materials and design
featured in the movie “Top Gun,” enjoyed a long choices. However, we should mention in this
service life. introduction how new materials are literally
shaping the future of aviation. These new
These aircraft were mechanized structures that materials are the so-called advanced composites.
included heavy wing pivots and motors or The picture below, taken from a website shows
hydraulics to move the wing in flight. The the essentials of composite construction
systems payoff for this weight addition was the
ability to operate efficiently over a wide speed When the Wright Brothers built their airplane
range. and for 20 years afterward, wood was the choice
of airplane builders for a number of reasons
In 1958 the U.S. Air Force began design studies previously cited. Wood is a blend of cellulose
for a supersonic bomber. The result was North fibers embedded in a matrix binder. As such it
American Aviation XB-70 supersonic bomber. is very strong for its weight. Plywood, an
artificial wood is composed of many different
layers glued together and was used for the first
monocoque fuselages and as wing covers in
stressed skin construction.
Chapter 1 Page 19
matrix that held them together as a unit. These this difference is the advanced composite
matrices were polymers, metals or resins. . fuselage of the Beech Premier Business Jet.
This fuselage is a three piece, filament wound
The individual fibers are very strong due to the component. An outer and inner layer of
fact that there are few imperfections in their composite material forms a sandwich structure
structure. They are also made of materials such that resists internal pressure and other fuselage
as carbon that has a very low molecular weight loads. The wing and tail use multi-spar metallic
compared to aluminum. structures.
Not only are composites strong and stiff That reticence on the part of commercial
compared to metals, but the manufacturing transport manufacturers is changing rapidly and
processes are different. An excellent example of
Chapter 1 Page 20
Introduction
Chapter 1 Page 21
mold lines”) right. For a wing, the structural There are thousands of such “load sets.”
effort includes the choice of the number of spars Understanding the origins of these loads, the
and ribs and the location of critical components. vocabulary associated with these loads - terms
The activities in Block 2 are characterized by such as “ultimate load,” “limit load” and “load
low level structural analysis using simple beam factor” - and how they are computed is one
models or other simple idealizations for element purpose of AAE352.
sizing and response estimation. Seldom will the
efforts in Block 2 revisit decisions made in
Block 1, particularly those involving planform Design loads
Environment and
shapes. discrete events Static strength
Damage tolerance
and safe life
Block 4, the certification activity, requires
extensive stress analysis. Stress analysis is labor Figure 1-3 - Structural design requirements
intensive and involves large, analytical models
developed from the drawings generated in Block To illustrate the diverse sources of these loads,
3. Next, Block 5 includes full scale and consider just a few from this list:10
component tests including fatigue and static
loads tests. Block 6 is the beginning of Air loads: lift, drag, gusts
production and delivery to the customer.
Acoustic loads
Thermal loads
Note that Blocks 2 through 6 involve the
opportunity and need to re-visit decisions made Landing loads: touchdown, arresting
in earlier Blocks. Analytical modeling furnishes Take-off loads: runway taxi, catapult,
additional information that allows discovery and assisted (e.g. jet or rocket) take-off
knowledge generation as the effort progresses. Power plant loads: thrust, engine torque,
On the other hand, you cannot analyze what you gyroscopic effects
have not defined, so surprises may visit the Special loads: towing, refueling,
development team often. pressurization, cargo, weapons recoil.
There is another way of looking at the structural Continuing clockwise to the right on the wheel
design activity. The structural design process in Figure 1-3, we have the requirement for static
begins with very general customer requirements strength. Static strength or “Not breaking” is a
that lead to clearly stated engineering design fundamental requirement. Static strength
criteria, complete with numbers or “metrics” that includes the requirement that the structure must
must be satisfied. These metrics include failure support: 1) structural limit loads, the largest
stresses and maximum allowable deflections. A loads likely to be encountered in practice during
summary of these general design criteria is the lifetime of the structure, without permanent,
shown in Figure 1-3. detrimental deformation; and, 2) the ultimate
Chapter 1 Page 22
Introduction
loads, defined as the limit loads times a factor of airplane flies, the more significant are
safety, for at least three seconds without aeroelastic interactions.
catastrophic failure or instability.
Chapter 1 Page 23
Crashworthiness is next on the wheel. Crash review statics, basic beam theory and truss
worthy designs protect occupants by absorbing analysis, as well as review structural failure
large amounts of energy to cushion passengers. criteria like buckling. To these topics we will
Protecting fuel tanks from rupture is also add a study of loads and how to calculate them,
important and can be done by using breakaway knowledge of materials sufficient to make
nacelle or flap attachments or by designing the choices, and new information about advanced
keel structure to absorb vertical impact and analysis, such as finite element methods.
scraping action from a wheels-up emergency
landing. Our primary goal is to develop a perspective of
what is required for creative structural design
Next on the structures wheel are the “abilities.” and reinforce your ability to function as a
Producibility is defined as the “ease of productive aerospace product team member.
manufacturing an item (or a group of items) in You should also be able to make or support
large enough quantities.” Producibility depends decisions in which structural weight, materials
strongly on design features that enable selection, component selection, structural form
economical fabrication, assembly, and and manufacturing processes are traded against
inspection or testing. each other.
Chapter 1 Page 24
Introduction
Chapter 1 Page 25
The book demonstrates that design and 2. Miller, R. and Sawers, D., The Technical
development of modern spacecraft structures Development of Modern Aviation, Praeger
is much more than structural analysis.) Publishers, New York, London, 1970.
22. Sechler, E.E. and Dunn, L.G., Airplane 3. Trimble, W.F., Jerome C. Hunsaker and the
Structural Analysis and Design, John Wiley Rise of American Aeronautics, Smithsonian
& Sons, Inc., 1942. Institution Press, Washington, D.C., 2002.
23. Shanley, F.R., Weight-Strength Analysis of 4. Loftin, L.K., Quest for Performance: The
Aircraft Structures, Dover Publications, Evolution of Modern Aircraft, NASA, 1985.
1952. (This is a classic book by a 5. Hallion, R.P., Taking flight : inventing the
practitioner. It is a “must-have” book.) aerial age from antiquity through the First
24. Steinbacher, F.R. and Gerard, G., Aircraft World War, New York : Oxford University
Structural Mechanics, Pitman Aeronautical Press, 2003.
Publications, 1952.
25. Sun, C.T., Mechanics of Aircraft Structures, 1.6 - Problems to consider
John Wiley & Sons, 1998. 1. The FAA is in charge of airplane
26. Williams, D., An Introduction to the Theory certification. What else does the FAA
of Aircraft Structures, 1960. do? Look up information about the
27. Young, W.C., Roark‟s Formulas for Stress Federal Aviation Regulations.
and Strain, 1989. 2. How maintainable is your cell phone?
What would you do if it has a problem?
The history or aircraft development is a Does anyone repair cell phones? Why
fascinating subject. Students, young and old, or why not? What is the trade-off?
can learn a great deal by observing the flow of 3. How does the design environment for an
history and the men and women who have automobile differ if it is operated in
contributed to aviation development. These Minnesota compared to south Texas?
include engineers, pilots and industrial What additional design requirements
developers. Some contributors were eccentric might you need to allow a car, originally
while others were quiet and u all were important. designed only for operation in Texas, to
operate in Minnesota?
The list that follows represents a small portion 4. How would design of an Indy Racing
of the material available today. These Car differ from the design of an ordinary
references are some of my favorites, but are not automobile? What design requirements
intended to be the end of the story. for maintainability would you add for
the Indy Racing Car?
1. Gibbs-Smith, C.H., Aviation, An Historical 5. Is a racing bicycle structural design fail-
Survey Form Its Origins to the End of World safe? Why or why not?
War II, Her Majesty‟s Stationery Office,
London, 1970.
Chapter 1 Page 26
Mechanics of truss structures
One of the essential skills required for aerospace You have already spent considerable time
structural design is the knowledge of drawing these types of diagrams in your
engineering mechanics and how to apply introductory courses, but we need to make sure
mechanics principles to the design process. we are all on the same page so we‟ll take some
Mechanics is an applied science that deals with time to review basic theory. You may need to
forces, moments and motions of isolated systems go back to the textbooks you used for statics and
or parts of systems. The simplest case occurs dynamics or strength of materials. However, in
when there is no system motion; we simply call this section we‟ll go heavy on review.
this situation “statics.” The concept of static or
dynamic equilibrium is fundamental to We will use truss structure design as an example
mechanics; all of the forces, moments and of how to apply mechanics to structures
motions are in balance. problems. Truss analysis was one of the first
analytical mechanics methods used by airplane
Central to the application of engineering designers. While truss construction is minimal
mechanics is the ability to draw clear diagrams in modern aircraft structures, as we saw in
that summarize the loads acting on the structural Chapter 1, truss structures are used for
unit and to draw isolated pieces of the structure spacecraft designs, such as that shown in Figure
where external loads and the internal reactive 2-1.
loads are shown in balance. The term used for
these diagrams is a “free body diagram.” The
ability to create these diagrams is the first
essential step in the process of structural
analysis.1
1
The first textbook describing the design and
construction of aircraft, “Elementary Principles of
Aeroplane Design,” by Arthur W. Judge, appeared in
England in 1916. Judge was an accomplished
automotive engineer at the time. His book
emphasized graphics and mechanics. Theoretical
principles of structural analysis applicable to bi- Figure 2-1– A truss design to support a solar array
planes of the time were well-known to railway and
bridge engineers. A book published later, in 1935,
Trusses are defined as structures composed of
had the title “Mechanics of Aircraft Structures” to
also emphasize this importance. Its author, John
slender members that resist loads along their
Younger, was a former senior engineer at Wright axes. These “slender” members include wires.
Field, an accomplished practioner, and later a Trusses are used extensively in bridge design
Professor at the University of Maryland. such as that shown in Figure 2-2. Trusses have
special properties and assumptions. First of all,
Chapter 2 Page 27
Truss structures
external loads are applied only at joints or considered to be the easiest solution to early
“nodes.” In addition, the components are aircraft structural problems. This type of
assumed to be pinned together, making them construction created a stiff, lightweight structure
“two-force” members that resist loading along fabricated from simple “two-force” members
an axis that passes through the node points. that could be easily replaced when damaged.
Seldom are the joints really pinned together, but Reliable methods for analysis and design of
the internal loads determined using this beams and trusses had been available for
assumption dominate over other “secondary” decades.2
loads.
Chapter 2 Page 28
Mechanics of truss structures
Analytical modeling is
an art; as your career
progresses you will
become better at this art.
Modeling is so important
that we‟ll have more to
say about this in a later
section of the book. The
Figure 2-5 – The Wright Flyer – 1903 – flight most important thing to An interstage structure
mechanics combined with structural mechanics remember is the obvious, connecting first and
in the form of an open space structure poor models and poor second stages of the
assumptions lead to poor, Zenit-2 vehicle.
pilot/passenger. Their structural approach was Copyright © 2001 by
“cut and try” rather than “analyze, cut and try.” sometimes fatal results.
In structures the most Anatoly Zak
Like other early designers, they got lucky.
important first step of idealization is the
While Newton provided three equations or construction of the Free Body Diagram. If you
“laws” for the motion of a particle in space, can‟t do this in your sleep, you are not going to
when we consider the motion of three- be a happy engineer.
dimensional bodies such as an aircraft or a
spacecraft, there are six “balance” equations that 2.2 Applied, Internal and Reactive forces
apply to our work: three equations for The end products of structural analysis include
equilibrium of forces; and, three equations for structural stresses and deflections. These two
static equilibrium of moments. In structures we items cannot be computed unless we know the
have special terms for these moments such as internal resultant forces and moments such as
“bending moment” or “torsional moment” to tell those we have just computed for the example
us what the moment does to the structure. structures. Determining how these forces and
moments are distributed over the internal cross-
For instance, in a two-dimensional x-y system sections is the subject of stress analysis and
we have three independent equations of static “strength of materials.”
equilibrium. These are:
Determining the size and directions of loads that
will be applied to our structure is a crucial first
x
Fx 0 step towards structural design. Applied forces or
loads act directly on the structure and occur as
y
Fy 0 the result of operation of the structure in its
environment. A second set of forces are applied
to the structure, but these are called reactions or
z
Mz 0 “reactive forces.” They are created by the
actions of one component interacting with
Forces sum and balance, but forces on what? another or with a “boundary.” This interaction
There is another feature of mechanics that is may be the reaction between a structure and a
important; this is the concept of a behavioral support for instance. It may also be the result of
“model.” To keep the problem within bounds, the connection between two elements which are
we need to look at the actual physical structure pulling or pushing on each other.
in an idealized way, eliminating non-essential
Central to modeling is the construction of
equilibrium diagrams or what we call “free body
Chapter 2 Page 29
Truss structures
diagrams” or FBD‟s. A free body diagram is a moments that the surroundings, such as
sketch of the system of interest and includes all supports, impose on the system. The free body
of the forces that act on the system. These diagram shows what it would take to
diagrams identify the unknowns in the structures mechanically fool the system into remaining at
problem and help us determine the magnitude of rest if it was cut free of its support.
the unknowns. Unknowns in structures
problems include reactions at supports or The motion of the system should be totally
internal forces such as shear forces in beams, unchanged if it were cut free and the forces
bending moments and shear forces. From these shown on the free body diagram were applied as
internal forces we can then determine stresses. a replacement for all external interactions. In
particular, the applied and reaction forces and
These stresses, also the result of behavioral moments must be shown on the free body
theories that include theory of elasticity, allow diagram at the points where they are applied. For
us to estimate how the resultant forces and reactions, these places are where you made
moments are distributed throughout the structure „cuts‟ to free the body.
and what deflections are caused by the external
loads. These stress distributions and deflections Internal reactions are the result of Newton‟s 3rd
allow us to assess, and perhaps change, essential Law which states that “for every action there is
design features of the structure. These features an equal and opposite reaction.” When we look
include component thickness or component at structures we have cut to expose internal
geometry such as we showed in Chapter 1. reactions, we notice that the forces come in
We‟ll show other examples of this later. pairs, equal, but opposite in direction.
A free body diagram is essential because it In a structure there are many pairs of “action-
precisely defines the system to which you are reaction” forces. We can see these pairs as we
applying Newton‟s mechanics equations and use FBD‟s to disassemble the structure. These
identifies all of the important forces to be action-reaction pairs, when plotted or displayed
considered. The FBD shows the system isolated on a diagram, show “load paths” from the
or “freed” from its environment. We use two applied loads to the supports. An example is the
kinds of FBD‟s: 1) a FBD showing the applied path that aerodynamic forces take from the wing
loads and the reactions; 2) a series of FBD‟s in surface to the fuselage attachment or “support”
which the structure is cut into pieces to show for the wing.
internal reactions.
2.3 Equivalent forces and moments
Drawing either of these FBD‟s is not always Aerodynamic forces acting on a wing are
easy. It takes practice. The free body diagram distributed as pressures over the wing. For some
of a structural system must shows the forces and calculations, we really don‟t need a precise
description of the
pressures. For instance, if
we are interested in the
acceleration of the airplane
we just want to know the
total or “integrated” force
in certain directions. We
want to know a single
Figure 2-6 – Examples of equivalent force systems. The effect of the four
loadings on the reaction at the wall is identical. Note that if the problem is to
quantity such as a vector
find the internal loads between points A and B the four systems give different that is equivalent to the
results. For internal loads we need to ask for the real applied loads. (Source: A. pressure distribution as far
Ruina, P. Rudr, Introduction to Statics and Dynamics, Oxford University as its effect on acceleration
Press, 2002 is concerned. The action
of the single vector on the
Chapter 2 Page 30
Mechanics of truss structures
center of mass of the airplane is equivalent to the must place a reaction moment or torque at that
action of the pressure distribution. In fact, all point. In case (b) we have used one of the
that the equations of motion know about the equations of equilibrium to eliminate “by
situation is the effect of the inspection” the
integrated force and moment horizontal force.
about some point. Examples of The concept of resultant or equivalent While this is
equivalent force and moment forces is central to all structural correct, it is also
systems are shown in Figure analysis. Two force systems are dangerous because
2-6 equivalent if they have the same sum we may forget that
(the same resultant) and the same net there is a boundary
This concept is important for moment about any arbitrary point in condition at this
structural analysis since it is space. point.
often more convenient to
replace a set of forces and Many of these
pressures as a single load. There are some rules bridge design used truss construction. Truss
about how to do this and we will deal with them elements were steel, iron or even wood. Trusses
as the situation requires. have attachments between members that result
in only a small transfer of moments at these
2.4 Drawing Free Body Diagrams joints. As a result we usually refer to the
We know enough about the representation of attachments as “pinned joints.” The presence of
forces as vector quantities and how to resolve or pinned joints makes the truss elements special in
“break up” these vectors into components. This that they are “two-force” elements. In this case
is a form of “static equivalence.” This term is the “line of action” of the internal forces in each
used to mean that, from the viewpoint of motion truss member must lie along the truss member
of a rigid body, two different systems, although itself. This is in contrast to beam elements
they look different (a single vector replaced by which can transfer bending moments and thus
two vectors), produce exactly the same permit shear forces to be carried within the
translational or rotational effects. When we do element.
static analysis we do not have translation or
rotation so the feature we are most interested in
is the calculation of the reactions at supports.
Chapter 2 Page 31
Truss structures
At joint B
and
Chapter 2 Page 32
Mechanics of truss structures
and
---------------------------------------------------------------------------------------------------------------------------
Chapter 2 Page 33
Truss structures
Let‟s first consider how many equations we have then the truss is externally determinate, but
at our disposal. All trusses considered in this internally statically indeterminate.
chapter lie in a single plane and are called
“planar” trusses. We have assumed that truss In general, when r is greater than 3, then, if
element ends are joined by smooth frictionless m+r>2*j we automatically have a statically
pins that do not transmit bending moments and indeterminate structure. The degree of
that the truss elements are loaded only at the indeterminancy DoI is given by the difference
joints or nodes. We have also made sure that DoI=(m+r)-2*j. Our review and the homework
enough truss elements are connected so that the problems in this chapter will only consider
truss cannot behave like a mechanism and statically determinate truss structures.
collapse.
2.7 The Method of Sections
The example in Figure 2-8 has 7 elements, 5 What happens if I don‟t really want to know all
joints or nodes and has 3 unknown external of the truss element forces? Instead of drawing
reactions. Let‟s assume that the more general seven free body diagrams, we can cut the
truss is composed of “j” joints, connecting “m” structure into two pieces, exposing the force that
bar elements with “r” reaction forces at the truss we want. Consider the problem of finding the
boundaries. No matter how many elements force in the element BC. The Free Body
there are in a truss, there are only three Diagram for this calculation is shown in
equations of static equilibrium that are useful to Figure 2-10. This FBD is the result of
solve for external reactions. These three are the sectioning the structure into two pieces, hence
two force summation equations and one moment the name “Method of Sections.” This method is
equation. We can always use two moment also known as the Method of Moments or the
equations and one force equation, but the Method of Shears.
number of useful equations can never be more
than three.
Chapter 2 Page 34
FBC. If we want to know the other internal force 20 P
in CE we can use the vertical force equilibrium M B 0 R2 h 20 P or R2
h
equation (the so-called “shear equation.” z
h 1
the two-force element BC. Although the h
reaction forces R2 and R3 are drawn so that they
2
appear to be independent, they are not. Like R1 R 20 Ph 20
these two reactions must have a resultant that is FAB 2 1
cos 20h h
in-line with member AB, as indicated in Figure
2
2-13. 20
FAB P 1
h
20 P
FBC R1
h
Figure 2-12 – Free-body diagram for truss
structure with supports removed We‟ll return to this example later when we show
how to calculate weights and optimize the
Let‟s solve for the reaction forces. Summing structural topology.
moments about point C we have:
In general, structures composed of two force
members make the analysis job easier, but you
Chapter 2 Page 35
Truss structures
Chapter 2 Page 36
or FCF = 1000 lb.
F
The static equilibrium equations at this node are:
x 0 3000 FCD FDF cos 45o
F
x
0 2000 FAG cos 45o
F
y
y 0 1000 FDF sin 45o
x
x
F
x
x 0 FFG 3000 FDF cos 45o
Figure 2-19 – FBD, Node A
or FFG = 2000 lb. Figure 2-19 shows a Free Body Diagram in the
vicinity of Node A. Note that we know that the
The vertical force equilibrium equation is value of FAG is 2828 lb. The FBD shows shear
force, V, and bending moment, M, at a position
F
y
y 0 FCF FDF sin 45o a distance x to the right of Node A.
Chapter 2 Page 37
Truss structures
F
x
x 0 P FAG cos 45o
so P=-2000 lb.
F
y
y 0 FAG sin 45o 1000 V
M
z
z 0 1000( x) FAG ( x) M
Figure 2-20 - Final results for truss element
internal forces
so that M=1000x in-lb. (0< x < 20) This
bending moment equation is valid as long as x is 2.10 Design terminology and weight
less than 20 inches. calculation
Although we are reviewing mechanics and
To the right of Node B, the Free Body Diagram applying it to simple trusses, this is a good time
changes since x>20 inches. (Draw the FBD to to introduce some new material. In particular
confirm this). The equilibrium equations are: we need to make sure that we understand some
important terms like “Factor of Safety.” We
F
also need to reinforce the fact that weight is
x 0 2000 P
important in structural design and also
x
understand the importance of putting together
F
y
y 0 1000 V good structural topology.
Chapter 2 Page 38
impossible. For instance, most aircraft (and product of a limit load times a factor of safety.
passengers) go through their entire service lives Applying an ultimate load to a structure will
without being subjected to ultra-violent result in failure, defined as permanent
maneuvers. However, loads from such deformation of some parts, fracture or rupture of
maneuvers might happen to one or two aircraft others or buckling that severely limits the
in the fleet over a period of decades. If applied, operation of the structure; the overall structure is
a limit load must be resisted by the structure not allowed to fail or collapse catastrophically.
without developing stresses so large that
permanent deformation would result when and if The FAR 25 requirements specify that the
the load is removed. Permanent deformation ultimate load be supported by the structure for at
stress levels, how high the stress levels need to least three seconds.4 The factor of safety is
be, depend on the materials used. usually equal to 1.5 for transport or military
airplanes with human occupants. For missiles or
Limit loads are computed or determined for for unmanned vehicles, this factor of safety
numerous operational conditions, including might be as low as 1.25.
symmetrical and unsymmetrical maneuvers.3
These load conditions include flight at different Originally, the factor of safety was intended to
altitudes and airspeeds with different weight account for unknowns and uncertainties in
payloads, locations and fuel weight. FAR materials quality, manufacturing and processing,
Section 25.321 states: assembly or operations. In the early days, this
ratio was 2.0 but was later reduced to its present
“Flight load factors represent the ratio of the value of 1.5. This factor is multiplied by the
aerodynamic force component (acting normal to limit loads to obtain the ultimate load factor.
the assumed longitudinal axis of the airplane) to
the weight of the airplane. A positive load factor Since we need to preclude failure, the
is one in which the aerodynamic force acts establishment of limit loads and ultimate loads
upward with respect to the airplane. naturally leads to tying them to allowable
stresses in the structure. How much stress is
Considering compressibility effects at each allowable if we need to prevent permanent
speed, compliance with the flight load deformation? How much is allowable if the
requirements of this subpart must be shown structural members are not allowed to buckle?
At each critical altitude within the range of
altitudes selected by the applicant; To understand the answers to these questions,
let‟s return to one of our truss examples.
At each weight from the design minimum weight
to the design maximum weight appropriate to
each particular flight load condition; and
For each required altitude and weight, for any 4
“The structure must be able to support ultimate
practicable distribution of disposable load
loads without failure for at least 3 seconds. However,
within the operating limitations recorded in the when proof of strength is shown by dynamic tests
Airplane Flight Manual.” simulating actual load conditions, the 3-second limit
does not apply. Static tests conducted to ultimate
Accidental exceedance of the limit load cannot load must include the ultimate deflections and
be excluded as a possibility so we are very ultimate deformation induced by the loading. When
careful to consider ultimate loads. An ultimate analytical methods are used to show compliance with
load for a structure is computed to be the the ultimate load strength requirements, it must be
shown that: (1) the effects of deformation are not
significant; (2) the deformations involved are fully
3
T. Lomax, Structural Loads Analysis for accounted for in the analysis; or (3) the methods and
Commercial Transport Aircraft: Theory and Practice, assumptions used are sufficient to cover the effects of
AIAA Education Series, Reston, Va., 1996. these deformations.”
Chapter 2 Page 39
Truss structures
Fi Fi
Ai
allowable allowable
The cross-sectional area of the compressive
member is now computed to be:
As a result, the weight formula becomes
Fi L 12 FBC
wi i Li ABC BC
allowable E
Now, the question becomes, “where do I find the The weight of member BC is
stress allowables?” The two members in the
example resist different loads. Member AB has
a tensile internal force while BC is compressed. 12 FBC 2
wBC LBC
2
The failure mode for AB will be tensile failure,
E
either due to exceeding a prescribed stress like
the yield stress or exceeding the rupture stress. with
Chapter 2 Page 40
element internal load, the weight of the
compressive element varies as the square root of
with the internal compressive load.
For compressive members, the material strength wspruce = AB, 0.149 square inches
parameter such as yield stress or ultimate stress, BC, 0.612 square inches
does not appear in the weight equation. Instead,
In this case, spruce is a better choice on the basis
the weight to stiffness ratio appears. For of element weight. However, we have not
compressive loads, the element axial stiffness factored in the joining mechanisms or the fact
becomes important; the larger the stiffness, the that “wooden wire” or very thin toothpick-like
smaller the cross-sectional area required. tension members might not be a good idea, from
a manufacturing standpoint. Note also that this
Notice also that very long members subjected to comparison is only for a single applied load and
compression will have very large weights since there are many other loads and considerations to
weight varies as the square of the element be taken into account. We‟ll have more to say
length. On the other hand, while weight of about this when we discuss materials in Chapter
tensile elements is linearly proportional to the 5.
Chapter 2 Page 41
Truss structures
If we plot the total weight of the two elements as computed using the simple strength of materials
a function of h, we find that, for spruce, the bending stress formula = My/I.
minimum weight is reached when the height h is
approximately 45 inches. The difference in The beam cross-section in (a) is square and the
weight between that found for h=20 inches and cross-sectional area is constant from the beam
h=45 inches is about 12%. root to the beam tip. The dimensions of the
cross-section are sufficient to make sure that the
When the two-bar truss is made of aluminum, stress at the beam support does not exceed the
the minimum weight occurs when h is near h=75 100 N/mm2 maximum allowable stress. In this
inches. The difference in the weight with h=20 case the stress is considerably smaller at other
inches compared to h=75 inches is about 28%. points on the beam so material there is wasted.
The optimum structural topology depends on the The weight of this design is 12 kg.
construction material, as well as the loading.
Simple models such as this example can help In (b) the cross-section is rectangular, but also
determine the most efficient structural constant along the beam itself. Once again, the
arrangement. Let‟s consider another example to large bending moment at the root determines the
make this point. height of the cross-section at every point on the
beam. The weight of this beam is 6.1 kg.
Figure 2-21 shows five different structural
configurations. Each structure transfers a load P In (c) the width is the same from beam root to
= 1000 Newtons to a wall 1 meter to the left of tip, but the cross-sectional height, h, varies so
the load application. The two basic design that the maximum stress at any cross-section is
choices are beam configurations (a, b, and c in 100 N/mm2. In this case we have a “fully-
Figure 2-21) or a truss configuration (d and e in stressed” beam design. The weight of this
Figure 2-21). Three different beam cross- design is 4.1 kg.
sections and two different truss geometries are
considered. The effect of using greater beam depth and the
effect of creating a fully stressed design are
The material for each of the five candidate apparent in the weight values. On the other
designs is steel; the maximum allowable tensile hand, the manufacturing cost of (c) will
stress is 100 N/mm2. The beam stresses are probably be greater than either (a) or (b).
Figure 2-21 – Geometry of structural choices to transfer a 1000 N. load to a wall support (Jorgen Kepler,
“Structural optimization as a design and styling tool-with emphasis on truss structures,”) The dimensions
given on the figures are given in millimeters.
Chapter 2 Page 42
Design (b) might be susceptible to lateral “law” states that energy cannot be created or
buckling instability. destroyed; it can only be changed from one form
to another. In the structures application, external
Cases (d) and (e) are truss alternatives to the work is converted to internal deformation or
beam design. In this case, hollow round tubes “strain energy.”
are used as compressive truss elements. The
heavy lines identify truss elements subjected to Strain energy is a useful concept in structural
compression while the thinner lines are truss analysis. Deflection calculations begin with
elements in tension. For truss elements, we need strain energy.
to consider not only the maximum tensile stress,
but also guard against buckling. The cross- Consider the concept of strain energy density for
sectional areas of the compressive members a linear elastic material, a material whose stress,
must be large enough to prevent buckling while , and strain, , are related through the
the areas of the tensile elements are large relationship where the constant e, is the
enough to keep the stresses below the maximum Young‟s modulus. If there is only a single stress
of 100 N/mm2. involved, this density is written as:
Chapter 2 Page 43
Truss structures
For the example two bar truss, the total strain Our discussion was limited to the application of
energy is the sum of the energies stored in each a single load to a structure and to statically
member. This means that we have determinate structures. In the usual case where
we have loads applied, removed and then
U total U AB U BC replaced with others, such as the case shown in
Figure 2-21 the situation becomes more
complicated.
20 2
1
1 400 P 2
U total P LAB
1 2 h
EAAB 2 EABC h 2
* 20
2
The relationship between the deflection where
the load P is applied, in the direction of P is
found by differentiating the strain energy with
respect to the applied load, P. This gives us:
20 2
1
U total h
PLAB 20 400 * P Figure 2-22- Two truss designs with two loads,
P EAAB EABC h 2 separately applied, at different times.
or
When the load P1 is removed and P2 is applied
20 2 the two left members carry no load. Inserting an
1 LAB
h 8000 additional element, such as shown in the truss
P diagram to the right in Figure 2-22 will tie in
EAAB EABC h 2 these members and cause them to participate in
the load carrying arrangement. Insertion of this
additional element also makes the structure
statically indeterminate so the analysis is more
This is the result of application of Castigliano‟s
complicated.
theorem.
We also briefly considered a fundamental
2.13 Summary
theorem of structural mechanics, Castigiliano‟s
We started Chapter 2 with the goal of reviewing
Theorem, to show how these relationships help
mechanics, including statics and Free Body
us to find deflections. Later we will develop an
Diagrams. Our vehicle for doing this was truss
analysis method, the finite element method; that
construction even though this design approach
will make truss analysis easier and automated.
sees limited use in aerospace structural design.
However, the fundamentals we have reviewed
We expanded upon this review by considering
here are still as important as they were 100 years
the estimation of structural weights and showed
ago when airplane enthusiasts took to the air.
how much of a penalty we incur when we load
truss members in compression instead of
2.14 – Historical notes
tension.
The strong influence of railway engineers and
civil engineers such as Octave Chanute on early
Configuration or “topology” is important for
aircraft structures has already been mentioned in
minimal weight structural design. So is the
Chapter 1. Some famous engineers and
choice of materials. We saw that we got
scientists were contributors to the engineering
different optimal geometries depending upon
mechanics theories that advanced bridge
which materials were used (our case considered
construction and, ultimately, aircraft structures.
aluminum and wood).
William J.M. Rankine‟s name is usually
Chapter 2 Page 44
associated with thermodynamics and the
Rankine cycle. However, he was an early
contributor to engineering mechanics. His book
Manual of Applied Mechanics was published in
1858 in England and was for years one of the
primary sources for structural theory and what
we know today as strength of materials.
Rankine, whose name we associate with
temperature measurement, penned a poem about
engineering units, one stanza of which reads:
5
Elementary theory and calculation of iron bridges
and roofs, Prof. Georg August D. Ritter, Professor
Polytechnique School Aix-La Chapelle, (Aachen)
translated by Henry Riall Sankey, Lt. Royal
Engineers, 1879.
6
S.T.G. Andrews and S.F. Benson, The Theory and
Practice of Aeroplane Design, E.P. Dutton &
Company, New York, 1920.
Chapter 2 Page 45
Truss structures
2.16 Homework Problems for Chapter Two Problem 2-4: The structure shown in the figure
Some of these problems will be assigned. represents an idealization of an engine mount.
Solutions should contain Free Body Diagrams Find the internal forces in the two-force
and equilibrium equations. Answers should be members. The supports on the right are pinned.
identified by outlining them. A summary
diagram of the answers similar to Figure 2-20
should appear for each problem solution.
Chapter 2 Page 46
The two idealizations shown in Problem Figure
2-6(b) and (c) and represent the math models of
the two competing designs. Complete the
following tasks.
c. The summation of the absolute value of the Problem 2-6 (c) - Competing spacecraft truss
internal forces times the lengths of each
members is written as and a) Compute the reaction forces and the internal
has the units of work. For the structure in 2- loads in members BC, BG and FG.
6(b), n=4 while for the structure in 2-6(c), b) Three members carry no loads. Identify
n=2. Compute the value of J for each those members by inspection and logic.
structure. Which has the smaller value of J?
(Answers: (b) J=400,000 in-lb; (c) Problem 2-8: The seven-element truss shown
J=450,000 in-lb.) below is fixed to two supports. The common
formula m=2j-3 used to determine static
indeterminacy does not apply to this structure.
Explain why.
Chapter 2 Page 47
Truss structures
Problem 2-9
Determine which truss structures shown below
are statically indeterminate and their degree of
indeterminancy.
(a) (b)
(c)
Chapter 2 Page 48
Inertia loads
A half a world away, a commercial airliner There are thousands of different aircraft loading
encounters turbulence. Overhead bins open and conditions. Identifying which of these
flight attendants and some passengers are thousands of loads are critical is important.
thrown to the ceiling while other passengers “Critical” means “might cause failure unless we
strain against their seatbelts. Still further away, design the component correctly.” An applied
fifty miles above the Earth, a space shuttle is load may be critical for one part like a bulkhead
returning from space, encountering re-entry in the front end of the airplane, but create no
forces that press its occupants firmly into their problem for other components else.
seats. Their vision blurs.
Airplanes loads come from diverse sources.
These three events recognize the importance of These include1
inertia loads in the design of aircraft and
spacecraft. What exactly are inertia loads? Component and payload weight
Air loads: lift, drag, gusts
As indicated in Chapter 1, aerospace structural Acoustic loads
design has a set of logical steps with one major Thermal loads
goal: produce a minimal (but not necessarily Landing loads: touchdown, arresting
minimum) weight structure that is reliable, cost- Take-off loads: runway taxi, catapult,
effective and, above all, safe to operate under a assisted (e.g. jet or rocket) take-off
variety of conditions. This task requires seven Power plant loads: thrust, engine
steps: torque, gyroscopic effects
special loads: towing, refueling,
1. Define top level vehicle system
pressurization, cargo, weapons recoil,
requirements
bomb release, blast effects.
2. Define airframe structural system design
requirements, including metrics for success
Most of the loads resisted by structural elements
3. Define load conditions, identify the critical
are due to aircraft weight (actually the aircraft‟s
loads
mass) and the fact that the airplane accelerates
4. Choose materials and layout the structural
during flight.
geometry and components
5. Perform analyses for sizing and to check for
1
compliance, compute factors of safety and The Federal Airworthiness Regulations (FAR) lists
design margins the loads that must be considered before a
6. Calculate system weights commercial aircraft can be certified as airworthy
Chapter 3 Page 49
Inertia loads
“Inertia loads” are a primary source of loads for object into a static load is called d‟Alembert‟s
aircraft and arise because of accelerations Principle.
created by engine thrust, maneuvers or gust
conditions. The concept of the “load factor” or Inertia forces have special properties. They are
“g-force” is used to describe inertia loads vectors that:
created by acceleration and resisted by the
structural components to keep it acting as a unit. act through the object‟s center of mass
have a direction opposite to the direction of
3.3 D’Alembert’s Principle – inertia loads acceleration
inside the vehicle are proportional to the mass of the object.
One of the most important loads on aircraft and
spacecraft is called the inertia load. How do we To illustrate this concept let‟s use a simple
compute inertia loads? Newton‟s Second Law is model with a weight, W, suspended from a wire,
written as F=ma; we compute airplane as shown in Figure 3-1. The weight is subjected
accelerations once we know the airplane mass, to gravitational force, W = mg, and to a force T
m, and the net forces, F. We include all external (as in tension) in the wire. The mass accelerates
forces and gravity acting on the airplane in our upward by an amount “a.” Newton‟s Second
calculation of F. This calculation allows us to
compute vehicle trajectories, but structural T
analysts must compute forces inside the a = (T - W)/m
accelerating airplane structure in a non-inertial T=ma+W
reference system.
Computing inertia forces is simple. While we As an observer with feet planted firmly on the
write the dynamics equation F=ma so that we ground, we see that the force T accelerates the
can compute the acceleration a=F/m, we can weight, W, but, if we are moving with the wire
also write the equation as F-ma=0. we don’t see the acceleration, we feel it. The
tension in the wire seems to counteract an
Written this way, this equation looks like a static acceleration “force” ma or (W/g)a.
equilibrium equation so it is called the “dynamic
equilibrium equation.” The term (–ma) is called D'Alembert‟s Principle transforms an
the “inertia force.” accelerating system into an equivalent static
system by adding so-called “inertia forces” and
Inertia forces are also called “D‟Alembert "inertia moments" to the forces acting on the
forces” to honor the French mathematician Jean system. D‟Alembert‟s Principle is just another
le Rond d‟Alembert (1717-1783) who first way of writing Newton‟s Second Law, but
formulated the concept and the process of written in a form more useful to structures
converting the mass times the acceleration of an analysts.
Chapter 3 Page 50
Figure 3-3 – Aircraft launch loads;
Source – David Peery
Chapter 3 Page 51
Inertia loads
We now divide both sides by the vehicle weight, In the y-direction we have a “1-g” condition.
W, to get a number called the load factor in the All mass items inside the airplane exert their
x-direction, nx. normal weight as forces downward.
Chapter 3 Page 52
Summing external forces (excluding gravity and direction is determined as follows: positive load
the inertia force) in the positive x-direction, we factors produce loads in the negative direction;
find that the total is: negative load factors produce loads in the
positive direction. Let‟s illustrate this with the
Fx
48,750
example shown in Figure 3-5.
nx external
3.25
W 15,000
Suppose we are in the cockpit of this airplane. An airplane is maneuvering so that the following
A pencil weighing one ounce when the airplane load factors are computed: nx = -1.5; nz = 2. An
is stationary will have two “weight” components item weighing 10 pounds is suspended from a
during launch, one in the x-direction and one in wire as shown in Figure 3-5. If the load factors
the y-direction. The magnitude of these weights remain constant for a long enough time, find the
and their directions are important for loads angle that the wire will make with respect to the
calculation. z-axis and find the tension in the wire.
Chapter 3 Page 53
Inertia loads
1 4 ft
z
m ax
4 ft T
1 3 ft
R1 W R2
0 .3 R 2
7 0 ft
3.7 Example – aircraft touchdown Next, we sum moments about the airplane center
The airplane shown in Figure 3-8 weighs W = of gravity.
120,000 lb. During touchdown its engines apply
a reverse thrust of 8,000 lb. and the rear wheels
M c. g . 0 56R1 14R2 4T 13(0.3* R2 )
exert 30% of their normal force during braking.
external
7 0 ft R1 29,500 lb.
Figure 3-8- Example problem, aircraft landing
(Source-W.L. Hallauer) To find ax we use the x force equilibrium
equation:
Find the:
vertical reactions at nose and main
landing gears
acceleration in the x direction
Chapter 3 Page 54
(a) Determine the vertical airplane acceleration
The negative sign tells us the aircraft is and load factor at the instant that the airplane
accelerating in the negative x direction (acting to touches down.
the right).
The positive z-direction is upward. The
To compute the load factor in the x-direction we unknown is the acceleration so we use the
use: dynamics equation:
W
F 8000 0.3* R2
Fx 45,000 45,000 30,000 ax
g
x
W 120,000 W
35,150 ax 60,000 lb.
0.293 g
120,000 ax
2
g
Anything with mass will exert a force to the
right (negative x-direction) equal to its weight
times the load factor (-0.293) due to its Now, compute the load factor, nz:
deceleration. As passengers in this airplane we
will feel ourselves being pushed forward into the F x 45,000 45,000 nzW 30,000* nz
seat. nz = 3
3.8 Calculating internal reaction forces and (b) Compute the compressive force and bending
moments moment in the left
When vehicles accelerate, how do we calculate landing gear
the structural component shear forces, bending strut.
moments, and other internal actions caused by
the combination of both static and dynamic Figure 3-11
loads? Let‟s answer this question with an shows the FBD of
example. the left wheel and
strut assembly
The 30,000 lb. airplane shown in Figure 3-10 is with the inertia
landing so that the ground reaction on each main force included.
wheel is 45,000 lb.2 Aerodynamic lift is The center of the
negligible during this landing impact. tire footprint is
offset from the
oleo strut Figure 3-11- Wheel
centerline by 6 assembly
inches. The
wheel assembly weighs 500 lbs. We will solve
for the oleo strut compression force C and the
bending moment, M, by summing forces and
moments.
C=43,500 lb.
M
2
This example is adapted from Aircraft Structures
cut 6* nz *500 6* 45,000 M 0
by David J. Peery, (1950) pages 50-52.
Chapter 3 Page 55
Inertia loads
M = 261,000 lb-in.
Each section has a c.g. location xi measured with
(c) Calculate the shear force and bending respect to the vehicle c.g. and a weight, wi. The
moment in the wing cross section 240" inboard relationship between the component weights,
from the wing tip (which is outboard of the c.g. locations and the vehicle center of mass is
landing gear/engine mount). This wing section such that
weighs 1,500 lb and has its center of gravity at 4
tail
Summing moments at the cut section we have Io nose
mx 2 dx
We could have used the inertia load of 30,000 To calculate the load factor on the vehicle in the
lb. acting downward, together with the actual z-direction we use
15,000 lb. weight to get the same results.
Fz Wnz
3.9 Adding in the angular acceleration
“forces” where W= missile total weight.
Figure 3-13 shows a missile acted upon by
external forces that produce a net force in the z The sum of the external forces times their
direction as well as a counter-clockwise pitching distance from the vehicle c.g. is equal to
moment about the missile c.g. The vehicle sub- M o I o .
divided into four pieces for analysis purposes.
These pieces may be actual modules or arbitrary
cut sections chosen where we wish to calculate
internal loads.
Chapter 3 Page 56
The vehicle has an angular pitching acceleration 2. An extra inertial load due to angular
as the result of the net clockwise moment. As acceleration about the vehicle c.g., located at
a result, two effects must be taken into account the component c.g. written as
when computing the inertia loads. but in a direction opposite to the
acceleration itself
First of all, because we are rotating the missile 3. An inertia moment equal to M Ri Ii in
in pitch there are accelerations with respect to the direction opposite to the angular
the vehicle c.g. not accounted for in the load acceleration. Although we usually draw this
factor computation. As indicated in Figure 3-13, moment as acting around the component c.g.
the translational acceleration at any position x it is a free vector.
from the missile c.g. is computed from the
formula The final inertial load results are shown in
Figure 3-15.
ari x
Chapter 3 Page 57
Inertia loads
Initially the cable tension is 20,000 lb. and its Rotational acceleration generates additional
line of action passes through the airplane center inertia forces in addition to those given by the
of gravity. The inclination of the cable is 20 load factors. These loads are shown in Figure
degrees with respect to the horizontal. This 3-17.
produces a 30,000 lb. reaction force in the main
landing gear located 20 inches aft of the
c.g.
Chapter 3 Page 58
Figure 3-18-Internal forces generated by one pound
components located at nose, c.g. and tail of an
aircraft due to translational and rotational
accelerations. The component shown in light green
is the result of a common load factor while the
yellow triangles are the result of angular
acceleration about the aircraft c.g..
Figure 3-20-Partial FBD of aft end of missile
3.10 A second example - missile joint bending
moments in a missile with rotational Data for the aft section shown in Figure 3-20 is
acceleration as follows: 1) mass moment of inertia, Io, about
Inertia moments are important. Consider this its own c.g. is Io = 16,670 slug-ft,2; 2) aft
problem. A missile in flight, depicted in component weight = 3000 lb.
Figure 3-19, is subjected to loads during
turbulence. Resultant forces due to the The first step is to convert this information into
distributed external side loads are shown in the inertia loads. This includes an inertia moment.
figure. As a result, a load factor nz = 1 and a These forces and moments are shown in Figure
counter-clockwise angular acceleration = 1.08 3-21, together with the shear force, V, at the
rad/sec.2 are created. missile joint and the internal bending moment,
MR .
An assembly joint is located 30 inches aft of the
missile c.g. Compute the bending moment and The c.g. of the component is subject to two
shear force in the field joint. A Free Body inertia forces. The first component is due to the
Diagram is shown in vehicle load factor and is 3000*nz. This load
Figure 3-20. acts at the missile component c.g in the negative
z-direction as shown in Figure 3-21.
Chapter 3 Page 59
Inertia loads
that the missile component is rotating about the This gives the result that MR = 1,225,000 in-lb.
missile c.g. with angular acceleration 1.08 in the clockwise direction.
rad/sec.2 This angular acceleration produces a
translational acceleration equal to aR 130 The shear force V is found by summing forces.
inches per sec2 in the positive z direction. It also
creates an inertia force located at the component V 1090 3000 20,000 3000 0
c.g of (3000/g)*aR directed to the left, in the V 18,910 lb.
negative z-direction.
The negative sign indicates that the shear force
In addition to these inertia forces, an inertia direction is opposite to that originally chosen. If
moment must be added to account for the fact the joint cannot furnish these internal reactions
that the missile component has rotational then the missile will break apart.
acceleration. Like the d‟Alembert inertia force,
the inertia moment is computed by multiplying 3.11 Example – another idealized missile in
the mass moment of inertia of the missile flight
component, Io = 16,670 slug-ft.2, by the missile Our structural dynamics group has chosen a
angular acceleration to produce the inertia lumped parameter or “stick” model shown in
moment that acts on the missile component. Figure 3- to analyze a missile in flight. .
or
M inertia 216,000 in lb.
Chapter 3 Page 60
c) Compute the internal axial force between aR1 22.5 in.
masses 5 and 6. sec.2
aR 2 12.5 in.
First, let‟s solve let‟s locate the missile c.g. sec.2
We‟ll use the temporary coordinate system in aR 3 2.5 in.
Figure 3-. sec.2
aR 4 7.5 in.
sec.2
aR 5 17.5 in.
sec.2
aR 6 27.5 in.
sec.2
The inertia loads due to these accelerations are
computed from the formula
w
FRi i aRi
g
Figure 3-24-Lumped mass model cg. locations
The direction of these inertia forces is opposite
to the accelerations, as indicated in Figure 3-26.
w X i i X c. g . w
i
Figure 3-25- Angular acceleration induced Figure 3-27- Final result, inertia forces due to
translational accelerations discrete masses, units are pounds.
The z-direction translational accelerations due to The bending moment at the missile c.g. due to
the angular acceleration are as follows: the inertia forces depends only on the forces in
Chapter 3 Page 61
Inertia loads
the z-direction. The FBD for this calculation is To obtain the final bending moments and axial
shown in Figure 3-28. forces inside the missile we would need to know
the distribution of external aerodynamic forces
over the missile and superimpose these forces
and moments onto the inertial effects we have
just computed.
The axial force is computed by drawing the FBD Figure 3-31-Smokestack structural model showing
shown in figure 3-30. pinned support at the bottom
Chapter 3 Page 62
wL2
M 0 w sin xdx
L
pin
2
3 g
sin
2 L
w 3 x
f ( x) maR x w sin
g 2 L
Figure 3-32- Gravity and inertia loads are
distributed along the smokestack. The gravity load Note that this expression depends only on the
producing stack compression is shown in the figure. smokestack angle.
The angular acceleration also produces an inertia We have drawn the beam model in Figure 3-33
force distribution, written as f(x) and shown in oriented in a horizontal direction to show all
Figure 3-. The units for this distributed force are forces, including the inertia load and
lb/ft. This transverse load produces a shear and components of weight in the transverse direction
bending moments in the smokestack. Unlike our (wsin).
earlier examples, the inertial reference point for
angular acceleration is the fixed pin at the
smokestack base, not the smokestack c.g.
F f ( x) wsin dx
L
where the mass moment of inertia about the pin y 0 Ry
0
is defined as or
L L3 w L3 3
I pin mx 2 dx m
0 3 g 3
Ry wL sin wL sin 1 wL sin
4 4
The bending moment for this problem is found
The sum of the moments about the pin is: by drawing the FBD shown in Figure 3-35. This
FBD shows the smokestack cut at an arbitrary
Chapter 3 Page 63
Inertia loads
Chapter 3 Page 64
total
Mo a w L x cos
a 4 a2
4
or
wL cos
L tan x 3
L 1
2
total a L 2 L
x 2 x
a2
We cannot develop this formula further without Load factors can act in any
knowing more smokestack geometry details, in direction
particular the ratio L/a. On the other hand, we
can understand that this formula should predict a
negative result (compression) unless we want
the material to fail. However, we expect the
smokestack will begin to fail in the region near
x/L = 0.35. Now, go watch then topple a
smokestack and see if we are right!
3.13 Summary
This chapter has discussed the concept of inertia
loads, where they come from and how to
calculate them. We defined and used the load
factor used by loads groups to define loads due
to acceleration of the vehicle. These loads,
combined with aerodynamic loads, can place
significant stress on the airframe in flight and Computing and interpreting load
during landing. factors
Chapter 3 Page 65
Inertia loads
3.14 Chapter 3 - Homework problems time a jet vane malfunction causes a force
F=20,000 lb. to be applied as shown.
Problem 3-1
Chapter 3 Page 66
(a) Compute the tension in the cable at this
acceleration.
(b) Compute the internal loads in each of the
seven elements due to the 60 pound hoist
and the accelerating 1000 pound mass.
Problem 3-5
Problem 3-7
Problem 3-5 Figure
Chapter 3 Page 67
Inertia loads
Chapter 3 Page 68
Problem 3-11
An aircraft is flying with a trimmed lift, L, equal
to the weight of the aircraft, W when a gust with
lift L = 1.5W strikes the aircraft. The resultant
of this gust produced lift is 90 inches aft of the
c.g. Ipitch for the entire airplane = 2.25x107 lb-
insec.2 computed about the aircraft c.g.
Chapter 3 Page 69
Matrix method analysis
4.1 Review - Statically indeterminate structures “nodes.” The two pieces must fit together –
Let‟s review what we know about statically before and after the loads are applied. This
indeterminate structures by considering a simple gives the solution its “uniqueness.” There are
example. The bar (also called a “rod”) shown in three nodes in this example. Two nodes are at
Figure 4-1 is subjected to a single load, Po. As the bar assembly boundaries and the third node
indicated, the bar consists of two components of is at the connection point.
different sizes.
Chapter 4 Page 70
Matrix methods
displacements of each bar. We have two From these results we conclude that the force in
equations: bar 2, the right bar, F2, is compressive and that
the force in bar 1, the left bar, F1, is tensile.
and
u2 – u3 = 0 We also see that the reaction forces at the wall
attachments are also equal and opposite to the
We also have two other equations that will allow forces F1 and F2. We also see that if the cross-
us to solve the problem. These are the bar sectional area of bar 2 decreases to nearly zero
displacement equations which read: then the external load is primarily resisted by bar
1.
F1 L1
u 2 u1 From these results we have two conclusions.
EA1 First of all, when we have a statically
and determinate structure, the internal loads depend
F2 L2 on the size of the individual components, in
u3 u 4 particular, the stiffness of the elements. In this
EA2
example, the bar stiffness is given as the
FL quantity
This gives us u2 1 1
EA1 EA
k1 1
and L1
F2 L2 Second, we need to take into account the support
u3 conditions and also enforce “compatibility”
EA2
conditions to make sure that the structure
deforms as a unit and doesn‟t come apart as it
Combining these equations, we have: deforms.
Chapter 4 Page 71
Matrix methods
element is heated it will try to change its length. This equation gives us the ability to solve for F1.
This change in length would be:
L T L
Chapter 4 Page 72
Matrix methods
element methods are used to solve structural In fact, structural analysis methods range from
analysis problems with incredible geometric the simple to the complex. At its simplest, the
complexity. Our discussion is confined to truss analysis of truss structures, with all of the
and bar element problems such as those we fundamental mechanics it entails, has been a tool
considered in Chapter 2 and in this section of for airplane structural design since the Wright
Chapter 4. Brothers first took to the air. The complexity of
structural designs, the high speeds involved and
4.4 Finite element methods the need to meet schedules were fortunately
In practice we are faced with the need for the accompanied by technology developments such
analysis of a wide variety of structural parts with as the high-speed digital computer. This gave
different, challenging geometries such as that rise to a methodology called matrix methods as a
shown in Figure 4-6. In addition, the way to model structures for a variety of loads,
information we need ranges from stresses due to including vibratory loads. Chief among these
loads and vibration dynamic response. matrix methods is the finite element method.
Providing this information is a challenge.
Finite element methods are a matrix method
approach to formulating and solving structural
problems. During the past 40 years advances in
this area have been phenomenal. The central
reason for this has been the development and
expansion of computers.
Chapter 4 Page 73
Matrix methods
One type of matrix method treats forces as the Deriving the element stiffness matrix
unknowns and this is called the “force method.” We will develop a matrix relationship between
The other type of analysis treats displacements the end forces at the nodes and the nodal
as the unknowns and this is called the displacements for all possible combinations of
“displacement” or “stiffness” method. We will forces and displacements for the rod element
consider only the stiffness method, a method shown in Figure 4-8. The restriction is that the
that has come to dominate finite element relationship is linear, that is
analysis.
F11 a11 b12
Before we go any further with definitions and
descriptions, let‟s go back to our examples in F12 c11 d 12
part one of this chapter and reformulate the
problem. How do we get the constants a, b, c and d? First
take the rod element and fix the right hand end
Our problem consisted of two connected bar so that 12 0 . To be in equilibrium, the two
elements. As indicated in Figure 4-7, if we end forces must be equal and opposite
apply end loads to either of these elements then
the deflection of the element is given as: F12 F11
F kx
This means that
EA
k
l F11 a11 F12 c12
or c = -a.
Chapter 4 Page 74
Matrix methods
EA EA
F11 11 12
L L
EA EA Figure 4-9- Two element rod arrangement with
F12 11 12 external loads
L L
which is written as
Figure 4-10 – Left hand bar element with element
F11 EA 1 1 11 k k 11
forces
F12 L 1 1 12 k k 12
At the left end, the reaction force R1 is in
or equilibrium with the element force F11 so that
F11 = R1. Similarly, at the right hand end of the
F11 k k 11
assembly the reaction force R2 is equal to the
F12 k k 12 elemental force F22, or F22=R2. Now the
elemental matrices are written as:
This seems like an overly complicated way to
write a simple relationship, but this expression R1 k1 k1 u1
accounts for every possible combination of
F12 k1 k1 u2
internal forces and displacements and after all it
is a matrix equation. Let‟s write the relationship
for the second bar element. And
F21 k2 k2 u2
F21 k2 k2 21 R3 k2 k2 u3
F22 k2 k2 22
Finally, we have the force equilibrium
For our first step in solution, we simply relationship at the center node, node 2. This
guarantee compatible deflections by observing relationship is:
(and writing): u1 = 12; 22 = 21 = u2; and, u =
22. We will do this a little more simply later.
This relationship requires that we write the
This means that we have two element matrices
element forces at node 3 in an expanded form as:
that are written as:
u1
F11 k1 k1 u1
F12 k1 k1 0u 2
F12 k1 k1 u2 u
and 3
F21 k2 k2 u2 u1
F22 k2 k2 u3 F21 0 k 2
k 2 u 2
u
3
Step 1 is a coordinate transformation and it
guarantees that there will be compatible
deformations – the element nodes will stick Then, our equation for force equilibrium
together where necessary. For Step 2, we need becomes:
to make sure that there is force equilibrium at
each of the three nodes. u1
Po F12 F21 k1 k1 k 2 k 2 u 2
u
3
Chapter 4 Page 75
Matrix methods
R1 k1 k1 0 u1
Po k1 k1 k 2 k 2 u 2
R 0 k2 k 2
2 u 3
The minus sign on the result for R2 means that
or our initial assumption about its direction was
incorrect. In fact, the R2 force vector goes from
right to left, opposing the applied force Po.
R1 k1 k1 u1
and then we solve for the reactions.
F12 k1 k1 u 2
Chapter 4 Page 76
Matrix methods
or
This gives us the result:
u1
F11 k1 k1 0
u2
F12 k1 k1 0 u
Compare this result to the result we found earlier 3
in this chapter and you will see that they are the
same. The elemental stiffness matrix is a square matrix
and so is the system stiffness matrix so there is
4.6 Assembly of matrix equations one more operation we need to do here. The
The building block elements of our bar analysis matrix of zeroes and ones is called the
is the 2 x 2 elemental stiffness matrix, written in coordinate transformation matrix because it
local deflection coordinates. This relationship transforms local displacement coordinates into
computes elemental forces which must be added system (also called global) coordinates. We
together to get the applied forces or the need to do the same thing to the elemental
reactions. The two key operations are to forces. This is done by multiplying the elemental
transform the elemental deflections (the ‟s) into forces by the transpose of the displacement
“system” deflections (the u‟s) and to add the transformation matrix.
appropriate elemental forces together to get the
applied loads. Let‟s see how this can be done. 1 0 1 0 u1
0 1 F11 0 1 k1 k1 0 u
First of all, let‟s consider the left hand bar F k
k1 0 2
0 0 12
0 0 1
element, element 1. We can use a matrix u3
relationship to relate its elemental coordinates
(also called “local” coordinates) to the system F12 k1 k1 0 u1
structural coordinates. The relationship reads:
F12 k1 k1 0 u2
u1 0 0 0 0 u3
11 1 0 0
u 2
12 0 1 0 u This last matrix relationship is simply the
3
elemental stiffness matrix arranged to be
consistent and compatible with the system
Similarly, for element 2, the right hand element, matrix coordinates. You should satisfy yourself
we have: that when we do the same operations with the
right hand element we get the following:
u1
21 0 1 0
u 2 0 0 0 0 u1
22 0 0 1 u
3 F21 0 k2 k2 u2
F 0 k k2 u3
Now we can substitute these results into the 22 2
original elemental stiffness relationship.
Now, when we add forces we have
u
F11 k1 k1 0 1
u F11 0 R1
F12 k1 k1 0 2
u3 F12 F21 Po
0 F R
or 22 2
u
F11 k1 k1 11 k1 k1 1 0 0 1
u
k1 12 k1 k1 0 1 0 2
and
F12 k1
u3
Chapter 4 Page 77
Matrix methods
R1 k1 k1 0 u1 0 0 0 u1
P0 k1 k1 0 u2 0 k2 k2 u2
R 0 0 0 u3 0 k2 k2 u3
2 Figure 4-11-Two element bar/rod assembly
Chapter 4 Page 78
Matrix methods
1 k k k1
1 Figure 4-13 shows the rod element at an
k1 k1
k 1 2 orientation with respect to the horizontal.
1 k1 k2 k1 k1
We want to write our stiffness relationship as a
The determinant is defined as: function of the ui coordinates due to loads acting
in the ui directions. The first thing we will do is
determinant k1 k1 k2 k12 k1k2 to define a relationship between the two i
coordinates that lie in-line with the rod and the
“global” displacements that will lie along the
The determinant is non-zero and positive. This system coordinate directions. This relationship
is a feature of all correctly restrained structural
is a transformation matrix, defined as:
stiffness matrices. The solution then becomes:
1 1
i ij u j
1
u1 k2 k1 k2 P1
Where the transformation matrix has the
u2 1 1 P2 following elements:
k2 k2
cos sin 0 0
ij
k2 0 0 cos sin
R3 k2u2 P1 P2 k1
And k1k2
We must also transform the forces according to
R3 P1 P2 the relationship:
Chapter 4 Page 79
Matrix methods
T k k u
F F
T
global
ij 1
F1 1 0 1 0 u1
k
i ij i ij
k u2 F 0 0 0
EA 0 u2
F (1)
2
F3 L 1 0 1 0 u3
i
or
F4 0 0 0 u4
k cos 2 k sin cos k cos 2 k sin cos
0
k sin cos k sin 2 k sin cos k sin 2
F
i
global
k cos 2 k sin cos k cos 2 k sin cos
u j The second element has nodes 1 and 3 at the left
and right ends of the element so that = 135o.
k sin cos k sin 2 k sin cos k sin 2
1 1
In this case, cos , sin and the
The transformation matrix operation is a 2 2
powerful way to construct global/system
element length is 2L
matrices. These operational steps are required:
1 1 1 1
1) Write the transformation matrix as a F3 2 2 2 2 u
3
relationship between “new” and “old” F 1 1 1 1
2 2 4
F EA 2 2 u
(2)
4
displacements. In our case here the new i
F
1
1 1 1 u5
coordinates were the u‟s while the old 5 2 L
2 2 2 2
F6 1 u6
coordinates were the ‟s. 2
1
2
1
2
1
2
2) Pre-multiply the “old” stiffness matrix (the
matrix to be transformed) by the transpose
of the transformation matrix. Finally, the third element has nodes 1 and 3
3) Post-multiply the old matrix by the
located at the left and right ends. In this case,
transformation matrix.
= 90o and while .
4) Add the resulting matrices together to get
the system stiffness matrix.
1
This example is taken from R.E. Sennett, Matrix
Analysis of Structures, Waveland Press, Long Grove,
Illinois, 1994, pp.36-39.
Chapter 4 Page 80
Matrix methods
F1 0 0 0 0 u1 1 0 1 0 0 0
0 1
F 1
EA 0 0 1 u2
1P 0 0 0 u1
(3) 1 P u2
Fi 2 2
1 0 (1 1 ) 1 1 1
2 2
F5 L 0 0 0 0 u5 2 2 2 2 2 2
P3 EA
1 1 1 1 u3
0 1 0 1 u 2 2 u4
0 0
F6 6 P4 L 2 2 2 2 2 2
P5 0 0 1 1 1 1 u5
2 2 2 2 2 2 2 2
P6 6
u
1 1 1 (1 1
These three elemental stiffness matrices have 0 1
2 2 2 2 2 2
)
2 2
dimensions 4x4 and need one more operation
before they can be assembled into the 6x6 global
stiffness matrix. The transformation matrices to
do these operations are:
u1
u1 1 0 0 0 0 0 u2 When we restrain the support displacements to
u 0 0 u3
2 1 0 0 0 be zero (u1 = u2 = u5 = u6 = 0) the reduced
stiffness matrix becomes:
u3 0 0 1 0 0 0 u4
u4 0 0 0 1 0
0 u5
1
P3 EA 1 2 2
1
u6 2 2 u3
P4 L 1 1 u4
u1 2 2 2 2
u3 0 0 1 0 0 0 u2
u 0 0 u3
The solutions are:
4 0 0 1 0
L
u3 P3 P4
u5 0 0 0 0 1 0 u4
EA
u6 0 0 0 0 0
1 u5
and
u6 L
u4 P3 1 2 2 P4
EA
u1
u1 1 0 0 0 0 0 u2 The reaction forces are found to be:
u 0 0 u3
2 1 0 0 0
R1 = -P3 – P4
u5 0 0 0 0 1 0 u4
u6 0 0 0 0 0
1 u5 R2 = 0
u6 R5 = P4
One step remains. We pre- and post-multiply R6 = -P4
the three elemental stiffness matrices by their
respective transformation matrices and then add The final results are shown in Figure 4-15 and
the results to get the unrestrained system global satisfy static equilibrium.
stiffness matrix. This matrix reads as follows:
Chapter 4 Page 81
Matrix methods
F1(2) 2 P4
(2)
F2 2 P4
Chapter 4 Page 82
Matrix methods
Chapter 4 Page 83
Matrix methods
On the other hand, if we totally restrict the F12 EA2 2 T k2 k2 1
k2 2
movement of the ends then the total amount of
expansion must zero. To make the rod element F21 EA2 2 T k2
expansion free we add a compressive load as k k2 u2
2
shown in Figure 4-17. The total rod expansion
k2 k2 u3
is then the sum of the thermal expansion minus
the contraction due to the compressive load.
This is given by the relationship From this relationship we see that we have
P included the equivalent thermal force Pthermal at
L2 0 2 (T ) L2 thermal L2 either end.
EA2
Where Pthermal is the compressive load required to How does this get into the analysis? First of all,
hold the expansion to zero. Solving for Pthermal recall that the end forces for each element are
we find: summed to be equal to the applied loads. When
Pthermal EA2 2 (T ) we take the model shown in Figure 4-18 and
superimpose it on the rest of the structure, we
The only way to totally restrict the rod get some forces in addition to those mechanical
expansion is to place it between two immovable loads we apply. The best way to illustrate this is
walls and then heat it (by the way, if we cool the with an example.
rod then T is negative and the end force is
negative or tensile). When the heated rod is
placed as a component within a flexible
structure, the end load will have a value that
ranges between zero and Pthermal. This value will
depend on how much the rest of the structure Figure 4-18-Heated element with applied loads
resists the expansion.
Shown in Figure 4-19 is a two element structure
How do we tell the structure that there is an for which we have heated the right hand
element that wants to expand (or contract)? We element. When we assemble the elements, the
must put information into the problem at node equivalent thermal loads appear at the nodes to
points. Here is how we approach the problem. which they are attached.
Referring again to Figure 4-17, we heat the rod
element, but we also place a fictitious set of
compressive forces at either end to hold the
expansion to zero. Now the rod is heated but we
have forces that do not belong to the problem.
We remove these forces by applying equal and
opposite forces as shown. These forces are now
added to the structural model as shown in Figure
4-18.
Chapter 4 Page 84
Matrix methods
Since u1 and u3 are both zero, the solution to this restraint forces and the temperature. The
problem is thermal restraint forces are equal and opposite to
P EA2 2 T the forces in the diagram shown in Figure 4-19.
u2 thermal
k1 k2 k1 k2
The subtraction of the equivalent force
or effectively adds back the temperature. Now the
EA T equivalent forces are gone, the correct
u2 2 2
EA1 EA2 displacements are present, the structure is hot (or
cold) and the mechanical reactions are present.
L1 L2
Chapter 4 Page 85
Matrix methods
result of actual loading that produces strain We begin by generating the element matrices in
within an element; and, thermal loads that are global coordinates. For this we need the data for
the result of thermal expansions being restrained element transformation matrices. The table
by other structural components or boundary below shows the transformation angles used to
conditions. generate the element stiffness matrices for this
structure. (Remember that when we measure
Thermal loads are different than mechanical angles we hold a left end fixed and rotate
loads. Mechanical loads produce deformation counter-clock-wise, angles are plus in this
and stresses while thermal loads are due to direction.)
restricted thermal expansion. Thermal stresses
are never produced in statically determinate Element Angle (degrees)
structures. 1-4 90
1-2 -45
Let‟s consider a final example that will illustrate 1-3 -135
how to assemble a two-dimensional truss and
apply “thermal loads.” This example
The three elements are identified by the
configuration is shown in Figure 4-20. numbering of the left end and right ends.
Element 1-4 has node 1 at the left end and node
4 at the right end. The upper rod (1-4) is heated
an amount Theat (this is a temperature increase
but we have eliminated the from the term)
while the rod 1-3 is cooled an amount Tcool.
Both of these symbols represent changes in
temperature, not the absolute temperature.
Chapter 4 Page 86
Matrix methods
Chapter 4 Page 87
Matrix methods
R3
R4 0.5 0.5
0.5 0.5
EA Tcool
5R
2 EA 0.5 0.5 u1
R EA Tcool L 0.5 0.5 u2
6 2 0 0
R7 0 1
R 8 EA Theat
R3 0.25 0.1768
R 0.25 0.1768 Figure 4-25-Applied and external forces
4
R5 0.25 0.1768 The three elements are shown with their local
EATheat EATcool
R6 0.25 0.1768 displacements in Figure 4-26.
R7 0 0
R8 0.5 0.3536
Chapter 4 Page 88
Matrix methods
F1 EA cos sin u1
F2 L cos sin u2
F1 EATheat EA 0 1 u1
F2 EATheat L 0 1 u2
1 1
F1 EA 2 2 u1
u
F2 L 1 1 2
2 2
and
1 2
F1 EATcool EATheat
4 4
and Figure 4-28 - Internal forces for element 1-3
showing effects of temperature increase and
temperature decrease separately.
Finally, for element 1-3 the relationship is:
1 1
F1 EA Tcool EA 2 2 u1
u
F2 EA Tcool L 1 1 2
2 2
Substituting we find:
1 2
F1 EATcool EATheat
4 4
and, again
F2 F1
Figure 4-29 Internal forces for element 1-2 showing
effects of temperature increase and temperature
These results are shown in Figure 4-27, 4-28 and decrease separately.
4-29. These diagrams should always be the final
step in a matrix analysis to make sure that all 4.11 Tricks of the trade-matrix partitioning
forces balance. Matrix methods are extremely powerful and
useful structural analysis tools. This approach
Chapter 4 Page 89
Matrix methods
u1
Within this matrix we have blocked out four u
matrix grouping. Let‟s re-write the matrix 2
u3
equation for equilibrium symbolically as:
from this u 4
Pa ua K aa K ab ua u5
Kij u6
Rb ub Kba Kbb ub
u 7
u
This matrix is written in “partitioned form.” The 8
partitioning requires the classification of degrees
of freedom. The “a” displacements are free to to a re-ordered set of system displacements:
move (they are caused by the applied load set u5
Pa) while the “b” displacements are restrained to u
be zero (they are associated with the reaction 6
u3
forces, Rb). In the previous example there were
two “a” displacements and four “b” u 4
displacements. The [Kaa] matrix was found by u1
simply crossing out the “b” rows and columns in u 2
the unrestrained stiffness matrix.
u 7
u
When you compare the elements outlined in the 8
stiffness matrix with the elements in the
partitioned equilibrium we can see that the This can be done by forming the following
unrestrained displacements are given as: relationship between the current displacements
(current means the ones we are using now) and
Chapter 4 Page 90
Matrix methods
u1 0 0 0 0 1 0 0 0 u5
u 0 0 0 0 0 1 0 0 u6
2
u3 0 0 1 0 0 0 0 0 u3
u4 0 0 0 1 0 0 0 0 u4
u5 1 0 0 0 0 0 0 0 u1
u6 0 1 0 0 0 0 0 0 u2
u7 0 0 0 0 0 0 1 0 u7
u 0 0 0 0 0 0 0 1 u8
8 Figure 4-30-Six element truss example
Chapter 4 Page 91
Matrix methods
P
1
Ri Kba Kaa j
Chapter 4 Page 92
Matrix methods
Chapter 4 Page 93
Matrix methods
4.13 Chapter 4 - Homework Problems one-half the cross-sectional area of the other
elements.
Problem 4-1
(a) Solve for the unrestrained 3 x 3
system stiffness matrix. This matrix
will have the factor EA1/L1, where
A1 is the area of element 1 and L1 is
the length of element 1.
(b) Solve for the system stiffness matrix
for the 2 x 2 restrained system.
(c) Find the deflections u2 and u3 and
Problem 4-1-Figure the total reaction force at node 1
where the three elements are
a) Develop the system stiffness matrix for the attached to the wall. These will be a
three element bar shown. function of the parameter L1/ EA1.
b) Use the applied loads and boundary (d) Solve for the internal forces in each
condition to develop the restrained system of the 4 elements in terms of L1, E
stiffness matrix with the applied loads and A1.
c) Find the deflections at each node
d) Solve for the internal forces in each rod Problem 4-3
element
e) Draw the applied forces in their correct
directions for each of the three elements
f) Draw the free body diagrams for each of the
three nodes showing the loads on the nodes
and the loads in the elements
g) Compute the strain energy in the structure.
Problem 4-2
Problem 4-2 Figure (a) Solve for the 4 x 4 system stiffness for
this structure.
A configuration has four rod elements connected (b) Solve for the 3 x 3 restrained stiffness
to three node points shown circled. All elements matrix.
have the same Young‟s modulus. Node 1 is (c) A 500 lb. load is placed at node 4; find
fixed (u1 = 0). The nodes move only in the the displacement of node 4 and the
horizontal direction. Node 2 has an applied load internal forces in each of the four
of 1000 lb. Node 2 has only elements 2 and 4 elements.
attached to it. Node 3 has three of the four
elements directly attached to it. Element 4 has Problem 4-4
Two loads are applied to a three element truss.
Chapter 4 Page 94
Matrix methods
a) Develop the 4x4 stiffness matrices for each b) If a single force Po=100 lb. is applied at
of the three elements using the nodal node 1, upward along the y-axis, solve for
numbering provided. Degrees of freedom the deflections.
are numbered odd in the x direction and c) Solve for the element forces in each spring.
even in the y-direction. Be sure to indicate which are in tension or
b) Develop the three transformation matrices compression.
required to expand 4x4 matrices to the three
8x8 matrices required to generate the system
stiffness matrix.
c) Write the three element stiffness matrices in
terms of 8x8 matrices in system coordinates
d) Superimpose (add together) the three
element stiffness matrices to find the
unrestrained system stiffness matrix.
e) Identify which degrees of freedom are
restrained.
f) Apply the loads and solve for the truss
displacements
g) Solve for the nodal reaction forces. Draw
them on a diagram (and check to see if
equilibrium is present).
h) Solve for the internal forces and draw them
on a diagram. Be sure to identify tension Problem 4-5 Figure
and compression.
Problem 4-6
Five elements are used to form the truss shown.
All members have equal areas. Nodes 3 and 4
are fully restrained against movement. The truss
members are made of aluminum (E=10.5 x 106
psi) and have cross-sectional areas equal to 0.5
in.2 Write and use a MATLAB program to do
the following:
Problem 4-5
Four springs are assembled as shown. Nodes 2,
3, 4 and 5 are restrained to be zero. The angle
is 45o. Two different axis systems are shown.
This problem will use the x-y system oriented
horizontally and vertically. Use the global
system numbering shown.
a) Solve for the 2x2 restrained system stiffness
matrix in terms of the x-y system when k1 = Problem 4-6-Figure
2k2.
Chapter 4 Page 95
Matrix methods
Problem 4-7
Develop the restrained stiffness matrix for the
redundant truss shown in the figure in terms of
the two nodal displacements shown.
Problem 4-8-Figure
Problem 4-7-Figure
Chapter 4 Page 96
Matrix methods
Problem 4-9-Figure
Chapter 4 Page 97
Materials and processes
Properties
INTRINSIC ATTRIBUTIVE
Bulk
mechanical Price &
properties availabilty
Figure 5-2-Mechanical properties are important for
aircraft operation (Source:Ashby)
Bulk non- Production
mechanical DESIGN properties thousands or even tens of thousands. How
properties
do we organize materials selection with such
a huge menu of choices? How do we match
Surface Aesthetic
properties properties our choices to respond to system
Chapter 5 Page 98
Materials
requirements? What are the metrics we not be of greater weight per square foot and
must use? must have reasonable resistance to tearing.”
Materials have an important effect on what The fact that this statement seems so
airplanes look like and the challenges of primitive to even the uneducated layman is a
aircraft and spacecraft development tribute to how far teams of materials
determine what materials are developed. engineers and structural engineers have
Figure 5-3 shows the major types of brought our aeronautical structures
materials used today in aerospace design. technology in less than a century.
These materials are largely metallic, but this
was not always so. In an article published in Materials also progress through a series of
1922 this observation was made.24 stages of development as they progress from
a scientific curiosity to a full-blown, reliable
“So long as wood is used, the making of product.
good aeroplanes will be a matter of art
rather than science. In the covering from New materials often make improved
aeroplane wings, nothing has yet been structural concepts feasible, but, as noted in
produced to compete seriously with linen or a National Research Council report,
cotton fabric treated with suitable dope and developing and “characterizing” a new
Chapter 5 Page 99
Materials and processes
repairability must be fully understood and About 80% of commercial aircraft airframes
acceptable. are composed of aluminum. While this is
changing – consider the new Boeing
We can summarize six requirements for a Dreamliner with its extensive use of
material to be considered for use in advanced composite carbon fiber materials –
aerospace products, these are: aluminum holds a prominent place in
airframe construction. Aluminum is the
Stabilized material and material third most abundant metal in the Earth‟s
processes crust. – after iron and silicon – but it usually
Producibility occurs as an oxide so it is expensive to
Characterized mechanical properties extract from ore. A century and a half ago,
Predictability of structural performance aluminum was so rare and expensive that the
Supportability – the ability to repair a French Emperor Napoleon III had a set of
structure in the field and inspect during cutlery made from it and used it while
manufacture and service entertaining the King of Siam. In 1857,
Availability and affordability Charles Dickens wrote this about aluminum:
We can also define “maturity levels” as “Within the course of the last two years… a
follows to assess whether or not the treasure has been divined, unearthed and
materials we are considering are risky for a brought to light ... what do you think of a
particular design. These maturity levels are metal as white as silver, as unalterable as
defined by a numbering scheme as follows: gold, as easily melted as copper, as tough as
iron, which is malleable, ductile, and with
1-2 For a ranking of 1 or 2, materials are not the singular quality of being lighter that
available and need development, no glass? Such a metal does exist and that in
complex parts are built, no test data considerable quantities on the surface of the
available globe.
3-4 For this range of numbers properties are
not reproducible, testing is not complete, no The advantages to be derived from a metal
sub-components or full scale components endowed with such qualities are easy to be
built, no cost model validation with full- understood. Its future place as a raw
scale parts material in all sorts of industrial
5-7 Limited data with large error bars, sub- applications is undoubted, and we may
expect soon to see it, in some shape or other,
components built but testing is incomplete,
full-scale parts built but not tested to design in the hands of the civilized world at large.”
loads, cost model validation is incomplete
8-10 Materials and processes are well- In 1884 aluminum was chosen to be the
characterized, solid data base for different metal for a pyramidal cap to be placed at the
conditions, reproducible sub-components top of the Washington Monument as part of
and full-scale components built and tested to a lightning protection system. At that time,
design loads, cost models validated the price of aluminum was $16 per pound
and the pyramid (whose replica is shown in
The development of aluminum is a prime Figure 5-4), was fabricated using 97.75%
example of a material advancing from a pure aluminum, weighed a little over 6
rarity to a standard material for aerospace pounds (2.85 kg.).
use. Aluminum is the primary material used
to construct aircraft today. There are many The only U.S. producer of aluminum at that
different aluminum alloys to choose from time was William Frishmuth of
today, but it was not always the case. Philadelphia. The total production of
aluminum by his foundry for 1884 was 51
advantages attending it, but from the these loads. Some materials must be
pressure of a world shortage of the limited transparent to light or radar.
supplies of that class of timber most suitable
for light structural purposes”. Materials selection also requires strong
consideration of how the part will be
North was referring to a situation that manufactured. Will the part be cast,
occurred during the latter part of the World stamped or machined? It is not unusual to
War I when a shortage of aircraft quality have a one pound part milled and cut from a
wood seriously threatened to disrupt British 50 pound forging. How do materials
aircraft production. Brought on by Britain's choices affect the size of the part? Can the
reliance on imported supplies of aircraft part be manufactured as a single piece or
woods, one would expect to find a strong must it be assembled? Do we use bolts or
interest on the part of British engineers in do we weld components together?
advancing metal construction. Such a In my opinion, the best book to learn the
movement failed to develop and a period of fundamentals of materials and the selection
British aeronautical stagnation followed process that we should follow is Michael
World War I. British engineers ignored the Ashby‟s textbook.28 I will make frequent
progress in metal construction being made reference to this text and use his PC based
by their counterparts in Europe and in the design code (CES EduPack). This code is
United States, and adopted a conservative accessible from the Purdue Libraries web
and unprogressive approach to design. site and the ECN Network. When you have
the chance, buy a copy of Ashby‟s book for
In 1924, as England moved to revitalize your own use.
procurement of military aircraft, the Air
Ministry ruled that all vital parts of future 5.2 Materials families and classification
service aircraft were to be made of metal. Let‟s start at the beginning. There are 103
Disappointingly short of a positive decision, naturally occurring elements. Everything
which might have promoted development of that ever has been or ever will be is
all-metal aircraft designed to take advantage composed of pieces from the periodic table
of the superior weight/strength properties of of elements. Aerospace materials are sub-
metals, the ruling did little more than foster divided into six different “families” or
continuation of fabric-covered biplanes with categories: 1) metals; 2) polymers; 3)
high alloy steel frames.27 elastomers; 4) ceramics; 5) glasses; and, 6)
hybrids such as advanced composites
While we usually think of choosing obtained by combining two or more
materials on the basis of strength there are different materials. Families have common
other important features of materials that features. These features include mechanical
come into play. properties, applications and processing.
Not all design requirements involve As shown in Figure 5-5, families are sub-
mechanical or thermal loads. The aerospace divided into classes, sub-classes and
operational environment can be harsh; for members. Each family member has
instance, naval aircraft operate in a highly
corrosive salt water environment. In 28
addition, materials must not only resist M. Ashby, H. Shecliff and D. Cebon,
loads, but may be called upon to store, Materials: Engineering, Science, Processing and
Design, Elsevier, 2007. This is a “must have”
transmit, dissipate or diffuse energy from
book. The majority of the figures used in this
chapter are taken from this textbook. We will
27
Charles H. Gibbs-Smith, Aviation, Her also use the CES/Edupack program for
Majesty's Stationery Office, London, 1970. homework.
This list uses terms that need further fighter aircraft. The United States imports
definition and we will do so in the next 100% of its manganese and 78% of its
section. Some terms, such as availability, cobalt and chromium.
are important but are not technically
measurable. Availability is sometimes 5.3 Material properties
affected by the political climate. For Materials scientists and engineers have their
instance, strategic materials are those that own language. Words and phrases must be
may be readily available in some countries used carefully. Unfortunately some of the
but not others. For instance, alloys might terms, like “strain,” “elastic” or “plastic” are
have certain elements such as chromium but also used imprecisely in common everyday
these elements might be in short supply in use to mean something different or less
the United States. precise.
Obtaining these elements during a crisis or Material properties fall into four general
during wartime can be a problem. In categories: 1) mechanical properties;
particular, the "Big 4" chromium, cobalt, thermal properties; electrical, magnetic and
manganese, and platinum are the most optical properties; and chemical properties.
critical. Because of the obvious importance of
strength, weight, deformation and durability
In addition to airframe development, to structural design, mechanical properties
materials development plays an essential are at the top of the list for designers.
role in new engine development. Materials
innovations such as titanium alloy Material properties change with
development, directional and single crystal temperature, a major problem for high speed
blade production, thermal barrier coatings, designs or designs that operate in space.
and composites development have opened Seldom does an increase or severe decrease
the design space to implement new engine in temperature improve a material property.
design concepts. As temperatures increase, materials such as
aluminum will “creep.” They slowly
deform or sag over time. Because of this
there is a maximum service temperature.
Some materials such as stainless steel have a
high service temperature – around 800o C –
while polymers have low service
temperatures – around 150o C.
families of materials. With this this column before it buckles is given by the
representation each family of materials Euler buckling formula:
occupies a particular area of the plot –
metals near the upper right, foams at the
b
2 2 4
lower left, polymers low central, and so on. P Pcrit EI E
L L 12
In terms of logarithms the specific stiffness EA
2 2
Consider the flat plate shown in Figure 5-10. First of all, the bending stress due to the
The load, P, the maximum deflection, max, applied load is:
and the width, b, are specified; the only free
design dimension is the plate thickness, h.
My
FL t 2 FL
6 3
4
The relationship between P and the center I t t
deflection is: 12
t4
EI Ebh3 where I
P 3 max
3 max
12
L 12 L
2
The thickness dimension, h, must be FL 3
The area is A t 6
2
y
1
12P 3
h L
Eb max The weight is w t 2 L
materials; stress is a relatively new concept where A is the area of the test coupon. For
in structural mechanics. Gordon30 credits test data purposes this area is assumed not to
Cauchy, who we generally associate more change as loading proceeds.
with mathematics than structures, with the
origins of the modern concept of stress and Materials data such as this are available in
strain, although Hooke is credited with the Metallic Material Properties
recognizing the linear relationship between Development and Standardization
load and deflection.31 (MMPDS) Handbook. This handbook
replaced another government sponsored
Stress causes strain. Strain is a measure of handbook, MIL-HNDBK-5, as the official
relative deformation; strain is dimensionless. source of materials properties in 2003. The
Normal strain is created by tensile and MMPDS Handbook is the only government
compressive stresses. recognized source of design allowable
materials properties for commercial and
military aircraft certification and
airworthiness. The MMPDS Handbook
contains data such as tensile and ultimate
We can also measure strain in terms of yield stresses for the most widely used
percent strain. In that case we multiply the metallic materials used in aerospace
number above by 100. construction. It also contains information on
fracture toughness and Young‟s modulus.
A stress-strain diagram for metals, such as
the one shown in Figure 5-16 provides an
illustration of many of the important
properties of structural materials. Such
stress-strain diagrams are developed by
loading small test specimens or “coupons”
in tension (or compression) with an axial
load, P, and measuring elongations, , under
loading. These elongations are converted
into strain measurements while the loading
is converted to stress using the relationship
=P/A.
30
J.E. Gordon, “Structures-or Why Things
Don‟t Fall Down.” This book and its
companion, “The New Science of Strong Figure 5-16- Typical stress-strain diagram for metals.
Materials” are “must haves.” Source: Ashby, Shercliff, Cebon
31
Cauchy derived the equations of equilibrium
for in infinitesimal element within the elastic Figure 5-16 shows that when first stretched,
body. He showed that on any plane within an most materials, such as the metals used in
elastic material the stress had three components aerospace construction, have a linear
(two shears and one normal). Timoshenko relationship between stress and strain so that
attributes Cauchy‟s interests in the theory of = E. The initial slope of the stress strain
elasticity to a paper by Navier (of Navier-Stokes curve is the Young’s modulus E and the
equation fame). Cauchy was a renowned stress-strain relationship is said to follow
mathematician who pioneered the study of real
and complex analysis, convergence and
Hooke‟s Law.32 When first loaded, most
divergence of infinite series, differential
32
equations, determinants, probability and Young‟s modulus is named after Thomas
mathematical physics. Young (1773-1829) who has been referred to as
materials are also elastic. Elastic Above the yield stress, when we unload a
deformation means that when a load is material the unloading will occur along a
applied to a material it deforms, but when straight line (shown as a dashed line in
the load is removed, the material will return Figure 5-16) whose slope is given by the
its original form; when an elastic material is modulus of elasticity or Young‟s modulus.
loaded above a stress called the “yield Note that this line is parallel to the original
stress” or “yield strength” and then loading line.
unloaded it will have a permanent
deformation. Above the yield stress the This unloading line will intersect the
material is inelastic or “plastic.” horizontal strain axis at some point. By
agreement among engineers, when this
Plasticity is the ability of a material to straight line intersects the strain of 0.002
undergo permanent changes in shape under (note that strain is dimensionless since it is
stress. When materials become plastic they defined as the change in length divided by
are easier to deform, as indicated by the the original length) we define the yield
reduction in the slope of the stress-strain stress. Thus the yield stress is often called
diagram. the 0.2% “offset yield” stress (this is the
strain times 100). The percentage unit is
Materials can also be nonlinear elastic, but used in the definition and can be confusing,
the key concept of material elasticity is that but the important thing is that this strain is
when stress is removed, the material returns extremely small.
to its original state with no permanent
deformation. Figure 5-17 shows four typical stress-strain
curves for polymers. Two of the four stress-
As indicated in Figure 5-16 ductile materials strain curves are for polymers with extreme
with may have no clear yield point. The plasticity before failure. On the other hand,
material is linear elastic up to a point and one of the polymers is brittle and shows no
then becomes nonlinear elastic for a small plasticity before failure.
range of stress, then becomes plastic.
Plasticity is a desirable feature and is
Yield strength is a transition point between related to ductility. Ductility is the ability of
elasticity and plasticity. Materials such as a material to deform substantially (of the
titanium have high yield strength while order of 10% increase in length) before
others like lead have low yield strength. failure. Ductility is characteristic of metals
like steel or aluminum alloys. At a large
strain, of the order of 10% change in length,
these materials reach their “ultimate” stress
as indicated in Figure 5-16. Further loading
causes fracture.
“England‟s Galileo.” Young developed his idea Ductility has several definitions but “the
without stress and strain as we use it today. He property of being capable of sustaining large
was truly multi-dimensional in his contributions plastic deformations without fracture (able
to science. His development of the double slit to be drawn into a wire)” is probably the
optical experiment helped to establish the wave best.
theory of light. A memorial plaque in
Westminster Abbey describes him as a “man Mechanical properties should not depend on
alike eminent in almost every department of
the size of the test specimen. However,
human learning.” Young‟s biography by
Figure 5-17-Stress-strain
Andrew behavior
Robinson is "The for polymers
Last Man Who Knew
once the material becomes plastic, the cross-
showing strain to
Everything" (Pi Press, 2006).failure sectional area changes substantially; the
stresses are no longer uniformly distributed
pressurized and depressurized half a dozen called the mode 1 stress intensity factor.
times in a day of operation. As a result, a The symbol for this quantity, including the
crack like those shown in Figure 5-19 and parameter Y, is:
Figure 5-21 may form and then grow. When
the crack reaches a critical size it will
suddenly lengthen, creating a breach in the
skin and the skin material will fracture. The The stress intensity factor K1 has units of
ability to resist the growth of cracks is called stress times the square root of length.
fracture toughness and is measured by the
material property K1c. A crack will propagate if the crack length
exceeds a value that will cause the stress
intensity factor to exceed a critical value.
This critical value is the fracture toughness.
Fracture toughness is a material property. It
is independent of how we measure it; it can
be used for design.
m
Figure 5-24- The endurance limit; fatigue strength at dc
107 cycles. (Source ASC) B( )m ( c) 2
dN
Integrating this equation, we find that In Australia, an Ansett Airlines Fokker F-28
turboprop short range airliner was retired
from service in 1988 after 29 years and 9
months of service. The flight time record for
1 1 1 1
N
the F-28 was 70,000 hours in 1989.82
B 2 m 1 m c m 2 1 c m 2 1
m
2 initial critical
Although the age of an airliner is not as
important to structural integrity as the
The design life of an airliner is defined by number of flights and flight hours, age
the aircraft manufacturer after discussions becomes important if the aircraft is
with the airlines. The design life is based operating in a wet or salty environment. This
upon an economic analysis of procurement will accelerate corrosion unless the airframe
and operating costs in combination with is properly protected, inspected and
forecasts of obsolescence with the maintained. The combination of age,
availability of the next generation of more corrosion, and improper maintenance may
economic aircraft. result in a significant reduction in the
damage tolerance of the structure.
After the First World War most airliners
were modified bombers and their life was 5.9 Special problems in aerospace
only a few hundred hours. By the late 1920s, structures
most of the aircraft in airline service were Aerospace structures face most of the
specially designed for airline operations. problems considered by mechanical and
The economical service life of these aircraft civil engineers plus a few others that make
in service with the major airlines was only aerospace structural design challenging.
four to five years and between 5,000 to Aerospace includes not only the design of
10,000 flight hours. airframes for use from low to extremely
high speeds (including re-entry) but also the
By the mid- 1950s with the introduction of design of gas turbines and rockets. Where
turboprop and turbojet airliners, the design these designs involve human passengers and
economic life had been increased from 10 to operators safety requires high degrees of
20 years. Design flight time was increased reliability.
from 10,000 hours to between 60,000 and
140,000 hours, depending on the aircraft 5.10 Homework – Chapter 5
range. The airframes of short range airliners Problem 5-1 – We want to design a
were being designed for more than 100,000 centrifuge disk (flywheel) that can rotate as
flights. fast as possible. The maximum stress
cannot exceed the yield stress of the
By 1997 a large number of the jet airliners centrifuge disk. The disc has radius, R, and
which were produced in the 1950s and thickness, t. The energy stored in the disc is
1960s were still in airline service and had
exceeded the original economic design life U R 4 2t while the disc mass is
4
of 20 years. By the end of 1988 eight
hundred and eight Boeing 707, 720, 727 and m R 2t . The quantity to be maximized
737 aircraft had been in service for more
U 1 2 2 1
is R . Find the material index to E2
m 4 M (neglect self weight). (Figure
aid in the choice of a material for the
centrifuge so that it can rotate at maximum d). Begin by computing the beam mass.
angular velocity. The maximum stress in
1
the disk is given as y 2 R 2 .
2
Problem 5-2 – A panel has a fixed width, b,
and length, L. The thickness is h. The
moment of inertia for calculating the
bending stress is I = bh3/12. The maximum
bending moment applied to the beam is
Mmax. If the beam stress cannot exceed the
yield stress, derive the material index to help
us choose the lightest weight beam material.
1 y
2
1 1 b
U A 2t h b I zz h3t 1 3
2 2 E 6 h
h t
Find the material index that defines a 3 0.10
minimum weight spring. b b
Problem 5-7
The expression for the bending shape factor
is
12 I
B 2
A
Problem 5-8
a) Calculate the bending shape factor B for
2 t
A 4bt I zz b3t 0.125 a thin-wall circular cross-section with outer
3 b radius 40 mm. and inner radius 36 mm.
b) Calculate the cross-sectional area of a
solid square cross-section with the same
value of cross-sectional moment of inertia as
the thin-wall circular tube in part (a).
c) Determine the relative weights of two
beams with equal length L, one with the
thin-wall circular cross-section and the other
with the solid square cross-section. Both
b
A a b t I zz a 3t 1 3 beams are constructed from the same
4 a material.
a t d) Calculate the shape factor for a hollow
2 0.10 square cross-section with b/t=5.
b a e) Calculate the shape factor for an I-beam
with h/b = 3, t/b = 0.125.
f) Which of the cross-sections in parts (d)
and (e) will provide the lightest weight beam
if the beam material used is the same?
Why?
Problem 5-9
The bending stress shape factor is defined
Z 6Z
as Bs 3 where the area and Z
Z square A 2
are those for the particular cross-section.
Problem 5-10
The I-beam cross-section in Problem 5-7
with h/b=4 and t/b=0.125 is to be designed
to carry a bending moment Mo.
a) Solve for the dimension h required
to resist this bending moment
without yielding. Use general
material properties such as y and
to develop an expression for the
weight (per unit length) of the beam
cross-section.
b) Develop the same relationship for
the weight per unit length for a thin-
wall elliptical cross-section with the
same material as part (a) with a/b =
3 and t/b=0.125.
Figure 6-1- An example of beam loads and Beam bending moments are converted to
structural behavior stresses by using the “flexure formula.” The
flexure formula provides the bending stress
(also called a “normal stress”) at a point y
“Strength of materials” deals with structures on the beam cross-section. The y-position is
such as that shown in Figure 6-1, loaded by measured from the neutral axis, an axis that
a 100 pound transverse load that creates passes through the beam cross-section area
support reactions as well as internal shear centroid. For a symmetrical beam cross-
forces and bending moments. The internal section with a bending moment M:
shear force is labeled as V while the internal My
I
M y dA
Figure 6-4 - Normal stress distribution area
pattern over a solid symmetrical cross-section
The minus sign appears because, as shown
In both of these figures the resultant bending in the diagrams in Figure 6- and 6-4,
moment (positive as shown) creates negative positive M produces a compressive
(compressive) normal stress above the (negative) normal stress when the location
neutral axis. is above the cross-section neutral axis.
The development of the seemingly simple Even though any normal stress distribution
flexure formula has a long interesting that produces the resultant bending moment
history. In Medieval times timbers were M is theoretically possible, the stress
used for bridge and building construction; it distribution must produce compatible
was important to understand how large a deformations. Compatible means that when
timber member should be to carry loads. we consider how the beam deforms under
stress we have to have something that makes
Galileo began the beam theory development, sense and is close to that observed in the
but used some basic wrong assumptions that actual system. As a result, there is only be
gave him the wrong answer. DaVinci also one stress distribution that satisfies
tried his hand at predicting beam stresses but equilibrium and produces the correct
was also was unsuccessful. He was deformation.
followed by Jacob Bernoulli, who assumed
Early researchers noticed that a beam axis for the flexure formula must be located
subjected only to bending moments at the cross-sectional area centroid. The
deformed into an arc of a circle, as shown in “neutral” surface with zero stress is located
Figure 6-2. This led to the classical at the cross-sectional area centroid. The
assumption that “plane sections remain neutral axis (y=0) is a transition point
plane.” between tensile and compressive stress.
area area
Figure 6-2– “Plane sections remain plane” is moment of inertia” and is written as
I NA y dA since it is computed about
a fundamental assumption about the 2
deflection of beam cross-sections loaded by area
bending moments. Bending produces stresses
the cross-section neutral axis.
but no longitudinal resultant forces
The deformation pattern in which cross- The expression for bending moment now
sections rotate but do not otherwise distort becomes:
means that the strain on internal cross-
sections varies as follows: M E I NA
P dA E y dA 0
area area
E y dA 0
area
The condition that bending moments create Figure 6-3-Example T-section used to
no in-plane force resultant means that the y- illustrate how we find the cross-sectional
neutral axis and area moment of inertia
3*13
I total 3* 3.5 2.5
2
12
1*33
3* 1.5 2.5
2
12
0.25 3 2.25 3 8.5 in 4
My
P dA
area area
I
dA
the centroids of the three areas with respect This integral is computed over the area of
to a temporary axis with an origin at the element 1.
base of Figure 6-6.
M
The centroid of the assembled area is P1
I y dA 0
Y 0.476 inches above the section base. Y area 1
in column 4. The contribution of the The resultant forces on each of the three
parallel axis term is shown in column 7. elemental areas are computed from these
These terms are summed to create the total three formulas:
area moment of inertia of the entire cross- M
section about its own area centroid. P1 Q1
I
M
Izz = 0.04623 in.4 P2 Q2
I
Now we can use the flexure formula to find M
P3 Q3
the bending stress. Let‟s look at a load I
carrying feature of each of the three cross-
My 2 M
Mi
area i
ydA
area
I
dA
I
area i
y 2 dA
y dA
Figure 6-10- A side view of a beam cross- 2
section FBD showing force resultants and the
Mi Ii
Ai
bending moment in force and moment
equilibrium M I I
section. The sum of the three forces equals The fraction of the bending moment resisted
zero and their net moment (the sum of the by each cross-section elemental area
products of the three elemental forces times moment of inertia (computed with respect to
their distances from the neutral axis) is equal the sectional neutral axis) is proportional to
to the bending moment, 1000 in-lb. the elemental area moment of inertia
Elements lying further away from the (computed with respect to the sectional
neutral axis (and having large areas) carry neutral axis), divided by the centroidal area
the highest loads. moment of inertia of the entire section.
Figure 6-11 shows the same side view of the For the cross-section in Figure 6-6 we get
beam element, but with all stresses shown in the following results:
terms of resultant normal forces on each of
the two beam “faces.” On both the left and
The beam cross-section is subdivided into From Figure 6-13 moment equilibrium
three elemental areas to expose internal requires that the bending moments are
shear forces. The shear force resultants are different on the left and right faces when a
shown as horizontal vectors in Figure 6-12 shear force V is present. Summing moments
acting parallel to the horizontal faces. about a point on the left “face” of the
element (in the clockwise direction), we
have
moments Va M M M 0
so M Va
Va Va
P1 ydA y1 A1
I area 1 I
Figure 6-14 – Bending stress distributions on left Va
and right faces of a beam element located a Q1
I
distance ”a” from each other when a shear force
Va
resultant V is present. Note that directions are P2 Q2
shown with positive values. I
Va
Now, subtract the distribution My/I from P3 Q3
I
each side of Figure 6-14 and redraw the
FBD showing only the “net stresses” = If the shear force V=100 lb. in Figure 6-1 is
My/I. This distribution is shown in applied, the net forces on each cross-
Figure 6-15. sectional element are:
The resultant forces in the horizontal
P1 87.8* a lb.
direction are computed by integrating the
normal stresses over the cross-sectional P2 4.23* a lb.
areas on the right hand or “back” face. P3 92.1* a lb.
Vay
dA
P These net forces are shown in Figure 6-16;
dA
area area
I their sum is zero. They act at the area
centroids of the individual area elements so
If the area of integration is over the entire the P‟s and the shear forces that create
cross-sectional face area then the resultant them are in moment equilibrium.
force P is zero. However, if we simply
compute this integral over any one of the Figure 6-17 shows the Free Body Diagrams
elemental areas then we get a non-zero we need to find the horizontal shear forces,
result. The results for each of the three labeled as F32 and F21.
elemental areas in Figure 6-9 (p. 105) are as
follows. To compute F32 we use the equation of
horizontal force equilibrium on the upper
area element, element 3 in Figure 6-6 (p.
105).
F32 87.8* a 0
Neither the shear flows nor the shear This means that 1 2 .
stresses are functions of the arbitrary
dimension “a.” This is an important result because once we
have found the shear stress on one surface of
We have computed the shear stress on a small element we know it on the other
longitudinal (horizontal) surfaces, three surfaces. This is called the “Principle
perpendicular to the front and back faces on of Complementary Shear.”
which the shear resultants act, but we want
the shear stresses on the left and right Figure 6-19 shows shear stresses at the point
vertical faces. Do we need an additional between the upper element and the middle
computation? The answer is “No.” element.
Figure 6-18 – A very small element with Figure 6-19 – The Principle of
shear stresses acting on all four sides. Complementary Shear – all shear stresses
on an infinitesimally small element must be
equal.
Figure 6-18 shows a small structural element
with dimensions x by y by z with shear The Principle of Complementary Shear
stresses acting on four sides. The allows us to use FBD‟s to compute the shear
dimensions x, y and z are extremely small so stresses and shear flows on longitudinal
that the shear stresses are constant over each sections and then conclude that the shear
surface. stresses are identical on the front and rear
faces of the beam section.
First of all, the two vertical shear stresses
(2) act on surfaces with the same areas Remember that the shear stress changes
(A2=zy) so they must be equal and opposite when we move upward or down ward on the
or this element would not be in static
equilibrium. The same holds true of the
shear stresses in the horizontal direction (1)
since they act on identical size surfaces
(A1=zx).
cross-section. These results are shown in The term Q is the first area moment of the
Figure 6-20. The shear stress changes its area above the point at which the cross-
value depending on where the cross-section section is cut to expose the shear stress, that
is cut. At the lower junction of element 2 is, where we wish to compute the shear
the shear stress is 921 psi. Note that the stress.
shear stresses are in the same direction as
the resultant shear force V=100 lb. that
created them.
y VaQ
P dA Va I dA
area area
I
Q2 upper = +5 in3.
Figure 6-22 – Example problem
The thickness dimension for this cut is t =
First, cut the cross-section at the upper glue 1.667” so that
joint. The thickness of this cut is t = 1.”
The (red) area above the glue joint is A1 = VQ2 1000*5
62.64 psi
1x4 = 4 in4. The centroidal coordinate of It I *1.33
this area is 2.5 inches above the neutral axis.
The value of Q1 is 2.5x4 = 10 in3. Using the Using the Principle of Complementary
simple formula for the shear stress we find Shear, this stress also occurs in the lower
that left longitudinal glue joint. A similar
VQ 1000*10 computation for the lower right glue joint
upper 1 208.8 psi
It I *1 will give the same answer.
To find the shear stress in the lower glue Consider one final example. Compute the
joint we cut the joint, being careful to cut shear stresses in the center element of an I-
away only the lower left element. We have beam cross-section shown in Figure 6-23.
two choices for computing Q. First of all,
our formula was derived considering the A shear resultant of 40,000 lb. acts on this
area above the cut. What do we mean cross-section. The moment of inertia for
“above the cut?” this section about its neutral axis is 43.3 in4.
We will compute the value of Q at a general
There are three areas in the cross-section position located a distance from the
above the cut for the lower left glue joint. In
general, the Q for areas “above the cut” will
be positive. The calculation for the value of
Q then seems to require consideration of
three elemental areas. On the other hand we
can simplify this calculation if we remember
that the sum of all of the Q‟s for the cross-
section is zero.
is the sum of the contributions from the two carrying passengers, cargo and fuel.
areas above the point y=: Suppose that a portion of the wing "box" is
made of aluminum alloy 7075-T6 and has
Q y1 A1 y2 A2 the hollow, rectangular, doubly-symmetric
cross section shown (not to scale):
2
Q 1 2 1 3 2.5
2 A bending moment Mz = 2000 kip-in is
2 applied to this beam.
Q 9.5
2 (-2<<2)
Calculate the maximum tensile and
The shear stress is: compressive bending normal stresses.
(Answer: 51,700 psi) The theory leading to
2
9.5 the flexure formula is not valid if the
VQ 40,000 2
predicted stresses exceed the material's yield
It 43.3 1 strength. Is the theory valid?
8780 462.5 2 Problem 6-2:
(-2<<2)
Problem 6-2
Problem 6-1
Problem 6-4: The symmetrical beam Problem 6-6: The symmetrical beam cross-
section shown has a total area moment of section shown is composed of three different
inertia Izz (about the horizontal neutral axis) types of element shapes. Element 1 has
equal to 1.504 x 109 mm.4 ( Do not dimensions 8x2 inches; element 2 has
recalculate this number) A vertical shear dimensions 2x2 inches; element 3 has
force resultant of 6 kN is applied. dimensions 20x1 inches.
**
Section 6.4 is adapted from notes written by
Professor W. L. Hallauer, U.S. Air Force
Academy.
Iyy is the area moment about the y axis Figure 6-27 – Example 6-1 – rectangular cross-
because z in the integral is distance of a section showing integration areas for the three
moments of inertia.
point away from the y axis.
For Izz we use the area defined in the upper
In general, a cross-section has three area
left figure to get a familiar formula for the
second moments, defined as:
area moment of inertia of a rectangular area.
I yy z 2 dA h/2
A I zz y 2 dA y 2 bdy
h/2
I zz y 2 dA
A
h/2
h/2 y3
b y dy b
A 2
I YY Z2 dA
I yz yzdA yzdydz A
A A I ZZ Y 2 dA
b /2 h /2 A
b /2 h /2
ydyzdz I YZ YZdA
A
b /2
h /2 y
2
h /2
zdz 0
2 b /2
Substituting the relationships between the
centroidal coordinates (Yo,Zo) and the
centroidal axis coordinates (y,z) we have:
For beam analysis we require the 2nd area
moments relative to the "global" centroid of I YY (Zo z ) 2 dA
the beam cross section. For a simple section A
Zo 2 dA Zo 2 A
A
I YY I yy Z o A
2
I YZ I yz Yo Z o A
Y
2
The integral dAi is the 2nd area
Ai
moment of subsection area Ai about the
cross-sectional or “global” centroid C,
which is generally not the centroid of Figure 6-30 – Example 6-3 beam cross-section
subsection area Ai , so we apply the Parallel
Axis Theorem: First determine the area centroid of the
cross-section. This position is shown in
Figure 6-31. We then apply the general
Y dAi I ZZi I zzi Yi Ai
2 2
parallel axis formulas.
Ai
Then our equation becomes
n
I ZZ I zzi Yi Ai
2
i 1
n
I YY I yyi Z i Ai
2
i 1
Figure 6-31 – Example 6-3 elements showing
n
I ZZ I zzi Yi Ai
2 cross-sectional centroid location
i 1
Starting with the z-axis we write:
I YZ I yzi Yi Z i Ai
n
i 1 I zz I yyi z
i
2
Ai I yy1 z12 A1 I yy2
z2 2 A 2 I yy3 z32 A3
It is essential that you understand that
Yi and Z i are the coordinates of the ith local This gives us the following result:
centroid relative to the global centroid and (3)(1)3
that Ai is the cross-sectional area of the ith I zz (3)(1)(-2.36) 2
12
sub-section
(1)(3)3 (5)(1)3
(1)(3)(-0.36) 2
6.4.3 Example 6-3 12 12
(5)(1)(1.64) 2
Find the second area moments of the
unsymmetrical beam cross-section shown in 33.46 in 4
Figure 6-30. To begin we break this cross- The area moment of inertia about the y-axis
section into three elemental areas. is found by using the following relationship.
5"
I yy I yyi Z i
2
Ai 1"
y
I yy1 I yy2
2 .5 8 "
Z12 A1
z
Z2 2 A 2 I yy3 Z32 A3 4"
3 .9 2 "
or 1"
(1)(3)3 (3)(1)3
I yy + (1)(3)(0.18) 2 + 1"
12 12
3 3"
(1)(5)
+ (3)(1)(1.18) 2 + + (1)(5)(-0.82) 2 Problem 6-5
12
= 20.56 in 4
I yz I yzi Y Z A i i i
6.4.4 PROBLEMS
6-4. Calculate Iyy, Izz, and Iyz for the cross
section shown. Answers: 207, 535, and 0 2"
in4
6"
2"
y
Z 1"
6"
z
Problem 6-6
2" 4 .2 7 "
2"
10"
Problem 6-4
7.1 - Introduction
In this section, we consider beam-like Bechereau's Deperdussin was a streamlined,
aircraft structures, and the simplified externally braced, mid-wing monoplane,
idealization of those structures to permit with a shell-like or “monocoque” fuselage of
approximate structural analysis for molded plywood. Contemporary aircraft
conceptual and preliminary design. We will structures used the skin only as a covering
also discuss the main characteristics of and relied on trusses and frames to carry the
reinforced, thin-walled structural design primary loads; monocoque construction
applicable to aircraft with high-aspect-ratio, derived its strength solely from the load-
relatively aerodynamic surfaces. carrying capacity of the skin.
Redistribute stresses around structural (2) The skin panels resist only shear stress;
discontinuities, such as skin cutouts this shear stress is assumed to be constant
(windows and bomb bay doors for across the thin wall. The shear stress is also
instance). constant along the arc length of a skin panel
Provide edge restraint for skin panels, to between stringers. Constancy of shear stress
increase buckling stress. along the skin panel arc length is equivalent
Provide end seals for pressurized to the wall thickness being so small that the
fuselages. axial force due to normal stress in the skin
panel due to bending is negligible.
7.2 The skin-stringer idealization -
fundamental assumptions A third assumption is that the transverse
Modern computers did not exist when semi- members are rigid within their own planes.
monocoque design evolved after World War The cross-section shape is unchanged during
I, so engineers developed a relatively simple loading. This condition is usually enforced
method of structural analysis that was by ribs and bulkheads.
appropriate for preliminary numerical
evaluation early in a design program. When the cross-sectional dimensions of
longitudinals are very small compared to the
cross-sectional dimensions of the entire
semi-monocoque beam, assumption (1)
above produces little error. However, skin
panels in such a beam may actually resist
significant normal stresses as well as
shearing stresses, so assumption (2) above
may not be reasonable. Therefore, it is
necessary to construct an idealized math
Figure 7-11 - Idealized semi-monocoque model of the structure which recognizes the
structures turn shells into beams normal-stress (bending stress) carrying
capability of some skin elements.
For this method, the structure is idealized, as
shown in Figure 7-11, to consist of two This is done by combining the areas of the
principal types of load carrying elements: 1) skin panels adjacent to a longitudinal
thin webs and skin, called skin panels; and, element with the area of the longitudinal
2) longitudinal stiffeners, called stringers. into a total effective stringer area. One
simple and common method of finding this
Skin and stringer elements resist many quantity is use of the one-half and one-sixth
different types of loads, but the main stress- “rules” that we will develop in the next
carrying purpose of each of these elements section.
is different. As a result, the skin-stringer
idealized models assume that each structural The fact that cross-sectional dimensions of
type carries its own type of stress most longitudinals are small compared with
exclusively. Specifically, this means that: those of the entire semi-monocoque beam
makes it possible to assume without serious
(1) The stringers will resist only normal error that the area of the effective stringer is
stresses such as those from bending. In concentrated at a point on the midline of the
addition, we assume that the load carrying skin where it joins the longitudinal. The
capability is “concentrated” at a point in locations of these idealized stringers are
space. As a result the stress is constant over indicated by small circles as shown in
the cross-section of a stringer. Figure 7-12.
The main objective in developing a skin- Figure 7-13 – Wing box cross-section and
stringer model is usually to estimate equivalent skin-stringer model (ignoring the
accurately the area moments of inertia so stringers)
important to computing bending stress. In
the last section we found that these area The most concern is accurately modeling the
moments of inertia were written as: moment of inertia about the z-axis since this
moment of inertia determines the bending
stresses due to airloads. The one-half and
I zz (I zz i Ai Yi )
2
one-sixth rules are intended to preserve as
closely as possible the true Izz in a skin-
I yy (I yyi Ai Z i )
2
stringer model (the two models ought to
I yz (I yzi Ai Yi Z i ) have approximately the same resistance to
bending even though we lump some areas
together).
Here the area elements represent the
contributions of stringers and spar caps to In the skin stringer model we do not care
the moments of inertia. what the actual shape of the longitudinal
member is. This means that areas such as
The 2nd area moments and products of the four corner areas shown in Figure 7-13
longitudinal members about their own are simply idealized circular cross-sections
elemental centroids are very small compared with no elemental moments of inertia about
to the other contributions in the above their own area centroids.
equations, so they are neglected. Therefore,
the appropriate equations to use for any For vertical thin web elements we develop
skin-stringer model of a cross section are: the one-half and one-sixth rules for a section
such as the wing box cross section drawn in
Izz A y i i
2
I yy A z i i
2
Figure 7-13.
I yz A y z i i i
First, let's develop the equivalent stringer
area A that preserves the Izz of the vertical
web shown in Figure 7-14 so it can get
1 2 For A3 we have:
8"
3 4 A3 = 0.75 + 1/2(.05)(48) + 1/6(.05)(8)
= 2.017 in2
Figure 7-17 - Example 7-1 cross-section Area A4 is:
The actual areas and skin thickness are A4 = 0.75 + 1/2(.05)(48) + 1/6(.05)(8)
2 2 = 2.017 in2
A1,2= 0.5 in A3,4= 0.75 in t = 0.05 in
The area centroid of this idealization is:
a) Idealize the cross section as a skin -
stringer model.
Yi Ai 8 A3 8 A4 Mz
My
Yo x
Ai A y z
I zz
I yy
8(2.017) 8(2.017)
500,000 100,000
1.767 1.767 2.017 2.017 0 (4.26) (24)
32.27 120.8 4,359
4.26"
7.568 x 171,000 psi
(I Ai yi 2 0 1.767 4.26
2
Izz zzi
20" 30" 20"
1 .7 5 0 "
The front spar caps, shown in figure 7-20, The stringers shown in Figure 7-22 have the
have the following dimensions and following dimensions and properties:
properties: 2 4
Area = 0.482 in Izz = 0.247 in
Total area = 0.502 in2 4 4
Iyy = 0.027 in Iyz = -0.055 in
Izz about its own centroidal axis Now we‟ll model the beam cross-section as
= 0.124 in4 eight idealized elements (stringers) that
carry only normal stress and skin panels that
Iyy = 0.070 in4 carry only shear stress. This idealized
model is shown in Figure 7-23.
Iyz = 0 y
0 .8 7 5 " 2 3 4 5
10 "
y 1 8 7 6
z
20" 30" 20"
1 .5 0 0 " z
Figure 7-23 – Skin stringer model for Example
0 .1 5 6 " 7-2
A3 = 0.458 + 1/2(.050)(50)
2
0 .6 0 " = 1.708 in
A4 = 0.458 + 1/2(.050)(50)
0 .1 5 6 " 2
y = 1.708 in
2 .0 "
z
A5 = 0.482 + 1/2(.050)(20)
0 .9 4 0 " 2
+ 1/6(.125)(10) = 1.190 in
A7 = 0.458 + 1/2(.032)(50)
2
I yy A z i i
2
= 1.258 in
A8 = 0.458 + 1/2(.032)(50)
I yz A i y i zi
2
= 1.258 in 2 2
Izz = 0.955(-5.62) + 1.135(4.38)
The sum of the areas is: 2 2
+ 1.708(4.38) + 1.708(4.38)
2 2
Atotal= A1+ A2+ . . . + A8 = 10.22 in2 + 1.190(4.38) + 1.010(-5.62)
2 2
+ 1.258(-5.62) + 1.258(-5.62)
The final results are shown in Figure 7-24 4
= 252 in
Now find the centroid of the idealized skin-
stringer cross section: 2 2
Iyy = 0.955(35.4) + 1.135(35.4)
2 2
y + 1.708(15.4) + 1.708(-14.6)
1 .1 3 5 1 .7 0 8 1 .7 0 8 1 .1 9 0
2 3 4 5 2 2
z + 1.190(-34.6) + 1.010(-34.6)
10 "
5 .6 2 " 2 2
1 8 7 6 + 1.258(-14.6) + 1.258(15.4)
0 .9 5 5 1 .2 5 8 1 .2 5 8 1 .0 1 0 4
= 6,590 in
20" 30" 20"
1 2 y 3
y z
5"
z
4 5 6
5 .4 3 "
Consider the example shown in Figure 7-26. Note that the sum of these three areas is
The cross-section consists of two webs and zero. Substituting the appropriate values of
three stringer elements loaded with a 10,000 Q into the formula we find that:
lb. shear resultant. We want to find the
shear flows in the two webs. VQ1 10, 000
q1 2 (7) 648.15 lb / in.
I zz 108
The first step is to compute the moment of
inertia about the cross-sectional neutral axis. V Q1 Q2 10, 000
q23 (7 3)
First set up a temporary coordinate at the I zz 108
bottom of the section and compute the
925.93 lb / in.
position of the area centroid. This is found
0*2 8*1 12*1
to be y or y 5 in. While these are the correct answers, we have
2 11 no clue about their direction. Let‟s try an
The moment of inertia Izz about the neutral alternative, better way of computing the
axis passing through the area centroid of the shear flows.
stringer areas is:
First of all, in Chapter 6 we cut a small
I zz 2* 5 1*3 1*7 horizontal section of length “a” from the
2 2 2
P1 P2
q2 925.9 lb / in.
a a
dFx q * dx
dFy q * dy
The resultant force Fx due to the shear flow
q, is
Fx qdx q dx qxo
and Figure 7-32-Example problem
Fy qdy q dy qyo
The line of action of the resultant force is
parallel to a line joining the web ends, but is
The dimensions xo and yo are lengths shown not necessarily colinear with this line. The
in Figure 7-30. The distance xo is the direction or sense of the resultant is
projected line length on the x-axis, determined from the shear flow direction.
constructed by joining the ends of the web.
A similar interpretation is valid for yo. For the shear flows q1 and q2 we have the
Notice that xo and yo are independent of the two different lengths:
shape of the web. xo 5"
For the shear flow q1
yo 0
xo 0
For the shear flow q2
yo 4"
This gives the following resultant forces The area A is the area formed by joining the
ends of the web to the point “o” with
Fx1 5q1 Fy1 0 straight lines, then
Fx 2 0 Fy 2 4q2
M o 2qA
While this appears to be a very messy Figure 7-34-Shear flow resultant moments
problem of analytic geometry, it is very depend on the point we are summing moments
simple, since, from geometry, a triangle of about.
perpendicular height r and base ds is formed
by joining point o with the ends of the For point “a” we have Aa=5*10=50 in2.
element ds.
r * ds M a 2*100*50 10,000 in lb.
The area of this triangle is dA while (clockwise)
2
the moment of the elemental force q*ds
For point “b” we have Ab= ½ 5*10 =25 in2.
about point “o” is
M b 2*100* 25 5000 in lb.
q * ds * r qr * ds 2q * dA (clockwise)
The resultant moment of q taken point “o”is The direction of Ma and Mb are determined
found by integration. by inspection.
Figure 7-35 – Closed cell section with constant The shear flow has no net resultant force and
shear flow is therefore the result of a force couple or
twisting moment. If point a lies outside the
Answer: This section is closed; the cell, some of the areas dA would be negative
perimeter is continuous. A torque causes no because qr*ds would be opposite in
net forces in stringers so the shear flow will direction or counter-clockwise for some
be constant around the perimeter, even if points on the web. The algebraic sum of all
there are stringers present (there are none these areas would still be A, the area
here). enclosed by the cell.
The area “A” swept out by integration Consider an additional example. A 100-
involves the total enclosed area of the web, inch-long, thin-walled hollow beam is fixed
50 square inches; it does not matter where to a wall. The cross section at the free end is
we choose to measure this moment. shown in the figure (all dimensions are
measured at the midline).
Aenclosed=50 in2.
Mo = 2*q*50=100q in-lb.
T=2qA=100q
Two 600-lb forces produce a torsional
This means that
couple (torque). The left side (semi-circle)
of the cross section has thickness t = 0.05"
q=T/2A=T/100
while the right side (three straight members)
has thickness t = 0.1 in. We will calculate
To summarize, the application of a torque to the shear flow and shear stresses in both
a closed section creates a shear flow that is
sides (neglecting possible stress
constant around the perimeter. The concentrations).
relationship between the shear flow and the
resultant torque is
First of all, the torque T is:
T=2qAenclosed
T=600*25=15,000 in-lb.
T 15, 000
q 53.85 lb / in
2 Aenclosed 2*139.27
Problem 7-4
The shear stress in the semi-circular section Problem 7-5: An idealized fuselage cross-
is section indicated in the figure is
q 53.85 symmetrical about both the y and z-axes.
1077 psi
t 0.05 The skin thickness is not given but is
assumed to be so small that we can ignore
The shear stress in the rectangular section is: it in our computation of the moments of
inertia and Q. The so-called flange areas
area also so small that their moments of
q 53.85
538.5 psi inertia about their own area centroids can
t 0.1 be neglected.
Observation: If a thin wall cross-section a) Compute both Izz and Iyy about the
carrying constant shear flow q is closed, the section centroidal axes. Data for the
shear flow q has a resultant twisting moment component areas is given in the table
T = 2qA about every point in the plane, below.
regardless of its position. The resultant
force is zero. Element Area z(ft.) y
(in2) (ft.)
7.7.1 –Homework problems A 1 0 4
Problem 7-3 – An unsymmetrical wing B&J 1.5 1,-1 3
cross-section has a shear resultant due to lift C&I 2 1.5,-1.5 1
of L=4500 lb. applied in the location D&H 2 1.5,-1.5 -1
indicated in the sketch. Compute the shear E&G 1.5 1,-1 -3
flow in the semi-circular leading edge F 1 0 -4
(radius = 3 inches).
b) A 1000 lb. downward shear force
resultant is applied along the y-axis
as the result of elevator and tail trim.
Compute the shear flows (and their
directions) in each of webs
connecting the flange elements.
Assume that the thin skin carries
only shear stress and no bending
stress. Also make use of the fact that
Problem 7-3
this section has two axes of
Problem 7-7
Problem 7-5
Figure 7-36 - Open section resultant shear
Problem 7-6: The cross-section idealization force and shear flows produce the same result.
shown is symmetrical about a vertical axis.
The five upper stringers have areas equal to Our problem is to find the shear flows in the
0.4 in.2 The five bottom stringers have webs of the cross-section in Figure 7-36.
areas of 0.8 in.2 The four stringers have equal areas of 0.50
in2. The webs carry no bending stress, only
shear. We also wish to locate the shear
center for the cross-section. Note that this
cross-section will resist bending in two
directions, but has only one axis of
symmetry.
I zz 4(0.5)(52 ) 50 in4
Problem 7-6
Problem7-7 - Find the y-z coordinates of The net forces on the right hand face of each
the shear center of the idealized open skin- stringer is:
stringer cross-section shown in the figure.
This y-z system origin is located mid-way VQi 10,000
Pi a Qi a 200Q1a
between the top and bottom stringers and I zz 50
mid-way between the left and right hand
stringers, but the y-z system is not where P > 0 implies compression on the
centroidal. back face. The values of Qi are either 2.5
in3. (for the two top elements) or -2.5 in3 (for shear force V must be exactly the same
the two bottom elements). This means that as the moment due to the shear flows.
for the top two elements This must be true because the two
representations (total resultant or shear
Pi 500a flow) are equivalent ways of
representing the same cross-sectional
For the bottom two elements stress system. We know the positions
and directions of the shear flows,
Pi 500a therefore we can compute their resultant
moment, MS about point o in Figure
7-37.
The negative sign indicates that the net
forces on the bottom face are tensile.
M S (500 lb / in)(4 in.)(10 in.)
As before, the shear flow analysis proceeds
by cutting each web to isolate it and then M S 20,000 in lb (clockwise)
using longitudinal static equilibrium to
determine the shear flows. These shear 2. We know that the resultant shear force,
flows are shown in Figure 7-37; their correct V=10,000 lb. must be acting through the
direction is important to the computations shear center. The distance “e” shown in
that follow. Figure 7-37 measures the distance from
the point 0 to the shear center. It is an
unknown. The clockwise moment,
about 0, due to shear force V=10,000 lb
is MV = 10,000e
MV = M S
The shear flows in this cross-section have Thus the shear center is located somewhere
been calculated assuming that they were the along a line 2 in. to the left of point o.
result of a cross-sectional shear resultant that Remember these three points about shear
produces only bending stress, but not center computations.
torsion. By definition, the shear center is
that point where a shear force causes 1. The shear center is a point in space and,
bending deformation, but no torsional in general, two calculations must be
deformation or shear stress. The shear flows performed: 1) apply an arbitrary shear
we have calculated therefore assume that the force to a cross-section and compute the
shear force is applied at the shear center. shear flows; 2) place the shear force at a
reference point in space located a
To find the shear center location we use the distance “e” from that point. Equate the
following procedure: moment due to the shear force to the
moment due to the shear flows.
1. At an arbitrary point in space (of our 2. The magnitude of the shear force V does
choosing), the moment (measured in the not make any difference in the final
plane of the page) due to the resultant result. If we had chosen V = 20,000 lb.
the magnitude of MV would have Find, e, the position of the shear center
doubled, but MS would also have with respect to the lower left stringer if
doubled, because the shear flows are b=2 inches.
linear functions of V and likewise would
have doubled. Find the dimension b required if the
3. The shear center lies along an axis of shear center distance is e=0 so that the
symmetry. Some shear centers may be shear center is located at the center web.
found from inspection or a combination
of inspection and computation. Problem 7-10
A thin-wall open cross-section is
Problem 7-8: The open, symmetrical thin- constructed with 6 one-inch stringers and 5
wall idealized cross-section is loaded with a thin webs. An upward shear resultant of
shear force resultant at the shear center, as 1000 lb. is applied. Find the shear flows in
indicated. Find the shear flows and the each of the 5 webs. Why is the shear flow in
position of the shear center. Two the angled webs zero?
coordinates are required.
Problem 7-10
force with a line of action located 8 in. to the point. To compute the moment due to the
left of the right hand corner. shear flows we use the relationship
I zz 4(0.5)(52 ) 2(1.0)(52 )
The areas in this relationship refer to the
100 in 4 area created by beginning with the reference
To begin the analysis, let's cut the webs on point, drawing a straight line to the
each side of the stringers as shown in Figure beginning of the web, tracing the web and
7-39. Call the shear flow in the leading edge then drawing a straight line to connect the
skin qo. end of the web to the reference point. Three
of these webs (webs 3, 4, 5) have zero areas
The net force on the back face of each associated with them.
stringer is
VQ 10,000
Pi i a Qi a 100Q1a
I zz 100
To compute torsional moments we choose Now, the shear flows in each web are
the lower right stringer as the reference calculated to be:
1 P 2 L3
2 2
1 P2 1 P 1 P Lh
0 L x dx
L 1
U U L ds
2
2 EI 2 EI 3 2 h Gt 2 h Gt
1
2
This simple example shows that the shear
Uo flexibility introduced by a thin web can
2 G
substantially increase beam deflection. In
For a single web such as that shown in the this case we see that an important geometric
variable influencing this increase is the ratio
example configuration, the shear stress is
constant across the web and along the web. h/t.
The total shear strain energy is:
Thin-wall web shear is also important for
1
2
the calculation of twisting deflection for
U U o dxdydz tLds
beam 2 G closed thin-wall beams. Let‟s compute a
useful relationship for thin-wall twist by
1 t 1 q
2 2
At any point on the cross-section, the shear The shear flows in the second integral are
flow will be qo +qi. The shear flow qo = proportional to the resultant shear force, V,
To/2A and the shear flow qi is the shear flow so the expression can be written as:
in the web due to the shear force resultant. To
The strain energy expression per unit length qi
1
ds
2A
is:
2 A Gt
qo qi
2
1
U' ds T 1
o2
V qi
2 Gt ds ds
2 4 A Gt 2 A Gt
To
qi where
1 2A ds
2 Gt qi
qi
V
300 60 60 60
360 Gal 0.04 360 Gst 0.07
or
GJ eff 3.05x1011 lb in.2
q0 1190 lb / in.
So that
q1 qo 1200 10 lb / in.
q2 qo 1600 410 lb / in.
q3 qo 1200 10 lb / in.
Figure 7-46
The free body diagrams of the four stringers Notice that the shear flows are negative,
produce the following equilibrium meaning that their directions are reversed
relationships for the shear flows. from our original assumptions. We‟ll leave
these in this form for now.
q1 qo 1200
Now we can solve for the position of the
q2 qo 1600 16,000 lb. shear resultant that creates these
q3 qo 1200 shear flows.
q0 qo
If the position e were specified then we
would simply use the moment equivalence
relationship to compute the 4th equation to
solve for the shear flows. Instead we will
use the twist rate equation:
1 N panels qi si
2 Aenclosed i 1 Gi ti
Figure 7-47
Since we have assumed that the 16,000 lb.
load is applied at the shear center then we Summing moments due to the shear flows
should have no twist angle or twist rate. about the lower right stringer we have:
This means that
M S 2 *1190 * Ao 2 * (10) * A1
1 N panels qi si 2 *1190 * 200 2 *10 * 200
0 472,000 in lb.
2 Aenclosed i 1 Gi ti
4 panels
qi si
or
i 1 Gi ti
0 The 16,000 lb. shear resultant has a moment
given as:
MV 16000* e
Using the three relationships for the shear
flows we have: Equating the two expressions we have:
Wbh
I and
8 allowable
E Wbh
EI
allowable 8
W E h
bc
8 allowable c
modulus, G, is the same for all three with correct directions and positive
webs. magnitudes.
b) Find the torsional stiffness of the
cross-section when the shear d) Compute the twist rate due to this 1000
modulus for the upper web is 1.5 lb. load.
times that of the other two webs.
7-14
c) Solve for the twist rate when the Compute the shear flows when the 1000 lb.
1000 lb. shear resultant is applied 6 force is applied through the shear center.
inches to the left of the lower right
stringer and the shear modulus is the
same for all webs.
Problem 7-14
Problem 7-12
7-15
A 100-inch-long, thin-walled hollow beam
is fixed to a wall. The cross section at the
free end is shown in the figure (all
7-13
dimensions are measured at the midline).
Two 600-lb forces produce a torsional
couple (torque). The left side (semi-circle)
of the cross section has thickness t = 0.05"
and G = 12x106 psi, while the right side
(three straight members) has thickness t =
0.1" and G = 5 x 106 psi.
600 lb
Problem 7-13
b) If all webs have the same thickness, t= Calculate the shear flows, and shear stresses
0.08 inches, and G=4x106 psi, compute the in both sides (neglecting possible stress
torsional stiffness. concentrations), and the angle of twist of the
free end.
c) Compute the shear flows if a 1000 lb. (Answers: q= 53.85 lb/in, left = 1,077 psi,
load is applied in the vertical direction 10 right = 538.5 psi, = 1.666e-3 rad =
inches from the left of the lower right
0.0955)
stringer. Show your results on a diagram
Problem 7-18
a) Compute the shear flows in the webs.
All skin stiffeners shown have cross-
sectional area equal to 0.1 in.2 The skin
elements have thickness 0.40 inches
with G=3.9 x 106 psi. Consider using
Excel.
b) Compute the rotation angle if this
7-17 section is an end section of a cantilever
Find the shear flows (magnitude and that is 100 inches long.
direction) in all of the webs for the single
cell cross-section shown if the 6000 lb. 7-19 - A 30“ long tube with a rectangular,
vertical shear resultant is applied at the shear thin-wall cross-section is subjected to a
center. The section is made entirely of torque, T. The tube is constructed by gluing
aluminum with web thicknesses all equal. two U-shaped sections together as indicated.
All six stringer areas are 1 square inch. The If the tube is pulled longitudinally, shear
nose web is closed, but is not a semi-circle. stress in the joint develops.
The developed length of the nose web is 14
inches. Its area is shown on the diagram.
Problem 7-17
7-21
A single cell thin-wall cross-section is
idealized as consisting of 3 webs with equal
thicknesses t=0.064 inches and 3 stringers
equal areas 0.50 in.2 Two shear resultants
are applied as shown. The enclosed cell
area is computed to be 128.54 in.2 The 3
unknown shear flows are labeled in the
figure below.
A mechanism is idealized as a weightless, rigid device shaped like a “T.” It is acted upon by two
equal forces, P. These forces always remain horizontal, no matter how the T-section rotates.
Find the value of P that makes this T-section neutrally stable.