万有引力的推导作者:老司机摘要很多中学生以为已知开普勒三大定律就能推导出牛顿万有引力定律,其实并不是如此.仅仅依靠开普勒三定律是没有办法推导出牛顿万有引力方程的.为了解开这个误区,今天就让我们一起来探索牛顿万有引力方程是如何推导出来的吧.233关键词:万有引力定律,比耐公式,开普勒定律,理论力学,物理拔高,毁梗用的1基本物理概念首先我们要引入角动量L与力矩M这个概念.定义角动量的表达式:L=r×p定义力矩的表达式:M=r×F力矩和角动量之间有如下关系:L=dM dt证明:我们已知牛顿第二定律F=ma=m dvdt,两边同时对牛顿第二定律叉乘r.得:r×F=mr×dv dt我们知道乘积求导有这样的一个关系:(ab) =a b+ab ,d(ab)dt=dadtb+adbdt 1把这个关系带入之前的式子:r×F=mr×dvdt=md(r×v)dt−mv×drdt我们知道v=drdt ,所以上式最后一项是mv×drdt=mv×v=0这样我们就得到了关系式L=r×F=d(r×mv)dt =dMdt角动量定理:L=dM dt我们现在知道了角动量和力矩的概念后我们就可以开始去探索如何推导出牛顿万有引力方程了.2比耐公式首先我们先看看在一般的中心力场中的规律.所谓中心力场,就是满足F=F(r) rr这样的力场,比如万有引力只与距离r有关,我们把这种只与距离有关而与其他无关(比如角度θ,φ)的力场叫做中心力场.中心力场中运动的物体一定是在一个平面内的轨迹.图1:中心力场(力心在o点处)根据牛顿第二定律:F(r) rr=ma=md2 rdt2写成x,y分量形式.m d2xdt2=F(r)xrm d2ydt2=F(r)yr2把直角坐标和极坐标互化公式x=rcosθ,y=rsinθ带进去.得:m(d2rdt2−r(dθdt)2)=F(r)(1)m(r d2θdt2+2drdtdθdt)=0(2)把(2)式凑全微分,得m1r dt(r2dθdt)=0,所以:r2dθdt=Constant(3) mr2dθdt=Constant(4)这样我们得到了中心力场的基本方程组m(d2rdt2−r(dθdt)2)=F(r)(5) r2dθdt=C(6)现在我们知道了在中心力场中的运动规律满足上面两个式子,那么我们对上面两个式子消去时间变量t就得到任意中心力场F(r)下的运动轨迹方程r= r(θ)(这是个极坐标方程).我们已知r2dθdt =C,做变量代换,以u=1r代换.得dθdt=Cu2dr dt =drdθdθdt=ddθ(1u)dθdt=−1u2dudθdθdt=−Cdudθd2r dt2=ddtdrdt=ddr(−Cdudθ)=ddθ(−Cdudθ)dθdt=−C2u2d2udθ2把以上三式带入(5)式,得:C2u2(d2udθ2+u)=−Fm这就是所要求的轨道微分方程,通常叫做比耐公式,引力时,F为负号,斥力时F为正号.由这个方程可知,若我们已知中心力场具体轨道形式,便可以求出该中心力场的力的形式.33万有引力定律好了,咱们回到推导万有引力定律来.1609年,开普勒发表了他的三大定律.开普勒第一定律,行星绕太阳做椭圆运动,太阳位于椭圆的一个焦点上.开普勒第二定律,行星和太阳之间的连线,在相等时间内扫过的面积相等.开普勒第三定律,行星公转的周期的平方和轨道半长轴的立方成正比.首先看开普勒第二定律.图2:开普勒第二定律设A 是矢径扫过的面积,由开普勒第二定律,知道单位面积内,矢径所扫过的面积相等,即dAdt=constant 现在来求dA dt 的表达式,P 1,P 2分别是行星沿着它轨道运动时的两个相邻位置,对太阳张开的角度为dθ,从P 1到P 2的时间是∆t ,在这一段时间内扫过的面积∆A 为OP 1P 2.当∆t →0的时候,P 2→P 1,∆A 就近似的等于∆OP 1P 2的面积,即12r (r ∆θ),所以.dA dt =lim ∆t →0∆A ∆t =lim ∆t →012r 2∆θ∆t =12r 2dθdt或2dA dt =r 2dθdt4这个就是开普勒第二定律的数学表达式.现在我们考虑开普勒第一定律,我们已知行星轨道是椭圆轨道.椭圆的极坐标方程为:r=p1+ecosθ令u=1ru=1p+epcosθ把这个式子带入比耐公式F=−mC2u2(d2udθ2+u)=−mC2u2p=−C2pmr2由于每个行星的C,p是不同的,所以我们还需要利用开普勒第三定律来确定C,p两边对2dAdt =r2dθdt积分.积分范围是全周期.2A=CT=2πabT=2πabC带入开普勒第三定律T2 a3=4π2b2C2a而b2 a =1a(a2−c2)=a(1−c2a2)=a(1−e2)=p即可以得到C2 p =4π2a3T2=k由开普勒第三定律知k是一个常数.与行星位置、质量无关F=−mk2 r2令k2=GMF(r)=−GMm r2牛顿万有引力公式得到证明.5。