REVIEW 4 major objections 6 minor 2 cited by
Conjugate momentum based thruster force estimate in dynamic multimodal robot
T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read A conjugate-momentum observer estimates thruster forces on a legged-aerial biped to within 5–19 percent error.
desk verdict Standard momentum observer, honestly applied to a new platform; simulation-only and no sensitivity analysis, but a plausible niche contribution. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the conjugate momentum observer, a standard momentum observer built on the generalized momentum $p = M\dot q$. The estimated thruster force $r$ evolves as $\dot r = K_0(B_t u_t - r)$, making $r$ a low-pass filter of the true thruster force with gain $K_0$; no joint accelerations or inertia-matrix inversions are needed. Because the observer equation needs ground reaction forces, the paper pairs it with a contact-constraint model that enforces zero stance-foot acceleration, $J_c \ddot q = \dot J_c \dot q$, to compute the ground reaction force via a Moore-Penrose pseudo-inverse. The controller used for the simulated walking is built on a variable-length inverted pendulum (VLIP) reduced-order model, with thrusters stabilizing roll and yaw, while the full simulation uses an Euler-Lagrangian model with compliant ground contact.
What would settle it
Measure the actual thruster force with a thrust stand while Harpy walks with thrusters on hardware and compare it with the observer's output: if the normalized RMSE grows well beyond the 0.05–0.19 range under realistic model mismatch, joint friction, or foot slip, the central claim fails. Alternatively, instrument the foot with a force plate to check whether the stance-foot acceleration is truly zero while the constraint model is active.
Extended reading notes
Core claim
Using the generalized momentum $p = M\dot q$, the observer defines an estimated generalized thruster force $r$ driven by $\dot r = K_0(B_t u_t - r)$, so that $r$ is a low-pass filtered version of the actual thruster force. Under ideal conditions, $M$ and $\beta$ are assumed exactly known. The body-frame thruster force is recovered with the Jacobian pseudo-inverse $\hat u_t = [J_t^\top J_t]^\dagger J_t^\top r$. When ground reaction forces are taken from the compliant ground model, the estimated generalized forces and torques match the actual values closely. When they are instead estimated through the contact constraint $J_c \ddot q = \dot J_c \dot q$, the estimates still track but with larger error, and the paper reports NRMSE values of 0.1156 for $F_x$, 0.0492 for $F_y$, 0.1946 for $F_z$, 0.1102 for $\tau_x$, 0.1342 for $\tau_y$, and 0.1239 for $\tau_z$. The vertical force $F_z$ has the largest error, which the authors attribute to the constraint model.
Load-bearing premise
The estimator assumes the robot's mass and bias terms are known exactly, and that the stance foot never accelerates or slips while the contact constraint is enforced; onboard hardware will violate both assumptions to some degree.
Editorial extensions
If this is right
- Controllers such as MPC and QP can use the estimated thruster forces directly, replacing thrust-stand calibration that misses battery voltage and other working conditions.
- A contact-constraint-based ground reaction force estimate removes the need for foot force sensors, at the cost of higher estimation error, with $F_z$ the least accurate.
- Because the estimator is a low-pass filter, its accuracy is governed by the quality of the ground reaction force information; better terrain knowledge yields better thrust estimates.
- The double-support phase makes $J_c M^{-1} J_c^\top$ rank-deficient, so ground force estimates are inherently less reliable during that phase.
- A second-order filter is the paper's stated next step to reduce the filter lag visible in the y-direction torque estimate.
Reading between the lines
- On hardware, the ideal-model assumption will be violated by joint friction and model error; the paper gives no sensitivity analysis, so the first test is whether the observer gain can absorb these errors without amplifying noise.
- The rank deficiency in double support suggests a practical deployment would gate the estimator by contact state or blend the constraint-model ground reaction force with any available foot force sensing.
- The same observer structure could be transferred to other thruster-augmented legged robots, since it only requires the generalized momentum, contact Jacobians, and a source of ground reaction force; the equations do not depend on Harpy's specific kinematics.
- A direct experimental check would be to compare the observer's thrust estimate against a thrust stand while battery voltage drops during flight, quantifying how much of the 5–19% error persists outside simulation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a generalized-momentum (conjugate momentum) observer to estimate thruster forces on the Harpy bipedal robot, a platform that combines legged locomotion with thruster actuation. The observer is derived in Section III and is evaluated in simulation during thruster-assisted walking. Two variants are presented: one that uses ground reaction forces from the compliant ground model directly, and one that estimates ground reaction forces from a contact constraint model. The reported results are normalized RMSE values between 0.05 and 0.19 for the generalized force/torque components in the constraint-model case (Table I). The paper claims that the observer accurately tracks thruster forces and that the method is suitable for onboard estimation where thrust-stand characterization is insufficient.
Significance. If the result holds, the paper offers a sensorless thruster-force estimation scheme for a multimodal legged-aerial robot, which is relevant for downstream MPC/QP control. A strength is that the paper builds on the well-established momentum-observer framework, avoiding acceleration estimation and matrix inversion. Another strength is that the simulation includes full-dynamics evaluation and reports quantitative errors, not just qualitative plots. However, the contribution is incremental relative to the existing momentum-observer literature, and the reported accuracy is obtained under strong ideal-model assumptions. The paper does not yet demonstrate robustness to model error, foot slip, or state-estimation delay, which are central to the stated motivation of onboard estimation under real operating conditions.
major comments (4)
- [III, Eq. (13)] The observer equation is not self-contained: the equality rdot = K0(Bt ut - r) = K0(pdot - pdot_hat) introduces pdot_hat without a definition or an explicit relation to the state and known inputs. In the standard momentum-observer derivation, one obtains rdot = K0(pdot - beta - r - Bg lambda - Bj uj), and the equivalence to K0(pdot - pdot_hat) must be shown. Please define pdot_hat and state the exact residual dynamics used in the simulator.
- [III-A] The zero-foot-acceleration constraint is stated as Jc qddot = Jdot qdot, but differentiating Jc qdot = 0 gives Jc qddot = -Jdot qdot. Equation (16) contains the correct minus sign, so this is an inconsistency in the text rather than in the implemented formula; please correct the sentence to avoid confusion.
- [IV-B] The constraint-model estimator is a feedback loop: lambda in Eq. (16) is computed from the residual r, and r is integrated in Eq. (14) using that lambda. No passivity, small-gain, or error-propagation analysis is given for this loop, and the only mitigation offered is raising K0 from 25 to 3000. Because Table I is generated under the ideal assumptions M-hat = M, beta-hat = beta, and no foot slip, the reported 0.05-0.19 NRMSE is not yet demonstrated as robust to the model error and terrain variation expected on hardware. Please add a sensitivity analysis (e.g., parameter perturbations, friction/slip, state-estimation delay) or a bounded-error argument.
- [IV, Table I] The stated contribution is a comparison of estimation with and without terrain knowledge, but Table I reports NRMSE only for the constraint-model case. The ground-model case (Fig. 6) is only qualitative. Please report the same NRMSE metric for both cases so the comparison is quantitative.
minor comments (6)
- [Abstract and Section I] There are numerous grammatical errors and typos, such as 'such our state-of-the-art Harpy platform' and 'we can characterize thruster force using a thrust stand but it generally does not account for working conditions.' Please proofread the manuscript.
- [II] Harpy's height is stated as 600 cm, which is presumably a typo for 60 cm. Please correct the unit.
- [II-C, Eq. (9)] The composition of u_t from u_t,c and u_t,L/u_t,R is unclear; specify the dimensions and the frames in which the components are expressed.
- [IV-B] The observer gains are reported for the two cases, but no tuning procedure or criterion is given. State how the gains were selected and whether the results are sensitive to their values.
- [IV-B, Table I] The 'Normalized RMSE' metric is not defined. Specify the normalization denominator (e.g., range, mean, standard deviation) so that the values in Table I are interpretable.
- [III] The term 'conjugate momentum' is used, but the derivation is the standard generalized momentum observer. Consider clarifying the terminology or providing a reference that uses this name.
Circularity Check
No significant circularity: the thruster-force estimator is a standard momentum-based observer, and its validation against simulated ground truth does not reduce to its own inputs by construction.
full rationale
The estimation scheme is the classical generalized-momentum observer: equation (13) defines the residual dynamics as dot-r = K0 (Bt ut - r) under the stated ideal-model assumption M-hat = M and beta-hat = beta, so the estimated thruster force is a low-pass filtered version of the actual thruster force rather than an algebraic renaming of it. The observer gains are hand-tuned and are not fitted to the NRMSE values reported in Table I, and the reported errors are comparisons against simulated ground truth, not a fitted parameter renamed as a prediction. The constraint-model branch (equation (16)) does feed the estimated residual r back into the GRF estimate, and that GRF estimate is then used in equation (14); this creates a coupled observer loop, but it is a genuine feedback loop with error dynamics, not a definitional identity, and the paper explicitly acknowledges that errors in r propagate into the GRF estimate and that the observer gain must be increased to mitigate this. The remaining limitations (exact model knowledge, zero stance-foot acceleration, no slip) are stated assumptions and affect external validity, but they are not circular reductions. Self-citations to prior Harpy hardware and control work are contextual and do not carry the load-bearing derivation, which rests on the standard momentum-observer literature and the paper's own dynamics equations.
Assumptions & free parameters
free parameters (3)
- Observer gain K0 =
K0 = diag(1,1,1,1,25,25,25,25,25,25) for ground-model case; K0 = diag(1,1,1,1,800,1200,60,3000,800,500) for…
- Ground contact model parameters =
µs=0.8, µc=0.64, µv=0.8, kg,p=8000, kg,d=268
- Robot mass, inertia, and link dimensions =
mB=2, mH=mK=0.5, IB=1e-3, IH=IK=1e-4; l1=[0,0.1,-0.1], l2=[0,0.5,0], l3=[0,0,-0.3], l4=[0,0.1,0]
assumptions (3)
- domain assumption The observer dynamic model is exact: M-hat = M and beta-hat = beta.
- domain assumption The stance foot acceleration is zero and there is no slip during contact.
- domain assumption Thrusters can produce force in any direction and their dynamics are ignored.
Cite this review
Pith. "Pith review of Conjugate momentum based thruster force estimate in dynamic multimodal robot." pith.science (2026). https://pith.science/paper/VKT33ZQZ
@misc{pith2026241114596,
author = {Pith},
title = {Pith review of: Conjugate momentum based thruster force estimate in dynamic multimodal robot},
year = {2026},
howpublished = {\url{https://pith.science/paper/VKT33ZQZ}},
note = {Machine review of arXiv:2411.14596}
}
read the original abstract
In a multi-modal system which combines thruster and legged locomotion such our state-of-the-art Harpy platform to perform dynamic locomotion. Therefore, it is very important to have a proper estimate of Thruster force. Harpy is a bipedal robot capable of legged-aerial locomotion using its legs and thrusters attached to its main frame. we can characterize thruster force using a thrust stand but it generally does not account for working conditions such as battery voltage. In this study, we present a momentum-based thruster force estimator. One of the key information required to estimate is terrain information. we show estimation results with and without terrain knowledge. In this work, we derive a conjugate momentum thruster force estimator and implement it on a numerical simulator that uses thruster force to perform thruster-assisted walking.
Figures
Figures from the paper (3 more)
Forward citations
Cited by 2 Pith papers
-
Estimation of Aerodynamics Forces in Dynamic Morphing Wing Flight
A physics-based observer and an MLP regression both estimate aerodynamic forces on a morphing-wing robot within about 0.13 N of load-cell measurements.
-
Analysis of Harpy's Constrained Trotting and Jumping Maneuver
The provided manuscript text does not contain the claimed analysis of the Harpy robot, making the abstract's conclusions unverifiable from this document.
Reference graph
Works this paper leans on
-
[1]
Robot Collisions: A Survey on Detection, Isolation, and Identification,
S. Haddadin, A. De Luca, and A. Albu-Sch ¨affer, “Robot Collisions: A Survey on Detection, Isolation, and Identification,” IEEE Trans- actions on Robotics , vol. 33, no. 6, pp. 1292–1312, Dec. 2017
work page 2017
-
[2]
Sensorless Robot Collision Detection and Hybrid Force/Motion Control,
A. de Luca and R. Mattone, “Sensorless Robot Collision Detection and Hybrid Force/Motion Control,” in Proceedings of the 2005 IEEE International Conference on Robotics and Automation , Apr. 2005, pp. 999–1004
work page 2005
-
[3]
Actuator failure detection and isolation using generalized momenta,
A. De Luca and R. Mattone, “Actuator failure detection and isolation using generalized momenta,” in 2003 IEEE International Conference on Robotics and Automation (Cat. No.03CH37422) , vol. 1, Sep. 2003, 634–639 vol.1
work page 2003
-
[4]
Predefined-Time External Force Esti- mation for Legged Robots,
P. Cai, D. Liu, and L. Zhu, “Predefined-Time External Force Esti- mation for Legged Robots,” in Intelligent Robotics and Applications, Singapore: Springer Nature, 2023, pp. 542–552
work page 2023
-
[5]
Sliding Mode Momentum Observers for Estimation of External Torques and Joint Acceleration,
G. Garofalo, N. Mansfeld, J. Jankowski, and C. Ott, “Sliding Mode Momentum Observers for Estimation of External Torques and Joint Acceleration,” in 2019 International Conference on Robotics and Automation (ICRA), May 2019, pp. 6117–6123
work page 2019
-
[6]
On making robots understand safety: Embedding injury knowledge into control - Sami Haddadin, Simon Haddadin, Augusto Khoury, Tim Rokahr, Sven Parusel, Rainer Burgkart, Antonio Bicchi, Alin Albu-Sch ¨affer, 2012 . [Online]. Available: https : / / journals . sagepub . com / doi / abs / 10 . 1177 / 0278364912462256 ? casa _ token = zT9wDmGIPrMAAAAA : gr3D3cW...
work page 2012
-
[7]
Sensorless Ground Reaction Force Observation With Disturbance Compensation in Heavy-Legged Robots,
S. Liu, Z. Pan, S. Zhou, Z. Niu, and R. Wang, “Sensorless Ground Reaction Force Observation With Disturbance Compensation in Heavy-Legged Robots,” IEEE/ASME Transactions on Mechatronics, pp. 1–12, 2024
work page 2024
-
[8]
Contact Model Fusion for Event-Based Locomotion in Unstructured Terrains,
G. Bledt, P. M. Wensing, S. Ingersoll, and S. Kim, “Contact Model Fusion for Event-Based Locomotion in Unstructured Terrains,” in 2018 IEEE International Conference on Robotics and Automation (ICRA), May 2018, pp. 4399–4406
work page 2018
Show all 25 references
-
[9]
[Online]
Residual-based contacts estimation for humanoid robots — IEEE Conference Publication — IEEE Xplore . [Online]. Available: https : / / ieeexplore . ieee . org / abstract / document/7803308?casa_token=2RhORD0rn4IAAAAA: 3fe9evpFqIEVpAHm23ly55x8MvR9ff _ HwCPHTNpoT10zm5e7xUHYW38kAI...
2024
-
[10]
External Wrench Estimation, Collision Detection, and Reflex Reaction for Flying Robots,
T. Tomi ´c, C. Ott, and S. Haddadin, “External Wrench Estimation, Collision Detection, and Reflex Reaction for Flying Robots,” IEEE Transactions on Robotics, vol. 33, no. 6, pp. 1467–1482, Dec. 2017
2017
-
[11]
Multi-Modal Mobility Morphobot (M4) with appendage repurpos- ing for locomotion plasticity enhancement,
E. Sihite, A. Kalantari, R. Nemovi, A. Ramezani, and M. Gharib, “Multi-Modal Mobility Morphobot (M4) with appendage repurpos- ing for locomotion plasticity enhancement,” Nature Communications, vol. 14, no. 1, p. 3323, Jun. 2023
2023
-
[12]
Efficient Path Planning and Tracking for Multi-Modal Legged- Aerial Locomotion Using Integrated Probabilistic Road Maps (PRM) and Reference Governors (RG),
E. Sihite, B. Mottis, P. Ghanem, A. Ramezani, and M. Gharib, “Efficient Path Planning and Tracking for Multi-Modal Legged- Aerial Locomotion Using Integrated Probabilistic Road Maps (PRM) and Reference Governors (RG),” in 2022 IEEE 61st Conference on Decision and Control (CDC)...
2022
-
[13]
Minimum Time Trajectory Generation for Bounding Flight: Combining Posture Control and Thrust Vectoring,
I. Mandralis, E. Sihite, A. Ramezani, and M. Gharib, “Minimum Time Trajectory Generation for Bounding Flight: Combining Posture Control and Thrust Vectoring,” in2023 European Control Conference (ECC), Jun. 2023, pp. 1–7
2023
-
[14]
A bipedal walking robot that can fly, slackline, and skateboard,
K. Kim, P. Spieler, E.-S. Lupu, A. Ramezani, and S.-J. Chung, “A bipedal walking robot that can fly, slackline, and skateboard,”Science Robotics, vol. 6, no. 59, eabf8136, Oct. 2021
2021
-
[15]
Control of Thruster-Assisted, Bipedal Legged Locomotion of the Harpy Robot,
P. Dangol, E. Sihite, and A. Ramezani, “Control of Thruster-Assisted, Bipedal Legged Locomotion of the Harpy Robot,” Frontiers in Robotics and AI , vol. 8, Dec. 2021
2021
-
[16]
Rough-Terrain Locomotion and Unilateral Contact Force Regula- tions With a Multi-Modal Legged Robot,
K. Liang, E. Sihite, P. Dangol, A. Lessieur, and A. Ramezani, “Rough-Terrain Locomotion and Unilateral Contact Force Regula- tions With a Multi-Modal Legged Robot,” in 2021 American Control Conference (ACC), May 2021, pp. 1762–1769
2021
-
[17]
Optimization-free Ground Contact Force Constraint Satisfaction in Quadrupedal Locomotion,
E. Sihite, P. Dangol, and A. Ramezani, “Optimization-free Ground Contact Force Constraint Satisfaction in Quadrupedal Locomotion,” in 2021 60th IEEE Conference on Decision and Control (CDC), Dec. 2021, pp. 713–719
2021
-
[18]
Generative Design of NU’s Husky Carbon, A Morpho-Functional, Legged Robot,
A. Ramezani, P. Dangol, E. Sihite, A. Lessieur, and P. Kelly, “Generative Design of NU’s Husky Carbon, A Morpho-Functional, Legged Robot,” in 2021 IEEE International Conference on Robotics and Automation (ICRA) , May 2021, pp. 4040–4046
2021
-
[19]
Feedback design for Harpy: A test bed to inspect thruster-assisted legged locomotion,
P. Dangol and A. Ramezani, “Feedback design for Harpy: A test bed to inspect thruster-assisted legged locomotion,” in Unmanned Systems Technology XXII, vol. 11425, SPIE, May 2020, pp. 49–55
2020
-
[20]
A multilayer control for multirotor UA Vs equipped with a servo robot arm,
F. Ruggiero, M. Trujillo, R. Cano, et al. , “A multilayer control for multirotor UA Vs equipped with a servo robot arm,” in 2015 IEEE International Conference on Robotics and Automation (ICRA) , May 2015, pp. 4014–4020
2015
-
[21]
Momentum- Based Extended Kalman Filter for Thrust Estimation on Flying Multibody Robots,
H. A. O. Mohamed, G. Nava, G. L’Erario, et al. , “Momentum- Based Extended Kalman Filter for Thrust Estimation on Flying Multibody Robots,” IEEE Robotics and Automation Letters , vol. 7, no. 1, pp. 526–533, Jan. 2022
2022
-
[22]
Dynamic multimodal locomotion: A quick overview of hardware and control.,
S. Pitroda, “Dynamic multimodal locomotion: A quick overview of hardware and control.,” 2023
2023
-
[23]
Control of Thruster-Assisted, Bipedal Legged Locomotion of the Harpy Robot,
P. Dangol, E. Sihite, and A. Ramezani, “Control of Thruster-Assisted, Bipedal Legged Locomotion of the Harpy Robot,” Frontiers in Robotics and AI , vol. 8, 2021
2021
-
[24]
Capture Point Control in Thruster-Assisted Bipedal Locomotion,
S. Pitroda, A. Bondada, K. Venkatesh, et al., “Capture Point Control in Thruster-Assisted Bipedal Locomotion,” in 2024 IEEE Interna- tional Conference on Advanced Intelligent Mechatronics (AIM) , Jul. 2024, pp. 1139–1144
2024
-
[25]
Col- lision Detection and Safe Reaction with the DLR-III Lightweight Manipulator Arm,
A. De Luca, A. Albu-Schaffer, S. Haddadin, and G. Hirzinger, “Col- lision Detection and Safe Reaction with the DLR-III Lightweight Manipulator Arm,” in 2006 IEEE/RSJ International Conference on Intelligent Robots and Systems , Oct. 2006, pp. 1623–1630
2006
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.