three_body_01.py 1.8 KB
Newer Older
M
march3 已提交
1
# -*- coding:utf-8 -*-
M
march3 已提交
2 3
# title           :三体场景模拟01
# description     :三体场景模拟(3个太阳、1个地球)
M
march3 已提交
4
# author          :Python超人
M
march3 已提交
5 6
# date            :2023-02-11
# link            :https://gitcode.net/pythoncr/
M
march3 已提交
7 8
# python_version  :3.8
# ==============================================================================
M
march3 已提交
9
from bodies import Sun, Earth
三月三net's avatar
三月三net 已提交
10
from common.consts import SECONDS_PER_WEEK, SECONDS_PER_YEAR, SECONDS_PER_MONTH, AU
三月三net's avatar
三月三net 已提交
11
from scenes.func import mayavi_run, ursina_run
M
march3 已提交
12 13

if __name__ == '__main__':
M
march3 已提交
14 15 16 17 18 19 20
    """
    3个太阳、1个地球(效果1)
    可以修改影响效果的参数为: 
    1、三个方向的初始位置 init_position[x, y, z]
    2、三个方向的初始速度 init_velocity[x, y, z]
    3、天体质量 mass    
    """
M
march3 已提交
21 22
    bodies = [
        Sun(mass=1.5e30, init_position=[849597870.700, 0, 0], init_velocity=[0, 7.0, 0],
三月三net's avatar
三月三net 已提交
23
            size_scale=5e1, texture="sun2.jpg"),  # 太阳放大 100 倍
M
march3 已提交
24
        Sun(mass=2e30, init_position=[0, 0, 0], init_velocity=[0, -8.0, 0],
M
march3 已提交
25
            size_scale=5e1, texture="sun2.jpg"),  # 太阳放大 100 倍
M
march3 已提交
26
        Sun(mass=2.5e30, init_position=[0, -849597870.700, 0], init_velocity=[18.0, 0, 0],
M
march3 已提交
27
            size_scale=5e1, texture="sun2.jpg"),  # 太阳放大 100 倍
M
march3 已提交
28
        Earth(init_position=[0, -349597870.700, 0], init_velocity=[15.50, 0, 0],
M
march3 已提交
29
              size_scale=4e3, distance_scale=1),  # 地球放大 4000 倍,距离保持不变
M
march3 已提交
30
    ]
三月三net's avatar
三月三net 已提交
31
    # 使用 mayavi 查看的运行效果
三月三net's avatar
三月三net 已提交
32 33 34
    # mayavi_run(bodies, SECONDS_PER_WEEK, view_azimuth=0)

    # 使用 ursina 查看的运行效果
三月三net's avatar
三月三net 已提交
35
    # 常用快捷键: P:运行和暂停  O:重新开始  I:显示天体轨迹
三月三net's avatar
三月三net 已提交
36
    # position = 左-右+、上+下-、前+后-
三月三net's avatar
三月三net 已提交
37
    ursina_run(bodies, SECONDS_PER_YEAR, position=(3 * AU, AU, -4 * AU), show_trail=True)