Моделирование солнечной системы - где мои вычисления ошибочны? - PullRequest
Новая
       24

Моделирование солнечной системы - где мои вычисления ошибочны?

0 голосов
/ 16 мая 2019

Я делаю 2D-моделирование солнечной системы для проекта с использованием физических уравнений, координат x и y и словаря с использованием 2D-списков. Я использую холст tkinter для создания анимации.

Планета Земля, кажется, исчезла с экрана сначала с очень медленным ускорением, а затем с огромным ускорением вскоре после этого. Я не вижу проблемы.

Ранее у меня была только земля, вращающаяся вокруг Солнца, которая была успешной, используя уравнения аналогичным образом и перемещая землю на холсте с помощью компонентов x и y вычисленного смещения. Словарь по порядку {(-body):[[x, y], [x velocity, y velocity], mass, [change in x displacement, change in y displacement]}. Рассчитанные значения сохраняются или добавляются к некоторым значениям в этом словаре при необходимости. Вот код, который я использую: 

G = 6.67384 * 10 ** -11                                              
scale = 10 ** 13                                                                             
speed = 1                 
global user_status
screen = Tk()
screen.title('Solar System' + ' - ' + user_status)
screen.geometry('1300x700')

ani = Canvas(screen, width=1300, height=700)
ani.pack()
ani.create_rectangle(0, 0, 1300, 700, fill='Black')

sunx = 636
suny = 343
sun = ani.create_oval(sunx-10, suny-10, sunx+10, suny+10, fill='yellow')

earthx = 746
earthy = 343

moonx = 747
moony = 343

bodies = {'sun':  [[sunx, suny],     [0, 0],                         1.989 * 10 ** 30 * speed / scale,      [0, 0]],
          'earth':[[earthx, earthy], [0, 347.3833062 * 1.033447099], 5.972 * 10 ** 24 * speed / scale,      [0, 0]],
          'moon': [[moonx, moony],   [0, 360],                       7.34767309 * 10 ** 22 * speed / scale, [0, 0]]
          }

body_names = []
for Ω in bodies.keys():
    body_names.append(Ω)

moon = ani.create_oval(moonx - 4, moony - 4, moonx + 4, moony + 4, fill='grey70')
earth = ani.create_oval(earthx-6, earthy-6, earthx+6, earthy+6, fill='sky blue')

timestep = 0.0001

while True:
    for i in range(len(body_names)):
        body1 = body_names[i]
        x1 = bodies[body1][0][0]
        y1 = bodies[body1][0][1]
        total_Fx = 0
        total_Fy = 0
        body1_mass = bodies[body1][2]
        for j in range(len(body_names)):
            body2 = body_names[j]
            if body1 != body2:
                x2 = bodies[body2][0][0]
                y2 = bodies[body2][0][1]
                body2_mass = bodies[body2][2]

                r = sqrt(((x1 - x2) ** 2) + ((y1 - y2) ** 2))
                rx = (x1 - x2)
                angle = (acos(rx/r))

                F = (G * body1_mass * body2_mass) / (r ** 2)
                Fx = F * cos(angle)
                Fy = F * sin(angle)

                total_Fx += Fx
                total_Fy += Fy

        ax = (total_Fx / body1_mass)
        ay = (total_Fy / body1_mass)

        ux = bodies[body1][1][0]
        uy = bodies[body1][1][1]

        vx = ux - (ax * timestep)
        if x1 <= sunx:
            vy = uy + (ay * timestep)
        else:
            vy = uy - (ay * timestep)

        sx = vx * timestep * speed
        sy = vy * timestep * speed

        bodies[body1][3][0] = sx
        bodies[body1][3][1] = sy

        bodies[body1][1][0] += vx
        bodies[body1][1][1] += vy

        bodies[body1][0][0] = x1 + sx
        bodies[body1][0][1] = y1 + sy
        print(bodies[body1][1], body1)

    move_e_x = bodies['earth'][3][0]
    move_e_y = bodies['earth'][3][1]
    ani.move(earth, move_e_x, move_e_y)

    move_m_x = bodies['moon'][3][0]
    move_m_y = bodies['moon'][3][1]
    ani.move(moon, move_m_x, move_m_y)

    ani.update()

Я хочу, чтобы объекты одновременно вращались вокруг Солнца и друг друга, но я получаю это сообщение об ошибке после того, как тела далеко ушли от экрана: 

Traceback (most recent call last):
  File "/Users/apple/Documents/School/Computer Science/NEA/NEA Programming/Solar System Simulator.py", line 380, in <module>
    simulate()                                                                               # Calls the 'simulate' function
  File "/Users/apple/Documents/School/Computer Science/NEA/NEA Programming/Solar System Simulator.py", line 290, in simulate
    r = sqrt(((x1 - x2) ** 2) + ((y1 - y2) ** 2))
OverflowError: (34, 'Result too large')

Я знаю, что это может быть не очень эффективный кусок кода, но все, что мне нужно, это некоторая помощь в том, как это сделать таким образом. Проблема, кажется, возникает в вычислениях смещения или скорости. Есть идеи?

1 Ответ

0 голосов
/ 07 июня 2019

Начальная скорость получается дважды: в каждом цикле вызывается и используется начальная скорость, поэтому сама и новая скорость складываются вместе, каждый раз она увеличивается более чем в два раза.

Новаядобавленная скорость теперь вычисляется, а затем добавляется только один раз к начальной скорости.Планеты сейчас находятся на орбите.В настоящее время вокруг Солнца вращаются 9 тел.

...