diff --git a/sim_scenes/science/jupiter_moon_protects_earth.py b/sim_scenes/science/jupiter_moon_protects_earth.py new file mode 100644 index 0000000000000000000000000000000000000000..5877f248e709aa542d5883d74b2bd7504ec7c336 --- /dev/null +++ b/sim_scenes/science/jupiter_moon_protects_earth.py @@ -0,0 +1,217 @@ +# -*- coding:utf-8 -*- +# title :木星、月球保护地球模拟 +# description :木星、月球保护地球模拟 +# author :Python超人 +# date :2023-07-11 +# link :https://gitcode.net/pythoncr/ +# python_version :3.8 +# ============================================================================== +from bodies import Sun, Jupiter, Earth, Moon +from common.func import calculate_distance +from objs import RockSnow, Rock, create_rock +from common.consts import SECONDS_PER_MONTH, SECONDS_PER_YEAR, AU +from sim_scenes.func import ursina_run, camera_look_at, two_bodies_colliding, create_text_panel +from simulators.ursina.entities.body_timer import TimeData, BodyTimer +from simulators.ursina.entities.entity_utils import create_directional_light +from simulators.ursina.ursina_event import UrsinaEvent +import random +import math +import numpy as np + + +class JupiterProtectsEarthSim: + + def __init__(self, comet_num=20): + """ + + @param comet_num: 随机生成 comet_num 个石头 + """ + self.earth_moon_d = 30000000 # 因为地球放大 2000 倍，为了较好的效果，地月距离要比实际大才行 + # 分别保存太阳碰撞次数、木星碰撞次数、地球碰撞次数、月球碰撞次数、木星保护次数、月球保护次数 + self.colliding_count = [0, 0, 0, 0, 0, 0] + self.sun = Sun(name="太阳", size_scale=0.8e2) # 太阳放大 80 倍，距离保持不变 + self.jupiter = Jupiter(name="木星", size_scale=0.6e3) # 木星放大 600 倍，距离保持不变 + self.earth = Earth(name="地球", size_scale=2e3) # 地球放大 2000 倍，距离保持不变 + self.moon = Moon(name="月球", size_scale=3e3, # 月球球放大 3000 倍，为了较好的效果，地月距离要比实际大 + init_position=[self.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) + + self.bodies = [self.sun, self.jupiter, self.earth, self.moon] + + for body in self.bodies: + body.rotation_speed /= 10 # 自转速度降低10倍 + + self.comets = [] + + for i in range(comet_num): + # 随机生成 comet_num 个石头 + comet = self.create_comet(i, gravity_only_for=[self.sun, self.jupiter, self.earth]) + self.bodies.append(comet) + self.comets.append(comet) + + # 显示板信息模板 + self.colliding_info = "太阳碰撞：%s次（%s）\n\n木星碰撞：%s次（%s）\n\n地球碰撞：%s次（%s）\n\n月球碰撞：%s次（%s）\n\n木星保护：%s次\n\n月球保护：%s次\n\n" + + def random_pos_vel(self): + # 随机生成石头的位置和初始速度信息 + radius = calculate_distance(self.jupiter.position, self.sun.position) * 1.5 + x = radius * math.cos(random.uniform(0, 2 * math.pi)) * (random.randint(100, 150) / 100) + z = radius * math.sin(random.uniform(0, 2 * math.pi)) * (random.randint(100, 150) / 100) + sun_pos = self.sun.position + pos = [x + sun_pos[0], 0 + sun_pos[1], z + sun_pos[2]] + + # 随机速度 + vel = [-random.randint(90, 200) / 300, 0, 0] + return pos, vel + + def create_comet(self, index, gravity_only_for): + """ + + 随机生成石头（随机位置、随机初始速度、随机大小、随机旋转） + @param index: 索引号 + @param gravity_only_for: 指定多个天体，石头只对该天体引力有效 + @return: + """ + # 随机生成石头的位置和初始速度信息 + pos, vel = self.random_pos_vel() + + # vel = [0, 0, 0] + # 石头随机大小 + size_scale = random.randint(600, 1200) * 1.5e4 + # 随机创建石头 + rock = create_rock( + no=index % 8 + 1, name=f'岩石{index + 1}', mass=size_scale / 1000, + size_scale=size_scale, color=(255, 200, 0), + init_position=pos, init_velocity=vel, gravity_only_for=gravity_only_for + ) + + # 给石头一个随机旋转的方向和值 + rock.rotation = [0, 0, 0] + rock.rotation[random.randint(0, 2)] = random.randint(90, 200) / 100 + return rock + + def set_moon_position(self, time_data): + # 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 = self.earth_moon_d * np.cos(angle) + pz = self.earth_moon_d * np.sin(angle) + self.moon.position = [self.earth.position[0] + px, 0, self.earth.position[2] + pz] + + def on_timer_changed(self, time_data: TimeData): + self.set_moon_position(time_data) + # 运行中，每时每刻都会触发 + for comet in self.comets: + if comet.planet.enabled: + collided = False + # 如果是否可见，则旋转石头 + comet.planet.rotation += comet.rotation + # 循环判断每个石头与木星是否相碰撞，如果相碰撞就爆炸 + if not hasattr(comet, "jupiter_collided") and two_bodies_colliding(comet, self.jupiter): + # 碰撞到木星，该石头不要爆炸，尝试看看是否会碰撞到地球，如果碰撞了地球，则“木星保护”统计加一 + self.colliding_count[1] += 1 + # 加入标记，说明该石头已经碰撞到木星 + setattr(comet, "jupiter_collided", True) + # collided = True + elif not hasattr(comet, "moon_collided") and two_bodies_colliding(comet, self.moon): + # 碰撞到月球，该石头不要爆炸，尝试看看是否会碰撞到地球，如果碰撞了地球，则“月球保护”统计加一 + self.colliding_count[3] += 1 + # 加入标记，说明该石头已经碰撞到月球 + setattr(comet, "moon_collided", True) + elif two_bodies_colliding(comet, self.sun): + # 将石头隐藏、设置引力无效后，展示爆炸效果 + comet.explode(self.sun) + if not hasattr(comet, "jupiter_collided") and not hasattr(comet, "moon_collided"): + # 该石头没有碰撞到木星和月球，才算一次 + self.colliding_count[0] += 1 + collided = True + elif two_bodies_colliding(comet, self.earth): + if hasattr(comet, "jupiter_collided") or hasattr(comet, "moon_collided"): + if hasattr(comet, "jupiter_collided"): + # 说明该石头已经碰撞到木星,也就是说木星保护了地球一次 + self.colliding_count[4] += 1 + else: + # 说明该石头已经碰撞到月球,也就是说月球保护了地球一次 + self.colliding_count[5] += 1 + else: + # 该石头没有碰撞到木星和月球，才算一次 + self.colliding_count[2] += 1 + + comet.explode(self.earth) + collided = True + + if collided: + # 如果碰撞了，则该石头重复再利用，这样才保证有无限个石头可用 + pos, vel = self.random_pos_vel() + comet.init_position = pos + comet.init_velocity = vel + comet.set_visible(True) + comet.ignore_mass = False + comet.planet.enabled = True + if hasattr(comet, "jupiter_collided"): + delattr(comet, "jupiter_collided") + if hasattr(comet, "moon_collided"): + delattr(comet, "moon_collided") + + print("Sun:%s Jupiter:%s Earth:%s Moon:%s" % (self.colliding_count[0], + self.colliding_count[1], + self.colliding_count[2], + self.colliding_count[3])) + sun_cnt, jupiter_cnt, earth_cnt, moon_cnt, j_protected_cnt, m_protected_cnt = self.colliding_count + total_cnt = sun_cnt + jupiter_cnt + earth_cnt + moon_cnt + if total_cnt == 0: + colliding_info = self.colliding_info % (0, "0.0%", + 0, "0.0%", + 0, "0.0%", + 0, "0.0%", + 0, 0) + else: + colliding_info = self.colliding_info % (sun_cnt, str(round(sun_cnt * 100 / total_cnt, 2)) + "%", + jupiter_cnt, str(round(jupiter_cnt * 100 / total_cnt, 2)) + "%", + earth_cnt, str(round(earth_cnt * 100 / total_cnt, 2)) + "%", + moon_cnt, str(round(moon_cnt * 100 / total_cnt, 2)) + "%", + j_protected_cnt, m_protected_cnt) + + self.text_panel.text = colliding_info + + def on_ready(self): + # 运行前触发 + self.text_panel = create_text_panel() + self.text_panel.text = self.colliding_info % (0, "0.0%", + 0, "0.0%", + 0, "0.0%", + 0, "0.0%", + 0, 0) + + +if __name__ == '__main__': + """ + 木星、月球保护地球模拟 + """ + + # 设置计时器的最小时间单位为年 + BodyTimer().min_unit = BodyTimer.MIN_UNIT_YEARS + + sim = JupiterProtectsEarthSim(comet_num=20) + + # 运行前会触发 on_ready + UrsinaEvent.on_ready_subscription(sim.on_ready) + # 运行中，每时每刻都会触发 on_timer_changed + UrsinaEvent.on_timer_changed_subscription(sim.on_timer_changed) + + # 使用 ursina 查看的运行效果 + # 常用快捷键： P：运行和暂停 O：重新开始 I：显示天体轨迹 + # position = 左-右+、上+下-、前+后- + ursina_run(sim.bodies, SECONDS_PER_MONTH, + position=(30000000, 0, -3000000000), + show_timer=True) diff --git a/sim_scenes/science/jupiter_protects_earth.py b/sim_scenes/science/jupiter_protects_earth.py index d5d66b8825c3a7134dc26957778db4c3e0d615c7..ace4c1d2257298902dd70f2c406a64af5e63f0df 100644 --- a/sim_scenes/science/jupiter_protects_earth.py +++ b/sim_scenes/science/jupiter_protects_earth.py @@ -1,6 +1,6 @@ # -*- coding:utf-8 -*- -# title :彗木保护地球模拟 -# description :彗木保护地球模拟 +# title :木星保护地球模拟 +# description :木星保护地球模拟 # author :Python超人 # date :2023-06-27 # link :https://gitcode.net/pythoncr/ @@ -157,7 +157,7 @@ class JupiterProtectsEarthSim: if __name__ == '__main__': """ - 彗木保护地球模拟 + 木星保护地球模拟 """ # 设置计时器的最小时间单位为年