地球引力模型程序代碼剖析

1、模型地球內(nèi)部引力圖地球模型計算的引力外地核處引力最大 ，外地核向內(nèi)或向外均在 減小，從地核核心區(qū)處逃離地球引力需要能量為 W=3.822MgR, 由于不均勻分布，勢能增加了 0.5MgR實際引力為分層球體的引力標(biāo)準(zhǔn)線疊加的結(jié)果。2. 96g2. 74g2. 500g2. 494gL 148g2. 84g2. 487g2. 484g0. 99gR3. 556g289g如果地球模型不錯我們完全可以根據(jù)分層密度計算它的內(nèi)引力圖半徑R10.994820.8461780.5448120.2622820.194318R/N11.0052071.1817841.8354963.812695.146204密度

2、2.753.45.3510.712.312.5分層密度2.750.651.955.351.60.2疊加起來，就是模型引力圖了。如果模型不錯，那么地球的實際引把密度確定以后，我們再把密度分層，先按最小密度繪制大圓圖，再 把第二層密度減去第一層后繪圖，如此循環(huán)，畫到地核，再把所有值力就會和所畫圖一樣。還必須用至V內(nèi)部引力公式g(a)二警1 2a2R22 a4 a2 二3g=3g0 1 2a-2 2(a2)3(1- a2)2 +1)2a41 - a232 + 1Ja二sinRsin為質(zhì)點到球心的距離,g0為地表重力加速度AUTOcadvba模型程序代碼Sub 地球引力模型 ()Dim ml As O

3、bjectDim a As Variant, a1 As Variant, a2 As Variant, a3 As Variant, a4 As VariantDim v As Variant, v1 As Variant, v2 As Variant, v3 As Variant, v4 As Variant, v5 As VariantDim U As Variant, U1 As Variant, U2 As Variant, U3 As Variant, U4 As VariantDim co As Variantco = 3Dim p(0 To 20001) As DoubleDi

4、m i As IntegerDim J As Varianta = 1.00520697211556: 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: v

5、1 = v * a: v2 = v * a1: v3 = v * a2: v4 = v * a3: v5 = v * a4p(i) = iIf v1 1 And v1 a ThenU3 = 0End IfIf v2 1 And v2 a1 ThenU3 = 0End IfIf v3 1 And v3 a2 ThenU2 = 0End IfIf v4 1 And v4 a3 ThenU3 = 0End IfIf v5 1 And v5 a4 ThenU4 = 0End IfIf v = 1 Thenp(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 + b4 * a4 八-3) / v 八 2 End IfNext iSe