张量与散度

张量的几个物理场景

应力状态

材料内部受力状况就是应力状态，应力在不同方向的大小不同，并且与感应面有关。如图所示材料某点的应力状态，可用矩阵表示
$$\sigma =\left(\begin{array}{ccc}{\sigma }_x& {\tau }_{xy}& {\tau }_{xz}\\ {\tau }_{yx}& {\sigma }_y& {\tau }_{yz}\\ {\tau }_{zx}& {\tau }_{zy}& {\sigma }_z\end{array}\right)$$
这是个对称矩阵，即${\tau }_{xy}={\tau }_{yx},{\tau }_{xz}={\tau }_{zx},{\tau }_{yz}={\tau }_{zy}$
现在已知某任意平面，那么怎么求该平面上的应力状态呢?如下图所示，以平面示意
红色的箭头是该平面的总的应力，该总应力可以分解成垂直于平面的正应力和平行于平面的切应力。若已知平面的法矢为$n=(n_x,n_y,n_z)$，其总应力计算如下
$${\sigma }^{\prime }=n\cdot \sigma =(n_x,n_y,n_z)\left(\begin{array}{ccc}{\sigma }_x& {\tau }_{xy}& {\tau }_{xz}\\ {\tau }_{yx}& {\sigma }_y& {\tau }_{yz}\\ {\tau }_{zx}& {\tau }_{zy}& {\sigma }_z\end{array}\right)$$
那么正应力为
$$\begin{gathered} {\sigma }_n^{\prime }={\sigma }^{\prime }/n\cdot n={\sigma }^{\prime }\cdot n\cdot n \\ =(n_x,n_y,n_z)\left(\begin{array}{ccc}{\sigma }_x& {\tau }_{xy}& {\tau }_{xz}\\ {\tau }_{yx}& {\sigma }_y& {\tau }_{yz}\\ {\tau }_{zx}& {\tau }_{zy}& {\sigma }_z\end{array}\right)\left(\begin{array}{c}{n}_x\\ n_y\\ n_z\end{array}\right)(n_x,n_y,n_z) \end{gathered}$$
切应力的为
$${\sigma }_{\tau }^{\prime }={\sigma }^{\prime }÷n={\sigma }^{\prime }-{\sigma }_n^{\prime }$$

转动惯量

一个物体的转动，和平移一样，也存在惯性。转动惯量就是衡量这个惯性的大小。对于规则的物体，其转动比较容易求得，如下图所示
该立方块绕$x$轴的转动惯量为$J_x$，绕$y$轴的转动惯量为$J_y$，绕$z$轴的转动惯量为$J_z$。那么如何求得其绕任意轴$n$的转动惯量$J^{\prime }$，如图红色所示。转动惯量也是一种张量，其数学形式如下
$$J=\left(\begin{array}{ccc}{J}_x& 0& 0\\ 0& J_y& 0\\ 0& 0& J_z\end{array}\right)$$
那么绕$n$的总转动惯量计算如下
$$J^{\prime }=n\cdot J=(n_x,n_y,n_z)\left(\begin{array}{ccc}{J}_x& 0& 0\\ 0& J_y& 0\\ 0& 0& J_z\end{array}\right)$$
总的转动惯量到转轴的投影就是绕这个轴的转动惯量
$$J_n^{\prime }=J^{\prime }/n\cdot n$$
总的转动惯量垂直于这个轴的投影为
$$J_{\tau }^{\prime }=J^{\prime }÷n$$
有个疑问，为什么会有这个垂直的转动惯量。这是因为，在物体转动的时候，两端会甩开远离轴，或者聚拢越来越接近轴，而这个转动的转动惯量就是总转动惯量垂直于这个轴的分量。

各向异性的感应电场

我们知道，一个电荷在真空中形成的电场是向各个方向发散，并且每个方向等强度，即等势面是球形。但如果在物质内部呢?这个电荷在物质内部形成的是感应电场(也常称为电位移矢量)，并且由于很多物质都是各向异性的，也就是说电荷在不同方向的感应效果不一样，所以等势面试是椭球，如下图所示

其实所有三维张量都可以用椭球表示，二维张量可以用椭圆表示。
电荷在真空发出的电场为
$$E=\frac{Q}{4\pi r^2}\left(\begin{array}{ccc}ϵ& 0& 0\\ 0& ϵ& 0\\ 0& 0& ϵ\end{array}\right)$$
电荷在各向异性物质内感应的电场为
$$D=\frac{Q}{4\pi r^2}\left(\begin{array}{ccc}{ϵ}_{xx}& {ϵ}_{xy}& {ϵ}_{xz}\\ {ϵ}_{yx}& {ϵ}_{yy}& {ϵ}_{yz}\\ {ϵ}_{zx}& {ϵ}_{zy}& {ϵ}_{zz}\end{array}\right)$$
该式中的矩阵为对称矩阵，即${ϵ}_{xy}={ϵ}_{yx},{ϵ}_{xz}={ϵ}_{zx},{ϵ}_{yz}={ϵ}_{zy}$。矩阵中${ϵ}_{ij}={ϵ}_{ij}^rϵ$，其中${ϵ}_{ij}^r$为方向$ij$的相对介电系数，所以任意方向的感应电场强度为
$$D^{\prime }=n\cdot D=(n_x,n_y,n_z)\left(\begin{array}{ccc}{ϵ}_{xx}& {ϵ}_{xy}& {ϵ}_{xz}\\ {ϵ}_{yx}& {ϵ}_{yy}& {ϵ}_{yz}\\ {ϵ}_{zx}& {ϵ}_{zy}& {ϵ}_{zz}\end{array}\right)$$
其大小为
$$|D^{\prime }|=D^{\prime }/n=D^{\prime }\cdot n=(n_x,n_y,n_z)\left(\begin{array}{ccc}{ϵ}_{xx}& {ϵ}_{xy}& {ϵ}_{xz}\\ {ϵ}_{yx}& {ϵ}_{yy}& {ϵ}_{yz}\\ {ϵ}_{zx}& {ϵ}_{zy}& {ϵ}_{zz}\end{array}\right)\left(\begin{array}{c}{n}_x\\ n_y\\ n_z\end{array}\right)$$
这就是一个二次型，其图像就是椭球，所以说张量可以用椭球表示。

张量的散度

我们知道张量可以表示感应电场，其实他可以表示任意种不同方向不同大小和方向的场，也就是张量场。既然有场，那就有散度的概念。已知下面的张量
$$M=\left(\begin{array}{ccc}{a}_{11}& a_{12}& a_{13}\\ a_{21}& a_{22}& a_{23}\\ a_{31}& a_{32}& a_{33}\end{array}\right)$$
那么该场为
$$f=(x,y,z)\left(\begin{array}{ccc}{a}_{11}& a_{12}& a_{13}\\ a_{21}& a_{22}& a_{23}\\ a_{31}& a_{32}& a_{33}\end{array}\right)$$
散度有两种计算方式，一种是极限的计算方式
$$divf=\underset{v\to 0}{\lim }{\oint }_{s}fds/v$$
还有一种方式为$$divf=\nabla \cdot f=\frac{\partial f}{\partial x}+\frac{\partial f}{\partial y}+\frac{\partial f}{\partial z}$$
采用第二种比较容易计算，那么张量场的散度为
$$divf=\nabla \cdot f=a_{11}+a_{22}+a_{33}$$
也就是说张量的散度就是张量的迹，可知一个方阵无论经过什么样的相似变换，它的迹不变；也就是说场无论经过什么样的坐标变换，它的散度不变。