西安科技大学

安全工程专业课程

力学的数学基础


金洪伟
西安科技大学安全科学与工程学院

返回目录⇡

如何浏览?

  1. 从浏览器地址栏打开 https://zimo.net/math/
  2. 点击章节列表中的任一链接,打开相应的演示文稿;
  3. 点击链接打开演示文稿,使用空格键或方向键导航;
  4. f键进入全屏播放,再按Esc键退出全屏;
  5. Alt键同时点击鼠标左键进行局部缩放;
  6. Esco键进入幻灯片浏览视图。

请使用最新版本浏览器访问此演示文稿以获得更好体验。

向量和向量分析

目 录

  1. 向量
  2. 向量代数
  3. 向量的微分运算
  4. 向量的积分运算

第1节 向量

1.1 力学中的“量”

在数学、物理和力学中,经常涉及以下三种基本的物理量:

标量(Scalar)
又称纯量,是只有大小,没有方向的量。像质量、时间、密度、温度、速率、长度、面积、体积、静水压力、能量、电阻、弹性模量等都是标量。
向量(Vector)
又称矢量,是同时具有大小和方向,并能按照向量代数和其他向量相加的几何对象。像位移、速度、力、动量、动量矩、力矩、电矩、磁矩、电场强度、磁场强度、温度梯度、压力梯度等都是向量。
张量(Tensor)
是描述向量、标量和其他张量之间线性关系的几何对象。张量特别适宜于表示各向异性介质的特征。像应力、应变、应变率、惯性矩、渗透率、弹性系数张量、位移梯度、应变梯度、应力梯度、电极化率等都是(二阶)张量。

1.1 力学中的“量”

在三维欧几里得空间中,可以用 1 个实数表示一个标量,用 3 个实数构成的数组表示一个向量,用 9 个实数构成的方阵表示一个(二阶)张量。在表示三种量的分量时,标量不需要下标(自由标),向量需要 1 个下标,表示成方阵的张量需要 2 个下标。

有时将标量、向量和张量统称为张量,其分量下标的个数称为张量的阶数(和矩阵秩的概念不一样)。从而标量为 0 阶张量,向量为 1 阶张量,而像应力、应变之类的张量为 2 阶张量。

1.2 向量的表示

在三维欧几里得空间中,可以使用三个数标量表示向量v,可以多种形式书写:

从点 AB 的向量

$$\boldsymbol{v}=\left( v_1,v_2,v_3 \right)$$

$$\boldsymbol{v}=\left( v_x,v_y,v_z \right)$$

$$\boldsymbol{v}=\boldsymbol{v}_1+\boldsymbol{v}_2+\boldsymbol{v}_3=v_1\boldsymbol{\hat{e}}_1+v_2\boldsymbol{\hat{e}}_2+v_3\boldsymbol{\hat{e}}_3$$

上式中,$\boldsymbol{\hat{e}}_1$, $\boldsymbol{\hat{e}}_2$和$\boldsymbol{\hat{e}}_3$分别是与三个坐标轴同方向的单位向量:

$$\boldsymbol{\hat{e}}_1=\left( 1,0,0 \right) ,\,\boldsymbol{\hat{e}}_2=\left( 0,1,0 \right) ,\,\boldsymbol{\hat{e}}_3=\left( 0,0,1 \right)$$

1.2 向量的表示

为了便于使向量参与矩阵运算,也常将向量定义为如下的列向量(只有在少数时候看作是行向量):

$$\boldsymbol{v}=\left[ \begin{array}{c} v_1\\ v_2\\ v_3\\ \end{array} \right] =\left[ \begin{matrix} v_1& v_2& v_3\\ \end{matrix} \right] ^{\text{T}}$$

以上描述向量的三个数 v1v2v3 (或 vxvy vz)称为向量在该坐标系统三个坐标轴方向上的分量

第2节 向量代数

2.1 向量加减

对于如下两个向量:

$$\boldsymbol{a}=a_1\boldsymbol{\hat{e}}_1+a_2\boldsymbol{\hat{e}}_2+a_3\boldsymbol{\hat{e}}_3,\,\boldsymbol{b}=b_1\boldsymbol{\hat{e}}_1+b_2\boldsymbol{\hat{e}}_2+b_3\boldsymbol{\hat{e}}_3$$

向量加法

$$\boldsymbol{a}+\boldsymbol{b}=\left( a_1+b_1 \right) \boldsymbol{\hat{e}}_1+\left( a_2+b_2 \right) \boldsymbol{\hat{e}}_2+\left( a_3+b_3 \right) \boldsymbol{\hat{e}}_3$$

向量减法

$$\boldsymbol{a}-\boldsymbol{b}=\left( a_1-b_1 \right) \boldsymbol{\hat{e}}_1+\left( a_2-b_2 \right) \boldsymbol{\hat{e}}_2+\left( a_3-b_3 \right) \boldsymbol{\hat{e}}_3$$

2.1 向量的加法和减法

这种求两向量之和的方法称为向量相加的平行四边形法则

向量加法
向量减法

应用:从点 $\boldsymbol{a}$ 到点 $\boldsymbol{b}$ 的位移为$\boldsymbol{b}-\boldsymbol{a}$。

2.2 向量的标量乘法

向量可以与一个标量(实数)相乘,其乘积也是一个向量

$$r\boldsymbol{a}=\left( ra_1 \right) \boldsymbol{\hat{e}}_1+\left( ra_2 \right) \boldsymbol{\hat{e}}_2+\left( ra_3 \right) \boldsymbol{\hat{e}}_3$$

向量乘以 -1 和 2

向量的加减法与标量乘法统称为向量的线性运算

2.2 向量的标量乘法

牛顿第二运动定律:

施加于物体的外力等于此物体的质量与加速度的乘积。

可用向量法表示

$$\boldsymbol{F}=m\boldsymbol{a}$$

其中

$$\boldsymbol{a}=\frac{\text{d}\boldsymbol{v}}{\text{d}t}=\left( \frac{\text{d}}{\text{d}t} \right) \left( \frac{\text{d}\boldsymbol{r}}{\text{d}t} \right) =\frac{\text{d}^2\boldsymbol{r}}{\text{d}t^2}$$

2.3 向量的长度

向量的长度又称为向量的,用$\lVert \boldsymbol{a} \rVert$或$\left| \boldsymbol{a} \right|$(容易和绝对值运算混淆)表示。

$$\lVert \boldsymbol{a} \rVert =\sqrt{a_{1}^{2}+a_{2}^{2}+a_{3}^{2}}$$

任何长度为 1 的向量称为单位向量,常用单位向量表示方向。常在向量的上方加一个帽子表示单位向量,如$\boldsymbol{\hat{a}}$。求任意向量对应的单位向量的过程称为向量的标准化,其计算方法为

$$ \boldsymbol{\hat{a}}=\frac{\boldsymbol{a}}{\lVert \boldsymbol{a} \rVert}=\frac{a_1}{\lVert \boldsymbol{a} \rVert}\boldsymbol{\hat{e}}_1+\frac{a_2}{\lVert \boldsymbol{a} \rVert}\boldsymbol{\hat{e}}_2+\frac{a_3}{\lVert \boldsymbol{a} \rVert}\boldsymbol{\hat{e}}_3=\cos \alpha \boldsymbol{\hat{e}}_1+\cos \beta \boldsymbol{\hat{e}}_2+\cos \gamma \boldsymbol{\hat{e}}_3 $$

$\alpha$、$\beta$ 和 $\gamma$ 分别为向量 $\boldsymbol{a}$ 与三个坐标轴的夹角。

长度为 0 的向量称为零向量,表示为 $\boldsymbol{0}$。

2.4 向量的标量积

向量的的标量积又称内积点积数量积,记作$\boldsymbol{a}\cdot \boldsymbol{b}$。标量积是一个标量

$$\boldsymbol{a}\cdot \boldsymbol{b}=\lVert \boldsymbol{a} \rVert \lVert \boldsymbol{b} \rVert \cos \theta =a_1b_1+a_2b_2+a_3b_3$$

式中,$\theta$ 是两向量的夹角。

其几何意义相等于 $\boldsymbol{b}$ 在 $\boldsymbol{a}$ 上的投影再乘以 $\boldsymbol{a}$ 的长度。

应用

  • 功 $W$ 是力向量 $\boldsymbol{F}$ 和位移向量 $\boldsymbol{x}$ 的标量积,即 $W=\boldsymbol{F}\cdot \boldsymbol{x}$。
  • 功率 $P$ 是力向量 $\boldsymbol{F}$ 和速度向量 $\boldsymbol{v}$ 的标量积,即 $P=\boldsymbol{F}\cdot \boldsymbol{v}$。

2.4 向量的向量积

向量的的向量积又称外积叉积,记作$\boldsymbol{a}\times \boldsymbol{b}$。向量积是一个向量

$$ \begin{aligned} \boldsymbol{a}\times \boldsymbol{b}&=\lVert \boldsymbol{a} \rVert \lVert \boldsymbol{b} \rVert \sin \left( \theta \right) \boldsymbol{n}\\ &=\left( a_2b_3-a_3b_2 \right) \boldsymbol{\hat{e}}_1+\left( a_3b_1-a_1b_3 \right) \boldsymbol{\hat{e}}_2+\left( a_1b_2-a_2b_1 \right) \boldsymbol{\hat{e}}_3\\ \end{aligned} $$

按右手法则确定向量积方向

式中,$\theta$ 是两向量的夹角;$\boldsymbol{n}$ 是按右手法则垂直于 $\boldsymbol{a}$ 和 $\boldsymbol{b}$ 所在平面的单位向量(如右图)。

所得的向量与 $\boldsymbol{a}$ 和 $\boldsymbol{b}$ 所在平面垂直,方向由右手法则规定,大小是两个被乘向量张成的平行四边形的面积。

2.4 向量的向量积

应用

  • 力矩 $\boldsymbol{M}$ 等于径向向量 $\boldsymbol{r}$ 与作用力 $\boldsymbol{F}$ 的向量积

    $$\boldsymbol{M}=\boldsymbol{r}\times \boldsymbol{F}$$

  • 角动量 $\boldsymbol{L}$ 等于径向向量 $\boldsymbol{r}$ 与动量向量 $\boldsymbol{p}$ 的向量积

    $$\boldsymbol{M}=\boldsymbol{r}\times \boldsymbol{F}=\boldsymbol{r}\times \boldsymbol{p}=\boldsymbol{r}\times \left( m\boldsymbol{v} \right) =\boldsymbol{r}\times \left( \boldsymbol{\omega }\times \left( m\boldsymbol{r} \right) \right) =mr^2\boldsymbol{\omega }=I\boldsymbol{\omega }$$

    式中,$I$ 表示质点的转动惯量, $\boldsymbol {\omega }$ 是角速度向量。

2.5 向量的标量三重积

由三个向量定义的平行六面体

向量的的标量三重积又称混合积框积三重标量积,记作$\boldsymbol{a}\cdot \left( \boldsymbol{b}\times \boldsymbol{c} \right)$或 $\left[ \boldsymbol{abc} \right]$。

$$\boldsymbol{a}\cdot \left( \boldsymbol{b}\times \boldsymbol{c} \right) =\det \left[ \begin{matrix} a_1& a_2& a_3\\ b_1& b_2& b_3\\ c_1& c_2& c_3\\ \end{matrix} \right] $$

其几何意义是三向量所构成的平行六面体的体积。三个向量的混合积具有循环轮换性,即

$$\left[ \boldsymbol{abc} \right] =\boldsymbol{a}\cdot \left( \boldsymbol{b}\times \boldsymbol{c} \right) =\boldsymbol{b}\cdot \left( \boldsymbol{c}\times \boldsymbol{a} \right) =\boldsymbol{c}\cdot \left( \boldsymbol{a}\times \boldsymbol{b} \right)$$

2.5 向量的向量三重积

向量的的向量三重积记作$\boldsymbol{a}\times \left( \boldsymbol{b}\times \boldsymbol{c} \right)$。对于任意向量,总有以下恒等式

$$\boldsymbol{a}\times \left( \boldsymbol{b}\times \boldsymbol{c} \right) =\left( \boldsymbol{a}\cdot \boldsymbol{c} \right) \boldsymbol{b}-\left( \boldsymbol{a}\cdot \boldsymbol{b} \right) \boldsymbol{c}$$

上式称作三重积展开拉格朗日公式。根据该式可得雅可比恒等式

$$\boldsymbol{a}\times \left( \boldsymbol{b}\times \boldsymbol{c} \right) +\boldsymbol{b}\times \left( \boldsymbol{c}\times \boldsymbol{a} \right) +\boldsymbol{c}\times \left( \boldsymbol{a}\times \boldsymbol{b} \right) =0$$

另外一个非常有用的式子

$$\left( \boldsymbol{a}\times \boldsymbol{b} \right) \times \boldsymbol{c}=\boldsymbol{a}\times \left( \boldsymbol{b}\times \boldsymbol{c} \right) -\boldsymbol{b}\times \left( \boldsymbol{a}\times \boldsymbol{c} \right)$$

以及一个和向量分析中梯度相关的恒等式

$$\nabla \times \left( \nabla \times \boldsymbol{f} \right) =\nabla \left( \nabla \cdot \boldsymbol{f} \right) -\left( \nabla \cdot \nabla \right) \boldsymbol{f}$$

2.6 多个笛卡尔坐标系之间的变换

设向量 $\boldsymbol{a}$ 在以 $\left\{ \boldsymbol{\hat{e}}_1,\boldsymbol{\hat{e}}_2,\boldsymbol{\hat{e}}_3 \right\}$ 为基的笛卡尔坐标系下表示为

$$\boldsymbol{a}=p\boldsymbol{\hat{e}}_1+q\boldsymbol{\hat{e}}_2+r\boldsymbol{\hat{e}}_3$$

同样的 $\boldsymbol{a}$ 在以 $\left\{ \boldsymbol{\hat{n}}_1,\boldsymbol{\hat{n}}_2,\boldsymbol{\hat{n}}_3 \right\}$ 为基的笛卡尔坐标系下表示为

$$\boldsymbol{a}=u\boldsymbol{\hat{n}}_1+v\boldsymbol{\hat{n}}_2+w\boldsymbol{\hat{n}}_3$$

则有如下关系成立

$$\left[ \begin{array}{c} u\\ v\\ w\\ \end{array} \right] =\left[ \begin{matrix} c_{11}& c_{12}& c_{13}\\ c_{21}& c_{22}& c_{23}\\ c_{31}& c_{32}& c_{33}\\ \end{matrix} \right] \left[ \begin{array}{c} p\\ q\\ r\\ \end{array} \right]$$

2.6 多个笛卡尔坐标系之间的变换

上页式中

$$C=\left[ \begin{matrix} c_{11}& c_{12}& c_{13}\\ c_{21}& c_{22}& c_{23}\\ c_{31}& c_{32}& c_{33}\\ \end{matrix} \right] =\left[ \begin{matrix} \boldsymbol{\hat{n}}_1\cdot \boldsymbol{\hat{e}}_1& \boldsymbol{\hat{n}}_1\cdot \boldsymbol{\hat{e}}_2& \boldsymbol{\hat{n}}_1\cdot \boldsymbol{\hat{e}}_3\\ \boldsymbol{\hat{n}}_2\cdot \boldsymbol{\hat{e}}_1& \boldsymbol{\hat{n}}_2\cdot \boldsymbol{\hat{e}}_2& \boldsymbol{\hat{n}}_2\cdot \boldsymbol{\hat{e}}_3\\ \boldsymbol{\hat{n}}_3\cdot \boldsymbol{\hat{e}}_1& \boldsymbol{\hat{n}}_3\cdot \boldsymbol{\hat{e}}_2& \boldsymbol{\hat{n}}_3\cdot \boldsymbol{\hat{e}}_3\\ \end{matrix} \right]$$

称 $C$ 为从 $e$ 到 $n$ 的变换矩阵旋转矩阵方向余弦矩阵

向量的分量值虽然随着坐标系的变化而变换,但这种变换符合一定的规则,是可逆的,保证不同坐标系下的向量仍为同一向量。

2.6 多个笛卡尔坐标系之间的变换

在二维坐标中,若将坐标系 $\left\{ \boldsymbol{\hat{e}}_1,\boldsymbol{\hat{e}}_2 \right\}$ 逆时针旋转 $\theta$ 变为坐标系 $\left\{ \boldsymbol{\hat{n}}_1,\boldsymbol{\hat{n}}_2 \right\}$,则有

$$C=\left[ \begin{matrix} \cos \theta& \sin \theta\\ -\sin \theta& \cos \theta\\ \end{matrix} \right]$$

第3节 向量的微分运算

3.1 场和场函数

在物理中,我们所关注的物理量往往是在待求解的空间区域分布的,这个和特定物理量相关的区域就是一个。如温度场、速度场、应变场等等。相应地,确定场的物理量称为场量,定义位置点(有时还要加上时间)与场量关系的函数称为场函数

3.1 场和场函数

1. 标量场

假如一个空间中的每一点的属性都可以以一个标量来代表的话,那么这个场就是一个标量场。平面内的标量场通常用不同颜色表示大小,有时候会画出等势线。

标量场云图
电势场等势线图

一个简单的电势场,两边各有一个点电荷

3.1 场和场函数

2. 向量场

假如一个空间中的每一点的属性都可以以一个向量代表的话,那么这个场就是一个向量场。平面内的向量场通常用具有不同长度和方向的箭头表示。

向量场$\left( \sin y,\sin x \right)$
磁铁的磁力线 (磁偶极子)

3.1 场和场函数

3. 场函数表述方法

标量场函数(或称标量函数)可表示为

$$f=f\left( \boldsymbol{x} \right) $$

向量场函数(或称向量值函数)可表示为

$$\boldsymbol{f}=\boldsymbol{f}\left( \boldsymbol{x} \right)$$

:向量场函数的自变量、因变量都是向量。

3.2 梯度

1. 什么是梯度?

设有这样一个定义在平面笛卡尔坐标系下的标量场函数

$$f = f\left( {x,y} \right)$$

该二元函数的图形是一个空间曲面。我们在高等数学中已经学习过该函数的偏导数全微分方向导数梯度的定义和几何意义,现在先来回顾一下。

1. 什么是梯度?

设一平面与该曲面在某一点\({P_0}\left( {{x_0},{y_0}} \right)\)相切,该切平面的方程为

$$f-f_0=f_x\left( x_0,y_0 \right) \left( x-x_0 \right) +f_y\left( x_0,y_0 \right) \left( y-y_0 \right)$$

函数的全微分

$$\text{d}f=\frac{\partial f}{\partial x}\text{d}x+\frac{\partial f}{\partial y}\text{d}y$$

表示该曲面在点 ${P_0}\left( {{x_0},{y_0}} \right)$ 处的切平面上点的竖坐标 f 的增量,即 f 的变化主要是由 xy 的变化引起的。

1. 什么是梯度?

函数的方向导数

$$\left. \frac{\partial f}{\partial l} \right|_{\left( x_0,y_0 \right)}=f_x\left( x_0,y_0 \right) \cos \alpha +f_y\left( x_0,y_0 \right) \cos \beta$$

它表示函数在点 ${P_0}\left( {{x_0},{y_0}} \right)$ 沿着处于水平面上的某方向 $\boldsymbol{l}$ 的变化率。其中 $\cos \alpha $, $\cos \beta $ 是方向 $\boldsymbol{l}$ 与坐标轴夹角的余弦。

偏导数 $\frac{\partial f}{\partial x}$ 和 $\frac{\partial f}{\partial y}$ 反映的是函数沿坐标轴方向的变化率。

1. 什么是梯度?

函数的梯度为

$$\boldsymbol{grad}f\left( x,y \right) =\nabla f=\frac{\partial f}{\partial x}\boldsymbol{\hat{e}}_x+\frac{\partial f}{\partial y}\boldsymbol{\hat{e}}_y$$

梯度是一个矢量,它与方向导数的关系为

$$\frac{\partial f}{\partial l}=\boldsymbol{grad}f\left( x,y \right) \cdot \left( \cos \alpha \boldsymbol{\hat{e}}_x+\sin \alpha \boldsymbol{\hat{e}}_y \right)$$

即方向导数 $\frac{\partial f}{\partial l}$ 是梯度 $\boldsymbol{grad}f\left( x,y \right) $ 在方向 $\boldsymbol{l}$ 上的投影。从而可知,梯度总是大于等于方向导数。进一步来说,函数在某点的梯度是这样一个矢量,它的方向与取得最大方向导数的方向一致,而它的模为方向导数的最大值。

2. $\nabla$ 算符

梯度表达式 $\nabla f$ 中的 $\nabla$ 叫做哈密顿(Hamilton)算符(运算符,算子),或称为Nabla符号,有时可直接读作 “del”,或直接根据其含义说成为“梯度”或 “grad”。这个符号是一个矢性微分算符,其定义为

$$\nabla =\frac{\partial}{\partial x_1}\boldsymbol{\hat{e}}_1+\frac{\partial}{\partial x_2}\boldsymbol{\hat{e}}_2+\frac{\partial}{\partial x_3}\boldsymbol{\hat{e}}_3$$

2. $\nabla$ 算符

$\nabla$ 算符本身并无意义,它同微分算符 $\text{d}f$、不定积分算符 $\int{f}$ 及偏微分算符 $\frac{\partial f}{\partial x}$ 一样,只是一种微分运算符号。但这是一个神奇的符号,它在运算中具有矢量和微分算符的双重性质,分别可与数量场和矢量场发生作用。

把任何标量场函数放到 $\nabla$ 的右边,就求得函数的梯度,如前面得出的梯度 $\nabla f$。

将函数 $f\left( x,y \right) =−\left( \cos 2x+\cos 2y \right) ^2$ 的梯度描绘为在底面上投影的向量场

3. 标量场的梯度

梯度是一个从标量场中得出的向量(或称向量值函数)。多元可微标量值函数 $f$ 在点 $P$ 处的梯度,是以 $f$ 在 $P$ 处的偏导数为分量的向量。

如果多元函数在点 $P$ 处的梯度不是零向量,它的方向是这个函数在 $P$ 处最大增长的方向,而它的模是在这个方向上的增长率。

$$\boldsymbol{grad}f=\nabla f=\frac{\partial f}{\partial x_1}\boldsymbol{\hat{e}}_1+\frac{\partial f}{\partial x_2}\boldsymbol{\hat{e}}_2+\frac{\partial f}{\partial x_3}\boldsymbol{\hat{e}}_3$$

4. 向量场的导数

向量场函数 $\boldsymbol{F}\left( \boldsymbol{x} \right)$ 的导数为

$$\boldsymbol{F}'\left( \boldsymbol{x} \right) =\left[ \begin{matrix} \frac{\partial F_1}{\partial x_1}& \frac{\partial F_1}{\partial x_2}& \frac{\partial F_1}{\partial x_3}\\ \frac{\partial F_2}{\partial x_1}& \frac{\partial F_2}{\partial x_2}& \frac{\partial F_2}{\partial x_3}\\ \frac{\partial F_3}{\partial x_1}& \frac{\partial F_3}{\partial x_2}& \frac{\partial F_3}{\partial x_3}\\ \end{matrix} \right]$$

该矩阵称作雅可比矩阵,常用符号 $\boldsymbol{J}$ 表示,其分量可表示为

$$\boldsymbol{J}_{ij}=\frac{\partial F_i}{\partial x_j}$$

5. 向量场的梯度

向量场函数 $\boldsymbol{F}\left( \boldsymbol{x} \right)$ 的梯度为

$$\nabla \boldsymbol{F}=\boldsymbol{J}'=\left[ \begin{matrix} \frac{\partial F_1}{\partial x_1}& \frac{\partial F_2}{\partial x_1}& \frac{\partial F_3}{\partial x_1}\\ \frac{\partial F_1}{\partial x_2}& \frac{\partial F_2}{\partial x_2}& \frac{\partial F_3}{\partial x_2}\\ \frac{\partial F_1}{\partial x_3}& \frac{\partial F_2}{\partial x_3}& \frac{\partial F_3}{\partial x_3}\\ \end{matrix} \right]$$

即向量场函数的梯度等于其导数(雅可比矩阵)的转置。

3.3 散度

1. 散度的定义

设 $\boldsymbol{F}$ 是向量场函数,将 $\nabla$ 看作向量,根据向量的标量积的计算方法可得

$$\nabla \cdot \boldsymbol{F}=\frac{\partial F_1}{\partial x_1}+\frac{\partial F_2}{\partial x_2}+\frac{\partial F_3}{\partial x_3}$$

$\nabla \cdot \boldsymbol{F}$ 在物理学中被称作向量场 $\boldsymbol{F}$ 的 散度,即

$$\text{div}\boldsymbol{F}=\nabla \cdot \boldsymbol{F}$$

散度的定义

如右图所示,在向量场 $\boldsymbol{F}\left( \boldsymbol{x} \right)$ 中有任意点 $\boldsymbol{x}$,$S_1$, $S_2$, $S_3$, … 是一系列包括点 $\boldsymbol{x}$ 的封闭曲面,其体积记作 $V_i$,则在点 $\boldsymbol{x}$ 处的向量场的散度为

$$\left. \text{div}\boldsymbol{F} \right|_{\boldsymbol{x}}=\lim_{V_i\rightarrow 0} \frac{\varPhi \left( S_i \right)}{|V_i|}$$

$$\left. \text{div}\boldsymbol{F} \right|_{\boldsymbol{x}}=\lim_{V\rightarrow 0} \frac{1}{|V|}\oiint_{S\left( V \right)}{S\left( V \right) \boldsymbol{F}\cdot \boldsymbol{\hat{n}}\,dS}$$

$\varPhi \left( S_i \right)$ 是通过曲面 $S_i$ 的通量,$\boldsymbol{\hat{n}}$ 是封闭曲面上任意面元 $\text{d}a$ 的法向单位向量。

2. 散度的物理意义

向量场的散度

向量场的散度表示向量场通量在特定点表现得像的程度,即对局部位置“流出”的度量——从一个无限小的空间区域流出的场向量比流入的场向量多出的程度。当一个点的散度为正时,该点为;当散度为负时,该点为。通过包围给定点的小表面的场通量越大,该点的散度绝对值就越大。通过封闭表面的通量为零的点具有零散度。

2. 散度的物理意义

例如,用气体的速度场来解释向量场的散度。气体的速度场是一个向量场。如果气体被加热,它会膨胀,这将导致气体分子向四周净运动。对场中任一个封闭的表面,气体将通过该曲面向外净流动,从而所有位置的速度场的散度为正。反之,若气体冷却,它将收缩,所有位置的速度场的散度为负。若空气温度不变,气体不膨胀和收缩,虽然气体仍然会流动,但流入任何封闭表面的气体流量速率等于流出的流量速率,其净流量为零,所有位置的的散度为零(这时的场称为螺线向量场)。

2. 散度的物理意义

如果气体仅在局部位置被加热,或者在一个点引入一个提供额外气体流量的小管,则那里的气体分子将向四周流出,这将在整个流场中产生一个以此点为中心的向外速度场。任何包围该点的封闭表面都会有气体从里面净流出,所以在该点具有正散度。当流场为稳定流时,任何没有包围该点的封闭表面内部都将有恒定的流速,流入和流出该区域的体积相等,这些位置的散度为零。

3.4 旋度

1. 旋度的定义

设 $\boldsymbol{F}$ 是向量场函数,将 $\nabla$ 看作向量,根据向量的向量积的计算方法可得

$$\nabla \times \boldsymbol{F}=\left| \begin{matrix} \boldsymbol{\hat{e}}_1& \boldsymbol{\hat{e}}_2& \boldsymbol{\hat{e}}_3\\ \frac{\partial}{\partial x}& \frac{\partial}{\partial y}& \frac{\partial}{\partial z}\\ F_x& F_y& F_z\\ \end{matrix} \right|$$

$\nabla \times \boldsymbol{F}$ 在物理学中被称作向量场 $\boldsymbol{F}$ 的 旋度,即

$$\text{curl} \, \boldsymbol{F}=\text{rot} \boldsymbol{F}=\nabla \times \boldsymbol{F}$$

1. 旋度的定义

向量场旋度的展开形式为

$$\begin{aligned} \nabla \times \boldsymbol{F}&=\left( \frac{\partial F_z}{\partial y}-\frac{\partial F_y}{\partial z} \right) \boldsymbol{\hat{e}}_1+\left( \frac{\partial F_x}{\partial z}-\frac{\partial F_z}{\partial x} \right) \boldsymbol{\hat{e}}_2+\left( \frac{\partial F_y}{\partial x}-\frac{\partial F_x}{\partial y} \right) \boldsymbol{\hat{e}}_3\\ &=\left[ \begin{array}{c} \frac{\partial F_z}{\partial y}-\frac{\partial F_y}{\partial z}\\ \frac{\partial F_x}{\partial z}-\frac{\partial F_z}{\partial x}\\ \frac{\partial F_y}{\partial x}-\frac{\partial F_x}{\partial y}\\ \end{array} \right]\\ \end{aligned}$$

旋度的定义

如右图所示,在向量场 $\boldsymbol{F}\left( \boldsymbol{x} \right)$ 中的任意点 $\boldsymbol{x}$ 的周围有一闭合有向曲线 $C$(其正方向规定为使得闭合曲线 $C$ 所包围的面积在它的左侧),$\boldsymbol{F}$ 沿着该曲线正方向的闭合曲线积分 $\oint_C{\boldsymbol{F}}\cdot \text{d}\boldsymbol{r}$ 称为向量场沿着该曲线的环量(或称旋涡量)。当闭合曲线所包围的面积趋于零时,环量与该面积微元的比值的极限值 $\lim_{A\rightarrow 0} \left( \frac{1}{|A|}\oint_C{\boldsymbol{F}}\cdot \text{d}\boldsymbol{r} \right)$ 称为向量场的环量面密度(或称环量强度)。

旋度的定义

对于不同方向的面元 $A$,所得环量面密度也不同。旋度是一个统一表征向量场中某一点旋转程度的向量,其所指方向对应的面元(与旋度方向垂直)的环量面密度最大(为旋度的模),它在一个方向上的投影的大小表示了在这个方向上的环量面密度的大小。即

$$\left( \nabla \times \boldsymbol{F} \right) \cdot \boldsymbol{\hat{n}}=\lim_{A\rightarrow 0} \left( \frac{1}{|A|}\oint_C{\boldsymbol{F}}\cdot \text{d}\boldsymbol{r} \right) $$

式中,$\boldsymbol{\hat{n}}$ 是面元 $A$ 的法向单位向量。

2. 旋度的物理意义

假设在流动的流体中放置一个表面粗糙的小球,球的中心固定在某一点(不可移动,可以绕中心向任意方向旋转),则流体的流动会使小球旋转。小球的旋转轴(根据右手法则确定)与小球中心位置的速度场旋度同向,旋转角速度是该点旋度大小的一半。

如果一个向量场中处处的旋度都是零,则称这个场为无旋场保守场

3.5 向量场的二阶微商

1. 梯度的旋度

针对如下二阶微商运算

$$\text{curl}\left( \text{grad}\,f \right) =\nabla \times \left( \nabla f \right) =\left( \nabla \times \nabla \right) f=0$$

从上式可得出两条结论:

  • 一个梯度的旋度为零。
  • 相反,如果 $\text{curl}\boldsymbol{F}$ 等于零,则 $\boldsymbol{F}$ 总是某标量场的梯度。即存在一个标量场 $\varPhi$(常称作标量势),使得 $\boldsymbol{F}=\nabla \varPhi$。

2. 旋度的散度

针对如下二阶微商运算(因为 $\nabla \times \boldsymbol{F}$ 垂直于 $\nabla$,所以 $\nabla$ 方向上没有 $\nabla \times \boldsymbol{F}$ 的分量)

$$\text{div}\left( \text{curl}\,\boldsymbol{F} \right) =\nabla \cdot \left( \nabla \times \boldsymbol{F} \right) =0$$

从上式可得出两条结论:

  • 一个旋度的散度为零。
  • 相反,如果 $\text{div}\boldsymbol{F}$ 等于零,则 $\boldsymbol{F}$ 总是某向量场的旋度。即存在一个向量场 $\boldsymbol{C}$,使得 $\boldsymbol{F}=\nabla \times \boldsymbol{C}$。

3. 拉普拉斯算符(作用于标量场)

针对如下二阶微商运算(因为 $\nabla \times \boldsymbol{F}$ 垂直于 $\nabla$,所以 $\nabla$ 方向上没有 $\nabla \times \boldsymbol{F}$ 的分量)

$$\begin{aligned} \nabla \cdot \left( \nabla f \right) &=\nabla _{x_1}\left( \nabla _{x_1}f \right) +\nabla _{x_2}\left( \nabla _{x_2}f \right) +\nabla _{x_3}\left( \nabla _{x_3} \right)\\ &=\frac{\partial ^2f}{\partial x_{1}^{2}}+\frac{\partial ^2f}{\partial x_{2}^{2}}+\frac{\partial ^2f}{\partial x_{3}^{2}}\\ \end{aligned}$$

上式可简写为

$$\nabla \cdot \left( \nabla f \right) =\nabla \cdot \nabla f=\left( \nabla \cdot \nabla \right) f=\nabla ^2f$$

3. 拉普拉斯算符(作用于标量场)

这里我们把 $\nabla ^2$ 看作一个新的算符。这是一个标量算符。由于它经常出现在物理学中,因而被赋予一个专用名称:拉普拉斯算符(或称拉普拉斯算子)。

$$\nabla ^2=\Delta =\frac{\partial ^2}{\partial x_{1}^{2}}+\frac{\partial ^2}{\partial x_{2}^{2}}+\frac{\partial ^2}{\partial x_{3}^{2}}$$

4. 拉普拉斯算符(作用于向量场)

虽然拉普拉斯算子是一个标量运算符,但同样可以对向量进行运算,只是要对每一个分量进行运算:

$$\nabla \left( \nabla \cdot \boldsymbol{F} \right) =\nabla ^2\boldsymbol{F}=\left( \nabla ^2F_{x_1},\nabla ^2F_{x_2},\nabla ^2F_{x_3} \right)$$

5. 旋度的旋度

对于下式,使用前面学到的向量三重积运算,可得

$$\begin{aligned} \nabla \times \left( \nabla \times \boldsymbol{F} \right) &=\nabla \left( \nabla \cdot \boldsymbol{F} \right) -\left( \nabla \cdot \nabla \right) \boldsymbol{F}\\ &=\nabla \left( \nabla \cdot \boldsymbol{F} \right) -\nabla ^2\boldsymbol{F}\\ \end{aligned}$$

3.6 向量的微分运算小结

向量的微分运算
运算 记法 描述 近似概念 范围
梯度 $$\text{grad}\left( f \right) =\nabla f$$ 表征标量场速率或方向的改变 标量乘法 将标量场映射为向量场
散度 $$\text{div}\left( \boldsymbol{F} \right) =\nabla \cdot \boldsymbol{F}$$ 表征向量场中给定点的是源还是汇 标量积 将向量场映射为标量场
旋度 $$\text{curl}\left( \boldsymbol{F} \right) =\nabla \times \boldsymbol{F}$$ 表征向量场中围绕一点旋转的趋势 向量积 将向量场映射为(伪)向量场

3.6 向量的微分运算小结

许多物理问题都可以用向量场的微分运算形式表示,如麦克斯韦方程组就可以表述为

$$\nabla \cdot \boldsymbol{E}=\frac{\rho}{\varepsilon _0} \\ \nabla \cdot \boldsymbol{B}=0 \\ \nabla \times \boldsymbol{E}=-\frac{\partial \boldsymbol{B}}{\partial t} \\ \nabla \times \boldsymbol{B}=\mu _0\left( \boldsymbol{J}+\varepsilon _0\frac{\partial \boldsymbol{E}}{\partial t} \right) $$

可见,采用向量微分运算符 $\nabla$ 后,许多物理问题都可以用优美而又简洁的公式描述。

第4节 向量的积分运算

上一节得出的三种基本的微分运算符都有对应的积分定理,从而使微积分基本定理推广到更高维度。

  • 梯度定理
  • 高斯散度定理
  • 旋度定理开尔文–斯托克斯定理)。

在二维空间内,散度和旋度定理变为格林公式

4.1 梯度定理

梯度定理也称线积分基本定理,可表述为:

标量场梯度沿曲线的积分可用标量场在该曲线两端的值之差来计算。
标量场梯度演曲线积分

设有一个标量场函数 $\varphi \left( \boldsymbol{r} \right)$(在三维直角坐标系中相当于 $\varphi \left( x_1,x_2,x_3 \right)$),则其梯度向量场函数为 $\nabla \varphi$。对于连接连接空间中两点 $\boldsymbol{p}$ 和 $\boldsymbol{q}$ 的任意曲线 $\varGamma$,沿此曲线从点 $\boldsymbol{p}$ 到 $\boldsymbol{q}$ 的积分符合如下关系:

4.1 梯度定理

$$\int_{\varGamma \left( \boldsymbol{p},\boldsymbol{q} \right)}{\nabla \varphi}\cdot \text{d}\boldsymbol{r}=\varphi \left( \boldsymbol{q} \right) -\varphi \left( \boldsymbol{p} \right)$$

梯度定理把微积分基本定理从直线数轴推广到平面、空间,乃至一般的 n 维空间中的曲线

梯度定理表明梯度场的曲线积分是路径无关的,这是物理学中保守力的定义方式之一。如果 $\varphi$ 是位势,则 $\nabla\varphi$ 就是保守向量场。上面的公式表明:保守力做功只和物体运动路径的端点有关,而与路径本身无关。

4.1 梯度定理

梯度定理的逆定理

如果 $\boldsymbol{F}$ 是个路径无关的向量场,则它是某个标量函数 $\varphi$ 的梯度。

当且仅当向量场沿任意闭曲线的积分为零时,该向量场是路径无关的。因此该逆定理也可表述如下:

如果向量场 $\boldsymbol{F}$ 沿定义域中每条闭曲线的积分为零,那么 $\boldsymbol{F}$ 是某个标量值函数 $\varphi$ 的梯度。

4.2 散度定理

散度定理又称高斯定理奥斯特罗格拉德斯基定理高斯-奥斯特罗格拉德斯基定理,这是一个把向量场通过闭合曲面的流动(即通量)与曲面内部的向量场的变化联系起来的定理,可表述为:

向量场穿过曲面的通量,等于散度在曲面围起来的体积上的积分。

或者用更直观的表述:

所有源点的和减去所有汇点的和,就是流出这区域的净流量。

4.2 散度定理

标量场梯度演曲线积分

假设空间中有一闭区域 V,该区域被分片光滑的闭合曲面 S 包围,$\boldsymbol{F}$ 是定义在 V 的邻域上的连续可微的向量场,则有

$$\iiint_V{\left( \nabla \cdot \boldsymbol{F} \right)}\,\mathrm{d}v=\oiint_S{\left( \boldsymbol{F}\cdot \boldsymbol{n} \right) \,\mathrm{d}s}$$

其中 $\boldsymbol{n}$ 是曲面 S 上指向区域外部的单位法向量。

4.2 散度定理

可以用液体的的流动来解释散度定理。运动的液体在每一点的速度构成一个向量场。假想流体内部有一个闭合曲面 S,在该闭合曲面内部形成了一个体积区域。液体流出该体积的流量等于穿过该表面的流体的体积速率,即表面上速度的表面积分。

  • 由于流体是不可压缩的,封闭区域内的液体量是恒定的。如果体积内没有源和汇,那么流出 S 的液体净流量必为零。
  • 如果体积内有源,如具有通向该体积的进水管道,则在源处额外出现的液体将对周围液体施加压力,导致液体在所有方向上向外流动。这将导致通过表面 S 的净向外流。通过表面 S 的向外流等于流体从管道流入 S 内的体积流量。

4.2 散度定理

  • 同样,如果体积内有汇,相当于有通向该体积的排水管,液体通过表面 S 向内流动的体积速率等于液体被排水管排出的速率。
  • 如果 S 内部有多个源和汇,液体流过源或汇的体积率(流过汇的流量给定一个负号)等于管口处速度场的散度,所以液体在由 S 包围的整个体积内的散度之和(积分)等于流过 S 的体积率。

散度定理适用于任何守恒定律:所有汇和源的总体积,即散度的体积积分,等于穿过体积边界的净流量。

4.2 散度定理

在三维条件下,散度定理可用高斯公式表述:

设空间闭区域 V 是由分片光滑的闭曲面 S 所围成,函数 $P\left( x,y,z \right)$、$Q\left( x,y,z \right)$ 和 $R\left( x,y,z \right)$ 在 V 上具有一阶连续偏导数,则有

$$\iiint_V{\left( \frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z} \right)}\mathrm{d}v=\oiint_S{\left( P\mathrm{d}y\mathrm{d}z+Q\mathrm{d}z\mathrm{d}x+R\mathrm{d}x\mathrm{d}y \right)}$$

$$\iiint_V{\left( \frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z} \right)}\mathrm{d}v=\oiint_S{\left( P\cos \alpha +Q\cos \beta +R\cos \gamma \right) \mathrm{d}S}$$

式中,$\cos \alpha$、$\cos \beta$ 和 $\cos \gamma$ 是 $S$ 在点 $\left( x,y,z \right)$ 处的法向量的方向余弦。

4.3 旋度定理

旋度定理又称开尔文-斯托克斯定理斯托克斯定理旋度基本定理,这是一个把向量场的旋度在曲面上的积分与向量场围绕该曲面的边界线的积分(即环流)联系起来的定理,可表述为:

向量场的旋度在一个闭曲面上的积分,等于向量场沿该曲面闭合的边界线的线积分。

向量场中闭曲面的旋度的通量,等于向量场围绕该曲面闭合的边界线的环流。

4.3 旋度定理

旋度定理示意图

假设空间中有一分段光滑的有向闭曲线 ΓS 是以 Γ 为边界的分片光滑的有向曲面,$\boldsymbol{F}$ 是定义在 S 的邻域上的连续可微的向量场,则有

$$\iint_S{\nabla \times \boldsymbol{F}\cdot \text{d}\boldsymbol{s}}=\oint_{\varGamma}{\boldsymbol{F}}\cdot \text{d}\boldsymbol{r}$$

4.3 旋度定理

在三维条件下,旋度定理可用斯托克斯公式表述:

Γ 为分段光滑的空间有向闭曲线,S 是以 Γ 为边界的分片光滑的有向曲面,Γ 的正方向与 S 的侧符合右手法则,函数 $P\left( x,y,z \right)$、$Q\left( x,y,z \right)$ 和 $R\left( x,y,z \right)$ 在 S (联通边界 Γ)上具有一阶连续偏导数,则有

$$\iint_S{\left( \frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z} \right) \text{d}y\text{d}z+\left( \frac{\partial P}{\partial z}-\frac{\partial R}{\partial x} \right) \text{d}z\text{d}x+\left( \frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y} \right) \text{d}x\text{d}y} \\ =\oint_{\varGamma}{\left( P\text{d}x+Q\text{d}y+R\text{d}z \right)}$$

4.4 格林公式

平面复连通区域

格林公式是旋度定律在二维空间中的表述。格林公式可表述为:

设闭区域 D 由分段光滑的曲线 L 围成,函数 P(x, y) 及 Q(x, y) 在 D 上具有一阶连续偏导数,则有

$$\iint_D{\left( \frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y} \right) \mathrm{d}x\mathrm{d}y}=\int_L{\left( P\,\mathrm{d}x+Q\,\mathrm{d}y \right)}$$

其中 LD 取正方向的边界线。

4.5 向量的积分运算小结

向量的积分运算
定理 表述 描述
梯度定理 $$\int_{\varGamma \left( \boldsymbol{p},\boldsymbol{q} \right)}{\nabla \varphi}\cdot \text{d}\boldsymbol{r}=\varphi \left( \boldsymbol{q} \right) -\varphi \left( \boldsymbol{p} \right)$$ 曲线 $L$ 上标量场梯度的曲线积分等于
曲线端点 $p$ 和 $q$ 之间标量场的变化。
散度定理 $$\iiint_V{\left( \nabla \cdot \boldsymbol{F} \right)}\,\mathrm{d}v=\oiint_S{\left( \boldsymbol{F}\cdot \boldsymbol{n} \right) \,\mathrm{d}s}$$ 一个向量场在一个 n 维固体上的散度积分等于
向量场通过该固体的 (n-1) 维封闭边界表面的通量。
旋度定理 $$\iint_S{\nabla \times \boldsymbol{F}\cdot \text{d}\boldsymbol{s}}=\oint_{\varGamma}{\boldsymbol{F}}\cdot \text{d}\boldsymbol{r}$$ 曲面 S 上的向量场的旋度积分等于包围
该曲面的封闭曲线 Γ 的向量场的环量。

参考文献

  谢谢!

返回目录
返回首页