《《地球引力模型程序代码》》由会员分享，可在线阅读，更多相关《《地球引力模型程序代码》（6页珍藏版）》请在金锄头文库上搜索。
1、L 148g模型地球部引力图地球模型计算的引力外地核处引力最大,外地核向或向外均在减 小，从地核核心区处逃离地球引力需要能量为 W=3.822MgR,由于 不均匀分布，势能增加了 0.5MgR 实际引力为分层球体的引力标准线叠加的结果。4.472g 12. 96g2 74g2. 500g2. 494g如果地球模型不错我们完全可以根据分层密度计算它的引力图半径R密度12.750.994823.40.8461785.350.54481210.70.26228212.30.19431812.5R/N分层密度12.751.0052070.651.1817841.951.8354965.353.8126
2、91.65.1462040.2把密度确定以后，我们再把密度分层，先按最小密度绘制大圆图，再把第二层密度减去第一层后绘图，如此循环，画到地核，再把所有值叠加起来，就是模型引力图了。如果模型不错，那么地球的实际引力就会和所画图一样。当然还2a4a33(1 a2)2 11 a22a4a33(1a2)21a2a sin , Rsin为质点到球心的距离为地表里力加速度AUTOcadvba模型程序代码Sub地球引力模型()Dim ml As ObjectDim a As Variant, al As Variant, a2 As Variant, a3 As Variant, a4 As VariantD
3、im 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 DoubleDim i As IntegerDim J As Varianta=1.00520697211556:al=1.18178444724396: a2=1.8354955
4、470878: 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 And v1 a ThenU3 = 0End
5、 IfIf v2 1 And v2 a1 ThenU3 = 0End IfIf v3 1 And v3 = a2 ThenU2 = 500 * b2 * a2 a -3 / v a 2U2 = 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 a 2 + (v a 4) / 3) / (1 + (1 - v a 2) a1.5) + U + U1 + U2 + U3 + U4Elsep(i + 1) = 500 * (1 + b * a A -3 + bl * al A -3 + b2 * a2 A -3 + b3 *a3 a -3 + b4 * a4 a -3) / v a 2End IfNext iSet ml = ThisDrawing.ModelSpace.AddLightWeightPolyline(p)co = co + 1ml.color = coEnd Sub
