Advanced Mechanics of

Materials and
Applied Elasticity

Sixth Edition
This page intentionally left blank
 Advanced Mechanics of
Materials and
Applied Elasticity

Sixth Edition

Ansel C. Ugural

Saul K. Fenster

Boston • Columbus • New York • San Francisco • Amsterdam • Cape Town

Dubai • London • Madrid • Milan • Munich • Paris • Montreal • Toronto • Delhi • Mexico City
São Paulo • Sydney • Hong Kong • Seoul • Singapore • Taipei • Tokyo
Many of the designations used by manufacturers and sellers to distinguish their products are claimed as trade-
marks. Where those designations appear in this book, and the publisher was aware of a trademark claim, the
designations have been printed with initial capital letters or in all capitals.

The authors and publisher have taken care in the preparation of this book, but make no expressed or implied
warranty of any kind and assume no responsibility for errors or omissions. No liability is assumed for incidental
or consequential damages in connection with or arising out of the use of the information or programs contained
herein.

For information about buying this title in bulk quantities, or for special sales opportunities (which may include
electronic versions; custom cover designs; and content particular to your business, training goals, marketing
focus, or branding interests), please contact our corporate sales department at [email protected] or
(800) 382-3419.

For government sales inquiries, please contact [email protected].

For questions about sales outside the U.S., please contact [email protected].

Visit us on the Web: informit.com

Library of Congress Control Number: 2019932748

Copyright © 2020 Pearson Education, Inc.

Cover illustration: Yulia Grigoryeva/Shutterstock

All rights reserved. This publication is protected by copyright, and permission must be obtained from the
publisher prior to any prohibited reproduction, storage in a retrieval system, or transmission in any form or
by any means, electronic, mechanical, photocopying, recording, or likewise. For information regarding permis-
sions, request forms and the appropriate contacts within the Pearson Education Global Rights & Permissions
Department, please visit www.pearsoned.com/permissions/.

ISBN-13: 978-0-13-485928-6
ISBN-10: 0-13-485928-6

ScoutAutomatedPrintCode
 Contents

Preface xvii
Acknowledgments xx
About the Authors xxi
List of Symbols xxii

Chapter 1  Analysis of Stress 1

1.1 Introduction 1
1.1.1 Mechanics of Materials and Theory of Elasticity 1
1.1.2 Historical Development 2
1.2 Scope of the Book 3
1.3 Analysis and Design 4
1.3.1 Role of Analysis in Design 6
1.3.2 Selection of Factor of Safety 6
1.3.3 Case Studies 7
1.4 Conditions of Equilibrium 8
1.5 Definition and Components of Stress 9
1.5.1 Sign Convention 11
1.5.2 Equality of Shearing Stresses 12
1.5.3 Some Special Cases of Stress 12
1.6 Internal Force Resultant and Stress Relations 13
1.6.1 Basic Formulas for Stress 15
1.6.2 Combined Stresses 17
1.7 Stresses on Inclined Sections 17
1.7.1 Axially Loaded Members 18
1.8 Variation of Stress within a Body 20
1.8.1 Equations of Equilibrium 20
1.9 Plane-Stress Transformation 23
1.9.1 Stress Tensor 25
1.9.2 Polar Representations of State of Plane Stress 25
1.9.3 Cartesian Representation of State of Plane Stress 25

 v
1.10 Principal Stresses and Maximum In-Plane Shear Stress 26
1.11 Mohr’s Circle for Two-Dimensional Stress 28
1.12 Three-Dimensional Stress Transformation 35
1.13 Principal Stresses in Three Dimensions 38
1.13.1 Invariants for Three-Dimensional Stress 40
1.14 Normal and Shear Stresses on an Oblique Plane 42
1.14.1 Octahedral Stresses 44
1.15 Mohr’s Circles in Three Dimensions 45
1.15.1 Absolute Maximum Shear Stress 46
1.15.2 Equations of Three Mohr’s Circles for Stress 47
1.16 Boundary Conditions in Terms of Surface Forces 49
1.17 Indicial Notation 50
References 51
Problems 51

Chapter 2  Strain and Material Properties 68

2.1 Introduction 68
2.2 Deformation 69
2.2.1 Superposition 69
2.3 Strain Defined 70
2.3.1 Plane Strain 70
2.3.2 Three-Dimensional Strain 72
2.3.3 Eulerian and Lagrangian Coordinates 73
2.3.4 Large Strains 74
2.4 Equations of Compatibility 75
2.5 State of Strain at a Point 76
2.5.1 Transformation of Two-Dimensional Strain 76
2.5.2 Transformation of Three-Dimensional Strain 78
2.5.3 Invariants in Three-Dimensional Strain 79
2.5.4 Mohr’s Circle for Plane Strain 80
2.6 Engineering Materials 83
2.6.1 General Properties of Some Common Materials 84
2.7 Stress-Strain Diagrams 86
2.7.1 Ductile Materials in Tension 86
2.7.2 Geometry Change of Specimen 87
2.7.3 True Stress and True Strain 88
2.7.4 Brittle Materials in Tension 89
2.7.5 Materials in Compression 89
2.7.6 Materials in Shear 90
2.7.7 Short-Time Effects of Temperature on Stress-Strain Properties 90
2.8 Elastic versus Plastic Behavior 91
2.9 Hooke’s Law and Poisson’s Ratio 92
2.9.1 Volume Change 93
2.9.2 Deflection of Axially Loaded Members 93

viContents
2.10 Generalized Hooke’s Law 96
2.11 Orthotropic Materials 101
2.11.1 Generalized Hook’s Law for Orthotropic Material 102
2.12 Measurement of Strain: Strain Gage 103
2.12.1 Strain Rosette of Three Gages 104
2.12.2 Rectangular and Delta Strain Rosettes 106
2.13 Strain Energy 107
2.13.1 Strain Energy Density for Normal and Shear Stresses 108
2.13.2 Strain Energy Density for Three-Dimensional Stresses 110
2.14 Strain Energy in Common Structural Members 111
2.14.1 Strain Energy for Axially Loaded Bars 111
2.14.2 Strain Energy of Circular Bars in Torsion 112
2.14.3 Strain Energy for Beams in Bending 113
2.15 Components of Strain Energy 113
2.16 Saint-Venant’s Principle 115
2.16.1 Confirmation of Saint-Venant’s Rule 116
References 117
Problems 118

Chapter 3  Problems in Elasticity 133

3.1 Introduction 133
3.2 Fundamental Principles of Analysis 134
3.2.1 Three-Dimensional Problems 134
3.2.2 Two-Dimensional Problems 134
Part A: Formulation and Methods of Solution 135
3.3 Plane Strain Problems 135
3.4 Plane Stress Problems 138
3.4.1 Stress–Strain Relations for Orthotropic Materials 139
3.5 Comparison of Two-Dimensional Isotropic Problems 140
3.6 Airy’s Stress Function 141
3.6.1 Generalized Plane Strain Problems 142
3.6.2 Antiplane Shear Deformations 142
3.7 Solution of Elasticity Problems 143
3.7.1 Polynomial Solutions 144
3.8 Thermal Stresses 149
3.8.1 Equations of Thermoelasticity 149
3.9 Basic Relations in Polar Coordinates 152
3.9.1 Equations of Equilibrium 153
3.9.2 Stress Function 153
3.9.3 Strain-Displacement Relations 154
3.9.4 Hooke’s Law 155
3.9.5 Transformation Equations 155
3.9.6 Compatibility Equation 156

Contents vii
Part B: Stress Concentrations 157
3.10 Stresses Due to Concentrated Loads 157
3.10.1 Compression of a Wedge (Fig. 3.10a) 157
3.10.2 Bending of a Wedge (Fig. 3.10b) 159
3.10.3 Concentrated Load on a Straight Boundary (Fig. 3.11a) 160
3.11 Stress Distribution Near a Concentrated Load Acting on a Beam 161
3.11.1 Accuracy of Results 163
3.12 Stress Concentration Factors 163
3.12.1 Circular Hole in a Large Plate in Simple Tension 165
3.12.2 Circular Hole in a Large Plate in Biaxial Tension 167
3.12.3 Elliptic Hole in a Large Plate in Tension 167
3.12.4 Graphs for Stress Concentration Factors 168
Part C: Contact Mechanics 169
3.13 Contact Stresses and Deflections 169
3.13.1 Hertz Theory 170
3.13.2 Johnson–Kendall–Roberts Theory 170
3.14 Spherical and Cylindrical Contacts 171
3.14.1 Two Spheres in Contact 171
3.14.2 Two Parallel Cylinders in Contact 173
3.15 Contact Stress Distribution 174
3.15.1 Two Spheres in Contact (Figure 3.18a) 174
3.15.2 Two Parallel Cylinders in Contact (Figure 3.20a) 174
3.16 General Contact 178
References 181
Problems 182

Chapter 4  Failure Criteria 192

4.1 Introduction192
4.1.1 Failure 192
Part A: Static Loading 193
4.2 Failure by Yielding 193
4.2.1 Creep: Time-Dependent Deformation 194
4.3 Failure by Fracture 195
4.3.1 Types of Fracture in Tension 196
4.4 Yield and Fracture Criteria 197
4.5 Maximum Shearing Stress Theory 198
4.6 Maximum Distortion Energy Theory 199
4.6.1 Yield Surfaces for Triaxial Stress 200
4.7 Octahedral Shearing Stress Theory 200
4.8 Comparison of the Yielding Theories 204
4.9 Maximum Principal Stress Theory 205
4.10 Mohr’s Theory 206
4.11 Coulomb–Mohr Theory 207
4.12 Introduction to Fracture Mechanics 210
4.12.1 Stress-Intensity Factors 211

viiiContents
4.13 Fracture Toughness 213
Part B: Repeated and Dynamic Loadings 216
4.14 Fatigue: Progressive Fracture 216
4.14.1 Fatigue Tests 216
4.14.2 Estimating the Endurance Limit and Fatigue Strength 217
4.15 Failure Criteria for Metal Fatigue 217
4.15.1 Uniaxial State of Stress 218
4.15.2 Comparison of Fatigue Failure Criteria 219
4.15.3 Design for Uniaxial Stress 219
4.15.4 Combined State of Stress 221
4.16 Fatigue Life 223
4.17 Impact Loads 225
4.17.1 Strain Rate 226
4.17.2 Basic Assumptions of Impact Analysis 227
4.18 Longitudinal and Bending Impact 227
4.18.1 Freely Falling Weight 227
4.18.2 Horizontally Moving Weight 228
4.19 Ductile–Brittle Transition 230
References 232
Problems 233

Chapter 5  Bending of Beams 242

5.1 Introduction 242
Part A: Exact Solutions 243
5.2 Pure Bending of Beams of Symmetrical Cross Section 243
5.2.1 Kinematic Relationships 244
5.2.2 Timoshenko Beam Theory 246
5.3 Pure Bending of Beams of Asymmetrical Cross Section 246
5.3.1 Stress Distribution 248
5.3.2 Transformation of Inertia Moments 248
5.4 Bending of a Cantilever of Narrow Section 251
5.4.1 Comparison of the Results with the Elementary Theory Results 253
5.5 Bending of a Simply Supported Narrow Beam 254
5.5.1 Use of Stress Functions 255
5.5.2 Comparison of the Results with the Elementary Theory Results 256
Part B: Approximate Solutions 256
5.6 Elementary Theory of Bending 256
5.6.1 Assumptions of Elementary Theory 257
5.6.2 Method of Integration 258
5.7 Normal and Shear Stresses 260
5.7.1 Rectangular Cross Section 262
5.7.2 Various Cross Sections 262
5.7.3 Beam of Constant Strength 267
5.8 Effect of Transverse Normal Stress 268

Contents ix
5.9 Composite Beams 270
5.9.1 Transformed Section Method 270
5.9.2 Equation of Neutral Axis 271
5.9.3 Stresses in the Transformed Beam 272
5.9.4 Composite Beams of Multi Materials 272
5.10 Shear Center 276
5.10.1 Thin-Walled Open Cross Sections 277
5.10.2 Arbitrary Solid Cross Sections 281
5.11 Statically Indeterminate Systems 281
5.11.1 The Method of Superposition 282
5.12 Energy Method for Deflections 284
5.12.1 Form Factor for Shear 285
Part C: Curved Beams 286
5.13 Elasticity Theory 286
5.13.1 Equations of Equilibrium and Compatibility 286
5.13.2 Boundary Conditions 287
5.13.3 Stress Distribution 288
5.13.4 Deflections 289
5.14 Curved Beam Formula 289
5.14.1 Basic Assumptions 289
5.14.2 Location of the Neutral Axis 290
5.14.3 Tangential Stress 291
5.14.4 Winkler’s Formula 293
5.15 Comparison of the Results of Various Theories 293
5.15.1 Correction of σ θ for Beams with Thin-Walled Cross Sections 294
5.16 Combined Tangential and Normal Stresses 296
References 300
Problems 300

Chapter 6  Torsion of Prismatic Bars 315

6.1 Introduction 315
6.2 Elementary Theory of Torsion of Circular Bars 316
6.2.1 Shearing Stress 317
6.2.2 Angle of Twist 317
6.2.3 Axial and Transverse Shear Stresses 320
6.3 Stresses on Inclined Planes 321
6.3.1 Stress Transformation 321
6.3.2 Transmission of Power by Shafts 323
6.4 General Solution of the Torsion Problem 324
6.4.1 Geometry of Deformation 324
6.4.2 Equations of Equilibrium 325
6.4.3 Equations of Compatibility 325
6.5 Prandtl’s Stress Function 326
6.5.1 Boundary Conditions 326

xContents
 6.5.2 Force and Moments over the Ends 327
6.5.3 Circular Cross Section 331
6.6 Prandtl’s Membrane Analogy 333
6.6.1 Equation of Equilibrium 333
6.6.2 Shearing Stress and Angle of Twist 335
6.7 Torsion of Narrow Rectangular Cross Section 338
6.7.1 Thin-Walled Open Cross Sections 339
6.8 Torsion of Multiply Connected Thin-Walled Sections 340
6.8.1 Shearing Stress 340
6.8.2 Angle of Twist 341
6.9 Fluid Flow Analogy and Stress Concentration 344
6.10 Torsion of Restrained Thin-Walled Members of Open Cross Section 346
6.10.1 Torsional and Lateral Shears 347
6.10.2 Boundary Conditions 348
6.10.3 Long Beams Under Torsion 348
6.10.4 Angle of Twist 348
6.11 Torsion Bar Springs 350
6.12 Curved Circular Bars 351
6.12.1 Helical Springs 352
References 354
Problems 355

Chapter 7  Numerical Methods 364

7.1 Introduction 364
Part A: Finite Difference Analysis 365
7.2 Finite Differences 365
7.2.1 Central Differences 366
7.3 Finite Difference Equations 368
7.4 Curved Boundaries 370
7.5 Boundary Conditions 373
Part B: Finite Element Analysis 377
7.6 Fundamentals 377
7.7 The Bar Element 379
7.7.1 Equilibrium Method 379
7.7.2 Energy Method 379
7.8 Arbitrarily Oriented Bar Element 380
7.8.1 Coordinate Transformation 380
7.8.2 Force Transformation 381
7.8.3 Displacement Transformation 383
7.8.4 Governing Equations 383
7.9 Axial Force Equation 384
7.10 Force-Displacement Relations for a Truss 386
7.10.1 The Assembly Process 386
7.11 Beam Element 393

Contents xi
7.12 Properties of Two-Dimensional Elements 399
7.12.1 Displacement Matrix 399
7.12.2 Strain, Stress, and Elasticity Matrices 401
7.13 General Formulation of the Finite Element Method 402
7.13.1 Outline of General Finite Element Analysis 403
7.14 Triangular Finite Element 407
7.14.1 Element Nodal Forces 410
7.15 Case Studies in Plane Stress 414
7.16 Computational Tools 423
References 423
Problems 424

Chapter 8  Thick-Walled Cylinders and Rotating Disks 434

8.1 Introduction 434
8.1.1 Basic Relations 434
8.2 Thick-Walled Cylinders Under Pressure 435
8.2.1 Special Cases 438
8.2.2 Closed-Ended Cylinder 440
8.3 Maximum Tangential Stress 441
8.4 Application of Failure Theories 442
8.5 Compound Cylinders: Press or Shrink Fits 443
8.6 Rotating Disks of Constant Thickness 446
8.6.1 Annular Disk 447
8.6.2 Solid Disk 448
8.7 Disk Flywheels 449
8.7.1 Design Factors 450
8.7.2 Stresses and Displacement 450
8.8 Rotating Disks of Variable Thickness 453
8.9 Rotating Disks of Uniform Stress 456
8.10 Thermal Stresses in Thin Disks 458
8.10.1 Annular Disk 459
8.10.2 Solid Disk 459
8.11 Thermal Stress in Long Circular Cylinders 460
8.11.1 Solid Cylinder 460
8.11.2 Cylinder with a Central Circular Hole 461
8.11.3 Special Case 464
8.12 Finite Element Solution 464
8.12.1 Axisymmetric Element 464
References 466
Problems 466

Chapter 9  Beams on Elastic Foundations 473

9.1 Introduction 473
9.2 General Theory 473

xiiContents
9.3 Infinite Beams 475
9.4 Semi-Infinite Beams 480
9.5 Finite Beams 483
9.6 Classification of Beams 484
9.7 Beams Supported by Equally Spaced Elastic Elements 485
9.8 Simplified Solutions for Relatively Stiff Beams 486
9.9 Solution by Finite Differences 488
9.10 Applications 490
9.10.1 Grid Configurations of Beams 490
References 492
Problems 493

Chapter 10  Applications of Energy Methods 496

10.1 Introduction 496
Part A: Energy Principles  497
10.2 Work Done in Deformation 497
10.3 Reciprocity Theorem 498
10.4 Castigliano’s Theorem 499
10.4.1 Application to Bars and Beams 500
10.4.2 Application to Trusses 500
10.4.3 Use of a Fictitious Load 501
10.5 Unit- or Dummy-Load Method 506
10.6 Crotti–Engesser Theorem 508
10.7 Statically Indeterminate Systems 510
Part B: Variational Methods 514
10.8 Principle of Virtual Work 514
10.8.1 Variation in Strain Energy 514
10.8.2 Virtual Work Done by Forces 515
10.9 Principle of Minimum Potential Energy 515
10.10 Deflections by Trigonometric Series 517
10.10.1 Strain Energy 518
10.10.2 Virtual Work 518
10.11 Rayleigh–Ritz Method 522
References 524
Problems 525

Chapter 11  Stability of Columns 534

11.1 Introduction 534
11.2 Critical Load 534
11.2.1 Equilibrium Method 535
11.2.2 Energy Method 536
11.3 Buckling of Pin-Ended Columns 536
11.3.1 Modes of Buckling 538
11.4 Deflection Response of Columns 539

Contents xiii
 11.4.1 Effects of Large Deflections 539
11.4.2 Effects of Imperfections 540
11.4.3 Effects of Inelastic Behavior 540
11.5 Columns with Different End Conditions 540
11.6 Critical Stress: Classification of Columns 543
11.6.1 Long Columns 543
11.6.2 Short Columns 544
11.6.3 Intermediate Columns: Inelastic Buckling 544
11.7 Design Formulas for Columns 548
11.8 Imperfections in Columns 550
11.9 Local Buckling of Columns 552
11.10 Eccentrically Loaded Columns: Secant Formula 552
11.10.1 Simplified Formula for Short Columns 554
11.11 Energy Methods Applied to Buckling 554
11.12 Solution by Finite Differences 562
11.13 Finite Difference Solution for Unevenly Spaced Nodes 567
References 568
Problems 569

Chapter 12  Plastic Behavior of Materials 578

12.1 Introduction 578
12.2 Plastic Deformation 579
12.2.1 Slip Action: Dislocation 579
12.3 Idealized Stress–Strain Diagrams 580
12.3.1 True Stress–True Strain Relationships 580
12.4 Instability in Simple Tension 582
12.5 Plastic Axial Deformation and Residual Stress 585
12.6 Plastic Deflection of Beams 588
12.7 Analysis of Perfectly Plastic Beams 590
12.7.1 Shape Factor 593
12.7.2 Plastic Hinge 593
12.8 Collapse Load of Structures: Limit Design 600
12.8.1 Collapse Mechanism 600
12.8.2 Ultimate Load by the Energy Method 601
12.9 Elastic–Plastic Torsion of Circular Shafts 605
12.9.1 Yield Torque 606
12.9.2 Elastic–Plastic Torque 606
12.9.3 Ultimate Torque 607
12.9.4 Residual Rotation and Stress 608
12.10 Plastic Torsion: Membrane Analogy 610
12.10.1 Membrane–Roof Analogy 610
12.10.2 Sand Hill Analogy 611
12.11 Elastic–Plastic Stresses in Rotating Disks 612
12.11.1 Initial Yielding 612

xivContents
 12.11.2 Partial Yielding 612
12.11.3 Complete Yielding 614
12.12 Plastic Stress–Strain Relations 614
12.13 Plastic Stress–Strain Increment Relations 620
12.14 Stresses in Perfectly Plastic Thick-Walled Cylinders 623
12.14.1 Complete Yielding 624
12.14.2 Partial Yielding 626
References 627
Problems 628

Chapter 13  Stresses in Plates and Shells 635

13.1 Introduction 635
Part A: Bending of Thin Plates 635
13.2 Basic Assumptions 635
13.3 Strain–Curvature Relations 636
13.4 Stress, Curvature, and Moment Relations 638
13.5 Governing Equations of Plate Deflection 640
13.6 Boundary Conditions 642
13.7 Simply Supported Rectangular Plates 644
13.8 Axisymmetrically Loaded Circular Plates 648
13.9 Deflections of Rectangular Plates by the Strain-Energy Method 650
13.10 Sandwich Plates 652
13.10.1 Design of Sandwich Beams and Plates 653
13.11 Finite Element Solution 654
13.11.1 Strain, Stress, and Elasticity Matrices 655
13.11.2 Displacement Function 655
13.11.3 Stiffness Matrix 657
13.11.4 External Nodal Forces 657
Part B: Membrane Stresses in Thin Shells 657
13.12 Theories and Behavior of Shells 657
13.13 Simple Membrane Action 658
13.14 Symmetrically Loaded Shells of Revolution 660
13.14.1 Equations of Equilibrium 661
13.14.2 Conditions of Compatibility 662
13.15 Some Typical Cases of Shells of Revolution 662
13.15.1 Spherical Shell 662
13.15.2 Conical Shell 663
13.15.3 Circular Cylindrical Shell 664
13.16 Thermal Stresses in Compound Cylinders 668
13.17 Cylindrical Shells of General Shape 670
References 673
Problems 673

Contents xv
Appendix A  Problem Formulation and Solution 679
A.1 Basic Method 679
A.1.1 Numerical Accuracy 680
A.1.2 Daily Planning 680

Appendix B  Solution of the Stress Cubic Equation 682

B.1 Principal Stresses 682
B.1.1 Direction Cosines 683

Appendix C  Moments of Composite Areas 687

C.1 Centroid 687
C.2 Moments of Inertia 690
C.2.1 Parallel Axis Theorem 690
C.2.2 Principal Moments of Inertia 692

Appendix D  Tables and Charts 699

D.1 Charts of Stress Concentration Factors 705

Appendix E  Introduction to MATLAB 710

Answers to Selected Problems 713

Index 722

xviContents
 Preface

INTRODUCTION

Advanced Mechanics of Materials and Applied Elasticity, Sixth Edition, is an outgrowth

of classroom notes prepared in connection with advanced undergraduate and first-year
graduate courses in the mechanics of solids and elasticity. It is designed to satisfy the
requirements of courses subsequent to an elementary treatment of the strength of materi-
als. In addition to its applicability to aeronautical, civil, and mechanical engineering and
to engineering mechanics curricula, the text is useful to practicing engineers. Emphasis is
given to numerical techniques (which lend themselves to computerization) in the solution
of problems resisting analytical treatment. The attention devoted to numerical solutions is
not intended to deny the value of classical analysis, which is given a rather full treatment.
Instead, the coverage provided here seeks to fill what we believe to be a void in the world
of textbooks.
We have attempted to present a balance between the theory necessary to gain insight
into the mechanics, but which can often offer no more than crude approximations to real
problems because of simplifications related to geometry and conditions of loading, and
numerical solutions, which are so useful in presenting stress analysis in a more realistic
setting. This text emphasizes those aspects of theory and application that prepare a student
for more advanced study or for professional practice in design and analysis.
The theory of elasticity plays three important roles in the text. First, it provides exact
solutions where the configurations of loading and boundary are relatively simple. Second,
it provides a check on the limitations of the mechanics of materials approach. Third, it
serves as the basis of approximate solutions employing numerical analysis.
To make the text as clear as possible, the fundamentals of the mechanics of materi-
als are addressed as necessary. The physical significance of the solutions and practical
applications are also emphasized. In addition, we have made a special effort to illus-
trate important principles and applications with numerical examples. Consistent with
announced national policy, problems are included in the text in which the physical quan-
tities are expressed in the International System of Units (SI). All important quantities
are defined in both SI and U.S. Customary System (USCS) of units. A sign convention,
consistent with vector mechanics, is employed throughout for loads, internal forces, and
stresses. This convention conforms to that used in most classical strength of materials and

 xvii
elasticity texts, as well as to that most often employed in the numerical analysis of com-
plex structures.

ORGANIZATION OF THE TEXT

Because of its extensive subdivision into a variety of topics and use of alternative meth-
ods of analysis, this text provides great flexibility for instructors when choosing assign-
ments to cover courses of varying length and content. Most chapters are substantially
self-­contained, so the order of presentation can be smoothly altered to meet an instructor’s
preference. Ideally, Chapters 1 and 2, which address the analysis of basic concepts, should
be studied first. The emphasis placed on the treatment of two-dimensional problems in
elasticity (Chapter 3) may then differ according to the scope of the course.
This sixth edition of Advanced Mechanics of Materials and Applied Elasticity seeks to
preserve the objectives and emphases of the previous editions. Every effort has been made
to provide a more complete and current text through the inclusion of new material dealing
with the fundamental principles of stress analysis and design: stress concentrations, contact
stresses, failure criteria, fracture mechanics, compound cylinders, finite element analysis
(FEA), energy and variational methods, buckling of stepped columns, common shell types,
case studies in analysis and design, and MATLAB solutions. The entire text has been reex-
amined, and many improvements have been made throughout by a process of elimination
and rearrangement. Some sections have been expanded to improve on previous expositions.
The references (identified in brackets), which are provided as an aid to those students
who wish to pursue certain aspects of a subject in further depth, have been updated and
listed at the end of each chapter. We have resisted the temptation to increase the mate-
rial covered except where absolutely necessary. Nevertheless, we have added a number
of illustrative examples and problems important in engineering practice and design. Extra
care has been taken in the presentation and solution of the sample problems. All the prob-
lem sets have been reviewed and checked to ensure both their clarity and their numerical
accuracy. Most changes in subject-matter coverage were prompted by the suggestions of
faculty familiar with earlier editions.
In this sixth edition, we have maintained the previous editions’ clarity of presentation,
simplicity as the subject permits, unpretentious depth, an effort to encourage intuitive
understanding, and a shunning of the irrelevant. In this context, as throughout, emphasis
is placed on the use of fundamentals to help build students’ understanding and ability to
solve the more complex problems.

SUPPLEMENTS

The book is accompanied by a comprehensive instructor’s Solutions Manual. Written and

class tested, it features complete solutions to all problems in the text. Answers to selected
problems are given at the end of the book. The password-protected Solutions Manual is
available for adopters at the Pearson Instructor Resource Center, pearsonhighered.com/irc.

xviiiPreface
 Optional Material is also available from the Pearson Resource Center, pearsonhigh-
ered.com/irc. This material includes PowerPoint slides of figures and tables, and solutions
using MATLAB for a variety of sample problems of practical importance. The book, how-
ever, is independent of any software package.

Register your copy of Advanced Mechanics of Materials and Applied Elasticity,

Sixth Edition, on the InformIT site for convenient access to updates and corrections
as they become available. To start the registration process, go to informit.com/reg-
ister and log in or create an account. Enter the product ISBN (9780134859286) and
click Submit. Look on the Registered Products tab for an Access Bonus Content link
next to this product, and follow that link to access any available bonus materials. If
you would like to be notified of exclusive offers on new editions and updates, please
check the box to receive email from us.

Preface xix
 Acknowledgments

It is a particular pleasure to acknowledge the contributions of those who assisted in the

evolution of the text. Thanks, of course, are due to the many readers who have contributed
general ideas and to reviewers who have made detailed comments on previous editions.
These notably include the following: F. Freudenstein, Columbia University; R. A. Scott,
University of Michigan; M. W. Wilcox and Y. Chan Jian, Southern Methodist University;
C. T. Sun, University of Florida; B. Koplik, H. Kountouras, K. A. Narh, R. Sodhi, and
C. E. Wilson, New Jersey Institute of Technology; H. Smith, Jr., South Dakota School
of Mines and Technology; B. P. Gupta, Gannon University; S. Bang, University of Notre
Dame; B. Koo, University of Toledo; J. T. Easley, University of Kansas; J. A. Bailey,
North Carolina State University; W. F. Wright, Vanderbilt University; R. Burks, SUNY
Maritime College; G. E. O. Widera, University of Illinois; R. H. Koebke, University of
South Carolina; B. M. Kwak, University of Iowa; G. Nadig, Widener University; R. L.
Brown, Montana State University; S. H. Advani, West Virginia University; E. Nassat,
Illinois Institute of Technology; R. I. Sann, Stevens Institute of Technology; C. O. Smith,
University of Nebraska; J. Kempner, Polytechnic University of New York; and P. C. Prister,
North Dakota State University; R. Wetherhold, University of Buffalo, SUNY; and Shaofan
Li, University of California at Berkeley.
Accuracy checking of the problems, proofreading, and typing of the Solutions
Manual were done expertly by my former student, Dr. Youngjin Chung. Also contribut-
ing considerably to this volume with typing new inserts, scanning some figures, limited
proofreading, and cover design was Errol A. Ugural. Their hard work is much appreciated.
I am deeply indebted to my colleagues who have found the text useful through the years,
as well as to Laura Lewin, executive editor at Pearson, who encouraged development of
this edition. Lastly, I am very thankful for the support and understanding of my wife Nora,
daughter Aileen, and son Errol during preparation of this book.

Ansel C. Ugural
Holmdel, New Jersey

xx
 About the Authors

Ansel C. Ugural, Ph.D., has been a research and visiting professor of mechanical and
civil engineering at the New Jersey Institute of Technology. He has taught in the engineer-
ing mechanics department at the University of Wisconsin. Dr. Ugural also served as chair-
man and tenured professor of mechanical engineering at Fairleigh Dickinson University
for twenty years. He has considerable and diverse industrial experience in both full-time
and consulting capacities as a design, development, and research engineer.
Professor Ugural earned his MS in mechanical engineering and PhD in engineer-
ing mechanics from the University of Wisconsin–Madison. Dr. Ugural was a National
Science Foundation (NSF) fellow. He has been a member of several professional societies,
including the American Society of Mechanical Engineers and the American Society of
Engineering Education. He is also listed in Who’s Who in Engineering.
Dr. Ugural is the author of several books, including Mechanics of Materials; Mechanical
Design: An Integrated Approach; Mechanical Design of Machine Components; Stresses
in Beams, Plates, and Shells; and Plates and Shells: Theory and Analysis. Most of these
texts have been translated into Korean, Chinese, and Portuguese. In addition, Professor
Ugural has published numerous articles in trade and professional journals.

Saul K. Fenster, Ph.D., served as president and tenured professor at New Jersey
Institute of Technology for more than two decades. In addition, he has held varied posi-
tions at Fairleigh Dickinson University and taught at the City University of New York.
His experience includes membership on a number of corporate boards and economic
development commissions. Fenster is a fellow of the American Society of Mechanical
Engineers, the American Society for Engineering Education, and the American Society
for Manufacturing Engineers. He is coauthor of a text on mechanics.

 xxi
 Symbols

Roman Letters
A area
B width
C carryover factor, torsional rigidity
c distance from neutral axis to outer fiber
D distribution factor, flexural rigidity of plate
[D] elasticity matrix
d diameter, distance
E modulus of elasticity in tension or compression
Es modulus of plasticity or secant modulus
Et tangent modulus
e dilatation, distance, eccentricity
{F} nodal force matrix of bar and beam finite elements
F body force per unit volume, concentrated force
f coefficient of friction
{   f} displacement function of finite element
G modulus of elasticity in shear or modulus of rigidity
g acceleration of gravity (≈9.81 m/s2)
h depth of beam, height, membrane deflection, mesh width
I moment of inertia of area, stress invariant
J polar moment of inertia of area, strain invariant
K bulk modulus, spring constant of an elastic support, stiffness ­factor, thermal con-
ductivity, fatigue factor, strength coefficient, stress concentration factor
[K  ] stiffness matrix of whole structure
k constant, modulus of elastic foundation, spring constant
[k] stiffness matrix of finite element
L length, span
l, m, n direction cosines
M moment
Mxy twisting moment in plates
m moment caused by unit load

xxii
N fatigue life (cycles), force
n factor of safety, number, strain hardening index
P concentrated force
p distributed load per unit length or area, pressure, stress resultant
Q first moment of area, heat flow per unit length, shearing force
{Q} nodal force matrix of two-dimensional finite element
R radius, reaction
r radius, radius of gyration
r, θ polar coordinates
S elastic section modulus, shear center
s distance along a line or a curve
T temperature, twisting couple or torque
t thickness
U strain energy
Uo strain energy per unit volume
U* complementary energy
u, u, w components of displacement
V shearing force, volume
u velocity
W weight, work
x, y, z rectangular coordinates
Z plastic section modulus

Greek Letters
α angle, coefficient of thermal expansion, form factor for shear
b numerical factor, angle
γ shear strain, weight per unit volume or specific weight, angle
δ deflection, finite difference operator, variational symbol, displacement
{δ} nodal displacement matrix of finite element
∆
change of a function
ε normal strain
θ angle, angle of twist per unit length, slope
ν Poisson’s ratio
λ axial load factor, Lamé constant
Π potential energy
ρ density (mass per unit volume), radius
σ normal stress
τ shear stress
φ total angle of twist
F stress function
ω angular velocity
ψ stream function

Symbols xxiii
This page intentionally left blank
 CHAPTER 5

Bending of Beams

5.1 INTRODUCTION

In this chapter we are concerned with the bending of straight as well as curved beams—
that is, structural elements possessing one dimension significantly greater than the other
two, usually loaded in a direction normal to the longitudinal axis. We first examine the
elasticity or “exact” solutions of beams that are straight and made of homogeneous, lin-
early elastic materials. Then, we consider solutions for straight beams using mechanics
of materials or elementary theory, special cases involving members made of composite
materials, and the shear center. The deflections and stresses in beams caused by pure
bending as well as those due to lateral loading are discussed. We analyze stresses in
curved beams using both exact and elementary methods, and compare the results of the
various theories.
Except in the case of very simple shapes and loading systems, the theory of elastic-
ity yields beam solutions only with considerable difficulty. Practical considerations often
lead to assumptions about stress and deformation that result in mechanics of materials or
elementary theory solutions. The theory of elasticity can sometimes be applied to test the
validity of such assumptions. This theory has three roles in these problems: It can serve to
place limitations on the use of the elementary theory, it can be used as the basis of approx-
imate solutions through numerical analysis, and it can provide exact solutions for simple
configurations of loading and shape.

242
Part A: Exact Solutions

5.2  PURE BENDING OF BEAMS OF SYMMETRICAL CROSS SECTION

The simplest case of pure bending is that of a beam possessing a vertical axis of sym-
metry, subjected to equal and opposite end couples (Fig. 5.1a). The semi-inverse method
is now applied to analyze this problem. The moment M z shown in Fig. 5.1a is defined as
positive, because it acts on a positive (negative) face with its vector in the positive (nega-
tive) coordinate direction. This sign convention agrees with that of stress (Section 1.5).
We will assume that the normal stress over the cross section varies linearly with y and that
the remaining stress components are zero:

σ x = ky, σ y = σ z = τ xy = τ xz = τ yz = 0 (5.1)

Here k is a constant, and y = 0 contains the neutral surface—that is, the surface along
which σ x = 0. The intersection of the neutral surface and the cross section locates the neu-
tral axis (abbreviated NA). Figure 5.1b shows the linear stress field in a section located an
arbitrary distance a from the left end.
Since Eqs. (5.1) indicate that the lateral surfaces are free of stress, we need only be
assured that the stresses are consistent with the boundary conditions at the ends. These
conditions of equilibrium require that the resultant of the internal forces be zero and that
the moments of the internal forces about the neutral axis equal the applied moment

∫Aσ x dA = 0, − ∫ yσ x dA = M z (5.2)
A

where A is the cross-sectional area. Note that the zero stress components τ xy ,τ xz in
Eqs. (5.1) satisfy the conditions that no y- and z-directed forces exist at the end faces.
Moreover, because of the y symmetry of the section, σ x = ky produces no moment about
the y axis. The negative sign in the second expression implies that a positive moment M z
is one that results in compressive (negative) stress at points of positive y. Substituting Eqs.
(5.1) into Eqs. (5.2) yields

k ∫ y dA = 0, − k ∫ y 2 dA = M z (5.3a, b)
A A

y
Mz
x x z
Mz N.A.
Mz σx
a Mz z dA
y a
y y
(a) (b)

Figure 5.1.  ( a) Beam of singly symmetric cross section in pure bending; (b) stress distribution
across cross section of the beam.

5.2   Pure Bending of Beams of Symmetrical Cross Section 243

 Since k ≠ 0, Eq. (5.3a) indicates that the first moment of cross-sectional area about the
neutral axis is zero. This requires that the neutral and centroidal axes of the cross section
coincide. Neglecting body forces, it is clear that the equations of equilibrium (3.4), are
satisfied by Eqs. (5.1). It can also readily be verified that Eqs. (5.1) together with Hooke’s
law fulfill the compatibility conditions, Eq. (2.12). Thus, Eqs. (5.1) represent an exact
solution.
The integral in Eq. (5.3b) defines the moment of inertia I z of the cross section about
the z axis of the beam cross section (Appendix C); therefore,
Mz
k=− (a)
Iz

An expression for normal stress can now be written by combining Eqs. (5.1) and (a):
Mz y
σx = − (5.4)
Iz

This is the familiar elastic flexure formula applicable to straight beams.

Since, at a given section, M and I are constant, the maximum stress is obtained from
Eq. (5.4) by taking | y|max = c:

Mc M M
σ max = = =
I I /c S (5.5)

where S is the elastic section modulus. Equation (5.5) is widely employed in practice
because of its simplicity. To facilitate its use, section moduli for numerous common sec-
tions are tabulated in various handbooks. A fictitious stress in extreme fibers, computed
from Eq. (5.5) for the experimentally obtained ultimate bending moment (Section 12.7), is
termed the modulus of rupture of the material in bending. This quantity, σ max = M u /S , is
frequently used as a measure of the bending strength of materials.

5.2.1  Kinematic Relationships

To gain further insight into the beam problem, we now consider the geometry of
­deformation—that is, beam kinematics. Fundamental to this discussion is the hypothesis
that sections originally plane remain so subsequent to bending. For a beam of symmetrical
cross section, Hooke’s law and Eq. (5.4) lead to

Mz y M y
εx = − , ε y = εz = ν z
EI z EI z (5.6)
γ xy = γ xz = γ yz = 0

where EI z is the flexural rigidity.

Let us examine the deflection of the beam axis, whose axial deformation is zero.
Figure 5.2a shows an element of an initially straight beam, now in a deformed state.
Because the beam is subjected to pure bending, uniform throughout, each element of

244 Chapter  5   Bending of Beams

 O2

rz

x
v εx ds
y m −y
+Mz m +Mz Beam
n ds n
+Mz axis
+Mz
rx
dθ rx Beam
axis
dθ
O1 O1

(a) (b)

Figure 5.2.  (a) Segment of a bent beam; (b) geometry of deformation.

infinitesimal length experiences identical deformation, with the result that the beam cur-
vature is everywhere the same. The deflected axis of the beam or the deflection curve is
thus shown deformed, with radius of curvature rx . The curvature of the beam axis in the xy
plane in terms of the y deflection υ is

1 d 2υ/dx 2 d 2υ
= 2 3/2 ≈ (5.7)
rx [1 + ( dυ/dx ) ] dx 2

where the approximate form is valid for small deformations (dυ /dx  1). The sign con-
vention for curvature of the beam axis is such that this sign is positive when the beam is
bent concave downward, as shown in the figure.
As shown by the geometry in Fig. 5.2b, the shaded sectors are similar. Hence, the
radius of curvature and the strain are related as follows:
ds ε ds
dθ = = − x (5.8)
rx y
where ds is the arc length mn along the longitudinal axis of the beam. For a small dis-
placement, ds ≈ dx and θ represents the slope dυ/dx of the beam axis. Clearly, for the
positive curvature shown in Fig. 5.2a, θ increases as we move from left to right along the
beam axis. On the basis of Eqs. (5.6) and (5.8),
1 ε M
= − x = z (5.9a)
rx y EI z

Following a similar procedure and noting that ε z ≈ −νε x , we may also obtain the curvature
in the yz plane as
1 ε vM z
=– z =– (5.9b)
rz y EI z

5.2   Pure Bending of Beams of Symmetrical Cross Section 245

 The basic equation of the deflection curve of a beam is obtained by combining Eqs.
(5.7) and (5.9a) as follows:

d 2υ M z
= (5.10)
dx 2 EI z

This expression, relating the beam curvature to the bending moment, is known as the
Bernoulli–Euler law of elementary bending theory. It is observed from Fig. 5.2 and
Eq. (5.10) that a positive moment produces a positive curvature. If the sign convention
adopted in this section for either moment or deflection (and curvature) is reversed, the
plus sign in Eq. (5.10) should likewise be reversed.
Reference to Fig. 5.2a reveals that the top and bottom lateral surfaces have been
deformed into saddle-shaped or anticlastic surfaces of curvature 1/rz . The vertical sides
have been simultaneously rotated as a result of bending. Examining Eq. (5.9b) suggests a
method for determining Poisson’s ratio [Ref. 5.1]. For a given beam and bending moment,
a measurement of 1/rz leads directly to v. The effect of anticlastic curvature is small when
the beam depth is comparable to its width.

5.2.2  Timoshenko Beam Theory

The Timoshenko theory of beams, developed by S. P. Timoshenko at the beginning
of the twentieth century, constitutes an improvement over the Euler–Bernoulli the-
ory. In the static case, the difference between the two hypotheses is that the former
includes the effect of shear stresses on the deformation by assuming a constant shear
over the beam height, whereas the latter ignores the influence of transverse shear on
beam deformation. The Timoshenko theory is also said to be an extension of the ordi-
nary beam theory that allows for the effect of the transverse shear deformation while
­relaxing the assumption that plane sections remain plane and normal to the deformed
beam axis.
The Timoshenko beam theory is well suited to describing the behavior of short
beams and sandwich composite beams. In the dynamic case, the theory incorporates
shear deformation as well as rotational inertia effects, and it will be more accurate for
not very slender beams. By effectively taking into account the mechanism of deforma-
tion, Timoshenko’s theory lowers the stiffness of the beam, with the result being a larger
deflection under static load and lower predicted fundamental frequencies of vibration for a
prescribed set of boundary conditions.

5.3  PURE BENDING OF BEAMS OF ASYMMETRICAL CROSS SECTION

In this section, we extend the discussion in Section 5.2 to the more general case in which
a beam of arbitrary cross section is subjected to end couples M y and M z about the y and
z axes, respectively (Fig. 5.3). Following a procedure similar to that described in Section
5.2, plane sections are again taken to remain plane. Assume that the normal stress σ x

246 Chapter  5   Bending of Beams

 Figure 5.3.  P
 ure bending of beams
z′ of asymmetrical cross
section.
C Mz
N.A. z
φ
y
z
My dA

θ
y y′

acting at a point within dA is a linear function of the y and z coordinates of the point;
assume further that the remaining stresses are zero. Then the stress field is

σ x = c1 + c2 y + c3 z (5.11)
σ y = σ z = τ xy = τ xz = τ yz = 0
where c1 , c2 , c3 are constants to be evaluated.
The equilibrium conditions at the beam ends, as before, relate to the force and bend-
ing moment:

∫Aσ x dA = 0 (a)
∫A zσ x dA = M y , − ∫ yσ x dA = M z
A
(b, c)

Carrying σ x , as given by Eq. (5.11), into Eqs. (a), (b), and (c) results in the following
expressions:
(d)
c1 ∫ dA + c2 ∫ y dA + c3 ∫ z dA = 0
A A A

(e)
c1 ∫ z dA + c2 ∫ yz dA + c3 ∫ z 2 dA = M y
A A A

(f)
c1 ∫ y dA + c2 ∫ y 2 dA + c3 ∫ yz dA = − M z
A A A

For the origin of the y and z axes to be coincident with the centroid of the section, it is
required that
(g)
∫A y dA = ∫A z dA = 0
Based on Eq. (d), we conclude that c1 = 0 ; based on Eqs. (5.11), we conclude that σ x = 0
at the origin. Thus, the neutral axis passes through the centroid, as in the beam of sym-
metrical section. In addition, the field of stress described by Eqs. (5.11) satisfies the
equations of equilibrium and compatibility and the lateral surfaces are free of stress. Now
consider the defining relationships

I y = ∫ z 2 dA, I z = ∫ y 2 dA, I yz = ∫ yz dA (5.12)

A A A

5.3   Pure Bending of Beams of Asymmetrical Cross Section 247

 The quantities I y and I z are the moments of inertia about the y and z axes, respectively,
and I yz is the product of inertia about the y and z axes. From Eqs. (e) and (f), together with
Eqs. (5.12), we obtain expressions for c2 and c3.

5.3.1  Stress Distribution

Substitution of the constants into Eqs. (5.11) results in the following generalized flexure
formula:
( M y I z + M z I yz ) z − ( M y I yz + M z I y ) y
σx = 2 (5.13)
I y I z − I yz

The equation of the neutral axis is found by equating this expression to zero:
( M y I z + M z I yz ) z − ( M y I yz + M z I y ) y = 0 (5.14)

This is an inclined line through the centroid C. The angle φ between the neutral axis and
the z axis is determined as follows:
y M y I z + M z I yz
tan θ = = (5.15)
z M y I yz + M z I y
The angle φ (measured from the z axis) is positive in the clockwise direction, as shown
in Fig. 5.3. The highest bending stress occurs at a point located farthest from the
­neutral axis.
There is a specific orientation of the y and z axes for which the product of inertia I yz
vanishes. Labeling the axes so oriented as y ′ and z ′ , we have I y ′z ′ = 0. The flexure for-
mula under these circumstances becomes
M y′ z ′ M z′ y ′
σx= − (5.16)
Iy I z′

The y′ and z′ axes now coincide with the principal axes of inertia of the cross section,
and we can find the stresses at any point by applying Eq. (5.13) or Eq. (5.16).
The kinematic relationships discussed in Section 5.2 are valid for beams of asym-
metrical section provided that y and z represent the principal axes.

5.3.2  Transformation of Inertia Moments

Recall that the two-dimensional stress (or strain) and the moment of inertia of an area
are second-order tensors (Section 1.17). Thus, the transformation equations for stress
and moment of inertia are analogous (Section C.2.2). In turn, the Mohr’s circle analysis
and all conclusions drawn for stress apply to the moment of inertia. With reference to the
coordinate axes shown in Fig. 5.3, applying Eq. (C.12a), the moment of inertia about the
y ′ axis is
I y + Iz I y − Iz
I y′ = + cos 2θ − I yz sin 2θ (5.17)
2 2

248 Chapter  5   Bending of Beams

From Eq. (C.13), the orientation of the principal axes is given by
2 I yz
tan 2θ p = − (5.18)
I y − Iz

The principal moments of inertia, I1 and I 2 , from Eq. (C.14) are

2
Ix + Iy  Ix − Iy  2
I1,2 = ±  2  + I xy (5.19)
2

where the subscripts 1 and 2 refer to the maximum and minimum values, respectively.
Determination of the moments of inertia and stresses in an asymmetrical section is
illustrated in Example 5.1.

EXAMPLE 5.1  Analysis of an Angle in Pure Bending

A 150- by 150-mm slender angle of 20-mm thickness is subjected to oppo-
sitely directed end couples M z = 11 kN ⋅ m at the centroid of the cross section.
What bending stresses exist at points A and B on a section away from the ends
(Fig. 5.4a)? Determine the orientation of the neutral axis.
Solution  Equations (5.13) and (5.16) are applied to obtain the normal stress.
This requires first determining a number of section properties through the use
of familiar expressions of mechanics given in Appendix C. The computer pro-
gram presented in Table C.2 provides a check of the numerical values obtained
here for the area characteristics and may be extended to compute the stresses.
Location of the Centroid C. Let y and z represent the distances from C to
arbitrary reference lines (denoted as Z and Y ):
Σ Ai zi A1 z1 + A2 z2 130 × 20 × 10 + 150 × 20 × 75
z= = = = 45 mm
Σ Ai A1 + A2 130 × 20 + 150 × 20

150 mm
z = 45 mm yB′ z′
C 75 mm
B 2 Z B
y
C A2 45° 33.72°
10 mm C
z z
C1
150 mm y N.A.

′
zA
20 mm 65 mm
A A1 yA′ y′
y
Y 10 mm A

(a) (b)

Figure 5.4.  Example 5.1. An equal-leg-angle cross section of beam.

5.3   Pure Bending of Beams of Asymmetrical Cross Section 249

 where zi represents the z distance from the Y reference line to the centroid of
each subarea (A1 and A2) composing the total cross section. Since the section
is symmetrical, z = y.
Moments and Products of Inertia. For a rectangular section of depth h and
width b, the moment of inertia about the neutral z axis is I z = bh3 /12 (Table
C.1). We now use the yz axes as reference axes through C. Representing the
distances from C to the centroids of each subarea by d y1 , d y 2, dz1, and dz2 we
obtain the moments of inertia with respect to these axes using the parallel-axis
theorem. Applying Eq. (C.9),

I z = ∑ ( I z + Ad y2 ) = I z 1 + A1d y21 + I z 2 + A2 d y22

Thus, referring to Fig. 5.4a,

I z = I y = 121 × 20 × (130)3 + 130 × 20 × (40)2

+ 121 × 150 × (20)3 +150 × 20 × (35)2
= 11.596 × 106 mm 4

The transfer formula, Eq. (C.11), for a product of inertia yields

I yz = ∑ ( I yz + Ad y d z )
= 0 + 130 × 20 × 40 × ( −35) + 0 + 150 × 20 × ( −35) × 30
= −6.79 × 106 mm 4

Stresses Using Eq. (5.13). We have y A = 0.105 m, y B = −0.045 m, z A = −0.045 m, z B = −0.045 m, and
A = 0.105 m, y B = −0.045 m, z A = −0.045 m, z B = −0.045 m, and M y = 0. Hence,

M z ( I yz z A − I y y A )
(σ x ) A = 2
I y I z − I yz

11(103 ) [ ( − 6.79)( − 0.045) − (11.596)(0.105) ]

= = − 114 MPa (h)
(11.596)2 − ( −6.79)2  10−6

Similarly,

11(103 ) [ ( −6.79)( −0.045) − (11.596)( −0.045) ]

(σ x ) B = = 103 MPa
(11.596)2 − ( −679)2 10−6

Alternatively, these stresses may be calculated as described next.

250 Chapter  5   Bending of Beams

 Directions of the Principal Axes and the Principal Moments of Inertia.
Employing Eq. (5.18), we have

−2( −6.79)
tan 2θ p = = ∞, 2θ p = 90° or 270°
11.596 − 11.596

Therefore, the two values of θ p are 45° and 135°. Substituting the first of these
values into Eq. (5.17), we obtain I y ′ = [11.596 + 6.79sin 90°]. Since the princi-
pal moments of inertia are, by application of Eq. (5.19),

I1,2 = [11.596 ± 0 + 6.792 ]106 = [11.596 ± 6.79 ]106

we see that I1 = I y ′ = 18.386 × 106 mm 4 and I 2 = I z ′ = 4.806 × 106 mm 4 . The

principal axes are indicated in Fig. 5.4b as the y′ z′ axes.
Stresses Using Eq. (5.16). The components of bending moment about the
principal axes are

M y ′ = 11(103 ) sin 45° = 7778 N ⋅ m

M z ′ = 11(103 ) cos 45° = 7778 N ⋅ m

Equation (5.16) is now applied, referring to Fig 5.4b, with y A′ = 0.043 m,

z A′ = −0.106 m, y B′ = −0.0636 m, and z B′ = 0 determined from geometric
considerations:
7778( −0.106) 7778(0.043)
(σ x ) A = − = −114 MPa
18.386 × 10−6 4806 × 10−6

7778(−0.0636)
(σ x ) B = 0 − = 103 MPa
4.806 × 10−6
as before.
Direction of the Neutral Axis. From Eq. (5.15), with M y = 0,
I yz −6.79
tan φ = or φ = arctan = −30.4°
Iy 11.596

The negative sign indicates that the neutral is located counterclockwise from
the z axis (Fig. 5.4b).

5.4  BENDING OF A CANTILEVER OF NARROW SECTION

Consider a narrow cantilever beam of rectangular cross section, loaded at its free end by
a concentrated force of magnitude such that the beam weight may be neglected (Fig. 5.5).
This situation may be regarded as a case of plane stress provided that the beam thickness t

5.4   Bending of a Cantilever of Narrow Section 251

 Figure 5.5.  D
 eflections of an end- L
loaded cantilever beam.
2h P x

is small relative to the beam depth 2h. The distribution of stress in the beam, as we found
in Example 3.1, is given by

Px P
σ x = −   y, σ y = 0, τ xy = − (h2 − y 2 ) (3.21)
 I  2I
To derive expressions for the beam displacement, we must relate stress, described by
Eq. (3.21), to strain. This is accomplished through the use of the strain-displacement rela-
tions and Hooke’s law:
∂u Pxy ∂υ ν Pxy
=− , = (a, b)
∂x EI ∂y EI
∂u ∂υ 2(1 + ν )τ xy (1 + ν ) P 2
+ = =− ( h − y 2 ) (c)
∂ y ∂x E EI

Integration of Eqs. (a) and (b) yields

Px 2 y
u=− + u1 ( y ) (d)
2 EI
ν Pxy 2
υ= + υ1 ( x ) (e)
2 EI
Differentiating Eqs. (d) and (e) with respect to y and x, respectively, and substituting into
Eq. (c), we have
du1 P dυ P 2 (1 + ν ) Ph2
− (2 + ν ) y 2 = − 1 + x −
dy 2 EI dx 2 EI EI
In this expression, note that the left and right sides depend only on y and x, respec-
tively. These variables are independent of each other, so we conclude that the equation can
be valid only if each side is equal to the same constant:
du1 P dυ1 Px 2 (1 + ν ) Ph2
− (2 + ν ) y 2 = a1 , − + = − a1
dy 2 EI dx 2 EI EI
These are integrated to yield
P
u1 ( y ) = (2 + ν ) y 3 + a1 y + a2
6 EI
Px 3 (1 + ν ) Pxh2
υ1 ( x ) = − − a1 x + a3
6 EI EI

252 Chapter  5   Bending of Beams

in which a2 and a3 are constants of integration. The displacements may now be written
Px 2 y P
u=− + (2 + ν ) y 3 + a1 y + a2 (5.20a)
2 EI 6 EI
ν Pxy 2 Px 3 (1 + ν ) Pxh2
υ= + − − a1 x + a2 (5.20b)
2 EI 6 EI EI
The constants a1, a2 , and a3 depend on known conditions. If, for example, the situa-
tion at the fixed end is such that
∂u
= 0, υ = u = 0 at x = L, y = 0
∂x
then, from Eqs. (5.20),
PL2 PL3 PLh2 (1 + ν )
a1 = , a2 = 0, a3 = +
2 EI 3EI EI
The beam displacement is therefore
( 2 + v ) Py 3
u=
p
2 EI
(
L2 − x 2 y + ) 6 EI
(5.21)

P  x 3 L3 x 
υ= + + (ν y 2 − L2 ) + h2 (1 + ν )( L − x )  (5.22)
EI  6 3 2 
On examining these equations, it becomes that u and υ do not obey a simple linear
­relationship with y and x. We conclude that plane sections do not, as assumed in elemen-
tary theory, remain plane subsequent to bending.

5.4.1  Comparison of the Results with the Elementary Theory Results

The vertical displacement of the beam axis is obtained by substituting y = 0 into Eq. (5.22):

Px 3 PL2 x PL3 Ph2 (1 + ν )

(υ ) y =0 = − + + (L − x ) (5.23)
6EI 2EI 3EI EI
Introducing this relation into Eq. (5.7), the radius of curvature is given by
1 Px M
≈ =
rx EI EI

provided that dυ /dx is a small quantity. Once again, we obtain Eq. (5.9a), the beam curvature–
moment relationship of elementary bending theory.
It is also a simple matter to compare the total vertical deflection at the free end ( x = 0)
with the deflection derived in elementary theory. Substituting x = 0 into Eq. (5.23), the
total deflection is
PL3 Ph2 (1 + ν ) L PL3 Ph2 L
(υ ) x = y =0 = + = + (5.24)
3EI EI 3EI 2GI

5.4   Bending of a Cantilever of Narrow Section 253

 where the deflection associated with shear is clearly Ph2 L/2GI = 3PL/2GA. The ratio
of the shear deflection to the bending deflection at x = 0 provides a measure of beam
slenderness:
2 2
Ph2 L /2GI 3 h2 E 3  2h  ≈  2h 
= = (1+ν )    
PL3 /3EI 2 L2G 4 L L
If, for example, L = 10(2h), the preceding quotient is only 100
1
. For a slender beam, 2h  L
and the deflection is mainly due to bending. In contrast, in cases involving vibration at
higher modes, and in wave propagation, the effect of shear is of great importance in slen-
der as well as in other beams.
In the case of wide beams (t  2h), Eq. (5.24) must be modified by replacing E and ν
as indicated in Table 3.1.

5.5  BENDING OF A SIMPLY SUPPORTED NARROW BEAM

In this section, we consider the stress distribution in a narrow beam of thickness t and
depth 2h subjected to a uniformly distributed loading (Fig. 5.6). The situation described
here is one of plane stress, subject to the following boundary conditions, consistent with
the origin of an x, y coordinate system located at midspan and midheight of the beam, as
shown:

(τ xy ) y =± h = 0, (σ y ) y =+ h = 0, (σ y ) y =− h = − p/t (a)

Since no longitudinal load is applied at the ends, it would appear reasonable to state that
σ x = 0 at x = ± L. However, this boundary condition leads to a complicated solution, and a
less severe statement is used instead:
h (b)
∫− hσ x t dy = 0
The corresponding condition for bending couples at x = ± L is
h
∫− hσ x ty dy = 0 (c)
For y equilibrium, it is required that
h
∫− h τ xy t dy = ± pL x = ± L (d)

Figure 5.6.  B
 ending of a simply supported p
beam with a uniform load.
2h x

L L
y

254 Chapter  5   Bending of Beams

5.5.1  Use of Stress Functions
The problem is treated by superimposing the solutions Φ2 , Φ3, and Φ5 (Section 3.6) with
c2 = b2 = a3 = c3 = a5 = b5 = c5 = e5 = 0
We then have
a b d d 2d
Φ = Φ 2 + Φ 3 + Φ5 = 2 x 2 + 3 x 2 y + 3 y 3 + 5 x 2 y 3 − 5 y5
2 2 6 6 60
The stresses are
σ x = d3 y + d5 ( x 2 y − 23 y 3 )
d
σ y = a2 + b3 y + 5 y 3 (e)
3
τ xy = −b3X − d5 xy 2
The conditions (a) are
−b3 − d5 h2 = 0
d5 3
a2 + b3h + h =0
3
d p
a2 − b3h − 5 h3 = −
3 t
and the solution is
p 3p 3p
a2 = − , b3 = , d5 = −
2t 4th 4th3
The constant d3 is obtained from condition (c) as follows:
h   2 3 

∫− h  d3 y + d5  L y − 3 y yt dy = 0
2
 
or
2 3p  L2 2 
d3 = − d5  L2 − h2  = −
 5  4th  h2 5 
Expressions (e), together with the values obtained for the constants, also fulfill conditions
(b) and (d).
The state of stress is thus represented by

py 2 py  y 2 h2 
σx = (L − x2 ) + −  (5.25a)
2I I  3 5

− p  y3 2 h3 
σy = − h 2
y + (5.25b)
2 I  3 3 
− px 2
τ xy =
2I
(
h − y2 ) (5.25c)

5.5   Bending of a Simply Supported Narrow Beam 255

 where I = 23 th3 is the area moment of inertia taken about a line through the centroid, par-
allel to the z axis. Although the solutions given by Eqs. (5.25) satisfy the equations of
elasticity and the boundary conditions, they are not exact. We can see this by substitut-
ing x = ± L into Eq. (5.25a) to obtain the following expression for the normal distributed
forces per unit area at the ends:

py  y 2 h2 
px = − 
I  3 5

This state cannot exist, as no forces act at the ends. From Saint-Venant’s principle, how-
ever, we may conclude that the solutions do predict the correct stresses throughout the
beam, except near the supports.

5.5.2  Comparison of the Results with the Elementary Theory Results

Recall that the longitudinal normal stress derived from elementary beam theory is
σ x = −My/ I ; this is equivalent to the first term of Eq. (5.25a). The second term is then the
difference between the longitudinal stress results given by the two approaches. To gauge
the magnitude of the deviation, consider the ratio of the second term of Eq. (5.25a) to the
result of elementary theory at x = 0. At this point, the bending moment is a maximum.
Substituting y = h for the condition of maximum stress, we obtain

∆σ x ( ph/ I )( h2 /3 − h2 /5) 4  h  2
= =
(σ x )elem. theory phL2 /2 I 15  L 

For a beam of length 10 times its depth, this ratio is small, 1500 1 . For beams of ordinary

proportions, we can conclude that elementary theory provides a result of sufficient accu-
racy for σ x . On the one hand, for σ y, this stress is not found in the elementary theory. On
the other hand, the result for τ xy is the same as that yielded by elementary beam theory.
The displacement of the beam may be determined in a manner similar to that
described for a cantilever beam (Section 5.4).

Part B: Approximate Solutions

5.6  ELEMENTARY THEORY OF BENDING

We may conclude, on the basis of the previous sections, that exact solutions are difficult
to obtain. We also observed that for a slender beam, the results of the exact theory do not
differ markedly from those found with the mechanics of materials or elementary approach
provided that solutions close to the ends are not required. The bending deflection is very
much larger than the shear deflection, so the stress associated with the former predomi-
nates. As a consequence, the normal strain ε y resulting from transverse loading may be
neglected.

256 Chapter  5   Bending of Beams

 Because it is more easily applied, the elementary approach is usually preferred in
engineering practice. The exact and elementary theories should be regarded as comple-
mentary—rather than competitive—approaches, enabling the analyst to obtain the degree
of accuracy required in the context of the specific problem at hand.

5.6.1  Assumptions of Elementary Theory

The basic presuppositions of the elementary theory [Ref. 5.2], for a slender beam whose
cross section is symmetrical about the vertical plane of loading, are
∂υ ∂u ∂υ
εy = = 0, γ xy = + =0
∂y ∂y ∂x
(5.26)
σx
ε x =   (independent of z)
E
ε z = 0, γ yz = γ xz = 0 (5.27)
The first equation of Eqs. (5.26) is equivalent to the assertion υ = υ ( x ). Thus, all points
in a beam at a given longitudinal location x experience identical deformation. The second
equation of Eqs. (5.26), together with υ = υ ( x ), yields, after integration,
dυ
u = −y + u0 (x ) (a)
dx
The third equation of Eqs. (5.26) and Eqs. (5.27) imply that the beam is considered nar-
row, and we have a case of plane stress.
At y = 0, the bending deformation should vanish. Referring to Eq. (a), it is clear that
uo ( x ) must represent axial deformation. The term dυ/dx is the slope θ  of the beam axis, as
shown in Fig. 5.7a, and is very much smaller than unity. Therefore,
dυ
u = −y = − yθ
dx
The slope is positive when clockwise, provided that the x and y axes have the directions
shown. Since u is a linear function of y, this equation restates the kinematic hypothesis of
the elementary theory of bending: Plane sections perpendicular to the longitudinal axis
of the beam remain plane subsequent to bending. This assumption is confirmed by the
exact theory only in the case of pure bending.
p
x

v V V + dV
x
y dv
=θ
y dx M M + dM
dv
dx
u = −y y
dx
Beam axis
(a) (b)

Figure 5.7.  (a) Longitudinal displacements in a beam due to rotation of a

plane section; (b) element between adjoining sections of a beam.

5.6   Elementary Theory of Bending 257

 5.6.2  Method of Integration
In the next section, we obtain the stress distribution in a beam according to the elemen-
tary theory. We now derive some useful relations involving the shear force V, the bending
moment M, the load per unit length p, the slope θ , and the deflection. Consider a beam
element of length dx subjected to a distributed loading (Fig. 5.7b). Since dx is small, we
omit the variation in the load per unit length p. In the free-body diagram, all the forces
and the moments are positive. The shear force obeys the sign convention discussed in
Section 1.4, while the bending moment is in agreement with the convention adopted in
Section 5.2.
In general, the shear force and bending moment vary with the distance x, such that
these quantities will have different values on each face of the element. The increments in
shear force and bending moment are denoted by dV and dM, respectively. The equilibrium
of forces in the vertical direction is governed by V − (V + dV ) − p dx = 0 or
dV
= − p (5.28)
dx
That is, the rate of change of shear force with respect to x is equal to the algebraic value of
the distributed loading. Equilibrium of the moments about a z axis through the left end of
the element, neglecting the higher-order infinitesimals, leads to
dM
= − V (5.29)
dx
This relation states that the rate of change of bending moment is equal to the algebraic
value of the shear force—a relation that is valid only if a distributed load or no load acts
on the beam segment. Combining Eqs. (5.28) and (5.29), we have
d2M
= p (5.30)
dx 2
The basic equation of bending of a beam, Eq. (5.10), combined with Eq. (5.30), may now
be written as
d 2  d 2υ 
EI = p (5.31)
dx 2  dx 2 
For a beam of constant flexural rigidity EI, the beam equations derived here may be
expressed as
d 4υ
EI = EIυ IV = p
dx 4
d 3υ
EI 3 = EIυ ′′′ = −V
dx (5.32)
d 2υ
EI 2 = EIυ ′′ = M
dx
dυ
EI = EIυ ′ = ∫ Mdx
dx

258 Chapter  5   Bending of Beams

These relationships also apply to wide beams provided that we substitute E/(1 − ν 2 ) for E
(Table 3.1).
In many problems of practical importance, the deflection due to transverse loading of
a beam may be obtained through successive integration of the beam equation:

EIυ IV = p
x
EIυ ′′′ = ∫ p dx + c1
0
x x
EIυ ′′ = ∫ dx ∫ p dx + c1 x + c2 (5.33)
0 0
x x x
EIυ ′ = ∫ dx ∫ dx ∫ p dx + 12 c1 x 2 + c2 x + c3
0 0 0
x x x x
EIυ = ∫ dx ∫ dx ∫ dx ∫ p dx + 61 c1 x 3 + 12 c2 x 2 + c3 x + c4
0 0 0 0

Alternatively, we could begin with EIυ ′′ = M ( x ) and integrate twice to obtain

x x
EIυ = ∫ dx ∫ M dx + c3 x + c4 (5.34)
0 0

In either case, the constants c1 , c2 , c3 and c4 , which correspond to the homogeneous

solution of the differential equations, may be evaluated from the boundary conditions.
The constants c1 , c2 , c3 / EI , and c4 / EI represent the values at the origin of V, M, θ , and
υ, respectively. In the method of successive integration, there is no need to distinguish
between statically determinate and statically indeterminate systems (Section 5.11),
because the equilibrium equations represent only two of the boundary conditions (on the
first two integrals), and because the total number of boundary conditions is always equal
to the total number of unknowns.

EXAMPLE 5.2  Displacements of a Cantilever Beam

A cantilever beam AB of length L and constant flexural rigidity EI carries a
moment M o at its free end A (Fig. 5.8a). Derive the equation of the deflection
curve and determine the slope and deflection at A.
Solution  From the free-body diagram of Fig. 5.8b, observe that the bending
moment is +M o throughout the beam. Thus, the third of Eqs. (5.32) becomes
EIυ′′ = M o

A B A
x
Mo Mo O M
L x
y
(a) (b)

Figure 5.8.  E
 xample 5.2. (a) A cantilever beam is subjected to moment at its
free end; (b) free-body diagram of part AO.

5.6   Elementary Theory of Bending 259

 Integration yields
EIυ ′ = M o x + c1
The constant of integration c1 can be found from the condition that the slope
is zero at the support; therefore, we have υ ′( L) = 0 from which c1 = − M o L.
The slope is then
Mo
υ′ = ( x − L) (5.35)
EI
Integrating, we obtain
Mo 2
υ= ( x − 2Lx ) + c2
2EI
The boundary condition on the deflection at the support is υ ( L ) = 0, which
yields c2 = M o L2 /2EI . The equation of the deflection curve is thus a parabola:
Mo 2
υ= ( L + x 2 − 2Lx ) (5.36)
2EI
However, every element of the beam experiences the same moments and defor-
mation. The deflection curve should, therefore, be part of a circle. This inconsis-
tency results from the use of an approximation for the curvature, Eq. (5.7). The
error is very small, however, when the deformation υ is small [Ref. 5.1].
The slope and deflection at A are readily found by letting x = 0 in Eqs.
(5.35) and (5.36):

MoL M L2
θA = − , υA = o (5.37)
EI 2EI

The minus sign indicates that the angle of rotation is counterclockwise

(Fig. 5.8a).

5.7  NORMAL AND SHEAR STRESSES

When a beam is bent by transverse loads, usually both a bending moment M and a shear
force V act on each cross section. The distribution of the normal stress associated with the
bending moment is given by the flexure formula, Eq. (5.4):
My
σx = − (5.38)
I
where M and I are taken with respect to the z axis (Fig. 5.7).
In accordance with the assumptions of elementary bending, Eqs. (5.26) and (5.27),
we omit the contribution of the shear strains to beam deformation in these calculations.
However, shear stresses do exist, and the shearing forces are the resultant of the stresses.

260 Chapter  5   Bending of Beams

The shearing stress τ xy acting at section mn, which is assumed to be uniformly distributed
over the area b ⋅ dx, can be determined on the basis of equilibrium of forces acting on the
shaded part of the beam element (Fig. 5.9). Here b is the width of the beam a distance
y1 from the neutral axis, and dx is the length of the element. The distribution of normal
stresses produced by M and M + dM is indicated in the figure. The normal force distrib-
uted over the left face mr on the shaded area A* is equal to
My
dA
b /2 h1
∫− b /2 ∫ y σ x dy dz = ∫A −
1
*
I
Similarly, an expression for the normal force on the right face ns may be written in
terms of M + dM. The equilibrium of x-directed forces acting on the beam element is gov-
erned by
( M + dM ) y My
−∫ dA − ∫ * − dA = τ xy b dx
*
A I A I
from which we have
1 dM
τ xy = −
Ib ∫ A*
dx
y dA

After substituting in Eq. (5.29), we obtain the shear formula (also called the shear
stress formula) for beams:
V VQ
τ xy =
Ib ∫ A*
y dA =
Ib
(5.39)

The integral represented by Q is the first moment of the shaded area A* with respect to the
neutral axis z:
Q = ∫ * y dA = A* y (5.40)
A

By definition, y is the distance from the neutral axis to the centroid of A*. In the case
of sections of regular geometry, A* y provides a convenient means of calculating Q. The
shear force acting across the width of the beam per unit length
VQ
q = τ xy b = (a)
I
is called the shear flow.

h2
C
τxy x z
M M +dM y y1
m n h1 A*
r s
y dx b

Figure 5.9.  (a) Beam segment for analyzing shear stress; (b) cross
section of beam.

5.7   Normal and Shear Stresses 261

 5.7.1  Rectangular Cross Section
In the case of a rectangular cross section of width b and depth 2h, the shear stress at y1 is

V b /2 h
Ib ∫− b /2 ∫ y1
τ xy = y dy dz =
V 2
2I
h − y12 (5.41) ( )
This equation shows that the shear stress varies parabolically with y1. It is zero when
y1 = ± h, and has its maximum value at the neutral axis:
Vh2 3 V
τ max =
= (5.42)
2I 2A
where A = 2bh is the area of the rectangular cross section. Note that the maximum shear
stress (either horizontal or vertical: τ xy = τ yx) is 1.5 times larger than the average shear
stress V/A. As observed in Section 5.4, for a thin rectangular beam, Eq. (5.42) is the exact
distribution of shear stress. More generally, for wide rectangular sections and for other
sections, Eq. (5.39) yields only approximate values of the shearing stress.

5.7.2  Various Cross Sections

Because the shear formula for beams is based on the flexure formula, the limitations of the
bending formula apply when it is used. Problems involving various types of cross sections
can be solved by following procedures identical to that for rectangular sections. Table 5.1

Table 5.1.  Maximum Shearing Stress for Some Typical Beam Cross-Sectional Forms
Cross Section Maximum Shearing Stress Location
Rectangle V 3V NA
τ max =
2A
NA

Circle V 4V NA
τ max =
3A
NA

Hollow Circle V V NA
τ max = 2
A
NA

Triangle V 3V Halfway between top and

τ max = bottom
2A
NA

Diamond V 9V At h/8 above and below

τ max = the NA
8A
h
NA

Notes: A, cross-sectional area; V, transverse shear force; NA, the neutral axis.

262 Chapter  5   Bending of Beams

shows some typical cases. Observe that shear stress can always be expressed as a constant
times the average shear stress (P/A), where the constant is a function of the cross-sectional
form. Nevertheless, the maximum shear stress does not always occur at the neutral axis. For
instance, in the case of a cross section having nonparallel sides, such as a triangular section,
the maximum value of Q/b (and thus τ xy) occurs at midheight, h/2, while the neutral axis is
located at a distance h/3 from the base.
The following examples illustrate the application of the normal and shear stress
formulas.

EXAMPLE 5.3  Stresses in a Beam of T-Shaped Cross Section

A simply supported beam of length L carries a concentrated load P
(Fig. 5.10a). Find: (a) The maximum shear stress, the shear flow qj, and the
shear stress T j in the joint between the flange and the web; (b) the maximum
bending stress. Given: P = 5 kN and L = 4 m.
Solution  The distance y from the Z axis to the centroid is obtained as fol-
lows (Fig. 5.10a):

A1 y1 + A2 y2 20 ( 60 ) 70 + 60 ( 20 ) 30
y= = = 50 mm
A1 + A2 20 ( 60 ) + 60 ( 20 )

The moment of inertia I about the NA is determined by the parallel axis theorem:

I = 121 ( 60 )( 20 )3 + 20 ( 60 )( 20 )2 + 121 ( 20 )( 60 )3 + 20 ( 60 )( 20 )2 = 136 × 104 mm 4

The shear and moment diagrams (Figs. 5.10b and c) are sketched by applying
the method of sections.
a. The maximum shearing stress in the beam takes place at the NA on the
cross section supporting the largest shear force V. Consequently, QNA = 50 ( 20 ) 25 = 25 × 103 mm 3
QNA = 50 ( 20 ) 25 = 25 × 103 mm 3. The shear force equals 2.5 kN on all cross sec-
tions of the beam (Fig. 5.10b), Thus,

τ max =
3
(
Vmax QNA 2.5 × 10 25 × 10
=
−6
)
= 2.3 MPa
lb 136 × 10−8 ( 0.02 )

The first moment of the area of the flange about the NA is

Q f = 20 ( 60 ) 20 = 24 × 103 mm 3 .

Shear flow and shear stress in the joint are

qj =
VQ f
=
(
2.5 × 103 24 × 10−6 ) = 44.1 kN/m
−8
I 136 × 10

and t j = qj/b = 44.1/0.02 = 2.205 MPa.

5.7   Normal and Shear Stresses 263

 60 mm
y
y P
20 mm A1
N. A.
A z
B x
C 60 mm
A2 y = 50 mm
P L L P
2
2 2 2
Z
20 mm
(a)

P
2
= 2.5 kN

x
P
2

(b)

M PL
= 5 kN·m
4

x
(c)

Figure 5.10.  E
 xample 5.3. (a) Load diagram and beam cross section; (b) shear dia-
gram; (c) moment diagram.

b. The largest moment takes place at midspan (Fig. 5.10c). Equation (5.39) is
therefore

Mc 5 × 103 ( 0.05)
σ max = = = 183.8 MPa
I 136 × 10−8

EXAMPLE 5.4  Shear Stresses in a Flanged Beam

A cantilever wide-flange beam is loaded by a force P at the free end act-
ing through the centroid of the section. The beam is of constant thickness t
(Fig. 5.11a). Determine the shear stress distribution in the section.
Solution  The vertical shear force at every section is P. It is assumed that the
shear stress τ xy is uniformly distributed over the web thickness. Then, in the
web, for 0 ≤ y1 ≤ h1 and applying Eq. (5.39),

V * P h − h1   h1 − y1  
τ xy = A y = b( h − h1 )  h1 +  + t ( h1 − y1 )  y1 + 
Ib It   2  2  

264 Chapter  5   Bending of Beams

 b τmax Figure 5.11.  E
 xample 5.4. (a)
Cross ­ section of a
wide-flange beam;
h τmin
h1 (b) shearing stress
P z
distribution.
C y1
h h1
t

y τavg
(a) (b)

This equation may be written as

τ xy =
P
( ) ( )
b h2 − h12 + t h12 − y12  (b)
2It  

The shearing stress thus varies parabolically in the web (Fig. 5.11b). The
extreme values of τ xy found at y1 = 0 and y1 = ± h1 are, from Eq. (b), as
follows:
P Pb 2
τ max = (bh2 − bh12 + th12 ), τ min = ( h − h12 )
2It 2It

Usually t  b, so that the maximum and minimum stresses do not differ

appreciably, as is seen in the figure. Similarly, the shear stress in the flange,
for h1 < y1 ≤ h, is
P h − y1   P 2
τ xz = b( h − y1 )  y1 + = ( h − y12 ) (c)
Ib   2   2I

This is the parabolic equation for the variation of stress in the flange, shown
by the dashed lines in Fig. 5.11b.
Comments  Clearly, for a thin flange, the shear stress is very small as
compared with the shear stress in the web. As a consequence, the approxi-
mate average value of shear stress in the beam may be found by dividing
P by the web cross section, with the web height assumed to be equal to the
beam’s overall height: τ avg = P/2. This result is indicated by the dotted lines
in Fig. 5.11b.
The distribution of stress given by Eq. (c) is fictitious, because the inner
planes of the flanges must be free of shearing stress, as they are load-free
boundaries of the beam. This contradiction cannot be resolved by the elemen-
tary theory; instead, the theory of elasticity must be applied to obtain the cor-
rect solution. Fortunately, this defect of the shearing stress formula does not
lead to serious error since, as pointed out previously, the web carries almost
all the shear force. To reduce the stress concentration at the juncture of the
web and the flange, the sharp corners should be rounded.

5.7   Normal and Shear Stresses 265

 EXAMPLE 5.5  Beam of Circular Cross Section
A cantilever beam of circular cross section supports a concentrated load P at
its free end (Fig. 5.12a). The shear force V in this beam is constant and equal
to the magnitude of the load P = V. Find the maximum shearing stresses (a) in
a solid cross section and (b) in a hollow cross section.
Assumptions: All shear stresses do not act parallel to the y axis. At a point
such as a or b on the boundary of the cross section, the shear stress τ must act
parallel to the boundary. The shear stresses at line ab across the cross section
are not parallel to the y axis and cannot be determined by the shear formula,
τ = VQ/Ib. The maximum shear stresses occur along the neutral axis z, are
uniformly distributed, and act parallel to the y axis. These stresses are within
approximately 5% of their true value [Ref. 5.1].
Solution
a. Solid Cross Section (Fig. 5.12b). The shear formula may be used to cal-
culate, with reasonable accuracy, the shear stresses at the neutral axis. The
area properties for a circular cross section of radius c (see Table C.1) are
π c4 π c 2  4c  2 c 2
I= Q = A* y = = (d)
4 2  3π  3
and b = 2c. The maximum shear stress is thus
VQ 4V 4V
τ max == 2 = (e)
Ib 3π c 3A
where A is the cross-sectional area of the beam.
Comment  This result shows that the largest shear stress in a circular beam is
4/3 times the average shear stress τ avg = V /A.
b. Hollow Circular Cross Section (Fig. 5.11c). Equation (5.43) applies with
equal rigor to circular tubes, since the same assumptions stated in the fore-
going are valid. But in this case, by Eq. (C.3), we have
Q = 23 ( c23 − c13 ), b = 2( c2 − c1 ), A = π ( c22 − c12 )

y y y

a b
τ
P r
z d e z
z N.A. c1 c2
τ max
x
2c
2c
(a) (b) (c)

Figure 5.12.  E
 xample 5.5. (a) A cantilever beam under a load P; (b) shear stress distribution on
a circular cross section; (c) hollow circular cross section.

266 Chapter  5   Bending of Beams

 and
π 4 4

4
I= (
c2 − c1 )
So, the maximum shear stress may be written in the following form:
VQ 4V c22 + c2 c1 + c12
τ max = = (f)
Ib 3 A c22 + c12
Comment  Note that for c1 = 0, Eq. (e) reduces to Eq. (5.43) for a solid circular
beam, as expected. In the special case of a thin-walled tube, we have r/t > 10,
where r and t represent the mean radius and thickness, respectively. As a theoreti-
cal limiting case, setting c2 = c1 = r means that Eq. (e) results in τ max = 2V / A.

5.7.3  Beam of Constant Strength

When a beam is stressed to a constant permissible stress, σ all throughout, then clearly the
beam material is used to its greatest capacity. For a given material, such a design is of
minimum weight. At any cross section, the required section modulus S is defined as
M
s= (5.43)
σ all
where M presents the bending moment on an arbitrary section. Tapered beams designed in
this way are called beams of constant strength. Ultimately, the shear stress at those beam
locations where the moment is small controls the design.
Examples of beams with uniform strength include leaf springs and certain cast
machine elements. For a structural member, fabrication and design constraints make it
impractical to produce a beam of constant stress. Hence, welded cover plates are often
used for parts of prismatic beams where the moment is large—for instance, in a bridge
girder. When the angle between the sides of a tapered beam is small, the flexure formula
allows for little error. On the contrary, the results found by applying the shear stress for-
mula may not be accurate enough for nonprismatic beams. Often, a modified form of
this formula is used for design purposes. The exact distribution in a rectangular wedge is
found by applying the theory of elasticity (Section 3.10).

EXAMPLE 5.6  Design of a Uniform Strength Beam

A cantilever beam of constant strength and rectangular cross section is to
carry a concentrated load P at the free end (Fig. 5.13a). Find the required
cross-sectional area for two cases: (a) the width b is constant; (b) the height h
is constant.
Solution
a. At a distance x from A, M = Px and S = bh2/6. By applying Eq. (5.43), we have
bh2 Px
= (g)
6 σ all

5.7   Normal and Shear Stresses 267

 Figure 5.13.  E
 xample 5.6. (a) Constant- P
strength cantilever; (b) side
view; (c) top view.

L
(a)

h
A h1

x B
(b)

P b b1

x
(c)

Likewise, at a fixed end (x = L and h = h1),

bh12 PL
=
6 σ all
Dividing Eq. (g) by the preceding relationship leads to
x
h = h1 (h)
L
Thus, the depth of the beam varies parabolically from the free end (Fig. 5.13b).
b. Equation (g) now gives
 6P  b
b= 2  x = 1 x (i)
 h σ all  L
Comments  In Eq. (i), the term in parentheses represents a constant and is
set equal to b1/L; hence, when x = L, the width is b1 (Fig. 5.13c). Clearly,
the cross section of the beam near end A must be designed to resist the shear
force, as depicted by the dashed lines in the figure.

5.8  EFFECT OF TRANSVERSE NORMAL STRESS

When a beam is subjected to a transverse load, a transverse normal stress is created.

According to Eq. (5.26), this stress is not related to the normal strain ε y , so it cannot be
determined using Hooke’s law. However, an expression for the average transverse normal
stress can be obtained from the equilibrium requirement of force balance along the axis of

268 Chapter  5   Bending of Beams

 p
p b

h σx σx + dσx σy
2h z x
x y
h ∂τxy
dx y τxy τxy + dx
L ∂x
y dx dx
y
(a) (b) (c)

Figure 5.14.  ( a) Uniformly loaded cantilever beam of rectangular cross section; (b) free-body
diagram of a segment; (c) stresses in a beam element.

the beam. For this purpose, a procedure is used similar to that employed for determining
the shear stress in Section 5.7.
Consider, for example, a rectangular cantilever beam of width b and depth 2h subject
to a uniform load of intensity p (Fig. 5.14a). The free-body diagram of an isolated beam
segment of length dx is shown in Fig. 5.14b. Passing a horizontal plane through this seg-
ment results in the free-body diagram of Fig. 5.14c, for which the condition of statics
ΣFy = 0 yields
b/2 h ∂τ xy h ∂τ xy
dx dy (a)
σ y ⋅ bdx = ∫− b / 2 ∫ y ∂x
dx ⋅ dy dz = b ∫
y ∂x
Here, the shear stress is defined by Eq. (5.41) as

V 2 3V   y  2 
τ xy = (h − y2 ) = 1 − 
2I 4bh   h   (b)
After substituting Eqs. (5.28) and (b) into Eq. (a), we have
h 3p   y  2 
σy =∫ − 1 −  dy
y 4bh   h  
Integration yields the transverse normal stress in the form
p 1 3 y 1  y3
σy =−  − + 
b  2 4  h  4  h   (5.44a)
This stress varies as a cubic parabola, ranging from − plb at the surface ( y = − h) where
the load acts, to zero at the opposite surface ( y = h).
The distribution of the bending and the shear stresses in a uniformly loaded cantilever
beam (Fig. 5.12a) is determined from Eqs. (5.38) and (b):
My 3p
σx = − =− ( L − x )2 y
I 4bh3
(5.44b)
3p  y 2
τ xy = ( L − x ) 1 −   
4bh   h 

5.8   Effect of Transverse Normal Stress 269

 The largest values of σ x ,τ xy , and σ y given by Eqs. (5.44) are

3pL2 3pL p
σ x , max = ± , τ max = , σ y , max = − (c)
4bh2 4bh b
To compare the magnitudes of the maximum stresses, consider the following ratios:

τ max h σ y , max 4  h  2
= , = (d, e)
σ x , max L σ x , max 3  L 

Because L is much greater than h in most beams ( L ≥ 20h), the shear and the transverse
normal stresses will usually be orders of magnitude smaller than the bending stresses.
This justification is the rationale for assuming γ xy = 0 and ε y = 0 in the technical theory
of bending. Note that Eq. (e) results in even smaller values than Eq. (d). Therefore, in
practice, it is reasonable to neglect σ y.
The foregoing conclusion applies, in most cases, to beams of a variety of cross-­
sectional shapes and under various load configurations. Clearly, the factor of proportion-
ality in Eqs. (d) and (e) will differ for beams of different sectional forms and for different
loadings of a given beam.

5.9  COMPOSITE BEAMS

Beams constructed of two or more materials having different moduli of elasticity are
referred to as composite beams. Examples include multilayer beams made by bonding
together multiple sheets, sandwich beams consisting of high-strength material faces sepa-
rated by a relatively thick layer of low-strength material such as plastic foam, and rein-
forced concrete beams. The assumptions of the technical theory for a homogeneous beam
(Section 5.6) are valid for a beam composed of more than one material.

5.9.1  Transformed Section Method

To analyze composite beams, we will use the common transformed-section method. In
this technique, the cross sections of several materials are transformed into an equivalent
cross section of one material on which the resisting forces and the neutral axis are the
same as on the original section. The usual flexure formula is then applied to the new sec-
tion. To illustrate this method, we will use a frequently encountered example: a beam with
a symmetrical cross section built of two different materials (Fig. 5.15a).
The cross sections of the beam remain plane during bending. Hence, the condition
of geometric compatibility of deformation is satisfied. It follows that the normal strain
ε x varies linearly with the distance y from the neutral axis of the section; that is, ε x = ky
(Figs. 5.15a and b). The location of the neutral axis is yet to be determined. Both materials
composing the beam are assumed to obey Hooke’s law, and their moduli of elasticity are
designated as E1 and E2 . Then, the stress–strain relation gives

σ x1 = E1ε x = E1ky, σ x 2 = E2 ε 2 = E2 ky (5.45a, b)

270 Chapter  5   Bending of Beams

 E1, A1 E1, A1
Z
y1 1 1
εx σx1
Y y y
y2
N.A. C N.A.
z x x z
2 σx2 2
E2, A2
M

y y E1, nA2 y
y
b
nb
(a) (b) (c) (d)

Figure 5.15.  B
 eam composed of two materials: (a) composite cross section; (b) strain dis-
tribution; (d) transformed cross section.

This result is sketched in Fig. 5.13c for the assumption that E2 > E1 . We introduce the
notation
E
n = 2 (5.46)
E1

where n is called the modular ratio. Note that n > 1 in Eq. (5.46). However, this choice is
arbitrary; the technique applies as well for n < 1.

5.9.2  Equation of Neutral Axis

Referring to the cross section (Figs. 5.15a and c), the equilibrium equations Σ Fx = 0 and
Σ M z = 0 lead to

∫Aσ x dA = ∫A σ x1 dA + ∫A σ x 2 dA = 0 (a)
1 2

− ∫ yσ x dA = − ∫ yσ x1 dA − yσ x 2 dA = M (b)
A A 1

where A1 and A2 denote the cross-sectional areas for materials 1 and 2, respectively.
Substituting σ x1 , σ x 2 and n, as given by Eqs. (5.45) and (5.46), into Eq. (a) results in

∫A y dA + n ∫A y dA = 0 (5.47)
1 2

Using the top of the section as a reference (Fig. 5.15a), from Eq. (5.47) with

∫A (Y − y ) dA + n∫A (Y − y ) dA = 0
1 2

or, setting

∫A Y dA = y1 A1
1
and ∫A Y dA = y2 A2
2

we have
A1 y1 − A1 y + nA2 y2 − nA2 y = 0

5.9   Composite Beams 271

 This expression yields an alternative form of Eq. (5.47):
A1 y1 + nA2 y2
y= (5.47′)
A1 + nA2
Equations (5.47) and (5.47′) can be used to locate the neutral axis for a beam of two
materials. These equations show that the transformed section will have the same neutral
axis as the original beam, provided the width of area 2 is changed by a factor n and area 1
remains the same (Fig. 5.15d). Clearly, this widening must be effected in a direction par-
allel to the neutral axis, since the distance y2 to the centroid of area 2 remains unchanged.
The new section constructed in this way represents the cross section of a beam made of a
homogeneous material with a modulus of elasticity E1 and with a neutral axis that passes
through its centroid, as shown in Fig. 5.15d.

5.9.3  Stresses in the Transformed Beam

Similarly, condition (b) together with Eqs. (5.45) and (5.46) leads to

M = − kE1 (∫ A1
2
y dA + n ∫A y
2
2
)
dA
or
M = − kE1 ( I1 + nI 2 ) = − kE1 I t (5.48)
where I1 and I 2 are the moments of inertia about the neutral axis of the cross-sectional
areas 1 and 2, respectively. Note that
I t = I1 + nI 2 (5.49)
is the moment of inertia of the entire transformed area about the neutral axis. From Eq.
(5.48), we have
M
k=−
E1 I t
The flexure formulas for a composite beam are obtained by introducing this relation
into Eqs. (5.45):
My nMy
σ x1 = − , σ x2 = − (5.50)
It It
where σ x1 and σ x2 are the stresses in materials 1 and 2, respectively. Note that when
E1 = E2 = E, Eqs. (5.50) reduce to the flexure formula for a beam of homogeneous mate-
rial, as expected.

5.9.4  Composite Beams of Multi Materials

The preceding discussion may be extended to include composite beams consisting of
more than two materials. It is readily shown that for m different materials, Eqs. (5.47′),
(5.49), and (5.50) take the forms
A1 y1 + Σ ni Ai yi Ei
y= , ni = (5.51)
A1 + Σ ni Ai E1

272 Chapter  5   Bending of Beams

 I t = I1 + Σni I i (5.52)

My ni My
σ x1 = − , σ xi = − (5.53)
It It

where i = 2, 3, . . . , m denotes the ith material.

The use of the formulas developed in this section is demonstrated in the solutions of
the numerical problems in Examples 5.7 and 5.8.

EXAMPLE 5.7  Aluminum-Reinforced Wood Beam

A wood beam with Ew = 8.75 GPa, 100 mm wide by 220 mm deep, has an
aluminum plate Ea = 70 GPa with a net section 80 mm by 20 mm securely
fastened to its bottom face, as shown in Fig. 5.16a. Dimensions are given in
millimeters. The beam is subjected to a bending moment of 20 kN ⋅ m around
a horizontal axis. Calculate the maximum stresses in both materials (a)
using a transformed section of wood and (b) using a transformed section of
aluminum.

Solution
a. The modular ratio is n = Ea / Ew = 8. The centroid and the moment of inertia
about the neutral axis of the transformed section (Fig. 5.16b) are

100(220)(110) + 20(640)(230)
y= = 154.1 mm
100(220) + 20(640)
I t = 121 (100)(220)3 + 100(220)(44.1)2 + 121 (640)(20)3 + 640(20)(75.9)2
= 205.6 × 106 mm 4

12.5
100 100
Z
Y 154.1
y
220 N.A. C N.A. N.A.
C z
z z 85.9

y y y
20 80 × 8 = 640 80
80
(a) (b) (c)

Figure 5.16.  E
 xample 5.7. (a) Composite cross section; (b) equivalent wood cross
­section; (c) ­equivalent aluminum cross section.

5.9   Composite Beams 273

 In turn, the maximum stresses in the wood and aluminum portions are

Mc 20 × 103 (0.1541)
σ w , max = = = 15 MPa
It 2056(10−6 )

nMc 8(20 × 103 )(0.0859)
σ a, max = = = 66.8 MPa
It 205.6(10−6 )

At the juncture of the two parts,

Mc 20 × 103 (0.0659)
σ w , min = = = 6.4 MPa
It 2056(10−6 )
σ a, min = n(σ w , min ) = 8(6.4) = 51.2 MPa

b. For this case, the modular ratio is n = Ew / Ea = 1/8 and the transformed area
is shown in Fig. 5.16c. We now have

I t = 121 (12.5)(220)3 + 12.5(220)(44.1)2 + 121 (80)(20)3 + 80(20)(75.9)2

= 25.7 × 106 mm 4
Then

Mc 20 × 103 (0.0859)
σ a, max = = = 66.8 MPa
It 25.7(10−6 )

nMc 20 × 103 (0.1541)
σ w , max = = = 15 MPa
It 8(25.7 × 10−6 )

as have already been found in part (a).

EXAMPLE 5.8  Steel-Reinforced-Concrete Beam

A concrete beam of width b = 250 mm and effective depth d = 400 mm
is reinforced with three steel bars, providing a total cross-sectional area
As = 1000 mm 2 (Fig. 5.17a). Dimensions are given in millimeters. Note that
it is usual for an approximate allowance a = 50 mm to be used to protect
the steel from corrosion and fire. Let n = Es / Ec = 10. Calculate the maxi-
mum stresses in the materials produced by a negative bending moment of
M = 60 kN ⋅ m.
Solution  Concrete is very weak in tension but strong in compression. Thus,
only the portion of the cross section located a distance kd above the neutral
axis is used in the transformed section (Fig. 5.17b); the concrete is assumed to
take no tension. The transformed area of the steel nAs is identified by a single
dimension from the neutral axis to its centroid. The compressive stress in the
concrete is assumed to vary linearly from the neutral axis, and the steel is
assumed to be uniformly stressed.

274 Chapter  5   Bending of Beams

 σc, max
b = 250
b
M kd
kd = 143.3 3
z N.A. x C
d = 400 z d (1 − k
)
d (1 − k) = 256.7 3

a = 50 T
y σs
y As y 3
nAs = 10(103) mm2
(a) (b) (c)

Figure 5.17.  E
 xample 5.8. (a) Reinforced-concrete cross section; (b) equivalent con-
crete cross section; (c) compressive force C in concrete and tensile force
T in the steel rods.

The condition that the first moment of the transformed section with
respect to the neutral axis be zero is satisfied by
kd
b( kd ) − nAs ( d − kd ) = 0
2
or
2n 2n
( kd )2 + ( kd ) As − dAs = 0 (5.54)
b b
By solving this quadratic expression for kd, we can obtain the position of the
neutral axis.
Introducing the data given, Eq. (5.54) reduces to

( kd )2 + 80( kd ) − 32 × 103 = 0
from which
kd = 143.3 mm and hence k = 0.358 (c)
The moment of inertia of the transformed cross section about the neutral axis is

I t = 121 (0.25) + (0.1433)3 + 0.25(0.1433)(0.0717)2 + 0 + 10 × 10−3 (0.2567)2

= 904.4 × 10−6 m 4

Thus, the peak compressive stress in the concrete and the tensile stress in the
steel are

Mc 60 × 103 (0.1433)
σ c , max = = = 9.5 MPa
It 904.4 × 10−6

nMc 10(60 × 103 )(0.2567)
σs = = = 170 MPa
It 904.4 × 10−6

5.9   Composite Beams 275

 These stresses act as shown in Fig. 5.15c.
An alternative method of solution is to obtain σ c, max and σ s from a free-
body diagram of the portion of the beam (Fig. 5.17c) without computing I t .
The first equilibrium condition, ΣFx = 0, gives C = T , where

C = 12 σ c , max (b ⋅ kd ), T = σ s ( As ) (d)

are the compressive and tensile stress resultants, respectively. From the sec-
ond requirement of statics, ΣM z = 0, we have

M = Cd (1 − kl 3) = Td (1 − kl 3) (e)

Equations (d) and (e) result in

2M M
σ c , max = 2 , σs = (5.55)
bd k (1 − k /3) As d (1 − k /3)

Substituting the data given and Eq. (c) into Eq. (5.55) yields

2(60 × 103 )
σ c , max = = 9.5 MPa
0.25(0.4 2 )(0.358)(1 − 0.358/3)

60 × 103
σs = = 170 MPa
1000 × 10−6 (0.4)(1 − 0.358/3)

as before.

5.10  SHEAR CENTER

Given any cross-sectional configuration, one point may be found in the plane of the
cross section through which the resultant of the transverse shearing stresses passes.
A  ­transverse  load applied on the beam must act through this point, called the shear
­center or flexural center, if no twisting is to occur [Ref. 5.3]. The center of shear
is sometimes defined as the point in the end section of a cantilever beam at which an
applied load results in bending only. When the load does not act through the shear ­center,
in addition to bending, a twisting action results (Section 6.1).
The location of the shear center is independent of the direction and magnitude of
the transverse forces. For singly symmetrical sections, the shear center lies on the axis of
symmetry, while for a beam with two axes of symmetry, the shear center coincides with
their point of intersection (also the centroid). It is not necessary, in general, for the shear
center to lie on a principal axis, and it may be located outside the cross section of the
beam.

276 Chapter  5   Bending of Beams

 τ
Vy
C
C A A, S z
z e τxz
S e
Vz

τxy
y y
(a) (b) (c)

Figure 5.18.  Shear centers.

5.10.1  Thin-Walled Open Cross Sections

For thin-walled sections, the shearing stresses are taken to be distributed uniformly over
the thickness of the wall and directed so as to parallel the boundary of the cross section.
If the shear center S for the typical section of Fig. 5.18a is required, we begin by calcu-
lating the shear stresses by means of Eq. (5.39). The moment M x of these stresses about
arbitrary point A is then obtained. Since the external moment attributable to V y about A is
V y e , the distance between A and the shear center is given by
Mx
e= (5.56)
Vy
If the force is parallel to the z axis rather than the y axis, the position of the line of
action may be established in the manner discussed previously. If both V y and Vz exist, the
shear center is located at the intersection of the two lines of action.
The determination of M x is simplified by propitious selection of point A, such as in
Fig. 5.18b. There, the moment M x of the shear forces about A is zero; point A is also the
shear center. For all sections consisting of two intersecting rectangular elements, the same
situation exists.
For thin-walled box beams (with boxlike cross section), the point or points in the wall
where the shear flow q = 0 (or τ xy = 0) is unknown. Here, shear flow is represented by the
superposition of transverse and torsional flow (see Section 6.8). Hence, the unit angle of
twist equation, Eq. (6.23), along with q = VQ/ I , is required to find the shear flow for a
cross section of a box beam. The analysis procedure is as follows: First, introduce a free
edge by cutting the section open; second, close it again by obtaining the shear flow that
makes the angle of twist in the beam zero [Refs. 5.4 through 5.6].
In the following examples, the shear center of an open, thin-walled section is deter-
mined for two typical situations. In the first, the section has only one axis of symmetry; in
the second, there is an asymmetrical section.

EXAMPLE 5.9  Shearing Stress Distribution in a Channel Section

Locate the shear center of the channel section loaded as a cantilever
(Fig.  5.19a). Assume that the flange thicknesses are small when compared
with the depth and width of the section.

5.10   Shear Center 277

 N + dN

dx
τ1
τ t1 • dx
τ2
P
F1
N C
(b) z
S
P F2 = P
s t1 A
n F3 = F1
m e
n
C y
t2 (c)
z 2h
(d)

(a) y
t1
b

Figure 5.19.  E
 xample 5.9. (a) Cantilever beam with a concentrated load at the
free end; (b) an element of the upper flange; (c) shear distribution;
(d) location of the shear center S.

Solution  The shearing stress in the upper flange at any section nn will be
found first. This section is located a distance s from the free edge m, as shown
in the figure. At m, the shearing stress is zero. The first moment of area st1
about the z axis is Qz = st1h. The shear stress at nn, from Eq. (5.39), is thus
V y Qz sh
τ xz = =P (a)
I zb Iz

The direction of τ along the flange can be determined from the equilibrium
of the forces acting on an element of length dx and width s (Fig. 5.19b). Here
the normal force N = t1sσ x , owing to the bending of the beam, increases with
dx by dN. Hence, the x equilibrium of the element requires that τ t1 ⋅ dx must
be directed as shown in the figure. This flange force is directed to the left,
because the shear forces must intersect at the corner of the element.
The distribution of the shear stress τ xz on the flange, as Eq. (a) indicates,
is linear with s. Its maximum value occurs at
bh
τ1 = P (b)
Iz

Similarly, the value of stress τ xy at the top of the web is

bt1h
τ2 = P (c)
t2 I 2
The stress varies parabolically over the web, and its maximum value is found
at the neutral axis. Figure 5.19c sketches the shear stress distribution in the

278 Chapter  5   Bending of Beams

 channel. As the shear stress is linearly distributed across the flange length,
from Eq. (b), the flange force is expressed by

b 2 ht1
F1 = 12 τ 1bt1 = P (d)
2I z

Symmetry of the section dictates that F1 = F3 (Fig. 5.19d). We will

assume that the web force F2 = P, since the vertical shearing force transmitted
by the flange is negligibly small, as shown in Example 5.3. The shearing force
components acting in the section must be statically equivalent to the resultant
shear load P. Thus, the principle of the moments for the system of forces in
Fig. 5.19d or Eq. (5.56), applied at A, yields M x = Pe = 2F1h. Substituting F1
from Eq. (d) into this expression, we obtain

b 2 h2 t1
e=
Iz

Since for the usual channel section t1 is small in comparison to b or h, the sim-
plified moment of inertia has the following form:
I z = 23 t2 h3 + 2bt1h2

The shear center is, in turn, located by the expression

b 2 t1
e= 3
2 (e)
ht2 + 3bt1
Comments  Note that e depends on only the section dimensions. Exami-
nation reveals that e may vary from a minimum of zero to a maximum of
b/2. A zero or near-zero value of e corresponds to either a flangeless beam
(b = 0, e = 0) or an especially deep beam ( h  b ). The extreme case, e = b/2,
is obtained for an infinitely wide beam.

EXAMPLE 5.10  Shear Flow in an Asymmetrical Channel Section

Locate the shear center S for the asymmetrical channel section shown in
Fig. 5.20a. All dimensions are in millimeters. Assume that the beam thickness
t = 125 mm is constant.
Solution  The centroid C of the section is located by y and z with respect to
nonprincipal axes z and y. By performing the procedure given in Example 5.1, we
obtain y = 15.63 mm, z = 5.21 mm,  I y = 4765.62 mm 4 , I z = 21,054.69 mm 4,
and I yz = 3984.37 mm 4 . Equation (5.18) then yields the direction of the prin-
cipal axis x ′, y ′ as θ p = 13.05°, and Eq. (5.19) gives the principal moments of
inertia as I y ′ = 3828.12 mm 4 , I z ′ = 21,953.12 mm 4 (Fig. 5.20a).
Let us now assume that a shear load V y ′ is applied in the y ′, z ′ plane (Fig.
5.20b). This force may be considered the resultant of force components F1, F2 ,

5.10   Shear Center 279

 12.5 s s

1.25 F1 F1
19.55 12.05
F2 z′ F2 z′
z′ C
37.5 z z Vy ′ Vy ′
Vz ′ S
C 13.05 A A F3 F4
y ey ′
F3
ez ′ ez ′
25 y′ y′
y′
y

(a) (b) (c)

Figure 5.20.  E
 xample 5.10. (a) Portion of a beam with a channel cross section;
(b) shear flow; (c) location of the shear center S.

and F3 acting in the flanges and web in the directions indicated in the figure.
The algebra will be minimized if we choose point A, where F2 and F3 inter-
sect, in finding the line of action of V y ′ by applying the principle of moments.
To do so, we need to determine the value of F1 acting in the upper flange. The
shear stress τ xz in this flange, from Eq. (5.39), is

V y ′ Qz ′ Vy′
τ xz = =  st (19.55 + 12 s sin13.05° )  (f )
I z′b I z′t 

where s is measured from right to left along the flange. Note that Qz ′ the
bracketed expression, is the first moment of the shaded flange element area
with respect to the z ′ axis. The constant 19.55 is obtained from the geometry
of the section. After substituting the numerical values and integrating Eq. (f),
the total shear force in the upper flange is found to be

s V y ′t 12.5
F1 = ∫ τ xz t ds =
I z ′ ∫0
s(19.55 + 12 s sin13.05°) ds = 0.0912V y ′ (g)
0

Application of the principle of moments at A gives V y ′ ez ′ = 37.5F1 .

Introducing F1 from Eq. (g) into this equation, the distance ez ′ , which locates
the line of action of V y ′ from A, is

ez ′ = 3.42 mm (h)

Next, assume that the shear loading Vz ′ acts on the beam (Fig. 5.20c). The
distance e y ′ may be obtained as in the situation just described. Because of Vz ′

280 Chapter  5   Bending of Beams

 the force components F1 to F4 will be produced in the section. The shear stress
in the upper flange is given by

Vz ′ Q y ′ Vz ′
τ xz = = [st (12.05 − 12 s cos13.05°)] (i)
I y′b I y′t

The quantity Q y ′ represents the first moment of the flange segment area with
respect to the y ′ axis, and 12.05 is found from the geometry of the section.
The total force F1 in the flange is

Vz ′ 12.5
F1 =
I y′ ∫0 st (12.05 − 12 s cos13.05°) ds = 0.204Vz ′

The principle of moments applied at A, Vz ′ e y ′ = 37.5F1 = 7.65Vz ′ leads to

e y ′ = 7.65 mm ( j)

Thus, the intersection of the lines of action of V y ′ and Vz ′ , and ez ′ and e y ′ ,

locates the shear center S of the asymmetrical channel section.

5.10.2  Arbitrary Solid Cross Sections

The preceding considerations can be extended to beams of arbitrary solid cross section,
in which the shearing stress varies with both cross-sectional coordinates y and z. For these
sections, the exact theory can, in some cases, be successfully applied to locate the shear
center. Examine the section of Fig. 5.18c subjected to the shear force Vz , which produces
the stresses indicated. Denote y and z as the principal directions. The moment about the x
axis is then

M x = ∫∫ (τ xy z − τ xz y ) dz dy (5.57)

Vz must be located a distance e from the z axis, where e = M x /Vz .

5.11  STATICALLY INDETERMINATE SYSTEMS

A large class of problems of considerable practical interest relate to structural systems for
which the equations of statics are not sufficient (but are necessary) for determination of
the reactions or other unknown forces. Such systems are called statically indeterminate,
and they require supplementary information for their solution. Additional equations usu-
ally describe certain geometric conditions associated with displacement or strain. These
equations of compatibility state that the strain owing to deflection or rotation must pre-
serve continuity.
With this additional information, the solution proceeds in essentially the same man-
ner as for statically determinate systems. The number of reactions in excess of the number

5.11   Statically Indeterminate Systems 281

 Figure 5.21.  S uperposition of displacements P
in a continuous beam. (a)
A
R1 R R2

P
(b)

(vA)P A
A
(c)
R
(vA)R

of equilibrium equations is called the degree of statical indeterminacy. Any reaction in

excess of that which can be obtained by statics alone is said to be redundant. Thus, the
number of redundants is the same as the degree of indeterminacy.
Several methods are available to analyze statically indeterminate structures. The prin-
ciple of superposition, briefly discussed next, is an effective approach for many cases. In
Section 5.6 and in Chapters 7 and 10, a number of commonly employed methods are dis-
cussed for the solution of the indeterminate beam, frame, and truss problems.

5.11.1  The Method of Superposition

In the event of complicated load configurations, the method of superposition may be used
to good advantage to simplify the analysis. Consider, for example, the continuous beam
shown in Fig. 5.21a, which is then replaced by the beams shown in Fig. 5.21b and c.
At point A, the beam now experiences the deflections (υ A ) P and (υ A ) R due to P and R,
respectively. Subject to the restrictions imposed by small deformation theory and a mate-
rial obeying Hooke’s law, the deflections and stresses are linear functions of transverse
loadings, and superposition is valid:

υ A = (υ A ) P + (υ A ) R

σ A = (σ A ) P + (σ A ) R

This procedure may, in principle, be extended to situations involving any degree of

indeterminacy.

EXAMPLE 5.11  Displacements of a Propped Cantilever Beam

Figure 5.22 shows a propped cantilever beam AB subject to a uniform load
of intensity p. Determine (a) the reactions, (b) the equation of the deflection
curve, and (c) the slope at A.
Solution Reactions RA , RB and M B are statically indeterminate because
there are only two equilibrium conditions ( Σ Fy = 0, Σ M z = 0); thus, the beam

282 Chapter  5   Bending of Beams

 p Figure 5.22.  E
 xample 5.11. A propped
MB beam under uniform load.

A x
B

RA L RB

is statically indeterminate to the first degree. With the origin of coordinates

taken at the left support, the equation for the beam moment is
M = − R x + 1 px 2
A 2
The third of Eqs. (5.32) then becomes
EIυ ′′ = − RA x + 12 px 2
and successive integrations yield
EIυ ′ = − 12 RA x 2 + 61 px 3 + c1
(a)
EIυ = − 61 RA x 3 + 241 px 4 + c1 x + c2
There are three unknown quantities in these equations ( c1, c2 , and RA) and
three boundary conditions:

υ (0) = 0, υ ′( L) = 0, υ ( L) = 0 (b)
a. Introducing Eqs. (b) into the preceding expressions, we obtain c2 = 0,, c1 = pL3 /48
c2 = 0, c1 = pL3 /48 and
RA = 83 pL (5.58a)
We can now determine the remaining reactions from the equations of
equilibrium:
RB = 85 pL, M B = 81 pL2 (5.58b, c)

b. By substituting for RA , c1 and c2 in Eq. (a), we obtain the equation of the

deflection curve:
p
υ= (2x 4 − 3Lx 3 + L3 x ) (5.59)
48EI
c. Differentiating Eq. (5.59) with respect to x, we obtain the equation of the
angle of rotation:
p
θ= (8x 3 − 9Lx 2 + L3 ) (5.60)
48EI
Setting x = 0, we have the slope at A:
pL3
θA = (5.61)
48EI

5.11   Statically Indeterminate Systems 283

 (vA)R

p p
A A B
= +
B (vA)p B RA
RA A
(a) (b) (c)

 xample 5.12. Method of superposition: (a) reaction RA is selected as redundant;

Figure 5.23.  E
(b) deflection at end A due to load P; (c) deflection at end A due to reaction RA.

EXAMPLE 5.12  Reactions of a Propped Cantilever Beam

Consider again the statically indeterminate beam shown in Fig. 5.22. Determine
the reactions using the method of superposition.
Solution Reaction RA is selected as redundant; it becomes an unknown
load when we eliminate the support at A (Fig. 5.23a). The loading is resolved
into those shown in Figs. 5.23b and 5.23c. The solution for each case is (see
Table D.4)

pL4 RA L3
(υ A ) P = , (υ A ) R = −
8EI 3EI
The compatibility condition for the original beam requires that
pL4 RA L3
υA = − =0
8EI 3EI
from which RA = 3pL /8. Reaction RB and moment M B can now be found
from the equilibrium requirements. The results correspond to those of
Example 5.11.

5.12  ENERGY METHOD FOR DEFLECTIONS

Strain energy methods are frequently employed to analyze the deflections of beams and
other structural elements. Of the many approaches available, Castigliano’s second theo-
rem is one of the most widely used. To apply this theory, the strain energy must be rep-
resented as a function of loading. Detailed discussions of energy techniques are found
in Chapter 10. In this section, we limit ourselves to a simple example to illustrate how
the strain energy in a beam is evaluated and how the deflection is obtained by the use of
Castigliano’s theorem (Section 10.4).
The strain energy stored in a beam under bending stress σ x only, substituting
M = EI ( d 2υ/dx 2 ) into Eq. (2.63), is expressed in the form
2
M 2 dx EI  d 2υ 
Ub = ∫ =∫ dx (5.62)
2EI 2  dx 2 

284 Chapter  5   Bending of Beams

where the integrations are carried out over the beam length. We next determine the strain
energy stored in a beam that is only due to the shear loading V. As described in Section
5.7, this force produces shear stress τ xy at every point in the beam. The strain energy den-
sity is, from Eq. (2.50), U o = τ xy /2G . Substituting τ xy as expressed by Eq. (5.39), we have
U o = V 2Q 2 /2GI 2b 2. Integrating this expression over the volume of the beam of cross-­
sectional area A, we obtain

V 2  Q2 
Us = ∫
2GI 2  ∫ b 2
dA dx (a)


5.12.1  Form Factor for Shear

Let us denote

A Q2
I 2 ∫ b2
α= dA (5.63)

This value is termed the shape or form factor for shear. When it is substituted in Eq. (a),
we obtain

αV 2 dx
Us = ∫ (5.64)
2 AG

where the integration is carried over the beam length. The form factor is a dimensionless
quantity specific to a given cross-section geometry.
For example, for a rectangular cross section of width b and height 2h, the first
moment Q, from Eq. (5.41), is Q = (b/2)( h2 − y12 ). Because A/ I 2 = 9/2bh5 , Eq. (5.63) pro-
vides the following result:

9 h
α= ∫− h 14 ( h
2
− y12 )2 b dy1 = 65 (b)
2bh5

In a like manner, the form factor for other cross sections can be determined. Table 5.2 lists
several typical cases. Following the determination of α , the strain energy is evaluated by
applying Eq. (5.64).

Table 5.2.  Form Factor for Shear for Various Beam Cross Sections

Cross Section Form Factor α

A. Rectangle 6/5
B. I-section, box section, or channelsa A/Aweb
C. Circle 10/9
D. Thin-walled circular 2
a
A = area of the entire section, Aweb = area of the web ht, where h is the beam depth and t is the web
thickness.

5.12   Energy Method for Deflections 285

 For a linearly elastic beam, Castigliano’s theorem, from Eq. (10.3), is expressed by

∂U
δ= (c)
∂P
where P is a load acting on the beam and δ is the displacement of the point of application
in the direction of P. Note that the strain energy U = U b + U s is expressed as a function of
the externally applied forces (or moments).
As an illustration, consider the bending of a cantilever beam of rectangular cross
­section and length L, subjected to a concentrated force P at the free end (Fig. 5.5). The
bending moment at any section is M = Px and the shear force V is equal in magnitude to
P. Substituting these together with α = 65 into Eqs. (5.62) and (5.64) and integrating, we
find the strain energy stored in the cantilever to be

P 2 L3 3P 2 L
U= +
6EI 5AG

The displacement of the free end owing to bending and shear is, by application of
Castigliano’s theorem, therefore

PL3 6PL
δ =υ = +
3EI 5AG

The exact solution is given by Eq. (5.24).

Part C: Curved Beams

5.13  ELASTICITY THEORY

Our treatment of stresses and deflections caused by the bending has been restricted so far
to straight members. In real-world applications, many members—such as crane hooks,
chain links, C-lamps, and punch-press frames—are curved and loaded as beams. Part C
deals with the stresses caused by the bending of bars that are initially curved.
A curved bar or beam is a structural element for which the locus of the centroids of
the cross sections is a curved line. This section focuses on an application of the theory of
elasticity to a bar characterized by a constant narrow rectangular cross section and a cir-
cular axis. The axis of symmetry of the cross section lies in a single plane throughout the
length of the member.

5.13.1  Equations of Equilibrium and Compatibility

Consider a beam subjected to equal end couples M such that bending takes place in the
plane of curvature, as shown in Fig. 5.24a. Given that the bending moment remains con-
stant along the length of the bar, the stress distribution should be identical in any radial

286 Chapter  5   Bending of Beams

 n n
+
σθ σr
h
σθ
m m
M t
r a
b (b)
dθ
M
(a)

Figure 5.24.  Pure bending of a curved beam of rectangular cross section.

cross section. Stated differently, we seek a distribution of stress displaying θ indepen-

dence. It is clear that the appropriate expression of equilibrium is Eq. (8.2):

dσ r σ r − σ θ
+ = 0 (a)
dr r
and that the condition of compatibility for plane stress, Eq. (3.41),

d 2 (σ r + σ θ ) 1 d (σ r + σ θ )
+ = 0
dr 2 r dr
must also be satisfied. The latter is an equidimensional equation, which can be reduced
to a second-order equation with constant coefficients by substituting r = e t or t = ln r.
Direct integration then leads to σ r + σ θ = c ′′ + c ′ ln r , which may be written in the form
σ r + σ θ = c ′′′ + c ′ In ( r/a ). Solving this expression together with Eq. (a) results in the fol-
lowing equations for the radial and tangential stress:
r c
σ r = c1 + c2 1n + 23
a r
(5.65)
 r  c3
σ θ = c1 + c2  1 + 1n  − 2
 a r

5.13.2  Boundary Conditions

To evaluate the constants of integration, the boundary conditions are applied as follows:
1. No normal forces act along the curved boundaries at r = a and r = b , so that
(σ r )r = a = (σ r )r = b = 0 (b)

2. Because there is no force acting at the ends, the normal stresses acting at the straight
edges of the bar must be distributed to yield a zero resultant:
b
t ∫ σ θ dr = 0 (c)
a

where t represents the beam thickness.

5.13   Elasticity Theory 287

 3. The normal stresses at the ends must produce a couple M:
b
t ∫ rσ θ dr = M (d)
a

The conditions (c) and (d) apply not just at the ends; that is, because of σ θ indepen-
dence, they apply at any θ . In addition, shearing stresses are assumed to be zero through-
out the beam, so τ rθ = 0 is satisfied at the boundaries, where no tangential forces exist.
Combining the first equation of Eqs. (5.65) with the condition (b), we find that

 a2  b
c3 = − a 2 c1, c1  2 − 1 = c2 ln
b  a

These constants together with the second of Eqs. (5.65) satisfy condition (c). Thus, we have

b 2 ln(b /a ) a 2b 2 ln(b /a )
c1 = c2 , c3 = c2 (e)
a2 − b2 b2 − a2
Finally, substitution of the second of Eqs. (5.65) and (e) into (d) provides
M 4(b 2 − a 2 )
c2 = (f)
N tb 4
where
2
 a2  a2 b
N =  1 − 2  − 4 2 ln 2 (5.66)
 b  b a

5.13.3  Stress Distribution

When the expressions for constants c1, c2 , and c3 are inserted into Eq. (5.65), the follow-
ing equations are obtained for the radial stress and tangential stress:

4M  a2  r  a2  b 
σr =  1 − 2  ln −  1 − 2  ln 
tb 2 N  b  a  r  a
(5.67)
4M  a2   r  a2  b 
σθ = 2  1 −  1 + ln  − 1 + ln 
tb N  b 2   a   r 2  a 

If the end moments are applied so that the force couples producing them are distrib-
uted in the manner indicated by Eq. (5.67), then these equations are applicable throughout
the bar. If the distribution of applied stress (to produce M) differs from Eq. (5.67), the
results may be regarded as valid in regions away from the ends, in accordance with Saint-
Venant’s principle.
These results, when applied to a beam with radius a, large relative to its depth h, yield
an interesting comparison between straight and curved beam theory. For slender beams
with h << a, the radial stress σ r in Eq. (5.67) becomes negligible, and the tangential stress
σ θ is approximately the same as that obtained from My/I. Note that the radial stresses
developed in nonslender curved beams made of isotropic materials are small enough that
they can be neglected in analysis and design.

288 Chapter  5   Bending of Beams

 The bending moment is taken as positive when it tends to decrease the radius of
c­ urvature of the beam, as in Fig. 5.24a. Employing this sign convention, σ r as determined
from Eq. (5.67) is always negative, indicating that this stress is compressive. Similarly,
when σ θ is positive, the stress is tensile; otherwise, it is compressive. Figure 5.24b plots
the stresses at section mn. Note that the maximum stress magnitude is found at the
extreme fiber of the concave side.

5.13.4 Deflections
Substitution of σ r and σ θ from Eq. (5.67) into Hooke’s law provides expressions for the
strains εθ , ε r , and γ rθ . The displacements u and υ then follow, upon integration, from the
strain-displacement relationships, described by Eqs. (3.33). The resulting displacements
indicate that plane sections of the curved beam subjected to pure bending remain plane
subsequent to bending. Castigliano’s theorem (Section 5.12) is a particularly attractive
method for determining the deflection of curved members.
For beams in which the depth of the member is small relative to the radius of curva-
ture or, as is usually assumed, r /c > 4, the initial curvature may be neglected in evaluating
the strain energy. In such a case, r represents the radius to the centroid, and c is the dis-
tance from the centroid to the extreme fiber on the concave side. Thus, the strain energy
due to the bending of a straight beam [Eq. (5.62)] is also a good approximation for curved,
slender beams.

5.14  CURVED BEAM FORMULA

The approach to curved beams explored in this section was developed by E. Winkler
(1835–1888). As an extension of the elementary theory of straight beams, Winkler’s the-
ory assumes that all conditions required to make the straight-beam formula applicable are
satisfied except that the beam is initially curved.
Consider the pure bending of a curved beam as illustrated in Fig. 5.25a. The distance
from the center of curvature to the centroidal axis is r . The positive y coordinate is mea-
sured toward the center of curvature O from the neutral axis (Fig. 5.25b). The outer and
inner fibers are at distances of ro and ri from the center of curvature, respectively.

5.14.1  Basic Assumptions

Derivation of the stress in the beam is again based on the three principles of solid mechan-
ics and the familiar presuppositions:
1. All cross sections possess a vertical axis of symmetry lying in the plane of the centroi-
dal axis passing through C.
2. The beam is subjected to end couples M. The bending moment vector is everywhere
normal to the plane of symmetry of the beam.
3. Sections originally plane and perpendicular to the centroidal beam axis remain so sub-
sequent to bending. (The influence of transverse shear on beam deformation is not
taken into account.)

5.14   Curved Beam Formula 289

 Centroidal
axis
b f b
a dθ
c Curved
x e z beam
θ y
R c
e c
d M Straight beam
θ + dθ N.A. y dA
(b) (c)
r O′
M
ri θ
R
ro r
(a)
O O

Figure 5.25.  ( a) Curved beam in pure bending with a cross-sectional vertical (y) axis of symmetry;
(b) cross section; (c) stress distributions over the cross section.

Referring to assumption (3), note the relationship in Fig. 5.25a between lines bc and ef
representing the plane sections before and after the bending of an initially curved beam.
Note also that the initial length of a beam fiber such as gh depends on the distance r from
the center of curvature O. On the basis of plane sections remaining plane, we can state
that the total deformation of a beam fiber obeys a linear law, as the beam element rotates
through small angle.

5.14.2  Location of the Neutral Axis

In Fig. 5.25a, it is clear that the initial length of any arbitrary fiber gh of the beam
depends on the distance r from the center of curvature O. Thus, the total deformation of
a beam fiber obeys a linear law, as the beam element rotates through a small angle dθ . 
Conversely, the normal or tangential strain εθ does not follow a linear relationship. The
contraction of fiber gh equals − ( R − r )dθ , where R is the distance from O to the neutral
axis (yet to be determined) and its initial length is rθ . So, the normal strain of this fiber is
given by εθ = −( R − r )dθ /rθ . For convenience, we denote λ = dθ /θ , which is constant for
any element.
The tangential normal stress, acting on an area dA of the cross section, can now be
obtained through the use of Hooke’s law, σ θ = Eεθ . It follows that

R−r
σ θ = − Eλ (a)
r
The equations of equilibrium, ∑ Fx = 0 and ∑ M z = 0 are, respectively,

∫ σ θ dA = 0 (b)

∫ σ θ ( R − r ) dA = M (c)

290 Chapter  5   Bending of Beams

When the tangential stress of Eq. (a) is inserted into Eq. (b), we obtain

 R − r  dA = 0
∫A − Eλ  r 
 (d)

Since Eλ and R are constants, they may be moved outside the integral sign, as follows:

dA
Eλ  R ∫ − dA = 0
 A r ∫A 

The radius of the neutral axis R is then written in the form

A
R= (5.68)
dA
∫A r
where A is the cross-sectional area of the beam. The integral in Eq. (5.68) may be evalu-
ated for various cross-sectional shapes (see Example 5.13 and Probs. 5.46 through 5.48).
For reference, Table 5.3 lists explicit formulas for R and A for some commonly used cases.
The distance e between the centroidal axis and the neutral axis (y = 0) of the cross
section of a curved beam (Fig. 5.25b) is equal to
e = r − R (5.69)
Thus, in a curved member, the neutral axis does not coincide with the centroidal axis. This
differs from the case involving straight elastic beams.

5.14.3  Tangential Stress

Once we know the location of the neutral axis, we can obtain the equation for the stress
distribution by introducing Eq. (a) into Eq. (c). Therefore,

( R − r )2
M = Eλ ∫ dA
A r
Expanding this equation, we have

dA
M = Eλ  R 2 ∫ − 2R ∫ dA + ∫ r dA
 A r A A 

Here, the first integral is equivalent to A/R, as determined by Eq. (5.68), and the second
integral equals the cross-sectional area A. The third integral, by definition, represents rA,
where r is the radius of the centroidal axis. Therefore, M = Eλ A( r − R) = Eλ Ae .
We now introduce E from Eq. (a) into the discussion and solve for σ θ from the result-
ing expression. Then, the tangential stress in a curved beam, subject to pure bending at a
distance r from the center of curvature, is expressed in the following form:
M (R − r)
σθ = (5.70)
Aer

5.14   Curved Beam Formula 291

 Table 5.3.  Properties for Various Cross-Sectional Shapes
Cross Section Radius of Neutral Surface R
A. Rectangle h
R=
r
h ri ln o
r1
b O A = bh
ro

B. Circle A
R=
2π ( r + r 2 − c 2 )
A = π c2
O

r–

C. Ellipse A
R=
2π b
2a (r − r 2 − a2 )
a
A = π ab
2b O

r–

D. Triangle A
R=
h ri bro  ro 
ln −b
h  ri 
O
b A = 12 bh
ro

E. Trapezoid A
R=
r
h ri 1
h [(b1ro − b2 ri ) ln o − h(b1 − b2 )]
ri
A = 12 (b1 + b2 ) h
b2 O
b1

ro

292 Chapter  5   Bending of Beams

where e is defined by Eq. (5.69). Alternatively, substituting y = R − r or r = R − y
(Fig. 5.25a) into Eq. (5.70) yields
My
σθ = − (5.71)
Ae ( R − y )

5.14.4  Winkler’s Formula

Equations (5.70) and (5.71) represent two forms of the curved-beam formula. Another
alternative form of these equations is often referred to as Winkler’s formula. The variation
of stress over the cross section is hyperbolic, as sketched in Fig. 5.25c. The sign convention
applied to bending moment is the same as that used in Section 5.13—namely, the bending
moment is positive when directed toward the concave side of the beam, as shown in the fig-
ure. If Eq. (5.70) or Eq. (5.71) results in a positive value, a tensile stress is present.

5.15  COMPARISON OF THE RESULTS OF VARIOUS THEORIES

We now examine the solutions obtained in Sections 5.13 and 5.14 with results determined
using the flexure formula for straight beams. To do so, we consider a curved beam of
rectangular cross section and unit thickness experiencing pure bending. The tangential
stress predicted by the elementary theory (based on a linear distribution of stress) is My/I.
The Winkler approach, which leads to a hyperbolic distribution, is given by Eq. (5.70) or
Eq. (5.71), while the exact theory results in Eqs. (5.67). In each case, the maximum and
minimum values of stress are expressed by

M
σθ = m (5.72)
a2

Table 5.4 lists values of m as a function of b/a for the four cases cited [Ref. 5.1], in
which b = ro and a = rj (see Figs. 5.24 and 5.25). As can be seen, there is good agree-
ment between the exact and Winkler results. On this basis, as well as from more exten-
sive comparisons, we may conclude that the Winkler approach is adequate for practical
applications. Its advantage lies in the relative ease with which it may be applied to any
symmetric section.

Table 5.4.  The Values of m for Typical Ratios of Outer Radius b to Inner Radius a
Curved Beam Formula Elasticity Theory
b/a Flexure Formula r=a r=b r=a r=b
1.3 ±66.67 -72.980 61.270 -73.050 61.350
1.5 ±24.00 -26.971 20.647 -27.858 21.275
2.0 ±6.00 -7.725 4.863 -7.755 4.917
3.0 ±1.50 -2.285 1.095 -2.292 1.130

5.15   Comparison of the Results of Various Theories 293

 The agreement between the Winkler approach and the exact analyses is not as good
in situations involved combined loading as it is for the case of pure bending. As might be
expected, for beams of only slight curvature, the simple flexure formula provides good
results while requiring only simple computation. The linear and hyperbolic stress distri-
butions are approximately the same for b/a = 1.1. As the curvature of the beam increases
(b/a > 1.3), the stress on the concave side rapidly increases over the one given by the flex-
ure formula.

5.15.1  Correction of σ θ for Beams with Thin-Walled Cross Sections

Where I-beams, T-beams, or thin-walled tubular curved beams are involved, the
approaches developed in this chapter will not accurately predict the stresses in the system.
The error in such cases is attributable to the high stresses existing in certain sections such
as the flanges, which cause significant beam distortion. A modified Winkler’s equation
can be applied in such situations if more accurate results are required [Ref. 5.6]. The dis-
tortion, and thus the error in σ θ , is reduced when the flange thickness is increased. Given
that material yielding is highly localized, its effect is not of concern unless the curved
beam is under fatigue loading.

EXAMPLE 5.13  Maximum Stress in a Curved Rectangular Bar

A rectangular aluminum bar having mean radius r carries end moments M,
as illustrated in Fig. 5.26. Calculate the stresses in the member using (a) the
flexure formula and (b) the curved beam formula.
Given: M = 1.2 kN ⋅ m, b = 30 mm, h = 50 mm, and r = 125 mm.
Solution  The subscripts i and o refer to the quantities of the inside and out-
side fibers, respectively.

dr dr
h

M r M ro
ro ri r
ri

O O
(a) (b)

Figure 5.26.  E
 xample 5.13. (a) Rectangular curved beam in pure bending;
(b) cross section.

294 Chapter  5   Bending of Beams

 a. Applying the flexure formula, Eq. (5.38) with y = hl2, we obtain
My 1200(0.025)
σ o = −σ i = = 1 3 = 96 MPa
I 12 (0.03)(0.05)

This is the result we would get for a straight beam.

b. We first derive the expression for the radius R of the neutral axis. From
Fig. 5.26, A = bh and dA = bdr. Integration of Eq. (5.68) between the limits
ri and ro results in

A bh h
R= = =
dA bdr
ro ro dr
∫A r ∫ri r ∫ri r
or
h
R= (5.73)
r
In o
ri
The given data lead to
A = bh = (30)(50) = 1500 mm 2

ri = r − 12 h = 125 − 25 = 100 mm
ro = r + 12 h = 125 + 25 = 150 mm

Then, Eqs. (5.73) and (5.69) yield, respectively,

h 50
R= = = 123.3152 mm
ro ln 3
ln 2
ri
e = r − R = 125 − 123.3152 = 1.6848 mm

Note that the radius of the neutral axis R must be calculated with five
­significant figures.
The maximum compressive and tensile stresses are calculated by using
Eq. (5.70) as follows:

M ( R − ri ) 1.2(123.3152 − 100)
σi = =
Aeri 1.5(10−3 )(1.6848(10−3 )(0.1))
= −110.7 MPa
M ( R − ro ) 1.2(123.3152 − 150)
σo = =−
Aero 1.5(10−3 )(1.6848(10−3 )(0.15))
= 84.5 MPa

The negative sign means a compressive stress.

5.15   Comparison of the Results of Various Theories 295

 Comment  The maximum stress 96 MPa obtained in part (a) by the flexure
formula represents an error of about 13% from the more accurate value for
the maximum stress (110.7 MPa) found in part (b).

5.16  COMBINED TANGENTIAL AND NORMAL STRESSES

Curved beams are often loaded so that there is both an axial force and a moment on the
cross section. The tangential stress given by Eq. (5.70) may then be algebraically added to
the stress due to an axial force P acting through the centroid of cross-sectional area A. For
this simple case of superposition, the total stress at a point located at distance r from the
center of curvature O may be expressed as follows:
P M (R − r)
σθ = − (5.74)
A Aer
As before, a negative sign would be associated with a compressive load P.
Of course, the theory developed in this section applies only to the elastic stress dis-
tribution in curved beams. Stresses in straight members under various combined loads are
discussed in detail throughout this text.
The following problems illustrate the application of the formulas developed to stati-
cally determinate and statically indeterminate beams under combined loadings. In the lat-
ter case, the energy method (Section 10.4) facilitates the determination of the unknown,
redundant moment in the member.

CASE STUDY 5.1  Stresses in a Steel Crane Hook by Various Methods

A load P is applied to the simple steel hook having a rectangular cross section,
as illustrated in Fig. 5.27a. Find the tangential stresses at points A and B, using
(a) the curved beam formula, (b) the flexure formula, and (c) elasticity theory.
Given: P = 6 kN, r = 50 mm, b = 25 mm, and h = 32 mm.
Solution
a. Curved Beam Formula. For the given numerical values, we obtain
(Fig. 5.27b):

A = bh = (25)(32) = 800 mm 2
ri = r − 12 h = 50 − 16 = 34 mm
ro = r + 12 h = 50 + 16 = 66 mm

Then, Eqs. (5.73) and (5.69) result in

h 32
R= = = 48.2441 mm
ro 66
ln ln
ri 34
e = r − R = 50 − 48.2441 = 1.7559 mm

296 Chapter  5   Bending of Beams

 P

h = 32 mm P
r = 50 mm
M
C
B O O B A O
A B A
b = 25 mm
r = 50 mm

50 mm
P P
(a) (b) (c)

Figure 5.27.  Case Study 5.1. A crane hook of rectangular cross section.

To maintain applied force P in equilibrium, there must be an axial tensile

force P and a moment M = -Pr at the centroid of the section (Fig.  5.27c).
Then, by Eq. (5.74), the stress at the inner edge (r = ri) of the section A–B is

P ( − Pr )( R − ri ) P  r ( R − ri ) 
(σ θ ) A = − = 1 +  (5.75a)
A Aeri A eri 

6000  50(48.2441 − 34) 
= 1+ = 97 MPa
0.0008  (1.7559)(34) 

Likewise, the stress at the outer edge ( r = r ) is

P ( − Pr )( R − ro ) P  r ( R − ro ) 
(σ θ ) B = − = 1 +  (5.75b)
A Aero A ero 

6000  50(48.2441 − 66) 
= 1+ = −50 MPa
0.0008  (1.7559)(66) 

The negative sign of (σ θ ) B means a compressive stress is present. The

maximum tensile stress is at A and equals 97 MPa.

5.16   Combined Tangential and Normal Stresses 297

 Comment  The stress due to the axial force,
P 6000
= = 7.5 MPa
A 0.0008
is negligibly small compared to the combined stresses at points A and B of the
cross section.
b. Flexure Formula. Equation (5.5), with M = Pr = 6(50) = 300N ⋅ m, gives

My 3.00(0.01.6)
(σ θ ) B = −(σ θ ) A = = 1 3 = 70.3 MPa
I 12 (0025)(0032)

c. Elasticity Theory. Using Eq. (5.66) with a = ri = 34 mm and b = ro = 66 mm,

we find
2
 34 2  34 2 66
N = 1 −    − 4   ln 2   = 0.0726
  
66   66   34 

Superposition of -P/A and the second of Eqs. (5.67) with t = 25 mm at

r = a leads to

6000 4(300)  34 2  66 
(σ θ ) A = − +  1 − 2  (1 + 0) − (1 + 1) ln 
2
0.0008 (0.025)(0.066) (0.0726)  66  34 
= −7.5 − 89.85 = −97.4 MPa

Similarly, at r = b, we find (σ θ ) B = −7.5 + 58.1 = 50.6 MPa.

Comments  The results obtained with the curved-beam formula and with
elasticity theory are in good agreement. In contrast, the flexure formula pro-
vides a result with unacceptable accuracy for the tangential stress in this non-
slender curved beam.

EXAMPLE 5.14  Ring with a Diametral Bar

A steel ring of 350-mm mean diameter and of uniform rectangular section
60  mm wide and 12 mm thick is shown in Fig. 5.28a. A rigid bar is fitted
across diameter AB, and a tensile force P is applied to the ring as shown.
Assuming an allowable stress of 140 MPa, determine the maximum tensile
force that can be carried by the ring.
Solution  Let the thrust induced in bar AB be denoted by 2F. The moment at
any section a–a (Fig. 5.28b) is then

Pr
Mθ = − Fr sin θ + M B + (1 − cosθ ) (a)
2

298 Chapter  5   Bending of Beams

 P
C C
ds
dθ a
175 mm a
A B θ B
F
MB

12 mm P/2
r

D (b)
P
(a)

Figure 5.28.  E
 xample 5.14. (a) Ring with a bar AB is
subjected to a concentrated load P; (b)
moment at a section.
Note that before and after deformation, the relative slope between B and C
remains unchanged. Therefore, the relative angular rotation between B and C
is zero. Applying Eq. (5.32), we obtain
C π /2
EIθ = 0 = ∫ Mθ dx = r ∫ Mθ dθ
B 0

where dx = ds = rdθ is the length of beam segment corresponding to dθ .

Substituting in Eq. (a), this becomes, after integrating,
π
2 πMB
1
+ 12 Pr  − 1 − Fr = 0 (b)
2 

This expression involves two unknowns, M B and F. Another expression in

terms of M B and F is found by recognizing that the deflection at B is zero. By
applying Castigliano’s theorem,

∂U 1 π /2 ∂ Mθ
δB = =
∂ F EI ∫0 Mθ
∂F
( r dθ ) = 0

where U is the strain energy of the segment. This expression, upon introduc-
tion of Eq. (a), takes the form

∫0
π /2
{ − F r sin θ + M B +
Pr
2 }
(1 − cosθ ) sin θ dθ = 0

After integration,
− 14 π F r + 14 Pr + M B = 0 (c)
Solution of Eqs. (b) and (c) yields
M B = 0.1106 Pr

5.16   Combined Tangential and Normal Stresses 299

 and FR = 0.4591Pr . Substituting Eq. (a) gives

M C = − Fr + M B + 12 Pr = 0.1515Pr

Thus, M C > M B .
Since r /c = 0.175/0.006 = 29, the simple flexure formula offers the most
efficient means of computation. The maximum stress is found at points
A and B:
P /2 M B c
(σ θ ) A, B = + = 694P + 13,441P = 14,135P
A I
Similarly, at C and D,
Mcc
(σ θ )C , D = = 18,411P
I

Hence σ θ C > σ θ B , and we have σ max = 140 × 106 = 18,411P. The maximum
tensile load is therefore P = 7.604 kN.

REFERENCES

5.1. Timoshenko, S. P., and Goodier, J. N. Theory of Elasticity, 3rd ed. New York:
McGraw-Hill, 1970.
5.2. Ugural, A. C. Mechanics of Materials. Hoboken, NJ: Wiley, 2008, Chapter 8.
5.3. Sokolnikoff, I. S. Mathematical Theory of Elasticity, 2nd ed. Melbourne, FL:
Krieger, 1986.
5.4. Cook, R. D., and Young, W. C. Advanced Mechanics of Materials. New York:
Macmillan, 1985.
5.5. Budynas, R. G. Advanced Strength and Applied Stress Analysis, 2nd ed. New York:
McGraw-Hill, 1999.
5.6. Boresi, A. P., and Schmidt, R. J. Advanced Mechanics of Materials, 6th ed.
Hoboken, NJ: Wiley, 2003.

P RO B L E M S
Sections 5.1 through 5.5
5.1. A simply supported beam constructed of a 0.15 × 0.015 m angle is loaded by
concentrated force P = 22.5 kN at its midspan (Fig. P5.1). Calculate stress σ x at
A and the orientation of the neutral axis. Neglect the effect of shear in bending
and assume that beam twisting is prevented.
5.2. A wood cantilever beam with cross section as shown in Fig. P5.2 is sub-
jected to an inclined load P at its free end. Determine (a) the orientation
of the neutral axis and (b) the maximum bending stress. Given: P = 1 kN,
α = 30°, b = 80 mm, h = 150 mm, and length L = 1.2 m.

300 Chapter  5   Bending of Beams

 b

P
P 0.9 m 0.9 m
h C z
1.8 m
x z
A
A
y α
y P y

Figure P5.1. Figure P5.2.

30
C C
h z 8 z
Mz
30 My
α y
y Mo 8 60 8

Figure P5.3. Figure P5.4.

24 kN 24 kN 24 kN
2m 2m
4m
D E

D, E z
C
y

4m

Figure P5.5.

5.3. A moment M o is applied to a beam of the cross section shown in Fig. P5.3 with
its vector forming an angle α . Use b = 100 mm, h = 40 mm, M o = 800N ⋅ m, and
α = 25°. Calculate (a) the orientation of the neutral axis and (b) the maximum
bending stress.
5.4. Couples M y = M o and M z = 1.5M o are applied to a beam of cross section shown
in Fig. P5.4. Determine the largest allowable value of M o for the maximum stress
not to exceed 80 MPa. All dimensions are in millimeters.
5.5. For the simply supported beam shown in Fig. P5.5, determine the bending stress
at points D and E. The cross section is a 0.15 × 0.15 × 0.02 m angle (Fig. 5.4).
5.6. A concentrated load P acts on a cantilever, as shown in Fig. P5.6. The
beam is constructed of a 2024-T4 aluminum alloy having a yield strength

Problems 301
 L α
P t
P
x b z
y
t
y c

Figure P5.6.

P
L h
P
t h
x
C z
h

y
y

Figure P5.8.

A Mo

α
80 mm C z
D

B y
90 mm

Figure P5.9.

σ yp = 290 MPa, L = 1.5 m, t = 20 mm, c = 60 mm, and b = 80 mm. Based on a

factor of safety n = 1.2 against initiation of yielding, calculate the magnitude of P
for (a) α = 0° and (b) α = 15°. Neglect the effect of shear in bending and assume
that beam twisting is prevented.
5.7. Re-solve Prob. 5.6 for α = 30°. Assume the remaining data are unchanged.
3 3
5.8. A cantilever beam has a Z section of uniform thickness for which I y = 23 th , I z = 83 th
3 8 3
I y = 3 th , I z = 3 th and I yz = −th . Determine the maximum bending stress in the beam sub-
2 3

jected to a load P at its free end (Fig. P5.8).

5.9. A beam with cross section as shown in Fig. P5.9 is acted on by a moment
M o = 3kN ⋅ m, with its vector forming an angle α = 20°. Determine (a) the orien-
tation of the neutral axis and (b) the maximum bending stress.
5.10 and 5.11. As shown in the cross section in Figs. P5.10 and P5.11, a beam carries a
moment M, with its vector forming an angle α with the horizontal axis. Find the
stresses at points A, B, and D.

302 Chapter  5   Bending of Beams

 M = 10 kN·m
D

α = 40°

60 mm
80 mm
C

A B
100 mm

120 mm

Figure P5.10.
y

M = 4 kN·m
75 mm
A
α = 35°
z C 50 mm

75 mm
B

100 mm

50 mm

Figure P5.11.
po t

x z h

L
y
y

Figure P5.12.

5.12. For the thin cantilever shown in Fig. P5.12, the stress function is given by

x3 x3 y xy 3 x3 y3 xy 5
Φ = − c1 xy + c2 − c3 − c4 − c5 − c6
6 6 6 9 20
a. Determine the stresses σ x ,σ y, and τ xy by using the elasticity method.
b. Determine the stress σ x by using the elementary method.
c. Compare the values of maximum stress obtained by the preceding approaches
for L = 10 h.

Problems 303
 5.13. Consider a cantilever beam of constant unit thickness subjected to a uniform load
of p = 2000 kN per unit length (Fig. P5.13). Determine the maximum stress in
the beam:
a. Based on a stress function

p  2 y 
Φ= − x + xy + ( x 2 + y 2 )  0.78 − tan −1  
0.43   x
b. Based on the elementary theory. Compare the results of (a) and (b).

Sections 5.6 through 5.11

5.14. A bending moment acting about the z axis is applied to a T-beam, as shown in
Fig. P5.14. Take the thickness t = 15 mm and depth h = 90 mm. Determine the
width b of the flange needed so that the stresses at the bottom and top of the
beam will be in the ratio 3:1, respectively.
5.15. A wooden, simply supported beam of length L is subjected to a uniform load
p. Determine the beam length and the loading necessary to develop simultane-
ously σ max = 8.4 MPa and τ max = 0.7 MPa. Take thickness t = 0.05 m and depth
h = 0.15 m.
5.16. A box beam supports the loading shown in Fig. P5.16. Determine the maximum
value of P such that a flexural stress σ = 7 MPa or a shearing stress τ = 0.7 will
not be exceeded.

a
2 MN/m

x b
t
y C
a z
h

t
y

Figure P5.13. Figure P5.14.

0.2
0.15
P 4.5 kN

x 0.25 0.2
z
A B

1.5 m y
1.5 m 1.5 m
y

Figure P5.16.

304 Chapter  5   Bending of Beams

 y

x 2b
A
B
L
b

Figure P5.17.

A B h
h h0

x b
L L
2 2

Figure P5.18.

y
t

z h h
M C

t
t

Figure P5.19.

5.17. Design a rectangular cantilever beam of constant strength and width b, to carry a
uniformly distributed load of intensity w (Fig. P5.17). Assumption: Only the nor-
mal stresses due to the bending need be taken into account; the permissible stress
equals σ all.
5.18. Design a simply supported rectangular beam of constant strength and width b,
supporting a uniformly distributed load of intensity w (Fig. P5.18). Assump-
tion: Only the normal stresses due to the bending need be taken into account; the
allowable stress is σ all.
5.19. A steel beam of the tubular cross section seen in Fig. P5.19 is subjected to
the bending moment M about the z axis. Determine (a) the bending moment
M and (b) the radius of curvature rx of the beam. Given: σ all = 150 MPa, E = 70 GPa, b = 120 mm, h =
σ all = 150 MPa, E = 70 GPa, b = 120 mm, h = 170 mm, E = 70 GPa, b = 120 mm, h = 170 mm,
and t = 10 mm.

Problems 305
 x
y

z
M
C

B
A
d
D

Figure P5.20.

P 400 mm
t y

A B z C 150 mm
t y
t = 25 mm
1.2 m 1.2 m

Figure P5.21.

h
A B
C
L/2
L b

Figure P5.22.

5.20. An aluminum alloy beam of hollow circular cross section is subjected to a bend-
ing moment M about the z axis (Fig. P5.20). Determine (a) the normal stress at
point A, (b) the normal stress at point B, and (c) the radius of curvature rx of
the beam of a transverse cross section. Given: M = 600 N ⋅ m, D = 60 mm, d =
40 mm, E = 70 GPa, and v = 0.29.
5.21. A simply supported beam AB of the channel cross section carries a concentrated
load P at midpoint (Fig. P5.21). Find the maximum allowable load P based on an
allowable normal stress of σ a11 = 60 MPa in the beam.
5.22. A uniformly loaded, simply supported rectangular beam has two 15-mm deep
vertical grooves opposite each other on the edges at midspan, as illustrated in
Fig. P5.22. Find the smallest permissible radius of the grooves for the case in which
the normal stress is limited to σ max = 95 MPa . Given: p = 12 kN/m, L = 3 m, b = 80 mm,
p = 12 kN/m, L = 3 m, b = 80 mm, and h = 120 mm.

306 Chapter  5   Bending of Beams

 p

A B h

b
L

Figure P5.23.

t b
y

h
z

Figure P5.25.
5.23. A simple wooden beam is under a uniform load of intensity p, as illustrated
in Fig. P5.23. (a) Find the ratio of the maximum shearing stress to the larg-
est bending stress in terms of the depth h and length L of the beam. (b) Using
σ a11 = 9 MPa,τ a11 = 1.4 MPa, b = 50 mm, and h = 160 mm, calculate the maximum
permissible length L and the largest permissible distributed load of intensity p.
5.24. A composite cantilever beam 140 mm wide, 300 mm deep, and 3 m long is fab-
ricated by fastening two timber planks, 60 mm × 300 mm, to the sides of a steel
plate ( Es = 200 GPa), 20 mm wide by 300 mm deep. Note that the 300-mm
dimension is vertical. The allowable stresses in bending for timber and steel are 7
and 120 MPa, respectively. Calculate the maximum vertical load P that the beam
can carry at its free end.
5.25. A simple beam pan length 3 m supports a uniformly distributed load of 40 kN/m.
Find the required thickness t of the steel plates. Given: The cross section of
the beam is a hollow box with wood flanges (Ew = 10.5 GPa) and steel (Es =
210 GPa), as seen in Fig. P5.25. Let a = 62.5 mm, b = 75 mm, and h = 225 mm.
Assumptions: The permissible stresses are 140 MPa for the steel and 10 MPa for
the wood.
5.26. A 180-mm-wide by 300-mm-deep wood beam ( Ew = 10 GPa), 4 m long, is rein-
forced with 180-mm-wide and 10-mm-deep aluminum plates ( Ea = 70 GPa) on
the top and bottom faces. The beam is simply supported and subject to a uniform
load of intensity 25 kN/m over its entire length. Calculate the maximum stresses
in each material.

Problems 307
 y s

D
Vy α
Vy
c=h
z
S
z
S C R
h t

e
e b t y

Figure P5.29. Figure P5.30.

e1 e2
P
t2
t1 C
b1 b2
z
S
y
h

Figure P5.31.

5.27. Referring to the reinforced concrete beam of Fig. 5.17a, assume b = 300, d = 450 mm, As = 1200 m
b = 300, d = 450 mm, As = 1200 m 2 , and n = 10. Given allowable stresses in steel and con-
crete of 150 and 12 MPa, respectively, calculate the maximum bending moment
that the section can carry.
5.28. Referring to the reinforced concrete beam of Fig. 5.17a, assume b = 300 mm,
d = 500 mm, and n = 8. Given the actual maximum stresses developed to be
σ s = 80 MPa and σ c = 5MPa, calculate the applied bending moment and the steel
area required.
5.29. A channel section of uniform thickness is loaded as shown in Fig. P5.29. Find (a)
the distance e to the shear center, (b) the shearing stress at D, and (c) the maxi-
mum shearing stress. Given: b = 100 mm, h = 90 mm, t = 4 mm, Vy = 5 kN.
5.30. A beam is constructed of half a hollow tube of mean radius R and wall thickness t
(Fig. P5.30). Assuming t  R, locate the shear center S. The moment of inertia of
the section about the z axis is I z = π R 3tl 2.
5.31. An H-section beam with unequal flanges is subjected to a vertical load P
(Fig. P5.31). The following assumptions are applicable:
1. The total resisting shear occurs in the flanges.
2. The rotation of a plane section during bending occurs about the symmetry axis
so that the of curvature of both flanges are equal.
Find the location of the shear center S.

308 Chapter  5   Bending of Beams

 75

15.5

9.5 250

po
z
C
x
y A B
9.5 L
125 y

Figure P5.32. Figure P5.33.

pL3/96El
p

L
y
(a) (b)

Figure P5.34.
Mo
A B
x

L
y

Figure P5.36.

5.32. Determine the shear center S of the section shown in Fig. P5.32. All dimensions
are in millimeters.
5.33. A cantilever beam AB supports a triangularly distributed load of maximum inten-
sity po (Fig. P5.33). Determine (a) the equation of the deflection curve, (b) the
deflection at the free end, and (c) the slope at the free end.
5.34. The slope at the wall of a built-in beam (Fig. P5.34a) is as shown in Fig. P5.34b
and is given by pL3 /96EI . Determine the force acting at the simple support,
expressed in terms of p and L.
5.35. A fixed-ended beam of length L is subjected to a concentrated force P at a dis-
tance c away from the left end. Derive the equations of the elastic curve.
5.36. A propped cantilever beam AB is subjected to a couple M o acting at support B,
as shown in Fig. P5.36. Derive the equation of the deflection curve and determine
the reaction at the roller support.

Problems 309
 po

A B
C
L L
2 2

Figure P5.37.

1 1

2h

2 h

Figure P5.38.
5.37. A clamped-ended beam AB carries a symmetric triangular load of maximum inten-
sity p0 (Fig. P5.37). Find all reactions, the equation of the elastic curve, and the
maximum deflection, using the second-order differential equation of the deflection.
5.38. A welded bimetallic strip (Fig. P5.38) is initially straight. A temperature incre-
ment ∆T causes the element to curve. The coefficients of thermal expansion of
the constituent metals are α 1 and α 2 . Assuming elastic deformation and α 2 > α 1
determine (a) the radius of curvature to which the strip bends, (b) the maximum
stress occurring at the interface, and (c) the temperature increase that would
result in the simultaneous yielding of both elements.

Sections 5.12 through 5.16

5.39. Verify the values of α for cases B, C, and D of Table 5.2.
5.40. Consider a curved bar subjected to pure bending (Fig. 5.24). Assume the stress
function
Φ = A ln r + Br 2 ln r + Cr 2 + D
Re-derive the stress field in the bar given by Eqs. (5.67).
5.41. The allowable stress in tension and compression for the clamp body shown in
Fig. P5.41 is 80 MPa. Calculate the maximum permissible load that the member
can resist. Dimensions are in millimeters.

46

a a 30
P
120
30 6 6

Section a-a

Figure P5.41.

310 Chapter  5   Bending of Beams

 P

A O C
B b O
r B A

Figure P5.42.

b
B 10 mm 10 mm

60 mm

20 mm
140 mm
A
60 mm Section A-B
M M
O

Figure P5.43.

B B

d Section A-B

A d A
O
b
P

Figure P5.44.

5.42. A curved frame of rectangular cross section is loaded as shown in Fig. P5.42.
Determine the maximum tangential stress (a) by using the second of Eqs. (5.67)
together with the method of superposition and (b) by applying Eq. (5.73). Given:
h = 100 mm, r = 150 mm, and P = 70 kN.
5.43. A curved frame having a channel-shaped cross section is subjected to bending by
end moments M, as illustrated in Fig. P5.43. Determine the dimension b required
if the tangential stresses at points A and B of the beam are equal in magnitude.
5.44. A curved beam of a circular cross section of diameter d is fixed at one end and
subjected to a concentrated load P at the free end (Fig. P5.44). Calculate (a) the
tangential stress at point A and (b) the tangential stress at point B. Given: P =
800 N, d = 20 mm, a = 25 mm, and b = 15 mm.

Problems 311
 P

ri
b2 B A b1 B A O

Section A-B ro

Figure P5.45.

ro

h ri

b O

Figure P5.46.

5.45. A circular steel frame has a cross section approximated by the trapezoidal form
shown in Fig. P5.45. Calculate (a) the tangential stress at point A and (b) the
tangential stress at point B. Given: ri = 100 mm, b = 75 mm, b = 50 mm, and
P = 50 kN.
5.46. The triangular cross section of a curved beam is shown in Fig. P5.46. Derive the
expression for the radius R along the neutral axis. Compare the result with that
given for Fig. D in Table 5.3.
5.47. The circular cross section of a curved beam is illustrated in Fig. P5.47. Derive the
expression for the radius R along the neutral axis. Compare the result with that
given for Fig. B in Table 5.3.
5.48. The trapezoidal cross section of a curved beam is depicted in Fig. P5.48. Derive
the expression for the radius R along the neutral axis. Compare the result with
that given for Fig. E in Table 5.3.
5.49. A machine component of channel cross-sectional area is loaded as shown in Fig.
P5.49. Calculate the tangential stress at points A and B. All dimensions are in
millimeters.
5.50. A load P is applied to an eye bar with rigid insert for the purpose of pulling
(Fig. P5.50). Determine the tangential stress at points A and B (a) by the elastic-
ity theory, (b) by Winkler’s theory, and (c) by the elementary theory. Compare the
results obtained in each case.

312 Chapter  5   Bending of Beams

 r–
h ri
ri

c w
b2 b1 O
O
C

dr
r
ro ro

Figure P5.47. Figure P5.48.

P 60

40
40
A B
10
30 5
Section A-B

P = 3 kN

Figure P5.49.

P P

b b t
2

Figure P5.50.

5.51. A ring of mean radius r and constant rectangular section is subjected to a concen-
trated load (Fig. P5.51). You may omit the effect of shear in bending. Derive the
following general expression for the tangential stress at any section of the ring:

( P /2)cosθ Mθ  R − r  (P5.51)
σθ = − +  
A A er 

Problems 313
 P

h
r
θ
t

Figure P5.51.
where
Mθ = 0.182Pr − 12 Pr (1 − cosθ )

Use Castigliano’s theorem.

5.52. The ring shown in Fig. P5.51 has the following dimensions: r = 150 mm,
t = 50 mm, and h = 100 mm. Taking E = 25 G , determine (a) the tangential stress
on the inner fiber at θ = π /4 and (b) the deflection along the line of action of the
load P, considering the effects of the normal and shear forces, as well as bending
moment (Section 10.4).

314 Chapter  5   Bending of Beams

 Index

A Bar springs, 348–351 general theory for, 473–474

Bars, torsion. see Torsion, grid configurations for,
Absolute maximum shear stress, 46
prismatic bars 490–492
Adhesion, JKR theory, 170–171
Bathe, K.I., 423, 673 infinite, 475–480
Advanced material, 85
Baumeister III, T., 117, 711 semi-infinite, 480–483
Airy’s stress function, elasticity
Beam column, 554, 556 solution by finite differences,
problems, 141–143
Beams, bending 488–490
Allowable stress, 7
analysis of perfectly plastic, solutions for relatively stiff,
Almroth, B.O., 569
590–600 486–488
Aluminium, 85
applying Castigliano’s theorem supported by elastic elements,
Angle of twist, 317, 324, 335–337,
to, 500, 510–511 485–486
348–349
asymmetrical section of, Becker, S. J., 466
Angular speed, rotating disks, 446
246–251 Bending
Anisotropy, 84, 96
cantilever, 145–146, 251–254, beams. see Beams, bending
Annular disk, rotating, 447, 459
269, 278 charts of bars/shafts in, 707–710
Antiplane strain, 142–143
composite, 270–276 elementary bending theory. see
Antiplastic surfaces, 246
of constant strength, 267–268 Elementary bending theory
Asymmetrical cross section,
curvature of, 245 impact, failure criteria, 227–230
beams, 246–251
curved, 286–299 pure bending. see Pure bending
Autofrettage, 623
elasticity theory applied to, of thin plates, 636, 640
Avallone, E. A., 117, 711
286–289 Bernoulli-Euler law, 246
Axial
kinematic relations of, 244–246 Betti, E., 499
biaxial tension, 621
normal and sheer stresses on, Bifurcation point, 535
loading, 227
260–268 Biharmonic equation, 141
rigidity, 94
numerical methods for, 393–399 Body forces, 2
triaxial stress, 12–13, 42, 45,
plastic deflection, 588–590 Boley, B. A., 181
98, 200
pure bending, 289 Bonded gages, foil and wire,
uniaxial stress, 13, 92–93, 218–221
shear center and, 276–281 103–104
uniaxial tension, 108
strain energy in, 113, 284–286 Boresi, A. P., 51, 117, 181, 300, 423
Axially loaded bars, 111–112,
stresses in sandwich, 653–654 Boundary conditions
387–389
table of deflections and, applied to plates, 642–644
Axis of symmetry, 689, 692
705–706 elasticity theory and, 287–288
Axisymmetric elements
tangential and normal stresses numerical methods for, 373–377
finite element solution for,
on, 296–300 three-dimensional problems
464–465
Timoshenko beam theory, 246 and, 49–50
stresses in circular plates, 648–650
under torsion, 348–349, in torsion problems, 326–327,
Azar, J .J., 568
373–377 348
transverse normal stress on, in two-dimensional problems,
B 268–270 136, 138
Bar elements Beams, elastic foundations Boundary value problem, elasticity,
Castigliano’s theorem for, 500 classification of, 484–485 133
numerical method for, 379–383 finite, 483–484 Brady, G. S., 117

722
Bredt’s formulas, 341 Classification, of columns, 543–544 Conversion, SI Unit/U.S. Unit,
Bridget, F. J., 569 Clauser, H. R., 117 701–704
Brittle material, 86, 88 Cold working, 92, 194 Cook, R. D., 300
Broek, D., 232 Collapse load, plastic, 600–605 Coordinate transformation, bar
Brush, D. O., 569 Columns element, 380–381
Buckling buckling of pin-ended, Coulomb, C. A., 198
of columns, 549–550, 552 536–539 Coulomb-Mohr theory, failure
energy methods applied to, critical load of, 534–535 criteria, 207–210
554–562 deflection response of, 539–540 Cozzarelli, F. A., 118
by finite differences, 562–567 design formulas for, 548–550 Crack, fracture mechanics, 210–213
intermediate columns and with different end conditions, Creep curve, 194–195
inelastic, 544 540–543 Critical load, columns, 534–535,
load, 535 eccentrically loaded, 552–554 560–562
modes of, 538–539 energy methods applied to Critical stress, classification of
of pin-ended columns, 536–538 buckling, 554–562 columns, 543–544
Buckling formula, Euler’s, 538 imperfections in, 550–552 Crotti-Engesser theorem, 508–509
Budynas, R. G., 181, 300, 354, 711 inelastic behavior and, 540 CST (constant strain triangular),
Bulk modulus, 98 Johnson’s buckling criterion for, 407
Burr, A. H., 524 546–548 Curvature, radius of, 245
local buckling of, 552 Curvature of plates, 636–639
long, short and intermediate, Curved beam formula, 289–293
C 543–544 Curved boundaries, numerical
CAD (computer-aided design) solution by finite differences, methods for, 370–373
software, 423 562–568 Cylinders
CAE (computer-aided tangent modulus theory and, filament-wound, 666–667
engineering), 423 544–546 hydraulic, 34–35, 445
Calculations, engineering, 680 Combined stress, 17 Cylinders, thick walled
Cantilever beams Compatibility equations. see compound, 443–446
bending of, 251–254, 259, 278 Equations of compatibility failure theories for, 442–443
Castigliano’s theorem for, Complete yielding, plastic, 614, finite element solution for,
510–511 624–626 464–466
plane stress in, 414–418 Composite areas, moments of maximum tangential stress,
solving deflection problems, centroid, 689–691 441–442
373–374, 519–522 inertia, 692–700 under pressure, 435–441
Cartesian representations, plane Composite beams stresses in perfectly plastic,
stress, 25–26 equation of neutral axis, 623–627
Case studies 271–272 thermal stress in long circular,
in plane stress, 414–422 of multi materials, 272–276 460–464
in analysis, 7–8 stresses in transformed, 272 Cylindrical contacts, elasticity,
Cast irons, 84–85 transformed section method for, 173–174
Castigliano’s theorem 270–271 Cylindrical shells, 668–673
applied to bars, beams and Composites, material, 85
trusses, 500–501 Compound cylinders, 443–446,
applied to strain energy, 668–670 D
284–286 Computational tools, numerical Da Vinci, Leonardo, 2
fictitious loads and, 501–506 methods, 423 Deflection curve, Rayleigh-Ritz
Central differences, numerical Computer-aided design (CAD) method, 522
methods, 366–368 software, 423 Deflections
Ceramics, 85 Computer-aided engineering of a beam column, 556–557
Cheatham, J. B., 524 (CAE), 423 of beams, 373–377
Chong, K. P., 51 Constant strain triangular (CST), for columns, 539–540
Chou, P. C., 51, 117 407 contact stresses and, 169–171
Circular Constantan, 104 elasticity theory and, 289
axisymmetrically loaded plates, Contact energy method for, 284–286
648–650 general, 178–181 equations of plate, 640–641
bars in torsion, 112–113, 316– pressure for compound perfectly plastic simple beam,
321, 351–354 cylinders, 443 599–600
contact, 171 spherical and cylindrical, 171–174 plastic deflection of beams,
cylinders and thermal stresses, stress distribution, 174–178 588–590
460–464 stresses and deflections, 169–171 of rectangular plates, 650–652

Index 723
Deflections (cont.) Elastic/elasticity for the bar element, 379–380
of a ring, 512 Airy’s stress function, 141–143 Castigliano’s theorem, 499–506
and slopes of beams, 705 basic relations in polar Crotti-Engesser theorem, 508–509
of statically indeterminate coordinates, 152–157 for deflection of beams,
beams, 706 beam foundations. see Beams, 284–286
by trigonometric series, 517–522 elastic foundations deflections by trigonometric
Deformation constants, 96 series, 517–522
analysis of stress and, 5 contact stress distribution, principle of minimum potential
beam kinematics and, 244–245 174–178 energy, 515–517
of plastic, 579–580 contact stresses and deflections, principle of virtual work,
plastic axial, 585–588 169–171 514–515
of a plate in bending, 636 elements, beams supported by, Rayleigh-Ritz method, 522–524
torsion problems and, 324–325 485–486 reciprocity theorem, 498–499
work done in, 497 general contact, 178–181 stability of columns and, 536
Deformational theory, 614 plane strain, 135–137 statically indeterminate systems,
Degree of statical indeterminacy, plane stress, 138–140 510–513
282. see also Statically versus plastic material, 91–92 unit-or dummy-load method,
indeterminate systems solutions for problems in, 506–508
Delta rosette, 106–107 143–148 work done in deformation, 497
Design, mechanics of solids, 4–8 spherical and cylindrical Engesser formula, 545
Design formulas, columns, 548–550 contacts, 171–174 Engineering calculations, 680
Deviator (distortional) stresses, strain, stress and matrices, Engineering materials, table, 702–703
114, 579 401–402 Equation of neutral axis, composite
Differential equations of stress concentration factors, beams, 271–272
equilibrium, 20 163–169 Equations. see also Formulas
Dilatation, 93, 98 stress distribution acting on a axial force, 384–386
Dilatational stress tensor, 113 beam, 161–163 cubic equation, 39–40, 682–686
Direction cosines, 44, 683–686 stresses due to concentrated for finite differences, 368–370
Disks, rotating loads, 157–161 Hencky’s equations, 616–619,
of constant thickness, 446–449 thermal stresses and, 149–152 620
disk flywheel, 449–453 three-dimensional problems, 134 stiffness matrix, 383
elastic-plastic stresses in, two-dimensional problems, of thermoelasticity, 149–152
612–614 134–135, 140–141 transformation equations,
thermal stresses in thin, 458–459 Elasticity theory, beams, 286–289 155–156, 248
of uniform stress, 456–458 Elastic-plastic Equations of compatibility
of variable thickness, 453–456 strain relations, 621 in beams, 286–287
Dislocation, 579–580 stresses in rotating disks, 612–614 in polar coordinates, 156–157
Displacement torque, 606 in solving torsion problem, 325
finite element analysis and, torsion of circular shafts, strain and, 75–76
377–378 605–609 in thermoelasticity, 149–152
finite element solution for thin Element nodal forces, 410–414 in two-dimensional problems,
plates, 655–657 Elementary bending theory 137, 141
matrix for two-dimensional comparison of results, 253–254, Equations of equilibrium
elements, 399–400 256 in beams, 286–287
Rayleigh-Ritz method, 522 conclusions about, 256–257 of plate deflection, 640–641
transformation for bar element, method of integration, 258–260 in polar coordinates, 153
383 Elementary formulas. see Prandtl’s membrane analogy,
Ductile material, 86–87, 89, 193, Mechanics of materials 333–335
230–232 Elementary theory of torsion, of shells, 661–662
Dummy load, energy method, circular bars, 316–321 in solving torsion problem, 325
506–508 Ellipsoid, stress, 44 in two-dimensional problems,
Dynamic loading, 225–226 Elliptical bars, under torsion, 136, 138
328–331, 372–373. see also Equidimensional equation, 455
Torsion, prismatic bars Equilateral hyperbola, 454
E Empirical formulas, columns, 548 Equilateral triangle bar, under
Eccentrically loaded columns, Endurance limit, 217 torsion, 331–332
552–554 Energy, strain, 107–110 Equilibrium
Edge dislocation, 579 Energy methods conditions of, 8–9
Effective stress, 199, 615 applied to buckling columns, equations. see Equations of
Eigenfunctions, 539 554–562 equilibrium

724Index
Equilibrium method Finite beams, elastic foundations, Newton’s interpolation formula,
applied to columns, 535 483–484 369
for the bar element, 379 Finite differences secant formula, 552–554, 615
Newton’s equilibrium method, 514 for beams on elastic transfer formula, 692–694
Equivalent alternating stress, 221 foundations, 488–490 Winkler’s formula, 293
Euler, Leonard, 2 boundary conditions, 373–377 Formulation, problem/solution,
Eulerian coordinates, 73–74 central differences, 366–368 679–681
Euler’s buckling formula, 541 columns and solution by, Fourier series, 517–518, 644
Euler’s curve, 544 562–568 Fracture
Exact theory, beams, 487 curved boundaries, 370–373 criteria, 197–198
Finite element analysis (FEA) failure by, 195–197
axial force equation, 384–386 mechanics, 210–213
F for the bar element, 379–383 progressive, 216–217
Faces, sandwich plates, 652 for a beam element, 393–399 toughness, 213–215
Factor of safety, 6–7, 213, 546–549 case studies in plane stress, Free-body diagram (FBD), 10, 679
Failure 414–422 Fundamental principals, of
comparison of the yielding computational tools, 423 elasticity, 134–135
theories, 204 for plates, 654–657
Coulomb-Mohr theory, 207–210 for a triangular finite element, G
criteria for metal fatigue, 217–223 407–414
Gages, strain, 103–105
ductile-brittle transition, for a truss, 386–393
Galilei, Galileo, 2
230–232 for two-dimensional elements,
General contact, elasticity,
fatigue life, 223–225 399–402
178–181
fatigue: progressive fracture, Finite element method (FEM), 3,
Generalized Hooke’s law, 96–100.
216–217 402–407
see also Hooke’s law
by fracture, 195–197 Fisher, F. E., 118, 181, 466, 524,
Generalized plane strain, 142
fracture toughness, 213–215 627, 673
Generator, shell, 670
impact loads, 225–227 Flexural
Geometry, of deformation. see also
introduction to fracture center, 276–277
Deformation
mechanics, 210–213 formula, 244
analysis of stress and, 5
longitude and bending impact, loading, 227
beam kinematics and, 244–245
227–230 rigidity, 244, 541, 639
torsion problems and, 324–325
maximum distortion energy Flugge, W., 492
Gerber criterion, fatigue failure,
theory, 199–200 Flywheel, disks, 449–453
219
maximum principal stress Foil gage, 103–104
Gere, J. M., 524, 568
theory, 204–206 Force transformation, bar element,
Goodier, J., 118, 181, 300, 354
maximum shearing stress 381–382
Goodman criterion, for uniaxial
theory, 198–199 Forces. see also by individual
stress, 219–220
Mohr’s theory, 206–207 types
Grid configuration, beams,
octahedral shearing stress axial force equation, 384–386
490–492
theory, 200–203 boundary conditions of surface
Griffith’s theory, failure by
theories. see Theories, failure forces, 49–50
fracture, 195
yield and fracture criteria, force transformation in bar
197–198 element, 381–382
by yielding, 193–195 force-displacement relations, H
Fatigue 386–393 Helical springs, 352–354
criterion table, 218 internal force resultants, 13–17 Hencky, H., 199
failure criteria for metal, 217–223 Kirchhoff ’s force, 642 Hencky’s equations, 616–619, 620
failure diagram, 219 nodal forces, 379, 410–414 Hertz, H., 169, 181
life, 223–225 Form factor, for shear, 285–286 Hertz Theory, contact stresses,
strength and endurance limits, Formulas. see also Equations 170–171
217 Bredt’s formulas, 341 Hetényi, M., 118, 492
tests, 216–217 column, 548–550, 554 Hinge, plastic, 593–601, 604
Faupel, J. H., 118, 181, 466, 524, curved beam, 289–293 Hodge, P. G., 627
627, 673 Engesser formula, 545 Hoffman, O., 627
Fictitious loads, Castigliano’s Euler’s buckling formula, 538, Hooke, Robert, 2
theorem, 501–506 541 Hooke’s law
Filament-wound cylinder, 666–667 finite element method, 402–407 for beams in bending, 244, 252
Filleted bars, stress concentration flexural, 244, 259, 272 for bending of thin plates, 638
charts, 707 Johnson’s buckling criterion, 546 generalized, 96–100

Index 725
Hooke’s law (cont.) J Loading rate, ductile-brittle
for orthotropic material, transition, 231
Johnson, K. L., 181
102–103 Local buckling, of columns, 552
Johnson-Kendall-Roberts (JKR)
for plane stress, 139–140 Logan, D. L., 423
theory, 170–171
in polar coordinates, 155 Logarithmic strain, 88
Johnson’s buckling criterion,
for uniaxial stress, 92–93 Long beams, 485
546–548
Huber, T. M., 199 Long columns, 544
Hydraulic cylinders, 34–35, Longitude impact, failure criteria,
445 K 227–230
Hydrodynamic analogy, 344–346 Kelvin, William Thomas, 2 Love, A. E. H., 2, 51
Hydrostatic stress, 44 Kendall, K., 181 LST (linear strain triangular)
Hyperbolic disk, 455–456 Kinematic relations. see Strain- elements, 407
displacement relations
Kirchhoff ’s force, 642 M
I
Knite, E., 423
Impact loads. see also Load Kotter, T., 118 Magnesium alloys, 85
basic assumptions, 227 Kreyzig, E., 423, 524 Marin, J., 232
bending, freely falling weight, Material
227–228 brittle. see Brittle material
bending, horizontally moving
L ductile. see Ductile material
weight, 228–230 Lagrange, Joseph-Louis, 2 elastic versus plastic, 91–92
strain rate and, 226 Lagrangian coordinates, 73–74 engineering materials table,
Imperfect column, 540 Lagrangian energy method, 514 702–703
Incremental theory, 620 Lame constants, 98 isotropic, 97
Indicial notation, analysis of stress, Large strain, 74 linearly elastic. see Linearly
50 Lateral shears, 347 elastic material
Inelastic behavior, columns, 540 Levy-Mises theory, 621–622 mechanics of. see Mechanics of
Inelastic buckling, 544 Limit design, plastic loads, 600–605 materials
Inertia moments Linear elasticity, 92 multi materials beams
parallel axis theorem, 692–694 Linear strain triangular (LST) (composite), 272–276
principal moments of inertia, elements, 407 orthotropic, 101–103, 139–140
694–700 Linearly elastic material plastic behavior. see Plastic/
Infinite beams, elastic foundations, beams, 473 plasticity
475–480 Castigliano’s theorem and, 499 MATLAB (Matrix Laboratory),
Inglis, C. E., 167 deformation and, 69 423, 712–714
Initial yielding, plastic, 612 unit load method and, 506 Matrices (strain, stress and
Instability phenomenon, plastic, Load elasticity), 655
582–583 axially loaded, 111–112, 387– Matrix methods, finite element
Interaction curves, elastic-plastic, 389, 585–588 analysis, 377–378
595–597 axisymmetrically loaded, Maugis, D., 181
Interface pressure, compound 648–650 Maximum
cylinders, 443 buckling, 535 distortion energy theory,
Intermediate beams, 485 collapse, 600–605 199–200
Internal force resultants, 13–17 critically loaded columns, principal stress theory, 204–206
Internal forces, 2 534–535, 560–562 radial stress, 448
Internal friction theory, 207–210 eccentrically loaded columns, shearing stress, 198–199
International System of Units (SI), 552–554 tangential stress, 441–442
2, 701–704 fictitious loads and Castigliano’s M-code, MATLAB, 712–713
Inverse method, 143 theorem, 501–506 Mean strain, 114
Iron, cast iron, 84–85 Mohr’s circle for torsional Mean stress, 218–219
Irwin, G. R., 232 loading, 322 Measurement, strain, 103–107
Isotropic repeated loading, 216 Mechanics of materials
bending of thin plates, 636 stiffness matrix for deflection problems in elasticity and, 133
engineering materials, 84 bar under combined loading, stress formulas, 15–17
materials and generalized 405–406 theory of elasticity and, 1–2
Hooke’s law, 97 symmetrically loaded shells, Mechanism of collapse, 600–601
sandwich plates, 654 660–662 Meguid, S. A., 232
two-dimensional problems, ultimate load, 601 Membrane action, shells, 658–660
140–141 unit-or dummy-load method, Membrane analogy. see Prandtl’s
Iyengar, K .T., 492 506–508 membrane analogy

726Index
Membrane-roof analogy, 610–611 Neutral Orthogonality function, 518–519
Mendelsson, A., 627 axes, 244 Orthotropic materials, strain,
Meridian curve, shells, 660 axis equation, composite beams, 101–103, 139–140
Mesh refinement, 378 271–272 Osgood, W. R., 581, 627
Messal, E. E., 686 equilibrium, 535 Out-of-plane principal strain, 139
Metal fatigue, failure criteria, Newton, Sir Isaac, 2
217–223 Newton’s equilibrium method, 514 P
Metals, 84–85 Newton’s interpolation formula, 369
Paez, T. L., 232
Method Nodal forces, 379, 410–414
Pagano, N. J., 51, 118
energy methods. see Energy Nodes, unevenly spaced, 567–568
Parallel-axis theorem, 692–694
methods Nonlinearly elastic, 91, 499,
Park, F. R., 466
equilibrium. see Equilibrium 508–509
Partial yielding, plastic, 612–613,
method Normal strain, 70
626–627
inverse and semi-inverse, 143 Normal stress
Particle mechanics, 514
numerical methods. see bending of beams and, 260–268
Pawlik, P. S., 51
Numerical methods on an oblique plane, 42–43
Peak compressive stress, 275
offset method, 86 shear stresses versus, 10
Pearson, K., 51, 354
of sections, 9–10, 270–276, 500 strain energy density for,
Perfectly plastic
superposition, 69–70, 282–284 108–109
beams, 590–600
Midsurface, shells, 658 tangential stresses combined
stresses in a flat disk, 612–614
Minimum potential energy, with, 296–300
stresses in thick-walled cylinder,
principle of, 515–517 transverse normal, 268–270
623–627
Mises, R. von, 199 transverse normal stress,
Permanent set/strain, 91
Mises-Hencky yield criterion, 268–270
Perry, D. J., 568
200–203, 610 Numerical accuracy, solving
Pierso, A. G., 232
Modes of failure, 192 problems, 680
Pilkey, W. D., 181, 711
Modes of force transmission, 14 Numerical methods
Pin-ended columns, buckling,
Modulus for an arbitrarily oriented bar
536–539, 541
bulk, 98 element, 380–383
Pivot points, 365
of elasticity, 90, 591 axial force equation, 384–386
Plane strain
of the foundation, 473 for a bar element, 379–380
elasticity, problems in, 135–137
of plasticity, 615 for a beam element, 393–399
generalized, 142
of rupture, 244 for boundary conditions,
Mohr’s circle for, 80–83
tangent modulus theory, 373–377
in a plate, 73–74
544–545 for case studies in plane stress,
transformation of, 23
of toughness/resilience, 109 414–422
two-dimensional, 70–72
Mohr’s circle central differences, 366–368
Plane stress
for curvature of plates, 638 computational tools, 423
bending of beams and, 257
for plane strain, 80–83 for curved boundaries, 370–373
in a cantilever beam, 414–418
three dimensional stress, 45–49 finite difference equations,
elasticity, problems in, 138–140
for torsional loading, 322 368–370
finite element applied to,
two-dimensional stress, 28–35 finite differences, 365
418–421
Mohr’s theory, failure criteria, finite element fundamentals,
on inclined planes, 321–324
206–207 377–378
introduction to, 13
Mollick, L., 466 finite element method formula,
plane-stress transformation,
Moment 402–407
23–26
bending, twisting, 14, 642 force-displacement relations,
in a thick walled cylinder, 436
force and moment over the ends 386–393
three-dimensional, 40
of a bar, 327–333 for a triangular finite element,
transformation of, 23–26
of inertia. see Inertia moments 407–414
Planes
plane stresses and, 638–639 for two-dimensional elements,
inclined planes, 321–324
plastic, 592 399–402
slip planes, 193
sign convention for, 243, 639
of symmetry, orthotropic
O materials, 101–103
N Oblique plane, stresses on, 42–45 Plastic/plasticity
Nadai, A., 232, 627 Octahedral shearing stress theory, axial deformation, 585–588
Navier, L., 2, 644 200–203 behavior of material, 578–579
Necking, 86 Octahedral stress, 44–45 collapse load of structures,
Neuber, H. P., 167, 181 Offset method, 86 600–605

Index 727
Plastic/plasticity (cont.) Prandtl’s stress function Rectangular bar, torsion of, 332–
deflection of beams, 588–590 boundary conditions of, 333, 338–340
deformation, 84, 579–580 326–327 Rectangular plates, stresses in,
elastic material versus, 91–92 force and moments over the 644–647, 650–652
elastic-plastic stresses in ends, 327–333 Rectangular rosette, 106–107
rotating disks, 612–614 Pressure Reddy, J. N., 524
elastic-plastic torsion, 605–609 contact, 443 Reismann, H., 51
hinge, 593–601, 604 cylinders under, 435–441 Repeated loading, 216
idealized stress-strain diagrams, Principal moments of inertia, Residual rotation, elastic-plastic,
580–582 694–700 608–609
instability in simple tension, Principal of superposition, 69. see Residual strain, 91
582–584 also Superposition Residual stress
perfectly plastic beams, Principal stresses elastic-plastic torsion and,
590–600 maximum in-plane shear stress 608–609
range, 88 and, 26–28 inelastically bent beam and, 592
stresses in perfectly plastic Mohr’s circle and, 29–30 perfectly plastic cylinders and,
cylinders, 623–627 solving cubic equations, 623
stress-strain relations, 614–623 682–686 plastic axial deformation and,
torsion: membrane analogy, in three dimensions, 38–42 585–588
610–611 Prismatic bars, torsion. see in rectangular beam, 598
types of materials, 85 Torsion, prismatic bars Revolution, shells of
Plate stress Problem formulation/solution, circular cylindrical, 664–667
axisymmetrically loaded 679–681 conical, 663
circular, 648–650 Progressive fracture, fatigue, spherical, 662–663
basic assumption of, 635–636 216–217 symmetrically loaded, 660–662
Plate stress (cont.) Proportional limit, 83 Right-hand screw rule, 14
boundary conditions of, Pure bending Rigid-plastic material, 580, 588,
642–644 615
of asymmetrical section of
with circular holes, 165–168 Robert Hooke, 92
beam, 246–251
curvature, and moment relations Roberts, A. D., 181
of inelastic beam, 588
of, 638–639 Rotating disks. see Disks, rotating
of perfectly plastic beams,
deflections of rectangular, Rotational inertia, 246
590–591
650–652 Rupture stress, 86
of plates, 638, 650–651
finite element solution for,
strain energy for beam in, 113
654–657
governing equations of plate of symmetrical section of beam, S
deflection, 640–641 242–246 Sachs, G., 627
sandwich, 652–654 Pure shear, 13, 92 Sadegh, A. M., 117, 181, 354, 711
simply supported rectangular, SAE relation, 219
644–647 Q Saint-Vanant’s principle, 115–117
strain-curvature relations of, Quinney, H., 232 Saint-Venant, semi-inverse method,
636–638 324–325
Poisson, S. D., 92–93 Sand hill analogy, 611
R Sandwich plates, stresses in,
Poisson’s ratio, 92–93, 579
Polar coordinates, basic relations Radial interference, 443 652–654
in, 152–157 Radial stress, 450 Schmidt, R. J., 117, 181, 300, 423
Polar representations, plane stress, Radians per second, rotating disks, Screw dislocation, 579
25 446 Secant formula, 552–554, 615
Polycrystalline structural metals, Radius Semi-infinite beams, elastic
failure, 193, 196 of curvature of beams, 245 foundations, 480–483
Polynomial solutions, 144–148 of curvature of shells, 658 Semi-inverse method, 143
Potential energy, 515–517 of gyration, 542, 692 Shaffer, B. W., 492
Power transmission, shafts in, Ramberg, W., 627 Shafts
323–324 Ramu, S. A., 492 stress concentration charts,
Prandtl’s membrane analogy Range stress, 218 708–710
equation of equilibrium, Rankine, W. J. M, 204 transmission of power by,
333–335 Ranov, T., 466 323–324
plastic torsion, 610–611 Rayleigh, Lord, 499, 522 Shames, I. H., 118
shearing stress and angle of Rayleigh-Ritz method, 522–524 Shape factor, 593
twist, 335–337 Reciprocity theorem, 498–499 Shear, form factor for, 285–286

728Index
Shear, pure, 13, 92 Spherical contacts, elasticity, material properties and, 68
Shear center, beams, 276–281 171–173 matrices, 401–402
Shear deformation, 109, 246 Spring constant, 227 measurement of, 103–107
Shear flow, 261, 279 Spring rate, 94 Mohr’s circle for plane strain,
Shear stress Springs 80–83
absolute maximum, 46 bar springs, 348–351 normal strain, 70
axial shear, 320–321 helical springs, 352–354 orthotropic materials, 101–103
bending of beams and, 260–268 Stable equilibrium, 515, 535 permanent and residual strain,
elementary theory of torsion Stainless steels, 85 91
and, 317 Statically indeterminate systems plane strain. see Plane strain
equality of, 12 applying energy methods to, rate, 226
fluid flow analogy and, 345 510–513 Saint-Vanant’s principle,
formula for, 261 axially loaded, 585–588 115–117
lateral, 347 bending of beams and, 281–284 state of strain at a point, 76
maximum in-plane, 26–28 deflections of beams, 706 stress-strain diagrams, 83–90,
normal stress versus, 10 infinite beams and, 474 580–582
oblique plane and, 42–45 statically determinate systems stress-strain relations, 614–623
Prandtl’s membrane analogy, versus, 9 superposition, 69–70
335–337 Steady-state creep, 195 thermal strain, 149–152
sign convention for, 640 Steels, 84–85 three-dimensional. see Three-
strain energy density for, Sternberg, E., 118 dimensional strain
108–109 Stiff beams, elastic foundations, transformation, 76–80
torsional, 335–337, 347 486–488 true and logarithmic strain, 88
transverse shear, 320–321 Stiffness matrix two-dimensional. see Two-
Shells, stress in the assembly process for, dimensional strain
cylindrical, of general shape, 386–387 Strain energy
670–673 for a beam element, 393–399 for axially loaded bars, 111–112
examples of revolution, for a cantilever beam, 415–416 in beams, 284–286
662–667 for deflection bar under for beams in bending, 113
simple membrane action of, combined loading, 405–406 Castigliano’s second theorem,
658–660 equilibrium method and the, 379 284–286
symmetrically loaded, 660–662 governing equations and the, of circular bars in torsion,
theories and behavior of, 383 112–113
657–658 nodal forces of a plate segment, in common structural members,
thermal stresses in compound 412–413 111
cylinders, 668–670 for a stepped axially loaded components of, 113–115
Shock loads, 225–226 bar, 388 deflections by stain energy
Short beams, 484 for thin plates, 657 method, 650–652
Short columns, 544, 554 for a three-bar truss, 390–391 deflections by trigonometric
Shrink-fit performance analysis, for triangular finite element, series, 518
disk flywheel, 451–452 410 for normal and shear stresses,
Shrinking allowance, 443 for a truss bar, 385 108–109
SI (International System of Units), Stiffness method, finite element for three-dimensional stresses,
2, 701–704 analysis, 377–378 110
Sign convention, stress, 11, 243 Strain variation in, 514–515
Slenderness ratio, 542 antiplane strain, 142–143 Strain-displacement relations
Slip action, 579 definition of, 70 bending of beams and, 244–246
Slip planes, failure by yielding, effective strain, 615 in large (finite) strains, 74
193 elastic versus plastic behavior, in plate bending, 636–638
Slope, 257 91–92 in polar coordinates, 154–155
Small deflection theory energy. see Strain energy in three-dimensional problems,
of bending, 636–637 equations of compatibility, 72, 78–80
of sandwich plates, 653 75–76 in two-dimensional problems,
of shells, 658 Eulerian and Lagrangian 76–78
Sokolnikoff, I. S., 51, 117, 181, coordinates, 73–74 Strength coefficient, 581
300, 524 hardening, 92, 194, 580–581 Stress
Solderberg criterion, 219, 222 Hooke’s law and Poisson’s ratio, bending of beams and, 260–268
Solid disks, rotating, 448–449, 92–100 boundary conditions of surface
459 large strain, 74 forces, 49–50
Solving problems, 679–681 LST and CST elements, 407 combined stress, 17

Index 729
Stress (cont.) residual stress. see Residual stress normal stress combined with,
concentration. see Stress resultants, 13–17, 24, 33 296–300
concentration in rotating disks, 612–614 thick walled cylinders under,
conditions of equilibrium, 8–9 shear stress. see Shear stress 441–442
contact stresses, 169–171 shells and. see Shells, stress in Tangential stress in cylinders,
critical stress in columns, sign convention, 11, 243 441–442
543–544 strain energy for normal and Taylor, G. L., 232
cubic equation, 39–40 shear stresses, 108–109 Taylor, R. I., 423, 466, 673
definitions and components, strain energy for three- Temperature
9–13 dimensional stresses, 110 ductile-brittle transition and,
dilatational stress, 113 stress-strain diagrams, 83–90, 230–231
distortional (deviator) stress, 580–582 stress-strain properties and
114, 579 stress-strain relations, 614–623 effects of, 90
distribution. see Stress tangential stress. see Tangential Tension, plastic, 582–584
distribution stress Theories
effective, 199, 615 tensor, 12 and behavior of shells, 657–658
elastic-plastic stresses in thermal stress. see Thermal stress elementary theory of torsion,
rotating disks, 612–614 three-dimensional 316–321
function, Airy’s, 141–143 transformation, 35–38 small deflection theory for
function, Prandtl’s. see Prandtl’s trajectories, 25 plates, 636
stress function in transformed composite theory of elasticity, 1–2, 133
historical development of stress beams, 272 Theories, bending of beams
analysis, 2–3 triaxial stress. see Triaxial stress comparison of, 293–296
idealized stress-strain diagrams, uniaxial stress. see Uniaxial stress elementary bending theory, 246,
580–582 variation within a body, 20–22 253–254, 256
on inclined planes, 17–19, Stress concentration Timoshenko beam theory, 246
321–324 for disk flywheel, 450 Theories, columns
indicial notation, 50 elasticity problems, 157–169 Johnson’s buckling criterion,
intensity factors, 211–213 fluid flow analogy and, 344–346 546–548
matrices, 401–402 Stress cubic equation, 682–686 tangent modulus theory,
maximum in-plane shear, 26–28 Stress distribution 544–546
maximum shearing stress charts of bars/shafts and, 707–710 Theories, failure
theory, 198–199 concentrated load, acting on a comparison of yielding, 204
mean stress, 218–219 beam, 161–163 Coulomb-Mohr, 207–210
Stress (cont.) contact stress distribution, maximum distortion energy,
mechanics of materials/theory 174–178 199–200, 625
of elasticity, 1–2 elasticity theory and, 288–289 maximum principal stress,
methods for design and, 4–8 pure bending of beams, 204–206
Mohr’s circle for two- asymmetrical, 248 maximum shear stress, 198, 625
dimensional, 28–35 Stress-strain diagrams, 83–90, Mohr’s, 206–207
Mohr’s circle in three 580–582 octahedral shearing stress,
dimensions, 45–49 Stress-strain relations, 139–140 200–203
normal stress. see Normal stress Sullivan, J. L., 232 for thick walled cylinders,
on an oblique plane, 42–45 Sundara, R., 492 442–443
octahedral shearing stress, Superposition, 69–70, 282–284, Thermal stress
200–203 450 in compound cylinders,
perfectly plastic thick-walled Surface forces, 2, 49–50 668–670
cylinders, 623–627 Symmetrical cross section, beams, elasticity, problems in, 149–152
plane stress. see Plane stress 242–246 in long circular cylinders,
plates and. see Plate stress 460–464
Symmetrically loaded shells,
principal stress theory, 204–206 in thin rotating disks, 458–459
660–662
principal stresses and cubic Thermoelasticity, 149–152
Systems of units, 2
equations, 682–686 Thermoplastic, 85
principal stresses in three Thermosets, engineering materials,
dimensions, 38–42 T 85
problems in stress analysis, Tangent modulus theory, columns, Thick walled cylinders. see
315–314 544–546 Cylinders, thick walled
pure shear, 13 Tangential stress Thin disks, thermal stresses,
range stress, 218 curved beams and, 291–293 458–459
relaxation, stress, 195 disk flywheel under, 450–452 Thin plates, bending of, 636

730Index
Thin-walled (plate and shell) strain energy for axially loaded V
structures, 8 bars, 111–112
Van Vlack, L. H., 232
Thin-walled cylinders, 435–436 of stresses on inclined planes,
Variational method, 496
Thin-walled sections 321–324
restrained, in torsion, 346–347 Virtual displacement, 514
of thin-walled sections, 340–
shear center in, 277–281 344, 346–349 Virtual work, 514–515, 518–522
torsion of multiply connected, of a three-cell tube, 343–344 Viscoelastic material, 83–84
340–344 Total strain theory, 614 Volume change, 93, 98–100
Thin-walled vessels, 15–16, 32–33 Transfer formula, 692–694 Von Mises stress, 615
Three-dimensional strain Transformation equations, 155– Von Mises theory, 199
for orthotropic materials, 156, 248
139–140 Transformed section method, W
overview of, 72 composite beams, 270–276 Wahl, A. M., 354
transformation of, 78–80 Transition, ductile-brittle, 230–232 Warping deformation, 315
Three-dimensional stress Transverse stress Wedge, concentrated loads,
invariants for, 40–42 normal, 268–270 157–160
Mohr’s circle and, 45–46 shear, 320–321 Weight
principal stresses in, 38–40 Tresca yield criterion, 198, 612 freely falling, 227–228
strain energy density for, 110 Triangular element, numerical horizontally moving,
transformation of, 35–38 method for, 407–414 228–230
Timoshenko, S. P., 3, 51, 118, Triaxial stress, 12–13, 42, 45, 98, Weiner, J. H., 181
181, 232, 300, 354, 524, 200, 231 Wheatstone bridge, 104
568, 673 Triaxiality, ductile-brittle Winkler model, 473
Timoshenko beam theory, 246 transition, 231 Winkler’s Formula, 293
Ting, B. Y., 492 Trigonometric functions, 680 Winkler’s theory, pure bending,
Titanic case study, 231–232 Trigonometric series, 483, 517–522 289
Todhunter, I., 51, 354 True factor of safety, 6–7 Wire gage, 103–104
Torque diagram, 317 True stress-strain, 88, 580–583, Woinowsky-Krieger, S., 673
Torque-displacement relations, 320 616, 623 Work hardening material, 580
Torsion Turbine application, 454 Work-strain energy, 497
charts of bars/shafts in, 707–709 Twisting moments, 14 Wriggers, P., 181
elementary theory of, 316–321 Two-dimensional elements,
formula, 112, 317 numerical method, 399–402
plastic behavior of materials, Two-dimensional strain
Y
605–611 elasticity problems, 134–135, Yang, T.Y., 423, 466, 673
torsional rigidity, 317, 330 140–141 Yield/yielding
torsional shear, 347 matrix, 401–402 comparison of theories, 204
Torsion, prismatic bars transformation of, 76–78 complete yielding in plastic,
bar springs, 350–351 614, 624–626
of bar with square cross-section, criteria of, 197–198
368–370
U
initial yielding in plastic, 612
basic formulas for stress, 15–16 Ugural, A. C., 51, 117, 181, 232, Mises-Hencky yield criterion,
of circular bars, 316–321 300, 354, 423, 466, 492, 524, 200–203, 610
of curved circular bars, 351–354 569, 627, 673, 681, 711 moment, 593
of elliptical bar, 328–331, Ultimate partial yielding in plastic,
372–373 load, 601 612–613, 626–627
of equilateral triangle bar, moment, 592 slip plane failure, 193–195
331–332 tensile stress, 86 surfaces for triaxial stress, 200
fluid flow analogy, 344–346 torque, 607 torque, 606
general solution to, 324–326 Uniaxial stress, 13, 92–93, Tresca yield criterion, 612
introduction to, 315–316 218–221 Young, Thomas, 2, 92
of narrow rectangular section, Uniaxial tension, 108 Young, W. C., 181, 300, 354, 700,
338–340 Uniform stress, rotating disks, 711
Prandtl’s membrane analogy, 456–458 Young’s modulus, 90, 92
333–337 Unit load, energy methods,
Prandtl’s stress function, 506–508
326–333 Units, System of, 701–704 Z
of a rectangular bar, 332–333 Unstable equilibrium, 535 Zhang, S., 673
of a rectangular tube, 342–343 U.S. Customary System of Units Zhao, D., 673
of a stepped bar, 336–337 (USCS), 701–704 Zienkiewicz, O. C., 423, 466, 673

Index 731