地球引力模型程序代码

模型地球内部引力图地球模型计算的引力外地核处引力最大,外地核向内或向外均在减小，从地核核心区处逃离地球引力需要能量为W=3.822MgR,由于不均匀分布，势能增加了0.5MgR实际引力为分层球体的引力标准线叠加的结果。如果地球模型不错我们完全可以根据分层密度计算它的内引力图半径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把密度确定以后，我们再把密度分层，先按最小密度绘制大圆图，再把第二层密度减去第一层后绘图，如此循环，画到地核，再把所有值叠加起来，就是模型引力图了。。如果模型不错，那么地球的实际引力就会和所画图一样。当然还必须要用到内部引力公式4422033222220113g333g()12=1222(1)1(1)1sin,singaaaaMGaaaRaaaR为质点到球心的距离,为地表重力加速度AUTOcadvba模型程序代码Sub地球引力模型()DimmlAsObjectDimaAsVariant,a1AsVariant,a2AsVariant,a3AsVariant,a4AsVariantDimvAsVariant,v1AsVariant,v2AsVariant,v3AsVariant,v4AsVariant,v5AsVariantDimUAsVariant,U1AsVariant,U2AsVariant,U3AsVariant,U4AsVariantDimcoAsVariantco=3Dimp(0To20001)AsDoubleDimiAsIntegerDimJAsVarianta=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-02Fori=0To10000Step2v=i/1000:v1=v*a:v2=v*a1:v3=v*a2:v4=v*a3:v5=v*a4p(i)=iIfv1=1ThenU=750*b*(1+v1*2-(1-v1^2+(v1^4)/3)/(1+(1-v1^2)^1.5))/aElseIfv11Andv1=aThenU=500*b*a^-3/v^2ElseIfv1aThenU3=0EndIfIfv2=1ThenU1=750*b1*(1+v2*2-(1-v2^2+(v2^4)/3)/(1+(1-v2^2)^1.5))/a1ElseIfv21Andv2=a1ThenU1=500*b1*a1^-3/v^2ElseIfv2a1ThenU3=0EndIfIfv3=1ThenU2=750*b2*(1+v3*2-(1-v3^2+(v3^4)/3)/(1+(1-v3^2)^1.5))/a2ElseIfv31Andv3=a2ThenU2=500*b2*a2^-3/v^2ElseIfv3a2ThenU2=0EndIfIfv4=1ThenU3=750*b3*(1+v4*2-(1-v4^2+(v4^4)/3)/(1+(1-v4^2)^1.5))/a3ElseIfv41Andv4=a3ThenU3=500*b3*a3^-3/v^2ElseIfv4a3ThenU3=0EndIfIfv5=1ThenU4=750*b4*(1+v5*2-(1-v5^2+(v5^4)/3)/(1+(1-v5^2)^1.5))/a4ElseIfv51Andv5=a4ThenU4=500*b4*a4^-3/v^2ElseIfv5a4ThenU4=0EndIfIfv=1Thenp(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+b1*a1^-3+b2*a2^-3+b3*a3^-3+b4*a4^-3)/v^2EndIfNextiSetml=ThisDrawing.ModelSpace.AddLightWeightPolyline(p)co=co+1ml.color=coEndSub