Re: astrodinamica

From: ernesto <ernesto.alto_at_libero.it>
Date: Wed, 12 Sep 2001 22:52:12 GMT

Nel caso che l'ipotesi fatta da alcuni astronomi in base alla
provenienza delle comete a lunga orbita sia vera, ovvero che la
presunta nube di Oort sia perturbata dalla presenza di un pianeta di
10 o pi� masse gioviane orbitante intorno al Sole a una distazna
compresa fra 1 e 2 anni luce, potrebbero essersi formati dei punti di
Lagrange equigravitazionali e quindi in quei punti essersi
"raggruppati" pi� "corpi" giustificando quelle che sembrano essere
provenienze privilegiate nel cielo?

Io ho buttato gi� un programmino in VB per calcolare i punti di
Lagrange tra due corpi di massa con rapporto tra M/m < 0,03 e quindi
con centro di massa nel corpo M.
L'ho usato con Sole-Giove e d� risultati accettabili anche se non ho
iterato le equazioni per raffinarli.
Lo si pu� usare anche per trovare eventuali punti di equilibrio
gravitazionale tra il supposto "corpo" orbitante fra 1 e 2 anni luce
dal Sole?

Forse gli eventuali punti di Lagrange del sistema Sole-Corpo X
avrebbero potuto raggrumare in qualche miliardo di anni una parte
delle comete?
Qualcuno � in grado di controllare e migliorare il programmino che
segue?
Grazie per l'attenzione.
Ernesto
-----------------------------------------------------
PROGRAMMA IN VB PER IL CALCOOLO DEI PUNTI DI LAGRANGE
-----------------------------------------------------
Option Explicit
Dim MassaCorpoGrande As Integer
Dim d As Integer
Dim kmL1 As Long
Dim kmL2 As Long
Dim kmL3 As Long
Dim kmL4 As Long
Dim kmL5 As Long
Dim kmL2daCorpoGrande As Long
Dim kmL1daCorpoGrande As Long
Dim kmL3daCorpoGrande As Long
Dim kmL3daCorpoPiccolo As Long
Dim kmL2daCorpoPiccolo As Long
Dim kmL1daCorpoPiccolo As Long
Dim a As Double
Dim b As Double
Dim m As Double
Dim x2 As Double
Dim x1 As Double
Dim x3 As Double
Dim L1 As Double
Dim L2 As Double
Dim L3 As Double
Dim L4 As Double
Dim L5 As Double

Private Sub Command1_Click()
MassaCorpoGrande = 1
m = InputBox("immetti la massa del corpo piccolo in millesimi rispetto
a quella del corpo grande", "Lagrange points")
m = m / 1000

d = InputBox("immetti la distanza media tra i due corpi in milioni di
km.", "Lagrange Points")

a = 1 / (1 + m) 'distanza di M dal centro di massa
b = m / (1 + m) 'distanza di m dal centro di massa
                  'si fa uguale a 1 la somma a+b

x2 = (b / 3 * (1 - b)) ^ (1 / 3) 'primo valore stimato valore per
                           ' eventuale iterazione
x1 = (b / 3 * (1 - b)) ^ (1 / 3) 'idem
x3 = 1 - (7 * b / 12) 'valore per eventuale iterazione per
                            aumentare la precisione

L2 = (b * (1 - x2) ^ 2) / (3 - 2 * b - x2 * (3 - b - x2))
L2 = Log(L2) * 1 / 3
L2 = Exp(L2)
L2 = Round(L2, 5)
kmL2 = d * L2 'trasformazione in km
kmL2daCorpoGrande = d - kmL2
kmL2daCorpoPiccolo = kmL2

L1 = (b * (1 + x1) ^ 2) / (3 - 2 * b + x1 * (3 - b + x1))
L1 = Log(L1) * 1 / 3
L1 = Exp(L1)
L1 = Round(L1, 5)
kmL1 = d * L1 'trasformazione in KM
kmL1daCorpoGrande = d + kmL1
kmL1daCorpoPiccolo = kmL1

L3 = ((1 - b) * (1 + x3) ^ 2) / (1 + 2 * b + x3 * (2 + b + x3))
L3 = Log(L3) * 1 / 3
L3 = Exp(L3)
L3 = Round(L3, 5)
kmL3 = d * L3 'trasformazione in KM
kmL3daCorpoGrande = kmL3
kmL3daCorpoPiccolo = 2 * d

kmL4 = d
kmL5 = d

lb1Risposta.Caption = " il valore di L2 � " & L2 & ", cio� di " &
kmL2daCorpoGrande & "milioni di km dal corpo M e di km " &
kmL2daCorpoPiccolo & " dal corpo m, sulla congiungente dei due corpi."
lb2Risposta.Caption = " il valore di L1 � " & L1 & ", cio� di " &
kmL1daCorpoGrande & "milioni di km dal corpo M e di km " &
kmL1daCorpoPiccolo & " oltre l'orbita del corpo m, sulla congiungente
dei due corpi"
lb3Risposta.Caption = " il valore di L3 � " & L3 & ", cio� di " &
kmL3daCorpoGrande & "milioni di km dal corpo M ma sulla parte opposta
dell'orbita di del corpo m, e di km " & kmL3daCorpoPiccolo & " dal
corpo m, sulla congiungente dei due corpi."
lb4Risposta.Caption = " il valore di L4 precede di 60� sull'orbita di
m a circa km " & kmL4
lb5Risposta.Caption = " il valore di L5 segue di 60� sull'orbita di m
a circa km " & kmL5

End Sub
-------------------------------------------------------------------

Il terrorismo spera si risponda alla barbarie con altra barbarie:� una trappola!
Received on Thu Sep 13 2001 - 00:52:12 CEST