SpaceEngineer, в скобках формулы (SemiMajorAxis = PericenterDist/(1.0 Eccentricity)) написан е=1.0, эта 1 для примера дана? Вместо неё вставляем е=N параметр так?
И ещё, не знал где спросить, спрошу тут: В Солнечной системе планеты вращвются вокруг Солнца. Но Солнце, физически не стоит на месте, оно ревгирует на движение каждой планеты. Я в Википедии недавно прочитал что Юпитер — единственная планета, у которой центр масс с Солнцем находится вне Солнца и отстоит от него примерно на 7 % солнечного радиуса.. Значит, у Солнца с каждой планетой есть Барицентр.
На такую мысль меня натолкнул следующий пример:
В двойной системе из Звёзд А и В с планетой С, которая совершает орбитальное движение вокруг звёзд А и В т.е. вокруг Барицентра (А+В), а это значит что между Барицентром (А+В) и планетой С существует Барицентр ((А+В)+С). Если тут так, значит и в Солнечной Системе тоже самое, только Барицентры планет земной группы находятся под поверхностью Солнца. Значит: У Солнца с Меркурием - Барицентр (С+М) (массы складываем); У Барицентра (С+М) с Венерой - Барицентр ((С+М)+В) ... и т.д. Из-за такой иерархической системы Барицентров, Солнце совершает сложные движения в пространстве... не так ли это?
Мои рассуждения верны? Если я неверно рассуждаю, объясните пожалуйста. Мне очень интересно... Спасибо.
P.S. - Если оффтоп, заранее извиняюсь... Windows 7 Ultimate x64 SP1 Intel Core i5 750 2.67 GHz 12 Gb 1333 MHz DDR3 RAM ASUS P7P55D PRO Sapphire HD7850 O.C. Edition 2 GB GDDR5 vRAM
SpaceEngineer, в скобках формулы (SemiMajorAxis = PericenterDist/(1.0 Eccentricity)) написан е=1.0, эта 1 для примера дана? Вместо неё вставляем е=N параметр так?
Там знак минуса должен быть: SemiMajorAxis = PericenterDist/(1.0 - Eccentricity)
Quote (NEMESIS)
В Солнечной системе планеты вращвются вокруг Солнца. Но Солнце, физически не стоит на месте, оно ревгирует на движение каждой планеты.
Да. Это естественно. Вспоминаем физику: центр масс системы (Солнечной) в отсутствие внешних сил движется линейно. В нашем случае влияет притяжение соседних звёзд и галактик, но этим можно пренебречь. Если рассматривать движение Солнца относительно центра масс системы, оно будет вот так болтаться. Раельно его не получится отписать иерархической системой барицентров, ведь все планеты влияют на Солнце одновременно, а плюс к этому влияют и друг на друга. Тут поможет только точное N-body моделирование.
Нет, я в какой-то толстой книге читал что солнце совершает незначительные сложные какие-то движения ...
Quote (SpaceEngineer)
Раельно его не получится отписать иерархической системой барицентров, ведь все планеты влияют на Солнце одновременно, а плюс к этому влияют и друг на друга.
SpaceEngineer, а если в SE свою систему описать иерархической системой Барицентров (чтоб как бы имитировать возмущения), движок поймёт? Хотелось бы посмотреть как будет болтаться звезда в ускоренном времени. Хотя при нажатии F2 будут видны одни барицентры... Но всётаки хотелось бы проэксперементировать... Windows 7 Ultimate x64 SP1 Intel Core i5 750 2.67 GHz 12 Gb 1333 MHz DDR3 RAM ASUS P7P55D PRO Sapphire HD7850 O.C. Edition 2 GB GDDR5 vRAM
Сообщение отредактировал NEMESIS - Четверг, 30.08.2012, 09:55
Хорошо, я попробую... только проблема с барицентром, который указан в каталоге звёзд. Он, получается, будет вращаться вокруг другого барицентра... В каталоге звёзд, могу ли я задать параметры орбиты этому Барицентру?
Параметры будут следующими:
StarBarycenter "Light Scream System/Nemesis System" { RA 17 00 05 // Right Ascention in hours, in decimal format or formatted as HH MM SS Dec 36 06 01 // Declination in degrees, in decimal format or formatted as DD MM SS Dist 1.66 // Distance from our Sun in parsecs
Orbit { RefPlane "Ecliptic" //Orbit reference plane Epoch 2451545.0 //Epoch, Ellian's epoch Period 0.36687512 //Orbital Period of Body, in Years Eccentricity 0.00351123 //Eccentricity of Body orbit (if e=0 is cyrcle, if 1>e>0 is ellipse, if e=1 is parabola, if e>1 is hyperbola) Inclination 9.49 //Incilaption of Orbit in degrees, 0<i<90 - is direct, 90<i<180 - is inverse AscendingNode -166.0 //Ascending Node in Degrees
//Next two parameters are interchangeable: SemiMajorAxis = PericenterDist/(1.0-Eccentricity) SemiMajorAxis 0.00000365 //Mass Ratio 583887.82289994:3.448123746 in AU PericenterDist 0.00000364 //Pericentric distanse in AU
//Next two parameters are interchangeable: ArgOfPericen = LongOfPericen - AscendingNode ArgOfPericen 289.2 //Argument of Pericenter, in Degrees. The ArgOfPericenter parameter of Body A must differ by 180 degrees from the ArgOfPericenter parameter of another Body B LongOfPericen 123.2 //Longitude of Pericenter, in Degrees
//Next two parameters are interchangeable: MeanAnomaly = MeanLongitude - LongOfPericen MeanAnomaly 107.0 //Meam anomaly, in Degrees MeanLongitude 230.2 //Mean longitude, in Degrees
AscNodePreces 101034 //Ascending node precession period, in Years ArgOfPeriPreces 123602 //Argument of pericenter precession period, in Years } }
Верно указал? Или всё таки надо описывать этот же Барицентр в каталоге планет?
Добавлено (30.08.2012, 14:12) --------------------------------------------- Я чуть не запутался... Нафиг тот скрипт... В новом я закомментировал другие планеты кроме одной, массивной... В Каталоге Звёзд:
Code
StarBarycenter "(Nemesis A-Nemesis B)-Goowlen" { RA 17 00 05 // Right Ascention in hours, in decimal format or formatted as HH MM SS Dec 36 06 01 // Declination in degrees, in decimal format or formatted as DD MM SS Dist 1.66 // Distance from our Sun in parsecs }
В каталоге Планет для Двойной звезды (Звезда-Звезда):
Code
Star "Nemesis A/Light Scream A" { ParentBody "Nemesis A-Nemesis B" //The Body, around wich this star orbits Class "F2V" //Spectral Class of a Star Lum 6.33 //Luminosity in Sun Luminosity's
Mass 1.71325519 //Mass of a Star, in Sun Masses Radius 1183155.125 //Radus 1/2 of Diameter, in Km Obliquity 0.00 //Obliquity of rotation axis, in Degrees //AppMagn 00.00 //Visual Star Magnitude RotationPeriod 705.01 //Rotation Period around axis in Hours RotationOffset 64.7 //In Degrees //AbsMagn 00.00 //Absolute Star Magnitude
Orbit { Period 0.106779 //Orbital Period of Body, in Years SemiMajorAxis 0.00767278 //Mass Ratio 1.71:0.40 in AU Eccentricity 0.39300467 //Eccentricity of Body orbit (if e=0 is cyrcle, if 1>e>0 is ellipse, if e=1 is parabola, if e>1 is hyperbola) Inclination 23.3 //Incilaption of Orbit in degrees, 0<i<90 - is direct, 90<i<180 - is inverse AscendingNode 16.8 //Ascending Node in Degrees MeanAnomaly -180.0 //Mean Anomaly in Degrees ArgOfPericen 90.0 //The ArgOfPericenter parameter of Body A must differ by 180 degrees from the ArgOfPericenter parameter of another Body B } }
Star "Nemesis B/Light Scream B" { ParentBody "Nemesis A-Nemesis B" ///The Body, around wich this star orbits Class "M4V" //Spectral Class of a Star Lum 0.04 //Luminosity in Sun Luminosity's
Mass 0.04025648 //Mass of a Star, in Sun Masses Radius 265430.5 //Radus 1/2 of Diameter, in Km Obliquity 0.00 //Obliquity of rotation axis, in Degrees //AppMagn 00.00 //Visual Star Magnitude RotationPeriod 111.39 //Rotation Period around axis in Hours RotationOffset 36.7 //In Degrees //AbsMagn 00.00 //Absolute Star Magnitude
Orbit { Period 0.106779 //Orbital Period of Body, in Years SemiMajorAxis 0.31880892 //Mass Ratio 1.71:0.40 in AU Eccentricity 0.39300467 //Eccentricity of Body orbit (if e=0 is cyrcle, if 1>e>0 is ellipse, if e=1 is parabola, if e>1 is hyperbola) Inclination 23.3 //Incilaption of Orbit in degrees, 0<i<90 - is direct, 90<i<180 - is inverse AscendingNode 16.8 //Ascending Node in Degrees MeanAnomaly -180.0 //Mean Anomaly in Degrees ArgOfPericen 270.0 //The ArgOfPericenter parameter of Body A must differ by 180 degrees from the ArgOfPericenter parameter of another Body B } }
В каталоге Планет для Барицентра ((Звезда+Звезда)+Планета):
Code
StarBarycenter "Nemesis A-Nemesis B" { ParentBody "(Nemesis A-Nemesis B)-Goowlen" //The Body, around wich this planet orbits
Orbit { RefPlane "Ecliptic" //Orbit reference plane Epoch 2451545.0 //Epoch, Ellian's epoch Period 6.29650347 //Orbital Period of Body, in Years Eccentricity 0.50679003 //Eccentricity of Body orbit (if e=0 is cyrcle, if 1>e>0 is ellipse, if e=1 is parabola, if e>1 is hyperbola) Inclination 23.0 //Incilaption of Orbit in degrees, 0<i<90 - is direct, 90<i<180 - is inverse AscendingNode 27.5 //Ascending Node in Degrees
//Next two parameters are interchangeable: SemiMajorAxis = PericenterDist/(1.0-Eccentricity) SemiMajorAxis 0.00329587 //Mass Ratio 583887.82289994:473.07032588 in AU PericenterDist 0.00000364 //Pericentric distanse in AU
//Next two parameters are interchangeable: ArgOfPericen = LongOfPericen - AscendingNode ArgOfPericen 90.0 //Argument of Pericenter, in Degrees. The ArgOfPericenter parameter of Body A must differ by 180 degrees from the ArgOfPericenter parameter of another Body B LongOfPericen 117.5 //Longitude of Pericenter, in Degrees
//Next two parameters are interchangeable: MeanAnomaly = MeanLongitude - LongOfPericen MeanAnomaly -19.1 //Meam anomaly, in Degrees MeanLongitude 98.4 //Mean longitude, in Degrees
AscNodePreces 129164 //Ascending node precession period, in Years ArgOfPeriPreces 260001 //Argument of pericenter precession period, in Years } }
В Каталоге Планет для Планеты:
Code
Planet "Goowlen" { ParentBody "(Nemesis A-Nemesis B)-Goowlen" //The Body, around wich this planet orbits Class "GasGiant" //Class of a Body
Radius 79099.34 //Radius, in Km Mass 473.07032588 //Mass of Body in Earth's masses Oblateness 0.0 //Oblateness of the Planet
Albedo 0.306 //Reflection Coefficient Color (1.000, 0.714, 0.940) //Color of particles
RotationPeriod 18.04277 //Rotation Period in Hours RotationOffset 33.421 //Rotation phase in Degrees RotationEpoch 2451545.0 //Epoch, Ellian's epoch Obliquity -4.396001 //Obliquity of rotation axis, in Degrees EqAscendNode 38.67 //Longitude of Ascending Node of rotation in Degrees Precession 22195.3 //Precession period, in Years
Surface { //These parameters describes a texture pack stored on the disk. If the planet uses a procedural generation of surface, comment out these parameters... //DiffMap "Cyclone B/Surface" //Color texture //DiffTileSize 512 //Color texture tile resolution //DiffTileBorder 1 //Color texture tile border width //BumpMap "Cyclone B/Bump" //Bump (heightmap) texture //BumpTileSize 512 //Bump texture tile resolution //BumpTileBorder 1 //Bump texture tile border width
BumpHeight 10.22 //Max Height of landscape, in Km BumpOffset 0.219 //Offset of Zero Height level, in Km DiffMapAlpha "Water" //Mode of specular effect (flecks of sunlight) SpecularBright 20.0 //Brightness of specular spot SpecularPower 150.0 //Size of specular spot DayAmbient 3.0 //Brightness of fake daytime ambient light Lommel 0.0 //Lambert to Lommel-Seeliger lighting mode ratio Exposure 1.99 //Global Surface Brightness Life false //false, if there is life on the planet - true
//These parameters for procedural surface generation... Style 0.04579 //Default color scheme Randomise (0.566, 0.490, -0.701) //Randomize vector colorDistMagn 1.6603 //Magnitude of detail textures distortion colorDistFreq 0.182 //Frequency of detail textures distortion //detailScale 0.0 //Scale of detail textures //colorConversion true //Perform textures color conversion //seaLevel 0.0 //Snow line level (0...1) //snowLevel 0.0 //Snow line level (0...1) //tropicWidth 0.0 //Width of tropics //tropicLatitude 0.0 //Latitude of tropics //icecapLatitude 0.0 //Latitude of ice caps boundaries //icecapHeight 0.0 // Height of ice caps //climatePole 0.0 //Climate of poles - index to color table (0...1) //climateTropic 0.0 //Climate of tropics - index to color table (0...1) //climateEquator 0.0 //Climate of equator - index to color table (0...1) mainFreq 0.0308 //Main noise frequency (oceans/continents) mainOctaves 8 //Number of octaves of Main Noise //venusFreq 0.0 //Venus-like global structures frequency //venusMagn 0.0 //Venus-like global structures magnitude //mareFreq 0.0 //Lunar mare frequency //mareDensity 0.0 //Lunar mare density //montesMagn 0.0 //Mountains magnitude (0...1) //montesFreq 0.0 //Mountains frequency //montesDensity 0.0 //Mountains density //dunesMagn 0.0 //Dunes magnitude (0...1) //dunesFreq 0.0 //Dunes frequency //dunesDensity 0.0 //Dunes density //hillsMagn 0.0 //Hills magnitude (0...1) //hillsFreq 0.0 //Hills frequency //hillsDensity 0.0 //Hills density //canyonMagn 0.0 //Canyons magnitude (0...1) //canyonFreq 0.0 //Canyons frequency //canyonDensity 0.0 //Canyons density //cracksMagn 0.0 //Ice cracks magnitude (0...1) //cracksFreq 0.0 //Ice cracks frequency //cracksOctaves 0 //Number of ice cracks octaves //craterMagn 0.0 //Craters magnitude (0...1) //craterFreq 0.0 //Craters frequency //craterDensity 0.0 //Craters density //craterOctaves 0.0 //Number of craters octaves //craterRayedFactor 0.0 //Number of rayed craters relative to usual craters cycloneMagn 14.165 //Cyclones magnitude (0...1) cycloneFreq 0.2793 //Cyclones frequency cycloneDensity 0.1037 //Cyclones density cycloneOctaves 1 //Number of cyckone octaves twistZones 2.119 //Number of Jupiter-like zones or strips twistMagn 0.617 //Strength of twist colorSea (0.630, 0.600, 0.550, 0.009) //Color of seas colorShelf (0.703, 0.706, 0.691, 0.010) //Color of shelf areas colorBeach (0.803, 0.870, 0.920, 0.000) //Color of beach lines colorDesert (0.650, 0.420, 0.270, 0.000) //Color of deserts colorLowland (0.730, 0.330, 0.270, 0.000) //Color of lower lands colorUpland (0.750, 0.346, 0.290, 0.000) //Color of upper lands colorRock (0.900, 0.550, 0.290, 0.000) //Color of rocks colorSnow (0.950, 0.330, 0.220, 0.016) //Color of snow colorLowPlants (0.730, 0.330, 0.270, 0.000) //Color of lower vegetation colorUpPlants (0.740, 0.330, 0.270, 0.000) //Color of upper vegetation }
Clouds { } Clouds { } Clouds { }
Atmosphere { Model "Jupiter" //Model of Atmosphere Height 770.12 //Height of the top boundary, in Km //Pressure 0.0 //Surface pressure, in Bars //Density 0.0 //Surface Dencity in Kg/m^3 //Greenhouse 0.0 //Greenhouse effect in Degrees Bright 10.0 //Brightness (rendering parameter) Opacity 1.0 //Opacity (rendering parameter) SkyLight 0.99 //Brightness of Skylight on terrain EclipseBright 250.0 //Brightness of solar eclipse shadow EclipseColor (1.000 0.679 0.204) //Color of solar eclipse shadow }
Rings { //These parameters describes a texture pack stored on the disk. If the planet uses a procedural generation of textures, comment out these parameters... InnerRadius 121347.07 //Inner Radius of Rings, in Km OuterRadius 222031.13 //Outer Radius of Rings, in Km RotationPeriod 44.33 //Rings Rotation Period, in Hours RotationOffset 24.93 //Rings Rotation phase, in Degrees FrontBright 1.0 //Direct lighting brightness BackBright 5.0 //Back lighting brightness (scattering) Density 1.0 //Density (Opacity) Exposure 2.0 //Global Brightness of Rings }
Orbit { RefPlane "Ecliptic" //Orbit reference plane Epoch 2451545.0 //Epoch, Ellian's epoch Period 6.29650347 // Orbital Period, in Years Eccentricity 0.50679003 //Eccentricity of Body orbit (if e=0 is cyrcle, if 1>e>0 is ellipse, if e=1 is parabola, if e>1 is hyperbola) Inclination 23.0 //Inclination of Orbit in degrees, 0<i<90 - is direct, 90<i<180 - is inverse AscendingNode 27.5 //Longitude of Ascending Node in Degrees
//Next two parameters are interchangeable: SemiMajorAxis = PericenterDist/(1.0 Eccentricity) SemiMajorAxis 4.06793137 //Mass Ratio 583887.82289994:473.07032588 in AU PericenterDist 2.00634431 //Pericentric distanse in AU
//Next two parameters are interchangeable: ArgOfPericen = LongOfPericen - AscendingNode ArgOfPericen 270.0 //Argument of Pericenter, in Degrees. The ArgOfPericenter parameter of Body A must differ by 180 degrees from the ArgOfPericenter parameter of another Body B LongOfPericen 297.5 //Longitude of Pericenter, in Degrees
//Next two parameters are interchangeable: MeanAnomaly = MeanLongitude - LongOfPericen MeanAnomaly -19.1 //Meam anomaly, in Degrees MeanLongitude 278.4 //Mean longitude, in Degrees
AscNodePreces 129164 //Ascending node precession period, in Years ArgOfPeriPreces 260001 //Argument of pericenter precession period, in Years } }
Дома проверю, как получится... Чем больше планет тем сложнее с Барицентрами париться... Получается, самый последнй Барицентр и будет в Каталоге Звёзд, потомучто последняя планета Газовый Гигант... Да? Windows 7 Ultimate x64 SP1 Intel Core i5 750 2.67 GHz 12 Gb 1333 MHz DDR3 RAM ASUS P7P55D PRO Sapphire HD7850 O.C. Edition 2 GB GDDR5 vRAM
Сообщение отредактировал NEMESIS - Четверг, 30.08.2012, 14:22
Именно используя это явление и открывают экзопланеты методом лучевых скоростей.
Да я это знаю но в тех случаях находятся планеты с большой массой и малым а.
Quote (SpaceEngineer)
Порядка миллиона км, а скорости - сотни м/с.
Я так понимаю для солнца это незаметное возмущение AMD Phenom II X2 555 3200 MHz; RAM 4GB 1333 Hz; ATI Radeon HD 5670 1024 Mb; ОС Win 7 tabulorasa edition 32 bit
У меня возникла проблемка с этой (T=2πR/v) формулой...
Короче...
Для начала, хотел расчитать орбитальную скорость между телами вручную (я знаю что SE автоматически вычисляет, но хочется научиться самому это делать). Решил поэксперементировать на двойной планете: 1-я = 473.07032588 массы Земли; 2-я = 24.06881274 массы Земли.
Общее расстояние между телами 4,620,541.17 Км. Расстояние 1-ой до Барицентра (по формуле вычислил) = 223,701.18 Км; Расстояние 2-ой до Барицентра (по формуле вычислил) = 4,396,839.99 Км.
Тут всё впорядке...
Исходя из вышеперечисленного, у меня не известна орбитальная скорость, и орбитальный период. Ор. период, ради примера, я взял как 36.498 дней... Неизвестным остаётся орбиталная скорость. Значит: v=2πR/T, верно?
Получается, v(1-ой Планеты)=0.094 км/с, а v(2-ой Планеты)=1.848 км/с. Скорость вообще какая-то нереальная получилась, не так ли?
Чтоб проверить это, для примера взял Нептун. В Википедии всё есть...
Большая Полуось = 30.10366151 а.е. Орбитальный Период = 164.79 лет В Вики Орбитальная Скорость Нептуна =5.43 км/с а у меня пучему-то 1.148 км/с получается. Где я ошибся? И вообще, правильно ли я вычисляю? Чувствую я где-то лохонулся...
Так нептун обращается вокруг солнца на огромном растоянии от него. А твои планеты вращаются друг вокруг друга. массы планет меньше солнца, значит и скорость будет меньше, ведь от массы зависит большая полуось. А от скорость их барицентра вокруг их солнца - другое уже дело AMD Phenom II X2 555 3200 MHz; RAM 4GB 1333 Hz; ATI Radeon HD 5670 1024 Mb; ОС Win 7 tabulorasa edition 32 bit
Если масса m не пренебрежимо мала по сравнению с M, то планета будет обращаться вокруг Солнца с той же скоростью и по той же орбите, что и материальная точка, обращающаяся вокруг массы M+m. При этом массу M в последней формуле нужно заменить на M+m:
Вообще, для чего?
M - это масса звезды, в массах Земли или Солнца? m - это масса планеты в массах Земли, как я понял; T - Период в Годах, да? G - это что? нельзя ли заменить его на 4π2? (заглянул Сюда) a - этр большая полуось...
Какую из них я могу использовать для системы Звезда-Планета а какую для Звезда-звезда? Объясните пожалуйста, очень хочется разобраться...
это усовершенствованый (или как его там называют) 3-тий закон кеплера. отношение квадратов периодов 2-х планет равно отношению кубов больших полуосей этих же планет. тоесть T/a=const.
Какую из них я могу использовать для системы Звезда-Планета а какую для Звезда-звезда?
обе разницы между звездой и планетой нет. они обе имеют массу и подчиняются законам физики AMD Phenom II X2 555 3200 MHz; RAM 4GB 1333 Hz; ATI Radeon HD 5670 1024 Mb; ОС Win 7 tabulorasa edition 32 bit
Две звезды: m1=1.21082119, sma=0.05915793 а.е. по формуле (Тело A: SemiMajorAxis = R * MassB / (MassA + MassB)) m2=0.74156327, sma=0.09659281 а.е. по формуле (Тело B: SemiMajorAxis = R * MassA / (MassA + MassB))
По формуле я заменяю буквы на соответствующие значения и фигня полная получается Судя по формуле мне необходимо найти квадрат периода? Мне период нужен а не квадратный период.... Переводя квадрат. пер. в пер. ответ получается очень мелким! Вычисляя формулу я сразу вывожу куб полуоси....
Я так вычислял:
Кв. Период = 4*(2*3.14159265)) / 6.67*(1.21082119+0.74156327) и это всё умножить на 0.09659281(m2) в кубе....
Вместо М я вставил массу m1, так как это главное тело системы...
У меня мозги взорвутся скоро... покажите хоть пример, а? Я хоть правильно делаю?
NEMESIS, скинь все даные про сис-тему и я посчитаю да и скрин тоже, что-бы я себе мог представить движение тел AMD Phenom II X2 555 3200 MHz; RAM 4GB 1333 Hz; ATI Radeon HD 5670 1024 Mb; ОС Win 7 tabulorasa edition 32 bit
Сообщение отредактировал GV_FiQst - Воскресенье, 09.09.2012, 20:36