這裡蒐索程式師資訊,查找有用的技術資料
当前位置:首页 » 编程语言 » 如何用c语言求万有引力
扩展阅读
webinf下怎么引入js 2023-08-31 21:54:13
堡垒机怎么打开web 2023-08-31 21:54:11

如何用c语言求万有引力

发布时间: 2022-05-03 07:33:38

A. 如何推出万有引力定律

zyg2506146, 看来也是物理系的,估计比我具有更扎实的数学基本功,但对这个问题的认识还不太清晰。

万有引力定律的确是“猜”出来的。从开普勒第三定律推导太阳和地球之间引力满足F=GMm/R^2是严格的数学结论,但并不能说明有质量的物体之间都有这样的引力。
牛顿发觉地面上,比如使“苹果落地”的力,和天体间的力,都是一种满足平方反比的力,很自然地(或许当时的历史条件下是很大胆地?)猜测,这是同一种力,并且世间万物,天体又或普通物体,都有符合F=GMm/R^2的引力。
“牛顿万有引力公式是有严格的公式推导!!”,但万有引力定律,是猜的,不是可以从其他理论推导出的。
==============================================================

当然,楼主的意思,就是进行万有引力公式的公式推导。如果还是高中生,把轨道当成圆,从开普勒第三定律出发就行了;如果严格些,按实际的椭圆轨道来推导,那是相当麻烦的,利用比耐公式,可以从轨道方程推出万有引力F(r)的形式。

=================================================================
还是先把圆形轨迹时的近似推导给出来吧。对于高中生够用了。
证明:
开普勒第三定律r^3/T^2=C(C是常数)
万有引力F,形式未知,但一定等于向心力F=mr(2π/T)^2
带入1/T^2=C/ r^3
F= mr 4π^2 *(C/ r^3)= C’* m/ r^2
因为引力的对称性F= C” * M/ r^2
所以F= GMm/ r^2 G是常数

B. c语言问题:我要用这个程序求万有引力,该怎么改这个程序有什么问题

const G=6.77*10^-11改为
const G=6.77*10E-11;加分号
const G=6.77*10E-11; r后面多了个逗号。

C. 物体内部某质点的万有引力怎么求

为了求解这个问题,首先必须引入几何中“立体角“这个概念

=========================================

假设存在一个无穷小面元dS并且规定参考点O,那么这个面元在垂直于与O点连线方向上的投影面积与面元到O点距离r的平方的比值就定义为立体角Ω

矢量式运算满足dS×r=Ωr^3


可以类比于平面几何中的弧度角定义

=========================================

证明这个问题,实际就是证明 ”均匀厚球壳在内部产生场强处处为0“ 这个命题!


假设存在一个薄球壳,球心为O,面质量密度为σ,已知P在壳内部

现在由P点任意引出两个 对顶 的立体角

现在考察两个立体角截出的面元产生的引力效应

F1=GmσS1/r1^2

F2=GmσS2/r2^2

由立体角定义得出S=Ω*(r^2)/cosa(此处为标量式)

又因为等腰三角形底角相等

因此式子又可化为

F1=GmσΩ/cosa

F2=GmσΩ/cosa


这说明,两个面元的引力抵消了

又因为此处立体角是任意选取的,所以整个球壳对于P点的引力都抵消了

对薄球壳尚且如此,那么厚球壳也是一样了


所以对于匀质实心球体中任意一点,其外部的那层球壳影响完全忽略!质点受到的万有引力完全是由于内部的球体造成的

得证

D. c语言----编程

#include "stdio.h"
#include "math.h"
#include "malloc.h"
#define G 6.67 /*万有引力常数(G:m^3/(kg*s^2))*/
main ()
{
FILE *fp;
int i;
double ra,rc,fia,fic;
double *Vxz,*Vzz;
double p,b,h,q; //double变量用%lf输入可以用%f或%lf输出

printf("\n***********A program of Gravitytabular***********\n");
printf("Remain Dentisity is p\n");
printf("Breadth is b\n");
printf("Burial Depth is h\n");
printf("Dip is q\n");
printf("Please enter the parameters:\n");
scanf("%lf%lf%lf%lf",&p,&b,&h,&q); //变量的顺序和上面提示的顺序不吻合
Vxz=(double *)calloc(50,sizeof(double));
Vzz=(double *)calloc(50,sizeof(double));
fp=fopen("answer2.txt","w");

for(i=0;i<=49;i++)
{
ra=sqrt(pow((i+b),2)+pow(h,2));
rc=sqrt(pow((i-b),2)+pow(h,2));
fia=atan((b+i)/h);
fic=atan((i-b)/h);
Vxz[i]=2*G*p*sin(q)*(sin(q)*log(rc/ra)+cos(q)*(fia-fic));
Vzz[i]=2*G*p*sin(q)*(sin(q)*(fia-fic)-cos(q)*log(rc/ra));
Vxz[i]=10*Vxz[i];
Vzz[i]=10000*Vzz[i];
fprintf(fp,"Vxz=%10.4f Vzz=%10.4f\n",Vxz[i],Vzz[i]);
}
fclose(fp);
printf("\n end of transform.\n");
}
//如果还有疑问我们可以交流交流
//祝你学习进步

E. C语言!!!!求两个质量分别为m1,m2,距离为r的物体间的万有引力,输出计算结果

# include <stdio.h>
# include <math.h>
main()
{
double F,m1,m2,r,g;
g=6.67*10E-11;
printf("请输入两个质量m1,m2和r:\n");
scanf("%d,%d,%d",&m1,&m2,&r);
F=(g*m1*m2)/(r*r);
printf("万有引力:%lf\n",F); //double类型对应为%lf
}

F. 求详细解答:关于万有引力公式F=Gm

C

G. 高一物理的万有引力对公式的应用技巧和计算技巧。

一、质点的运动(1)------直线运动
1)匀变速直线运动
1.平均速度V平=s/t(定义式) 2.有用推论Vt2-Vo2=2as
3.中间时刻速度Vt/2=V平=(Vt+Vo)/2 4.末速度Vt=Vo+at
5.中间位置速度Vs/2=[(Vo2+Vt2)/2]1/2 6.位移s=V平t=Vot+at2/2=Vt/2t
7.加速度a=(Vt-Vo)/t {以Vo为正方向,a与Vo同向(加速)a>0;反向则a<0}
8.实验用推论Δs=aT2 {Δs为连续相邻相等时间(T)内位移之差}
9.主要物理量及单位:初速度(Vo):m/s;加速度(a):m/s2;末速度(Vt):m/s;时间(t)秒(s);位移(s):米(m);路程:米;速度单位换算:1m/s=3.6km/h。
注:
(1)平均速度是矢量;
(2)物体速度大,加速度不一定大;
(3)a=(Vt-Vo)/t只是量度式,不是决定式;
(4)其它相关内容:质点、位移和路程、参考系、时间与时刻〔见第一册P19〕/s--t图、v--t图/速度与速率、瞬时速度〔见第一册P24〕。
2)自由落体运动
1.初速度Vo=0 2.末速度Vt=gt
3.下落高度h=gt2/2(从Vo位置向下计算) 4.推论Vt2=2gh
注:
(1)自由落体运动是初速度为零的匀加速直线运动,遵循匀变速直线运动规律;
(2)a=g=9.8m/s2≈10m/s2(重力加速度在赤道附近较小,在高山处比平地小,方向竖直向下)。
(3)竖直上抛运动
1.位移s=Vot-gt2/2 2.末速度Vt=Vo-gt (g=9.8m/s2≈10m/s2)
3.有用推论Vt2-Vo2=-2gs 4.上升最大高度Hm=Vo2/2g(抛出点算起)
5.往返时间t=2Vo/g (从抛出落回原位置的时间)
注:
(1)全过程处理:是匀减速直线运动,以向上为正方向,加速度取负值;
(2)分段处理:向上为匀减速直线运动,向下为自由落体运动,具有对称性;
(3)上升与下落过程具有对称性,如在同点速度等值反向等。
二、质点的运动(2)----曲线运动、万有引力
1)平抛运动
1.水平方向速度:Vx=Vo 2.竖直方向速度:Vy=gt
3.水平方向位移:x=Vot 4.竖直方向位移:y=gt2/2
5.运动时间t=(2y/g)1/2(通常又表示为(2h/g)1/2)
6.合速度Vt=(Vx2+Vy2)1/2=[Vo2+(gt)2]1/2
合速度方向与水平夹角β:tgβ=Vy/Vx=gt/V0
7.合位移:s=(x2+y2)1/2,
位移方向与水平夹角α:tgα=y/x=gt/2Vo
8.水平方向加速度:ax=0;竖直方向加速度:ay=g注:
(1)平抛运动是匀变速曲线运动,加速度为g,通常可看作是水平方向的匀速直线运与竖直方向的自由落体运动的合成;
(2)运动时间由下落高度h(y)决定与水平抛出速度无关;
(3)θ与β的关系为tgβ=2tgα;
(4)在平抛运动中时间t是解题关键;(5)做曲线运动的物体必有加速度,当速度方向与所受合力(加速度)方向不在同一直线上时,物体做曲线运动。
2)匀速圆周运动
1.线速度V=s/t=2πr/T 2.角速度ω=Φ/t=2π/T=2πf
3.向心加速度a=V2/r=ω2r=(2π/T)2r 4.向心力F心=mV2/r=mω2r=mr(2π/T)2=mωv=F合
5.周期与频率:T=1/f 6.角速度与线速度的关系:V=ωr
7.角速度与转速的关系ω=2πn(此处频率与转速意义相同)
8.主要物理量及单位:弧长(s):米(m);角度(Φ):弧度(rad);频率(f):赫(Hz);周期(T):秒(s);转速(n):r/s;半径(r):米(m);线速度(V):m/s;角速度(ω):rad/s;向心加速度:m/s2。
注:
(1)向心力可以由某个具体力提供,也可以由合力提供,还可以由分力提供,方向始终与速度方向垂直,指向圆心;
(2)做匀速圆周运动的物体,其向心力等于合力,并且向心力只改变速度的方向,不改变速度的大小,因此物体的动能保持不变,向心力不做功,但动量不断改变。
3)万有引力
1.开普勒第三定律:T2/R3=K(=4π2/GM){R:轨道半径,T:周期,K:常量(与行星质量无关,取决于中心天体的质量)}
2.万有引力定律:F=Gm1m2/r2 (G=6.67×10-11N?m2/kg2,方向在它们的连线上)
3.天体上的重力和重力加速度:GMm/R2=mg;g=GM/R2 {R:天体半径(m),M:天体质量(kg)}
4.卫星绕行速度、角速度、周期:V=(GM/r)1/2;ω=(GM/r3)1/2;T=2π(r3/GM)1/2{M:中心天体质量}
5.第一(二、三)宇宙速度V1=(g地r地)1/2=(GM/r地)1/2=7.9km/s;V2=11.2km/s;V3=16.7km/s
6.地球同步卫星GMm/(r地+h)2=m4π2(r地+h)/T2{h≈36000km,h:距地球表面的高度,r地:地球的半径}
注:
(1)天体运动所需的向心力由万有引力提供,F向=F万;
(2)应用万有引力定律可估算天体的质量密度等;
(3)地球同步卫星只能运行于赤道上空,运行周期和地球自转周期相同;
(4)卫星轨道半径变小时,势能变小、动能变大、速度变大、周期变小(一同三反);
(5)地球卫星的最大环绕速度和最小发射速度均为7.9km/s。
三、力(常见的力、力的合成与分解)
1)常见的力
1.重力G=mg (方向竖直向下,g=9.8m/s2≈10m/s2,作用点在重心,适用于地球表面附近)
2.胡克定律F=kx {方向沿恢复形变方向,k:劲度系数(N/m),x:形变量(m)}
3.滑动摩擦力F=μFN {与物体相对运动方向相反,μ:摩擦因数,FN:正压力(N)}
4.静摩擦力0≤f静≤fm (与物体相对运动趋势方向相反,fm为最大静摩擦力)
5.万有引力F=Gm1m2/r2 (G=6.67×10-11N?m2/kg2,方向在它们的连线上)
6.静电力F=kQ1Q2/r2 (k=9.0×109N?m2/C2,方向在它们的连线上)
7.电场力F=Eq (E:场强N/C,q:电量C,正电荷受的电场力与场强方向相同)
8.安培力F=BILsinθ (θ为B与L的夹角,当L⊥B时:F=BIL,B//L时:F=0)
9.洛仑兹力f=qVBsinθ (θ为B与V的夹角,当V⊥B时:f=qVB,V//B时:f=0)注:
(1)劲度系数k由弹簧自身决定;
(2)摩擦因数μ与压力大小及接触面积大小无关,由接触面材料特性与表面状况等决定;
(3)fm略大于μFN,一般视为fm≈μFN;
(4)其它相关内容:静摩擦力(大小、方向)〔见第一册P8〕;
(5)物理量符号及单位B:磁感强度(T),L:有效长度(m),I:电流强度(A),V:带电粒子速度(m/s),q:带电粒子(带电体)电量(C);
(6)安培力与洛仑兹力方向均用左手定则判定。
2)力的合成与分解
1.同一直线上力的合成同向:F=F1+F2, 反向:F=F1-F2 (F1>F2)
2.互成角度力的合成:
F=(F12+F22+2F1F2cosα)1/2(余弦定理) F1⊥F2时:F=(F12+F22)1/2
3.合力大小范围:|F1-F2|≤F≤|F1+F2|
4.力的正交分解:Fx=Fcosβ,Fy=Fsinβ(β为合力与x轴之间的夹角tgβ=Fy/Fx)
注:
(1)力(矢量)的合成与分解遵循平行四边形定则;
(2)合力与分力的关系是等效替代关系,可用合力替代分力的共同作用,反之也成立;
(3)除公式法外,也可用作图法求解,此时要选择标度,严格作图;
(4)F1与F2的值一定时,F1与F2的夹角(α角)越大,合力越小;
(5)同一直线上力的合成,可沿直线取正方向,用正负号表示力的方向,化简为代数运算。四、动力学(运动和力)
1.牛顿第一运动定律(惯性定律):物体具有惯性,总保持匀速直线运动状态或静止状态,直到有外力迫使它改变这种状态为止
2.牛顿第二运动定律:F合=ma或a=F合/ma{由合外力决定,与合外力方向一致}
3.牛顿第三运动定律:F=-F?{负号表示方向相反,F、F?各自作用在对方,平衡力与作用力反作用力区别,实际应用:反冲运动}
4.共点力的平衡F合=0,推广 {正交分解法、三力汇交原理}
5.超重:FN>G,失重:FN<G {加速度方向向下,均失重,加速度方向向上,均超重}
6.牛顿运动定律的适用条件:适用于解决低速运动问题,适用于宏观物体,不适用于处理高速问题,不适用于微观粒子〔见第一册P67〕
注:平衡状态是指物体处于静止或匀速直线状态,或者是匀速转动。
五、振动和波(机械振动与机械振动的传播)
1.简谐振动F=-kx {F:回复力,k:比例系数,x:位移,负号表示F的方向与x始终反向}
2.单摆周期T=2π(l/g)1/2 {l:摆长(m),g:当地重力加速度值,成立条件:摆角θ<100;l>>r}
3.受迫振动频率特点:f=f驱动力
4.发生共振条件:f驱动力=f固,A=max,共振的防止和应用〔见第一册P175〕
5.机械波、横波、纵波〔见第二册P2〕
6.波速v=s/t=λf=λ/T{波传播过程中,一个周期向前传播一个波长;波速大小由介质本身所决定}
7.声波的波速(在空气中)0℃:332m/s;20℃:344m/s;30℃:349m/s;(声波是纵波)
8.波发生明显衍射(波绕过障碍物或孔继续传播)条件:障碍物或孔的尺寸比波长小,或者相差不大
9.波的干涉条件:两列波频率相同(相差恒定、振幅相近、振动方向相同)
10.多普勒效应:由于波源与观测者间的相互运动,导致波源发射频率与接收频率不同{相互接近,接收频率增大,反之,减小〔见第二册P21〕}
注:
(1)物体的固有频率与振幅、驱动力频率无关,取决于振动系统本身;
(2)加强区是波峰与波峰或波谷与波谷相遇处,减弱区则是波峰与波谷相遇处;
(3)波只是传播了振动,介质本身不随波发生迁移,是传递能量的一种方式;
(4)干涉与衍射是波特有的;
(5)振动图象与波动图象;
(6)其它相关内容:超声波及其应用〔见第二册P22〕/振动中的能量转化〔见第一册P173〕。
六、冲量与动量(物体的受力与动量的变化)
1.动量:p=mv {p:动量(kg/s),m:质量(kg),v:速度(m/s),方向与速度方向相同}
3.冲量:I=Ft {I:冲量(N?s),F:恒力(N),t:力的作用时间(s),方向由F决定}
4.动量定理:I=Δp或Ft=mvt–mvo {Δp:动量变化Δp=mvt–mvo,是矢量式}
5.动量守恒定律:p前总=p后总或p=p’?也可以是m1v1+m2v2=m1v1?+m2v2?
6.弹性碰撞:Δp=0;ΔEk=0 {即系统的动量和动能均守恒}
7.非弹性碰撞Δp=0;0<ΔEK<ΔEKm {ΔEK:损失的动能,EKm:损失的最大动能}
8.完全非弹性碰撞Δp=0;ΔEK=ΔEKm {碰后连在一起成一整体}
9.物体m1以v1初速度与静止的物体m2发生弹性正碰:
v1?=(m1-m2)v1/(m1+m2) v2?=2m1v1/(m1+m2)
10.由9得的推论-----等质量弹性正碰时二者交换速度(动能守恒、动量守恒)
11.子弹m水平速度vo射入静止置于水平光滑地面的长木块M,并嵌入其中一起运动时的机械能损失
E损=mvo2/2-(M+m)vt2/2=fs相对 {vt:共同速度,f:阻力,s相对子弹相对长木块的位移}
注:
(1)正碰又叫对心碰撞,速度方向在它们“中心”的连线上;
(2)以上表达式除动能外均为矢量运算,在一维情况下可取正方向化为代数运算;
(3)系统动量守恒的条件:合外力为零或系统不受外力,则系统动量守恒(碰撞问题、爆炸问题、反冲问题等);
(4)碰撞过程(时间极短,发生碰撞的物体构成的系统)视为动量守恒,原子核衰变时动量守恒;
(5)爆炸过程视为动量守恒,这时化学能转化为动能,动能增加;(6)其它相关内容:反冲运动、火箭、航天技术的发展和宇宙航行〔见第一册P128〕。

H. 求大神高手把下面这个公式用C语言编写出来,是在是弄的头昏了,多谢!(在线等待)

delta g(t) = 2 * pi * G * σ * Δθ * (Ri+1 - Ri + sqrt(Ri*Ri + h*h) - sqrt(Ri+1 * Ri+1 + h*h) )
难道是要具体化出来?

I. 万有引力解析。。。。。

万有引力定律
编辑

万有引力即万有引力定律,更多含义,请参阅万有引力(多义词)。
万有引力定律(Law of universal gravitation)是艾萨克·牛顿在1687年于《自然哲学的数学原理》上发表的解释物体之间相互作用的引力的定律。定律内容如下:任意两个质点通过连心线方向上的力相互吸引。该引力的大小与它们的质量乘积成正比,与它们距离的平方成反比,与两物体的化学本质或物理状态以及中介物质无关,公式表示:F=G*M1M2/(R*R)(G=6.67×10^-11N·m^2/kg^2)。
牛顿在前人(开普勒、胡克、雷恩、哈雷)研究的基础上,凭借他超凡的数学能力证明,在1687年于《自然哲学的数学原理》上发表了万有引力定律。万有引力定律的发现,是17世纪自然科学最伟大的成果之一。
目录
1推理依据
2引力规律
牛顿的猜想
猜想的依据
检验的思想
检验的结果
3详细内容
4公式表示
适用范围
意义
G的数值
5其它关系
重力加速度
引力~重力
匀速圆周运动
6存在问题
简介
理论问题
观测问题
局限性
7参考理论
过往理论
最新理论
8应用
1推理依据

科学家牛顿
伽利略在1632年实际上已经提出离心力和向心力的初步想法。布里阿德在1645年提出了引力平方比关系的思想.牛顿在1665—1666年的手稿中,用自己的方式证明了离心力定律,但向心力这个词可能首先出现在《论运动》的第一个手稿中。一般人认为离心力定律是惠更斯在1673年发表的《摆钟》一书中提出来的。根据1684年8月-10月间牛顿写的《论回转物体的运动》一文手稿中,牛顿很可能在这个手稿中第一次提出向心力及其定义。
万有引力与相作用的物体的质量乘积成正比,是发现引力平方反比定律过渡到发现万有引力定律的必要阶段.·牛顿从1665年至1685年,花了整整20年的时间,才眼中离心力—向心力—重力—万有引力概念的演化顺序,终于提出“万有引力”这个概念和词汇。·牛顿在《原理》第三卷中写道:最后,如果由实验和天文学观测,普遍显示出地球周围的一切天体被地球重力所吸引,并且其重力与它们各自含有的物质之量成比例,则月球同样按照物质之量被地球重力所吸引。另一方面,它显示出,我们的海洋被月球重力所吸引;并且一切行星相互被重力所吸引,彗星同样被太阳的重力所吸引。由于这个规则,我们必须普遍承认,一切物体,不论是什么,都被赋与了相互的引力(gravitation)的原理。因为根据这个表象所得出的一切物体的万有引力(universal gravitation)的论证……".
牛顿在1665—1666年间只用离心力定律和开普勒第三定律,因而只能证明圆轨道上的而不是椭圆轨道上的引力平方反比关系。在1679年,他知道运用开普勒第二定律,但是在证明方法上没有突破,仍停留在1665—1666年的水平。只是到了1684年1月,哈雷、雷恩、胡可和牛顿都能够证明圆轨道上的引力平方反比关系,都已经知道椭圆轨道上遵守引力平方反比关系,但是最后可能只有牛顿才根据开普勒三个定律、从离心力定律演化出的向心力定律和数学上的极限概念或微积分概念,才用几何法证明了这个难题。
2引力规律
牛顿的猜想
地球与太阳之间的吸引力与地球对周围物体的引力可能是同一种力,遵循相同的规律。
猜想的依据
(1)行星与太阳之间的引力使行星不能飞离太阳,物体与地球之间的引力使物体不能离开地球;(2)在离地面很高的距离里,都不会发现重力有明显的减弱,那么这个力必然延伸到很远的地方。
检验的思想
如果猜想正确,月球在轨道上运动的向心加速度与地面重力加速度的比值,应该等于地球半径平方与月球轨道半径平方之比,即


检验的结果
地面物体所受地球的引力,与月球所受地球的引力是同一种力。
3详细内容
自然界中任何两个物体都是相互吸引的,引力的大小与两物体的质量的乘积成正比,与两物体间距离的平方成反比。
4公式表示



F: 两个物体之间的引力
G:万有引力常量
m1: 物体1的质量
m2: 物体2的质量
r: 两个物体之间的距离
依照国际单位制,F的单位为牛顿(N),m1和m2的单位为千克(kg),r 的单位为米(m),常数G近似地等于

(牛顿平方米每二次方千克)。
由此可知排斥力F一直都将不存在,这意味着净加速度的力是绝对的。(这个符号规约是为了与库仑定律相容而订立的,在库仑定律中绝对的力表示两个电子之间的作用力。)
适用范围
经典万有引力定律反映了一定历史阶段人类对引力的认识,在十九世纪末发现,水星在近日点的移动速度比理论值大,即发现水星轨道有旋紧,轨道旋紧的快慢的实际值为每世纪42.9″。这种现象用万有引力定律无法解释,而根据广义相对论计算的结果旋紧是每世纪43.0″,在观测误差允许的范围内。此外,广义相对论还能较好地解释谱线的红移和光线在太阳引力作用下的偏转等现象。这表明广义相对论的引力理论比经典的引力理论进了一步。
在法拉第和麦克斯韦之后,人们看到物理的实在除了粒子还有场。电磁场具有动量和能量且能传播电磁波。这使人们联想万有引力定律也是物理的实在,能传播引力波,也有许多人努力探测它,但尚无很好的结果。电磁波的传播可用光子解释,类似地,光子也导致引力子概念的引出。万有引力也不再是超距作用,而以引力子为媒介。但这些都是物理学家正在探索的领域。
经典力学的适用范围并引入普朗克常量和真空中光速来界定经典力学的领地。粗糙的说,经典的万有引力定律适用范围也可用一数量表示。现在引入引力半径

,G、m分别表示引力常量和产生引力场的球体的球体的质量,c为光速。用R表示产生力场球体之半径,若

,则可用牛顿引力定律。对于太阳,

,应用牛顿引力定律无问题;即使是对致密的白矮星,

,也仍然可用牛顿万有引力定律;至于中子星,

,这便有必要引用广义相对论。至于黑洞和宇宙大爆炸,应当是应用广义相对论的。
意义
万有引力定律的发现,是17世纪自然科学最伟大的成果之一。它把地面上物体运动的规律和天体运动的规律统一了起来,对以后物理学和天文学的发展具有深远的影响。它第一次解释了(自然界中四种相互作用之一)一种基本相互作用的规律,在人类认识自然的历史上树立了一座里程碑。
万有引力定律揭示了天体运动的规律,在天文学上和宇宙航行计算方面有着广泛的应用。它为实际的天文观测提供了一套计算方法,可以只凭少数观测资料,就能算出长周期运行的天体运动轨道,科学史上哈雷彗星、海王星、冥王星的发现,都是应用万有引力定律取得重大成就的例子。利用万有引力公式,开普勒第三定律等还可以计算太阳、地球等无法直接测量的天体的质量。牛顿还解释了月亮和太阳的万有引力引起的潮汐现象。他依据万有引力定律和其他力学定律,对地球两极呈扁平形状的原因和地轴复杂的运动,也成功的做了说明。推翻了古代人类认为的神之引力。
对文化发展有重大意义:使人们建立了有能力理解天地间的各种事物的信心,解放了人们的思想,在科学文化的发展史上起了积极的推动作用。
G的数值
牛顿在推出万有引力定律的同时,并没能得出引力常量G的具体值。G的数值于1789年由卡文迪许利用他所发明的扭秤得出。卡文迪许的扭秤试验,不仅以实践证明了万有引力定律,同时也让此定律有了更广泛的使用价值。
卡文迪许测出的

,与现在的公认值

极为接近;直到1969年G的测量精度还保持在卡文迪许的水平上。
5其它关系
重力加速度
令a1为事先已知质点的重力加速度。由牛顿第二定律知, 即。取代前面方程中的F
同理亦可得出a2.
依照国际单位制,重力加速度(同其他一般加速度)的单位被规定为米每平方秒 (m/s^2 或 m s^三)。非国际单位制的单位有伽利略、单位g(见后)以及 英尺每秒的平方。
请注意上述方程中的a1,质量m1的加速度,在实际上并不取决于m1的取值。因此可推论出对于任何物体,无论它们的质量为多少,它们都将按照同样的比率向地面坠落(忽略空气阻力)。
如果物体运动过程中r只有极微小的改变——譬如地面附近的自由落体运动——重力加速度将几乎保持不变(参看条目地心引力)。而对于一个庞大物体,由于r的变化导致的不同位点所受重力的变化,将会引起巨大而可观的潮汐力作用。
令m1为地球质量5.98*10^24kg,m2为1kg,R为地球半径6380000m,代入万有引力公式,计算出F=9.8N,这说明1kg的物体在地球表面受重力为9.8N。换句话说,等式两边同除以m2,结果就是重力加速度g。
具有空间广度的物体:
如果被讨论的物体具有空间广度(远大于理论上的质点),它们之间的万有引力可以以物体的各个等效质点所受万有引力之和来计算。在极限上,当组成质点趋近于“无限小”时,将需要求出两物体间的力(矢量式见下文)在空间范围上的积分。
从这里可以得出:如果物体的质量分布呈现均匀球状时,其对外界物体施加的万有引力吸引作用将同所有的质量集中在该物体的几何中心原理时的情况相同。(这不适用于非球状对称物体)。
矢量式:
地球附近空间内的重力示意图:在此数量级上地球表面的弯曲可被忽略不计,因此力线可以近似地相互平行并且指向地球的中心牛顿万有引力定律亦可通过矢量方程的形式进行表述而用以计算万有引力的方向和大小。在下列公式中,以粗体显示的量代表矢量。

地球的重力示意图
其中:
F12: 物体1对物体2的引力
G: 万有引力常量
m1与m2: 分别为物体1和物体2的质量
r21 = | r2 r1 |: 物体2和物体1之间的距离
r21= r1+r2 物体2和物体1之间的距离
: 物体1到物体2的单位矢量
可以看出矢量式方程的形式与之前给出的标量式方程相类似,区别仅在于在矢量式中的F是一个矢量,以及在矢量式方程的右端被乘上了相应的单位向量。而且,我们可以看出:F12 = F21.
同样,重力加速度的矢量式方程与其标量式方程相类似。
引力~重力
1.重力是由于地球的吸引而产生的,但能否说万有引力就是重力呢?分析这个问题应从地球自转入手。由于地球自转,地球上的物体随之做圆周运动,所受的向心力F1=mrw^2=mRw^2cosa,F1是引力F提供的,它是F的一个分力,cosa是引力F与赤道面的夹角的余弦值,F的另一个分力F2就是物体所受的重力,即F2=mg。
由此可见,地球对物体的万有引力是物体受到重力的原因,但重力不完全等于万有引力,这是因为物体随地球自转,需要有一部分万有引力来提供向心力。
2.重力与万有引力间的大小关系
(1)重力与纬度的关系
在赤道上满足mg=F-F向(物体受万有引力和地面对物体的支持力Fn的作用,其合力充当向心力,Fn的大小等于物体的重力的大小)。
在地球两极处,由于F向=0,即mg=F,在其他位置,mg、F与F向 间符合平行四边形定则。同一物体在赤道处重力最小,并随纬度的增加而增大。
(2)重力、重力加速度与高度的关系
在距地面高度为h的高处,若不考虑地球自转的影响时,则mg'=F=GMm/(R+h)^2;而在地面处mg=GMm/R^2.
距地面高为h处,其重力加速度g'=GM/(R+h)^2,在地面处g=GM/R^2.
在距地面高度为h的轨道上运行的宇宙飞船中,质量为m的物体的重力即为该处受到的万有引力,即mg'=GmM/(R+h)^2,但无法用测力计测出其重力。
匀速圆周运动
一个天体环绕另一个中心天体做匀速圆周运动。其向心力由万有引力提供。即F引=GMm/r^2≈mg=ma向,而a向=v^2/r=w^2r=vw=(4π^2/T^2)r=4π^2f^2r,因此应用万有引力定律解决天体的有关问题,主要有以下几个度量关系:F引=GMm/r^2(r为轨道半径)=mg=ma向=mv^2/r=mw^2r=m(4π^2/T^2)r=m4π^2f^2r.
重力场:
球状星团 M13 证明重力场的存在。重力场是用于描述在任意空间内某一点的物体每单位质量所受万有引力的矢量场。而在实际上等于该点物体所受的重力加速度。
以下是一个普适化的矢量式,可被应用于多于两个物体的情况(例如在地球与月球之间穿行的火箭)的计算。对于两个物体的情况(比如说物体1是火箭,物体2是地球)来说,我们可以用 替代并用m替代m1来将重力场表示为:
因此我们可以得到:
该公式不受产生重力场的物体的限制。重力场的单位为力除以质量的单位;在国际单位制上,被规定为N·kgㄢ(牛顿每千克)。
6存在问题
简介
尽管牛顿对重力的描述对于众多实践运用来说十分地精确,但它也具有几大理论问题且被证明是不完全正确的。
理论问题
没有任何征兆表明重力的传送媒介可以被识别出,牛顿自己也对这种无法说明的超距作用感到不满意(参看后文条目“局限性”)。
牛顿的理论需要定义重力可以瞬时传播。因此给出了古典自然时空观的假设,这样亦能使约翰内斯·开普勒所观测到的角动量守恒成立。但是,这与爱因斯坦的狭义相对论理论有直接的冲突,因为狭义相对论定义了速度的极限——真空中的光速——在此速度下信号可以被传送。
观测问题
牛顿的理论并不能完全地解释出水星在沿其轨道运动到近日点时出现的进动现象。牛顿学说的预言(由其它行星的重力拖曳产生)与实际观察到的进动相比每世纪会出现43弧秒的误差。
牛顿的理论预言的重力作用下光线的偏折只有实际观测结果的一半。广义相对论则与观察结果更为接近。
所有物体的重力质量与惯性质量相同的这一观测现象是牛顿的系统所不能解释的。广义相对论则将它作为一个基本条件。参看条目等效原理。
局限性
当牛顿非凡的工作使万有引力定律能够为数学公式所表示后,他仍然不满于公式中所隐含的“超距作用”观点。他从来没有在他的文字中“赋予产生这种能力的原因”。在其它情况下,他使用运动的现象来解释物体受到不同力的作用的原因,但是对于重力这种情况,他却无法用实验方法来确认运动产生了重力。此外,他甚至还拒绝对这个由地面产生的力的起因提出假设,而这一切都违背了科学证据的原则。
牛顿的经典力学只适用于低速、宏观、弱引力,而不适用于高速、微观与强引力。
牛顿对重力的发现埋葬了“哲学家至今仍在愚蠢地试图探索自然”(philosophers have hitherto attempted the search of nature in vain)这句所谓的真理,就同他深信着的“有各种因素”使得“各种迄今未知的原因”是所有“自然现象”的基础。这些基本的现象至今仍在研究中,而且,虽然存在着许多种的假设,最终答案仍然没有找出。 虽然爱因斯坦的假设的确比牛顿的假设更能精确地解释确定案例中万有引力的作用效果,但是他也从来没有在他的理论中为这种能力赋予一个原因。在爱因斯坦的方程式中,“物质告诉空间怎么扭曲,空间告诉物质怎么移动”(matter tells space how to curve, and space tells matter how to move),但是这个完全异于牛顿世界的新的思想,也不能使爱因斯坦所赋予“产生这种能力的原因”比万有引力定律使牛顿所赋予的原因更能使空间产生扭曲。牛顿自己说:
我还没有能力去从现象中发现产生这些重力特性的原因,而且我无法臆测……我所解释的定律和丰富的天体运动的计算已经足够于说明重力的确存在并能产生效果。一个物体可以不通过任何介质穿过真空间的距离对另一个物体产生作用,在此之上它们的活动和力可以传送自对方,这对于我来说简直就是一个天大的谬论。因此,我相信,任何有足够的哲学思维能力的人都不会沉溺于此。I have not yet been able to discover the cause of these properties of gravity from phenomena and I feign no hypotheses... It is enough that gravity does really exist and acts according to the laws I have explained, and that it abundantly serves to account for all the motions of celestial bodies. That one body may act upon another at a distance through a vacuum without the mediation of anything else, by and through which their action and force may be conveyed from one another, is to me so great an absurdity that, I believe, no man who has in philosophic matters a competent faculty of thinking could ever fall into it.
需注意的是,这里使用的单词“原因(cause)”并不是“起因(cause)和影响”或者“被告导致(cause)受害者死亡”中所表示的意义。何况,当牛顿使用单词“原因(cause)”时,他(明显地)意指为一种“解释”。或者说,像“牛顿学说的重力是行星运动的原因”这个短语的意思就是牛顿学说的重力解释了行星的运动。
7参考理论
过往理论
亚里士多德引力理论 亚里士多德认为,物体的运动速度和其所受外界的合力是成正比(或者是该物体所受的自己本身的引力),并且和物体运动介质的粘度成反比。
尼古拉·特斯拉(Nikola Tesla)宣布但是从未发表的引力动力学理论;部分原因是因为理论的细节(如果有的话)并没有透露,并没有得到物理学家们的重视。
感应引力(Inced Gravity),由安德烈·萨哈罗夫(Andrei Sakharov)提出,认为广义相对论可能起源于量子场论。
雷萨吉万有引力理论(Le Sage's Theory of Gravitation)(也叫做雷萨吉引力理论),由乔治-路易斯·雷萨吉(Georges-Louis Le Sage)提出,以一种充满整个宇宙轻的气体的流动来解释这种现象。
万有引力理论(Nordström's Theory of Gravitation),广义相对论的早期竞争者。
怀特黑德万有引力理论,(Whitehead's Theory of Gravitation)广义相对论的另一个早期竞争者。
最新理论
布兰斯-迪克(Brans-Dicke)有关引力的理论(Brans-Dicke theory)。
Rosen有关引力的理论(Rosen Bi-metric Theory)。
莫德采·米尔格若姆(Mordehai Milgrom)在修正牛顿引力理论(Modified Newtonian Dynamics, MOND)中,提出在微小加速运动上对牛顿第二定律的修正。
新近提出的且被高度争论的程序物理学(Process Physics)理论试图处理(address)引力问题。
自建宇宙理论(Self Creation Cosmology)将布兰斯-迪克理论修正为允许创造质量。
8应用
1.计算天体质量
(1)计算地球质量
若不考虑地球自转,地面上物体所受重力即地球对它的万有引力
mg=GmM/R^2 由此可得地球质量 M=gR^2/G
(2)计算太阳质量
测量地球绕太阳公转周期,公转轨道半径,将轨道看成圆,匀速圆周运动向心力就是万有引力
即 GMm/R^2=m(2π/T)^2 R 地球质量为m, 太阳质量 M=4π^2R^3/GT^2
运用类似方法已知人造卫星质量,卫星绕某天体运动的周期和轨道半径
可算出天体质量
2.估算天体密度
若设某天体半径R,卫星绕天体表面运行时,轨道半径为R,
又测得已知运行周期为T
设卫星质量为m 则 GMm/R^2=m(2π/T)^2R 天体质量M=4π^2R^3/GT^2
体积V=4πR^3/3 ρ=M/V=3π/GT^2

J. 如何用C语言计算公式

C语言计算公式通常需要#include math.h的数学库