|
|
| Line 1: |
Line 1: |
| [[File:Hookes-law-springs.png|thumb|Hooke's law: the force is proportional to the extension]]
| | Myrtle Benny is how I'm known as and I feel comfy when people use the full name. Bookkeeping is my occupation. The preferred hobby for my children and me is to play baseball and I'm trying to make it a occupation. Years ago we moved to Puerto Rico and my family members enjoys it.<br><br>Feel free to surf to my homepage; [http://Leqo.ru/dietdelivery10007 http://Leqo.ru/dietdelivery10007] |
| [[File:Manometer anim 02.gif|thumb|200px|[[Pressure_measurement|Manometers]] are based on Hooke's law. The force created by gas [[pressure]] inside the coiled metal tube at right unwinds it by an amount proportional to the pressure.]]
| |
| [[File:Balancier avec ressort spiral.png|thumb|200px|The balance wheel at the core of many mechanical clocks and watches depends on Hooke's law. Since the torque generated by the coiled spring is proportional to the angle turned by the wheel, its oscillations have a nearly constant period.]]
| |
| | |
| '''Hooke's law''' is a principle of [[physics]] that states that the [[force]] <math>F</math> needed to extend or compress a [[spring (mechanics)|spring]] by some distance <math>X</math> is proportional to that distance. That is,
| |
| :<math>F = k X,</math>
| |
| where <math>k</math> is a constant factor characteristic of the spring, its [[stiffness]].
| |
| | |
| Hooke's equation in fact holds (to some extent) in many other situations where an [[elasticity (physics)|elastic]] body is deformed, such as wind blowing on a tall building, a musician plucking a [[string (music)|string]] of a [[violin]], or the filling of a [[party balloon]]. An elastic body or material for which this equation can be assumed is said to be [[linear elasticity|linear-elastic]] or '''Hookean'''.
| |
| | |
| Hookes's law is only a [[Taylor series|first order linear approximation]] to the real response of springs and other elastic bodies to applied forces. It must eventually fail once the forces exceed some limit, since no material can be compressed beyond a certain minimum size, or stretched beyond a maximum size, without some permanent deformation or change of state. In fact, many materials will noticeably deviate from Hooke's law well before those [[elastic limit]]s are reached.
| |
| | |
| On the other hand, Hooke's law is an accurate approximation for most solid bodies, as long as the forces and deformations are small enough. For this reason, Hooke's law is extensively used in all branches of [[science]] and [[engineering]], and is the foundation of many disciplines such as [[seismology]], [[molecular mechanics]] and [[acoustics]]. It is also the fundamental principle behind the [[spring scale]], the [[manometer]], and the [[balance wheel]] of the [[mechanical clock]].
| |
| | |
| The modern [[theory of elasticity]] generalizes Hooke's law to say that the [[deformation (mechanics)|strain]] (deformation) of an elastic object or material is proportional to the [[stress (mechanics)|stress]] applied to it. However, since general stresses and strains may have multiple independent components, the "proportionality factor" may no longer be just a single real number, but rather a [[linear map]] (a [[tensor]]) that can be represented by a [[matrix (mathematics)|matrix]] of real numbers.
| |
| | |
| In this general form, Hooke's law and [[Isaac Newton|Newton]]'s laws of [[static equilibrium]] make it possible to deduce the relation between strain and stress for complex objects in terms of intrinsic properties of the materials it is made of. For example, one can deduce that a [[homogeneous]] rod with uniform [[cross section (geometry)|cross section]] will behave like a simple spring when stretched, with a stiffness <math>k</math> directly proportional to its cross-section area and inversely proportional to its length.
| |
| | |
| Hooke's law is named after the 17th century British physicist [[Robert Hooke]]. He first stated this law in 1660 as a [[Latin]] [[anagram]],<ref>The anagram was given in alphabetical order, ''ceiiinosssttuv'', representing ''Ut tensio, sic vis'' — "As the extension, so the force": {{cite book|last=Petroski|first=Henry|authorlink=Henry Petroski|title=Invention by Design: How Engineers Get from Thought to Thing|year=1996|publisher=Harvard University Press|location=Cambridge, MA|isbn=0674463684|page=11}}</ref><ref>See [http://www.lindahall.org/events_exhib/exhibit/exhibits/civil/design.shtml], where one can find also an anagram for the [[catenary]].</ref> whose solution he published in 1678 as ''Ut tensio, sic vis''; literally translated as: "As the extension, so the force" or a more common meaning is "The extension is proportional to the force".
| |
| | |
| ==Formal definition==
| |
| ===For linear springs===
| |
| [[File:HookesLawForSpring-English.png|thumb|250px|Plot of applied force ''F'' vs. elongation ''X'' for a helical spring according to Hooke's law (red line) and what the actual plot might look like (dashed line). At bottom, pictures of spring states corresponding to some points of the plot; the middle one is in the relaxed state (no force applied).]]
| |
| Consider a simple [[helix|helical]] spring that has one end attached to some fixed object, while the free end is being pulled by a force whose magnitude is <math>F</math>. Suppose that the spring has reached a state of [[mechanical equilibrium|equilibrium]], where its length is not changing anymore. Let <math>X</math> be the amount by which the free end of the spring was displaced from its "relaxed" position (when it is not being stretched). Hooke's law states that
| |
| :<math>F= k X\,</math>
| |
| or, equivalently,
| |
| :<math>X= \frac{1}{k}F\,</math>
| |
| where <math>k</math> is a positive real number, characteristic of the spring. Moreover, the same formula holds when the spring is compressed, with <math>F</math> and <math>X</math> both negative in that case. According to this formula, the [[graph (mathematics)|graph]] of the applied force <math>F</math> as a function of the displacement <math>X</math> will be a straight line passing through the [[Cartesian coordinates|origin]], whose [[slope]] is <math>k</math>.
| |
| | |
| Hooke's law for a spring is often stated under the convention that <math>F</math> is the [[reaction force|restoring (reaction) force]] exerted by the spring on whatever is pulling its free end. In that case the equation becomes
| |
| :<math>F= -k X\,</math>
| |
| since the direction of the restoring force is opposite to that of the displacement.
| |
| | |
| ===General "scalar" springs===
| |
| Hooke's spring law usually applies to any elastic object, of arbitrary complexity, as long as both the deformation and the stress can be expressed by a single number that can be both positive and negative.
| |
| | |
| For example, when a block of rubber attached to two parallel plates is deformed by [[simple shear|shear]]ing, rather than stretching or compression, the shearing force <math>F</math> and the sideways displacement of the plates <math>X</math> obey Hooke's law (for small enough deformations).
| |
| | |
| Hooke's law also applies when a straight [[steel]] bar or [[concrete]] beam, supported at both ends, is bent by a weight <math>F</math> placed at some intermediate point. The displacement <math>X</math> in this case is the deviation of the beam, measured in the transversal direction, relative to its unloaded shape.
| |
| | |
| The law also applies when a stretched steel wire is twisted by pulling on a lever attached to one end. In this case the stress <math>F</math> can be taken as the force applied to the lever, and <math>X</math> as the distance traveled by it along its circular path. Or, equivalently, one can let <math>F</math> be the [[torque]] applied by the lever to the end of the wire, and <math>X</math> be the angle by which that end turns. In either case <math>F</math> is proportional to <math>X</math> (although the constant <math>k</math> is different in each case.)
| |
| | |
| ===Vector formulation===
| |
| In the case of a helical spring that is stretched or compressed along its [[axis (mathematics)|axis]], the applied (or restoring) force and the resulting elongation or compression have the same direction (which is the direction of said axis). Therefore, if <math>F</math> and <math>X</math> are defined as [[vector (mathematics)|vector]]s, Hooke's [[equation]] still holds, and says that the force vector is the [[displacement (vector)|elongation vector]] multiplied by a fixed [[scalar (mathematics)|scalar]].
| |
| | |
| ===General tensor form===
| |
| Some elastic bodies will deform in one direction when subjected to a force with a different direction. One example is a horizontal wood beam with non-square rectangular cross section that is bent by a transverse load that is neither vertical nor horizontal. In such cases, the ''magnitude'' of the displacement <math>X</math> will be proportional to the magnitude of the force <math>F</math>, as long as the direction of the latter remains the same (and its value is not too large); so the scalar version of Hooke's law <math>F = k X</math> will hold. However, the force and displacement ''vectors'' will not be scalar multiples of each other, since they have different directions. Moreover, the ratio <math>k</math> between their magnitudes will depend of the direction of the vector <math>F</math>.
| |
| | |
| Yet, in such cases there is often a fixed [[linear map|linear relation]] between the force and deformation vectors, as long as they are small enough. Namely, there is a [[function (mathematics)|function]] <math>\kappa</math> from vectors to vectors, such that <math>F = \kappa(X)</math>, and <math>\kappa(\alpha X_1 + \beta X_2) = \alpha \kappa(X_1) + \beta \kappa(X_2)</math> for any real numbers <math>\alpha,\beta</math> and any displacement vectors <math>X_1,X_2</math>. Such a function is called a (second-order) [[tensor]].
| |
| | |
| With respect to an arbitrary [[Cartesian coordinates|Cartesian coordinate system]], the force and displacement vectors can be represented by 1×3 [[matrix (mathematics)|matrices]] of real numbers. Then the tensor <math>\kappa</math> connecting them can be represented by a 3×3 matrix <math>\kappa</math> of real coefficients,
| |
| that, when [[matrix product|multiplied]] by the displacement vector, gives the force vector:
| |
| :<math>
| |
| F \;=\;
| |
| \begin{bmatrix} F_1\\ F_2 \\ F_3 \end{bmatrix} \;=\;
| |
| \begin{bmatrix}
| |
| \kappa_{1\,1}& \kappa_{1\,2}& \kappa_{1\,3}\\
| |
| \kappa_{2\,1}& \kappa_{2\,2}& \kappa_{2\,3}\\
| |
| \kappa_{3\,1}& \kappa_{3\,2}& \kappa_{3\,3}
| |
| \end{bmatrix}
| |
| \begin{bmatrix} X_1\\ X_2 \\ X_3 \end{bmatrix}
| |
| \;=\; \kappa X
| |
| </math>
| |
| That is,
| |
| :<math>F_i \;=\; \kappa_{i\,1} X_1 + \kappa_{i\,2} X_2 + \kappa_{i\,3} X_3 \;=\; \sum_{j=1}^3 \kappa_{i\,j} X_j</math>
| |
| for <math>i</math> equal to 1,2, and 3. Therefore, Hooke's law <math>F = \kappa X</math> can be said to hold also when <math>X</math> and <math>F</math> are vectors with variable directions, except that the stiffness of the object is a tensor <math>\kappa</math>, rather than a single real number <math>k</math>.
| |
| | |
| === Hooke's law for continuous media ===
| |
| {{main|linear elasticity}}
| |
| The stresses and strains of the material inside a [[continuum mechanics|continuous]] elastic material (such as a block of [[rubber]], the wall of a [[boiler]], or a [[steel]] bar) are connected by a linear relationship that is mathematically similar to Hooke's spring law, and is often referred to by that name.
| |
| | |
| However, the strain state in a solid medium around some point cannot be described by a single vector. The same parcel of material, no matter how small, can be compressed, stretched, and sheared at the same time, along different directions. Likewise, the stresses in that parcel can be at once pushing, pulling, and shearing.
| |
| | |
| In order to capture this complexity, the relevant state of the medium around a point must be represented by two second-order tensors, the [[strain tensor]] <math>\epsilon</math> (in lieu of the displacement <math>X</math>) and the [[Cauchy stress tensor|stress tensor]] <math>\sigma</math> (replacing the restoring force <math>F</math>).The analogous of Hooke's spring law for continuous media is then
| |
| :<math>
| |
| \sigma = -c \epsilon,
| |
| </math>
| |
| where <math>c</math> is a fourth-order tensor (that is, a linear map between second-order tensors) usually called the [[stiffness tensor]] or [[elasticity tensor]]. One may also write it as
| |
| :<math>
| |
| \epsilon = - s \sigma,
| |
| </math>
| |
| where the tensor <math>s</math>, called the [[stiffness tensor|compliance tensor]], represents the inverse of said linear map.
| |
| | |
| In a Cartesian coordinate system, the stress and strain tensors can be represented by 3×3 matrices
| |
| :<math>
| |
| \epsilon =
| |
| \begin{bmatrix}
| |
| \epsilon_{1\,1}& \epsilon_{1\,2}& \epsilon_{1\,3}\\
| |
| \epsilon_{2\,1}& \epsilon_{2\,2}& \epsilon_{2\,3}\\
| |
| \epsilon_{3\,1}& \epsilon_{3\,2}& \epsilon_{3\,3}
| |
| \end{bmatrix}
| |
| \quad\quad\quad\quad
| |
| \sigma =
| |
| \begin{bmatrix}
| |
| \sigma_{1\,1}& \sigma_{1\,2}& \sigma_{1\,3}\\
| |
| \sigma_{2\,1}& \sigma_{2\,2}& \sigma_{2\,3}\\
| |
| \sigma_{3\,1}& \sigma_{3\,2}& \sigma_{3\,3}
| |
| \end{bmatrix}
| |
| </math>
| |
| Being a linear mapping between the nine numbers <math>\sigma_{i\,j}</math> and the nine numbers <math>\epsilon_{k\,\ell}</math>, the stiffness tensor <math>c</math> is represented by a matrix of 3×3×3×3 = 81 real numbers <math>c_{i\,j\,k\,\ell}</math>. Hooke's law then says that
| |
| :<math>\sigma_{i\,j} = - \sum_{k=1}^3 \sum_{\ell=1}^3 c_{i\,j\,k\,\ell} \epsilon_{k\,\ell}</math>
| |
| where <math>i</math> and <math>j</math> are 1, 2, or 3.
| |
| | |
| All three tensors generally vary from point to point inside the medium, and may vary with time as well. The strain tensor <math>\epsilon</math> merely specifies the displacement of the medium particles in the neighborhood of the point, while the stress tensor <math>\sigma</math> specifies the forces that neighboring parcels of the medium are exerting on each other. Therefore, they are independent of the composition and physical state of the material. The stiffness tensor <math>c</math>, on the other hand, is a property of the material, and often depends on physical state variables such as [[temperature]], [[pressure]], and [[microstructure]].
| |
| | |
| Due to the inherent symmetries of <math>\sigma</math>, <math>\epsilon</math>, and <math>c</math>, only 21 elastic coefficients of the latter are independent. For [[isotropic]] media (which have the same physical properties in any direction), <math>c</math> can be reduced to only two independent numbers, the [[bulk modulus]] <math>K</math> and the [[shear modulus]] <math>G</math>, that quantify the material's resistance to changes in volume and to shearing deformations, respectively.
| |
| | |
| == Analogous laws ==
| |
| Since Hooke's law is a simple proportionality between two quantities, its formulas and consequences are mathematically similar to those of many other physical laws, such as those describing the motion of [[fluid]]s, or the [[ionic polarization|polarization]] of a [[dielectric]] by an [[electric field]].
| |
| | |
| In particular, the tensor equation <math>\sigma = c \epsilon</math> relating elastic stresses to strains is entirely similar to the equation <math>\tau = \mu \dot\epsilon</math> relating the [[viscous stress tensor]] <math>\tau</math> and the [[strain rate tensor]] <math>\dot\epsilon</math> in flows of [[viscosity|viscous]] fluids; although the former pertains to [[statics|static]] stresses (related to ''amount'' of deformation) while the latter pertains to [[dynamics (physics)|dynamical]] stresses (related to the ''rate'' of deformation).
| |
| | |
| ==Units of measurement==
| |
| In [[International System of Units|SI units]], displacements are measured in [[metre]]s (m), and forces in [[newton (unit)|newton]]s (N or kg·m/s<sup>2</sup>). Therefore the spring constant <math>k</math>, and each element of the tensor <math>\kappa</math>, is measured in newtons per metre (N/m), or [[kilogram]]s per [[second]] squared (kg/s<sup>2</sup>).
| |
| | |
| For continuous media, each element of the stress tensor <math>\sigma</math> is a force divided by an area; it is therefore measured in units of pressure, namely [[pascal (unit)|pascal]]s (Pa, or N/m<sup>2</sup>, or kg/m/s<sup>2</sup>. The elements of the strain tensor <math>\epsilon</math> are [[dimensionless]] (displacements divided by distances). Therefore the entries of <math>c_{ijk\ell}</math> are also expressed in units of pressure.
| |
| | |
| == General application to elastic materials ==
| |
| [[Image:Stress v strain A36 2.svg|thumb|right|300px|[[Stress–strain curve]] for low-carbon steel. Hooke's law is only valid for the portion of the curve between the origin and the yield point(2).<br />
| |
| 1. [[Ultimate strength]]<br />
| |
| 2. [[Yield strength]] – corresponds to yield point<br />
| |
| 3. [[Rupture (engineering)|Rupture]]<br />
| |
| 4. [[Strain hardening]] region<br />
| |
| 5. [[Necking (engineering)|Necking]] region<br />
| |
| A: Engineering stress (F/A<sub>0</sub>)<br />
| |
| B: True stress (F/A)]]
| |
| Objects that quickly regain their original shape after being deformed by a force, with the molecules or atoms of their material returning to the initial state of stable equilibrium, often obey Hooke's law.
| |
| | |
| Hooke's law only holds for some materials under certain loading conditions. Steel exhibits linear-elastic behavior in most engineering applications; Hooke's law is valid for it throughout its '''elastic range''' (i.e., for stresses below the [[Yield (engineering)|yield strength]]). For some other materials, such as aluminium, Hooke's law is only valid for a portion of the elastic range. For these materials a [[proportional limit]] stress is defined, below which the errors associated with the linear approximation are negligible.
| |
| | |
| Rubber is generally regarded as a "non-hookean" material because its elasticity is stress dependent and sensitive to temperature and loading rate.
| |
| | |
| Generalizations of Hooke's law for the case of [[Finite strain theory|large deformations]] is provided by models of [[neo-Hookean solid]]s and [[Mooney-Rivlin solid]]s.
| |
| | |
| ==Derived formulas==
| |
| === Tensional stiffness of a uniform bar ===
| |
| We may view a rod of any [[Elasticity (physics)|elastic]] material as a linear [[Spring (device)|spring]]. The rod has length ''L'' and cross-sectional area ''A''. Its extension (strain) is linearly proportional to its [[tensile stress]] ''σ'' by a constant factor <math>\varepsilon</math>, the inverse of its [[modulus of elasticity]] ''E'', such that
| |
| | |
| <math>E \varepsilon = [constant] = \sigma</math>.
| |
| | |
| In turn,
| |
| | |
| <math>\varepsilon = \frac{\Delta L}{L}</math> (''i.e.'', [change in length] as a fraction or percentage of total length),
| |
| | |
| and because
| |
| | |
| <math>\sigma = \frac{F}{A} </math> ,
| |
| | |
| such that
| |
| | |
| <math> \varepsilon = \frac{\sigma}{E} = \frac{(\frac{F}{A})}{E} = \frac{F}{A E}</math> ,
| |
| | |
| | |
| this relationship may also be expressed as
| |
| | |
| :<math>\Delta L = \varepsilon L = \frac{\sigma}{E} L = \frac{F}{A E} L = \frac{F L}{A E}</math> .
| |
| | |
| === Spring energy ===
| |
| The potential energy stored in a spring is given by
| |
| :<math> \mathrm{PE}={1\over2}kx^2</math>
| |
| which comes from adding up the energy it takes to incrementally compress the spring. That is, the integral of force over displacement. Since the external force has the same general direction as the displacement, the potential energy of a spring is always non-negative.
| |
| | |
| This potential "U" can be visualized as a [parabola] on the ''U''-''x'' plane such that <math>U(x)= {1\over2}kx^2</math>. As the spring is stretched in the positive x-direction, the potential energy increases parabolically (the same thing happens as the spring is compressed). Since the change in potential energy changes at a constant rate: <math> \frac{d^2}{dx^2}U=k</math>. Note that the change in the change in U is constant even when the displacement and acceleration are zero.
| |
| | |
| === Relaxed force constants (generalized compliance constants) ===
| |
| Relaxed force constants (the inverse of Grunenberg`s generalized compliance constants) are uniquely defined for molecular systems, in contradistinction to the usual “rigid” force constants, and thus their use allows meaningful correlations to be made between force fields calculated for reactants, transition state, and product of a chemical reaction. Just as the potential energy can be written as a quadratic form in the internal coordinates, so it can also be written in terms of generalized forces. The resulting coefficients are termed compliance constants. The suitability of relaxed force constants (inverse compliance constants) as covalent bond strength descriptors was demonstrated as early as ~1980, at the University of Cambridge. Recently, the suitability as non-covalent bond strength descriptors was demonstrated, too.<ref>Alla G. Ponomareva, Yevgen P. Yurenko , Roman O. Zhurakivsky, Tanja van Mourik and Dmytro M. Hovorun Phys. Chem. Chem. Phys., 2012,14, 6787-6795</ref>
| |
| | |
| === Harmonic oscillator ===
| |
| {{see also|Harmonic oscillator}}
| |
| [[File:Mass-spring-system.png|thumb|upright|A mass suspended by a spring is the classical example of an harmonic oscillator]]
| |
| A mass ''m'' is attached to the end of a spring is the classic example of a [[harmonic oscillator]]. By pulling slightly on the mass and then releasing it, the system will be set in [[sine wave|sinusoid]] oscillating motion about the equilibrium position. To the extent that the spring obeys Hooke's law, and that one can neglect [[friction]] and the mass of the spring, the amplitude of the oscillation will remain constant; and its [[frequency]] <math>f</math> will be independent of its amplitude, determined only by the mass and the stiffness of the spring:
| |
| :<math>f = {1 \over 2 \pi} \sqrt{k \over m}</math>
| |
| This phenomenon made possible the construction of accurate [[mechanical clock]]s and [[watch]]es that could be carried on ships and people's pockets.
| |
| <!--TO BE MOVED TO [[harmonic oscillator]]: either as an [[angular frequency]]
| |
| :<math>\omega = \sqrt{k \over m}</math>
| |
| or as a [[natural frequency]]
| |
| :<math>
| |
| -->
| |
| <!-- The various lattice spring constants are to be included later.-->
| |
| | |
| == Linear elasticity theory for continuous media ==
| |
| {{Einstein_summation_convention}}
| |
| <!-- MERGE THIS SECTION INTO [[linear elasticity]]-->
| |
| === Isotropic materials ===<!--- NOTE: there is a link from "Linear elasticity" to this section --->
| |
| | |
| (see [[viscosity]] for an analogous development for viscous fluids.)
| |
| | |
| Isotropic materials are characterized by properties which are independent of direction in space. Physical equations involving isotropic materials must therefore be independent of the coordinate system chosen to represent them. The strain tensor is a symmetric tensor. Since the [[Trace (linear algebra)|trace]] of any tensor is independent of any coordinate system, the most complete coordinate-free decomposition of a symmetric tensor is to represent it as the sum of a constant tensor and a traceless symmetric tensor.<ref>{{cite book|author=Symon, Keith |title=Mechanics |publisher=Addison-Wesley |location=Reading, MA |year=1971 |isbn=0-201-07392-7 }}
| |
| </ref>{{rp|Ch. 10}} Thus in [[Ricci calculus|index notation]]:
| |
| :<math>
| |
| \varepsilon_{ij} = \left(\tfrac{1}{3}\varepsilon_{kk}\delta_{ij}\right) +
| |
| \left(\varepsilon_{ij}-\tfrac{1}{3}\varepsilon_{kk}\delta_{ij}\right)
| |
| </math>
| |
| where <math>\delta_{ij}</math> is the [[Kronecker delta]]. In direct tensor notation:
| |
| :<math>
| |
| \boldsymbol{\varepsilon} = \mathrm{vol}(\boldsymbol{\varepsilon}) +
| |
| \mathrm{dev}(\boldsymbol{\varepsilon}) ~;~~
| |
| \mathrm{vol}(\boldsymbol{\varepsilon}) := \tfrac{1}{3}~\mathrm{tr}(\boldsymbol{\varepsilon})~\mathbf{I} ~;~~
| |
| \mathrm{dev}(\boldsymbol{\varepsilon}) := \boldsymbol{\varepsilon} - \mathrm{vol}(\boldsymbol{\varepsilon})
| |
| </math>
| |
| where <math>\mathbf{I}</math> is the second-order identity tensor.
| |
| The first term on the right is the constant tensor, also known as the '''volumetric strain tensor''', and the second term is the traceless symmetric tensor, also known as the '''deviatoric strain tensor''' or [[shear tensor]].
| |
| | |
| The most general form of Hooke's law for isotropic materials may now be written as a linear combination of these two tensors:
| |
| :<math>
| |
| \sigma_{ij}=3K\left(\tfrac{1}{3}\varepsilon_{kk}\delta_{ij}\right)
| |
| +2G\left(\varepsilon_{ij}-\tfrac{1}{3}\varepsilon_{kk}\delta_{ij}\right)\,~;~~
| |
| \boldsymbol{\sigma} = 3K~\mathrm{vol}(\boldsymbol{\varepsilon}) + 2G~\mathrm{dev}(\boldsymbol{\varepsilon})
| |
| </math>
| |
| where ''K'' is the [[bulk modulus]] and ''G'' is the [[shear modulus]].
| |
| | |
| Using the relationships between the [[elastic modulus|elastic moduli]], these equations may also be expressed in various other ways. A common form of Hooke's law for isotropic materials, expressed in direct tensor notation, is
| |
| <ref name=Simo98>{{Citation |last=Simo |first=J. C. |last2=Hughes |first2=T. J. R. |year=1998 |title=Computational Inelasticity |publisher=Springer |isbn=9780387975207 }}</ref>
| |
| :<math>
| |
| \boldsymbol{\sigma} = \lambda~\mathrm{tr}(\boldsymbol{\varepsilon})~\mathbf{I} + 2\mu~\boldsymbol{\varepsilon}
| |
| = \mathsf{c}:\boldsymbol{\varepsilon} ~;~~ \mathsf{c} = \lambda~\mathbf{I}\otimes\mathbf{I} + 2\mu~\mathsf{I}
| |
| </math>
| |
| where <math>\lambda := K - 2/3 G</math> and <math>\mu := G</math> are the [[Lamé constants]], <math>\mathbf{I}</math> is the second-rank identity tensor, and <math>\mathsf{I}</math> is the symmetric part of the fourth-rank identity tensor. In index notation:
| |
| :<math>
| |
| \sigma_{ij} = \lambda~\varepsilon_{kk}~\delta_{ij} + 2\mu~\varepsilon_{ij} = c_{ijk\ell}~\varepsilon_{k\ell} ~;~~ c_{ijk\ell} = \lambda~\delta_{ij}~\delta_{k\ell} + \mu~(\delta_{ik}~\delta_{j\ell} + \delta_{i\ell}~\delta_{jk})
| |
| </math>
| |
| The inverse relationship is<ref name=Milton02>{{Citation |last=Milton |first=Graeme W. |year=2002 |title=The Theory of Composites |publisher=Cambridge University Press |isbn=9780521781251 |series=Cambridge Monographs on Applied and Computational Mathematics }}</ref>
| |
| :<math>
| |
| \boldsymbol{\varepsilon} = \tfrac{1}{2\mu}~\boldsymbol{\sigma} - \tfrac{\lambda}{2\mu(3\lambda+2\mu)}~\mathrm{tr}(\boldsymbol{\sigma})~\mathbf{I} = \tfrac{1}{2G}~\boldsymbol{\sigma} + \left(\tfrac{1}{9K} - \tfrac{1}{6G}\right)~\mathrm{tr}(\boldsymbol{\sigma})~\mathbf{I}
| |
| </math>
| |
| Therefore the compliance tensor in the relation <math>\boldsymbol{\varepsilon} = \mathsf{s}:\boldsymbol{\sigma}</math> is
| |
| :<math>
| |
| \mathsf{s} = - \tfrac{\lambda}{2\mu(3\lambda+2\mu)}~\mathbf{I}\otimes\mathbf{I} + \tfrac{1}{2\mu}~\mathsf{I}
| |
| = \left(\tfrac{1}{9K} - \tfrac{1}{6G}\right)~\mathbf{I}\otimes\mathbf{I} + \tfrac{1}{2G}~\mathsf{I}
| |
| </math>
| |
| In terms of [[Young's modulus]] and [[Poisson's ratio]], Hooke's law for isotropic materials can then be expressed as
| |
| :<math>
| |
| \varepsilon_{ij}=\tfrac{1}{E}(\sigma_{ij}-\nu[\sigma_{kk}\delta_{ij}-\sigma_{ij}]) ~;~~
| |
| \boldsymbol{\varepsilon} = \tfrac{1}{E}(\boldsymbol{\sigma} - \nu[\mathrm{tr}(\boldsymbol{\sigma})~\mathbf{I} - \boldsymbol{\sigma}])
| |
| </math>
| |
| This is the form in which the strain is expressed in terms of the stress tensor in engineering. The expression in expanded form is
| |
| :<math>
| |
| \begin{align}
| |
| \varepsilon_{11} & = \tfrac{1}{E}\left[ \sigma_{11} - \nu(\sigma_{22}+\sigma_{33}) \right] \\
| |
| \varepsilon_{22} & = \tfrac{1}{E}\left[\sigma_{22} - \nu(\sigma_{11}+\sigma_{33}) \right] \\
| |
| \varepsilon_{33} & = \tfrac{1}{E}\left[\sigma_{33} - \nu(\sigma_{11}+\sigma_{22}) \right] \\
| |
| \varepsilon_{12} & = \tfrac{1}{2G}~\sigma_{12} ~;~~
| |
| \varepsilon_{13} = \tfrac{1}{2G}~\sigma_{13} ~;~~
| |
| \varepsilon_{23} = \tfrac{1}{2G}~\sigma_{23}
| |
| \end{align}
| |
| </math>
| |
| | |
| where ''E'' is the [[Young's modulus]] and <math>\nu</math> is [[Poisson's ratio]]. (See [[3-D elasticity]]).
| |
| | |
| :{| class="toccolours collapsible collapsed" width="60%" style="text-align:left"
| |
| !Derivation of Hooke's law in 3D
| |
| |-
| |
| |The 3-D form of Hooke's law can be derived using Poisson's ratio and the 1-D form of Hooke's law as follows.
| |
| Consider the strain and stress relation as a superposition of two effects: stretching in direction of the load (1) and shrinking (caused by the load) in perpendicular directions (2 and 3),
| |
| | |
| :<math>\varepsilon_1' = \frac{1}{E}\sigma_1</math>,
| |
| :<math>\varepsilon_2' = -\frac{\nu}{E}\sigma_1</math>,
| |
| :<math>\varepsilon_3' = -\frac{\nu}{E}\sigma_1</math>,
| |
| | |
| where <math>\nu</math> is the [[Poisson's ratio]] and <math>E</math> the [[Young's modulus]].
| |
| We get similar equations to the loads in directions 2 and 3,
| |
| | |
| :<math>\varepsilon_1'' = -\frac{\nu}{E}\sigma_2</math>,
| |
| :<math>\varepsilon_2'' = \frac{1}{E}\sigma_2</math>,
| |
| :<math>\varepsilon_3'' = -\frac{\nu}{E}\sigma_2</math>,
| |
| | |
| and
| |
| | |
| :<math>\varepsilon_1''' = -\frac{\nu}{E}\sigma_3</math>,
| |
| :<math>\varepsilon_2''' = -\frac{\nu}{E}\sigma_3</math>,
| |
| :<math>\varepsilon_3''' = \frac{1}{E}\sigma_3</math>.
| |
| | |
| Summing the three cases together (<math>\varepsilon_i = \varepsilon_i' + \varepsilon_i'' +\varepsilon_i'''</math>) we get
| |
| | |
| :<math>\varepsilon_1 = \frac{1}{E}(\sigma_1-\nu(\sigma_2+\sigma_3))</math>
| |
| :<math>\varepsilon_2 = \frac{1}{E}(\sigma_2-\nu(\sigma_1+\sigma_3))</math>
| |
| :<math>\varepsilon_3 = \frac{1}{E}(\sigma_3-\nu(\sigma_1+\sigma_2))</math>
| |
| | |
| or by adding and subtracting one <math>\nu\sigma</math>
| |
| | |
| :<math>\varepsilon_1 = \frac{1}{E}((1+\nu)\sigma_1-\nu(\sigma_1+\sigma_2+\sigma_3))</math>
| |
| :<math>\varepsilon_2 = \frac{1}{E}((1+\nu)\sigma_2-\nu(\sigma_1+\sigma_2+\sigma_3))</math>
| |
| :<math>\varepsilon_3 = \frac{1}{E}((1+\nu)\sigma_3-\nu(\sigma_1+\sigma_2+\sigma_3))</math>
| |
| | |
| and further we get by solving <math>\sigma_1</math>
| |
| | |
| :<math>\sigma_1 = \frac{E}{1+\nu}\varepsilon_1 + \frac{\nu}{1+\nu}(\sigma_1+\sigma_2+\sigma_3)</math>.
| |
| | |
| Calculating the sum
| |
| | |
| :<math>\sum_{i=1,2,3}\varepsilon_i = \frac{1}{E}((1+\nu)\sum_{i=1,2,3}\sigma_i - 3\nu(\sum_{i=1,2,3}\sigma_i)) = \frac{1-2\nu}{E}\sum_{i=1,2,3}\sigma_i</math>
| |
| :<math> \sigma_1 +\sigma_2+\sigma_3 = \frac{E}{1-2\nu}(\varepsilon_1 + \varepsilon_2 +\varepsilon_3)</math>
| |
| | |
| and substituting it to the equation solved for <math>\sigma_1</math> gives
| |
| | |
| :<math>\sigma_1 = \frac{E}{1+\nu}\varepsilon_1 + \frac{E\nu}{(1+\nu)(1-2\nu)}(\varepsilon_1 + \varepsilon_2 +\varepsilon_3)</math>,
| |
| :<math>\sigma_1 = 2\mu\varepsilon_1 + \lambda(\varepsilon_1 + \varepsilon_2 +\varepsilon_3)</math>,
| |
| | |
| where <math>\mu</math> and <math>\lambda</math> are the [[Lamé parameters]].
| |
| Similar treatment of directions 2 and 3 gives the Hooke's law in three dimensions.
| |
| |}
| |
| In matrix form, Hooke's law for isotropic materials can be written as
| |
| :<math>
| |
| \begin{bmatrix}\varepsilon_{11} \\ \varepsilon_{22} \\ \varepsilon_{33} \\ 2\varepsilon_{23} \\ 2\varepsilon_{31} \\ 2\varepsilon_{12} \end{bmatrix} =
| |
| \begin{bmatrix}\varepsilon_{11} \\ \varepsilon_{22} \\ \varepsilon_{33} \\ \gamma_{23} \\ \gamma_{31} \\ \gamma_{12} \end{bmatrix} =
| |
| \cfrac{1}{E}
| |
| \begin{bmatrix} 1 & -\nu & -\nu & 0 & 0 & 0 \\
| |
| -\nu & 1 & -\nu & 0 & 0 & 0 \\
| |
| -\nu & -\nu & 1 & 0 & 0 & 0 \\
| |
| 0 & 0 & 0 & 2(1+\nu) & 0 & 0 \\
| |
| 0 & 0 & 0 & 0 & 2(1+\nu) & 0 \\
| |
| 0 & 0 & 0 & 0 & 0 & 2(1+\nu) \end{bmatrix}
| |
| \begin{bmatrix}\sigma_{11} \\ \sigma_{22} \\ \sigma_{33} \\ \sigma_{23} \\ \sigma_{31} \\ \sigma_{12} \end{bmatrix}
| |
| </math>
| |
| where <math>\gamma_{ij} := 2\varepsilon_{ij}</math> is the '''engineering shear strain'''.
| |
| The inverse relation may be written as
| |
| :<math>
| |
| \begin{bmatrix}\sigma_{11} \\ \sigma_{22} \\ \sigma_{33} \\ \sigma_{23} \\ \sigma_{31} \\ \sigma_{12} \end{bmatrix}
| |
| = \cfrac{E}{(1+\nu)(1-2\nu)}
| |
| \begin{bmatrix} 1-\nu & \nu & \nu & 0 & 0 & 0 \\
| |
| \nu & 1-\nu & \nu & 0 & 0 & 0 \\
| |
| \nu & \nu & 1-\nu & 0 & 0 & 0 \\
| |
| 0 & 0 & 0 & (1-2\nu)/2 & 0 & 0 \\
| |
| 0 & 0 & 0 & 0 & (1-2\nu)/2 & 0 \\
| |
| 0 & 0 & 0 & 0 & 0 & (1-2\nu)/2 \end{bmatrix}
| |
| \begin{bmatrix}\varepsilon_{11} \\ \varepsilon_{22} \\ \varepsilon_{33} \\ 2\varepsilon_{23} \\ 2\varepsilon_{31} \\ 2\varepsilon_{12} \end{bmatrix}
| |
| </math>
| |
| | |
| which can be simplified thanks to the Lamé constants :
| |
| :<math>
| |
| \begin{bmatrix}\sigma_{11} \\ \sigma_{22} \\ \sigma_{33} \\ \sigma_{23} \\ \sigma_{31} \\ \sigma_{12} \end{bmatrix}
| |
| =
| |
| \begin{bmatrix} 2\mu+\lambda & \lambda & \lambda & 0 & 0 & 0 \\
| |
| \lambda & 2\mu+\lambda & \lambda & 0 & 0 & 0 \\
| |
| \lambda & \lambda & 2\mu+\lambda & 0 & 0 & 0 \\
| |
| 0 & 0 & 0 & \mu & 0 & 0 \\
| |
| 0 & 0 & 0 & 0 & \mu & 0 \\
| |
| 0 & 0 & 0 & 0 & 0 & \mu \end{bmatrix}
| |
| \begin{bmatrix}\varepsilon_{11} \\ \varepsilon_{22} \\ \varepsilon_{33} \\ 2\varepsilon_{23} \\ 2\varepsilon_{31} \\ 2\varepsilon_{12} \end{bmatrix}
| |
| </math>
| |
| | |
| ==== Plane stress ====
| |
| Under [[Plane_stress#Plane_stress|plane stress]] conditions <math>\sigma_{31} = \sigma_{13} = \sigma_{32} = \sigma_{23} = \sigma_{33} = 0</math>. In that case Hooke's law takes the form
| |
| :<math>
| |
| \begin{bmatrix}\varepsilon_{11} \\ \varepsilon_{22} \\ 2\varepsilon_{12} \end{bmatrix} = \cfrac{1}{E}
| |
| \begin{bmatrix} 1 & -\nu & 0 \\
| |
| -\nu & 1 & 0 \\
| |
| 0 & 0 & 2(1+\nu) \end{bmatrix}
| |
| \begin{bmatrix}\sigma_{11} \\ \sigma_{22} \\ \sigma_{12} \end{bmatrix}
| |
| </math>
| |
| The inverse relation is usually written in the reduced form
| |
| :<math>
| |
| \begin{bmatrix}\sigma_{11} \\ \sigma_{22} \\ \sigma_{12} \end{bmatrix}
| |
| = \cfrac{E}{1-\nu^2}
| |
| \begin{bmatrix} 1 & \nu & 0 \\
| |
| \nu & 1 & 0 \\
| |
| 0 & 0 & \cfrac{1-\nu}{2} \end{bmatrix}
| |
| \begin{bmatrix}\varepsilon_{11} \\ \varepsilon_{22} \\ 2\varepsilon_{12} \end{bmatrix}
| |
| </math>
| |
| | |
| === Anisotropic materials ===
| |
| The symmetry of the [[stress (physics)|Cauchy stress tensor]] (<math>\sigma_{ij} = \sigma_{ji}\,</math>) and the generalized Hooke's laws (<math>\sigma_{ij} = c_{ijk\ell}~ \epsilon_{k\ell}</math>) implies that <math>c_{ijk\ell} = c_{jik\ell}\,</math>. Similarly, the symmetry of the [[infinitesimal strain theory|infinitesimal strain tensor]] implies that <math>c_{ijk\ell} = c_{ij\ell k}\,</math>. These symmetries are called the '''minor symmetries''' of the '''stiffness tensor''' (<math>\mathsf{c}</math>). This reduces the number of elastic constants from 81 to 36.
| |
| | |
| If in addition, since the displacement gradient and the Cauchy stress are work conjugate, the stress-strain relation can be derived from a strain energy density functional (<math>U</math>), then
| |
| :<math>
| |
| \sigma_{ij} = \cfrac{\partial U}{\partial \epsilon_{ij}} \quad \implies \quad
| |
| c_{ijk\ell} = \cfrac{\partial^2 U}{\partial \epsilon_{ij}\partial \epsilon_{k\ell}}~.
| |
| </math>
| |
| The arbitrariness of the order of differentiation implies that <math>c_{ijk\ell} = c_{k\ell ij}\,</math>. These are called the '''major symmetries''' of the stiffness tensor. This reduces the number of elastic constants to 21 from 36. The major and minor symmetries indicate that the stiffness tensor has only 21 independent components.
| |
| | |
| ==== Matrix representation (stiffness tensor) ====
| |
| It is often useful to express the anisotropic form of Hooke's law in matrix notation, also called [[Voigt notation]]. To do this we take advantage of the symmetry of the stress and strain tensors and express them as six-dimensional vectors in an orthonormal coordinate system (<math>\mathbf{e}_1,\mathbf{e}_2,\mathbf{e}_3</math>) as
| |
| :<math>
| |
| [\boldsymbol{\sigma}] = \begin{bmatrix}\sigma_{11}\\ \sigma_{22} \\ \sigma_{33} \\ \sigma_{23} \\ \sigma_{31} \\ \sigma_{12} \end{bmatrix} \equiv
| |
| \begin{bmatrix} \sigma_1 \\ \sigma_2 \\ \sigma_3 \\ \sigma_4 \\ \sigma_5 \\ \sigma_6 \end{bmatrix} ~;~~
| |
| [\boldsymbol{\epsilon}] = \begin{bmatrix}\epsilon_{11}\\ \epsilon_{22} \\ \epsilon_{33} \\ 2\epsilon_{23} \\ 2\epsilon_{31} \\ 2\epsilon_{12} \end{bmatrix} \equiv
| |
| \begin{bmatrix} \epsilon_1 \\ \epsilon_2 \\ \epsilon_3 \\ \epsilon_4 \\ \epsilon_5 \\ \epsilon_6 \end{bmatrix}
| |
| </math>
| |
| Then the stiffness tensor (<math>\mathsf{c}</math>) can be expressed as
| |
| :<math>
| |
| [\mathsf{C}] = \begin{bmatrix} c_{1111} & c_{1122} & c_{1133} & c_{1123} & c_{1131} & c_{1112} \\
| |
| c_{2211} & c_{2222} & c_{2233} & c_{2223} & c_{2231} & c_{2212} \\
| |
| c_{3311} & c_{3322} & c_{3333} & c_{3323} & c_{3331} & c_{3312} \\
| |
| c_{2311} & c_{2322} & c_{2333} & c_{2323} & c_{2331} & c_{2312} \\
| |
| c_{3111} & c_{3122} & c_{3133} & c_{3123} & c_{3131} & c_{3112} \\
| |
| c_{1211} & c_{1222} & c_{1233} & c_{1223} & c_{1231} & c_{1212}
| |
| \end{bmatrix} \equiv \begin{bmatrix}
| |
| C_{11} & C_{12} & C_{13} & C_{14} & C_{15} & C_{16} \\
| |
| C_{12} & C_{22} & C_{23} & C_{24} & C_{25} & C_{26} \\
| |
| C_{13} & C_{23} & C_{33} & C_{34} & C_{35} & C_{36} \\
| |
| C_{14} & C_{24} & C_{34} & C_{44} & C_{45} & C_{46} \\
| |
| C_{15} & C_{25} & C_{35} & C_{45} & C_{55} & C_{56} \\
| |
| C_{16} & C_{26} & C_{36} & C_{46} & C_{56} & C_{66} \end{bmatrix}
| |
| </math>
| |
| and Hooke's law is written as
| |
| :<math>
| |
| [\boldsymbol{\sigma}] = [\mathsf{C}][\boldsymbol{\epsilon}] \qquad \text{or} \qquad \sigma_i = C_{ij} \epsilon_j ~.
| |
| </math>
| |
| Similarly the compliance tensor (<math>\mathsf{s}</math>) can be written as
| |
| :<math>
| |
| [\mathsf{S}] = \begin{bmatrix}
| |
| s_{1111} & s_{1122} & s_{1133} & 2s_{1123} & 2s_{1131} & 2s_{1112} \\
| |
| s_{2211} & s_{2222} & s_{2233} & 2s_{2223} & 2s_{2231} & 2s_{2212} \\
| |
| s_{3311} & s_{3322} & s_{3333} & 2s_{3323} & 2s_{3331} & 2s_{3312} \\
| |
| 2s_{2311} & 2s_{2322} & 2s_{2333} & 4s_{2323} & 4s_{2331} & 4s_{2312} \\
| |
| 2s_{3111} & 2s_{3122} & 2s_{3133} & 4s_{3123} & 4s_{3131} & 4s_{3112} \\
| |
| 2s_{1211} & 2s_{1222} & 2s_{1233} & 4s_{1223} & 4s_{1231} & 4s_{1212}
| |
| \end{bmatrix} \equiv \begin{bmatrix}
| |
| S_{11} & S_{12} & S_{13} & S_{14} & S_{15} & S_{16} \\
| |
| S_{12} & S_{22} & S_{23} & S_{24} & S_{25} & S_{26} \\
| |
| S_{13} & S_{23} & S_{33} & S_{34} & S_{35} & S_{36} \\
| |
| S_{14} & S_{24} & S_{34} & S_{44} & S_{45} & S_{46} \\
| |
| S_{15} & S_{25} & S_{35} & S_{45} & S_{55} & S_{56} \\
| |
| S_{16} & S_{26} & S_{36} & S_{46} & S_{56} & S_{66} \end{bmatrix}
| |
| </math>
| |
| | |
| ==== Change of coordinate system ====
| |
| If a linear elastic material is rotated from a reference configuration to another, then the material is symmetric with respect to the rotation if the components of the stiffness tensor in the rotated configuration are related to the components in the reference configuration by the relation<ref name=Slaughter>{{Citation |last=Slaughter |first=William S. |year=2001 |title=The Linearized Theory of Elasticity |publisher=Birkhäuser |isbn= 978-0817641177 }}</ref>
| |
| :<math>
| |
| c_{pqrs} = l_{pi}~l_{qj}~l_{rk}~l_{s\ell}~c_{ijk\ell}
| |
| </math>
| |
| where <math>l_{ab}</math> are the components of an [[orthogonal matrix|orthogonal rotation matrix]] <math>[L]</math>. The same relation also holds for inversions.
| |
| | |
| In matrix notation, if the transformed basis (rotated or inverted) is related to the reference basis by
| |
| :<math>
| |
| [\mathbf{e}_i'] = [L][\mathbf{e}_i]
| |
| </math>
| |
| then
| |
| :<math>
| |
| C_{ij}~\epsilon_i~\epsilon_j = C_{ij}'~\epsilon'_i~\epsilon'_j ~.
| |
| </math>
| |
| In addition, if the material is symmetric with respect to the transformation <math>[L]</math> then
| |
| :<math>
| |
| C_{ij} = C'_{ij} \quad \implies \quad C_{ij}~(\epsilon_i~\epsilon_j - \epsilon'_i~\epsilon'_j) = 0 ~.
| |
| </math>
| |
| | |
| ==== Orthotropic materials ====
| |
| {{Main|Orthotropic material}}
| |
| [[Orthotropic material]]s have three [[orthogonal]] [[plane of symmetry|planes of symmetry]]. If the basis vectors (<math>\mathbf{e}_1,\mathbf{e}_2,\mathbf{e}_3</math>) are normals to the planes of symmetry then the coordinate transformation relations imply that
| |
| :<math>
| |
| \begin{bmatrix} \sigma_1 \\ \sigma_2 \\ \sigma_3 \\ \sigma_4 \\ \sigma_5 \\ \sigma_6 \end{bmatrix} =
| |
| \begin{bmatrix}
| |
| C_{11} & C_{12} & C_{13} & 0 & 0 & 0 \\
| |
| C_{12} & C_{22} & C_{23} & 0 & 0 & 0 \\
| |
| C_{13} & C_{23} & C_{33} & 0 & 0 & 0 \\
| |
| 0 & 0 & 0 & C_{44} & 0 & 0 \\
| |
| 0 & 0 & 0 & 0 & C_{55} & 0 \\
| |
| 0 & 0 & 0 & 0 & 0 & C_{66} \end{bmatrix}
| |
| \begin{bmatrix} \epsilon_1 \\ \epsilon_2 \\ \epsilon_3 \\ \epsilon_4 \\ \epsilon_5 \\ \epsilon_6 \end{bmatrix}
| |
| </math>
| |
| The inverse of this relation is commonly written as<ref name=Boresi>Boresi, A. P, Schmidt, R. J. and Sidebottom, O. M., 1993, ''Advanced Mechanics of Materials'', Wiley.</ref>
| |
| :<math>
| |
| \begin{bmatrix}
| |
| \epsilon_{{\rm xx}} \\ \epsilon_{\rm yy} \\ \epsilon_{\rm zz} \\ 2\epsilon_{\rm yz} \\ 2\epsilon_{\rm zx} \\ 2\epsilon_{\rm xy}
| |
| \end{bmatrix}
| |
| = \begin{bmatrix}
| |
| \tfrac{1}{E_{\rm x}} & - \tfrac{\nu_{\rm yx}}{E_{\rm y}} & - \tfrac{\nu_{\rm zx}}{E_{\rm z}} & 0 & 0 & 0 \\
| |
| -\tfrac{\nu_{\rm xy}}{E_{\rm x}} & \tfrac{1}{E_{\rm y}} & - \tfrac{\nu_{\rm zy}}{E_{\rm z}} & 0 & 0 & 0 \\
| |
| -\tfrac{\nu_{\rm xz}}{E_{\rm x}} & - \tfrac{\nu_{\rm yz}}{E_{\rm y}} & \tfrac{1}{E_{\rm z}} & 0 & 0 & 0 \\
| |
| 0 & 0 & 0 & \tfrac{1}{G_{\rm yz}} & 0 & 0 \\
| |
| 0 & 0 & 0 & 0 & \tfrac{1}{G_{\rm zx}} & 0 \\
| |
| 0 & 0 & 0 & 0 & 0 & \tfrac{1}{G_{\rm xy}} \\
| |
| \end{bmatrix}
| |
| \begin{bmatrix}
| |
| \sigma_{\rm xx} \\ \sigma_{\rm yy} \\ \sigma_{\rm zz} \\ \sigma_{\rm yz} \\ \sigma_{\rm zx} \\ \sigma_{\rm xy}
| |
| \end{bmatrix}
| |
| </math>
| |
| where
| |
| :<math>{E}_{\rm i}\,</math> is the [[Young's modulus]] along axis <math>i</math>
| |
| :<math>G_{\rm ij}\,</math> is the [[shear modulus]] in direction <math>j</math> on the plane whose normal is in direction <math>i</math>
| |
| :<math>\nu_{\rm ij}\,</math> is the Poisson's ratio that corresponds to a contraction in direction <math>j</math> when an extension is applied in direction <math>i</math>.
| |
| | |
| Under ''plane stress'' conditions, <math>\sigma_{zz} = \sigma_{zx} = \sigma_{yz} = 0</math>, Hooke's law for an orthotropic material takes the form
| |
| :<math>
| |
| \begin{bmatrix}\varepsilon_{\rm xx} \\ \varepsilon_{\rm yy} \\ 2\varepsilon_{\rm xy} \end{bmatrix} =
| |
| \begin{bmatrix} \frac{1}{E_{\rm x}} & -\frac{\nu_{\rm yx}}{E_{\rm y}} & 0 \\
| |
| -\frac{\nu_{\rm xy}}{E_{\rm x}} & \frac{1}{E_{\rm y}} & 0 \\
| |
| 0 & 0 & \frac{1}{G_{\rm xy}} \end{bmatrix}
| |
| \begin{bmatrix}\sigma_{\rm xx} \\ \sigma_{\rm yy} \\ \sigma_{\rm xy} \end{bmatrix} \,.
| |
| </math>
| |
| The inverse relation is
| |
| :<math>
| |
| \begin{bmatrix}\sigma_{\rm xx} \\ \sigma_{\rm yy} \\ \sigma_{\rm xy} \end{bmatrix}
| |
| = \cfrac{1}{1-\nu_{\rm xy}\nu_{\rm yx}}
| |
| \begin{bmatrix} E_{\rm x} & \nu_{\rm yx}E_{\rm x} & 0 \\
| |
| \nu_{\rm xy}E_{\rm y} & E_{\rm y} & 0 \\
| |
| 0 & 0 & G_{\rm xy}(1-\nu_{\rm xy}\nu_{\rm yx}) \end{bmatrix}
| |
| \begin{bmatrix}\varepsilon_{\rm xx} \\ \varepsilon_{\rm yy} \\ 2\varepsilon_{\rm xy} \end{bmatrix} \,.
| |
| </math>
| |
| The transposed form of the above stiffness matrix is also often used.
| |
| | |
| ==== Transversely isotropic materials ====
| |
| A [[transversely isotropic]] material is symmetric with respect to a rotation about an [[axis of symmetry]]. For such a material, if <math>\mathbf{e}_3</math> is the axis of symmetry, Hooke's law can be expressed as
| |
| :<math>
| |
| \begin{bmatrix} \sigma_1 \\ \sigma_2 \\ \sigma_3 \\ \sigma_4 \\ \sigma_5 \\ \sigma_6 \end{bmatrix} =
| |
| \begin{bmatrix}
| |
| C_{11} & C_{12} & C_{13} & 0 & 0 & 0 \\
| |
| C_{12} & C_{11} & C_{13} & 0 & 0 & 0 \\
| |
| C_{13} & C_{13} & C_{33} & 0 & 0 & 0 \\
| |
| 0 & 0 & 0 & C_{44} & 0 & 0 \\
| |
| 0 & 0 & 0 & 0 & C_{44} & 0 \\
| |
| 0 & 0 & 0 & 0 & 0 & \tfrac{1}{2}(C_{11}-C_{12}) \end{bmatrix}
| |
| \begin{bmatrix} \epsilon_1 \\ \epsilon_2 \\ \epsilon_3 \\ \epsilon_4 \\ \epsilon_5 \\ \epsilon_6 \end{bmatrix}
| |
| </math>
| |
| | |
| More frequently, the <math>x \equiv \mathbf{e}_1</math> axis is taken to be the axis of symmetry and the inverse Hooke's law is written as
| |
| <ref name=Tan>Tan, S. C., 1994, ''Stress Concentrations in Laminated Composites'', Technomic Publishing Company, Lancaster, PA.</ref>
| |
| :<math>
| |
| \begin{bmatrix}
| |
| \epsilon_{{\rm xx}} \\ \epsilon_{\rm yy} \\ \epsilon_{\rm zz} \\ 2\epsilon_{\rm yz} \\ 2\epsilon_{\rm zx} \\ 2\epsilon_{\rm xy}
| |
| \end{bmatrix}
| |
| = \begin{bmatrix}
| |
| \tfrac{1}{E_{\rm x}} & - \tfrac{\nu_{\rm yx}}{E_{\rm y}} & - \tfrac{\nu_{\rm yx}}{E_{\rm y}} & 0 & 0 & 0 \\
| |
| -\tfrac{\nu_{\rm xy}}{E_{\rm x}} & \tfrac{1}{E_{\rm y}} & - \tfrac{\nu_{\rm yz}}{E_{\rm y}} & 0 & 0 & 0 \\
| |
| -\tfrac{\nu_{\rm xy}}{E_{\rm x}} & - \tfrac{\nu_{\rm yz}}{E_{\rm y}} & \tfrac{1}{E_{\rm y}} & 0 & 0 & 0 \\
| |
| 0 & 0 & 0 & \tfrac{2(1+\nu_{\rm yz})}{E_{\rm y}} & 0 & 0 \\
| |
| 0 & 0 & 0 & 0 & \tfrac{1}{G_{\rm xy}} & 0 \\
| |
| 0 & 0 & 0 & 0 & 0 & \tfrac{1}{G_{\rm xy}} \\
| |
| \end{bmatrix}
| |
| \begin{bmatrix}
| |
| \sigma_{\rm xx} \\ \sigma_{\rm yy} \\ \sigma_{\rm zz} \\ \sigma_{\rm yz} \\ \sigma_{\rm zx} \\ \sigma_{\rm xy}
| |
| \end{bmatrix}
| |
| </math>
| |
| | |
| == Thermodynamic basis ==
| |
| <!-- MERGE THIS SECTION INTO [[linear elasticity]]-->
| |
| Linear deformations of elastic materials can be approximated as [[adiabatic]]. Under these conditions and for quasistatic processes the [[first law of thermodynamics]] for a deformed body can be expressed as
| |
| :<math>
| |
| \delta W = \delta U\,
| |
| </math>
| |
| where <math>\delta U</math> is the increase in [[internal energy]] and <math>\delta W</math> is the [[work (physics)|work]] done by external forces. The work can be split into two terms
| |
| :<math>
| |
| \delta W = \delta W_s + \delta W_b\,
| |
| </math>
| |
| where <math>\delta W_s</math> is the work done by [[surface force]]s while <math>\delta W_b</math> is the work done by [[body force]]s. If <math>\delta\mathbf{u}</math> is a [[Calculus of variations|variation]] of the displacement field <math>\mathbf{u}</math> in the body, then the two external work terms can be expressed as
| |
| :<math>
| |
| \delta W_s = \int_{\partial\Omega} \mathbf{t}\cdot\delta\mathbf{u}~{\rm dS} ~;~~
| |
| \delta W_b = \int_{\Omega} \mathbf{b}\cdot\delta\mathbf{u}~{\rm dV}
| |
| </math>
| |
| where <math>\mathbf{t}</math> is the surface [[Stress (mechanics)|traction]] vector, <math>\mathbf{b}</math> is the body force vector, <math>\Omega\,</math> represents the body and <math>\partial\Omega</math> represents its surface. Using the relation between the [[stress (mechanics)|Cauchy stress]] and the surface traction, <math>\mathbf{t} = \mathbf{n}\cdot\boldsymbol{\sigma}</math> (where <math>\mathbf{n}</math> is the unit outward normal to <math>\partial\Omega</math>), we have
| |
| :<math>
| |
| \delta W = \delta U = \int_{\partial\Omega} (\mathbf{n}\cdot\boldsymbol{\sigma})\cdot\delta\mathbf{u}~{\rm dS} + \int_{\Omega} \mathbf{b}\cdot\delta\mathbf{u}~{\rm dV}
| |
| </math>
| |
| Converting the surface integral into a volume integral via the [[divergence theorem]] gives
| |
| :<math>
| |
| \delta U = \int_{\Omega} [\boldsymbol{\nabla}\cdot(\boldsymbol{\sigma}\cdot\delta\mathbf{u}) + \mathbf{b}\cdot\delta\mathbf{u}]~ {\rm dV} ~.
| |
| </math>
| |
| Using the symmetry of the Cauchy stress and the identity
| |
| :<math>
| |
| \boldsymbol{\nabla}\cdot(\boldsymbol{A}\cdot\mathbf{b}) = (\boldsymbol{\nabla}\cdot\boldsymbol{A})\cdot\mathbf{b}+
| |
| \tfrac{1}{2}[\boldsymbol{A}^T:\boldsymbol{\nabla}\mathbf{b}+
| |
| \boldsymbol{A}:(\boldsymbol{\nabla}\mathbf{b})^T]
| |
| </math>
| |
| we have the following
| |
| :<math>
| |
| \delta U = \int_{\Omega} [\boldsymbol{\sigma}:
| |
| \tfrac{1}{2}\{\boldsymbol{\nabla}\delta\mathbf{u}+(\boldsymbol{\nabla}\delta\mathbf{u})^T\} + \{\boldsymbol{\nabla}\cdot\boldsymbol{\sigma}+\mathbf{b}\}\cdot\delta\mathbf{u}]~{\rm dV} ~.
| |
| </math>
| |
| From the definition of [[infinitesimal strain theory|strain]] and from the equations of [[linear elasticity|equilibrium]] we have
| |
| :<math>
| |
| \delta\boldsymbol{\epsilon} = \tfrac{1}{2}[\boldsymbol{\nabla}\delta\mathbf{u}+(\boldsymbol{\nabla}\delta\mathbf{u})^T] ~;~~
| |
| \boldsymbol{\nabla}\cdot\boldsymbol{\sigma}+\mathbf{b}=\mathbf{0} ~.
| |
| </math>
| |
| Hence we can write
| |
| :<math>
| |
| \delta U = \int_{\Omega} \boldsymbol{\sigma}:\delta\boldsymbol{\epsilon}~{\rm dV}
| |
| </math>
| |
| and therefore the variation in the [[internal energy]] density is given by
| |
| :<math>
| |
| \delta U_0 = \boldsymbol{\sigma}:\delta\boldsymbol{\epsilon} ~.
| |
| </math>
| |
| An [[Elasticity (physics)|elastic]] material is defined as one in which the total internal energy is equal to the [[potential energy]] of the internal forces (also called the '''elastic strain energy'''). Therefore the internal energy density is a function of the strains, <math>U_0 = U_0(\boldsymbol{\epsilon})</math> and the variation of the internal energy can be expressed as
| |
| :<math>
| |
| \delta U_0 = \cfrac{\partial U_0}{\partial\boldsymbol{\epsilon}}:\delta\boldsymbol{\epsilon} ~.
| |
| </math>
| |
| Since the variation of strain is arbitrary, the stress-strain relation of an elastic material is given by
| |
| :<math>
| |
| \boldsymbol{\sigma} = \cfrac{\partial U_0}{\partial\boldsymbol{\epsilon}}~.
| |
| </math>
| |
| For a linear elastic material, the quantity <math>\partial U_0/\partial\boldsymbol{\epsilon}</math> is a linear function of <math>\boldsymbol{\epsilon}</math>, and can therefore be expressed as
| |
| :<math>
| |
| \boldsymbol{\sigma} = \mathsf{c}:\boldsymbol{\epsilon}
| |
| </math>
| |
| where <math>\mathsf{c}</math> is a fourth-rank tensor of material constants, also called the '''stiffness tensor'''. We can see why <math>\mathsf{c}</math> must be a fourth-rank tensor by noting that, for a linear elastic material,
| |
| :<math>
| |
| \cfrac{\partial}{\partial\boldsymbol{\epsilon}}[\boldsymbol{\sigma}(\boldsymbol{\epsilon})] = \text{constant} = \mathsf{c} \,.
| |
| </math>
| |
| In index notation
| |
| :<math>
| |
| \cfrac{\partial\sigma_{ij}}{\partial\epsilon_{k\ell}} = \text{constant} = c_{ijk\ell} \,.
| |
| </math>
| |
| Clearly, the right-hand side constant requires four indices and is a fourth-rank quantity. We can also see that this quantity must be a tensor because it is a linear transformation that takes the strain tensor to the stress tensor. We can also show that the constant obeys the tensor transformation rules for fourth-rank tensors.
| |
| | |
| == See also ==
| |
| {{Continuum mechanics}}
| |
| * [[Elastic potential energy]]
| |
| * [[List of scientific laws named after people]]
| |
| * [[Quadratic form]]
| |
| * [[Spring system]]
| |
| * [[Simple harmonic motion#Mass on a spring|Simple harmonic motion of a mass on a spring]]
| |
| * [[Sine wave]]
| |
| * [[Solid mechanics]]
| |
| * [[Spring pendulum]]
| |
| | |
| {{clear}}
| |
| | |
| == Notes ==
| |
| {{Reflist|30em}}
| |
| | |
| == References ==
| |
| * A.C. Ugural, S.K. Fenster, ''Advanced Strength and Applied Elasticity'', 4th ed
| |
| * [[:File:hookeslaw.ogg|Walter Lewin explains Hooke's law]]. From {{cite video
| |
| | people = [[Walter Lewin]] | date = 1 October 1999
| |
| | title = Hooke's Law, Simple Harmonic Oscillator. MIT Course 8.01: Classical Mechanics, Lecture 10.
| |
| | url = http://ocw.mit.edu/courses/physics/8-01-physics-i-classical-mechanics-fall-1999/video-lectures/lecture-10/
| |
| | format = ogg | medium = videotape | language = English
| |
| | publisher=[[MIT OpenCourseWare|MIT OCW]] | location = Cambridge, MA USA
| |
| | accessdate = 23 December 2010 | time = 1:21–10:10 | ref = lewin1
| |
| | quote = ...arguably the most important equation in all of Physics.
| |
| }}
| |
| * [[:File:hookeslawexample.ogg|A test of Hooke's law]]. From {{cite video
| |
| | people = [[Walter Lewin]] | date = 1 October 1999
| |
| | title = Hooke's Law, Simple Harmonic Oscillator. MIT Course 8.01: Classical Mechanics, Lecture 10.
| |
| | url = http://ocw.mit.edu/courses/physics/8-01-physics-i-classical-mechanics-fall-1999/video-lectures/lecture-10/
| |
| | format = ogg | medium = videotape | language = English
| |
| | publisher=MIT OCW | location = Cambridge, MA USA
| |
| | accessdate = 23 December 2010 | time = 10:10–16:33 | ref = lewin2
| |
| }}
| |
| | |
| == External links ==
| |
| * [http://webphysics.davidson.edu/applets/animator4/demo_hook.html Java Applet demonstrating Hooke's Law in motion]
| |
| * [http://webphysics.davidson.edu/applets/animator4/demo_hook.html Hooke's Law for General Stress/Strain Calculations ]
| |
| {{Elastic moduli}}
| |
| | |
| [[Category:1660 in science]]
| |
| [[Category:Springs (mechanical)]]
| |
| [[Category:Elasticity (physics)]]
| |
| [[Category:Solid mechanics]]
| |