# -*- coding:utf-8 -*- # title :地月场景模拟(观看月相变化的过程) # description :地月场景模拟(观看月相变化的过程) # author :Python超人 # date :2023-02-11 # link :https://gitcode.net/pythoncr/ # python_version :3.8 # ============================================================================== from bodies import Sun, Earth, Moon from common.consts import SECONDS_PER_HOUR, SECONDS_PER_HALF_DAY, SECONDS_PER_DAY, SECONDS_PER_WEEK, SECONDS_PER_MONTH from sim_scenes.func import mayavi_run, ursina_run, camera_look_at, set_camera_parent from bodies.body import AU from simulators.ursina.entities.body_timer import TimeData from simulators.ursina.ursina_config import UrsinaConfig from simulators.ursina.ursina_event import UrsinaEvent import numpy as np if __name__ == '__main__': """ 地月场景模拟 """ # # 地球的Y方向初始速度 # EARTH_INIT_VELOCITY = 0 # 0km/s # sun = Sun(init_position=[0, 0, AU], size_scale=2e1) # 太阳放大 20 倍 # # 忽略质量的引力 # sun.ignore_mass = True # # # 观看月相变化的过程:分别是 新月、蛾眉月、上弦月、盈凸、满月、亏凸、下弦月、残月 # # 参考:images/moon/月相变化过程.jpeg # # TODO: 月球在摄像机的前方(从 “新月” 开始) # moon_pos, moon_vel = [0, 0, 384400], [-(EARTH_INIT_VELOCITY + 1.03), 0, 0] # # TODO: 月球在摄像机的右方(从 “下弦月” 开始),将会从右方出现 # # moon_pos, moon_vel = [384400, 0, 0], [-EARTH_INIT_VELOCITY, 0, 1.03] # # TODO: 月球在摄像机的左方(从 “上弦月” 开始) # # moon_pos, moon_vel = [-384400, 0, 0], [EARTH_INIT_VELOCITY, 0, -1.03] # moon = Moon(init_position=moon_pos, # 距地距离约: 363104 至 405696 km # init_velocity=moon_vel, # size_scale=2e1) # 月球放大 20 倍,距离保持不变 # bodies = [ # sun, # Earth(init_position=[0, 0, 0], # init_velocity=[0, EARTH_INIT_VELOCITY, 0], # size_scale=1e1), # 地球放大 10 倍,距离保持不变 # moon # ] sun = Sun(name="太阳", size_scale=1e2) # 太阳放大 100 倍,距离保持不变 earth = Earth(name="地球", size_scale=1.8e3) # 地球放大 1800 倍,距离保持不变 earth_moon_d = 20000000 # 因为地球放大 1800 倍,为了较好的效果,地月距离要比实际大才行 moon = Moon(name="月球", size_scale=3e3, # 月球球放大 3000 倍,为了较好的效果,地月距离要比实际大 init_position=[earth_moon_d, 0, AU], init_velocity=[0, 0, 0], ignore_mass=True, rotation_speed=0.4065, # gravity_only_for_earth=True )# .set_light_disable(True) bodies = [ sun, earth, moon ] def on_timer_changed(time_data: TimeData): # 显示轨迹,并设置轨迹长度 UrsinaConfig.show_trail = True UrsinaConfig.trail_length = 800 # 设置运行速度(加速10倍) UrsinaConfig.run_speed_factor = 10 # 1个月有29.5天 days_per_month = 29.5 # 1天多少角度 angle_per_day = 360 / days_per_month # 当前天数的角度(度) angle = time_data.total_days * angle_per_day # 当前天数的角度(弧度) angle = angle * np.pi / 180 # 计算月亮的坐标(这里没有用到万有引力) px = earth_moon_d * np.cos(angle) pz = earth_moon_d * np.sin(angle) moon.position = [earth.position[0] + px, 0, earth.position[2] + pz] def on_ready(): camera_look_at(sun) # 订阅事件后,上面2个函数功能才会起作用 # 运行前会触发 on_ready UrsinaEvent.on_ready_subscription(on_ready) # 运行中,每时每刻都会触发 on_timer_changed UrsinaEvent.on_timer_changed_subscription(on_timer_changed) # 使用 ursina 查看的运行效果 # 常用快捷键: P:运行和暂停 O:重新开始 I:显示天体轨迹 # position = 左-右+、上+下-、前+后- # position=(0, 0, 0) 的位置是站在地球视角,可以观看月相变化的过程 ursina_run(bodies, SECONDS_PER_DAY, position=(0, 4 * AU, 0), show_timer=True)