太阳辐射压 SRP 模拟器 — 航天器轨道摄动与帆推进

参数设置

航天器质量 m

受光有效面积 A

m²

垂直于太阳光的投影面积 (电池板+结构)

反射率 ρ

0=完全吸收 / 1=完全镜面反射

吸收系数 c

表面吸收的比例 (热分析的参考值)

太阳距离 d

水星 0.39 / 地球 1.0 / 火星 1.52 / 木星 5.2

任务周期

day

轨道类型

SRP 是否占主导地位取决于轨道

计算结果

—

SRP 压 (μN/m²)

—

SRP 力 (mN)

—

SRP 加速度 (μm/s²)

—

β 参数

—

累积 ΔV (m/s)

—

轨道影响

—

SRP 示意图 — 光子、航天器、摄动矢量

太阳辐射的光子流打到航天器受光面，SRP 力矢量 (橙色) 沿反向太阳方向作用。颜色表示 β 的大小。

SRP 力 vs 太阳距离 (AU)

反射率别 SRP 力对比

理论·主要公式

$$P_{SRP} = \frac{\Phi}{c},\quad F = (1+\rho)\,P\,A,\quad \beta = \frac{a_{SRP}}{a_{grav}} = \frac{(1+\rho)\,\Phi\,A}{c\,m\,a_g}$$

Φ = 太阳流量 (W/m²，1 AU 处为 1361)、c = 光速 2.998×10⁸ m/s、ρ = 反射率、A = 受光面积 (m²)、m = 航天器质量 (kg)、a_g = 太阳引力加速度。β > 1 时航天器被太阳推离。

$$\Phi(d) = \frac{\Phi_{1\,AU}}{d^{2}},\qquad \Delta V_{cum} = a_{SRP} \cdot t_{mission}$$

流量与距离平方成反比 (逆平方律)。累积 ΔV 是加速度乘以任务秒数的概算值。实际中受入射角、地球阴影 (eclipse)、自转影响，时间平均值会降低。

太阳辐射压 (SRP) 与航天器摄动

🙋

「太阳辐射压」就是太阳光推动航天器对吧？可是光没有质量，怎么会产生力呢？

🎓

好问题！光子虽然静止质量为零，但具有能量 E 和动量 p = E/c。当光子打到航天器表面进行反射和吸收时，根据动量守恒，必然将动量传递给航天器。在 1 AU 处，太阳的能量流量约为 1361 W/m²，用光速 c 除得到的压力约为 4.54 μN/m²。确实，每平方米只有几微牛顿，极其微弱。但在真空中没有摩擦，这个力会不断积累，经过数月乃至数年，就能产生数 m/s 甚至数十 m/s 的速度变化。这对 GEO 静止卫星的保持和深空探测任务的轨道确定都是不可忽视的摄动。

🙋

那普通的通信卫星这样的人工卫星，SRP 真的有影响吗？左边选「GEO」时显示「对 GEO 位置保持重要」。

🎓

对，GEO 静止卫星正是 SRP 影响明显的例子。GEO 卫星的太阳电池板很大 (常见的 50 m² 级)，相对质量而言受光面积比 A/m 很大。因此 SRP 加速度在白天和晚上会不断变化，导致轨道离心率增减变化。从年度看，这是东西向保持 (E-W station-keeping) 的主要干扰之一，需要消耗约 1〜3 m/s/年的 ΔV 来补正。这直接影响卫星的燃料寿命，所以轨道设计阶段必须考虑。用默认参数 (1000 kg、20 m²)，一年就会产生约 3.7 m/s 的累积 ΔV，就是这个原因。

🙋

什么是 β 参数？右边显示「β = 1.99e-5」，这个数字大的话会怎样？

🎓

β 是「SRP 加速度与太阳引力加速度的比值」，是无量纲量，用来衡量 SRP 相对重力有多强。定义为 β = a_SRP / a_grav = (1+ρ)·Φ·A / (c·m·a_g)。普通卫星的 β 约为 1e-5 〜 1e-4，非常小，SRP 被作为摄动处理。但太阳帆的目标是让 β = 0.01〜0.1 这个量级。当 β = 1 时，SRP 加速度等于太阳引力，航天器开始被太阳推离轨道。JAXA 的 IKAROS (2010) 是 14 m 见方、约 310 kg 的太阳帆，实现了 β ≈ 1.4e-4，是世界上首次完全用太阳帆进行行星间飞行的验证。试试把左边的轨道类型改为「太阳帆」，再把面积改大、质量改小，你会看到 β 急剧上升。

🙋

拉格朗日点任务中，SRP 的问题在哪里？我听说 JWST 在那里。

🎓

JWST、Gaia、SOHO 都在太阳-地球拉格朗日点 (L1 或 L2) 运行。拉格朗日点本质上是不稳定的「鞍点」，微小的摄动就会导致轨道漂移。JWST 在 L2，有巨大的遮阳板 (约 22 m × 12 m)，受光面积很大。更糟的是，遮阳板的反射特性不均匀，SRP 的合力矢量会偏离质量中心，还会产生转矩。为了补正这些，JWST 每隔几个月要进行一次保持喷射，这是消耗推进剂的主要原因之一，直接影响卫星的使用寿命。

🙋

太阳帆的情况下，SRP 就不再是「麻烦」而是「主推力」了。这样真的能到深空吗？

🎓

原理上是可以的。最大的优势是「推进剂质量为零」。常规化学火箭为了获得 ΔV，推进剂需求呈指数增长。但太阳帆只要时间足够，理论上可以获得无限的 ΔV。实际案例包括 IKAROS (金星方向)、NASA 的 NEA Scout、LightSail 2，还有即将进行的 Solar Cruiser 等。不过加速度极小 (β=1e-4 时约 5e-7 m/s²)，到达目标地点需要好几年。但在小行星勘测、太阳极地轨道等化学推进难以到达的任务中，太阳帆逐渐成为可行方案，相关研究在持续推进。

常见问题

光子具有能量 E 和动量 p = E/c，当它们打到航天器表面进行反射和吸收时，会传递动量。1 AU 处太阳辐射流量约 1361 W/m²，用光速 c 除得到的压力约 4.54 μN/m²，极其微小。但在没有大气的真空中没有摩擦力，经年累积可达数 m/s 至数十 m/s 的 ΔV，成为 GEO 静止卫星保持和深空任务轨道确定中不可忽视的摄动。

β 是 SRP 加速度与太阳引力加速度的比值，是无量纲量，定义为 β = a_SRP / a_grav = (1+ρ)·Φ·A / (c·m·a_g)。常规通信卫星的 β 约为 1e-5 〜 1e-4，SRP 作为摄动处理。而太阳帆的 β 以 0.01 〜 0.1 为目标，IKAROS 约为 1.4e-4，当 β 超过 1 时，航天器能被太阳辐射压推离太阳。

GEO 静止卫星相对质量而言太阳电池板面积大，SRP 加速度的日周变化会改变轨道离心率。典型通信卫星 (m=3000kg, A=50m²) 的 SRP 年度 ΔV 需求约 1〜3 m/s，是东西向保持 (E-W) 的主要干扰之一。与太阳月球三体摄动南北向保持 (N-S) 的 50 m/s/年相比较小，但直接关系到推进剂寿命，必须在轨道计划中考虑。

太阳帆用大面积薄膜 (典型为 7.5 μm 聚酰亚胺上蒸镀铝) 反射太阳光，以 SRP 为主推进源。不消耗推进剂，单位 ΔV 的质量增分为零，长期运行可到达深空。JAXA 的 IKAROS (2010) 为 14m 方形、约 310kg，实现 β≈1.4e-4，NASA NEA Scout、LightSail 2 也采用相同原理。通过改变帆面角度可产生切向力分量，能自由改变轨道半径。

实际应用

GEO 静止卫星的东西向保持：通信和广播卫星需要相对地球保持在同一经度。SRP 导致的日周变化加速度使轨道离心率逐渐增大，卫星位置在东西方向摇晃，因此需要每年进行几次推进器喷射进行补正。设计时会将 A/m 比 (面积质量比) 纳入 SRP 分析模型，据此确定推进剂装载量。在本工具中选择 GEO、增大面积、减小质量，会看到累积 ΔV 急剧增加。

拉格朗日点任务的轨道确定：太阳-地球系统的 L1 (SOHO、WIND、ACE) 和 L2 (Gaia、JWST、Euclid) 是本质上不稳定的平衡点，即使是微弱的 SRP 也会导致利萨如轨道漂移。实际任务中需要高精度分析 SRP 合力与质量中心偏离产生的转矩两个因素，每隔几个月进行一次保持 (SK) 喷射来维持利萨如轨道。JWST 的 SK 推进剂枯竭实际上决定了卫星的使用寿命上限。

太阳帆的深空探测：JAXA IKAROS (2010、金星方向、β≈1.4e-4)、Planetary Society LightSail 2 (2019、绕地球、验证了轨道高度上升)、NASA NEA Scout (2022、小行星接近)，以及正在进行中的 Solar Cruiser、Heliogyro 等中大型帆计划。对于太阳极地轨道、星际边界、长周期彗星追踪等化学推进难以到达的任务，以 SRP 为主推力的探测器正成为现实可行的选择。

精密轨道确定 (POD) 的高精度模型：GPS、Galileo 等定位卫星和地球观测卫星 (Sentinel、GRACE) 要求厘米级的轨道精度。通过高精度 SRP 模型 (Box-wing 法、CODE 模型等)，与地球反照和红外辐射压力一起纳入运动方程，实现了厘米级精密轨道重构。SRP 不仅是简单的摄动，而是决定测地、太空天气、暗物质探测等精度的关键技术。

常见误解与注意事项

最常见的误解是「SRP 微弱可以忽略」。确实 1 AU 处的压力仅为 4.54 μN/m²，极其微小。但在真空中没有摩擦和大气阻力，微小的加速度经过时间积分会产生显著的 ΔV。默认参数 (1000 kg、20 m²) 在一年内会累积约 3.7 m/s，五年运行达到 18 m/s。这虽然不及月球三体摄动的南北向 ΔV (50 m/s/年)，但推进剂计划中绝不能忽视。「微弱可忽视」的说法仅对短期估算 (时间到日级) 成立。

其次，「反射率 ρ 为 1 时 SRP 力翻倍」这个简化。本工具的公式 F = (1+ρ)·P·A 基于完美镜面反射的模型。实际航天器表面 (太阳电池、隔热层、白色涂层、铝箔等) 是镜面反射、漫反射、吸收的混合。详细分析使用「Box-wing 模型」，对每个表面贴片分别设定镜面反射率 ρ_s、漫反射率 ρ_d、吸收率 α，或用更细致的 ray-tracing 计算。本工具数值仅为概算，实际运用需要超越简单球形反射模型的方法。

最后，「地球阴影 (eclipse) 未考虑」。LEO 卫星每圈 (约 90 分钟) 中最多 36 分钟在地球阴影中，此时 SRP 为零。GEO 在春秋分前后 45 天也会出现最长 72 分钟的食。本工具假设 100% 日照，给出了上限值。LEO 实际的年累积 ΔV 约为本工具值的 60〜65%。另外卫星自身的自转、帆角度变化、地球反照 (约 0.3) 的影响也被省略。详细分析需配合 SPICE 等轨道历算库进行时间积分。

太阳辐射压 SRP 模拟器 — 航天器轨道摄动与帆推进

太阳辐射压 (SRP) 与航天器摄动

常见问题

实际应用

常见误解与注意事项

使用指南

具体计算例

实际运用要点