(1) 物理问题#

一固定长度的柔软、均匀细弦，两端拉紧（一般两端固定）。在垂直于弦线方向施加外加激励（可能在初始时刻施加，也可能是在运动过程中持续施加）让它离开原来平衡位置， 在此外加激励作用下弦将作微小的横向振动，需要找到弦上个每个点的运动方程。

(2) 问题简化#

题目中提到的几处字眼：

• 柔软：只能承受轴向的载荷（弦只能承受张力，而杆能承受压力和拉力），忽略弯曲刚度（或者说弯矩）和剪切的影响，轴向载荷的方向与弦变形后的切线方向相一致。
• 轻质：不考虑重力所产生的影响。
• 均匀：单位长度的质量密度为一常数。
• 微小的横向振动：不考虑在张力作用下弦(（绳）的伸长变形（沿轴向的变形），只考虑弦（绳）上每个点的横向运动。弦（绳）的横向位移相对于弦（绳）的长度来说，是个小量。图中弦变形后的角度$\alpha$非常小，$\alpha\ll1$。

下面是问题简化后的研究示意图：

图示比较夸张，以便于画图，实际角度很小。

(3) 分析的参数与变量#

用$x$表示弦（绳）的位置，$u(x,t)$表示$t$时刻$x$位置处的横向位移。由于是微小的横向振动，故每个点只发生横向移动，忽略水平方向运动的影响。$\rho$表示单位长度的质量，$\mathrm{d}s$表示所分析微段的长度。

(4) 问题求解#

微段水平方向（$x$轴）的平衡：

微段垂直方向（$u$轴）的平衡：

对于水平方向上的平衡方程，根据小变形假设：横向位移相对于弦（绳）的长度来说是小量，图中弦变形后的角度$\alpha$非常小。

从而：

不考虑轴向运动影响,那么弦内张力保持不变。

对于垂直方向上的平衡方程，根据小变形假设：横向位移相对于弦（绳）的长度来说是小量。

由于$T_1=T_2=T$，则：

忽略高阶小量：

即：

设$a^2=\frac{T}{\rho}$，则：

由此，得到了弦（绳）的自由振动方程。$a$的量纲为速度量纲，可以将$a$称为波速，即波传播的速度。

如果考虑有单位长度$F(x,t)$的外加横向力作用，即$f(x,t)=F(x,t)/\rho$，可以得到：

化简得出强迫振动方程：

(1) 物理问题#

一张均匀张紧的轻质薄膜，离开水平位置沿垂直于水平位置的微小振动。

(2) 分析的参数与变量#

用$(x,y)$表示膜的位置，$u(x,y,t)$表示在$t$时刻在$(x,y)$位置处膜的横向位移。和弦问题一样， 根据水平方向的力平衡（任意方向），可知膜上的张力处处相同，记单位长度膜上的张力（注意是法线方向）为$T$。

(3) 问题求解#

取平面上微小单元$\mathrm{d}x\mathrm{d}y$进行受力分析，沿$x$轴方向两对边上由张力引起的横向合力为：（泰勒展开并忽略高阶小量）

沿$y$轴方向两对边上由张力引起的横向合力为：（泰勒展开并忽略高阶小量）

使用$\rho$表示单位面积的质量，对微元体受力分析，应用牛顿第二定律，可以得到：

可以将上式改写为：

NOTE

其中$\Delta_2$为二维拉普拉斯算符，也可以直接记作$\Delta$：

$$\Delta_2=\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}$$

(1) 物理问题#

一根均匀的细长杆，寻找沿轴向的纵向位移所应该遵循的方程。

(2) 分析的参数与变量#

$x$表示某点的位置，$u(x,t)$表示$t$时刻$x$位置处的纵向位移。

(3) 问题简化#

横截面变形后还是横截面，横截面形状保持不变。用位移的梯度$\varepsilon=\frac{\mathrm{d}u}{\mathrm{d}x}$（应变）来表征变形的程度。用单位横截面上的平均载荷来表征受力情况。

(4) 问题求解#

$B$段左边的位移为$u(x,t)$，$B$段右边的位移为$u(x+\mathrm{d}x,t)$，所以$B$段的变形量为：

根据胡克定律：

其中$S$为杆的横截面积，$E$为材料的弹性模量（类同于弹簧的劲度系数），$N$为杆拉压中的内力。

$B$ 段左边由于变形受到$A$ 段对它的作用力为：$ESu_x(x,t)$；$B$段右边由于变形受到$C$段对它的作用力为：$ESu_x(x+\mathrm{d}x,t)$。这里$S$为杆的横截面积，$E$为材料的弹性模量。

根据牛顿第二定律及力的平衡关系：

得到：

其中：

$a$具有速度量纲，也是波速。

(1) 物理问题#

由于物体内部物质浓度的分布不匀，导致物质从浓度高的区域向浓度低的区域转移（我们称这种现象为扩散）。分析物体扩散所需要满足的方程。

扩散问题遵循菲克第一定律：通过单位面积例子数率（通常称为粒子流通量）与粒子浓度的梯度成正比，方向与梯度方向相反，其比例系数称为扩散系数。并遵循质量守恒定律。

(2) 分析的参数与变量#

取空间中任意一个微元体（一般是一个平行六面体），根据扩散定律（菲克第一定律）和质量守恒定律，来分析浓度的变化过程。

用$(x,y,z)$表示位置，$t$表示时间，$u(x,y,z,t)$表示$(x,y,z)$位置在$t$时刻的浓度，那么，在某一个位置处，浓度的扩散流强度$q$为：

比例系数$D$即为扩散系数，与物质的性态、温度等有关。

(3) 问题求解#

单位时间内，沿$x$轴左边面上的流入量$q_1|_x=\mathrm{d}y\mathrm{d}z$，沿$x$轴右边面上的流出量$q_1|_{x+\mathrm{d}x}\mathrm{d}y\mathrm{d}z$，所以沿着$x$轴方向的净流入量为：

同理获得沿着其他轴方向的净流入量，获得总的流入量：

进一步化简得到：

根据质量守恒定律，在没有源（产生物质）或者汇（吸收物质）的条件下，该区域内物质总量的变化量等于流入的物质总量。即：

如果扩散系数在空间中是一致的，那么三维的扩散方程为：

其中$a^2=D$。

NOTE

其中$\Delta_3$为三维拉普拉斯算符：

$$\Delta_2=\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2}$$

退化到一维问题（也就是只有沿一个方向的扩散），那么扩散方程为：

对于有源或者有汇的扩散问题，那么相应的扩散方程为：

(1) 物理问题#

由于物体内部温度分布的不匀，导致热量从温度高的区域向温度低的区域转移（我们称这种现象为热传导）。分析热传导所需要满足的方程。

热传导问题遵循热传导定律（又称为傅里叶定律）：在导热过程中，单位时间内通过给定截面的导热量，正比于垂直于该截面方向上的温度变化率和截面面积；热量传递的方向与温度升高的方向相反（和扩散问题的菲克第一定律类似）。并遵循能量守恒定律。

(2) 分析的参数与变量#

与扩散问题的分析方法相同。取任意的一个微元体（一般为一个平行六面体），根据热传导定律来分析流经六面体上六个面单位时间的热能，再根据能（热）量守恒定律：通过边界进入微元体的中热量导致微元体内部的温度变化，从而得到相应的温度变化规律。

用$u(x,y,z,t)$表示$t$时刻在$(x,y,z)$处的温度，并利用热传导定律：

其中$q$成为热流强度，即单位时间通过单位横截面积的热量，反应了热传导能力的强弱，$k$称为热传导系数，或者说是热传导率。

(3) 问题求解#

仿照扩散问题的推导过程，应用热传导定律和能（热）量守恒定律，对于没有热源（汇）的一维和三维问题，相应的热传导方程分别为：

其中$c$为比热容，$\rho$为物体的质量密度。

化为标准形式为：

其中$a=\frac{k}{c\rho}$。

(1) 物理问题#

在电荷周围存在电场。电荷和电荷之间有力的作用，这个作用依靠电场来进行传递。由静电荷所产生的电场，我们称之为静电场。这里研究静电场满足的方程。

(2) 分析的参数与变量#

静电场遵循的物理规律是高斯定理：

表示的是闭合曲面$\sum$包裹的空间$T$中电荷量和通过闭合曲面电通量之间的关系。其中$\rho$表示单位体积的电荷量，$\varepsilon_0$为真空介电常数。

(3) 问题求解#

数学上学过高斯公式：

两个公式联立得到：

由于区域$T$是任意选择的，所以被积函数是相等的：

由于电场是无旋场（电力线不相交，电场中任意一点的电场强度不可能有两个方向），$\nabla\times T=\boldsymbol{0}$，因此存在一个势函数，使得：

结合被积函数的相等关系，得到：

最终得到静电场的势函数满足泊松方程：

那么在无源静电场的势函数满足拉普拉斯方程：

(1) 第一类（狄利克雷条件）#

第一类条件为给定了边界上物理量所要满足的条件。

1. 细杆导热问题：比如端点$x=a$的温度$u(x=a,t)$按已知函数$f(t)$发生变化。
2. 细杆振动问题：比如端点$x=a$固定，$u(x=a,t)=0$。
3. 弦的振动问题：比如端点$a=a,x=b$同时固定，$u(x=a,t)=u(x=b,t)=0$。

(2) 第二类（纽曼条件）#

第二类条件为给定物理量沿边界上外法线方向（$\vec{n}$）的导数（边界的梯度）所要满足的条件。

1. 细杆纵向振动问题：

   

   如果在右端部$x=b$有一个沿着轴向的外力$f(t)$该外力作用在边界上，那么在右边界取一个非常小的微段，根据受力，由于$\mathrm{d}m$很小，所以惯性力可以忽略，即：

   $$ES\frac{\partial u}{\partial x}\big|_{x=b}=f(t)=ES\frac{\partial u}{\partial \vec{n}}$$

   同样地，在左端侧进行同样的处理：

   $$-ES\frac{\partial u}{\partial x}=f(t)=ES\frac{\partial u}{\partial \vec{n}}$$

   这里外力$f(t)$指向$x$轴正方向。合起来：

   $$ES\frac{\partial u}{\partial \vec{n}}=f(t)$$

   如果边界自由：

   $$ES\frac{\partial u}{\partial \vec{n}}=0$$

2. 热传导问题：

   根据傅里叶定理，在任意位置的热流强度为：

   $$q=-k\nabla u$$

   所以沿着边界外法线方向的热流强度为：

   $$q_n=n\cdot q=-\vec{n}\cdot k\nabla u=-k\frac{\partial u}{\partial \vec{n}}$$

   边界上沿外法线的热流强度$q_n$大于零，说明热量是从物体内部向外部传递。反之亦然。如果在边界上给定了沿外法线方向的热流强度$f(t)$，那么根据傅里叶定律得到的热流强度在边界上应该等于给定的热流强度（$f(t)$只是 一个数值）：

   $$q_n=f(t)$$

第三类（混合边界条件，罗宾边界条件）#

第三类条件为给定边界上物理量与物理量沿外法线方向的方向导数之间所应该满足的关系。

1. 细长杆的纵向振动问题：

   在右边端点$x=b$有一个固定的弹簧支撑，右边端点附近去一个微段，忽略惯性力，根据胡克定律可以得到：

   $$\left(ES\frac{\partial u}{\partial \vec{n}}+Ku\right)\big|_{x=b}=0,\left(ES\frac{\partial u}{\partial x}+Ku\right)\big|_{x=b}=0$$

   在左边端点$x=a$有：

   $$\left(ES\frac{\partial u}{\partial \vec{n}}+Ku\right)\big|_{x=a}=0,\left(ES\frac{\partial u}{\partial x}-Ku\right)\big|_{x=a}=0$$

2. 细杆自由冷却问题：

   端部放在一个确定温度的环境中自由冷却，端部与环境之间由于存在温差有热对流，此时既不能立即确定温度，也不能立即确定温度的梯度。根据牛顿冷却定理：

   $$\frac{\mathrm{d}T}{\mathrm{d}t}=h(T_e-T)$$

   其中$h$为热交换系数，$T_e$为环境温度。

   根据热能守恒定律，边界上的热传导引起的热量变化量应该等于热交换引起的热量变化量。

   

   如果左端部$x=a$放置在温度$\theta$的环境中自由冷却，环境的热交换系数为$h$。

   $$q_n=-k\frac{\partial u}{\partial \vec{n}}|_{x=a}=h\left(u|_{x=a}-\theta\right)$$ $$-k\frac{\partial u}{\partial x}|_{x=a}=h\left(\theta-u|_{x=a}\right)$$ $$\frac{\partial u}{\partial x}|_{x=a}=-\frac{\partial u}{\partial \vec{n}}|_{x=a}$$

   如果右端部$x=b$放置在温度$\theta$的环境中自由冷却，环境的热交换系数为$h$。

   $$q_n=-k\frac{\partial u}{\partial \vec{n}}|_{x=b}=h\left(u|_{x=b}-\theta\right)$$ $$-k\frac{\partial u}{\partial x}|_{x=b}=h\left(u|_{x=b}-\theta\right)$$ $$\frac{\partial u}{\partial x}|_{x=b}=\frac{\partial u}{\partial \vec{n}}|_{x=b}$$

   对于一般表面的自由冷却问题：

   $$q_n=-k\frac{\partial u}{\partial \vec{n}}|_{x=a}=h(u-\theta)$$

(1) 相关定义与性质#

如果系数 $a_ij,b_i,c$及$f$与$u$无关，只是$x_1,x_2,...,x_n$的函数或者是常数，那么该方程对于$u$及其各阶导数来说是线性的，因此微分方程也是线性的。

记微分算子$L$：

$L$是线性，可以得到：

线性二阶偏微分方程中，如果自由项$f\equiv0$，称方程为齐次的，否则称为非齐次的。

(2) 定解条件#

对于齐次方程：

对于非齐次方程：

(3) 叠加原理#

若泛定方程和定解条件都是线性的（可以是非齐次），可以把泛定方程和定解条件分别拆成多个线性方程的合理叠加。只要各部分的解分别满足各自的泛定方程和定解条件，那么叠加后的函数可以满足原来的方程和定解条件。也就是说原来问题的解可以看成各部分解的线性叠加。利用叠加原理可以将复杂方程的解拆分成几个比较简单问题的解。

在物理上，可以理解为：几种不同原因综合所产生的效果，等于这些不同原因单独产生效果的累加。

(1) 双曲型方程#

特征方程有两组不同的解，用与之对应的两个$z$的解作为变换函数，可以将$A_{11}$与$A_{22}$变为$0$；根据上文推导，知：

所以$z(x, y)$ 为：

于是得到了最开始的变量代换方法：

再作自变量变换：

从而得到用自变量 $\alpha$ 与 $\beta$ 所表示的双曲型方程标准形式：

前面的各种波动（或振动）方程就属于双曲型方程。

(2) 抛物型方程#

由于$\Delta =0$，所以

对应的特征变量变换为：

只取一组特征变量，令：

两边分别对 $x$ 和 $y$ 求偏导，得到：

把 $\xi$ 代入 $A_{11}$、$A_{12}$的表达式：

由于 $\frac{\xi_x}{\xi_y} = -\frac{a_{12}}{a_{11}}$，故有：

把 $\xi$ 代入$A_{22}$ 的表达式

即：

只要选择合适的 $\eta$，使得

即可保证 $A_{22} \neq 0$。从而得到新坐标下的控制方程，即抛物方程的标准形式：

前面的扩散方程、热传导方程等都属于抛物方程。

(3) 椭圆方程#

特征方程只有一组共轭复函数解 $\xi(x, y)$ 及 $\xi^*(x, y)$。由于 $A_{11} = 0,\quad A_{22} = 0$，此时方程为：

我们取 $\alpha$ 和 $\beta$ 分别为 $\xi(x, y)$ 函数的实部和虚部，即：

那么，方程变为：

这是椭圆方程的标准形式。前面的稳定浓度分布、稳定扩散分布及静电场势函数等都属于椭圆方程。

(1) 考虑无限长弦的波动问题#

对于足够长的弦，在有限时间内，远处的边界条件不会对结果有影响。 对于足够长的弦，我们可以把它近似为一个无限长的弦问题（没有边界条件）。

初始条件分别为位移的初始条件和速度的初始条件：

把两个初始条件代入通解表达式：

分别可以得到：

对第二个方程的两边积分

联立第一个初始条件方程，得到：

带入通解表达式，得到达朗贝尔公式：

其中包含两项：第一项表示初始形状带来的影响，后一项积分项表示初始速度带来的影响。下面将讨论该问题的两种特殊情况

1. 如果无限长弦初始是静止的，即初始条件为：

   $$u|_{t=0}=\varphi(x),u_t|_{t=0}=0\quad (-\infty < x < \infty)$$

   根据达朗贝尔公式，得到相应的解为：

   $$u(x,t)=\frac{1}{2}[\varphi(x+at)+\varphi(x-at)]$$

   只保留下了达朗贝尔公式的第一项，即该波动问题的解可以看成是两只波的叠加，其中一支波形$\varphi(x)/2$以波速$a$向左移动，另一支波形$\varphi(x)/2$以波速$a$向右移动。

   

2. 如果无限长的弦初始位置为零，并具有初始速度，即初始条件为：

   $$u|_{t=0}=0, u_t|_{t=0}=\psi(x)\quad (-\infty < x < \infty)$$

   根据达朗贝尔公式，相应的解可以写为：

   $$u(x,t)=\frac{1}{2a}\int_{x-at}^{x+at}\psi(\xi)\mathrm{d}\xi$$

   利用定积分性质写为：

   $$u(x,t)=\frac{1}{2a}\int_{-\infty}^{x+at}\psi(\xi)\mathrm{d}\xi-u(x,t)=\frac{1}{2a}\int_{-\infty}^{x-at}\psi(\xi)\mathrm{d}\xi$$

   记：

   $$\Psi(x)=\frac{1}{2a}\int_{-\infty}^{x}\psi(\xi)\mathrm{d}\xi$$

   那么解为：

   $$u(x,t)=\Psi(x+at)-\Psi(x-at)$$

(2) 半无限长弦的振动#

该问题与无限长弦振动问题的区别在于方程成立的定义域，其中一头是有边界的：

初始条件为：

在左端常见的有三种边界条件：

1. 端点固定

   

   $$u(0,t)=0$$

2. 张力的竖直方向的分离为$0$

   $$ u_x(0,t)=0 $$

3. 弹性支撑

   $$ u_x(0,t)-hu(0,t)=0 $$

需要根据左边的边界条件来考虑进行奇延拓还是偶延拓，以左端固定的半无限长弦振动问题举例，满足：

把初始条件进行奇延拓（因为函数$x=0$处值为$0$，如果一阶导为$0$，则使用偶延拓）：

这样该问题就能等价成一个无限长弦的振动问题：

$x>0$部分的解$v(x,t)$就是原来问题的解$u(x,t)$。根据达朗贝尔公式，写出无限长弦振动问题的解为：

我们关心的事$x>0$部分的解，所以将原函数回代时，要进行分段函数的转化：

对于$x>at$时：

对于$x<at$时：