无穷小应变理论 (infinitesimal strain theory)也称为无限小应变理论 ,是连续介质力学 中描述固体形变 的数学分析法,适用在其形变量远小于物体尺寸(无穷小量 )的情形,因此若是均质材料,可以假设材料每一点的结构性质(密度 及刚度 )都相等,不会随变形而不同。
在此假设下,连续介质力学的方程可以简化。此作法也称为是小形变理论 、小位移理论 或小位移梯度理论 。无穷小应变理论和有限应变理论 的假设恰好相反,后者假设形变量没有远小于物体尺寸。
无穷小应变理论常用在土木工程 及机械工程 中,其中会进行结构的应力分析 ,而材料是用强度较高的混凝土 及钢 制成,而结构设计的目标也是在一般结构荷重 下,希望其形变量可以降到最小。不过若分析的结构物是较细较薄,较容易变形的元件(例如杆、平板及薄壳),用无限小应变理论来分析就不可靠了[ 1] 。
在连续体 的无限小变形中(位移梯度张量 远小于1,也就是
‖
∇
u
‖
≪
1
{\displaystyle \|\nabla \mathbf {u} \|\ll 1}
),可以用有限应变理论中的任何一个有限应变张量(例如拉格朗日有限应变张量
E
{\displaystyle \mathbf {E} }
,或是尤拉有限应变张量
e
{\displaystyle \mathbf {e} }
)进行线性化。在线性化中,可以省略有限应变张量中的二次项或是非线性项,因此可得
E
=
1
2
(
∇
X
u
+
(
∇
X
u
)
T
+
(
∇
X
u
)
T
∇
X
u
)
≈
1
2
(
∇
X
u
+
(
∇
X
u
)
T
)
{\displaystyle \mathbf {E} ={\frac {1}{2}}\left(\nabla _{\mathbf {X} }\mathbf {u} +(\nabla _{\mathbf {X} }\mathbf {u} )^{T}+(\nabla _{\mathbf {X} }\mathbf {u} )^{T}\nabla _{\mathbf {X} }\mathbf {u} \right)\approx {\frac {1}{2}}\left(\nabla _{\mathbf {X} }\mathbf {u} +(\nabla _{\mathbf {X} }\mathbf {u} )^{T}\right)}
或
E
K
L
=
1
2
(
∂
U
K
∂
X
L
+
∂
U
L
∂
X
K
+
∂
U
M
∂
X
K
∂
U
M
∂
X
L
)
≈
1
2
(
∂
U
K
∂
X
L
+
∂
U
L
∂
X
K
)
{\displaystyle E_{KL}={\frac {1}{2}}\left({\frac {\partial U_{K}}{\partial X_{L}}}+{\frac {\partial U_{L}}{\partial X_{K}}}+{\frac {\partial U_{M}}{\partial X_{K}}}{\frac {\partial U_{M}}{\partial X_{L}}}\right)\approx {\frac {1}{2}}\left({\frac {\partial U_{K}}{\partial X_{L}}}+{\frac {\partial U_{L}}{\partial X_{K}}}\right)}
以及
e
=
1
2
(
∇
x
u
+
(
∇
x
u
)
T
−
∇
x
u
(
∇
x
u
)
T
)
≈
1
2
(
∇
x
u
+
(
∇
x
u
)
T
)
{\displaystyle \mathbf {e} ={\frac {1}{2}}\left(\nabla _{\mathbf {x} }\mathbf {u} +(\nabla _{\mathbf {x} }\mathbf {u} )^{T}-\nabla _{\mathbf {x} }\mathbf {u} (\nabla _{\mathbf {x} }\mathbf {u} )^{T}\right)\approx {\frac {1}{2}}\left(\nabla _{\mathbf {x} }\mathbf {u} +(\nabla _{\mathbf {x} }\mathbf {u} )^{T}\right)}
或
e
r
s
=
1
2
(
∂
u
r
∂
x
s
+
∂
u
s
∂
x
r
−
∂
u
k
∂
x
r
∂
u
k
∂
x
s
)
≈
1
2
(
∂
u
r
∂
x
s
+
∂
u
s
∂
x
r
)
{\displaystyle e_{rs}={\frac {1}{2}}\left({\frac {\partial u_{r}}{\partial x_{s}}}+{\frac {\partial u_{s}}{\partial x_{r}}}-{\frac {\partial u_{k}}{\partial x_{r}}}{\frac {\partial u_{k}}{\partial x_{s}}}\right)\approx {\frac {1}{2}}\left({\frac {\partial u_{r}}{\partial x_{s}}}+{\frac {\partial u_{s}}{\partial x_{r}}}\right)}
线性化意味着连续体中特定点的物质坐标(material coordinate)和空间坐标(spatial coordinate)差异很小,拉格朗日描述和尤拉描述近似相等。因此,物质位移梯度张量和空间位移梯度张量的分量也相近相等。可得
E
≈
e
≈
ε
=
1
2
(
(
∇
u
)
T
+
∇
u
)
{\displaystyle \mathbf {E} \approx \mathbf {e} \approx {\boldsymbol {\varepsilon }}={\frac {1}{2}}\left((\nabla \mathbf {u} )^{T}+\nabla \mathbf {u} \right)}
或
E
K
L
≈
e
r
s
≈
ε
i
j
=
1
2
(
u
i
,
j
+
u
j
,
i
)
{\displaystyle E_{KL}\approx e_{rs}\approx \varepsilon _{ij}={\frac {1}{2}}\left(u_{i,j}+u_{j,i}\right)}
其中
ε
i
j
{\displaystyle \varepsilon _{ij}}
是无穷小应变张量(也称为柯西应变张量、线性应变张量、小应变张量)的分量。
ε
i
j
=
1
2
(
u
i
,
j
+
u
j
,
i
)
=
[
ε
11
ε
12
ε
13
ε
21
ε
22
ε
23
ε
31
ε
32
ε
33
]
=
[
∂
u
1
∂
x
1
1
2
(
∂
u
1
∂
x
2
+
∂
u
2
∂
x
1
)
1
2
(
∂
u
1
∂
x
3
+
∂
u
3
∂
x
1
)
1
2
(
∂
u
2
∂
x
1
+
∂
u
1
∂
x
2
)
∂
u
2
∂
x
2
1
2
(
∂
u
2
∂
x
3
+
∂
u
3
∂
x
2
)
1
2
(
∂
u
3
∂
x
1
+
∂
u
1
∂
x
3
)
1
2
(
∂
u
3
∂
x
2
+
∂
u
2
∂
x
3
)
∂
u
3
∂
x
3
]
{\displaystyle {\begin{aligned}\varepsilon _{ij}&={\frac {1}{2}}\left(u_{i,j}+u_{j,i}\right)\\&={\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{21}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon _{31}&\varepsilon _{32}&\varepsilon _{33}\\\end{bmatrix}}\\&={\begin{bmatrix}{\frac {\partial u_{1}}{\partial x_{1}}}&{\frac {1}{2}}\left({\frac {\partial u_{1}}{\partial x_{2}}}+{\frac {\partial u_{2}}{\partial x_{1}}}\right)&{\frac {1}{2}}\left({\frac {\partial u_{1}}{\partial x_{3}}}+{\frac {\partial u_{3}}{\partial x_{1}}}\right)\\{\frac {1}{2}}\left({\frac {\partial u_{2}}{\partial x_{1}}}+{\frac {\partial u_{1}}{\partial x_{2}}}\right)&{\frac {\partial u_{2}}{\partial x_{2}}}&{\frac {1}{2}}\left({\frac {\partial u_{2}}{\partial x_{3}}}+{\frac {\partial u_{3}}{\partial x_{2}}}\right)\\{\frac {1}{2}}\left({\frac {\partial u_{3}}{\partial x_{1}}}+{\frac {\partial u_{1}}{\partial x_{3}}}\right)&{\frac {1}{2}}\left({\frac {\partial u_{3}}{\partial x_{2}}}+{\frac {\partial u_{2}}{\partial x_{3}}}\right)&{\frac {\partial u_{3}}{\partial x_{3}}}\\\end{bmatrix}}\end{aligned}}}
或者使用不同的表示方式:
[
ε
x
x
ε
x
y
ε
x
z
ε
y
x
ε
y
y
ε
y
z
ε
z
x
ε
z
y
ε
z
z
]
=
[
∂
u
x
∂
x
1
2
(
∂
u
x
∂
y
+
∂
u
y
∂
x
)
1
2
(
∂
u
x
∂
z
+
∂
u
z
∂
x
)
1
2
(
∂
u
y
∂
x
+
∂
u
x
∂
y
)
∂
u
y
∂
y
1
2
(
∂
u
y
∂
z
+
∂
u
z
∂
y
)
1
2
(
∂
u
z
∂
x
+
∂
u
x
∂
z
)
1
2
(
∂
u
z
∂
y
+
∂
u
y
∂
z
)
∂
u
z
∂
z
]
{\displaystyle {\begin{bmatrix}\varepsilon _{xx}&\varepsilon _{xy}&\varepsilon _{xz}\\\varepsilon _{yx}&\varepsilon _{yy}&\varepsilon _{yz}\\\varepsilon _{zx}&\varepsilon _{zy}&\varepsilon _{zz}\\\end{bmatrix}}={\begin{bmatrix}{\frac {\partial u_{x}}{\partial x}}&{\frac {1}{2}}\left({\frac {\partial u_{x}}{\partial y}}+{\frac {\partial u_{y}}{\partial x}}\right)&{\frac {1}{2}}\left({\frac {\partial u_{x}}{\partial z}}+{\frac {\partial u_{z}}{\partial x}}\right)\\{\frac {1}{2}}\left({\frac {\partial u_{y}}{\partial x}}+{\frac {\partial u_{x}}{\partial y}}\right)&{\frac {\partial u_{y}}{\partial y}}&{\frac {1}{2}}\left({\frac {\partial u_{y}}{\partial z}}+{\frac {\partial u_{z}}{\partial y}}\right)\\{\frac {1}{2}}\left({\frac {\partial u_{z}}{\partial x}}+{\frac {\partial u_{x}}{\partial z}}\right)&{\frac {1}{2}}\left({\frac {\partial u_{z}}{\partial y}}+{\frac {\partial u_{y}}{\partial z}}\right)&{\frac {\partial u_{z}}{\partial z}}\\\end{bmatrix}}}
进一步来说,因为形变梯度可以表示成
F
=
∇
u
+
I
{\displaystyle {\boldsymbol {F}}={\boldsymbol {\nabla }}\mathbf {u} +{\boldsymbol {I}}}
其中
I
{\displaystyle {\boldsymbol {I}}}
是二阶单位张量,可得
ε
=
1
2
(
F
T
+
F
)
−
I
{\displaystyle {\boldsymbol {\varepsilon }}={\frac {1}{2}}\left({\boldsymbol {F}}^{T}+{\boldsymbol {F}}\right)-{\boldsymbol {I}}}
另外,根据拉格朗日有限应变张量及尤拉有限应变张量的通用表示法,可得
E
(
m
)
=
1
2
m
(
U
2
m
−
I
)
=
1
2
m
[
(
F
T
F
)
m
−
I
]
≈
1
2
m
[
{
∇
u
+
(
∇
u
)
T
+
I
}
m
−
I
]
≈
ε
e
(
m
)
=
1
2
m
(
V
2
m
−
I
)
=
1
2
m
[
(
F
F
T
)
m
−
I
]
≈
ε
{\displaystyle {\begin{aligned}\mathbf {E} _{(m)}&={\frac {1}{2m}}(\mathbf {U} ^{2m}-{\boldsymbol {I}})={\frac {1}{2m}}[({\boldsymbol {F}}^{T}{\boldsymbol {F}})^{m}-{\boldsymbol {I}}]\approx {\frac {1}{2m}}[\{{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{T}+{\boldsymbol {I}}\}^{m}-{\boldsymbol {I}}]\approx {\boldsymbol {\varepsilon }}\\\mathbf {e} _{(m)}&={\frac {1}{2m}}(\mathbf {V} ^{2m}-{\boldsymbol {I}})={\frac {1}{2m}}[({\boldsymbol {F}}{\boldsymbol {F}}^{T})^{m}-{\boldsymbol {I}}]\approx {\boldsymbol {\varepsilon }}\end{aligned}}}
^ Boresi, Arthur P. (Arthur Peter), 1924-. Advanced mechanics of materials. Schmidt, Richard J. (Richard Joseph), 1954- 6th. New York: John Wiley & Sons. 2003: 62. ISBN 1601199228 . OCLC 430194205 .