地球引力模型程序代码讲稿.docVIP

下载本文档

8
0
约2.71千字
约 6页
2017-09-03 发布于湖北
举报
版权申诉

地球引力模型程序代码讲稿.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

模型地球内部引力图地球模型计算的引力外地核处引力最大,外地核向内或向外均在减小，从地核核心区处逃离地球引力需要能量为W=3.822MgR,由于不均匀分布，势能增加了0.5MgR 实际引力为分层球体的引力标准线叠加的结果。如果地球模型不错我们完全可以根据分层密度计算它的内引力图半径R 密度 1 2.75 0.99482 3.4 0.846178 5.35 0.544812 10.7 0.262282 12.3 0.194318 12.5 R/N 分层密度 1 2.75 1.005207 0.65 1.181784 1.95 1.835496 5.35 3.81269 1.6 5.146204 0.2 把密度确定以后，我们再把密度分层，先按最小密度绘制大圆图，再把第二层密度减去第一层后绘图，如此循环，画到地核，再把所有值叠加起来，就是模型引力图了。。如果模型不错，那么地球的实际引力就会和所画图一样。当然还必须要用到内部引力公式 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 Then U = 750 * b * (1 + v1 * 2 - (1 - v1 ^ 2 + (v1 ^ 4) / 3) / (1 + (1 - v1 ^ 2) ^ 1.5)) / a ElseIf v1 1 And v1 = a Then U = 500 * b * a ^ -3 / v ^ 2 ElseIf v1 a Then U3 = 0 End If If v2 = 1 Then U1 = 750 * b1 * (1 + v2 * 2 - (1 - v2 ^ 2 + (v2 ^ 4) / 3) / (1 + (1 - v2 ^ 2) ^ 1.5)) / a1 ElseIf v2 1 And v2 = a1 Then U1 = 500 * b1 * a1 ^ -3 / v ^ 2 ElseIf v2 a1 Then U3 = 0 End If If v3 = 1 Then U2 = 750 * b2 * (1 + v3 * 2 - (1 - v3 ^ 2 + (v3 ^ 4) / 3) / (1 + (1 - v3 ^ 2) ^ 1.5)) / a2 ElseIf v3 1 And v3 = a2 Then U2 = 500 * b2 * a2 ^ -3 / v ^ 2 ElseIf v3 a2 Then U2 = 0 End If If v4 = 1 Then U3 = 750 * b3 * (1 + v4 * 2 - (1 - v4 ^ 2 + (v4 ^ 4) / 3) / (1 + (1 - v4 ^ 2) ^ 1.5)) / a3 ElseIf v4 1 And