二维波动方程

                     

贡献者: 待更新

  • 本文处于草稿阶段。
预备知识 一维波动方程,梯度,散度

   假设我们有一个紧绷的橡皮薄膜,单位长度的张力为 $\lambda$。面密度为 $\sigma$,均为常数。取一小块微元进行受力分析,其质量为 $m = \sigma S$。要分析其受力,可以将面元的边界划分成许多小段,每小段给面元施加的竖直方向的力等于 $\lambda$ 乘以法相斜率 $ \boldsymbol\nabla u$ 再乘以该小段长度 $ \,\mathrm{d}{l} $。注意只有再假设其振动幅度很小时这才成立,因为斜率是 $\tan\theta$,而严格来说我们需要 $\sin\theta$,类比同一维波动方程的推导

\begin{equation} \frac{\partial^{2}{u}}{\partial{t}^{2}} = \frac{\lambda}{\sigma S}\oint [ \boldsymbol\nabla u( \boldsymbol{\mathbf{r}} )] \boldsymbol\cdot \hat{\boldsymbol{\mathbf{n}}} ( \boldsymbol{\mathbf{r}} ) \,\mathrm{d}{l} ~, \end{equation}
$ \hat{\boldsymbol{\mathbf{n}}} ( \boldsymbol{\mathbf{r}} )$ 为边界在 $ \boldsymbol{\mathbf{r}} $ 处向外的法向量。

   当 $S \to 0$,根据散度的定义

\begin{equation} \frac{1}{S} \oint [ \boldsymbol\nabla u( \boldsymbol{\mathbf{r}} )] \boldsymbol\cdot \hat{\boldsymbol{\mathbf{n}}} ( \boldsymbol{\mathbf{r}} ) \,\mathrm{d}{l} \to \boldsymbol{\nabla}\boldsymbol{\cdot} ( \boldsymbol\nabla u) = \boldsymbol{\nabla}^2 u~, \end{equation}
所以波动方程为
\begin{equation} \boldsymbol{\nabla}^2 u - \frac{\sigma}{\lambda} \frac{\partial^{2}{u}}{\partial{t}^{2}} = 0~, \end{equation}
其中 $v_0 = \sqrt{\lambda/\sigma}$。

   特殊地,当薄膜保持静止不动,我们就得到了拉普拉斯方程


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利