地球引力模型程序代碼

模型地球內(nèi)部引力圖地球模型計(jì)算的引力外地核處引力最大,外地核向內(nèi)或向外均在減小，從地核核心區(qū)處逃離地球引力需要能量為W=3.822MgR,由于不均勻分布，勢(shì)能增加了0.5MgR實(shí)際引力為分層球體的引力標(biāo)準(zhǔn)線疊加的結(jié)果。如果地球模型不錯(cuò)我們完全可以根據(jù)分層密度計(jì)算它的內(nèi)引力圖半徑R密度12.750.994823.40.8461785.350.54481210.70.26228212.30.19431812.5R/N分層密度12.751.0052070.651.1817841.951.8354965.353.812691.65.1462040.2把密度確定以后，我們?cè)侔衙芏确謱?，先按最小密度繪制大圓圖，再把第二層密度減去第一層后繪圖，如此循環(huán)，畫到地核，再把所有值疊加起來，就是模型引力圖了。如果模型不錯(cuò)，那么地球的實(shí)際引力就會(huì)和所畫圖一樣。當(dāng)然還必須要用到內(nèi)部引力公式AUTOcadvba模型程序代碼Sub 地球引力模型() Dim ml As Object Dim a As Variant, 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 Integer Dim J As Variant a = 1.00520697211556: a1 = 1.18178444724396: a2 = 1.8354955470878: a3 = 3.81269015792163: a4 = 5.14620364557066 b = 0.236363636363636: b1 = 0.709090909090909: 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 * 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 ml = T