ursina_simulator.py 13.4 KB
Newer Older
三月三net's avatar
三月三net 已提交
1 2 3 4 5 6 7 8 9
# -*- coding:utf-8 -*-
# title           :ursina天体运行模拟器
# description     :ursina天体运行模拟器
# author          :Python超人
# date            :2023-02-11
# link            :https://gitcode.net/pythoncr/
# python_version  :3.8
# ==============================================================================
# pip install -i http://pypi.douban.com/simple/ --trusted-host=pypi.douban.com ursina
三月三net's avatar
三月三net 已提交
10
from ursina import Ursina, window, Entity, Grid, Mesh, camera, Text, application, color, mouse, Vec2, Vec3, \
三月三net's avatar
三月三net 已提交
11
    load_texture, held_keys, distance
M
march3 已提交
12
from ursina.prefabs.first_person_controller import FirstPersonController
三月三net's avatar
三月三net 已提交
13
import itertools
三月三net's avatar
三月三net 已提交
14
from simulators.ursina.ursina_event import UrsinaEvent
三月三net's avatar
三月三net 已提交
15 16 17 18
# from simulators.ursina.ursina_ui import UrsinaUI
from simulators.ursina.ui.control_ui import ControlUI
from simulators.ursina.ui.control_handler import ControlHandler

M
march3 已提交
19
from simulators.views.ursina_view import UrsinaView, UrsinaPlayer
三月三net's avatar
三月三net 已提交
20
from simulators.ursina.ursina_config import UrsinaConfig
三月三net's avatar
三月三net 已提交
21 22
from simulators.simulator import Simulator
from common.system import System
M
march3 已提交
23 24
import time
import datetime
三月三net's avatar
三月三net 已提交
25
import math
三月三net's avatar
三月三net 已提交
26
import os
三月三net's avatar
三月三net 已提交
27
from ursina import EditorCamera, PointLight, SpotLight, AmbientLight, DirectionalLight
三月三net's avatar
三月三net 已提交
28
from scenes.func import ursina_run
三月三net's avatar
三月三net 已提交
29

M
march3 已提交
30

三月三net's avatar
三月三net 已提交
31 32 33 34
class WorldGrid(Entity):  # Entity # 定义构造方法
    def __init__(self):
        super().__init__()
        s = 100
三月三net's avatar
三月三net 已提交
35
        grid = Entity(model=Grid(s, s), scale=s * 20, color=color.rgba(255, 255, 255, 20), rotation_x=90,
三月三net's avatar
三月三net 已提交
36
                      position=(0, -80, 0))
三月三net's avatar
三月三net 已提交
37 38 39 40 41
        # 坐标轴
        # vertsx = ((0, 0, 0), (10, 0, 0))
        # Entity(model=Mesh(vertices=vertsx, mode='line', thickness=3), color=color.cyan).set_light_off()
        # vertsyz = [(0, 0, 0), (0, 10, 0), (0, 0, 0), (0, 0, 10)]
        # Entity(model=Mesh(vertices=vertsyz, mode='line', thickness=3), color=color.yellow).set_light_off()
三月三net's avatar
三月三net 已提交
42 43 44
        grid.set_light_off()


三月三net's avatar
三月三net 已提交
45
class UrsinaSimulator(Simulator):
三月三net's avatar
三月三net 已提交
46 47 48
    """
    Ursina官网: https://www.ursinaengine.org/
    """
三月三net's avatar
三月三net 已提交
49
    def __init__(self, bodies_sys: System):
三月三net's avatar
三月三net 已提交
50 51 52 53 54 55
        self.app = Ursina()
        # import os
        # os.environ['CUDA_VISIBLE_DEVICES'] = '1'  # 选择第二个GPU
        # self.app = Ursina(window_title='GPU模拟',
        #              window_kwargs={'vsync': True, 'fullscreen': False, 'borderless': False, 'show_ursina_splash': True,
        #                             'high_resolution': True})
三月三net's avatar
三月三net 已提交
56

M
march3 已提交
57
        self.ursina_views = []
三月三net's avatar
三月三net 已提交
58 59 60 61 62 63 64
        window.color = color.black

        super().__init__(bodies_sys, UrsinaView)

        # ps = ["sun", "mercury", "venus", "earth", "mars", "jupiter", "saturn", "uranus", "neptune"]
        # cp = [200, 15, 35, 42, 20, 160, 145, 90, 80]
        # x, y, z = 0, 0, 0
M
march3 已提交
65
        for i, view in enumerate(self.body_views):
三月三net's avatar
三月三net 已提交
66
            # pos = tuple(body.position)
M
march3 已提交
67 68
            # ursina_view = UrsinaView(body)
            view.update()
M
march3 已提交
69
            self.ursina_views.append(view)
三月三net's avatar
三月三net 已提交
70 71
            # planets.append(newPlanet)
            # x += cp[i] * 10
三月三net's avatar
三月三net 已提交
72
        self.adj_run_params()
三月三net's avatar
三月三net 已提交
73 74
        UrsinaEvent.on_searching_bodies_subscription(type(self).__name__, self.on_searching_bodies)

三月三net's avatar
三月三net 已提交
75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93
    # def get_bodies_max_distance(self, body_views):
    #     max_distance = 0
    #     for b1 in body_views:
    #         if b1.body.ignore_mass:
    #             continue
    #         for b2 in body_views:
    #             if (b1 is b2) or b2.body.ignore_mass:
    #                 continue
    #         d = distance(b1.planet, b2.planet)
    #         if d > max_distance:
    #             max_distance = d
    #     return max_distance

    def get_bodies_max_distance(self, body_views):
        """
        算法优化
        :param body_views:
        :return:
        """
三月三net's avatar
三月三net 已提交
94
        max_distance = 0
三月三net's avatar
三月三net 已提交
95 96 97
        for b1, b2 in itertools.combinations(body_views, 2):
            if b1.body.ignore_mass or b2.body.ignore_mass:
                continue
三月三net's avatar
三月三net 已提交
98 99 100
            d = distance(b1.planet, b2.planet)
            if d > max_distance:
                max_distance = d
三月三net's avatar
三月三net 已提交
101 102 103 104 105 106 107 108
        return max_distance

    def adj_run_params(self):
        """

        :return:
        """
        max_distance = self.get_bodies_max_distance(self.body_views)
三月三net's avatar
三月三net 已提交
109 110

        # UrsinaConfig.control_camera_speed = round(max_distance * 10, 2)
三月三net's avatar
三月三net 已提交
111
        time_scale = round(pow(max_distance, 1 / 4), 2)
三月三net's avatar
三月三net 已提交
112 113
        if time_scale < 0.01:
            time_scale = 0.01
三月三net's avatar
三月三net 已提交
114

三月三net's avatar
三月三net 已提交
115
        application.time_scale = time_scale
三月三net's avatar
三月三net 已提交
116
        # UrsinaConfig.auto_scale_factor = 1.0e-9
三月三net's avatar
三月三net 已提交
117

三月三net's avatar
三月三net 已提交
118 119 120
    def on_searching_bodies(self, **kwargs):
        views = []
        for view in self.body_views:
三月三net's avatar
三月三net 已提交
121 122
            # if view.appeared:
            views.append(view)
三月三net's avatar
三月三net 已提交
123 124
        return views

M
march3 已提交
125 126 127 128
    def check_elapsed_time(self):
        """检查时间间隔是否已过"""
        now = datetime.datetime.now()
        elapsed_time = now - self.last_time
M
march3 已提交
129 130 131 132
        value = elapsed_time >= self.interval
        if value:
            self.last_time = now
        return value
M
march3 已提交
133 134 135

    def check_and_evolve(self):
        if self.check_elapsed_time():
三月三net's avatar
三月三net 已提交
136
            run_speed_factor = UrsinaConfig.run_speed_factor
三月三net's avatar
三月三net 已提交
137 138 139 140 141
            if UrsinaConfig.seconds_per == 0:
                evolve_dt = self.evolve_dt * run_speed_factor
            else:
                evolve_dt = UrsinaConfig.seconds_per * run_speed_factor

三月三net's avatar
三月三net 已提交
142
            evolve_dt = evolve_dt * self.interval_fator
三月三net's avatar
三月三net 已提交
143
            super().evolve(evolve_dt)
M
march3 已提交
144

三月三net's avatar
三月三net 已提交
145
    def cosmic_background(self, texture='../textures/cosmic2.jpg'):
三月三net's avatar
三月三net 已提交
146 147 148 149 150
        """
        加入宇宙背景
        :param texture:
        :return:
        """
三月三net's avatar
三月三net 已提交
151 152
        # Add skybox
        from ursina import Sky
三月三net's avatar
三月三net 已提交
153 154
        Sky(texture=texture).scale = 10000

三月三net's avatar
三月三net 已提交
155 156 157 158 159
        # texture = load_texture(texture)
        # sky_dome = Entity(model='sky_dome', texture=texture, scale=10000,
        #                   color=color.white,
        #                   position=(0, 0, 0),
        #                   rotation=(0, 0, 0))
三月三net's avatar
三月三net 已提交
160

三月三net's avatar
三月三net 已提交
161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178
    def __add_glow(self, entity, intensity=2, light_color=color.white, attenuation=3):
        """
        未用,保留代码
        :param entity:
        :param intensity:
        :param light_color:
        :param attenuation:
        :return:
        """
        lights = []
        import math
        for i in range(5):
            glow_entity = Entity(parent=entity, model='sphere', color=color.rgba(1.0, 0.6, 0.2, 1),
                                 scale=math.pow(1.03, i), alpha=0.2)
            lights.append(glow_entity)
        # 创建一个新的 Entity 对象,作为光晕的容器
        # glow_entity = Entity(parent=entity, model='sphere', scale=entity.scale * 1.2)
        # 创建 PointLight 对象,并设置它的属性
三月三net's avatar
三月三net 已提交
179 180 181
        for i in range(2):
            light = PointLight(parent=lights[0], intensity=intensity, color=light_color, attenuation=attenuation)
            lights.append(light)
三月三net's avatar
三月三net 已提交
182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201

        # 把 Entity 对象放到星星的后面,使得光晕看起来像是从星星发出来的
        glow_entity.world_position = entity.world_position
        glow_entity.world_parent = entity.parent
        glow_entity.y += entity.scale_y * 0.1
        glow_entity.depth_test = False
        return lights

    def create_fixed_star_lights(self, entity):
        """
        创建恒星的发光的效果、并作为灯光源
        :param entity:
        :return:
        """

        # 如果是恒星(如:太阳),自身会发光,则需要关闭灯光
        entity.set_light_off()

        lights = []
        # 创建多个新的 Entity 对象,作为光晕的容器
三月三net's avatar
三月三net 已提交
202
        for i in range(10):
三月三net's avatar
三月三net 已提交
203
            glow_entity = Entity(parent=entity, model='sphere', color=color.rgba(1.0, 0.6, 0.2, 1),
三月三net's avatar
三月三net 已提交
204
                                 scale=math.pow(1.03, i), alpha=0.1)
三月三net's avatar
三月三net 已提交
205 206

            lights.append(glow_entity)
三月三net's avatar
三月三net 已提交
207 208 209 210
        for i in range(2):
            # 创建 PointLight 对象,作为恒星的灯光源
            light = PointLight(parent=entity, intensity=10, range=10, color=color.white)
            lights.append(light)
三月三net's avatar
三月三net 已提交
211 212 213 214 215

        # light = DirectionalLight(shadows=True, direction=Vec3(0, 0, 1), color=color.white)
        # light.look_at(Vec3(0, 0, -1))
        # light = SpotLight(parent=entity,shadows=True, direction=Vec3(1,1,1), color=color.white)

三月三net's avatar
三月三net 已提交
216 217
        return lights

三月三net's avatar
三月三net 已提交
218 219 220 221 222 223 224 225 226 227 228
    def switch_to_english_input_method(self):
        pass
        # os.system("powershell -Command \"Set-WinUserLanguageList -LanguageList 'en-US', 'zh-CN' -Force\"")

        # if os.name == 'nt':  # Windows 系统
        #     os.system('powershell -Command "Set-WinUserLanguageList en-US -Force"')
        # elif os.name == 'posix':  # macOS 或 Linux 系统
        #     os.system(
        #         'defaults write ~/Library/Preferences/com.apple.HIToolbox.plist AppleSelectedInputSources -array-add \'{ "InputSourceKind" = "Keyboard Layout"; "KeyboardLayout ID" = 252; "KeyboardLayout Name" = "ABC"; }\'')
        #     os.system('killall cfprefsd')  # 重启 CoreFoundation 进程使更改生效

三月三net's avatar
三月三net 已提交
229
    def run(self, dt, **kwargs):
三月三net's avatar
三月三net 已提交
230 231

        window.title = '宇宙模拟器'
三月三net's avatar
三月三net 已提交
232 233 234 235 236 237 238 239

        # 设置 camera 的裁剪面和位置
        # camera.clip_plane_near = 0.00000009
        # camera.fov = 120
        # camera.clip_plane_far = 1000
        # camera.position = (0, 10, -20)
        # camera.rotation_x = -30

三月三net's avatar
三月三net 已提交
240
        # 设定时间间隔为0.01秒
三月三net's avatar
三月三net 已提交
241 242 243
        self.interval_fator = 0.01
        self.evolve_dt = dt
        self.interval = datetime.timedelta(seconds=self.interval_fator)
三月三net's avatar
三月三net 已提交
244

M
march3 已提交
245
        self.last_time = datetime.datetime.now() - datetime.timedelta(seconds=2)
三月三net's avatar
三月三net 已提交
246 247
        if "light" in kwargs:
            if kwargs["light"]:
三月三net's avatar
三月三net 已提交
248 249
                for v in self.ursina_views:
                    if v.body.is_fixed_star:
三月三net's avatar
三月三net 已提交
250 251
                        # self.lights = self.create_fixed_star_lights(v.planet)
                        pass
三月三net's avatar
三月三net 已提交
252

三月三net's avatar
三月三net 已提交
253 254 255 256
        if "show_grid" in kwargs:
            if kwargs["show_grid"]:
                WorldGrid()

三月三net's avatar
三月三net 已提交
257 258
        if "cosmic_bg" in kwargs:
            cosmic_bg = kwargs["cosmic_bg"]
三月三net's avatar
三月三net 已提交
259
            if cosmic_bg is None:
三月三net's avatar
三月三net 已提交
260
                # cosmic_bg = '../textures/cosmic1.png'
三月三net's avatar
三月三net 已提交
261 262
                # cosmic_bg = '../textures/cosmic2.jpg'
                cosmic_bg = '../textures/cosmic3.jpg'
三月三net's avatar
三月三net 已提交
263 264 265 266
            import os
            if cosmic_bg is not None and os.path.exists(cosmic_bg):
                self.cosmic_background(cosmic_bg)

三月三net's avatar
三月三net 已提交
267 268 269
        # ui = UrsinaUI()
        ctl = ControlUI(ControlHandler(), position=(0.6, 0.5))

三月三net's avatar
三月三net 已提交
270
        EditorCamera(ignore_paused=True)
三月三net's avatar
三月三net 已提交
271 272
        # 防止打开中文输入法
        # self.switch_to_english_input_method()
三月三net's avatar
三月三net 已提交
273 274
        self.app.run()

M
march3 已提交
275

三月三net's avatar
三月三net 已提交
276
if __name__ == '__main__':
M
march3 已提交
277
    from bodies import Sun, Mercury, Venus, Earth, Mars, Jupiter, Saturn, Uranus, Neptune, Pluto, Moon
三月三net's avatar
三月三net 已提交
278
    from bodies.body import AU
三月三net's avatar
三月三net 已提交
279
    from common.consts import SECONDS_PER_WEEK, SECONDS_PER_DAY, SECONDS_PER_HALF_DAY
三月三net's avatar
三月三net 已提交
280 281 282 283 284

    """
    3个太阳、1个地球
    """
    bodies = [
M
march3 已提交
285
        Sun(name='太阳1', mass=1.5e30, init_position=[849597870.700, 0, 0], init_velocity=[0, 7.0, 0],
三月三net's avatar
三月三net 已提交
286
            size_scale=5e1, texture="sun.png"),  # 太阳放大 100 倍
M
march3 已提交
287
        Sun(name='太阳2', mass=2e30, init_position=[0, 0, 0], init_velocity=[0, -8.0, 0],
三月三net's avatar
三月三net 已提交
288
            size_scale=5e1, texture="sun.png"),  # 太阳放大 100 倍
M
march3 已提交
289
        Sun(name='太阳3', mass=2.5e30, init_position=[0, -849597870.700, 0], init_velocity=[18.0, 0, 0],
三月三net's avatar
三月三net 已提交
290
            size_scale=5e1, texture="sun.png"),  # 太阳放大 100 倍
M
march3 已提交
291
        Earth(name='地球', init_position=[0, -349597870.700, 0], init_velocity=[15.50, 0, 0],
三月三net's avatar
三月三net 已提交
292 293
              size_scale=4e3, texture="earth.png", distance_scale=1),  # 地球放大 4000 倍,距离保持不变
    ]
M
march3 已提交
294 295 296 297 298 299 300 301 302 303 304
    # bodies = [
    #
    #     Sun(name='太阳2', mass=1.5e30, init_position=[0, 0, 0], init_velocity=[0, -8.0, 0],
    #         size_scale=5e1, texture="sun.png"),  # 太阳放大 100 倍
    #     Sun(name='太阳2', mass=1.5e30, init_position=[849597870.700, 0, 0], init_velocity=[0, -8.0, 0],
    #         size_scale=5e1, texture="sun.png"),  # 太阳放大 100 倍
    #     Sun(name='太阳2', mass=1.5e30, init_position=[0, -849597870.700, 0], init_velocity=[0, -8.0, 0],
    #         size_scale=5e1, texture="sun.png"),  # 太阳放大 100 倍
    #     Earth(name='地球', mass=5.97237e24, init_position=[0, -349597870.700, 0], init_velocity=[15.50, 0, 0],
    #           size_scale=4e3, texture="earth.png", distance_scale=1),  # 地球放大 4000 倍,距离保持不变
    # ]
M
march3 已提交
305
    bodies = [
三月三net's avatar
三月三net 已提交
306 307 308
        Sun(size_scale=0.8e2),  # 太阳放大 80 倍
        Mercury(size_scale=4e3, distance_scale=1.3),  # 水星放大 4000 倍,距离放大 1.3 倍
        Venus(size_scale=4e3, distance_scale=1.3),  # 金星放大 4000 倍,距离放大 1.3 倍
三月三net's avatar
三月三net 已提交
309 310 311
        Earth(init_position=[1.12 * AU, 0, 0],
              init_velocity=[0, 29.79, 0], size_scale=4e3, distance_scale=1.3),  # 地球放大 4000 倍,距离放大 1.3 倍
        Moon(init_position=[363104 + 1.12 * AU, 0, 0],
三月三net's avatar
三月三net 已提交
312
             init_velocity=[-9, 29.79 + 1.023, 0], size_scale=4e3, distance_scale=1.3),
三月三net's avatar
三月三net 已提交
313
        Mars(size_scale=4e3, distance_scale=1.3),  # 火星放大 4000 倍,距离放大 1.3 倍
M
march3 已提交
314
        Jupiter(size_scale=0.68e3, distance_scale=0.65),  # 木星放大 680 倍,距离缩小到真实距离的 0.65
三月三net's avatar
三月三net 已提交
315 316 317 318
        Saturn(size_scale=0.68e3, distance_scale=0.52),  # 土星放大 680 倍,距离缩小到真实距离的 0.52
        Uranus(size_scale=0.8e3, distance_scale=0.36),  # 天王星放大 800 倍,距离缩小到真实距离的 0.36
        Neptune(size_scale=1e3, distance_scale=0.27),  # 海王星放大 1000 倍,距离缩小到真实距离的 0.27
        Pluto(size_scale=10e3, distance_scale=0.23),  # 冥王星放大 10000 倍,距离缩小到真实距离的 0.23(从太阳系的行星中排除)
M
march3 已提交
319
    ]
M
march3 已提交
320

三月三net's avatar
三月三net 已提交
321
    ursina_run(bodies, SECONDS_PER_DAY, position=(AU * 2, AU * 2, AU * 3), show_grid=True)