地球引力模型程序代码

1、L 148g模型地球部引力图地球模型计算的引力外地核处引力最大，外地核向或向外均在减 小，从地核核心区处逃离地球引力需要能量为 W=3.822MgR,由于 不均匀分布，势能增加了 0.5MgR 实际引力为分层球体的引力标准线叠加的结果。2. 96g2. 74g2. 500g2. 494g如果地球模型不错我们完全可以根据分层密度计算它的引力图半径R10.994820.8461780.5448120.2622820.194318R/N11.0052071.1817841.8354963.812695.146204密度2.753.45.3510.712.312.5分层密度2.750.651.955.

2、351.60.2把第二层密度减去第一层后绘图,如此循环，画到地核，再把所有值叠加起来，就是模型引力图了。如果模型不错，那么地球的实际引把密度确定以后，我们再把密度分层，先按最小密度绘制大圆图，再力就会和所画图一样。当然还必须要2a4a33(1 a2)211a2=3gc=22a4a33(1a2)21a2a sin , Rsin为质点到球心的距离,g。为地表重力加速度AUTOcadvba模型程序代码Sub地球引力模型()Dim ml As ObjectDim a As Variant, al As Variant, a2 As Variant, a3 As Variant, a4 AsVaria

3、ntDim v As Variant, v1 As Variant, v2 As Variant, v3 As Variant, v4 AsVaria nt, v5 As Varia ntDim U As Variant, U1 As Variant, U2 As Variant, U3 As Variant, U4 AsVaria ntDim co As Varia ntco = 3Dim p(0 To 20001) As DoubleDim i As In tegerDim J As Varia nta=1.00520697211556: a1=1.18178444724396:a2=1.

4、8354955470878: a3 = 3.81269015792163: a4 = 5.14620364557066b=0.236363636363636: b1=0.709090909090909:b2=1.94545454545455:b3 =0.581818181818182:b4 =7.27272727272727E-02For i = 0 To 10000 Step 2v = i / 1000: v1 = v * a: v2 = v * a1: v3 = v * a2: v4 = v * a3: v5 = v *a4P(i) = iIf v1 <= 1 The nU = 75

5、0 * b * (1 + v1 * 2 - (1 - v1 八 2 + (v1 八 4) / 3) / (1 + (1 - v1 八 2) 八 1.5)/ aElself v1 > 1 And v1 <= a The nU = 500 * b * a 八-3 / v 八 2Elself v1 > a The nU3 = 0End IfIf v2 <= 1 The nU1 = 750 * b1 * (1 + v2 * 2 - (1 - v2 八 2 + (v2 八 4) / 3) / (1 + (1 - v2 八2)八 1.5) / a1Elself v2 > 1

6、And v2 <= a1 The nU1 = 500 * b1 * a1 八-3 / v 八 2Elself v2 > a1 The nU3 = 0End IfIf v3 <= 1 The nU2 = 750 * b2 * (1 + v3 * 2 - (1 - v3 八 2 + (v3 八 4) / 3) / (1 + (1 - v3 八2)八 1.5) / a2ElseIf v3 > 1 And v3 <= a2 The nU2 = 500 * b2 * a2 八-3 / v 八 2ElseIf v3 > a2 The nU2 = 0End IfIf v4

7、 <= 1 The nU3 = 750 * b3 * (1 + v4 * 2 - (1 - v4 八 2 + (v4 八 4) / 3) / (1 + (1 - v4 八2)八 1.5)/ a3Elself v4 > 1 And v4 <= a3 The nU3 = 500 * b3 * a3 八-3 / v 八 2Elself v4 > a3 The nU3 = 0End IfIf v5 <= 1 The nU4 = 750 * b4 * (1 + v5 * 2 - (1 - v5 八 2 + (v5 八 4) / 3) / (1 + (1 -v5 八 2)八

8、1.5) / a4Elself v5 > 1 And v5 <= a4 The nU4 = 500 * b4 * a4 八-3 / v 八 2Elself v5 > a4 The nU4 = 0End IfIf v <= 1 The np(i + 1) = 750 * (1 + v * 2 - (1 - v 八 2 + (v 八 4) / 3) / (1 + (1 - v 八 2)八1.5) + U + U1 + U2 + U3 + U4Elsep(i + 1) = 500 * (1 + b * a 八-3 + bl * al 八-3 + b2 * a2 八-3 + b3a3 八-3