本文作者：天疆说
本站地址：https://cislunarspace.cn
轨道优化中的约束（trajectory constraints）是航天器轨道优化问题里，除动力学方程本身之外、轨迹必须满足的所有附加限制。它们把抽象的工程要求翻译成数学条件：航天器何时必须到达何处、沿途必须避开什么、哪些状态或控制量的边界不可逾越。一条同时满足动力学与全部约束的轨迹称为可行解 。优化算法的任务是在可行解集合中最小化性能指标[200] [260] 。
按作用的位置与方式，约束可分四类：
终端/边界约束（terminal / boundary）：作用在某段的初始或终端时刻，例如交会、入轨、拦截条件。 路径约束（path）：在某段全程每一时刻都必须成立的不等式，例如日食规避、热流密度、推力幅值上界。 状态约束（state）：仅含状态量、不含控制量的路径约束，例如近地高度下界、最大飞行距离。它们在间接法中触发协态跳跃。 事件/相位约束（event / phase）：作用在段内特定时刻或跨段联立变量上，例如系统间转移界面处的相位角匹配。 拼接点处的弧间连续性条件，在直接配点法里常被表达为缺陷约束 ，见 缺陷约束 。弧间匹配的理论框架见 多弧段最优控制 。
含四类约束的单段最优控制问题标准形式：
min u ( ⋅ ) , t f J = Φ [ x ( t 0 ) , x ( t f ) , t f ] + ∫ t 0 t f L d t , \min_{\boldsymbol{u}(\cdot),\,t_f} J = \Phi[\boldsymbol{x}(t_0),\boldsymbol{x}(t_f),t_f] + \int_{t_0}^{t_f} L\,\mathrm{d}t, u ( ⋅ ) , t f min J = Φ [ x ( t 0 ) , x ( t f ) , t f ] + ∫ t 0 t f L d t ,
约束为
x ˙ = f ( x , u , t ) , \dot{\boldsymbol{x}} = \boldsymbol{f}(\boldsymbol{x},\boldsymbol{u},t), x ˙ = f ( x , u , t ) ,
ψ [ x ( t 0 ) , x ( t f ) , t f ] = 0 （终端/边界） , \boldsymbol{\psi}[\boldsymbol{x}(t_0),\boldsymbol{x}(t_f),t_f] = \boldsymbol{0} \quad \text{（终端/边界）}, ψ [ x ( t 0 ) , x ( t f ) , t f ] = 0 （终端 / 边界） ,
g l ≤ g [ x ( t ) , u ( t ) , t ] ≤ g u （路径） , \boldsymbol{g}_l \leq \boldsymbol{g}[\boldsymbol{x}(t),\boldsymbol{u}(t),t] \leq \boldsymbol{g}_u \quad \text{（路径）}, g l ≤ g [ x ( t ) , u ( t ) , t ] ≤ g u （路径） ,
x l ≤ x ( t ) ≤ x u （状态量简单上下界） . \boldsymbol{x}_l \leq \boldsymbol{x}(t) \leq \boldsymbol{x}_u \quad \text{（状态量简单上下界）}. x l ≤ x ( t ) ≤ x u （状态量简单上下界） .
某分量取 g l = g u g_l = g_u g l = g u 积分函数 ∫ q d t \int q\,\mathrm{d}t ∫ q d t
在间接法中，终端约束通过乘子 ν \boldsymbol{\nu} ν λ ( t f ) = ∂ Φ / ∂ x ( t f ) + ( ∂ ψ / ∂ x ( t f ) ) T ν \boldsymbol{\lambda}(t_f) = \partial\Phi/\partial\boldsymbol{x}(t_f) + (\partial\boldsymbol{\psi}/\partial\boldsymbol{x}(t_f))^{T}\boldsymbol{\nu} λ ( t f ) = ∂ Φ/ ∂ x ( t f ) + ( ∂ ψ / ∂ x ( t f ) ) T ν λ ( t 0 ) \boldsymbol{\lambda}(t_0) λ ( t 0 ) 两点边值问题 。
当纯状态不等式 S ( x ) ≤ 0 S(\boldsymbol{x})\leq 0 S ( x ) ≤ 0 [ t 1 , t 2 ] [t_1,t_2] [ t 1 , t 2 ] 边界弧 ）时，标准一阶必要条件需引入非负乘子 η ( t ) ≥ 0 \eta(t)\geq 0 η ( t ) ≥ 0 t i t_i t i 跳跃条件
λ ( t i + ) = λ ( t i − ) − η ( t i ) ∇ S [ x ( t i ) ] . \boldsymbol{\lambda}(t_i^{+}) = \boldsymbol{\lambda}(t_i^{-}) - \eta(t_i)\,\nabla S[\boldsymbol{x}(t_i)]. λ ( t i + ) = λ ( t i − ) − η ( t i ) ∇ S [ x ( t i )] .
对 p p p S S S p p p u \boldsymbol{u} u λ ( p − 1 ) \boldsymbol{\lambda}^{(p-1)} λ ( p − 1 ) λ \boldsymbol{\lambda} λ [678] 。
多圈 Lambert 问题在固定 r 0 \boldsymbol{r}_0 r 0 r f \boldsymbol{r}_f r f T T T N N N [287] 证明，轨道动力学边界约束会重塑允许解流形随 N N N N N N ( r 0 , r f , T ) (\boldsymbol{r}_0,\boldsymbol{r}_f,T) ( r 0 , r f , T ) N N N 打靶法 选初值或为不确定 Lambert 问题建代理模型时必须显式建模。
太阳能电推进任务在遮挡天体的本影/半影锥内无法（或只能降级）推力。Sowell & Taheri[603] 把这一约束表达为推力加速度的平滑熄火因子
Δ prop = T max m δ δ ˉ e α ^ , δ ˉ e = 1 2 [ 1 + tanh ( S e / ρ e ) ] , \boldsymbol{\Delta}_{\text{prop}} = \frac{T_{\max}}{m}\,\delta\,\bar{\delta}_e\,\hat{\boldsymbol{\alpha}}, \quad \bar{\delta}_e = \tfrac{1}{2}\left[1 + \tanh(S_e/\rho_e)\right], Δ prop = m T m a x δ δ ˉ e α ^ , δ ˉ e = 2 1 [ 1 + tanh ( S e / ρ e ) ] ,
其中 S e = a D − a R S − a R B S_e = a_D - a_{RS} - a_{RB} S e = a D − a R S − a R B ρ e \rho_e ρ e δ ˉ e ∈ [ 0 , 1 ] \bar{\delta}_e\in[0,1] δ ˉ e ∈ [ 0 , 1 ] ρ e \rho_e ρ e 限制累计日食时长 ，只需把 ∫ ( 1 − δ ˉ e ) d t \int (1-\bar{\delta}_e)\,\mathrm{d}t ∫ ( 1 − δ ˉ e ) d t
日地与地月两套 CR3BP 系统之间的转移，在几何上由不稳定/稳定流形管在庞加莱截面上的交点定位。在星历模型中，这些交点只在特定的月球相位朝向下才存在。Howell & Kakoi[399] 用 3-1-3 欧拉角序列 ( α , i , β ) (\alpha, i, \beta) ( α , i , β ) i ≈ 5 ∘ i\approx 5^{\circ} i ≈ 5 ∘ ψ \psi ψ ( α , β , ψ ) (\alpha, \beta, \psi) ( α , β , ψ ) α \alpha α β \beta β α \alpha α β \beta β
直接转录 把每条约束翻成 NLP 行：每个配点处的缺陷方程（见 缺陷约束 ），加上节点/中点的路径约束估值。利用稀疏性的 SQP/IP 求解器（SNOPT、IPOPT）可处理 10 4 10^4 1 0 4 10 6 10^6 1 0 6
软约束 / 松弛 把硬等式 X P = X Q \boldsymbol{X}_P = \boldsymbol{X}_Q X P = X Q min ∥ X P − X Q ∥ \min\|\boldsymbol{X}_P - \boldsymbol{X}_Q\| min ∥ X P − X Q ∥ [95] 把此法用于地月 L 1 L_1 L 1 过约束 ，优化器无解。松弛为优化目标里的 Lambert Δ V \Delta V Δ V Δ V \Delta V Δ V
平滑化（双曲正切 / CSC） 把分段定义的约束（日食开关、bang-bang 控制）正则化，使间接积分器与 NLP 雅可比保持光滑。与庞特里亚金极值原理并用时常见 复合光滑控制 。
增广拉格朗日 / 罚函数 把活动不等式约束移入目标函数，见 增广拉格朗日法 。
终端约束主导 打靶法 雅可比的结构。终端盒越紧，收敛域越窄。 路径约束活动状态频繁切换（如低推力螺旋上的日食）是伪谱法网格加密的主要驱动。 多体问题中的段间联立约束，是把读者引向 多弧段最优控制 理论的桥梁。 