地球引力模型程序代码

1、模型地球内部引力图 地球模型计算的引力外地核处引力最大 ,外地核向内或向外均在 减小，从地核核心区处逃离地球引力需要能量为 W=3.822MgR, 由于不均匀分布，势能增加了 0.5MgR 实际引力为分层球体的引力标准线叠加的结果。 密度 2.75 3.4 5.35 10.7 12.3 12.5 分层密度 2.75 0.65 1.95 5.35 1.6 0.2 叠加起来，就是模型引力图了。 。如果模型不错，那么地球的实际引 力就会和所画图一样。 当然还必须要用到内部引力公式 4 2 a4 1 a2 2a 3 3 3 (1 a2) 2 1 g(a) 3MG 2R2 2a a2 (1 a2)2 3

2、 = 3g0 1 如果地球模型不错我们完全可以根据分层密度计算它的内 引力图 半径 R 1 0.99482 0.846178 0.544812 0.262282 0.194318 R/N 1 1.005207 1.181784 1.835496 3.81269 5.146204 把密度确定以后，我们再把密度分层，先按最小密度绘制大圆图，再 把第二层密度减去第一层后绘图，如此循环，画到地核，再把所有值 a sin ,Rsin 为质点到球心的距离 , g0为地表重力加速度 AUTOcadvba 模型程序代码 Sub 地球引力模型 () Dim ml As Object Dim a As Varia

3、nt, a1 As Variant, a2 As Variant, a3 As Variant, a4 As Variant Dim v As Variant, v1 As Variant, v2 As Variant, v3 As Variant, v4 As Variant, v5 As Variant Dim U As Variant, U1 As Variant, U2 As Variant, U3 As Variant, U4 As Variant Dim co As Variant co = 3 Dim p(0 To 20001) As Double Dim i As Intege

4、r Dim J As Variant a = 1.556: a1 = 1.18178444724396: a2 = 1.8354955470878: a3 = 3.81269015792163: a4 = 5.14620364557066 b = 0.236363636363636: b1 = 0.7909: b2 = 1.94545454545455: b3 = 0.581818181818182: b4 = 7.27272727272727E-02 For i = 0 To 10000 Step 2 v = i / 1000: v1 = v * a: v2 = v * a1: v3 = v

5、 * a2: v4 = v * a3: v5 = v * a4 p(i) = i If v1 1 And v1 a Then U3 = 0 End If If v2 1 And v2 a1 Then U3 = 0 End If If v3 1 And v3 a2 Then U2 = 0 End If If v4 1 And v4 a3 Then U3 = 0 End If If v5 1 And v5 a4 Then U4 = 0 End If If v = 1 Then p(i + 1) = 750 * (1 + v * 2 - (1 - v 2 + (v 4) / 3) / (1 + (1 - v 2) 1.5) + U + U1 + U2 + U3 + U4 Else p(i + 1) = 500 * (1 + b * a -3 + b1 * a1 -3 + b2 * a2 -3 + b3 * a3 -3 + b4 * a4 -3) / v 2 End If Next i Set m