From d4ce23a07e63b202d15c9b24872749c464c9d792 Mon Sep 17 00:00:00 2001 From: Ying-Li Niu <64801511@qq.com> Date: Thu, 18 Jun 2026 23:03:36 +0800 Subject: [PATCH] =?UTF-8?q?feat:=20case04=20=E4=BB=A5=E6=9C=88=E7=90=83?= =?UTF-8?q?=E8=B4=A8=E9=87=8F=3D1=20=E4=B8=BA=E5=9F=BA=E5=87=86=E6=9B=B4?= =?UTF-8?q?=E6=96=B0=E8=B4=A8=E9=87=8F=E6=AF=94=E5=92=8C=E8=B4=A8=E9=87=8F?= =?UTF-8?q?=E7=BC=A9=E6=94=BE?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit - coord.txt: 太阳 27000 / 地球 81 / 月球 1（月球=1 质量单位） - 地球绕日轨道速度 v=520（v = sqrt(g_strength * M_sun / r)） - 月球绕地速度 v=526（地球速度 + 6 绕地附加速度） - box_a: 20 → 30（容纳更大轨道半径） - Readme.html 同步更新质量表和初始构型表 - 已知局限说明太阳质量缩至 1/1000 的原因 --- examples/case04/Readme.html | 27 ++++++++++++++------------- examples/case04/input/coord.txt | 6 +++--- examples/case04/input/input.txt | 2 +- 3 files changed, 18 insertions(+), 17 deletions(-) diff --git a/examples/case04/Readme.html b/examples/case04/Readme.html index 22d9a9b..fded8aa 100644 --- a/examples/case04/Readme.html +++ b/examples/case04/Readme.html @@ -211,14 +211,15 @@ MathJax = {

2.1 质量

- - - - + + + +
天体质量 / kg与地球质量之比
太阳$1.989 \times 10^{30}$$333\,000$
地球$5.972 \times 10^{24}$$1$
月球$7.35 \times 10^{22}$$0.0123 \;(1/81.3)$
天体质量 / kg模拟质量（月球=1）
太阳$1.989 \times 10^{30}$$27\,000$
地球$5.972 \times 10^{24}$$81$
月球$7.35 \times 10^{22}$$1$
- 质量比：$M_\odot : M_\oplus : M_\text{moon} \approx 27\,100\,000 : 81.3 : 1$ + 质量比：$M_\odot : M_\oplus : M_\text{moon} \approx 27\,100\,000 : 81.3 : 1$
+ 模拟中采用缩放比例 $27\,000 : 81 : 1$（太阳质量缩至 $1/1000$ 以保持数值稳定）。

2.2 轨道参数

@@ -310,8 +311,8 @@ MathJax = { - - + + @@ -321,12 +322,12 @@ MathJax = {
参数值说明
gravity_strength100.0万有引力强度（已缩放）
box_a20.0盒子半边长
gravity_strength100.0万有引力强度
box_a30.0盒子半边长
methodleapfrog蛙跳法（辛积分器，能量守恒）
T_total10.0 s总模拟时间
NSTEP2抽帧间隔（密采帧）
- - - + + +
天体质量位置 $(x,y,z)$速度 $(v_x,v_y,v_z)$固定约束
太阳1$(0,0,0)$$(0,0,0)$全部固定
地球1$(4,0,0)$$(0,0,4)$无
月球0.1$(5,0,0)$$(0,0,6)$无
太阳$27\,000$$(0,0,0)$$(0,0,0)$全部固定
地球$81$$(10,0,0)$$(0,0,520)$无
月球$1$$(11.5,0,0)$$(0,0,526)$无
-

地球在 $z$ 方向获得初速度，产生绕太阳的轨道运动；月球在地球基础上获得额外的 $z$ 方向速度，形成绕地球的轨道。

+

地球在 $z$ 方向获得初速度 $v=520$，产生绕太阳的轨道运动；月球在地球基础上附加 $v=6$ 的绕地速度，形成绕地球的微小轨道。速度由圆形轨道公式 $v = \sqrt{G_{\text{eff}} M / r}$ 计算，其中 $G_{\text{eff}} = \text{gravity\_strength} = 100$。

五、使用方法

@@ -357,8 +358,8 @@ MathJax = {

七、已知局限

diff --git a/examples/case04/input/coord.txt b/examples/case04/input/coord.txt index 7b5d01a..036ade1 100644 --- a/examples/case04/input/coord.txt +++ b/examples/case04/input/coord.txt @@ -1,4 +1,4 @@ n mass radius x y z vx vy vz fix_x fix_y fix_z -1 1 0.28 0 0 0 0 0 0 1 1 1 -2 1 0.28 4 0 0 0 0 4 0 0 0 -3 0.1 0.18 5 0 0 0 0 6 0 0 0 \ No newline at end of file +1 27000 2.0 0 0 0 0 0 0 1 1 1 +2 81 0.5 10 0 0 0 0 520 0 0 0 +3 1 0.3 11.5 0 0 0 0 526 0 0 0 \ No newline at end of file diff --git a/examples/case04/input/input.txt b/examples/case04/input/input.txt index 6d70e56..1bf6fbd 100644 --- a/examples/case04/input/input.txt +++ b/examples/case04/input/input.txt @@ -23,7 +23,7 @@ force_calc: 0 # 强制重新计算：1=跳过缓存强算，0=自动使用 engine: python # 默认使用 Python 引擎 # ── 盒子 ────────────────────────────────────── -box_a: 20.0 # 立方体半边长，粒子被限制在 [-box_a, box_a]³ 内 +box_a: 30.0 # 立方体半边长，粒子被限制在 [-box_a, box_a]³ 内 # ── 初始构型 ────────────────────────────────── # 坐标文件格式：