Phase-Gated Constrained Reinforcement Learning for End-to-End Capture and Detumbling of Non-Cooperative Targets
Authors: Álvaro Belmonte-Baeza, Celia Redondo-Verdú, José L. Ramón, Miguel Cazorla, Leonard Felicetti, Jorge Pomares
Organizations: Department of Computer Science and Artificial Intelligence, University of Alicante, Spain. · Department of Physics, Systems Engineering and Signal Theory, University of Alicante, Spain. · Centre for Space Systems, Cranfield University, Cranfield, UK.
Active debris removal requires a servicer spacecraft to approach a non-cooperative tumbling object, capture it, and arrest its rotation. This sequence is conventionally divided among phase-specific controllers, each depending on knowledge of the target's inertia, a dedicated grappling fixture, or a hand-designed switch at the contact-rich seam where such pipelines are most fragile. We show that a single constrained reinforcement learning policy can execute the entire mission, driving a dual-arm free-flying servicer from free-space closing through friction-based fixture-free capture to full detumbling of a target whose mass and inertia are unknown and whose state is observed only through noisy estimates. The enabling mechanism is a smooth proximity gate applied at three loci of a constrained Markov decision process. In the reward it induces the phase structure rather than imposing it; in the action space it forecloses a degenerate grab-and-pull exploit as a hard guarantee; and in the constraint set it resolves a normalization failure of constraint-based learning, in which a constraint spanning phases that load it unequally is rendered nearly inert in the phase where precision matters most. Both the failure mode and the remedy apply to any multi-phase problem of this kind. Over 1000 randomized trials the policy succeeds in 99.5% while holding its contact-force, actuator, base-reaction and collision limits throughout, transfers zero-shot to unseen geometries and centre-of-mass offsets, and extrapolates across a tenfold range of target mass.
Figures & tables
Fig. 1 : The three phases of the ADR mission . The servicer closes the separation under thruster force fthr and reaction torque τrw with the arms stowed; the arms are then actuated by joint torques τ to press two disk-shaped end effectors on opposing faces of the tumbling target, and the resulting contact wrenches fc,L , fc,R dissipate its angular momentum ωt . The phases are named for exposition only: no phase index is supplied to the controller and no switching condition is evaluated at run time.
Quantity
Symbol
Value
Reference
Relative position
σp
1×10−2m
[ 26 ]
Relative attitude
σθ
1.75×10−2rad
[ 26 , 45 ]
Relative linear vel.
σv
5×10−3ms−1
[ 27 ]
Target angular vel.
σω
1.7×10−3rads−1
[ 27 ]
TABLE I : Relative-state estimation noise (standard deviations).
Fig. 2 : Method overview . Both gates are computed from the state, not from the policy: the proximity gate gp modulates the reward, the manipulator command and the constraint terms, while the contact gate gc cascades onto the stabilization rewards. The reward and constraint machinery acts only during training; at deployment the phase structure it induced is already internal to the policy, which requires no supervision beyond the action gating of Eq. 15 .
Group
Quantity
Symbol
Dim.
Target (noisy)
Relative position, base frame
p~tb
3
Relative orientation
φ(R~tb)
6
Relative linear velocity
v~tb
3
Angular velocity
ω~tb
3
Servicer
Base linear / angular velocity
vbb,ωbb
6
Joint positions
q
14
TABLE II : Observation space ( 93 dimensions). Target quantities are corrupted by the estimation noise of ( 7 ); the remainder are measured onboard. φ(⋅) denotes the 6 D rotation encoding [ 46 ] .
Term
Error x
σ
Gate
w
Approach
Closing distance
d−d⋆
0.25 m
—
5.0
Arm stillness
q−q0
0.3 rad
1−gp
5.0
Capture
EE pose tracking a
ep ; eR
0.3 ; 0.5
gp
10.0
EE centering
ei
0.10 m
gp
2.0
TABLE III : Reward terms . Each term is κ(x;σ) on the listed error unless marked; gp is the proximity gate ( 18 ) and gc the contact gate ( 22 ). All terms are scaled by the control timestep.
Constraint
Limit
Gate
pmax
Mission
Target linear velocity
0.01m/s
1c
0.25
Target angular velocity
0.02rad/s
1c
0.25
Base linear velocity (loose)
0.08m/s
—
0.25
Base angular velocity (loose)
0.10rad/s
—
0.25
Base linear velocity (tight)
0.01m/s
1c
0.50
TABLE IV : Constraint set . 1c denotes the end-effector contact indicator. pmax is the maximum termination probability of ( 11 ).
Stage
Rate [95% CI] (%)
Approach complete
100.0 [99.6, 100.0]
Grasp established
100.0 [99.6, 100.0]
Detumble complete (kinematic)
99.5 [98.8, 99.8]
Full success (grasp held)
99.5 [98.8, 99.8]
TABLE V : Phase completion rates over n=1000 randomized trials . Intervals are 95% Wilson score intervals.
Quantity
Value
Residual target linear velocity ∥vt∥
0.0020±0.0009 m/s
Residual target angular velocity ∥ωt∥
0.0036±0.0024 rad/s
Residual base linear velocity ∥vb∥
0.0028±0.0012 m/s
Residual base angular velocity ∥ωb∥
0.0040±0.0014 rad/s
Time to full stabilization (mean)
9.7±8.9 s
Time to full stabilization (median)
5.9 s
TABLE VI : Residual state, timing, and actuation effort at episode end , mean ± std over the n=995 successful trials.
Fig. 3 : Time to full stabilization against initial target angular velocity over 1000 trials ; the five failures are marked at the episode horizon. The side panel shows the distribution of the initial tumbling magnitude.
Constraint
Symbol
Violation (%)
Mission
Target linear velocity
cvt
2.25
Target angular velocity
cωt
3.55
Base linear velocity (loose)
cvb
0.00
Base angular velocity (loose)
cωb
0.01
Base linear velocity (tight)
cvbc
1.43
TABLE VII : Constraint satisfaction over the 1000 trials of section VI-B , as the percentage of control steps at which each limit is exceeded.
Sample ID
∥ωt∥0
min∥ωt∥
tmin
∥ωt∥f
mt
(rad/s)
(rad/s)
(s)
(rad/s)
(kg)
38
0.329
0.0041
28.9
0.0289
30.9
550
0.314
0.0016
12.7
0.0584
33.4
578
0.358
0.0018
8.2
0.0192
40.5
662
0.305
0.0008
2.9
0.0259
44.4
825
0.331
0.0009
6.4
0.0445
26.6
TABLE VIII : The five failed trials of section VI-B . min∥ωt∥ is the lowest angular velocity reached, at tmin ; ∥ωt∥f is the value at the end of the episode. The tolerance is εω=0.03 rad/s.
Fig. 4 : Emergent phase structure in one representative episode . Left column, top to bottom: servicer–target distance with the proximity gate gp ; target and base angular velocity with the success tolerance εω ; manipulator joint torques (solid: left arm; dashed: right arm). Right column, top to bottom: end-effector contact forces with the contact gate gc ; target and base linear velocity; commanded thruster force and reaction torque magnitudes. Vertical lines mark the behaviorally-defined phase boundaries — the onset of bilateral contact and the onset of sustained tumble reduction — rather than any gate value.
Performance
Mission constraints
Actuation
Safety
Scenario
Succ.
∥ωt∥f
∥vt∥f
tstab
cvt
cωt
cvb
cωb
cvbc
cωbc
cf
cτ
cq
cF
ccol
(%)
(rad/s)
(m/s)
(s)
(% of steps)
Nominal
99.5
0.0036
0.0020
9.7
2.25
3.55
0.00
0.01
1.43
2.99
0.12
0.06
0.02
0.27
0.00
Unknown geometry (zero-shot)
Cylinder
99.2
0.0036
0.0021
6.5
2.24
3.83
0.13
0.11
3.00
4.51
0.19
0.06
0.03
0.22
0.00
Cuboid
99.7
0.0039
0.0021
7.7
2.22
4.08
0.04
0.04
2.74
5.03
0.22
0.08
0.03
0.25
0.00
TABLE IX : Generalization and robustness . Left: success rate, residual target velocities, and settling time, the latter three as means over successful trials. Right: constraint satisfaction, as the percentage of control steps at which each limit is exceeded; symbols are defined in Table IV . The geometry and inertia scenarios are zero-shot; tumbling rates are per-axis bounds.
Fig. 5 : Zero-shot capture of unseen geometries . The policy is unchanged from the cube-trained one.
Fig. 6 : Success rate versus the initial tumbling magnitude . Error bars are 95% Wilson intervals; the shaded region marks magnitudes attainable during training.
R
A
C
Success
Grasp loss
∥ωt∥f
∥ωb∥f
Δv
(%)
(%)
(rad/s)
(rad/s)
(m/s)
✓
✓
dec.
99.5
0.0
0.0036
0.0040
0.045
✓
–
dec.
97.3
0.0
0.0087
0.0058
0.189
–
✓
dec.
92.1
0.0
0.0154
0.0115
0.727
–
–
dec.
54.6
10.8
0.0433
0.0296
0.707
✓
✓
int.
23.8
3.5
0.0311
0.0149
0.254
TABLE X : Ablations . R and A denote the reward and action-space gating of sections V-D and V-E ; C denotes the constraint formulation of section V-F , either the decomposed loose and tight pair ( dec. ) or the single constraint whose limit is interpolated between the same two values according to the contact signal ( int. ). Grasp loss is the percentage of trials in which contact is established and subsequently lost. Residuals are means over successful trials.
Method
Target state
Inertial params.
Phase structure
Constraints in loop
Fixture- free
Reported outcome
Friction detumbling [ 32 ]
full
required
switched
no
yes
20∘ /s halted in 30 s, base <0.15∘
Adaptive detumbling [ 33 ]
full
online ID
post-capture
no
yes
inertia uncertainty rejected
Robust SMC [ 13 ]
full
bounded
switched
no
no
chatter-free second-order SMC
Time-optimal [ 11 ]
full
bounded
plan+track
bounds
no
detumble and pose in one action
Optimal concurrent [ 34 ]
full
online ID
plan+track
actuator +motion
yes
finite-time, low chatter
Position/force [ 35 ]
partial
none
switched SMC+RL
joints
yes
PPO force loop, unknown target
TABLE XI : Comparison with representative model-based and hybrid capture/detumbling methods, along the assumptions each requires. Reported outcomes are quoted under each method’s own conditions and are not comparable across rows.
ETSI Aeronáutica y del Espacio, Universidad Politécnica de Madrid, Madrid, Spain · Indra Sistemas, Madrid, Spain · Department of Systems and Industrial Engineering, University of Arizona, Tucson, USA +2
School of Automation Science and Electrical Engineering, Beihang University, 100191, Beijing, China · School of Mechanical and Aerospace Engineering, Nanyang Technological University, 639798, Singapore