（六）速度梯度与加速度梯度

1. 速度梯度、变形率张量与物质旋律张量

对于速度场 $v=v(x,t)=v^i{g}_i$，其右梯度 (称作速度梯度) 可写作
$$\begin{array}{rl}& v▽=\frac{\partial v}{\partial x^j}\otimes g{}^j=v^i{|}_j{g}_i\otimes g^j\triangleq \mathbf{L}\\ & \\ & =\frac{\partial v}{\partial X^A}\frac{\partial X^A}{\partial x^j}\otimes g{}^j=\frac{\partial v}{\partial X^A}\otimes C^A=v^B|{|}_A{C}_B\otimes C^A={\stackrel{\bullet }{C}}_A\otimes C^A\\ & \\ & =({\stackrel{\bullet }{C}}_B\otimes G^B)\cdot (G_A\otimes C^A)=\stackrel{\bullet }{\mathbf{F}}\cdot \stackrel{-1}{\mathbf{F}}\end{array}$$
将速度梯度的对称部分称作变形率张量：
$$\mathbf{D}=\frac{1}{2}(\mathbf{L}+{\mathbf{L}}^T)=\frac{1}{2}(v▽+▽v)$$
将速度梯度的反对称部分称作物质旋率张量：
$$\mathbf{W}=\frac{1}{2}(\mathbf{L}-{\mathbf{L}}^T)=\frac{1}{2}(v▽-▽v)$$
其轴向量为：
$$\omega =\frac{1}{2}(▽\times v)$$
此外还可以表示为：
$$\omega =-\frac{1}{4}ϵ:({\stackrel{\bullet }{C}}_A\otimes C^A-C^A\otimes {\stackrel{\bullet }{C}}_A)=\frac{1}{2}{C}^A\times {\stackrel{\bullet }{C}}_A$$
定义 处处满足 $w=0$ 或 $\mathbf{W}=0$ 的运动称为无旋运动。

2. 加速度梯度

加速度为：
$$\begin{array}{rl}& a=\stackrel{\bullet }{v}=v^{\prime }+(v▽)\cdot v\\ & \\ & =v^{\prime }+\mathbf{L}\cdot v=v^{\prime }+(\mathbf{D}+\mathbf{W})\cdot v\\ & \\ & =v^{\prime }+\mathbf{D}\cdot v+\omega \times v\\ & \\ & =v^{\prime }+\frac{1}{2}▽(v\cdot v)+2\mathbf{W}\cdot v\end{array}$$
其中，最后一式是由于
$$\mathbf{W}\cdot v=w\times v=\frac{1}{2}(▽\times v)\times v=\frac{1}{2}(v▽-▽v)\cdot v=\frac{1}{2}(\mathbf{L}\cdot v-▽v\cdot v)=\frac{1}{2}[\mathbf{L}\cdot v-\frac{1}{2}▽(v\cdot v)]$$

加速度的 (右) 梯度 定义为：
$$\begin{array}{rl}& a▽=\frac{\partial a}{\partial x^i}\otimes g^i=\frac{\partial a}{\partial X^A}\otimes C^A=a^B|{|}_A{C}_B\otimes C^A\\ & \\ & ={\stackrel{\bullet \bullet }{C}}_A\otimes C^A=({\stackrel{\bullet \bullet }{C}}_A\otimes G^A)\cdot (G_B\otimes C^B)=\stackrel{\bullet \bullet }{\mathbf{F}}\cdot \stackrel{-1}{\mathbf{F}}\end{array}$$

继续讨论前，给出如下命题：

对于任意仿射量 $\Psi$，有：
$${\left(\stackrel{\bullet }{\Psi }\right)}^T=\stackrel{\bullet }{\left({\Psi }^T\right)}$$
证明：对任意向量 $a，b$ 满足：
$$a\cdot \mathbf{\Psi }\cdot b=b\cdot {\mathbf{\Psi }}^T\cdot a$$
那么，
$$\stackrel{\bullet }{a}\cdot \mathbf{\Psi }\cdot b+a\cdot \stackrel{\bullet }{\mathbf{\Psi }}\cdot b+a\cdot \mathbf{\Psi }\cdot \stackrel{\bullet }{b}=\stackrel{\bullet }{b}\cdot ({\mathbf{\Psi }}^T)\cdot a+b\cdot \stackrel{\bullet }{({\mathbf{\Psi }}^T)}\cdot a+b\cdot ({\mathbf{\Psi }}^T)\cdot \stackrel{\bullet }{a}$$
则
$$a\cdot \stackrel{\bullet }{\mathbf{\Psi }}\cdot b=b\cdot \stackrel{\bullet }{({\mathbf{\Psi }}^T)}\cdot a⟹b\cdot \left[{\left(\stackrel{\bullet }{\Psi }\right)}^T-\stackrel{\bullet }{({\mathbf{\Psi }}^T)}\right]\cdot a=0$$
由 $a，b$ 的任意性可知命题成立。故此二者可以不加区分运算的先后次序。另外还表明：对称仿射量的物质导数是对称仿射量，反对称仿射量的物质导数为反对称仿射量。

加速度梯度的反对称部分记为：
$$\mathbf{J}=\frac{1}{2}(a▽-▽a)$$

由于恒等式：
$$2{\mathbf{F}}^T\cdot \mathbf{W}\cdot \mathbf{F}={\mathbf{F}}^T\cdot (\stackrel{\bullet }{\mathbf{F}}\cdot {\mathbf{F}}^{-1}-{\mathbf{F}}^{-T}\cdot \stackrel{\bullet }{\mathbf{F}}{}^T)\cdot \mathbf{F}={\mathbf{F}}^T\cdot \stackrel{\bullet }{\mathbf{F}}-\stackrel{\bullet }{\mathbf{F}}{}^T\cdot \mathbf{F}$$
对左侧求物质导数有：
$$\begin{array}{rl}& 2[\stackrel{\bullet }{\mathbf{F}}{}^T\cdot \mathbf{W}\cdot \mathbf{F}+{\mathbf{F}}^T\cdot \stackrel{\bullet }{\mathbf{W}}\cdot \mathbf{F}+{\mathbf{F}}^T\cdot \mathbf{W}\cdot \stackrel{\bullet }{\mathbf{F}}]\\ & \\ & =2{\mathbf{F}}^T\cdot ({\mathbf{L}}^T\cdot \mathbf{W}+\stackrel{\bullet }{\mathbf{W}}+\mathbf{W}\cdot \mathbf{L})\cdot \mathbf{F}\\ & \\ & =2{\mathbf{F}}^T\cdot (\mathbf{D}\cdot \mathbf{W}+\stackrel{\bullet }{\mathbf{W}}+\mathbf{W}\cdot \mathbf{D})\cdot \mathbf{F}+2{\mathbf{F}}^T\cdot (\frac{{\mathbf{L}}^T-\mathbf{L}}{2}\cdot \mathbf{W}+\mathbf{W}\cdot \frac{\mathbf{L}-{\mathbf{L}}^T}{2})\cdot \mathbf{F}\\ & \\ & =2{\mathbf{F}}^T\cdot (\mathbf{D}\cdot \mathbf{W}+\stackrel{\bullet }{\mathbf{W}}+\mathbf{W}\cdot \mathbf{D})\cdot \mathbf{F}+2{\mathbf{F}}^T\cdot (-{\mathbf{W}}^2+{\mathbf{W}}^2)\cdot \mathbf{F}\\ & \\ & =2{\mathbf{F}}^T\cdot (\mathbf{D}\cdot \mathbf{W}+\stackrel{\bullet }{\mathbf{W}}+\mathbf{W}\cdot \mathbf{D})\cdot \mathbf{F}\end{array}$$
对右侧求物质导数：
$$\begin{array}{rl}& \stackrel{\bullet }{\mathbf{F}}{}^T\cdot \stackrel{\bullet }{\mathbf{F}}+{\mathbf{F}}^T\cdot \stackrel{\bullet \bullet }{\mathbf{F}}-\stackrel{\bullet \bullet }{\mathbf{F}}{}^T\cdot \mathbf{F}-\stackrel{\bullet }{\mathbf{F}}{}^T\cdot \stackrel{\bullet }{\mathbf{F}}\\ & \\ & ={\mathbf{F}}^T\cdot \stackrel{\bullet \bullet }{\mathbf{F}}-\stackrel{\bullet \bullet }{\mathbf{F}}{}^T\cdot \mathbf{F}\\ & \\ & ={\mathbf{F}}^T\cdot (\stackrel{\bullet \bullet }{\mathbf{F}}\cdot {\mathbf{F}}^{-1}-{\mathbf{F}}^{-T}\cdot \stackrel{\bullet \bullet }{\mathbf{F}}{}^T)\cdot \mathbf{F}\\ & \\ & ={\mathbf{F}}^T\cdot (a▽-▽a)\cdot \mathbf{F}\\ & \\ & =2{\mathbf{F}}^T\cdot \mathbf{J}\cdot \mathbf{F}\end{array}$$
比较两式得：
$$\mathbf{J}=\stackrel{\bullet }{\mathbf{W}}+\mathbf{D}\cdot \mathbf{W}+\mathbf{W}\cdot \mathbf{D}$$

3. Lagrange-Cauchy 定理

[Lagrange-Cauchy 定理] \quad 对于加速度为势的梯度的运动，若某一时刻运动是无旋的，则运动始终是无旋的。

证明：由于加速度为势的梯度，则
$$a=▽\alpha (x,t)$$
那么
$$\begin{array}{rl}& (▽\alpha )▽=\frac{\partial }{\partial x^j}\left(\frac{\partial \alpha }{\partial x^i}\right)g^i\otimes g^j\\ & \\ & ▽(▽\alpha )=\frac{\partial }{\partial x^j}\left(\frac{\partial \alpha }{\partial x^i}\right)g^j\otimes g^i=\frac{\partial }{\partial x^i}\left(\frac{\partial \alpha }{\partial x^j}\right)g^i\otimes g^j\\ & \\ & (▽\alpha )▽=▽(▽\alpha )\end{array}$$
则，加速度梯度的反对称部分为
$$\mathbf{J}=\frac{(▽\alpha )▽-▽(▽\alpha )}{2}=0$$
那么
$$\frac{D}{Dt}({\mathbf{F}}^T\cdot \mathbf{W}\cdot \mathbf{F})={\mathbf{F}}^T\cdot (\mathbf{D}\cdot \mathbf{W}+\stackrel{\bullet }{\mathbf{W}}+\mathbf{W}\cdot \mathbf{D})\cdot \mathbf{F}={\mathbf{F}}^T\cdot \mathbf{J}\cdot \mathbf{F}=0$$
若运动某一时刻无旋，即某一时刻 $\mathbf{W}=0$，则运动过程中始终有：
$${\mathbf{F}}^T\cdot \mathbf{W}\cdot \mathbf{F}=0$$
由于 $\mathbf{F}$ 正则，故在运动过程中始终有：
$$\mathbf{W}=0$$
即，运动始终无旋。

命题 \quad 若速度场为势的梯度，则加速度场也为势的梯度。

证明：由于速度有势，则
$$v=▽\phi (x,t)$$
那么
$$a=\stackrel{\bullet }{v}=\stackrel{\bullet }{▽\phi }=(▽\phi {)}^{\prime }+(v▽)\cdot v$$
由于
$$\begin{gathered} (▽\times v)\times v=(v▽-▽v)\cdot v=v▽\cdot v-▽v\cdot v \\ ▽v\cdot v=▽\left(\frac{v\cdot v}{2}\right) \\ ▽\times v=▽\times (▽\phi )=\frac{{\partial }^2\phi }{\partial x^i\partial x^j}{ϵ}^{ijk}{g}_k=0 \end{gathered}$$
故
$$a=▽\left({\phi }^{\prime }+\frac{v\cdot v}{2}\right)$$
证毕。

4. 涡旋传输定理

5. 变形率 $\mathbf{D}$ 与物质旋率 $\mathbf{W}$ 的几何意义

5.1. 线元的物质导数与相关的相对变化率

由于任意线元 $dx$ 的物质导数为：
$$\frac{D}{Dt}(dx)=\frac{D}{Dt}(\mathbf{F}\cdot dX)=\stackrel{\bullet }{\mathbf{F}}\cdot X=(\stackrel{\bullet }{\mathbf{F}}\cdot {\mathbf{F}}^{-1})\cdot (\mathbf{F}\cdot dX)=\mathbf{L}\cdot dx$$

5.1.1. 线元长度的相对变化率

任意线元 $dx$ 长度 $|dx|$ 的物质导数为：
$$\begin{array}{rl}& \frac{D}{Dt}(|dx|)=\frac{D}{Dt}\sqrt{dx\cdot dx}=\frac{1}{2|dx|}\left[\frac{D}{Dt}(dx)\cdot dx+dx\cdot \frac{D}{Dt}(dx)\right]\\ & \\ & =\frac{1}{|dx|}\left(dx\cdot \frac{{\mathbf{L}}^T+\mathbf{L}}{2}\cdot dx\right)=\frac{1}{|dx|}\left(dx\cdot \mathbf{D}\cdot dx\right)=l\cdot (\mathbf{D}\cdot dx)\end{array}$$
可理解为 $\mathbf{D}\cdot dx$ 在 $l$ 上的投影。故沿单位切向量 $l$ 的线元 $dx$ 的长度的相对变化率 $d_l$ 为：
$$d_l\triangleq \frac{\stackrel{\bullet }{|dx|}}{|dx|}=\frac{dx}{|dx|}\cdot \mathbf{D}\cdot \frac{dx}{|dx|}=l\cdot \mathbf{D}\cdot l$$

进一步讨论，在当前时刻什么方向上使得线元的长度相对变化率最大？即求解
$$\{\begin{array}{l}\underset{l}{\max }(l\cdot \mathbf{D}\cdot l)\\ \\ |l|=1\end{array}$$
采用 Lagrange 乘子法进行求解，使得问题化为如下无约束极值问题的必要条件：
$$\{\begin{array}{l}\frac{\partial }{\partial l^k}[l^i{l}^j{D}_{ij}-\eta (l^i{l}^j{g}_{ij}-1)]=0\\ \\ l^i{l}_i-1=0\end{array}$$
即
$$\{\begin{array}{l}(D_{ik}-\eta g_{ik})l^i=0\\ \\ l^i{l}^j{g}_{ij}-1=0\end{array}⟹\{\begin{array}{l}(D_{\bullet k}^i-\eta {\delta }_{\bullet k}^i)l_i=0\\ \\ l^i{l}^j{g}_{ij}-1=0\end{array}$$
其求解等价于求解 $\mathbf{D}$ 的特征值与单位特征向量问题：
$$(\mathbf{D}-\eta \mathbf{I})l=0$$
由于 $\mathbf{D}$ 是对称张量，故存在三个实特征值 $d_{\alpha }(\alpha =1,2,3)$ 并总可以找到两两垂直的单位特征方向 $v_{\alpha }(\alpha =1,2,3)$。且
$$d_{v_{\alpha }}=v_{\alpha }\cdot (\mathbf{D}\cdot v_{\alpha })=v_{\alpha }\cdot (d_{\alpha }{v}_{\alpha })=d_{\alpha }(\alpha =1,2,3)$$

上述讨论实际上说明：在各时刻，线元长度的相对变化率取极值（驻值）的方向为伸长率张量的主方向，长度相对变化率的大小为伸长率张量的对应的驻值。

5.1.2. 线元方向的变化率 (与 W 相关)

考虑线元 $dx$ 单位切向量 $l$ 的物质导数：
$$\begin{array}{rl}& \frac{Dl}{Dt}=\frac{D}{Dt}\left(\frac{dx}{|dx|}\right)=\frac{1}{|dx|}\frac{D}{Dt}(dx)-\frac{dx}{|dx{|}^2}\frac{D}{Dt}(|dx|)\\ & \\ & =\mathbf{L}\cdot \frac{dx}{|dx|}-\frac{dx}{|dx|}(l\cdot \mathbf{D}\cdot l)\\ & \\ & =(\mathbf{L}-d_l\mathbf{I})\cdot l\end{array}$$
其中，$d_l$ 为 $l$ 方向线元长度的相对变化率。又
$$\frac{Dl}{Dt}\cdot l=l\cdot (\mathbf{L}-d_l\mathbf{I})\cdot l=l\cdot (\mathbf{L}-\mathbf{D})\cdot l=l\cdot \mathbf{W}\cdot l=0$$
由于物质旋率张量是反对称张量。这也说明了单位切向量的物质导数与自身垂直的几何含义。

5.1.3. 线元夹角余弦的变化率

在当前时刻，存在两线元 $d{x}_{(1)}$ 与 $d{x}_{(2)}$，其单位切向量分别为：$l_{(1)}、l_{(2)}$，夹角为：$\theta$。那么
$$\begin{array}{rl}& \frac{D}{Dt}(cos\theta )=\frac{D}{Dt}(l_{(1)}\cdot l_{(2)})\\ & \\ & =l_{(1)}\cdot \frac{D{l}_{(2)}}{Dt}+\frac{D{l}_{(1)}}{Dt}\cdot l_{(2)}\\ & \\ & =l_{(1)}\cdot \left(\mathbf{L}-d_{l_{(2)}}\mathbf{I}\right)\cdot l_{(2)}+\left(\mathbf{L}-d_{l_{(1)}}\mathbf{I}\right)\cdot l_{(1)}\cdot l_{(2)}\\ & \\ & =l_{(1)}\cdot \left(\mathbf{L}-d_{l_{(2)}}\mathbf{I}\right)\cdot l_{(2)}+l_{(1)}\cdot \left({\mathbf{L}}^T-d_{l_{(1)}}\mathbf{I}\right)\cdot l_{(2)}\\ & \\ & =2l_{(1)}\cdot \mathbf{D}\cdot l_{(2)}-\left(d_{l_{(1)}}+d_{l_{(1)}}\right)l_{(1)}\cdot l_{(2)}\end{array}$$
特别地，若考虑 $\theta =90°$，则
D D t ( c o s θ ) ∣ θ = 90 ° = − D θ D t ∣ θ = 90 ° = 2 l ⃗ ( 1 ) ⋅ D ⋅ l ⃗ ( 2 ) \left.\dfrac{D}{Dt}(cos\theta)\right|_{\theta=90\degree}=-\left.\dfrac{D\theta}{Dt}\right|_{\theta=90\degree}=2\vec{l}_{(1)}\cdot\bold D\cdot\vec{l}_{(2)} DtD​(cosθ) ​θ=90°​=−DtDθ​ ​θ=90°​=2l (1)​⋅D⋅l (2)​

5.2. 体元的物质导数与与其相对变化率

5.2.1. 体元的物质导数与与其相对变化率

由于
$$\frac{D}{Dt}(dv)=\frac{D}{Dt}(\mathcal{J}d{v}_0)=\stackrel{\bullet }{\mathcal{J}}d{v}_0=▽\cdot v\mathcal{J}d{v}_0=▽\cdot vdv$$
又
$$▽\cdot v=v^i{|}_i=tr(v\otimes ▽)=tr(\mathbf{L})=tr(\mathbf{D})$$
故任意体元 $dv$ 的物质导数为：
$$\frac{D}{Dt}(dv)=tr(\mathbf{D})dv$$
则体元的相对变化率为：
$$\frac{\stackrel{\bullet }{\mathrm{dv}}}{dv}=tr(\mathbf{D})$$

5.2.2 体元物质导数的另一种表达

对任意仿射量 $\mathbf{F}、\mathbf{B}$ 与 实数 $\delta$, 记
$$\hat{\mathbf{F}}=\mathbf{F}+\delta \mathbf{B}$$
则
$$\hat{\mathcal{J}}\triangleq det(\hat{\mathbf{F}})=det(\mathbf{F}+\delta \mathbf{B})=det(\mathbf{F})det(\mathbf{I}+\delta \mathbf{B}{\mathbf{F}}^{-1})$$
进一步有：
D J ^ D δ ∣ δ = 0 = d e t ( F ) t r ( B F − 1 ) \left.\dfrac{D\hat{\mathscr J}}{D\delta}\right|_{\delta=0} =det(\bold F)tr(\bold B\bold F^{-1}) DδDJ^​​ ​δ=0​=det(F)tr(BF−1)
上式运用了特征多项式展开的相关性质。得到上述结论后，令上式中 $\mathbf{F}$ 代表变形梯度，且
$$\delta =t，\mathbf{B}=\stackrel{\bullet }{\mathbf{F}}$$
则有：
D   d e t ( F + t F ∙ ) D t ∣ t = 0 = J t r L = J t r D \left.\dfrac{D\ det(\bold F +t\overset{\bullet}{\bold F})}{Dt}\right|_{t=0}=\mathscr{J}tr{\bold L}=\mathscr{J}tr{\bold D} DtD det(F+tF∙)​ ​t=0​=JtrL=JtrD
则
[ D   d e t ( F + t F ∙ ) D t ∣ t = 0 ] d v 0 = t r ( D ) d v = d v ∙ \left[\left.\dfrac{D\ det(\bold F +t\overset{\bullet}{\bold F})}{Dt}\right|_{t=0}\right]dv_0 =tr(\bold D)dv =\overset{\bullet}{dv} ​DtD det(F+tF∙)​ ​t=0​ ​dv0​=tr(D)dv=dv∙

5.3. 有向面元的物质导数与其面积的相对变化率

由于
$$\begin{array}{rl}& \frac{D}{Dt}(NdS)=\frac{D}{Dt}(\mathcal{J}{\mathbf{F}}^{-T}\cdot {}_{0}Nd{S}_0)=\frac{D}{Dt}(\mathcal{J}{\mathbf{F}}^{-T})\cdot {}_{0}Nd{S}_0\\ & \\ & =[tr(\mathbf{D})(\mathcal{J}{\mathbf{F}}^{-T})+\mathcal{J}(\stackrel{\bullet }{\mathbf{F}}{}^{-1}{)}^T)]\cdot {}_{0}Nd{S}_0\end{array}$$
又
$$\begin{gathered} \frac{D}{Dt}(\mathbf{F}\cdot {\mathbf{F}}^{-1})=\stackrel{\bullet }{\mathbf{F}}\cdot {\mathbf{F}}^{-1}+\mathbf{F}\cdot \stackrel{\bullet }{\mathbf{F}}{}^{-1}=0 \\ \stackrel{\bullet }{\mathbf{F}}{}^{-1}=-{\mathbf{F}}^{-1}\cdot \stackrel{\bullet }{\mathbf{F}}\cdot {\mathbf{F}}^{-1}=-{\mathbf{F}}^{-1}\cdot \mathbf{L} \end{gathered}$$
则任意有向线元的物质导数为：
$$\frac{D}{Dt}(NdS)=[tr(\mathbf{D})(\mathcal{J}{\mathbf{F}}^{-T})-{\mathbf{L}}^T\cdot (\mathcal{J}{\mathbf{F}}^{-T})]\cdot {}_{0}Nd{S}_0=[tr(\mathbf{D})\mathbf{I}-{\mathbf{L}}^T]\cdot NdS$$
进一步，微元面积的物质导数为：
$$\begin{array}{rl}& \frac{DdS}{Dt}=\frac{D}{Dt}\sqrt{NdS\cdot NdS}\\ & \\ & =\frac{1}{2dS}\left\{[tr(\mathbf{D})\mathbf{I}-{\mathbf{L}}^T]\cdot NdS\cdot NdS+NdS\cdot [tr(\mathbf{D})\mathbf{I}-{\mathbf{L}}^T]\cdot NdS\right\}\\ & \\ & =\frac{1}{2dS}\left\{NdS\cdot [tr(\mathbf{D})\mathbf{I}-\mathbf{L}]\cdot NdS+NdS\cdot [tr(\mathbf{D})\mathbf{I}-{\mathbf{L}}^T]\cdot NdS\right\}\\ & \\ & =\frac{1}{dS}\left[tr(\mathbf{D})(dS{)}^2-NdS\cdot \frac{\mathbf{L}+{\mathbf{L}}^T}{2}\cdot NdS\right]\\ & \\ & =dS\left[tr(\mathbf{D})-N\cdot \mathbf{D}\cdot N\right]\end{array}$$
那么，任意有向面元 $NdS$ 面积的相对变化率为：
$$\frac{\stackrel{\bullet }{\mathrm{dS}}}{dS}=tr(\mathbf{D})-N\cdot \mathbf{D}\cdot N=tr(\mathbf{D})-d_N$$

5.4. 物质旋率张量 W 的几何意义

当前时刻，伸长率张量 $\mathbf{D}$ 的单位特征向量为： $v_{\alpha }(\alpha =1,2,3)$ ，对应的特征值分别为：$d_{\alpha }(\alpha =1,2,3)$，考虑沿当前 $v_{\alpha }$ 方向的线元的方向变化率：
$$\stackrel{\bullet }{l}(v_{\alpha })=\stackrel{\bullet }{l}{|}_{l=v_{\alpha }}$$
需要注意的是：
$$\stackrel{\bullet }{l}(v_{\alpha })\ne \stackrel{\bullet }{v_{\alpha }}$$
因为对线元的单位方向的物质导数表示研究的始终是同一线元，但前一时刻与下一时刻伸长率张量的特征向量可能与不同线元相重合。
l ⃗ ∙ ( v ⃗ α ) = ( L − d α I ) ⋅ v ⃗ α = ( L − D ) ⋅ v ⃗ α = W ⋅ v ⃗ α = ω ⃗ × v ⃗ α \overset{\bullet}{\vec{l}}(\vec{v}_\alpha) =(\bold L-d_{\alpha}\bold I)\cdot\vec{v}_{\alpha} =(\bold L-\bold D)\cdot\vec{v}_{\alpha} =\bold W\cdot\vec{v}_{\alpha} =\vec{\omega}\times\vec{v}_{\alpha} l ∙(v α​)=(L−dα​I)⋅v α​=(L−D)⋅v α​=W⋅v α​=