Freshmark
← 返回全部文章

微积分基础导论:微积分应用1

本文介绍微积分在物理中的应用,涵盖质心计算(圆锥、半球壳、半球体)、保守力势能与平衡稳定性判据、泰勒小量展开及其在电偶极子受力、相对论动能还原、非均匀浮力等经典问题中的应用。

朝花夕拾

alt text

例1

有一个质量均匀的圆锥,其母线与高的夹角为θ\theta,底面半径为rr,密度为ρ\rho.

(1)求总质量

(2)求重心位置

设原点为锥顶,圆锥的高为x轴

(1)

dm=ρSdh=πρ(htan⁡θ)2dhm=∫dm=πρtan⁡2θ∫0rcot⁡θh2dh=πρtan⁡2θ(13h3)∣0rcot⁡θ=13πρr3cot⁡θ\begin{gathered} dm=\rho Sdh\\ =\pi\rho (h\tan\theta)^2dh\\ m=\int dm\\ =\pi\rho\tan^2\theta\int_0^{r\cot\theta}h^2dh\\ =\pi\rho\tan^2\theta\left(\frac{1}{3}h^3)\right|_0^{r\cot\theta}\\ =\frac{1}{3}\pi\rho r^3\cot\theta \end{gathered}

(2)

由对称性知,重心在圆锥的高上.

xc=∫mixim=∫0rcot⁡θπρtan⁡2θx3dxm=πρtan⁡2θ∫0rcot⁡θx3dx13πρr3cot⁡θ=3tan⁡3θr3(rcot⁡θ)44=34rcot⁡θ\begin{gathered} x_c=\frac{\int m_ix_i}{m}\\ =\frac{\int_0^{r\cot\theta}\pi\rho\tan^2\theta x^3dx}{m}\\ =\frac{\pi\rho\tan^2\theta\int_0^{r\cot\theta}x^3dx}{\frac{1}{3}\pi\rho r^3\cot\theta}\\ =3\frac{\tan^3\theta}{r^3}\frac{(r\cot\theta)^4}{4}\\ =\frac{3}{4}r\cot\theta \end{gathered}

例3

求均匀半球壳的重心位置

以球心为原点,垂直截面方向为x轴

不妨设单位面积质量为1.

仿照求n维球的体积,我们考虑在θ,dθ\theta,d\theta处各切一刀。

与半球壳的本质区别

n维球体积半球壳面积
切片方式水平薄片,厚度 ∥dy∥=Rsin⁡θ dθ\|dy\| = R\sin\theta\,d\theta球面环形,弧长 dl=R dθdl = R\,d\theta
Rsin⁡θ dθR\sin\theta\,d\theta✅ 正确(竖直厚度)❌ 错误(应为弧长)

dm=(2πRsin⁡θ)Rdθxc=∫0π2Rcos⁡θdm2πR2=2πR3∫0π2sin⁡θcos⁡θdθ2πR2=2πR3∫0π2sin⁡θdsin⁡θ2πR2=2πR3∫01udu2πR2=2πR3×122πR2=12R\begin{gathered} dm=(2\pi R\sin\theta)Rd\theta\\ x_c=\frac{\int_{0}^{\frac{\pi}{2}} R\cos\theta dm}{2\pi R^2}\\ =\frac{2\pi R^3\int_{0}^{\frac{\pi}{2}}\sin\theta\cos\theta d\theta}{2\pi R^2}\\ =\frac{2\pi R^3\int_{0}^{\frac{\pi}{2}}\sin\theta d\sin\theta}{2\pi R^2}\\ =\frac{2\pi R^3\int_{0}^{1}udu}{2\pi R^2}\\ =\frac{2\pi R^3\times\frac{1}{2}}{2\pi R^2}\\ =\frac{1}{2}R \end{gathered}

例4

求均匀半球体的重心位置

以球心为原点,垂直截面方向为x轴

不妨设单位体积质量为1.

我们平行截面横切:

法一

dm=π(Rsin⁡θ)2sin⁡θRdθxc=∫Rcos⁡θdm23πR3=πR4∫0π2sin⁡3θcos⁡θdθ23πR3=32R∫0π2sin⁡3θdsin⁡θ=32R∫01u3du=38R\begin{gathered} dm=\pi (R\sin\theta)^2\sin\theta Rd\theta\\ x_c=\frac{\int R\cos\theta dm}{\frac{2}{3}\pi R^3}\\ =\frac{\pi R^4\int_{0}^{\frac{\pi}{2}}\sin^3\theta\cos\theta d\theta}{\frac{2}{3}\pi R^3}\\ =\frac{3}{2}R\int_{0}^{\frac{\pi}{2}}\sin^3\theta d\sin\theta\\ =\frac{3}{2}R\int_{0}^1u^3du\\ =\frac{3}{8}R \end{gathered}

法二

dm=π(R2−x2)dxxc=∫dmx23πR3=π∫0R(R2−x2)xdx23πR3=π(R22x2−14x4)∣0R23πR3=38R\begin{gathered} dm=\pi (R^2-x^2)dx\\ x_c=\frac{\int dmx}{\frac{2}{3}\pi R^3}\\ =\frac{\pi \int_{0}^R (R^2-x^2)xdx}{\frac{2}{3}\pi R^3}\\ =\frac{\pi \left(\frac{R^2}{2}x^2-\frac{1}{4}x^4\right)|_0^R}{\frac{2}{3}\pi R^3}\\ =\frac{3}{8}R \end{gathered}

平衡的稳定性

在保守力下,若势能处于极值时,则平衡:

  • 势能处在极大值:不稳定平衡
  • 势能处在极小值:稳定平衡

说文解字:保守力做功与路径无关,只与初末位置相关

重力,引力,静电力都是保守力

摩擦力是典型的非保守力

alt text

导数极值第二充分条件

定理内容

设函数 f(x)f(x) 在点 x0x_0 处具有二阶导数,且 f′(x0)=0f'(x_0) = 0(即 x0x_0 为驻点):

  1. 若 f′′(x0)<0f''(x_0) < 0,则 f(x)f(x) 在 x0x_0 处取得极大值。
  2. 若 f′′(x0)>0f''(x_0) > 0,则 f(x)f(x) 在 x0x_0 处取得极小值。
  3. 若 f′′(x0)=0f''(x_0) = 0,则该定理失效,无法判定是否取得极值(需使用第一充分条件或更高阶导数判断)。

应用示例

求函数 f(x)=x3−3xf(x) = x^3 - 3x 的极值。

解: 首先求一阶导数: f′(x)=3x2−3f'(x) = 3x^2 - 3

令 f′(x)=0f'(x) = 0,解得驻点为: x1=1,x2=−1x_1 = 1, \quad x_2 = -1

再求二阶导数: f′′(x)=6xf''(x) = 6x

将驻点代入二阶导数进行判断:

  • 对于 x1=1x_1 = 1: f′′(1)=6>0f''(1) = 6 > 0 所以 f(x)f(x) 在 x=1x = 1 处取得极小值,极小值为 f(1)=−2f(1) = -2。
  • 对于 x2=−1x_2 = -1: f′′(−1)=−6<0f''(-1) = -6 < 0 所以 f(x)f(x) 在 x=−1x = -1 处取得极大值,极大值为 f(−1)=2f(-1) = 2。

例5

有一个质量为mm,劲度系数为kk,原长为2πr02\pi r_0的橡皮筋,位于母线与高线夹角为θ\theta的圆锥上.

(1)求平衡位置

(2)求上下扰动稳定性

alt text

(1)

设平衡位置到锥顶的高度差为zz,平衡时橡皮筋伸长量为xx

x=2π(ztan⁡θ−r0)Ep=12kx2−mgz=2π2k(ztan⁡θ−r0)2−mgz=2π2ktan⁡2θz2−(mg+4π2kr0tan⁡θ)z+2π2kr02z0=mg+4π2kr0tan⁡θ4π2ktan⁡2θ\begin{gathered} x=2\pi (z\tan\theta-r_0)\\ E_p=\frac{1}{2}kx^2-mgz\\ =2\pi^2 k(z\tan\theta-r_0)^2-mgz\\ =2\pi^2k\tan^2\theta z^2-(mg+4\pi^2kr_0\tan\theta)z+2\pi^2kr_0^2\\ z_0=\frac{mg+4\pi^2kr_0\tan\theta}{4\pi^2k\tan^2\theta} \end{gathered}

(2)

Ep(z)E_p(z)在z=z0z=z_0处取得极小值,故平衡为稳定平衡.

例6

光滑圆环半径为RR,弹簧一端挂在圆环正上方,另一端拴着质量为mm的小球.圆环穿过小球.弹簧原长为0,劲度系数为kk

(1)求平衡位置

(2)求平衡的稳定性.

alt text

(1)

设弹簧与竖直方向夹角为θ\theta

Ep=12k(2Rcos⁡θ)2+mg(2Rsin⁡2θ)=2kR2cos⁡2θ+2mgR(1−cos⁡2θ)=2R(kR−mg)cos⁡2θ+2mgRdEpdθ=4R(mg−kR)cos⁡θsin⁡θ=0→sin⁡2θ=0θ=0 or π2\begin{gathered} E_p=\frac{1}{2}k(2R\cos\theta)^2+mg(2R\sin^2\theta)\\ =2kR^2\cos^2\theta+2mgR(1-\cos^2\theta)\\ =2R(kR-mg)\cos^2\theta+2mgR\\ \frac{dE_p}{d\theta}=4R(mg-kR)\cos\theta\sin\theta=0\\ \rightarrow \sin2\theta=0\\ \theta=0 \text{ or } \frac{\pi}{2} \end{gathered}

平衡位置为A点,B点

(2)

d2Epdθ2=4R(mg−kR)cos⁡2θ\begin{gathered} \frac{d^2E_p}{d\theta^2}=4R(mg-kR)\cos2\theta \end{gathered}

若mg<kRmg\lt kR,则θ=0\theta=0时,d2Epdθ2<0\frac{d^2E_p}{d\theta^2}\lt0

为势能极大值点,B点属于不稳定平衡.

θ=π2\theta=\frac{\pi}{2}时,d2Epdθ2>0\frac{d^2E_p}{d\theta^2}\gt0

为势能极小值点,A点属于稳定平衡

mg>kRmg\gt kR,情况恰相反

小量展开

例7

f(x)=xf(x)=\sqrt{x},估算f(4.01)f(4.01)

泰勒展开

鸭子原则:有一个长得像鸭子的东西,听声音也很像鸭子,吃肉也是鸭子味,那这就是一只鸭子

不断地获得更多的信息,即可逼近真相

f(x0+x)=f(x0)+f′(x0)x+f′′(x)2!x2+f′′′(x)3!x3+...\boxed{f(x_0+x)=f(x_0)+f'(x_0)x+\frac{f''(x)}{2!}x^2+\frac{f'''(x)}{3!}x^3+...}

例8

f(x)=xαf(x)=x^\alpha在1附近的展开:

f′(x)=αxα−1,f′(1)=αf′′(x)=α(α−1)xα−2,f′′(1)=α(α−1)...f(1+x)=1+αx+α(α−1)2x2+...其中,∣x∣<1\begin{gathered} f'(x)=\alpha x^{\alpha-1},f'(1)=\alpha\\ f''(x)=\alpha(\alpha-1)x^{\alpha-2},f''(1)=\alpha(\alpha-1)\\ ...\\ f(1+x)=1+\alpha x+\frac{\alpha(\alpha-1)}{2}x^2+...\\ \text{其中,}|x|\lt 1 \end{gathered}

同理得到其他常见泰勒展开

sin⁡θ=θ−13!θ3+15!θ5+...cos⁡θ=1−12!θ2+14!θ4ex=1+x+12!x2+13!x3+...ln⁡(1+x)=x−12x2+13x3−...\begin{gathered} \sin\theta=\theta-\frac{1}{3!}\theta^3+\frac{1}{5!}\theta^5+...\\ \cos\theta=1-\frac{1}{2!}\theta^2+\frac{1}{4!}\theta^4\\ e^x=1+x+\frac{1}{2!}x^2+\frac{1}{3!}x^3+...\\ \ln (1+x)=x-\frac{1}{2}x^2+\frac{1}{3}x^3-... \end{gathered}

例9

证明:1+eiπ=01+e^{i\pi}=0或eiθ=cos⁡θ+isin⁡θe^{i\theta}=\cos\theta+i\sin\theta,其中ii为虚数单位.

把iθi\theta带入exe^x泰勒展开易证.

例10

例题:电偶极子的受力

已知条件

库仑定律:

F=kQqr2F = \frac{kQq}{r^2}

  • 若 QQ、qq 同号,则 F>0F > 0,为斥力
  • 若 QQ、qq 异号,则 F<0F < 0,为引力

题目

电偶极子由一对等量异种电荷 +q+q 与 −q-q 组成,两者间距为 ll。

求:

(1) 电偶极子轴线延长线上,距偶极子中心距离为 rr 处的点电荷 QQ 所受的力。

条件:r≫lr \gg l,保留到最低阶小量。

(2) 电偶极子轴线与连线成 θ\theta 角、距偶极子中心距离为 rr 处的点电荷 QQ 所受的力(分解为径向分量 FrF_r 和切向分量 FθF_\theta)。

条件:r≫lr \gg l,保留到最低阶小量。

alt text

(1)

F=kQq(r−l2)2+kQ(−q)(r+l2)2=kQqr2[(1−l2r)−2−(1+l2r)−2]≈kQqr22lr=2kQqlr3\begin{gathered} F=\frac{kQq}{(r-\frac{l}{2})^2}+\frac{kQ(-q)}{(r+\frac{l}{2})^2}\\ =\frac{kQq}{r^2}[(1-\frac{l}{2r})^{-2}-(1+\frac{l}{2r})^{-2}]\\ \approx \frac{kQq}{r^2}{\frac{2l}{r}}\\ =\frac{2kQql}{r^3} \end{gathered}

其中p=qlp=ql,称为偶极矩

(2)

alt text

Fr=F1cos⁡Δθ1−F2cos⁡Δθ2=F1(1−12Δθ12)−F2(1−12Δθ22)=忽略二阶小量F1−F2=kQqr2+(l2)2−2rl2cos⁡θ−kQqr2+(l2)2+2rl2cos⁡θ=kQqr2(11−lrcos⁡θ+l24r2−11+lrcos⁡θ+l24r2)=泰勒展开kQqr2(2lcos⁡θr)=2kQqlcos⁡θr3\begin{gathered} F_r=F_1\cos\Delta\theta_1-F_2\cos\Delta\theta_2\\ =F_1(1-\frac{1}{2}\Delta\theta_1^2)-F_2(1-\frac{1}{2}\Delta\theta_2^2)\\ \xlongequal{\text{忽略二阶小量}}F_1-F_2\\ =\frac{kQq}{r^2+(\frac{l}{2})^2-2r\frac{l}{2}\cos\theta}-\frac{kQq}{r^2+(\frac{l}{2})^2+2r\frac{l}{2}\cos\theta}\\ =\frac{kQq}{r^2}(\frac{1}{1-\frac{l}{r}\cos\theta+\frac{l^2}{4r^2}}-\frac{1}{1+\frac{l}{r}\cos\theta+\frac{l^2}{4r^2}})\\ \xlongequal{\text{泰勒展开}}\frac{kQq}{r^2}(\frac{2l\cos\theta}{r})\\ =\frac{2kQql\cos\theta}{r^3} \end{gathered}

Fθ=kQqlsin⁡θr3\begin{gathered} F_\theta=\frac{kQql\sin\theta}{r^3} \end{gathered}

例11

在相对论中,物体的总能量E=mc2E=mc^2

其中mm为动质量:m=11−v2c2m0m=\frac{1}{\sqrt{1-\frac{v^2}{c^2}}}m_0

其中vv为物体速度,cc为光速,m0m_0为静质量

称E0=m0c2E_0=m_0c^2为静能量,有

Ek=E−E0\boxed{E_k=E-E_0}

而根据经典物理Ek=12mv2E_k=\frac{1}{2}mv^2

证明:在v<<cv\lt\lt c时,相对论中的EkE_k还原为经典物理学中的动能

Ek=11−v2c2m0c2−m0c2=m0c2[(1−v2c2)−12−1]=保留一阶小量12m0v2\begin{gathered} E_k=\frac{1}{\sqrt{1-\frac{v^2}{c^2}}}m_0c^2-m_0c^2\\ =m_0c^2[(1-\frac{v^2}{c^2})^{\frac{-1}{2}}-1]\\ \xlongequal{\text{保留一阶小量}}\frac{1}{2}m_0v^2 \end{gathered}

例12

有一杯糖水,其密度ρ(x)\rho(x)与深度相关,ρ(x)=ρ0exl0\rho(x)=\rho_0e^\frac{x}{l_0}

有一个密度为ρ1=1.5ρ0\rho_1=1.5\rho_0的均匀硬棒,截面积为SS,长度为l0l_0

求平衡时,顶端到水面的距离(在水面上还是水面下?). m0=ρ1Sl0=1.5ρ0Sl0m_0=\rho_1Sl_0=1.5\rho_0Sl_0 假设上端在水面上,浸入水中的长度x0<l0x_0\lt l_0

F=∫0x0ρ(x)gSdx=ρ0gS∫0x0exldx=ρ0gSl(ex0l−1)\begin{gathered} F=\int_{0}^{x_0}\rho(x)gSdx\\ =\rho_0gS\int_{0}^{x_0}e^\frac{x}{l}dx\\ =\rho_0gSl(e^{\frac{x_0}{l}}-1) \end{gathered}

要求F=m0g=ρ1Sl0=1.5ρ0gSl0F=m_0g=\rho_1Sl_0=1.5\rho_0gSl_0

即ex0l0=2.5,x0<l0e^\frac{x_0}{l_0}=2.5,x_0\lt l_0,假设成立

Sonnet 4.6总结

总结

本节围绕微积分在物理建模中的核心应用,串联了以下几条主线:

质心计算:通过微元积分,分别求解了均匀圆锥(34rcot⁡θ\frac{3}{4}r\cot\theta)、半球壳(12R\frac{1}{2}R)和半球体(38R\frac{3}{8}R)的重心位置,体会切片方式的选取对积分的影响。

平衡与稳定性:引入势能极值判据——势能极小对应稳定平衡,势能极大对应不稳定平衡。通过圆锥上的橡皮筋和圆环弹簧小球两个模型,将该判据具体化为对 EpE_p 二阶导数的符号分析。

泰勒展开与小量近似:掌握 (1+x)α(1+x)^\alpha、sin⁡θ\sin\theta、cos⁡θ\cos\theta、exe^x、ln⁡(1+x)\ln(1+x) 等常见展开式,并应用于:

  • 电偶极子受力(引出偶极矩 p=qlp = ql)
  • 相对论动能在 v≪cv \ll c 下还原为 12m0v2\frac{1}{2}m_0v^2
  • 非均匀密度糖水中硬棒的浮力平衡

核心思想:复杂问题的关键往往在于选取合适的微元与坐标系,再结合小量展开保留主要项,从而将繁琐的精确计算化为简洁的近似结果。

评论

评论

登录后参与讨论。注册时只需验证一次邮箱。

正在加载评论…

正在检查登录状态…