令
Ω
0
{\displaystyle \Omega _{0}}
為區域
Ω
(
t
)
{\displaystyle \Omega (t)}
的參考組態,令其運動及形變梯度為
x
=
φ
(
X
,
t
)
;
⟹
F
(
X
,
t
)
=
∇
∘
φ
.
{\displaystyle \mathbf {x} ={\boldsymbol {\varphi }}(\mathbf {X} ,t)~;\qquad \implies \qquad {\boldsymbol {F}}(\mathbf {X} ,t)={\boldsymbol {\nabla }}_{\circ }{\boldsymbol {\varphi }}~.}
令
J
(
X
,
t
)
=
det
[
F
(
X
,
t
)
]
{\displaystyle J(\mathbf {X} ,t)=\det[{\boldsymbol {F}}(\mathbf {X} ,t)]}
.
則目前組態及參考組態的積分有以下的關係
∫
Ω
(
t
)
f
(
x
,
t
)
dV
=
∫
Ω
0
f
[
φ
(
X
,
t
)
,
t
]
J
(
X
,
t
)
dV
0
=
∫
Ω
0
f
^
(
X
,
t
)
J
(
X
,
t
)
dV
0
.
{\displaystyle \int _{\Omega (t)}\mathbf {f} (\mathbf {x} ,t)~{\text{dV}}=\int _{\Omega _{0}}\mathbf {f} [{\boldsymbol {\varphi }}(\mathbf {X} ,t),t]~J(\mathbf {X} ,t)~{\text{dV}}_{0}=\int _{\Omega _{0}}{\hat {\mathbf {f} }}(\mathbf {X} ,t)~J(\mathbf {X} ,t)~{\text{dV}}_{0}~.}
That this derivation is for a material element is implicit in the time constancy of the reference configuration: it is constant in material coordinates. 針對體積積分的微分定義為
d
d
t
(
∫
Ω
(
t
)
f
(
x
,
t
)
dV
)
=
lim
Δ
t
→
0
1
Δ
t
(
∫
Ω
(
t
+
Δ
t
)
f
(
x
,
t
+
Δ
t
)
dV
−
∫
Ω
(
t
)
f
(
x
,
t
)
dV
)
.
{\displaystyle {\cfrac {\mathrm {d} }{\mathrm {d} t}}\left(\int _{\Omega (t)}\mathbf {f} (\mathbf {x} ,t)~{\text{dV}}\right)=\lim _{\Delta t\rightarrow 0}{\cfrac {1}{\Delta t}}\left(\int _{\Omega (t+\Delta t)}\mathbf {f} (\mathbf {x} ,t+\Delta t)~{\text{dV}}-\int _{\Omega (t)}\mathbf {f} (\mathbf {x} ,t)~{\text{dV}}\right)~.}
將上式轉換為對參考組態的積分,可得
d
d
t
(
∫
Ω
(
t
)
f
(
x
,
t
)
dV
)
=
lim
Δ
t
→
0
1
Δ
t
(
∫
Ω
0
f
^
(
X
,
t
+
Δ
t
)
J
(
X
,
t
+
Δ
t
)
dV
0
−
∫
Ω
0
f
^
(
X
,
t
)
J
(
X
,
t
)
dV
0
)
.
{\displaystyle {\cfrac {\mathrm {d} }{\mathrm {d} t}}\left(\int _{\Omega (t)}\mathbf {f} (\mathbf {x} ,t)~{\text{dV}}\right)=\lim _{\Delta t\rightarrow 0}{\cfrac {1}{\Delta t}}\left(\int _{\Omega _{0}}{\hat {\mathbf {f} }}(\mathbf {X} ,t+\Delta t)~J(\mathbf {X} ,t+\Delta t)~{\text{dV}}_{0}-\int _{\Omega _{0}}{\hat {\mathbf {f} }}(\mathbf {X} ,t)~J(\mathbf {X} ,t)~{\text{dV}}_{0}\right)~.}
因為
Ω
0
{\displaystyle \Omega _{0}}
和時間無關,可得
d
d
t
(
∫
Ω
(
t
)
f
(
x
,
t
)
dV
)
=
∫
Ω
0
[
lim
Δ
t
→
0
f
^
(
X
,
t
+
Δ
t
)
J
(
X
,
t
+
Δ
t
)
−
f
^
(
X
,
t
)
J
(
X
,
t
)
Δ
t
]
dV
0
=
∫
Ω
0
∂
∂
t
[
f
^
(
X
,
t
)
J
(
X
,
t
)
]
dV
0
=
∫
Ω
0
(
∂
∂
t
[
f
^
(
X
,
t
)
]
J
(
X
,
t
)
+
f
^
(
X
,
t
)
∂
∂
t
[
J
(
X
,
t
)
]
)
dV
0
{\displaystyle {\begin{aligned}{\cfrac {\mathrm {d} }{\mathrm {d} t}}\left(\int _{\Omega (t)}\mathbf {f} (\mathbf {x} ,t)~{\text{dV}}\right)&=\int _{\Omega _{0}}\left[\lim _{\Delta t\rightarrow 0}{\cfrac {{\hat {\mathbf {f} }}(\mathbf {X} ,t+\Delta t)~J(\mathbf {X} ,t+\Delta t)-{\hat {\mathbf {f} }}(\mathbf {X} ,t)~J(\mathbf {X} ,t)}{\Delta t}}\right]~{\text{dV}}_{0}\\&=\int _{\Omega _{0}}{\frac {\partial }{\partial t}}[{\hat {\mathbf {f} }}(\mathbf {X} ,t)~J(\mathbf {X} ,t)]~{\text{dV}}_{0}\\&=\int _{\Omega _{0}}\left({\frac {\partial }{\partial t}}[{\hat {\mathbf {f} }}(\mathbf {X} ,t)]~J(\mathbf {X} ,t)+{\hat {\mathbf {f} }}(\mathbf {X} ,t)~{\frac {\partial }{\partial t}}[J(\mathbf {X} ,t)]\right)~{\text{dV}}_{0}\end{aligned}}}
現在,
det
F
{\displaystyle \det {\boldsymbol {F}}}
的時間導數為
[ 6]
∂
J
(
X
,
t
)
∂
t
=
∂
∂
t
(
det
F
)
=
(
det
F
)
(
∇
⋅
v
)
=
J
(
X
,
t
)
∇
⋅
v
(
φ
(
X
,
t
)
,
t
)
=
J
(
X
,
t
)
∇
⋅
v
(
x
,
t
)
.
{\displaystyle {\frac {\partial J(\mathbf {X} ,t)}{\partial t}}={\frac {\partial }{\partial t}}(\det {\boldsymbol {F}})=(\det {\boldsymbol {F}})({\boldsymbol {\nabla }}\cdot \mathbf {v} )=J(\mathbf {X} ,t)~{\boldsymbol {\nabla }}\cdot \mathbf {v} ({\boldsymbol {\varphi }}(\mathbf {X} ,t),t)=J(\mathbf {X} ,t)~{\boldsymbol {\nabla }}\cdot \mathbf {v} (\mathbf {x} ,t)~.}
因此
d
d
t
(
∫
Ω
(
t
)
f
(
x
,
t
)
dV
)
=
∫
Ω
0
(
∂
∂
t
[
f
^
(
X
,
t
)
]
J
(
X
,
t
)
+
f
^
(
X
,
t
)
J
(
X
,
t
)
∇
⋅
v
(
x
,
t
)
)
dV
0
=
∫
Ω
0
(
∂
∂
t
[
f
^
(
X
,
t
)
]
+
f
^
(
X
,
t
)
∇
⋅
v
(
x
,
t
)
)
J
(
X
,
t
)
dV
0
=
∫
Ω
(
t
)
(
f
˙
(
x
,
t
)
+
f
(
x
,
t
)
∇
⋅
v
(
x
,
t
)
)
dV
{\displaystyle {\begin{aligned}{\cfrac {\mathrm {d} }{\mathrm {d} t}}\left(\int _{\Omega (t)}\mathbf {f} (\mathbf {x} ,t)~{\text{dV}}\right)&=\int _{\Omega _{0}}\left({\frac {\partial }{\partial t}}[{\hat {\mathbf {f} }}(\mathbf {X} ,t)]~J(\mathbf {X} ,t)+{\hat {\mathbf {f} }}(\mathbf {X} ,t)~J(\mathbf {X} ,t)~{\boldsymbol {\nabla }}\cdot \mathbf {v} (\mathbf {x} ,t)\right)~{\text{dV}}_{0}\\&=\int _{\Omega _{0}}\left({\frac {\partial }{\partial t}}[{\hat {\mathbf {f} }}(\mathbf {X} ,t)]+{\hat {\mathbf {f} }}(\mathbf {X} ,t)~{\boldsymbol {\nabla }}\cdot \mathbf {v} (\mathbf {x} ,t)\right)~J(\mathbf {X} ,t)~{\text{dV}}_{0}\\&=\int _{\Omega (t)}\left({\dot {\mathbf {f} }}(\mathbf {x} ,t)+\mathbf {f} (\mathbf {x} ,t)~{\boldsymbol {\nabla }}\cdot \mathbf {v} (\mathbf {x} ,t)\right)~{\text{dV}}\end{aligned}}}
其中
f
˙
{\displaystyle {\dot {\mathbf {f} }}}
為
f
{\displaystyle \mathbf {f} }
的材料導數 ,現在材料導數為
f
˙
(
x
,
t
)
=
∂
f
(
x
,
t
)
∂
t
+
[
∇
f
(
x
,
t
)
]
⋅
v
(
x
,
t
)
.
{\displaystyle {\dot {\mathbf {f} }}(\mathbf {x} ,t)={\frac {\partial \mathbf {f} (\mathbf {x} ,t)}{\partial t}}+[{\boldsymbol {\nabla }}\mathbf {f} (\mathbf {x} ,t)]\cdot \mathbf {v} (\mathbf {x} ,t)~.}
因此
d
d
t
(
∫
Ω
(
t
)
f
(
x
,
t
)
dV
)
=
∫
Ω
(
t
)
(
∂
f
(
x
,
t
)
∂
t
+
[
∇
f
(
x
,
t
)
]
⋅
v
(
x
,
t
)
+
f
(
x
,
t
)
∇
⋅
v
(
x
,
t
)
)
dV
{\displaystyle {\cfrac {\mathrm {d} }{\mathrm {d} t}}\left(\int _{\Omega (t)}\mathbf {f} (\mathbf {x} ,t)~{\text{dV}}\right)=\int _{\Omega (t)}\left({\frac {\partial \mathbf {f} (\mathbf {x} ,t)}{\partial t}}+[{\boldsymbol {\nabla }}\mathbf {f} (\mathbf {x} ,t)]\cdot \mathbf {v} (\mathbf {x} ,t)+\mathbf {f} (\mathbf {x} ,t)~{\boldsymbol {\nabla }}\cdot \mathbf {v} (\mathbf {x} ,t)\right)~{\text{dV}}}
或者
d
d
t
(
∫
Ω
(
t
)
f
dV
)
=
∫
Ω
(
t
)
(
∂
f
∂
t
+
∇
f
⋅
v
+
f
∇
⋅
v
)
dV
.
{\displaystyle {\cfrac {\mathrm {d} }{\mathrm {d} t}}\left(\int _{\Omega (t)}\mathbf {f} ~{\text{dV}}\right)=\int _{\Omega (t)}\left({\frac {\partial \mathbf {f} }{\partial t}}+{\boldsymbol {\nabla }}\mathbf {f} \cdot \mathbf {v} +\mathbf {f} ~{\boldsymbol {\nabla }}\cdot \mathbf {v} \right)~{\text{dV}}~.}
利用以下的恆等式
∇
⋅
(
v
⊗
w
)
=
v
(
∇
⋅
w
)
+
∇
v
⋅
w
{\displaystyle {\boldsymbol {\nabla }}\cdot (\mathbf {v} \otimes \mathbf {w} )=\mathbf {v} ({\boldsymbol {\nabla }}\cdot \mathbf {w} )+{\boldsymbol {\nabla }}\mathbf {v} \cdot \mathbf {w} }
可得
d
d
t
(
∫
Ω
(
t
)
f
dV
)
=
∫
Ω
(
t
)
(
∂
f
∂
t
+
∇
⋅
(
f
⊗
v
)
)
dV
.
{\displaystyle {\cfrac {\mathrm {d} }{\mathrm {d} t}}\left(\int _{\Omega (t)}\mathbf {f} ~{\text{dV}}\right)=\int _{\Omega (t)}\left({\frac {\partial \mathbf {f} }{\partial t}}+{\boldsymbol {\nabla }}\cdot (\mathbf {f} \otimes \mathbf {v} )\right)~{\text{dV}}~.}
利用高斯散度定理 及恆等式
(
a
⊗
b
)
⋅
n
=
(
b
⋅
n
)
a
{\displaystyle (\mathbf {a} \otimes \mathbf {b} )\cdot \mathbf {n} =(\mathbf {b} \cdot \mathbf {n} )\mathbf {a} }
,可得
d
d
t
(
∫
Ω
(
t
)
f
dV
)
=
∫
Ω
(
t
)
∂
f
∂
t
dV
+
∫
∂
Ω
(
t
)
(
f
⊗
v
)
⋅
n
dA
=
∫
Ω
(
t
)
∂
f
∂
t
dV
+
∫
∂
Ω
(
t
)
(
v
⋅
n
)
f
dA
◻
{\displaystyle {{\cfrac {\mathrm {d} }{\mathrm {d} t}}\left(\int _{\Omega (t)}\mathbf {f} ~{\text{dV}}\right)=\int _{\Omega (t)}{\frac {\partial \mathbf {f} }{\partial t}}~{\text{dV}}+\int _{\partial \Omega (t)}(\mathbf {f} \otimes \mathbf {v} )\cdot \mathbf {n} ~{\text{dA}}=\int _{\Omega (t)}{\frac {\partial \mathbf {f} }{\partial t}}~{\text{dV}}+\int _{\partial \Omega (t)}(\mathbf {v} \cdot \mathbf {n} )\mathbf {f} ~{\text{dA}}\qquad \square }}