Кейс: солнечный парус изменяет ориентацию и траекторию под действием солнечного ветра и фотонного давления — какие силы действуют на тонкую отражающую пластину и как оптимизировать манёвры
Кратко и по существу — какие силы действуют и как оптимизировать манёвры. 1) Основные силы на тонкую отражающую пластину - Давление фотонов (радиационное давление). Для идеального зеркала с площадью AAA, интенсивностью солнечного потока III и углом между нормалью пластины и направлением солнечного света α\alphaα: Prad=2Iccos2α
P_{\text{rad}}=\frac{2I}{c}\cos^2\alpha Prad=c2Icos2α
сила на всю панель: F⃗rad=2IAccos2α n^
\vec F_{\text{rad}}=\frac{2I A}{c}\cos^2\alpha\,\hat n Frad=c2IAcos2αn^
радиальная и поперечная компоненты: Fr=2IAccos3α,Ft=2IAccos2αsinα.
F_r=\frac{2I A}{c}\cos^3\alpha,\qquad F_t=\frac{2I A}{c}\cos^2\alpha\sin\alpha. Fr=c2IAcos3α,Ft=c2IAcos2αsinα.
Для частично поглощающей/диффузно отражающей поверхности формула усложняется (коэффициенты отражения/поглощения). - Солнечный ветер (динамическое давление частиц): порядок величины около 10−910^{-9}10−9–10−8 Па10^{-8}\ \text{Па}10−8Па на 1 AU, обычно на 2–3 порядка меньше фотонного давления и потому часто пренебрегаем. - Гравитационная сила от центрального тела: F⃗g=−μmr2r^\vec F_g=-\dfrac{\mu m}{r^2}\hat rFg=−r2μmr^. - Торки и моменты: при несоосности центра давления и центра масс возникает момент τ⃗=(r⃗CP−r⃗CM)×F⃗rad.
\vec\tau=(\vec r_{CP}-\vec r_{CM})\times\vec F_{\text{rad}}. τ=(rCP−rCM)×Frad. - Дополнительно: силы из-за термального перераспределения излучения (реэмиссия), электростатические и лоренцевы силы при зарядке в магнитном поле (обычно малы), атмосферный дрейф у низких орбит, удары микрометеоритов. 2) Удобные параметры и числовые ориентиры - Солнечная постоянная на 1 AU: I0≈1361 W/m2I_0\approx1361\ \text{W/m}^2I0≈1361W/m2. - Давление идеального зеркала на 1 AU: P0≈2I0c≈9.1×10−6 PaP_0\approx\frac{2I_0}{c}\approx9.1\times10^{-6}\ \text{Pa}P0≈c2I0≈9.1×10−6Pa (около 9 μN/m29\ \mu\text{N/m}^29μN/m2). - Характерное ускорение через массовую плотность (areal mass) σ=m/A\sigma=m/Aσ=m/A: a(α)=2Icσcos2α.
a(\alpha)=\frac{2I}{c\sigma}\cos^2\alpha. a(α)=cσ2Icos2α.
Пример: при σ=10 g/m2=0.01 kg/m2\sigma=10\ \text{g/m}^2=0.01\ \text{kg/m}^2σ=10g/m2=0.01kg/m2 максимальное a(0)≈9.1×10−4 m/s2a(0)\approx9.1\times10^{-4}\ \text{m/s}^2a(0)≈9.1×10−4m/s2 (~0.9 mm/s²). 3) Как оптимизировать манёвры (практические принципы) - Уменьшайте areal mass σ\sigmaσ и увеличивайте отражательность — основное средство повышения тяги. - Для увеличения орбитальной энергии (спиральный вывод наружу) нужен про-градуальный (поперечный) компонент тяги; для этого ставят платину под наклоном. Транзитный компонент максимален при угле, где функция f(α)=cos2αsinαf(\alpha)=\cos^2\alpha\sin\alphaf(α)=cos2αsinα достигает максимума: ddα(cos2αsinα)=0 ⇒ cosα=23 ⇒ α≈35.26∘.
\frac{d}{d\alpha}(\cos^2\alpha\sin\alpha)=0\ \Rightarrow\ \cos\alpha=\sqrt{\tfrac{2}{3}}\ \Rightarrow\ \alpha\approx35.26^\circ. dαd(cos2αsinα)=0⇒cosα=32⇒α≈35.26∘.
То есть для максимального поперечного ускорения держите наклон около 35∘35^\circ35∘. - Для быстрого увеличения радиальной компоненты (выброс наружу, коррекция перицентра) — максимально лицом к Солнцу (α≈0\alpha\approx0α≈0). - Для изменения наклонения орбиты нужен нормальный (вне орбитальной плоскости) компонент — наклоняйте плоскость так, чтобы часть силы была в нормальном направлении. - Управление ориентацией: реакционные колёса, управляющие планки (vanes), перемещение груза для смещения центра масс, локальная модуляция отражающей способности. Учитывайте устойчивость центра давления относительно центра масс. - Алгоритмы оптимизации траектории: оптимизация с постоянным конусным углом (simple fixed-cone law) даёт простую реализацию; для минимального времени/топлива применять формулировки минимального времени с непрерывной тягой (метод Понтрягина, численные оптимизаторы). Часто оптимально — комбинировать эпохи «лицом к Солнцу» и «наклон» (bang–bang или поконные переключения) в зависимости от цели. 4) Практические рекомендации - Проектируйте форму и расположение отражателей так, чтобы центр давления максимально близко совпадал с центром масс или предусматривать активное управление моментом. - Для межпланетных трансфертов использовать угол ~35° для эффективного изменения орбитального момента, лицом к Солнцу — для быстрого набора энергии. - Включайте модели затенения, деградации отражающей поверхности и случайные возмущения (солнечная активность) в оптимизатор траектории. Если нужно — могу дать формулы для частично поглощающей/диффузной поверхности, вывести уравнения Гаусса для малых возмущений орбит при заданных Fr,FtF_r,F_tFr,Ft или поставить задачу оптимизации (минимальное время/минимальная масса) для конкретных исходных данных.
1) Основные силы на тонкую отражающую пластину
- Давление фотонов (радиационное давление). Для идеального зеркала с площадью AAA, интенсивностью солнечного потока III и углом между нормалью пластины и направлением солнечного света α\alphaα:
Prad=2Iccos2α P_{\text{rad}}=\frac{2I}{c}\cos^2\alpha
Prad =c2I cos2α сила на всю панель:
F⃗rad=2IAccos2α n^ \vec F_{\text{rad}}=\frac{2I A}{c}\cos^2\alpha\,\hat n
Frad =c2IA cos2αn^ радиальная и поперечная компоненты:
Fr=2IAccos3α,Ft=2IAccos2αsinα. F_r=\frac{2I A}{c}\cos^3\alpha,\qquad F_t=\frac{2I A}{c}\cos^2\alpha\sin\alpha.
Fr =c2IA cos3α,Ft =c2IA cos2αsinα. Для частично поглощающей/диффузно отражающей поверхности формула усложняется (коэффициенты отражения/поглощения).
- Солнечный ветер (динамическое давление частиц): порядок величины около 10−910^{-9}10−9–10−8 Па10^{-8}\ \text{Па}10−8 Па на 1 AU, обычно на 2–3 порядка меньше фотонного давления и потому часто пренебрегаем.
- Гравитационная сила от центрального тела: F⃗g=−μmr2r^\vec F_g=-\dfrac{\mu m}{r^2}\hat rFg =−r2μm r^.
- Торки и моменты: при несоосности центра давления и центра масс возникает момент
τ⃗=(r⃗CP−r⃗CM)×F⃗rad. \vec\tau=(\vec r_{CP}-\vec r_{CM})\times\vec F_{\text{rad}}.
τ=(rCP −rCM )×Frad .
- Дополнительно: силы из-за термального перераспределения излучения (реэмиссия), электростатические и лоренцевы силы при зарядке в магнитном поле (обычно малы), атмосферный дрейф у низких орбит, удары микрометеоритов.
2) Удобные параметры и числовые ориентиры
- Солнечная постоянная на 1 AU: I0≈1361 W/m2I_0\approx1361\ \text{W/m}^2I0 ≈1361 W/m2.
- Давление идеального зеркала на 1 AU: P0≈2I0c≈9.1×10−6 PaP_0\approx\frac{2I_0}{c}\approx9.1\times10^{-6}\ \text{Pa}P0 ≈c2I0 ≈9.1×10−6 Pa (около 9 μN/m29\ \mu\text{N/m}^29 μN/m2).
- Характерное ускорение через массовую плотность (areal mass) σ=m/A\sigma=m/Aσ=m/A:
a(α)=2Icσcos2α. a(\alpha)=\frac{2I}{c\sigma}\cos^2\alpha.
a(α)=cσ2I cos2α. Пример: при σ=10 g/m2=0.01 kg/m2\sigma=10\ \text{g/m}^2=0.01\ \text{kg/m}^2σ=10 g/m2=0.01 kg/m2 максимальное a(0)≈9.1×10−4 m/s2a(0)\approx9.1\times10^{-4}\ \text{m/s}^2a(0)≈9.1×10−4 m/s2 (~0.9 mm/s²).
3) Как оптимизировать манёвры (практические принципы)
- Уменьшайте areal mass σ\sigmaσ и увеличивайте отражательность — основное средство повышения тяги.
- Для увеличения орбитальной энергии (спиральный вывод наружу) нужен про-градуальный (поперечный) компонент тяги; для этого ставят платину под наклоном. Транзитный компонент максимален при угле, где функция f(α)=cos2αsinαf(\alpha)=\cos^2\alpha\sin\alphaf(α)=cos2αsinα достигает максимума:
ddα(cos2αsinα)=0 ⇒ cosα=23 ⇒ α≈35.26∘. \frac{d}{d\alpha}(\cos^2\alpha\sin\alpha)=0\ \Rightarrow\ \cos\alpha=\sqrt{\tfrac{2}{3}}\ \Rightarrow\ \alpha\approx35.26^\circ.
dαd (cos2αsinα)=0 ⇒ cosα=32 ⇒ α≈35.26∘. То есть для максимального поперечного ускорения держите наклон около 35∘35^\circ35∘.
- Для быстрого увеличения радиальной компоненты (выброс наружу, коррекция перицентра) — максимально лицом к Солнцу (α≈0\alpha\approx0α≈0).
- Для изменения наклонения орбиты нужен нормальный (вне орбитальной плоскости) компонент — наклоняйте плоскость так, чтобы часть силы была в нормальном направлении.
- Управление ориентацией: реакционные колёса, управляющие планки (vanes), перемещение груза для смещения центра масс, локальная модуляция отражающей способности. Учитывайте устойчивость центра давления относительно центра масс.
- Алгоритмы оптимизации траектории: оптимизация с постоянным конусным углом (simple fixed-cone law) даёт простую реализацию; для минимального времени/топлива применять формулировки минимального времени с непрерывной тягой (метод Понтрягина, численные оптимизаторы). Часто оптимально — комбинировать эпохи «лицом к Солнцу» и «наклон» (bang–bang или поконные переключения) в зависимости от цели.
4) Практические рекомендации
- Проектируйте форму и расположение отражателей так, чтобы центр давления максимально близко совпадал с центром масс или предусматривать активное управление моментом.
- Для межпланетных трансфертов использовать угол ~35° для эффективного изменения орбитального момента, лицом к Солнцу — для быстрого набора энергии.
- Включайте модели затенения, деградации отражающей поверхности и случайные возмущения (солнечная активность) в оптимизатор траектории.
Если нужно — могу дать формулы для частично поглощающей/диффузной поверхности, вывести уравнения Гаусса для малых возмущений орбит при заданных Fr,FtF_r,F_tFr ,Ft или поставить задачу оптимизации (минимальное время/минимальная масса) для конкретных исходных данных.