地球引力模型程序代碼

1、L 148g模型地球部引力圖地球模型計算的引力外地核處引力最大，外地核向或向外均在減 小，從地核核心區(qū)處逃離地球引力需要能量為 W=3.822MgR,由于 不均勻分布，勢能增加了 0.5MgR 實際引力為分層球體的引力標準線疊加的結果。2. 96g2. 74g2. 500g2. 494g如果地球模型不錯我們完全可以根據(jù)分層密度計算它的引力圖半徑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把第二層密度減去第一層后繪圖,如此循環(huán)，畫到地核，再把所有值疊加起來，就是模型引力圖了。如果模型不錯，那么地球的實際引把密度確定以后，我們再把密度分層，先按最小密度繪制大圓圖，再Word文檔Word文檔力就會和所畫圖一樣。當然還必須要2a4a33(1 a2)211a2=3gc=22a4a33(1a2)21a2a sin , Rsin為質點到球心的距離,g。為地表重力加速度Word文檔Word文檔AUTOcadvba模型程序代碼Word文檔Sub地球引力模型()Dim ml As ObjectDim a As Variant, al As Variant, a2 As Varia

3、nt, a3 As Variant, a4 AsVaria 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.00520697211

4、556: a1=1.18178444724396:a2=1.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

5、) = iIf v1 <= 1 The nU = 750 * 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 八

6、2)八 1.5) / a1Elself v2 > 1 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

7、> a2 The nU2 = 0End IfIf v4 <= 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

8、八 4) / 3) / (1 + (1 -v5 八 2)八 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 + b