Stiffness and Deflection Analysis of Complex Structures

Total Page:16

File Type:pdf, Size:1020Kb

Load more

JOURNAL OF THE AERONAUTICAL SCIENCES VOLUME 23 SEPTEMBER, 1956 NUMBER 9 Stiffness and Deflection Analysis of Complex Structures M. J. TURNER,* R. W. CLOUGH,t H. C. MARTIN,t AND L. J. TOPP** ABSTRACT tion on static air loads, and theoretical analysis of aero- A method is developed for calculating stiffness influence co- elastic effects on stability and control. This is a prob- efficients of complex shell-type structures. The object is to pro- lem of exceptional difficulty when thin wings and tail vide a method that will jdeld structural data of sufficient accuracy surfaces of low aspect ratio, either swept or unswept, to be adequate for subsequent dynamic and aeroelastic analyses. are involved. Stiffness of the complete structure is obtained by summing It is recognized that camber bending (or rib bending) stiffnesses of individual units. Stiffnesses of typical structural components are derived in the paper. Basic conditions of con- is a significant feature of the vibration modes of the tinuity and equilibrium are established at selected points (nodes) newer configurations, even of the low-order modes; in the structure. Increasing the number of nodes increases the in order to encompass these characteristics it seems accuracy of results. Any physically possible support conditions likely that the load-deflection relations of a practical can be taken into account. Details in setting up the analysis can structure must be expressed in the form of either de- be performed by nonengineering trained personnel; calculations are conveniently carried out on automatic digital computing flection or stiffness influence coefficients. One ap- equipment. proach is to employ structural models and to determine Method is illustrated by application to a simple truss, a flat the influence coefficients experimentally; it is antici- plate, and a box beam. Due to shear lag and spar web deflection, pated that the experimental method will be employed the box beam has a 25 per cent greater deflection than predicted extensively in the future, either in lieu of or as a final from beam theory. It is shown that the proposed method cor- rectly accounts for these effects. check on the result of analysis. However, elaborate Considerable extension of the material presented in the paper models are expensive, they take a long time to build, is possible. and tend to become obsolete because of design changes ; for these reasons it is considered essential that a con- (I) INTRODUCTION tinuing research effort should be applied to the devel- RESENT CONFIGURATION TRENDS in the design of opment of analytical methods. It is to be expected that modern developments in high-speed digital com- Downloaded by VIRGINIA TECHNICAL UNIVERSITY on September 12, 2014 | http://arc.aiaa.org DOI: 10.2514/8.3664 Phigh-speed aircraft have created a number of difficult, fundamental structural problems for the puting machines will make possible a more fundamental worker in aeroelasticity and structural dynamics. The approach to the problems of structural analysis; we chief problem in this category is to predict, for a given shall expect to base our analysis on a more realistic elastic structure, a comprehensive set of load-deflection and detailed conceptual model of the real structure relations which can serve as structural basis for dynamic than has been used in the past. As indicated by the load calculations, theoretical vibration and flutter title, the present paper is exclusively concerned with analyses, estimation of the effects of structural deflec- methods of theoretical analysis; also it is our object to outline the development of a method that is well Received June 29, 1955. This paper is based on a paper adapted to the use of high-speed digital computing presented at the Aeroelasticity Session, Twenty-Second Annual machinery, Meeting, IAS, New York, January 25-29, 1954. * Structural Dynamics Unit Chief, Boeing Airplane Company, Seattle Division. (II) REVIEW OF EXISTING METHODS OF STRUCTURAL f Associate Professor of Civil Engineering, University of Cali- ANALYSIS fornia, Berkeley. J Professor of Aeronautical Engineering, University of Wash- (1) Elementary Theories of Flexure and Torsion ington, Seattle. ** Structures Engineer, Structural Dynamics Unit, Boeing Air- The limitations of these venerable theories are too plane Company, Wichita Division. well known to justify extensive comment. They are 805 806 JOURNAL OF THE AERONAUTICAL SCIENCES — S E P T E M B E R , 1956 lQ n adequate only for low-order modes of elongated struc- (5) Direct Stiffness Calculation: Levy, Schuerch > tures. When the loading is complex (as in the case In a recent paper Levy has presented a method of of inertia loading associated with a mode of high order) analysis for highly redundant structures which is par- refinements are required to account for secondary ticularly suited to the use of high-speed digital com- effects such as shear lag and torsion-bending. puting machines. The structure is regarded as an assemblage of beams (ribs and spars) and interspar (2) Wide Beam Theory: Schuerch1 torque cells. The stiffness matrix for the entire struc- ture is computed by simple summation of the stiff- Schuerch has devised a generalized theory of com- ness matrices of the elements of the structure. Fi- bined flexure and torsion which is applicable to multi- nally, the matrix of deflection influence coefficients is spar wide beams having essentially rigid ribs. Torsion- obtained by inversion of the stiffness matrix. Schuerch bending effects are included but not shear lag. It is has also presented a discussion of the problem from the expected that wide beam theory will be used extensively point of view of determining the stiffness coefficients. in the solution of static aeroelastic problems (effect of air-frame flexibility on steady air loads, stability, etc.). However, the rigid rib assumption appears to limit its (Ill) S O M E U N S O L V E D P R O B L E M S utility rather severely for vibration and flutter anal- At the present time, it is believed that the greatest ysis of thin low aspect ratio wings. need is to derive a numerical method of analysis for a class of structures intermediate between the thin (3) Method of Redundant Forces: Levy, Bisplinghoff and 2 stiffened shell and the solid plate. These are hollow Lang, Langefors, Rand, Wehle and Lansing ^ structures having a rather large share of the bending These writers have contributed the basic papers material located in the skin, which is relatively thick leading to the present widespread use of energy prin- but still thin enough so that we may safely neglect ciples, matrix algebra, and influence coefficients in the its plate bending stiffness. In order to cope with this solution of structural deflection problems. Redundant class of structures successfully, we must base our internal loads are determined by the principle of least analysis upon a structural idealization that is suffi- work, and deflections are obtained by application of ciently realistic to encompass a fairly general two- Castigliano's theorem. The method is, of course, dimensional stress distribution in the cover plates; perfectly general. However, the computational diffi- and our method of analysis must yield the load-deflec- culties become severe if the structure is highly re- tion relations associated with such stresses. It is char- dundant, and the method is not particularly well acteristic of these problems that the directions of prin- adapted to the use of high-speed computing machines. cipal stresses in certain critical parts of the structure Rand has suggested a method of solution for stresses cannot be determined by inspection. Hence, the in highly redundant structures which might also be familiar methods of structural analysis based upon the used for calculating deflections. Instead of using concepts of axial load carrying members, joined by member loads as redundants, he proposes to employ membranes carrying pure shear, are not satisfactory, systems of self-equilibrating internal stresses. These even if we employ effective width concepts to account redundant stresses may be regarded as perturbations for the bending resistance of the skin. We should like of a primary stress distribution that is in equilibrium to include shear lag, torsion-bending, and Poisson's with the external loads (but does not generally satisfy ratio effects to a sufficient approximation for reliable compatibility conditions). The number of properly prediction of vibration modes and natural frequencies chosen redundants required to obtain a satisfactory of moderate order. Also, we should like to avoid any solution may be considerably less than the "degree of assumptions of closely spaced rigid diaphragms or of Downloaded by VIRGINIA TECHNICAL UNIVERSITY on September 12, 2014 | http://arc.aiaa.org DOI: 10.2514/8.3664 redundancy." Successful application of this method orthotropic cover plates, which have been introduced requires a high degree of engineering judgment, and in many papers on advanced structural analysis. The the accuracy of the results is very difficult to evaluate. actual rib spacing and finite rib stiffnesses should be accounted for in a realistic fashion. In summary, what (4) Plate Methods: Fung, Reissner, Bens cot er, and is required is an approximate numerical method of MacNeaV* analysis which avoids drastic modification of the As the trend toward thinner sections approaches the geometry of the structure or artificial constraints of its ultimate limit, we enter first a regime of very thick elastic
Recommended publications
  • Phys 140A Lecture 8–Elastic Strain, Compliance, and Stiffness

    Phys 140A Lecture 8–Elastic Strain, Compliance, and Stiffness

    Lecture 8 Elastic strains, compliance, and stiffness Review for exam Stress: force applied to a unit area Strain: deformation resulting from stress Purpose of today’s derivations: Generalize Hooke’s law to a 3D body that may be subjected to any arbitrary force. We imagine 3 orthogonal vectors 풙̂, 풚̂, 풛̂ embedded in a solid before we have deformed it. After we have deformed the solid, these vectors might be of different length and they might be pointing in different directions. We describe these deformed vectors 풙′, 풚′, 풛′ in the following way: ′ 풙 = (1 + 휖푥푥)풙̂ + 휖푥푦풚̂ + 휖푥푧풛̂ ′ 풚 = 휖푦푥풙̂ + (1 + 휖푦푦)풚̂ + 휖푦푧풛̂ ′ 풛 = 휖푧푥풙̂ + 휖푧푦풚̂ + (1 + 휖푧푧)풛̂ The components 휖훼훽 define the deformation. They are dimensionless and have values much smaller than 1 in most instances in solid state physics. How one ‘reads’ these vectors is (for example) 휖푥푥force along x, deformation along x; 휖푥푦shear xy-plane along x or y direction (depending which unit vector it is next to) The new axes have new lengths given by (for example): ′ ′ 2 2 2 풙 ∙ 풙 = 1 + 2휖푥푥 + 휖푥푥 + 휖푥푦 + 휖푥푧 Usually we are dealing with tiny deformations so 2nd order terms are often dropped. Deformation also changes the volume of a solid. Considering our unit cube, originally it had a volume of 1. After distortion, it has volume 1 + 휖푥푥 휖푥푦 휖푥푧 ′ ′ ′ ′ 푉 = 풙 ∙ 풚 × 풛 = | 휖푦푥 1 + 휖푦푦 휖푦푧 | ≈ 1 + 휖푥푥 + 휖푦푦 + 휖푧푧 휖푧푥 휖푧푦 1 + 휖푧푧 Dilation (훿) is given by: 푉 − 푉′ 훿 ≡ ≅ 휖 + 휖 + 휖 푉 푥푥 푦푦 푧푧 (2nd order terms have been dropped because we are working in a regime of small distortion which is almost always the appropriate one for solid state physics.
  • 19.1 Attitude Determination and Control Systems Scott R. Starin

    19.1 Attitude Determination and Control Systems Scott R. Starin

    19.1 Attitude Determination and Control Systems Scott R. Starin, NASA Goddard Space Flight Center John Eterno, Southwest Research Institute In the year 1900, Galveston, Texas, was a bustling direct hit as Ike came ashore. Almost 200 people in the community of approximately 40,000 people. The Caribbean and the United States lost their lives; a former capital of the Republic of Texas remained a tragedy to be sure, but far less deadly than the 1900 trade center for the state and was one of the largest storm. This time, people were prepared, having cotton ports in the United States. On September 8 of received excellent warning from the GOES satellite that year, however, a powerful hurricane struck network. The Geostationary Operational Environmental Galveston island, tearing the Weather Bureau wind Satellites have been a continuous monitor of the gauge away as the winds exceeded 100 mph and world’s weather since 1975, and they have since been bringing a storm surge that flooded the entire city. The joined by other Earth-observing satellites. This weather worst natural disaster in United States’ history—even surveillance to which so many now owe their lives is today—the hurricane caused the deaths of between possible in part because of the ability to point 6000 and 8000 people. Critical in the events that led to accurately and steadily at the Earth below. The such a terrible loss of life was the lack of precise importance of accurately pointing spacecraft to our knowledge of the strength of the storm before it hit. daily lives is pervasive, yet somehow escapes the notice of most people.
  • Beam Element Stiffness Matrices

    Beam Element Stiffness Matrices

    Beam Element Stiffness Matrices CEE 421L. Matrix Structural Analysis Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall, 2014 Truss elements carry only axial forces. Beam elements carry shear forces and bending moments. Frame elements carry shear forces, bending moments, and axial forces. This document presents the development of beam element stiffness matrices in local coordinates. 1 A simply supported beam carrying end-moments Consider a simply supported beam resisting moments M1 and M2 applied at its ends. The flexibility relates the end rotations {θ1, θ2} to the end moments {M1,M2}: θ1 F11 F12 M1 = . θ2 F21 F22 M2 The flexibility coefficients, Fij, may be obtained from Castigliano’s 2nd Theo- ∗ rem, θi = ∂U (Mi)/∂Mi. First Column Second Column 2 CEE 421L. Matrix Structural Analysis – Duke University – Fall 2014 – H.P. Gavin The applied moments M1 and M2 are in equilibrium with the reactions forces V1 and V2; V1 = (M1 + M2)/L and V2 = −(M1 + M2)/L M + M x ! x V (x) = 1 2 M(x) = M − 1 + M L 1 L 2 L The total potential energy of a beam with these forces and moments is: 1Z L M 2 1Z L V 2 U = dx + dx 2 0 EI 2 0 G(A/α) By Castigliano’s Theorem, ∂U θ1 = ∂M1 ∂M(x) ∂V (x) Z L M(x) Z L V (x) = ∂M1 dx + ∂M1 dx 0 EI 0 G(A/α) x 2 x x Z L Z L Z L Z L L − 1 dx αdx L L − 1 dx αdx = + M1 + + M2 0 EI 0 GAL2 0 EI 0 GAL2 and ∂U θ2 = ∂M2 ∂M(x) ∂V (x) Z L M(x) Z L V (x) = ∂M2 dx + ∂M2 dx 0 EI 0 G(A/α) x x x 2 Z L Z L Z L Z L L L − 1 dx αdx L dx αdx = + M1 + + M2 0 EI 0 GAL2 0 EI 0 GAL2
  • Introduction to Stiffness Analysis Stiffness Analysis Procedure

    Introduction to Stiffness Analysis Stiffness Analysis Procedure

    Introduction to Stiffness Analysis displacements. The number of unknowns in the stiffness method of The stiffness method of analysis is analysis is known as the degree of the basis of all commercial kinematic indeterminacy, which structural analysis programs. refers to the number of node/joint Focus of this chapter will be displacements that are unknown development of stiffness equations and are needed to describe the that only take into account bending displaced shape of the structure. deformations, i.e., ignore axial One major advantage of the member, a.k.a. slope-deflection stiffness method of analysis is that method. the kinematic degrees of freedom In the stiffness method of analysis, are well-defined. we write equilibrium equations in 1 2 terms of unknown joint (node) Definitions and Terminology Stiffness Analysis Procedure Positive Sign Convention: Counterclockwise moments and The steps to be followed in rotations along with transverse forces and displacements in the performing a stiffness analysis can positive y-axis direction. be summarized as: 1. Determine the needed displace- Fixed-End Forces: Forces at the “fixed” supports of the kinema- ment unknowns at the nodes/ tically determinate structure. joints and label them d1, d2, …, d in sequence where n = the Member-End Forces: Calculated n forces at the end of each element/ number of displacement member resulting from the unknowns or degrees of applied loading and deformation freedom. of the structure. 3 4 1 2. Modify the structure such that it fixed-end forces are vectorially is kinematically determinate or added at the nodes/joints to restrained, i.e., the identified produce the equivalent fixed-end displacements in step 1 all structure forces, which are equal zero.
  • Force/Deflection Relationships

    Force/Deflection Relationships

    Introduction to Finite Element Dynamics ⎡ x x ⎤ 1. Element Stiffness Matrices u()x = ⎢1− ⎥ u = []C u (3) ⎣ L L ⎦ 1.1 Bar Element (iii) Derive Strain - Displacement Relationship Consider a bar element shown below. At the two by using Mechanics Theory2 ends of the bar, axially aligned forces [F1, F2] are The axial strain ε(x) is given by the following applied producing deflections [u ,u ]. What is the 1 2 relationship between applied force and deflection? du d ⎡ 1 1 ⎤ ε()x = = ()[]C u = []B u = ⎢− ⎥u (4) dx dx ⎣ L L ⎦ Deformed shape u2 u1 (iv) Derive Stress - Displacement Relationship by using Elasticity Theory The axial stresses σ (x) are given by the following assuming linear elasticity and homogeneity of F1 x F2 material throughout the bar element. Where E is Element Node the material modulus of elasticity. There are typically five key stages in the σ (x)(= Eε x)= E[]B u (5) 1 analysis . (v) Use principle of Virtual Work (i) Conjecture a displacement function There is equilibrium between WI , the internal The displacement function is usually an work done in deforming the bar, and WE , approximation, that is continuous and external work done by the movement of the differentiable, most usually a polynomial. u(x) is applied forces. The bar cross sectional area the axial deflections at any point x along the bar A = ∫∫dydz is assumed constant with respect to x. element. WI = ∫∫∫σ (x).ε (x)(dxdydz = ∫σ x).ε (x)dx∫∫dydz ⎡a1 ⎤ (6) u()x = a1 + a2 x = []1 x ⎢ ⎥ = [N] a (1) = A∫σ ()x .ε ()x dx ⎣a2 ⎦ ⎡F ⎤ W = []u u 1 = uT F (7) E 1 2 ⎢F ⎥ (ii) Express u()x in terms of nodal ⎣ 2 ⎦ displacements by using boundary conditions.
  • Deflection of Beams Introduction

    Deflection of Beams Introduction

    Deflection of Beams Introduction: In all practical engineering applications, when we use the different components, normally we have to operate them within the certain limits i.e. the constraints are placed on the performance and behavior of the components. For instance we say that the particular component is supposed to operate within this value of stress and the deflection of the component should not exceed beyond a particular value. In some problems the maximum stress however, may not be a strict or severe condition but there may be the deflection which is the more rigid condition under operation. It is obvious therefore to study the methods by which we can predict the deflection of members under lateral loads or transverse loads, since it is this form of loading which will generally produce the greatest deflection of beams. Assumption: The following assumptions are undertaken in order to derive a differential equation of elastic curve for the loaded beam 1. Stress is proportional to strain i.e. hooks law applies. Thus, the equation is valid only for beams that are not stressed beyond the elastic limit. 2. The curvature is always small. 3. Any deflection resulting from the shear deformation of the material or shear stresses is neglected. It can be shown that the deflections due to shear deformations are usually small and hence can be ignored. Consider a beam AB which is initially straight and horizontal when unloaded. If under the action of loads the beam deflect to a position A'B' under load or infact we say that the axis of the beam bends to a shape A'B'.
  • Mechanics of Materials Chapter 6 Deflection of Beams

    Mechanics of Materials Chapter 6 Deflection of Beams

    Mechanics of Materials Chapter 6 Deflection of Beams 6.1 Introduction Because the design of beams is frequently governed by rigidity rather than strength. For example, building codes specify limits on deflections as well as stresses. Excessive deflection of a beam not only is visually disturbing but also may cause damage to other parts of the building. For this reason, building codes limit the maximum deflection of a beam to about 1/360 th of its spans. A number of analytical methods are available for determining the deflections of beams. Their common basis is the differential equation that relates the deflection to the bending moment. The solution of this equation is complicated because the bending moment is usually a discontinuous function, so that the equations must be integrated in a piecewise fashion. Consider two such methods in this text: Method of double integration The primary advantage of the double- integration method is that it produces the equation for the deflection everywhere along the beams. Moment-area method The moment- area method is a semigraphical procedure that utilizes the properties of the area under the bending moment diagram. It is the quickest way to compute the deflection at a specific location if the bending moment diagram has a simple shape. The method of superposition, in which the applied loading is represented as a series of simple loads for which deflection formulas are available. Then the desired deflection is computed by adding the contributions of the component loads (principle of superposition). 6.2 Double- Integration Method Figure 6.1 (a) illustrates the bending deformation of a beam, the displacements and slopes are very small if the stresses are below the elastic limit.
  • Ch. 8 Deflections Due to Bending

    Ch. 8 Deflections Due to Bending

    446.201A (Solid Mechanics) Professor Youn, Byeng Dong CH. 8 DEFLECTIONS DUE TO BENDING Ch. 8 Deflections due to bending 1 / 27 446.201A (Solid Mechanics) Professor Youn, Byeng Dong 8.1 Introduction i) We consider the deflections of slender members which transmit bending moments. ii) We shall treat statically indeterminate beams which require simultaneous consideration of all three of the steps (2.1) iii) We study mechanisms of plastic collapse for statically indeterminate beams. iv) The calculation of the deflections is very important way to analyze statically indeterminate beams and confirm whether the deflections exceed the maximum allowance or not. 8.2 The Moment – Curvature Relation ▶ From Ch.7 à When a symmetrical, linearly elastic beam element is subjected to pure bending, as shown in Fig. 8.1, the curvature of the neutral axis is related to the applied bending moment by the equation. ∆ = = = = (8.1) ∆→ ∆ For simplification, → ▶ Simplification i) When is not a constant, the effect on the overall deflection by the shear force can be ignored. ii) Assume that although M is not a constant the expressions defined from pure bending can be applied. Ch. 8 Deflections due to bending 2 / 27 446.201A (Solid Mechanics) Professor Youn, Byeng Dong ▶ Differential equations between the curvature and the deflection 1▷ The case of the large deflection The slope of the neutral axis in Fig. 8.2 (a) is = Next, differentiation with respect to arc length s gives = ( ) ∴ = → = (a) From Fig. 8.2 (b) () = () + () → = 1 + Ch. 8 Deflections due to bending 3 / 27 446.201A (Solid Mechanics) Professor Youn, Byeng Dong → = (b) (/) & = = (c) [(/)]/ If substitutng (b) and (c) into the (a), / = = = (8.2) [(/)]/ [()]/ ∴ = = [()]/ When the slope angle shown in Fig.
  • Fundamentals of Biomechanics Duane Knudson

    Fundamentals of Biomechanics Duane Knudson

    Fundamentals of Biomechanics Duane Knudson Fundamentals of Biomechanics Second Edition Duane Knudson Department of Kinesiology California State University at Chico First & Normal Street Chico, CA 95929-0330 USA [email protected] Library of Congress Control Number: 2007925371 ISBN 978-0-387-49311-4 e-ISBN 978-0-387-49312-1 Printed on acid-free paper. © 2007 Springer Science+Business Media, LLC All rights reserved. This work may not be translated or copied in whole or in part without the written permission of the publisher (Springer Science+Business Media, LLC, 233 Spring Street, New York, NY 10013, USA), except for brief excerpts in connection with reviews or scholarly analysis. Use in connection with any form of information storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter developed is forbidden. The use in this publication of trade names, trademarks, service marks and similar terms, even if they are not identified as such, is not to be taken as an expression of opinion as to whether or not they are subject to proprietary rights. 987654321 springer.com Contents Preface ix NINE FUNDAMENTALS OF BIOMECHANICS 29 Principles and Laws 29 Acknowledgments xi Nine Principles for Application of Biomechanics 30 QUALITATIVE ANALYSIS 35 PART I SUMMARY 36 INTRODUCTION REVIEW QUESTIONS 36 CHAPTER 1 KEY TERMS 37 INTRODUCTION TO BIOMECHANICS SUGGESTED READING 37 OF UMAN OVEMENT H M WEB LINKS 37 WHAT IS BIOMECHANICS?3 PART II WHY STUDY BIOMECHANICS?5 BIOLOGICAL/STRUCTURAL BASES
  • Chapter 10: Elasticity and Oscillations

    Chapter 10: Elasticity and Oscillations

    Chapter 10 Lecture Outline 1 Copyright © The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Chapter 10: Elasticity and Oscillations •Elastic Deformations •Hooke’s Law •Stress and Strain •Shear Deformations •Volume Deformations •Simple Harmonic Motion •The Pendulum •Damped Oscillations, Forced Oscillations, and Resonance 2 §10.1 Elastic Deformation of Solids A deformation is the change in size or shape of an object. An elastic object is one that returns to its original size and shape after contact forces have been removed. If the forces acting on the object are too large, the object can be permanently distorted. 3 §10.2 Hooke’s Law F F Apply a force to both ends of a long wire. These forces will stretch the wire from length L to L+L. 4 Define: L The fractional strain L change in length F Force per unit cross- stress A sectional area 5 Hooke’s Law (Fx) can be written in terms of stress and strain (stress strain). F L Y A L YA The spring constant k is now k L Y is called Young’s modulus and is a measure of an object’s stiffness. Hooke’s Law holds for an object to a point called the proportional limit. 6 Example (text problem 10.1): A steel beam is placed vertically in the basement of a building to keep the floor above from sagging. The load on the beam is 5.8104 N and the length of the beam is 2.5 m, and the cross-sectional area of the beam is 7.5103 m2.
  • Glossary of Materials Engineering Terminology

    Glossary of Materials Engineering Terminology

    Glossary of Materials Engineering Terminology Adapted from: Callister, W. D.; Rethwisch, D. G. Materials Science and Engineering: An Introduction, 8th ed.; John Wiley & Sons, Inc.: Hoboken, NJ, 2010. McCrum, N. G.; Buckley, C. P.; Bucknall, C. B. Principles of Polymer Engineering, 2nd ed.; Oxford University Press: New York, NY, 1997. Brittle fracture: fracture that occurs by rapid crack formation and propagation through the material, without any appreciable deformation prior to failure. Crazing: a common response of plastics to an applied load, typically involving the formation of an opaque banded region within transparent plastic; at the microscale, the craze region is a collection of nanoscale, stress-induced voids and load-bearing fibrils within the material’s structure; craze regions commonly occur at or near a propagating crack in the material. Ductile fracture: a mode of material failure that is accompanied by extensive permanent deformation of the material. Ductility: a measure of a material’s ability to undergo appreciable permanent deformation before fracture; ductile materials (including many metals and plastics) typically display a greater amount of strain or total elongation before fracture compared to non-ductile materials (such as most ceramics). Elastic modulus: a measure of a material’s stiffness; quantified as a ratio of stress to strain prior to the yield point and reported in units of Pascals (Pa); for a material deformed in tension, this is referred to as a Young’s modulus. Engineering strain: the change in gauge length of a specimen in the direction of the applied load divided by its original gauge length; strain is typically unit-less and frequently reported as a percentage.
  • On the Relationship Between Indenation Hardness and Modulus, and the Damage Resistance of Biological Materials

    On the Relationship Between Indenation Hardness and Modulus, and the Damage Resistance of Biological Materials

    bioRxiv preprint doi: https://doi.org/10.1101/107284; this version posted February 9, 2017. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. On the relationship between indenation hardness and modulus, and the damage resistance of biological materials. David Labontea,∗, Anne-Kristin Lenzb, Michelle L. Oyena aThe NanoScience Centre, Department of Engineering, Cambridge, UK bUniversity of Applied Sciences, Bremen, Germany Abstract The remarkable mechanical performance of biological materials is based on intricate structure-function relation- ships. Nanoindentation has become the primary tool for characterising biological materials, as it allows to relate structural changes to variations in mechanical properties on small scales. However, the respective theoretical back- ground and associated interpretation of the parameters measured via indentation derives largely from research on ‘traditional’ engineering materials such as metals or ceramics. Here, we discuss the functional relevance of inden- tation hardness in biological materials by presenting a meta-analysis of its relationship with indentation modulus. Across seven orders of magnitude, indentation hardness was directly proportional to indentation modulus, illus- trating that hardness is not an independent material property. Using a lumped parameter model to deconvolute indentation hardness into components arising from reversible and irreversible deformation, we establish crite- ria which allow to interpret differences in indentation hardness across or within biological materials. The ratio between hardness and modulus arises as a key parameter, which is a proxy for the ratio between irreversible and reversible deformation during indentation, and the material’s yield strength.