قانون هوك: لزيادة القوة يزيد الامتداد.
مانومترات تعتمد في عملها على قانون هوك. القوة التي تشكلت بفعل
ضغط الغاز داخل الأنبوب المعدني الملفوف تتناسب مع الضغط.
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.
قانون هوك هو مبدأ في الفيزياء ينص على أن القوة
F
𝐹
{\displaystyle{\displaystyle F}}
التي يتغير بها الجسم (الإجهاد ) مرتبطة خطيًا بالقوة المسببة لهذا التغير (الشد). المواد التي ينطبق عليها قانون هوك تقريبًا هي مواد خطية المرونة .
سمى قانون هوك على اسم الفيزيائي الإنجليزي روبرت هوك الذي عاش في القرن السابع عشر. لقد ذكر هذا القانون في 1676 كبديل لاتيني, نشره في 1678 كجملة تعني :
"لزيادة القوة يزيد الامتداد"
من أجل الأنظمة التي يطبق عيها قانون هوك، الامتداد الناتج يتناسب مباشرة مع الحمل:
𝐅
=
k
𝐱
𝐅
𝑘
𝐱
{\displaystyle{\displaystyle\mathbf{F}=k\mathbf{x}\ }}
حيث :
𝐱
𝐱
{\displaystyle{\displaystyle\ \mathbf{x}}}
هي الفرق في المسافة بين موضع الجسم الجديد وموقعه الأصلي سواء كان مضغوطًا أو ممدودا"الازاحة الحاصلة" (عادة تقاس بالمتر )
𝐅
𝐅
{\displaystyle{\displaystyle\ \mathbf{F}}}
هي قوة الإعادة أو كما يطلق عليها القوة المشوهه للجسم اي معناها ان هذه القوة تغير من ابعاد الجسم ولو وصلت لحد معين قد تسبب تشوه للجسم اي لا يعود لشكله الاصلي قبل ان تؤثر عليه تلك القوة التي تمارسها المادة (عادة تقاس بالنيوتن )
و
𝐤
𝐤
{\displaystyle{\displaystyle\ \mathbf{k}}}
هو ثابت القوة ووحدته القوة إلى الطول (يقاس بالنيوتن لكل متر)
نظرة عامة
يتميز كثير من الأجسام، كالسلك الزنبركي او القضيب المعدني، بخاصية تسمى المرونة، فعندما يستطيل الجسم أو ينضغط تحت تأثير قوة مسلطة فإنه يميل إلى العودة إلى طوله الأصلي عند إزالة القوة. لنفرض مثلاً ان الزنبرك المبين بالشكل (1) طوله الأصلي L0 وانه قد استطال بمقدار LΔ تحت تأثير القوة المسلطة F. بدراسة هذا السلوك وجد روبرت هوك (1635 - 1703) أن الاستطالة تتضاعف مرتين إذا تضاعفت القوة المسلطة مرتين، بشرط ألا تكون الاستطالة كبيرة جداً، أي ان L α FΔ عموماً. وقد وضع هوك اكتشافاته هذه في صورة قاعدة تعرف الآن بقانون هوك:
عندما يتمدد جسم مرن أو يتشوه بأي صورة اخرى فإن مقدار التشوه يتناسب خطياً مع القوة المشوهة.
ولكن عند امتداد (استطالة) الزنبرك بمقدار كبير بحيث يتعدى ما يعرف بحد المرونة فإن ينحرف عن هذا التناسب الطردي بين LΔ و F وعلاوة على ذلك سنلاحظ أن الزنبرك لن يعود إلى طوله الأصلي عند إزالة القوة المسلطة.[1]
وعند استبدال الزنبرك المبين بالشكل (1) بقضيب مصمت سنجد أيضاً أن القضيب يتبع قانون هوك. وبالرغم من أن الاستطالة النسبية للقضيب أصغر كثيراً من قيمتها في حالة الزنبرك فإن القضيب يستطيل بانتظام بما يتفق مع قانون هولك ، ولكن قيم الاستطالة تكون أصغر مما في حالة الزنبرك؛ ويوضح الشكل (2) السلوك المشاهد عملياً في تجربة نموذجية من هذا النوع. لاحظ ان قانون هوك ينطبق في المنطقة المرنة فقط ، وسوف يفترض في المناقشة الآتية أن القوة والاستطالة صغيران بحيث لا يتعدى تشوه المادة حد مرونتها.
لاستخدام قانون هوك في وصف الخواص المرنة للجوامد سوف نستخدم مصطلحين هامين هما الإجهاد والانفعال ، وسنقوم بتعريف هاتين الكميتين بمساعدة تجربة الاستطالة ( او الشد) المبينة بالشكل (3). في هذه التجربة تؤثر القوة الشادة (المطيلة) F عمودياً على المساحة الطرفية A لقضيب طوله الأصلي L0 فيستطيل القضيب نتيجة لذلك بمقدار LΔ. يعرف الإجهاد الناتج عن F كالتالي:
ويعرف انفعال القضيب في الشكل 3)) كما يلي:
شكل 3: إجهاد الشد وإجهاد الضغط في حالة قضيب منتظم الإجهاد هو F/A والانفعال هو L / L0Δ.
وقد عرف الانفعال بالنسبة L / L0Δ، بدلا ً من LΔ، لأن أي جسم مرن يستطيع بمقدار يتناسب طردياً مع طوله الأصلي. وبقسمة LΔ على L0 نكون قد تخلصنا من تأثير طول الجسم على الاستطالة، وهو تأثير لا يمثل أي أهمية فيما يتعلق بخواص مادة القضيب ذاتها.
ونظراً لأن الانفعال نسبة بين طولين فإنه كمية ليست لها وحدات. وسنرى مؤخراً في هذا القسم أن هناك انواعاً اخرى من الانفعال ،وهذا يتوقف على الناحية الهندسية للموقف. اما في هذه الحالة الحالية فإننا نتحدث عن انفعال شد. ولكن إذا ضغط القضيب في اتجاه مواز لطوله فإن الانفعال، طبقاً للتعريف، سيكون أيضاً هو النسبة بين التغير في الطول والطول الاصلي.
الآن يمكننا إعادة صياغة قانون هوك. ذلك أن الإجهاد مقياس للقوة المشوهة والانفعال مقياس للتشوه. وعليه يمكن كتابة قانون هوك على الصورة:
(الانفعال) (ثابت) = الإجهاد
وبهذه الصورة يمكن تطبيق قانون هوك على مواقف كثيرة تختلف عن استطالة القضيب، وقد أثبتت تجارب هوك أن هذا القانون صالح للتطبيق في حالات استطالة وانحناء وفي العديد من الزنبركات والأجسام الأخرى. وكما أوضحنا سابقاً فإن قانون هولك ينطبق طبعاً في المنطقة المرنة من التشوهات فقط.
يعتمد ثابت التناسب في المعادلة (3) على طبيعة المادة ونوع التشوه الذي تعانيه، وهو يعرف بمعامل مرونة المادة. إذن ، طبقاً للتعريف:
الاجهاد/الانفاعل= معامل المرونة
وحيث أن الانفعال كمية ليس لها وحدات، فإن وحدات معامل المرونة هي نفس وحدات الإجهاد. لاحظ ان معامل المرونة يكون كبيراً عندما يسبب الإجهاد الكبير انفعالاً صغيراً فقط. وعليه فإن معامل المرونة مقياس لجسوءة المادة. وهناك، وفي الواقع، عدد انواع من معاملات المرونة ، وهذا يتوقف على تفاصيل الطريقة التي تستطيع بها المادة أو تنحني او تتشوه بأي طريقة أخرى من الطرق.
التعريف الرسمي
الزنبرك الخطي
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 helical spring that has one end attached to some fixed object, while the free end is being pulled by a force whose magnitude is
F
𝐹
{\displaystyle{\displaystyle F}}
. Suppose that the spring has reached a state of equilibrium , where its length is not changing anymore. Let
X
𝑋
{\displaystyle{\displaystyle X}}
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
F
=
k
X
𝐹
𝑘
𝑋
{\displaystyle{\displaystyle F=kX\,}}
or, equivalently,
X
=
1
k
F
𝑋
1
𝑘
𝐹
{\displaystyle{\displaystyle X={\frac{1}{k}}F\,}}
where
k
𝑘
{\displaystyle{\displaystyle k}}
is a positive real number, characteristic of the spring. Moreover, the same formula holds when the spring is compressed, with
F
𝐹
{\displaystyle{\displaystyle F}}
and
X
𝑋
{\displaystyle{\displaystyle X}}
both negative in that case. According to this formula, the graph of the applied force
F
𝐹
{\displaystyle{\displaystyle F}}
as a function of the displacement
X
𝑋
{\displaystyle{\displaystyle X}}
will be a straight line passing through the origin , whose slope is
k
𝑘
{\displaystyle{\displaystyle k}}
.
Hooke's law for a spring is often stated under the convention that
F
𝐹
{\displaystyle{\displaystyle F}}
is the restoring (reaction) force exerted by the spring on whatever is pulling its free end. في تلك الحالة تصبح المعادلة:
F
=
-
k
X
𝐹
𝑘
𝑋
{\displaystyle{\displaystyle F=-kX\,}}
since the direction of the restoring force is opposite to that of the displacement.
الزنبرك "العددي" العام
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.
صياغة المتجه
In the case of a helical spring that is stretched or compressed along its 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
F
𝐹
{\displaystyle{\displaystyle F}}
and
X
𝑋
{\displaystyle{\displaystyle X}}
are defined as vectors , Hooke's equation still holds, and says that the force vector is the elongation vector multiplied by a fixed scalar .
General tensor form
With respect to an arbitrary Cartesian coordinate system , the force and displacement vectors can be represented by 3×1 matrices of real numbers. Then the tensor
κ
𝜅
{\displaystyle{\displaystyle\kappa}}
connecting them can be represented by a 3×3 matrix
κ
𝜅
{\displaystyle{\displaystyle\kappa}}
of real coefficients,
that, when multiplied by the displacement vector, gives the force vector:
F
=
[
F
1
F
2
F
3
]
=
[
κ
1 1
κ
1 2
κ
1 3
κ
2 1
κ
2 2
κ
2 3
κ
3 1
κ
3 2
κ
3 3
]
[
X
1
X
2
X
3
]
=
κ
X
𝐹
subscript
𝐹
1
subscript
𝐹
2
subscript
𝐹
3
subscript
𝜅
1 1
subscript
𝜅
1 2
subscript
𝜅
1 3
subscript
𝜅
2 1
subscript
𝜅
2 2
subscript
𝜅
2 3
subscript
𝜅
3 1
subscript
𝜅
3 2
subscript
𝜅
3 3
subscript
𝑋
1
subscript
𝑋
2
subscript
𝑋
3
𝜅
𝑋
{\displaystyle{\displaystyle 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}}
That is,
F
i
=
κ
i
1
X
1
+
κ
i
2
X
2
+
κ
i
3
X
3
=
∑
j
=
1
3
κ
i
j
X
j
subscript
𝐹
𝑖
subscript
𝜅
𝑖
1
subscript
𝑋
1
subscript
𝜅
𝑖
2
subscript
𝑋
2
subscript
𝜅
𝑖
3
subscript
𝑋
3
superscript
subscript
𝑗
1
3
subscript
𝜅
𝑖
𝑗
subscript
𝑋
𝑗
{\displaystyle{\displaystyle 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}}}
for
i
𝑖
{\displaystyle{\displaystyle i}}
equal to 1,2, and 3. Therefore, Hooke's law
F
=
κ
X
𝐹
𝜅
𝑋
{\displaystyle{\displaystyle F=\kappa X}}
can be said to hold also when
X
𝑋
{\displaystyle{\displaystyle X}}
and
F
𝐹
{\displaystyle{\displaystyle F}}
are vectors with variable directions, except that the stiffness of the object is a tensor
κ
𝜅
{\displaystyle{\displaystyle\kappa}}
, rather than a single real number
k
𝑘
{\displaystyle{\displaystyle k}}
.
قانون هوك للوسائط المستمرة
σ
=
-
c
ϵ
,
𝜎
𝑐
italic-ϵ
{\displaystyle{\displaystyle\sigma=-c\epsilon,}}
where
c
𝑐
{\displaystyle{\displaystyle c}}
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
ϵ
=
-
s
σ
,
italic-ϵ
𝑠
𝜎
{\displaystyle{\displaystyle\epsilon=-s\sigma,}}
where the tensor
s
𝑠
{\displaystyle{\displaystyle s}}
, called the 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
ϵ
=
[
ϵ
1 1
ϵ
1 2
ϵ
1 3
ϵ
2 1
ϵ
2 2
ϵ
2 3
ϵ
3 1
ϵ
3 2
ϵ
3 3
]
σ
=
[
σ
1 1
σ
1 2
σ
1 3
σ
2 1
σ
2 2
σ
2 3
σ
3 1
σ
3 2
σ
3 3
]
formulae-sequence
italic-ϵ
subscript
italic-ϵ
1 1
subscript
italic-ϵ
1 2
subscript
italic-ϵ
1 3
subscript
italic-ϵ
2 1
subscript
italic-ϵ
2 2
subscript
italic-ϵ
2 3
subscript
italic-ϵ
3 1
subscript
italic-ϵ
3 2
subscript
italic-ϵ
3 3
𝜎
subscript
𝜎
1 1
subscript
𝜎
1 2
subscript
𝜎
1 3
subscript
𝜎
2 1
subscript
𝜎
2 2
subscript
𝜎
2 3
subscript
𝜎
3 1
subscript
𝜎
3 2
subscript
𝜎
3 3
{\displaystyle{\displaystyle\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}}}}
Being a linear mapping between the nine numbers
σ
i
j
subscript
𝜎
𝑖
𝑗
{\displaystyle{\displaystyle\sigma_{i\,j}}}
and the nine numbers
ϵ
k
ℓ
subscript
italic-ϵ
𝑘
ℓ
{\displaystyle{\displaystyle\epsilon_{k\,\ell}}}
, the stiffness tensor
c
𝑐
{\displaystyle{\displaystyle c}}
is represented by a matrix of 3×3×3×3 = 81 real numbers
c
i
j
k
ℓ
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
{\displaystyle{\displaystyle c_{i\,j\,k\,\ell}}}
. Hooke's law then says that
σ
i
j
=
-
∑
k
=
1
3
∑
ℓ
=
1
3
c
i
j
k
ℓ
ϵ
k
ℓ
subscript
𝜎
𝑖
𝑗
superscript
subscript
𝑘
1
3
superscript
subscript
ℓ
1
3
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
subscript
italic-ϵ
𝑘
ℓ
{\displaystyle{\displaystyle\sigma_{i\,j}=-\sum_{k=1}^{3}\sum_{\ell=1}^{3}c_{i%
\,j\,k\,\ell}\epsilon_{k\,\ell}}}
where
i
𝑖
{\displaystyle{\displaystyle i}}
and
j
𝑗
{\displaystyle{\displaystyle j}}
are 1, 2, or 3.
قوانين مماثلة
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 fluids , or the polarization of a dielectric by an electric field .
In particular, the tensor equation
σ
=
c
ϵ
𝜎
𝑐
italic-ϵ
{\displaystyle{\displaystyle\sigma=c\epsilon}}
relating elastic stresses to strains is entirely similar to the equation
τ
=
μ
ϵ
˙
𝜏
𝜇
˙
italic-ϵ
{\displaystyle{\displaystyle\tau=\mu{\dot{\epsilon}}}}
relating the viscous stress tensor
τ
𝜏
{\displaystyle{\displaystyle\tau}}
and the strain rate tensor
ϵ
˙
˙
italic-ϵ
{\displaystyle{\displaystyle{\dot{\epsilon}}}}
in flows of viscous fluids; although the former pertains to static stresses (related to amount of deformation) while the latter pertains to dynamical stresses (related to the rate of deformation).
وحدات القياس
In SI units , displacements are measured in metres (m), and forces in newtons (N or kg·m/s2 ). Therefore the spring constant
k
𝑘
{\displaystyle{\displaystyle k}}
, and each element of the tensor
κ
𝜅
{\displaystyle{\displaystyle\kappa}}
, is measured in newtons per metre (N/m), or kilograms per second squared (kg/s2 ).
For continuous media, each element of the stress tensor
σ
𝜎
{\displaystyle{\displaystyle\sigma}}
is a force divided by an area; it is therefore measured in units of pressure, namely pascals (Pa, or N/m2 , or kg/m/s2 . The elements of the strain tensor
ϵ
italic-ϵ
{\displaystyle{\displaystyle\epsilon}}
are dimensionless (displacements divided by distances). Therefore the entries of
c
i
j
k
ℓ
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
{\displaystyle{\displaystyle c_{ijk\ell}}}
are also expressed in units of pressure.
General application to elastic materials
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.
صيغ مشتقة
Tensional stiffness of a uniform bar
We may view a rod of any elastic material as a linear 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
ε
𝜀
{\displaystyle{\displaystyle\varepsilon}}
, the inverse of its modulus of elasticity E , such that
E
ε
=
[
c
o
n
s
t
a
n
t
]
=
σ
𝐸
𝜀
delimited-[]
𝑐
𝑜
𝑛
𝑠
𝑡
𝑎
𝑛
𝑡
𝜎
{\displaystyle{\displaystyle E\varepsilon=[constant]=\sigma}}
.
In turn,
ε
=
Δ
L
L
𝜀
Δ
𝐿
𝐿
{\displaystyle{\displaystyle\varepsilon={\frac{\Delta L}{L}}}}
(i.e. , [change in length] as a fraction or percentage of total length),
and because
σ
=
F
A
𝜎
𝐹
𝐴
{\displaystyle{\displaystyle\sigma={\frac{F}{A}}}}
,
such that
ε
=
σ
E
=
(
F
A
)
E
=
F
A
E
𝜀
𝜎
𝐸
𝐹
𝐴
𝐸
𝐹
𝐴
𝐸
{\displaystyle{\displaystyle\varepsilon={\frac{\sigma}{E}}={\frac{({\frac{F}{A%
}})}{E}}={\frac{F}{AE}}}}
,
this relationship may also be expressed as
Δ
L
=
ε
L
=
σ
E
L
=
F
A
E
L
=
F
L
A
E
Δ
𝐿
𝜀
𝐿
𝜎
𝐸
𝐿
𝐹
𝐴
𝐸
𝐿
𝐹
𝐿
𝐴
𝐸
{\displaystyle{\displaystyle\Delta L=\varepsilon L={\frac{\sigma}{E}}L={\frac{%
F}{AE}}L={\frac{FL}{AE}}}}
.
طاقة الزنبرك
The potential energy stored in a spring is given by
PE
=
1
2
k
x
2
PE
1
2
𝑘
superscript
𝑥
2
{\displaystyle{\displaystyle\mathrm{PE}={1\over 2}kx^{2}}}
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.
Harmonic oscillator
A mass suspended by a spring is the classical example of a harmonic oscillator
Rotation in Gravity-Free Space
If the mass m was attached to a spring with force constant k and rotating in free space, the spring tension (Ft ) would balance the required centripetal force (Fc ) as follows -
F
t
=
k
x
subscript
𝐹
𝑡
𝑘
𝑥
{\displaystyle{\displaystyle F_{t}=kx}}
F
c
=
m
ω
2
r
subscript
𝐹
𝑐
𝑚
superscript
𝜔
2
𝑟
{\displaystyle{\displaystyle F_{c}=m\omega^{2}r}}
Since Ft = Fc and x = r, therefore:
k
=
m
ω
2
𝑘
𝑚
superscript
𝜔
2
{\displaystyle{\displaystyle k=m\omega^{2}}}
Given that
ω
=
2
π
f
𝜔
2
𝜋
𝑓
{\displaystyle{\displaystyle\omega=2\pi f}}
, this leads to the same frequency equation as above -
f
=
1
2
π
k
m
𝑓
1
2
𝜋
𝑘
𝑚
{\displaystyle{\displaystyle f={1\over 2\pi}{\sqrt{k\over m}}}}
Linear elasticity theory for continuous media
Note: the Einstein summation convention of summing on repeated indices is used below.
Isotropic materials
(see viscosity for an analogous development for viscous fluids.)
Thus in index notation :
ε
i
j
=
(
1
3
ε
k
k
δ
i
j
)
+
(
ε
i
j
-
1
3
ε
k
k
δ
i
j
)
subscript
𝜀
𝑖
𝑗
1
3
subscript
𝜀
𝑘
𝑘
subscript
𝛿
𝑖
𝑗
subscript
𝜀
𝑖
𝑗
1
3
subscript
𝜀
𝑘
𝑘
subscript
𝛿
𝑖
𝑗
{\displaystyle{\displaystyle\varepsilon_{ij}=\left({\tfrac{1}{3}}\varepsilon_{%
kk}\delta_{ij}\right)+\left(\varepsilon_{ij}-{\tfrac{1}{3}}\varepsilon_{kk}%
\delta_{ij}\right)}}
where
δ
i
j
subscript
𝛿
𝑖
𝑗
{\displaystyle{\displaystyle\delta_{ij}}}
is the Kronecker delta . In direct tensor notation:
𝜺
=
vol
(
𝜺
)
+
dev
(
𝜺
)
;
vol
(
𝜺
)
:=
1
3
tr
(
𝜺
)
𝐈
;
dev
(
𝜺
)
:=
𝜺
-
vol
(
𝜺
)
formulae-sequence
𝜺
vol
𝜺
dev
𝜺
formulae-sequence
assign
vol
𝜺
1
3
tr
𝜺
𝐈
assign
dev
𝜺
𝜺
vol
𝜺
{\displaystyle{\displaystyle{\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}})}}
where
𝐈
𝐈
{\displaystyle{\displaystyle\mathbf{I}}}
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:
σ
i
j
=
3
K
(
1
3
ε
k
k
δ
i
j
)
+
2
G
(
ε
i
j
-
1
3
ε
k
k
δ
i
j
)
;
𝝈
=
3
K
vol
(
𝜺
)
+
2
G
dev
(
𝜺
)
formulae-sequence
subscript
𝜎
𝑖
𝑗
3
𝐾
1
3
subscript
𝜀
𝑘
𝑘
subscript
𝛿
𝑖
𝑗
2
𝐺
subscript
𝜀
𝑖
𝑗
1
3
subscript
𝜀
𝑘
𝑘
subscript
𝛿
𝑖
𝑗
𝝈
3
𝐾
vol
𝜺
2
𝐺
dev
𝜺
{\displaystyle{\displaystyle\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}})}}
where K is the bulk modulus and G is the shear modulus .
Using the relationships between the 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
[2]
𝝈
=
λ
tr
(
𝜺
)
𝐈
+
2
μ
𝜺
=
𝖼
:
𝜺
;
𝖼
=
λ
𝐈
⊗
𝐈
+
2
μ
𝖨
:
𝝈
𝜆
tr
𝜺
𝐈
2
𝜇
𝜺
𝖼
𝜺
𝖼
tensor-product
𝜆
𝐈
𝐈
2
𝜇
𝖨
{\displaystyle{\displaystyle{\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}}}}
where
λ
:=
K
-
2
/
3
G
assign
𝜆
𝐾
2
3
𝐺
{\displaystyle{\displaystyle\lambda:=K-2/3G}}
and
μ
:=
G
assign
𝜇
𝐺
{\displaystyle{\displaystyle\mu:=G}}
are the Lamé constants ,
𝐈
𝐈
{\displaystyle{\displaystyle\mathbf{I}}}
is the second-rank identity tensor, and
𝖨
𝖨
{\displaystyle{\displaystyle{\mathsf{I}}}}
is the symmetric part of the fourth-rank identity tensor. In index notation:
σ
i
j
=
λ
ε
k
k
δ
i
j
+
2
μ
ε
i
j
=
c
i
j
k
ℓ
ε
k
ℓ
;
c
i
j
k
ℓ
=
λ
δ
i
j
δ
k
ℓ
+
μ
(
δ
i
k
δ
j
ℓ
+
δ
i
ℓ
δ
j
k
)
formulae-sequence
subscript
𝜎
𝑖
𝑗
𝜆
subscript
𝜀
𝑘
𝑘
subscript
𝛿
𝑖
𝑗
2
𝜇
subscript
𝜀
𝑖
𝑗
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
subscript
𝜀
𝑘
ℓ
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
𝜆
subscript
𝛿
𝑖
𝑗
subscript
𝛿
𝑘
ℓ
𝜇
subscript
𝛿
𝑖
𝑘
subscript
𝛿
𝑗
ℓ
subscript
𝛿
𝑖
ℓ
subscript
𝛿
𝑗
𝑘
{\displaystyle{\displaystyle\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})}}
The inverse relationship is[3]
𝜺
=
1
2
μ
𝝈
-
λ
2
μ
(
3
λ
+
2
μ
)
tr
(
𝝈
)
𝐈
=
1
2
G
𝝈
+
(
1
9
K
-
1
6
G
)
tr
(
𝝈
)
𝐈
𝜺
1
2
𝜇
𝝈
𝜆
2
𝜇
3
𝜆
2
𝜇
tr
𝝈
𝐈
1
2
𝐺
𝝈
1
9
𝐾
1
6
𝐺
tr
𝝈
𝐈
{\displaystyle{\displaystyle{\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}}}
Therefore the compliance tensor in the relation
𝜺
=
𝗌
:
𝝈
:
𝜺
𝗌
𝝈
{\displaystyle{\displaystyle{\boldsymbol{\varepsilon}}={\mathsf{s}}:{%
\boldsymbol{\sigma}}}}
is
𝗌
=
-
λ
2
μ
(
3
λ
+
2
μ
)
𝐈
⊗
𝐈
+
1
2
μ
𝖨
=
(
1
9
K
-
1
6
G
)
𝐈
⊗
𝐈
+
1
2
G
𝖨
𝗌
tensor-product
𝜆
2
𝜇
3
𝜆
2
𝜇
𝐈
𝐈
1
2
𝜇
𝖨
tensor-product
1
9
𝐾
1
6
𝐺
𝐈
𝐈
1
2
𝐺
𝖨
{\displaystyle{\displaystyle{\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}}}}
In terms of Young's modulus and Poisson's ratio , Hooke's law for isotropic materials can then be expressed as
ε
i
j
=
1
E
(
σ
i
j
-
ν
[
σ
k
k
δ
i
j
-
σ
i
j
]
)
;
𝜺
=
1
E
(
𝝈
-
ν
[
tr
(
𝝈
)
𝐈
-
𝝈
]
)
formulae-sequence
subscript
𝜀
𝑖
𝑗
1
𝐸
subscript
𝜎
𝑖
𝑗
𝜈
delimited-[]
subscript
𝜎
𝑘
𝑘
subscript
𝛿
𝑖
𝑗
subscript
𝜎
𝑖
𝑗
𝜺
1
𝐸
𝝈
𝜈
delimited-[]
tr
𝝈
𝐈
𝝈
{\displaystyle{\displaystyle\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}}])}}
This is the form in which the strain is expressed in terms of the stress tensor in engineering. The expression in expanded form is
ε
11
=
1
E
[
σ
11
-
ν
(
σ
22
+
σ
33
)
]
ε
22
=
1
E
[
σ
22
-
ν
(
σ
11
+
σ
33
)
]
ε
33
=
1
E
[
σ
33
-
ν
(
σ
11
+
σ
22
)
]
ε
12
=
1
2
G
σ
12
;
ε
13
=
1
2
G
σ
13
;
ε
23
=
1
2
G
σ
23
subscript
𝜀
11
absent
1
𝐸
delimited-[]
subscript
𝜎
11
𝜈
subscript
𝜎
22
subscript
𝜎
33
subscript
𝜀
22
absent
1
𝐸
delimited-[]
subscript
𝜎
22
𝜈
subscript
𝜎
11
subscript
𝜎
33
subscript
𝜀
33
absent
1
𝐸
delimited-[]
subscript
𝜎
33
𝜈
subscript
𝜎
11
subscript
𝜎
22
subscript
𝜀
12
formulae-sequence
absent
1
2
𝐺
subscript
𝜎
12
formulae-sequence
subscript
𝜀
13
1
2
𝐺
subscript
𝜎
13
subscript
𝜀
23
1
2
𝐺
subscript
𝜎
23
{\displaystyle{\displaystyle{\begin{aligned} \displaystyle\varepsilon_{11}&%
\displaystyle={\tfrac{1}{E}}\left[\sigma_{11}-\nu(\sigma_{22}+\sigma_{33})%
\right]\\
\displaystyle\varepsilon_{22}&\displaystyle={\tfrac{1}{E}}\left[\sigma_{22}-%
\nu(\sigma_{11}+\sigma_{33})\right]\\
\displaystyle\varepsilon_{33}&\displaystyle={\tfrac{1}{E}}\left[\sigma_{33}-%
\nu(\sigma_{11}+\sigma_{22})\right]\\
\displaystyle\varepsilon_{12}&\displaystyle={\tfrac{1}{2G}}~{}\sigma_{12}~{};~%
{}~{}\varepsilon_{13}={\tfrac{1}{2G}}~{}\sigma_{13}~{};~{}~{}\varepsilon_{23}=%
{\tfrac{1}{2G}}~{}\sigma_{23}\end{aligned}}}}
where E is the Young's modulus and
ν
𝜈
{\displaystyle{\displaystyle\nu}}
is Poisson's ratio . (See 3-D elasticity ).
أظهر 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),
ε
1
′
=
1
E
σ
1
superscript
subscript
𝜀
1
′
1
𝐸
subscript
𝜎
1
{\displaystyle{\displaystyle\varepsilon_{1}^{\prime}={\frac{1}{E}}\sigma_{1}}}
,
ε
2
′
=
-
ν
E
σ
1
superscript
subscript
𝜀
2
′
𝜈
𝐸
subscript
𝜎
1
{\displaystyle{\displaystyle\varepsilon_{2}^{\prime}=-{\frac{\nu}{E}}\sigma_{1%
}}}
,
ε
3
′
=
-
ν
E
σ
1
superscript
subscript
𝜀
3
′
𝜈
𝐸
subscript
𝜎
1
{\displaystyle{\displaystyle\varepsilon_{3}^{\prime}=-{\frac{\nu}{E}}\sigma_{1%
}}}
,
where
ν
𝜈
{\displaystyle{\displaystyle\nu}}
is the Poisson's ratio and
E
𝐸
{\displaystyle{\displaystyle E}}
the Young's modulus .
We get similar equations to the loads in directions 2 and 3,
ε
1
′′
=
-
ν
E
σ
2
superscript
subscript
𝜀
1
′′
𝜈
𝐸
subscript
𝜎
2
{\displaystyle{\displaystyle\varepsilon_{1}^{\prime\prime}=-{\frac{\nu}{E}}%
\sigma_{2}}}
,
ε
2
′′
=
1
E
σ
2
superscript
subscript
𝜀
2
′′
1
𝐸
subscript
𝜎
2
{\displaystyle{\displaystyle\varepsilon_{2}^{\prime\prime}={\frac{1}{E}}\sigma%
_{2}}}
,
ε
3
′′
=
-
ν
E
σ
2
superscript
subscript
𝜀
3
′′
𝜈
𝐸
subscript
𝜎
2
{\displaystyle{\displaystyle\varepsilon_{3}^{\prime\prime}=-{\frac{\nu}{E}}%
\sigma_{2}}}
,
and
ε
1
′′′
=
-
ν
E
σ
3
superscript
subscript
𝜀
1
′′′
𝜈
𝐸
subscript
𝜎
3
{\displaystyle{\displaystyle\varepsilon_{1}^{\prime\prime\prime}=-{\frac{\nu}{%
E}}\sigma_{3}}}
,
ε
2
′′′
=
-
ν
E
σ
3
superscript
subscript
𝜀
2
′′′
𝜈
𝐸
subscript
𝜎
3
{\displaystyle{\displaystyle\varepsilon_{2}^{\prime\prime\prime}=-{\frac{\nu}{%
E}}\sigma_{3}}}
,
ε
3
′′′
=
1
E
σ
3
superscript
subscript
𝜀
3
′′′
1
𝐸
subscript
𝜎
3
{\displaystyle{\displaystyle\varepsilon_{3}^{\prime\prime\prime}={\frac{1}{E}}%
\sigma_{3}}}
.
Summing the three cases together (
ε
i
=
ε
i
′
+
ε
i
′′
+
ε
i
′′′
subscript
𝜀
𝑖
superscript
subscript
𝜀
𝑖
′
superscript
subscript
𝜀
𝑖
′′
superscript
subscript
𝜀
𝑖
′′′
{\displaystyle{\displaystyle\varepsilon_{i}=\varepsilon_{i}^{\prime}+%
\varepsilon_{i}^{\prime\prime}+\varepsilon_{i}^{\prime\prime\prime}}}
) we get
ε
1
=
1
E
(
σ
1
-
ν
(
σ
2
+
σ
3
)
)
subscript
𝜀
1
1
𝐸
subscript
𝜎
1
𝜈
subscript
𝜎
2
subscript
𝜎
3
{\displaystyle{\displaystyle\varepsilon_{1}={\frac{1}{E}}(\sigma_{1}-\nu(%
\sigma_{2}+\sigma_{3}))}}
ε
2
=
1
E
(
σ
2
-
ν
(
σ
1
+
σ
3
)
)
subscript
𝜀
2
1
𝐸
subscript
𝜎
2
𝜈
subscript
𝜎
1
subscript
𝜎
3
{\displaystyle{\displaystyle\varepsilon_{2}={\frac{1}{E}}(\sigma_{2}-\nu(%
\sigma_{1}+\sigma_{3}))}}
ε
3
=
1
E
(
σ
3
-
ν
(
σ
1
+
σ
2
)
)
subscript
𝜀
3
1
𝐸
subscript
𝜎
3
𝜈
subscript
𝜎
1
subscript
𝜎
2
{\displaystyle{\displaystyle\varepsilon_{3}={\frac{1}{E}}(\sigma_{3}-\nu(%
\sigma_{1}+\sigma_{2}))}}
or by adding and subtracting one
ν
σ
𝜈
𝜎
{\displaystyle{\displaystyle\nu\sigma}}
ε
1
=
1
E
(
(
1
+
ν
)
σ
1
-
ν
(
σ
1
+
σ
2
+
σ
3
)
)
subscript
𝜀
1
1
𝐸
1
𝜈
subscript
𝜎
1
𝜈
subscript
𝜎
1
subscript
𝜎
2
subscript
𝜎
3
{\displaystyle{\displaystyle\varepsilon_{1}={\frac{1}{E}}((1+\nu)\sigma_{1}-%
\nu(\sigma_{1}+\sigma_{2}+\sigma_{3}))}}
ε
2
=
1
E
(
(
1
+
ν
)
σ
2
-
ν
(
σ
1
+
σ
2
+
σ
3
)
)
subscript
𝜀
2
1
𝐸
1
𝜈
subscript
𝜎
2
𝜈
subscript
𝜎
1
subscript
𝜎
2
subscript
𝜎
3
{\displaystyle{\displaystyle\varepsilon_{2}={\frac{1}{E}}((1+\nu)\sigma_{2}-%
\nu(\sigma_{1}+\sigma_{2}+\sigma_{3}))}}
ε
3
=
1
E
(
(
1
+
ν
)
σ
3
-
ν
(
σ
1
+
σ
2
+
σ
3
)
)
subscript
𝜀
3
1
𝐸
1
𝜈
subscript
𝜎
3
𝜈
subscript
𝜎
1
subscript
𝜎
2
subscript
𝜎
3
{\displaystyle{\displaystyle\varepsilon_{3}={\frac{1}{E}}((1+\nu)\sigma_{3}-%
\nu(\sigma_{1}+\sigma_{2}+\sigma_{3}))}}
and further we get by solving
σ
1
subscript
𝜎
1
{\displaystyle{\displaystyle\sigma_{1}}}
σ
1
=
E
1
+
ν
ε
1
+
ν
1
+
ν
(
σ
1
+
σ
2
+
σ
3
)
subscript
𝜎
1
𝐸
1
𝜈
subscript
𝜀
1
𝜈
1
𝜈
subscript
𝜎
1
subscript
𝜎
2
subscript
𝜎
3
{\displaystyle{\displaystyle\sigma_{1}={\frac{E}{1+\nu}}\varepsilon_{1}+{\frac%
{\nu}{1+\nu}}(\sigma_{1}+\sigma_{2}+\sigma_{3})}}
.
Calculating the sum
∑
i
=
1
,
2
,
3
ε
i
=
1
E
(
(
1
+
ν
)
∑
i
=
1
,
2
,
3
σ
i
-
3
ν
(
∑
i
=
1
,
2
,
3
σ
i
)
)
=
1
-
2
ν
E
∑
i
=
1
,
2
,
3
σ
i
subscript
𝑖
1
2
3
subscript
𝜀
𝑖
1
𝐸
1
𝜈
subscript
𝑖
1
2
3
subscript
𝜎
𝑖
3
𝜈
subscript
𝑖
1
2
3
subscript
𝜎
𝑖
1
2
𝜈
𝐸
subscript
𝑖
1
2
3
subscript
𝜎
𝑖
{\displaystyle{\displaystyle\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}}}
σ
1
+
σ
2
+
σ
3
=
E
1
-
2
ν
(
ε
1
+
ε
2
+
ε
3
)
subscript
𝜎
1
subscript
𝜎
2
subscript
𝜎
3
𝐸
1
2
𝜈
subscript
𝜀
1
subscript
𝜀
2
subscript
𝜀
3
{\displaystyle{\displaystyle\sigma_{1}+\sigma_{2}+\sigma_{3}={\frac{E}{1-2\nu}%
}(\varepsilon_{1}+\varepsilon_{2}+\varepsilon_{3})}}
and substituting it to the equation solved for
σ
1
subscript
𝜎
1
{\displaystyle{\displaystyle\sigma_{1}}}
gives
σ
1
=
E
1
+
ν
ε
1
+
E
ν
(
1
+
ν
)
(
1
-
2
ν
)
(
ε
1
+
ε
2
+
ε
3
)
subscript
𝜎
1
𝐸
1
𝜈
subscript
𝜀
1
𝐸
𝜈
1
𝜈
1
2
𝜈
subscript
𝜀
1
subscript
𝜀
2
subscript
𝜀
3
{\displaystyle{\displaystyle\sigma_{1}={\frac{E}{1+\nu}}\varepsilon_{1}+{\frac%
{E\nu}{(1+\nu)(1-2\nu)}}(\varepsilon_{1}+\varepsilon_{2}+\varepsilon_{3})}}
,
σ
1
=
2
μ
ε
1
+
λ
(
ε
1
+
ε
2
+
ε
3
)
subscript
𝜎
1
2
𝜇
subscript
𝜀
1
𝜆
subscript
𝜀
1
subscript
𝜀
2
subscript
𝜀
3
{\displaystyle{\displaystyle\sigma_{1}=2\mu\varepsilon_{1}+\lambda(\varepsilon%
_{1}+\varepsilon_{2}+\varepsilon_{3})}}
,
where
μ
𝜇
{\displaystyle{\displaystyle\mu}}
and
λ
𝜆
{\displaystyle{\displaystyle\lambda}}
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
[
ε
11
ε
22
ε
33
2
ε
23
2
ε
13
2
ε
12
]
=
[
ε
11
ε
22
ε
33
γ
23
γ
13
γ
12
]
=
1
E
[
1
-
ν
-
ν
0
0
0
-
ν
1
-
ν
0
0
0
-
ν
-
ν
1
0
0
0
0
0
0
2
(
1
+
ν
)
0
0
0
0
0
0
2
(
1
+
ν
)
0
0
0
0
0
0
2
(
1
+
ν
)
]
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
subscript
𝜀
11
subscript
𝜀
22
subscript
𝜀
33
2
subscript
𝜀
23
2
subscript
𝜀
13
2
subscript
𝜀
12
subscript
𝜀
11
subscript
𝜀
22
subscript
𝜀
33
subscript
𝛾
23
subscript
𝛾
13
subscript
𝛾
12
continued-fraction
1
𝐸
1
𝜈
𝜈
0
0
0
𝜈
1
𝜈
0
0
0
𝜈
𝜈
1
0
0
0
0
0
0
2
1
𝜈
0
0
0
0
0
0
2
1
𝜈
0
0
0
0
0
0
2
1
𝜈
subscript
𝜎
11
subscript
𝜎
22
subscript
𝜎
33
subscript
𝜎
23
subscript
𝜎
13
subscript
𝜎
12
{\displaystyle{\displaystyle{\begin{bmatrix}\varepsilon_{11}\\
\varepsilon_{22}\\
\varepsilon_{33}\\
2\varepsilon_{23}\\
2\varepsilon_{13}\\
2\varepsilon_{12}\end{bmatrix}}={\begin{bmatrix}\varepsilon_{11}\\
\varepsilon_{22}\\
\varepsilon_{33}\\
\gamma_{23}\\
\gamma_{13}\\
\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_{13}\\
\sigma_{12}\end{bmatrix}}}}
where
γ
i
j
:=
2
ε
i
j
assign
subscript
𝛾
𝑖
𝑗
2
subscript
𝜀
𝑖
𝑗
{\displaystyle{\displaystyle\gamma_{ij}:=2\varepsilon_{ij}}}
is the engineering shear strain .
The inverse relation may be written as
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
=
E
(
1
+
ν
)
(
1
-
2
ν
)
[
1
-
ν
ν
ν
0
0
0
ν
1
-
ν
ν
0
0
0
ν
ν
1
-
ν
0
0
0
0
0
0
(
1
-
2
ν
)
/
2
0
0
0
0
0
0
(
1
-
2
ν
)
/
2
0
0
0
0
0
0
(
1
-
2
ν
)
/
2
]
[
ε
11
ε
22
ε
33
2
ε
23
2
ε
13
2
ε
12
]
subscript
𝜎
11
subscript
𝜎
22
subscript
𝜎
33
subscript
𝜎
23
subscript
𝜎
13
subscript
𝜎
12
continued-fraction
𝐸
1
𝜈
1
2
𝜈
1
𝜈
𝜈
𝜈
0
0
0
𝜈
1
𝜈
𝜈
0
0
0
𝜈
𝜈
1
𝜈
0
0
0
0
0
0
1
2
𝜈
2
0
0
0
0
0
0
1
2
𝜈
2
0
0
0
0
0
0
1
2
𝜈
2
subscript
𝜀
11
subscript
𝜀
22
subscript
𝜀
33
2
subscript
𝜀
23
2
subscript
𝜀
13
2
subscript
𝜀
12
{\displaystyle{\displaystyle{\begin{bmatrix}\sigma_{11}\\
\sigma_{22}\\
\sigma_{33}\\
\sigma_{23}\\
\sigma_{13}\\
\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_{13}\\
2\varepsilon_{12}\end{bmatrix}}}}
which can be simplified thanks to the Lamé constants :
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
=
[
2
μ
+
λ
λ
λ
0
0
0
λ
2
μ
+
λ
λ
0
0
0
λ
λ
2
μ
+
λ
0
0
0
0
0
0
μ
0
0
0
0
0
0
μ
0
0
0
0
0
0
μ
]
[
ε
11
ε
22
ε
33
2
ε
23
2
ε
13
2
ε
12
]
subscript
𝜎
11
subscript
𝜎
22
subscript
𝜎
33
subscript
𝜎
23
subscript
𝜎
13
subscript
𝜎
12
2
𝜇
𝜆
𝜆
𝜆
0
0
0
𝜆
2
𝜇
𝜆
𝜆
0
0
0
𝜆
𝜆
2
𝜇
𝜆
0
0
0
0
0
0
𝜇
0
0
0
0
0
0
𝜇
0
0
0
0
0
0
𝜇
subscript
𝜀
11
subscript
𝜀
22
subscript
𝜀
33
2
subscript
𝜀
23
2
subscript
𝜀
13
2
subscript
𝜀
12
{\displaystyle{\displaystyle{\begin{bmatrix}\sigma_{11}\\
\sigma_{22}\\
\sigma_{33}\\
\sigma_{23}\\
\sigma_{13}\\
\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_{13}\\
2\varepsilon_{12}\end{bmatrix}}}}
Plane stress
Under plane stress conditions
σ
31
=
σ
13
=
σ
32
=
σ
23
=
σ
33
=
0
subscript
𝜎
31
subscript
𝜎
13
subscript
𝜎
32
subscript
𝜎
23
subscript
𝜎
33
0
{\displaystyle{\displaystyle\sigma_{31}=\sigma_{13}=\sigma_{32}=\sigma_{23}=%
\sigma_{33}=0}}
. In that case Hooke's law takes the form
[
ε
11
ε
22
2
ε
12
]
=
1
E
[
1
-
ν
0
-
ν
1
0
0
0
2
(
1
+
ν
)
]
[
σ
11
σ
22
σ
12
]
subscript
𝜀
11
subscript
𝜀
22
2
subscript
𝜀
12
continued-fraction
1
𝐸
1
𝜈
0
𝜈
1
0
0
0
2
1
𝜈
subscript
𝜎
11
subscript
𝜎
22
subscript
𝜎
12
{\displaystyle{\displaystyle{\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}}}}
The inverse relation is usually written in the reduced form
[
σ
11
σ
22
σ
12
]
=
E
1
-
ν
2
[
1
ν
0
ν
1
0
0
0
1
-
ν
2
]
[
ε
11
ε
22
2
ε
12
]
subscript
𝜎
11
subscript
𝜎
22
subscript
𝜎
12
continued-fraction
𝐸
1
superscript
𝜈
2
1
𝜈
0
𝜈
1
0
0
0
continued-fraction
1
𝜈
2
subscript
𝜀
11
subscript
𝜀
22
2
subscript
𝜀
12
{\displaystyle{\displaystyle{\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}}}}
Anisotropic materials
σ
i
j
=
∂
U
∂
ϵ
i
j
⟹
c
i
j
k
ℓ
=
∂
2
U
∂
ϵ
i
j
∂
ϵ
k
ℓ
.
formulae-sequence
subscript
𝜎
𝑖
𝑗
continued-fraction
𝑈
subscript
italic-ϵ
𝑖
𝑗
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
continued-fraction
superscript
2
𝑈
subscript
italic-ϵ
𝑖
𝑗
subscript
italic-ϵ
𝑘
ℓ
{\displaystyle{\displaystyle\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}}}~{}.}}
The arbitrariness of the order of differentiation implies that
c
i
j
k
ℓ
=
c
k
ℓ
i
j
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
subscript
𝑐
𝑘
ℓ
𝑖
𝑗
{\displaystyle{\displaystyle c_{ijk\ell}=c_{k\ell ij}\,}}
. 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 (
𝐞
1
,
𝐞
2
,
𝐞
3
subscript
𝐞
1
subscript
𝐞
2
subscript
𝐞
3
{\displaystyle{\displaystyle\mathbf{e}_{1},\mathbf{e}_{2},\mathbf{e}_{3}}}
) as
[
𝝈
]
=
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
≡
[
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
]
;
[
ϵ
]
=
[
ϵ
11
ϵ
22
ϵ
33
2
ϵ
23
2
ϵ
13
2
ϵ
12
]
≡
[
ϵ
1
ϵ
2
ϵ
3
ϵ
4
ϵ
5
ϵ
6
]
formulae-sequence
delimited-[]
𝝈
subscript
𝜎
11
subscript
𝜎
22
subscript
𝜎
33
subscript
𝜎
23
subscript
𝜎
13
subscript
𝜎
12
subscript
𝜎
1
subscript
𝜎
2
subscript
𝜎
3
subscript
𝜎
4
subscript
𝜎
5
subscript
𝜎
6
delimited-[]
bold-italic-ϵ
subscript
italic-ϵ
11
subscript
italic-ϵ
22
subscript
italic-ϵ
33
2
subscript
italic-ϵ
23
2
subscript
italic-ϵ
13
2
subscript
italic-ϵ
12
subscript
italic-ϵ
1
subscript
italic-ϵ
2
subscript
italic-ϵ
3
subscript
italic-ϵ
4
subscript
italic-ϵ
5
subscript
italic-ϵ
6
{\displaystyle{\displaystyle[{\boldsymbol{\sigma}}]={\begin{bmatrix}\sigma_{11%
}\\
\sigma_{22}\\
\sigma_{33}\\
\sigma_{23}\\
\sigma_{13}\\
\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_{13}\\
2\epsilon_{12}\end{bmatrix}}\equiv{\begin{bmatrix}\epsilon_{1}\\
\epsilon_{2}\\
\epsilon_{3}\\
\epsilon_{4}\\
\epsilon_{5}\\
\epsilon_{6}\end{bmatrix}}}}
Then the stiffness tensor (
𝖼
𝖼
{\displaystyle{\displaystyle{\mathsf{c}}}}
) can be expressed as
[
𝖢
]
=
[
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
]
≡
[
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
]
delimited-[]
𝖢
subscript
𝑐
1111
subscript
𝑐
1122
subscript
𝑐
1133
subscript
𝑐
1123
subscript
𝑐
1131
subscript
𝑐
1112
subscript
𝑐
2211
subscript
𝑐
2222
subscript
𝑐
2233
subscript
𝑐
2223
subscript
𝑐
2231
subscript
𝑐
2212
subscript
𝑐
3311
subscript
𝑐
3322
subscript
𝑐
3333
subscript
𝑐
3323
subscript
𝑐
3331
subscript
𝑐
3312
subscript
𝑐
2311
subscript
𝑐
2322
subscript
𝑐
2333
subscript
𝑐
2323
subscript
𝑐
2331
subscript
𝑐
2312
subscript
𝑐
3111
subscript
𝑐
3122
subscript
𝑐
3133
subscript
𝑐
3123
subscript
𝑐
3131
subscript
𝑐
3112
subscript
𝑐
1211
subscript
𝑐
1222
subscript
𝑐
1233
subscript
𝑐
1223
subscript
𝑐
1231
subscript
𝑐
1212
subscript
𝐶
11
subscript
𝐶
12
subscript
𝐶
13
subscript
𝐶
14
subscript
𝐶
15
subscript
𝐶
16
subscript
𝐶
12
subscript
𝐶
22
subscript
𝐶
23
subscript
𝐶
24
subscript
𝐶
25
subscript
𝐶
26
subscript
𝐶
13
subscript
𝐶
23
subscript
𝐶
33
subscript
𝐶
34
subscript
𝐶
35
subscript
𝐶
36
subscript
𝐶
14
subscript
𝐶
24
subscript
𝐶
34
subscript
𝐶
44
subscript
𝐶
45
subscript
𝐶
46
subscript
𝐶
15
subscript
𝐶
25
subscript
𝐶
35
subscript
𝐶
45
subscript
𝐶
55
subscript
𝐶
56
subscript
𝐶
16
subscript
𝐶
26
subscript
𝐶
36
subscript
𝐶
46
subscript
𝐶
56
subscript
𝐶
66
{\displaystyle{\displaystyle[{\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}}}}
and Hooke's law is written as
[
𝝈
]
=
[
𝖢
]
[
ϵ
]
or
σ
i
=
C
i
j
ϵ
j
.
formulae-sequence
delimited-[]
𝝈
delimited-[]
𝖢
delimited-[]
bold-italic-ϵ
or
subscript
𝜎
𝑖
subscript
𝐶
𝑖
𝑗
subscript
italic-ϵ
𝑗
{\displaystyle{\displaystyle[{\boldsymbol{\sigma}}]=[{\mathsf{C}}][{%
\boldsymbol{\epsilon}}]\qquad{\text{or}}\qquad\sigma_{i}=C_{ij}\epsilon_{j}~{}%
.}}
Similarly the compliance tensor (
𝗌
𝗌
{\displaystyle{\displaystyle{\mathsf{s}}}}
) can be written as
[
𝖲
]
=
[
s
1111
s
1122
s
1133
2
s
1123
2
s
1131
2
s
1112
s
2211
s
2222
s
2233
2
s
2223
2
s
2231
2
s
2212
s
3311
s
3322
s
3333
2
s
3323
2
s
3331
2
s
3312
2
s
2311
2
s
2322
2
s
2333
4
s
2323
4
s
2331
4
s
2312
2
s
3111
2
s
3122
2
s
3133
4
s
3123
4
s
3131
4
s
3112
2
s
1211
2
s
1222
2
s
1233
4
s
1223
4
s
1231
4
s
1212
]
≡
[
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
]
delimited-[]
𝖲
subscript
𝑠
1111
subscript
𝑠
1122
subscript
𝑠
1133
2
subscript
𝑠
1123
2
subscript
𝑠
1131
2
subscript
𝑠
1112
subscript
𝑠
2211
subscript
𝑠
2222
subscript
𝑠
2233
2
subscript
𝑠
2223
2
subscript
𝑠
2231
2
subscript
𝑠
2212
subscript
𝑠
3311
subscript
𝑠
3322
subscript
𝑠
3333
2
subscript
𝑠
3323
2
subscript
𝑠
3331
2
subscript
𝑠
3312
2
subscript
𝑠
2311
2
subscript
𝑠
2322
2
subscript
𝑠
2333
4
subscript
𝑠
2323
4
subscript
𝑠
2331
4
subscript
𝑠
2312
2
subscript
𝑠
3111
2
subscript
𝑠
3122
2
subscript
𝑠
3133
4
subscript
𝑠
3123
4
subscript
𝑠
3131
4
subscript
𝑠
3112
2
subscript
𝑠
1211
2
subscript
𝑠
1222
2
subscript
𝑠
1233
4
subscript
𝑠
1223
4
subscript
𝑠
1231
4
subscript
𝑠
1212
subscript
𝑆
11
subscript
𝑆
12
subscript
𝑆
13
subscript
𝑆
14
subscript
𝑆
15
subscript
𝑆
16
subscript
𝑆
12
subscript
𝑆
22
subscript
𝑆
23
subscript
𝑆
24
subscript
𝑆
25
subscript
𝑆
26
subscript
𝑆
13
subscript
𝑆
23
subscript
𝑆
33
subscript
𝑆
34
subscript
𝑆
35
subscript
𝑆
36
subscript
𝑆
14
subscript
𝑆
24
subscript
𝑆
34
subscript
𝑆
44
subscript
𝑆
45
subscript
𝑆
46
subscript
𝑆
15
subscript
𝑆
25
subscript
𝑆
35
subscript
𝑆
45
subscript
𝑆
55
subscript
𝑆
56
subscript
𝑆
16
subscript
𝑆
26
subscript
𝑆
36
subscript
𝑆
46
subscript
𝑆
56
subscript
𝑆
66
{\displaystyle{\displaystyle[{\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}}}}
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[4]
c
p
q
r
s
=
l
p
i
l
q
j
l
r
k
l
s
ℓ
c
i
j
k
ℓ
subscript
𝑐
𝑝
𝑞
𝑟
𝑠
subscript
𝑙
𝑝
𝑖
subscript
𝑙
𝑞
𝑗
subscript
𝑙
𝑟
𝑘
subscript
𝑙
𝑠
ℓ
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
{\displaystyle{\displaystyle c_{pqrs}=l_{pi}~{}l_{qj}~{}l_{rk}~{}l_{s\ell}~{}c%
_{ijk\ell}}}
where
l
a
b
subscript
𝑙
𝑎
𝑏
{\displaystyle{\displaystyle l_{ab}}}
are the components of an orthogonal rotation matrix
[
L
]
delimited-[]
𝐿
{\displaystyle{\displaystyle[L]}}
. The same relation also holds for inversions.
In matrix notation, if the transformed basis (rotated or inverted) is related to the reference basis by
[
𝐞
i
′
]
=
[
L
]
[
𝐞
i
]
delimited-[]
superscript
subscript
𝐞
𝑖
′
delimited-[]
𝐿
delimited-[]
subscript
𝐞
𝑖
{\displaystyle{\displaystyle[\mathbf{e}_{i}^{\prime}]=[L][\mathbf{e}_{i}]}}
then
C
i
j
ϵ
i
ϵ
j
=
C
i
j
′
ϵ
i
′
ϵ
j
′
.
subscript
𝐶
𝑖
𝑗
subscript
italic-ϵ
𝑖
subscript
italic-ϵ
𝑗
superscript
subscript
𝐶
𝑖
𝑗
′
subscript
superscript
italic-ϵ
′
𝑖
subscript
superscript
italic-ϵ
′
𝑗
{\displaystyle{\displaystyle C_{ij}~{}\epsilon_{i}~{}\epsilon_{j}=C_{ij}^{%
\prime}~{}\epsilon^{\prime}_{i}~{}\epsilon^{\prime}_{j}~{}.}}
In addition, if the material is symmetric with respect to the transformation
[
L
]
delimited-[]
𝐿
{\displaystyle{\displaystyle[L]}}
then
C
i
j
=
C
i
j
′
⟹
C
i
j
(
ϵ
i
ϵ
j
-
ϵ
i
′
ϵ
j
′
)
=
0
.
formulae-sequence
subscript
𝐶
𝑖
𝑗
subscript
superscript
𝐶
′
𝑖
𝑗
subscript
𝐶
𝑖
𝑗
subscript
italic-ϵ
𝑖
subscript
italic-ϵ
𝑗
subscript
superscript
italic-ϵ
′
𝑖
subscript
superscript
italic-ϵ
′
𝑗
0
{\displaystyle{\displaystyle C_{ij}=C^{\prime}_{ij}\quad\implies\quad C_{ij}~{%
}(\epsilon_{i}~{}\epsilon_{j}-\epsilon^{\prime}_{i}~{}\epsilon^{\prime}_{j})=0%
~{}.}}
Orthotropic materials
Orthotropic materials have three orthogonal planes of symmetry . If the basis vectors (
𝐞
1
,
𝐞
2
,
𝐞
3
subscript
𝐞
1
subscript
𝐞
2
subscript
𝐞
3
{\displaystyle{\displaystyle\mathbf{e}_{1},\mathbf{e}_{2},\mathbf{e}_{3}}}
) are normals to the planes of symmetry then the coordinate transformation relations imply that
[
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
]
=
[
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
]
[
ϵ
1
ϵ
2
ϵ
3
ϵ
4
ϵ
5
ϵ
6
]
subscript
𝜎
1
subscript
𝜎
2
subscript
𝜎
3
subscript
𝜎
4
subscript
𝜎
5
subscript
𝜎
6
subscript
𝐶
11
subscript
𝐶
12
subscript
𝐶
13
0
0
0
subscript
𝐶
12
subscript
𝐶
22
subscript
𝐶
23
0
0
0
subscript
𝐶
13
subscript
𝐶
23
subscript
𝐶
33
0
0
0
0
0
0
subscript
𝐶
44
0
0
0
0
0
0
subscript
𝐶
55
0
0
0
0
0
0
subscript
𝐶
66
subscript
italic-ϵ
1
subscript
italic-ϵ
2
subscript
italic-ϵ
3
subscript
italic-ϵ
4
subscript
italic-ϵ
5
subscript
italic-ϵ
6
{\displaystyle{\displaystyle{\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}}}}
The inverse of this relation is commonly written as[5]
[
ϵ
xx
ϵ
yy
ϵ
zz
2
ϵ
yz
2
ϵ
zx
2
ϵ
xy
]
=
[
1
E
x
-
ν
yx
E
y
-
ν
zx
E
z
0
0
0
-
ν
xy
E
x
1
E
y
-
ν
zy
E
z
0
0
0
-
ν
xz
E
x
-
ν
yz
E
y
1
E
z
0
0
0
0
0
0
1
G
yz
0
0
0
0
0
0
1
G
zx
0
0
0
0
0
0
1
G
xy
]
[
σ
xx
σ
yy
σ
zz
σ
yz
σ
zx
σ
xy
]
subscript
italic-ϵ
xx
subscript
italic-ϵ
yy
subscript
italic-ϵ
zz
2
subscript
italic-ϵ
yz
2
subscript
italic-ϵ
zx
2
subscript
italic-ϵ
xy
1
subscript
𝐸
x
subscript
𝜈
yx
subscript
𝐸
y
subscript
𝜈
zx
subscript
𝐸
z
0
0
0
subscript
𝜈
xy
subscript
𝐸
x
1
subscript
𝐸
y
subscript
𝜈
zy
subscript
𝐸
z
0
0
0
subscript
𝜈
xz
subscript
𝐸
x
subscript
𝜈
yz
subscript
𝐸
y
1
subscript
𝐸
z
0
0
0
0
0
0
1
subscript
𝐺
yz
0
0
0
0
0
0
1
subscript
𝐺
zx
0
0
0
0
0
0
1
subscript
𝐺
xy
subscript
𝜎
xx
subscript
𝜎
yy
subscript
𝜎
zz
subscript
𝜎
yz
subscript
𝜎
zx
subscript
𝜎
xy
{\displaystyle{\displaystyle{\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}}}}
where
E
i
subscript
𝐸
i
{\displaystyle{\displaystyle{E}_{\rm{i}}\,}}
is the Young's modulus along axis
i
𝑖
{\displaystyle{\displaystyle i}}
G
ij
subscript
𝐺
ij
{\displaystyle{\displaystyle G_{\rm{ij}}\,}}
is the shear modulus in direction
j
𝑗
{\displaystyle{\displaystyle j}}
on the plane whose normal is in direction
i
𝑖
{\displaystyle{\displaystyle i}}
ν
ij
subscript
𝜈
ij
{\displaystyle{\displaystyle\nu_{\rm{ij}}\,}}
is the Poisson's ratio that corresponds to a contraction in direction
j
𝑗
{\displaystyle{\displaystyle j}}
when an extension is applied in direction
i
𝑖
{\displaystyle{\displaystyle i}}
.
Under plane stress conditions,
σ
z
z
=
σ
z
x
=
σ
y
z
=
0
subscript
𝜎
𝑧
𝑧
subscript
𝜎
𝑧
𝑥
subscript
𝜎
𝑦
𝑧
0
{\displaystyle{\displaystyle\sigma_{zz}=\sigma_{zx}=\sigma_{yz}=0}}
, Hooke's law for an orthotropic material takes the form
[
ε
xx
ε
yy
2
ε
xy
]
=
[
1
E
x
-
ν
yx
E
y
0
-
ν
xy
E
x
1
E
y
0
0
0
1
G
xy
]
[
σ
xx
σ
yy
σ
xy
]
.
subscript
𝜀
xx
subscript
𝜀
yy
2
subscript
𝜀
xy
1
subscript
𝐸
x
subscript
𝜈
yx
subscript
𝐸
y
0
subscript
𝜈
xy
subscript
𝐸
x
1
subscript
𝐸
y
0
0
0
1
subscript
𝐺
xy
subscript
𝜎
xx
subscript
𝜎
yy
subscript
𝜎
xy
{\displaystyle{\displaystyle{\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}}\,.}}
The inverse relation is
[
σ
xx
σ
yy
σ
xy
]
=
1
1
-
ν
xy
ν
yx
[
E
x
ν
yx
E
x
0
ν
xy
E
y
E
y
0
0
0
G
xy
(
1
-
ν
xy
ν
yx
)
]
[
ε
xx
ε
yy
2
ε
xy
]
.
subscript
𝜎
xx
subscript
𝜎
yy
subscript
𝜎
xy
continued-fraction
1
1
subscript
𝜈
xy
subscript
𝜈
yx
subscript
𝐸
x
subscript
𝜈
yx
subscript
𝐸
x
0
subscript
𝜈
xy
subscript
𝐸
y
subscript
𝐸
y
0
0
0
subscript
𝐺
xy
1
subscript
𝜈
xy
subscript
𝜈
yx
subscript
𝜀
xx
subscript
𝜀
yy
2
subscript
𝜀
xy
{\displaystyle{\displaystyle{\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}}\,.}}
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
𝐞
3
subscript
𝐞
3
{\displaystyle{\displaystyle\mathbf{e}_{3}}}
is the axis of symmetry, Hooke's law can be expressed as
[
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
]
=
[
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
1
2
(
C
11
-
C
12
)
]
[
ϵ
1
ϵ
2
ϵ
3
ϵ
4
ϵ
5
ϵ
6
]
subscript
𝜎
1
subscript
𝜎
2
subscript
𝜎
3
subscript
𝜎
4
subscript
𝜎
5
subscript
𝜎
6
subscript
𝐶
11
subscript
𝐶
12
subscript
𝐶
13
0
0
0
subscript
𝐶
12
subscript
𝐶
11
subscript
𝐶
13
0
0
0
subscript
𝐶
13
subscript
𝐶
13
subscript
𝐶
33
0
0
0
0
0
0
subscript
𝐶
44
0
0
0
0
0
0
subscript
𝐶
44
0
0
0
0
0
0
1
2
subscript
𝐶
11
subscript
𝐶
12
subscript
italic-ϵ
1
subscript
italic-ϵ
2
subscript
italic-ϵ
3
subscript
italic-ϵ
4
subscript
italic-ϵ
5
subscript
italic-ϵ
6
{\displaystyle{\displaystyle{\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}}}}
More frequently, the
x
≡
𝐞
1
𝑥
subscript
𝐞
1
{\displaystyle{\displaystyle x\equiv\mathbf{e}_{1}}}
axis is taken to be the axis of symmetry and the inverse Hooke's law is written as
[6]
[
ϵ
xx
ϵ
yy
ϵ
zz
2
ϵ
yz
2
ϵ
zx
2
ϵ
xy
]
=
[
1
E
x
-
ν
yx
E
y
-
ν
yx
E
y
0
0
0
-
ν
xy
E
x
1
E
y
-
ν
yz
E
y
0
0
0
-
ν
xy
E
x
-
ν
yz
E
y
1
E
y
0
0
0
0
0
0
2
(
1
+
ν
yz
)
E
y
0
0
0
0
0
0
1
G
xy
0
0
0
0
0
0
1
G
xy
]
[
σ
xx
σ
yy
σ
zz
σ
yz
σ
zx
σ
xy
]
subscript
italic-ϵ
xx
subscript
italic-ϵ
yy
subscript
italic-ϵ
zz
2
subscript
italic-ϵ
yz
2
subscript
italic-ϵ
zx
2
subscript
italic-ϵ
xy
1
subscript
𝐸
x
subscript
𝜈
yx
subscript
𝐸
y
subscript
𝜈
yx
subscript
𝐸
y
0
0
0
subscript
𝜈
xy
subscript
𝐸
x
1
subscript
𝐸
y
subscript
𝜈
yz
subscript
𝐸
y
0
0
0
subscript
𝜈
xy
subscript
𝐸
x
subscript
𝜈
yz
subscript
𝐸
y
1
subscript
𝐸
y
0
0
0
0
0
0
2
1
subscript
𝜈
yz
subscript
𝐸
y
0
0
0
0
0
0
1
subscript
𝐺
xy
0
0
0
0
0
0
1
subscript
𝐺
xy
subscript
𝜎
xx
subscript
𝜎
yy
subscript
𝜎
zz
subscript
𝜎
yz
subscript
𝜎
zx
subscript
𝜎
xy
{\displaystyle{\displaystyle{\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}}}}
Thermodynamic basis
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
δ
W
=
δ
U
𝛿
𝑊
𝛿
𝑈
{\displaystyle{\displaystyle\delta W=\delta U\,}}
where
δ
U
𝛿
𝑈
{\displaystyle{\displaystyle\delta U}}
is the increase in internal energy and
δ
W
𝛿
𝑊
{\displaystyle{\displaystyle\delta W}}
is the work done by external forces. The work can be split into two terms
δ
W
=
δ
W
s
+
δ
W
b
𝛿
𝑊
𝛿
subscript
𝑊
𝑠
𝛿
subscript
𝑊
𝑏
{\displaystyle{\displaystyle\delta W=\delta W_{s}+\delta W_{b}\,}}
where
δ
W
s
𝛿
subscript
𝑊
𝑠
{\displaystyle{\displaystyle\delta W_{s}}}
is the work done by surface forces while
δ
W
b
𝛿
subscript
𝑊
𝑏
{\displaystyle{\displaystyle\delta W_{b}}}
is the work done by body forces . If
δ
𝐮
𝛿
𝐮
{\displaystyle{\displaystyle\delta\mathbf{u}}}
is a variation of the displacement field
𝐮
𝐮
{\displaystyle{\displaystyle\mathbf{u}}}
in the body, then the two external work terms can be expressed as
δ
W
s
=
∫
∂
Ω
𝐭
⋅
δ
𝐮
dS
;
δ
W
b
=
∫
Ω
𝐛
⋅
δ
𝐮
dV
formulae-sequence
𝛿
subscript
𝑊
𝑠
subscript
Ω
⋅
𝐭
𝛿
𝐮
dS
𝛿
subscript
𝑊
𝑏
subscript
Ω
⋅
𝐛
𝛿
𝐮
dV
{\displaystyle{\displaystyle\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}}}}
where
𝐭
𝐭
{\displaystyle{\displaystyle\mathbf{t}}}
is the surface traction vector,
𝐛
𝐛
{\displaystyle{\displaystyle\mathbf{b}}}
is the body force vector,
Ω
Ω
{\displaystyle{\displaystyle\Omega\,}}
represents the body and
∂
Ω
Ω
{\displaystyle{\displaystyle\partial\Omega}}
represents its surface. Using the relation between the Cauchy stress and the surface traction,
𝐭
=
𝐧
⋅
𝝈
𝐭
⋅
𝐧
𝝈
{\displaystyle{\displaystyle\mathbf{t}=\mathbf{n}\cdot{\boldsymbol{\sigma}}}}
(where
𝐧
𝐧
{\displaystyle{\displaystyle\mathbf{n}}}
is the unit outward normal to
∂
Ω
Ω
{\displaystyle{\displaystyle\partial\Omega}}
), we have
δ
W
=
δ
U
=
∫
∂
Ω
(
𝐧
⋅
𝝈
)
⋅
δ
𝐮
dS
+
∫
Ω
𝐛
⋅
δ
𝐮
dV
𝛿
𝑊
𝛿
𝑈
subscript
Ω
⋅
⋅
𝐧
𝝈
𝛿
𝐮
dS
subscript
Ω
⋅
𝐛
𝛿
𝐮
dV
{\displaystyle{\displaystyle\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}}}}
Converting the surface integral into a volume integral via the divergence theorem gives
δ
U
=
∫
Ω
[
∇
⋅
(
𝝈
⋅
δ
𝐮
)
+
𝐛
⋅
δ
𝐮
]
dV
.
𝛿
𝑈
subscript
Ω
delimited-[]
⋅
bold-∇
⋅
𝝈
𝛿
𝐮
⋅
𝐛
𝛿
𝐮
dV
{\displaystyle{\displaystyle\delta U=\int_{\Omega}[{\boldsymbol{\nabla}}\cdot(%
{\boldsymbol{\sigma}}\cdot\delta\mathbf{u})+\mathbf{b}\cdot\delta\mathbf{u}]~{%
}{\rm{dV}}~{}.}}
Using the symmetry of the Cauchy stress and the identity
∇
⋅
(
𝑨
⋅
𝐛
)
=
(
∇
⋅
𝑨
)
⋅
𝐛
+
1
2
[
𝑨
T
:
∇
𝐛
+
𝑨
:
(
∇
𝐛
)
T
]
fragments
bold-∇
⋅
fragments
(
A
⋅
b
)
fragments
(
bold-∇
⋅
A
)
⋅
b
1
2
fragments
[
superscript
𝑨
𝑇
:
bold-∇
b
A
:
superscript
fragments
(
bold-∇
b
)
𝑇
]
{\displaystyle{\displaystyle{\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}]}}
we have the following
δ
U
=
∫
Ω
[
𝝈
:
1
2
{
∇
δ
𝐮
+
(
∇
δ
𝐮
)
T
}
+
{
∇
⋅
𝝈
+
𝐛
}
⋅
δ
𝐮
]
dV
.
fragments
δ
U
subscript
Ω
fragments
[
σ
:
1
2
fragments
{
bold-∇
δ
u
superscript
fragments
(
bold-∇
δ
u
)
𝑇
}
fragments
{
bold-∇
⋅
σ
b
}
⋅
δ
u
]
dV
.
{\displaystyle{\displaystyle\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}}~{}.}}
From the definition of strain and from the equations of equilibrium we have
δ
ϵ
=
1
2
[
∇
δ
𝐮
+
(
∇
δ
𝐮
)
T
]
;
∇
⋅
𝝈
+
𝐛
=
𝟎
.
formulae-sequence
𝛿
bold-italic-ϵ
1
2
delimited-[]
bold-∇
𝛿
𝐮
superscript
bold-∇
𝛿
𝐮
𝑇
⋅
bold-∇
𝝈
𝐛
0
{\displaystyle{\displaystyle\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}~{}.}}
Hence we can write
δ
U
=
∫
Ω
𝝈
:
δ
ϵ
dV
:
𝛿
𝑈
subscript
Ω
𝝈
𝛿
bold-italic-ϵ
dV
{\displaystyle{\displaystyle\delta U=\int_{\Omega}{\boldsymbol{\sigma}}:\delta%
{\boldsymbol{\epsilon}}~{}{\rm{dV}}}}
and therefore the variation in the internal energy density is given by
δ
U
0
=
𝝈
:
δ
ϵ
.
:
𝛿
subscript
𝑈
0
𝝈
𝛿
bold-italic-ϵ
{\displaystyle{\displaystyle\delta U_{0}={\boldsymbol{\sigma}}:\delta{%
\boldsymbol{\epsilon}}~{}.}}
An 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,
U
0
=
U
0
(
ϵ
)
subscript
𝑈
0
subscript
𝑈
0
bold-italic-ϵ
{\displaystyle{\displaystyle U_{0}=U_{0}({\boldsymbol{\epsilon}})}}
and the variation of the internal energy can be expressed as
δ
U
0
=
∂
U
0
∂
ϵ
:
δ
ϵ
.
:
𝛿
subscript
𝑈
0
continued-fraction
subscript
𝑈
0
bold-italic-ϵ
𝛿
bold-italic-ϵ
{\displaystyle{\displaystyle\delta U_{0}={\cfrac{\partial U_{0}}{\partial{%
\boldsymbol{\epsilon}}}}:\delta{\boldsymbol{\epsilon}}~{}.}}
Since the variation of strain is arbitrary, the stress–strain relation of an elastic material is given by
𝝈
=
∂
U
0
∂
ϵ
.
𝝈
continued-fraction
subscript
𝑈
0
bold-italic-ϵ
{\displaystyle{\displaystyle{\boldsymbol{\sigma}}={\cfrac{\partial U_{0}}{%
\partial{\boldsymbol{\epsilon}}}}~{}.}}
For a linear elastic material, the quantity
∂
U
0
/
∂
ϵ
subscript
𝑈
0
bold-italic-ϵ
{\displaystyle{\displaystyle\partial U_{0}/\partial{\boldsymbol{\epsilon}}}}
is a linear function of
ϵ
bold-italic-ϵ
{\displaystyle{\displaystyle{\boldsymbol{\epsilon}}}}
, and can therefore be expressed as
𝝈
=
𝖼
:
ϵ
:
𝝈
𝖼
bold-italic-ϵ
{\displaystyle{\displaystyle{\boldsymbol{\sigma}}={\mathsf{c}}:{\boldsymbol{%
\epsilon}}}}
where
𝖼
𝖼
{\displaystyle{\displaystyle{\mathsf{c}}}}
is a fourth-rank tensor of material constants, also called the stiffness tensor . We can see why
𝖼
𝖼
{\displaystyle{\displaystyle{\mathsf{c}}}}
must be a fourth-rank tensor by noting that, for a linear elastic material,
∂
∂
ϵ
[
𝝈
(
ϵ
)
]
=
constant
=
𝖼
.
continued-fraction
bold-italic-ϵ
delimited-[]
𝝈
bold-italic-ϵ
constant
𝖼
{\displaystyle{\displaystyle{\cfrac{\partial}{\partial{\boldsymbol{\epsilon}}}%
}[{\boldsymbol{\sigma}}({\boldsymbol{\epsilon}})]={\text{constant}}={\mathsf{c%
}}\,.}}
In index notation
∂
σ
i
j
∂
ϵ
k
ℓ
=
constant
=
c
i
j
k
ℓ
.
continued-fraction
subscript
𝜎
𝑖
𝑗
subscript
italic-ϵ
𝑘
ℓ
constant
subscript
𝑐
𝑖
𝑗
𝑘
ℓ
{\displaystyle{\displaystyle{\cfrac{\partial\sigma_{ij}}{\partial\epsilon_{k%
\ell}}}={\text{constant}}=c_{ijk\ell}\,.}}
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.
انظر أيضاً
الهوامش
^ قانون هوك، المرجع الالكتروني للمعلوماتية
^ Simo, J. C.; Hughes, T. J. R. (1998), Computational Inelasticity , Springer, ISBN 9780387975207
^ Milton, Graeme W. (2002), The Theory of Composites , Cambridge Monographs on Applied and Computational Mathematics, Cambridge University Press, ISBN 9780521781251
^ Slaughter, William S. (2001), The Linearized Theory of Elasticity , Birkhäuser, ISBN 978-0817641177
^ Boresi, A. P, Schmidt, R. J. and Sidebottom, O. M., 1993, Advanced Mechanics of Materials , Wiley.
^ Tan, S. C., 1994, Stress Concentrations in Laminated Composites , Technomic Publishing Company, Lancaster, PA.
المصادر
وصلات خارجية
أخف صيغ التحويل
Homogeneous isotropic linear elastic materials have their elastic properties uniquely determined by any two moduli among these; thus, given any two, any other of the elastic moduli can be calculated according to these formulas.
K
=
𝐾
absent
{\displaystyle{\displaystyle K=\,}}
E
=
𝐸
absent
{\displaystyle{\displaystyle E=\,}}
λ
=
𝜆
absent
{\displaystyle{\displaystyle\lambda=\,}}
G
=
𝐺
absent
{\displaystyle{\displaystyle G=\,}}
ν
=
𝜈
absent
{\displaystyle{\displaystyle\nu=\,}}
M
=
𝑀
absent
{\displaystyle{\displaystyle M=\,}}
Notes
(
K
,
E
)
𝐾
𝐸
{\displaystyle{\displaystyle(K,\,E)}}
K
𝐾
{\displaystyle{\displaystyle K}}
E
𝐸
{\displaystyle{\displaystyle E}}
3
K
(
3
K
-
E
)
9
K
-
E
3
𝐾
3
𝐾
𝐸
9
𝐾
𝐸
{\displaystyle{\displaystyle{\tfrac{3K(3K-E)}{9K-E}}}}
3
K
E
9
K
-
E
3
𝐾
𝐸
9
𝐾
𝐸
{\displaystyle{\displaystyle{\tfrac{3KE}{9K-E}}}}
3
K
-
E
6
K
3
𝐾
𝐸
6
𝐾
{\displaystyle{\displaystyle{\tfrac{3K-E}{6K}}}}
3
K
(
3
K
+
E
)
9
K
-
E
3
𝐾
3
𝐾
𝐸
9
𝐾
𝐸
{\displaystyle{\displaystyle{\tfrac{3K(3K+E)}{9K-E}}}}
(
K
,
λ
)
𝐾
𝜆
{\displaystyle{\displaystyle(K,\,\lambda)}}
K
𝐾
{\displaystyle{\displaystyle K}}
9
K
(
K
-
λ
)
3
K
-
λ
9
𝐾
𝐾
𝜆
3
𝐾
𝜆
{\displaystyle{\displaystyle{\tfrac{9K(K-\lambda)}{3K-\lambda}}}}
λ
𝜆
{\displaystyle{\displaystyle\lambda}}
3
(
K
-
λ
)
2
3
𝐾
𝜆
2
{\displaystyle{\displaystyle{\tfrac{3(K-\lambda)}{2}}}}
λ
3
K
-
λ
𝜆
3
𝐾
𝜆
{\displaystyle{\displaystyle{\tfrac{\lambda}{3K-\lambda}}}}
3
K
-
2
λ
3
𝐾
2
𝜆
{\displaystyle{\displaystyle 3K-2\lambda\,}}
(
K
,
G
)
𝐾
𝐺
{\displaystyle{\displaystyle(K,\,G)}}
K
𝐾
{\displaystyle{\displaystyle K}}
9
K
G
3
K
+
G
9
𝐾
𝐺
3
𝐾
𝐺
{\displaystyle{\displaystyle{\tfrac{9KG}{3K+G}}}}
K
-
2
G
3
𝐾
2
𝐺
3
{\displaystyle{\displaystyle K-{\tfrac{2G}{3}}}}
G
𝐺
{\displaystyle{\displaystyle G}}
3
K
-
2
G
2
(
3
K
+
G
)
3
𝐾
2
𝐺
2
3
𝐾
𝐺
{\displaystyle{\displaystyle{\tfrac{3K-2G}{2(3K+G)}}}}
K
+
4
G
3
𝐾
4
𝐺
3
{\displaystyle{\displaystyle K+{\tfrac{4G}{3}}}}
(
K
,
ν
)
𝐾
𝜈
{\displaystyle{\displaystyle(K,\,\nu)}}
K
𝐾
{\displaystyle{\displaystyle K}}
3
K
(
1
-
2
ν
)
3
𝐾
1
2
𝜈
{\displaystyle{\displaystyle 3K(1-2\nu)\,}}
3
K
ν
1
+
ν
3
𝐾
𝜈
1
𝜈
{\displaystyle{\displaystyle{\tfrac{3K\nu}{1+\nu}}}}
3
K
(
1
-
2
ν
)
2
(
1
+
ν
)
3
𝐾
1
2
𝜈
2
1
𝜈
{\displaystyle{\displaystyle{\tfrac{3K(1-2\nu)}{2(1+\nu)}}}}
ν
𝜈
{\displaystyle{\displaystyle\nu}}
3
K
(
1
-
ν
)
1
+
ν
3
𝐾
1
𝜈
1
𝜈
{\displaystyle{\displaystyle{\tfrac{3K(1-\nu)}{1+\nu}}}}
(
K
,
M
)
𝐾
𝑀
{\displaystyle{\displaystyle(K,\,M)}}
K
𝐾
{\displaystyle{\displaystyle K}}
9
K
(
M
-
K
)
3
K
+
M
9
𝐾
𝑀
𝐾
3
𝐾
𝑀
{\displaystyle{\displaystyle{\tfrac{9K(M-K)}{3K+M}}}}
3
K
-
M
2
3
𝐾
𝑀
2
{\displaystyle{\displaystyle{\tfrac{3K-M}{2}}}}
3
(
M
-
K
)
4
3
𝑀
𝐾
4
{\displaystyle{\displaystyle{\tfrac{3(M-K)}{4}}}}
3
K
-
M
3
K
+
M
3
𝐾
𝑀
3
𝐾
𝑀
{\displaystyle{\displaystyle{\tfrac{3K-M}{3K+M}}}}
M
𝑀
{\displaystyle{\displaystyle M}}
(
E
,
λ
)
𝐸
𝜆
{\displaystyle{\displaystyle(E,\,\lambda)}}
E
+
3
λ
+
R
6
𝐸
3
𝜆
𝑅
6
{\displaystyle{\displaystyle{\tfrac{E+3\lambda+R}{6}}}}
E
𝐸
{\displaystyle{\displaystyle E}}
λ
𝜆
{\displaystyle{\displaystyle\lambda}}
E
-
3
λ
+
R
4
𝐸
3
𝜆
𝑅
4
{\displaystyle{\displaystyle{\tfrac{E-3\lambda+R}{4}}}}
2
λ
E
+
λ
+
R
2
𝜆
𝐸
𝜆
𝑅
{\displaystyle{\displaystyle{\tfrac{2\lambda}{E+\lambda+R}}}}
E
-
λ
+
R
2
𝐸
𝜆
𝑅
2
{\displaystyle{\displaystyle{\tfrac{E-\lambda+R}{2}}}}
R
=
E
2
+
9
λ
2
+
2
E
λ
𝑅
superscript
𝐸
2
9
superscript
𝜆
2
2
𝐸
𝜆
{\displaystyle{\displaystyle R={\sqrt{E^{2}+9\lambda^{2}+2E\lambda}}}}
(
E
,
G
)
𝐸
𝐺
{\displaystyle{\displaystyle(E,\,G)}}
E
G
3
(
3
G
-
E
)
𝐸
𝐺
3
3
𝐺
𝐸
{\displaystyle{\displaystyle{\tfrac{EG}{3(3G-E)}}}}
E
𝐸
{\displaystyle{\displaystyle E}}
G
(
E
-
2
G
)
3
G
-
E
𝐺
𝐸
2
𝐺
3
𝐺
𝐸
{\displaystyle{\displaystyle{\tfrac{G(E-2G)}{3G-E}}}}
G
𝐺
{\displaystyle{\displaystyle G}}
E
2
G
-
1
𝐸
2
𝐺
1
{\displaystyle{\displaystyle{\tfrac{E}{2G}}-1}}
G
(
4
G
-
E
)
3
G
-
E
𝐺
4
𝐺
𝐸
3
𝐺
𝐸
{\displaystyle{\displaystyle{\tfrac{G(4G-E)}{3G-E}}}}
(
E
,
ν
)
𝐸
𝜈
{\displaystyle{\displaystyle(E,\,\nu)}}
E
3
(
1
-
2
ν
)
𝐸
3
1
2
𝜈
{\displaystyle{\displaystyle{\tfrac{E}{3(1-2\nu)}}}}
E
𝐸
{\displaystyle{\displaystyle E}}
E
ν
(
1
+
ν
)
(
1
-
2
ν
)
𝐸
𝜈
1
𝜈
1
2
𝜈
{\displaystyle{\displaystyle{\tfrac{E\nu}{(1+\nu)(1-2\nu)}}}}
E
2
(
1
+
ν
)
𝐸
2
1
𝜈
{\displaystyle{\displaystyle{\tfrac{E}{2(1+\nu)}}}}
ν
𝜈
{\displaystyle{\displaystyle\nu}}
E
(
1
-
ν
)
(
1
+
ν
)
(
1
-
2
ν
)
𝐸
1
𝜈
1
𝜈
1
2
𝜈
{\displaystyle{\displaystyle{\tfrac{E(1-\nu)}{(1+\nu)(1-2\nu)}}}}
(
E
,
M
)
𝐸
𝑀
{\displaystyle{\displaystyle(E,\,M)}}
3
M
-
E
+
S
6
3
𝑀
𝐸
𝑆
6
{\displaystyle{\displaystyle{\tfrac{3M-E+S}{6}}}}
E
𝐸
{\displaystyle{\displaystyle E}}
M
-
E
+
S
4
𝑀
𝐸
𝑆
4
{\displaystyle{\displaystyle{\tfrac{M-E+S}{4}}}}
3
M
+
E
-
S
8
3
𝑀
𝐸
𝑆
8
{\displaystyle{\displaystyle{\tfrac{3M+E-S}{8}}}}
E
-
M
+
S
4
M
𝐸
𝑀
𝑆
4
𝑀
{\displaystyle{\displaystyle{\tfrac{E-M+S}{4M}}}}
M
𝑀
{\displaystyle{\displaystyle M}}
S
=
±
E
2
+
9
M
2
-
10
E
M
𝑆
plus-or-minus
superscript
𝐸
2
9
superscript
𝑀
2
10
𝐸
𝑀
{\displaystyle{\displaystyle S=\pm{\sqrt{E^{2}+9M^{2}-10EM}}}}
There are two valid solutions.
The plus sign leads to
ν
≥
0
𝜈
0
{\displaystyle{\displaystyle\nu\geq 0}}
.
The minus sign leads to
ν
≤
0
𝜈
0
{\displaystyle{\displaystyle\nu\leq 0}}
.
(
λ
,
G
)
𝜆
𝐺
{\displaystyle{\displaystyle(\lambda,\,G)}}
λ
+
2
G
3
𝜆
2
𝐺
3
{\displaystyle{\displaystyle\lambda+{\tfrac{2G}{3}}}}
G
(
3
λ
+
2
G
)
λ
+
G
𝐺
3
𝜆
2
𝐺
𝜆
𝐺
{\displaystyle{\displaystyle{\tfrac{G(3\lambda+2G)}{\lambda+G}}}}
λ
𝜆
{\displaystyle{\displaystyle\lambda}}
G
𝐺
{\displaystyle{\displaystyle G}}
λ
2
(
λ
+
G
)
𝜆
2
𝜆
𝐺
{\displaystyle{\displaystyle{\tfrac{\lambda}{2(\lambda+G)}}}}
λ
+
2
G
𝜆
2
𝐺
{\displaystyle{\displaystyle\lambda+2G\,}}
(
λ
,
ν
)
𝜆
𝜈
{\displaystyle{\displaystyle(\lambda,\,\nu)}}
λ
(
1
+
ν
)
3
ν
𝜆
1
𝜈
3
𝜈
{\displaystyle{\displaystyle{\tfrac{\lambda(1+\nu)}{3\nu}}}}
λ
(
1
+
ν
)
(
1
-
2
ν
)
ν
𝜆
1
𝜈
1
2
𝜈
𝜈
{\displaystyle{\displaystyle{\tfrac{\lambda(1+\nu)(1-2\nu)}{\nu}}}}
λ
𝜆
{\displaystyle{\displaystyle\lambda}}
λ
(
1
-
2
ν
)
2
ν
𝜆
1
2
𝜈
2
𝜈
{\displaystyle{\displaystyle{\tfrac{\lambda(1-2\nu)}{2\nu}}}}
ν
𝜈
{\displaystyle{\displaystyle\nu}}
λ
(
1
-
ν
)
ν
𝜆
1
𝜈
𝜈
{\displaystyle{\displaystyle{\tfrac{\lambda(1-\nu)}{\nu}}}}
Cannot be used when
ν
=
0
⇔
λ
=
0
⇔
𝜈
0
𝜆
0
{\displaystyle{\displaystyle\nu=0\Leftrightarrow\lambda=0}}
(
λ
,
M
)
𝜆
𝑀
{\displaystyle{\displaystyle(\lambda,\,M)}}
M
+
2
λ
3
𝑀
2
𝜆
3
{\displaystyle{\displaystyle{\tfrac{M+2\lambda}{3}}}}
(
M
-
λ
)
(
M
+
2
λ
)
M
+
λ
𝑀
𝜆
𝑀
2
𝜆
𝑀
𝜆
{\displaystyle{\displaystyle{\tfrac{(M-\lambda)(M+2\lambda)}{M+\lambda}}}}
λ
𝜆
{\displaystyle{\displaystyle\lambda}}
M
-
λ
2
𝑀
𝜆
2
{\displaystyle{\displaystyle{\tfrac{M-\lambda}{2}}}}
λ
M
+
λ
𝜆
𝑀
𝜆
{\displaystyle{\displaystyle{\tfrac{\lambda}{M+\lambda}}}}
M
𝑀
{\displaystyle{\displaystyle M}}
(
G
,
ν
)
𝐺
𝜈
{\displaystyle{\displaystyle(G,\,\nu)}}
2
G
(
1
+
ν
)
3
(
1
-
2
ν
)
2
𝐺
1
𝜈
3
1
2
𝜈
{\displaystyle{\displaystyle{\tfrac{2G(1+\nu)}{3(1-2\nu)}}}}
2
G
(
1
+
ν
)
2
𝐺
1
𝜈
{\displaystyle{\displaystyle 2G(1+\nu)\,}}
2
G
ν
1
-
2
ν
2
𝐺
𝜈
1
2
𝜈
{\displaystyle{\displaystyle{\tfrac{2G\nu}{1-2\nu}}}}
G
𝐺
{\displaystyle{\displaystyle G}}
ν
𝜈
{\displaystyle{\displaystyle\nu}}
2
G
(
1
-
ν
)
1
-
2
ν
2
𝐺
1
𝜈
1
2
𝜈
{\displaystyle{\displaystyle{\tfrac{2G(1-\nu)}{1-2\nu}}}}
(
G
,
M
)
𝐺
𝑀
{\displaystyle{\displaystyle(G,\,M)}}
M
-
4
G
3
𝑀
4
𝐺
3
{\displaystyle{\displaystyle M-{\tfrac{4G}{3}}}}
G
(
3
M
-
4
G
)
M
-
G
𝐺
3
𝑀
4
𝐺
𝑀
𝐺
{\displaystyle{\displaystyle{\tfrac{G(3M-4G)}{M-G}}}}
M
-
2
G
𝑀
2
𝐺
{\displaystyle{\displaystyle M-2G\,}}
G
𝐺
{\displaystyle{\displaystyle G}}
M
-
2
G
2
M
-
2
G
𝑀
2
𝐺
2
𝑀
2
𝐺
{\displaystyle{\displaystyle{\tfrac{M-2G}{2M-2G}}}}
M
𝑀
{\displaystyle{\displaystyle M}}
(
ν
,
M
)
𝜈
𝑀
{\displaystyle{\displaystyle(\nu,\,M)}}
M
(
1
+
ν
)
3
(
1
-
ν
)
𝑀
1
𝜈
3
1
𝜈
{\displaystyle{\displaystyle{\tfrac{M(1+\nu)}{3(1-\nu)}}}}
M
(
1
+
ν
)
(
1
-
2
ν
)
1
-
ν
𝑀
1
𝜈
1
2
𝜈
1
𝜈
{\displaystyle{\displaystyle{\tfrac{M(1+\nu)(1-2\nu)}{1-\nu}}}}
M
ν
1
-
ν
𝑀
𝜈
1
𝜈
{\displaystyle{\displaystyle{\tfrac{M\nu}{1-\nu}}}}
M
(
1
-
2
ν
)
2
(
1
-
ν
)
𝑀
1
2
𝜈
2
1
𝜈
{\displaystyle{\displaystyle{\tfrac{M(1-2\nu)}{2(1-\nu)}}}}
ν
𝜈
{\displaystyle{\displaystyle\nu}}
M
𝑀
{\displaystyle{\displaystyle M}}