Researchers at Beijing Institute of Technology have now tackled this measurement gap head-on. In a study published in Space: Science & Technology, the team proposed an error-state extended Kalman filtering method built directly on multibody integrators, allowing the dynamics of a flexible cable net to be predicted and corrected without the cumbersome model transformations that conventional filtering demands. The work addresses a problem that has long frustrated the multibody dynamics community: the cable-net system is inherently a high-dimensional, nonlinear, strongly rigid-flexible coupled constrained multibody system, and its governing equations take the form of index-3 differential-algebraic equations, or DAEs. These equations mix differential equations with algebraic constraints, and traditional approaches force analysts to convert them into ordinary differential equations and then linearize them before a Kalman filter can touch them—a process that is complex, lacks generality, and often degrades accuracy.

The numerical obstacles run deeper still. Direct integration of index-3 DAEs is notoriously stiff, and conventional numerical integrators converge slowly on such problems. The Jacobian matrix becomes severely ill-conditioned when large time steps are used—its condition number, in fact, grows with the cube of the inverse of the time step—and the algebraic constraints that hold the system together are prone to violation, a phenomenon known as constraint drift. For a state estimator that must run fast enough to guide a capture in real time, these are not academic inconveniences; they are showstoppers. The Beijing Institute of Technology team therefore set out to develop a filtering method that requires no model transformation, can use the multibody DAE directly for state prediction, and can maintain constraint satisfaction throughout the estimation process.

The foundation of the new approach is a carefully constructed multibody model of the recovery system itself. The cable-net recovery sequence unfolds in three phases—following, entering, and capturing—during which trolleys drive the cable motion to track the rocket’s landing point. Controlling the position of the capture container center depends on accurate cable dynamic information, which is exactly what the model must supply. Each arresting cable is discretized into variable-length cable elements using the Arbitrary Lagrangian–Eulerian formulation of the Absolute Nodal Coordinate Formulation, a technique well suited to cables whose length changes as they are paid out and reeled in. The complete system comprises four cables, each driven at both ends by trolleys and subjected to distributed loads such as plume impingement, and the resulting dynamic equations are high-dimensional, nonlinear, and strongly coupled between rigid and flexible modes.

To tame the index-3 DAEs, the researchers applied two complementary remedies. First, they performed a scaling transformation on the equations with respect to dimensionless time and the magnitude of the constraint equations. Second, they combined this with an augmented Lagrangian formulation for index reduction. The combined effect is striking: the Jacobian matrix condition number and the error propagation magnitude become independent of the time step, and the numerical behavior of the constrained system reaches stability comparable to that of ordinary differential equations. This means the integrator no longer punishes the user for choosing the larger time steps that real-time operation demands—a property that proves decisive later in the validation.

On top of this stabilized model sits the filtering framework itself. The method adopts an indirect filtering strategy that decouples state prediction from error covariance update. In the prediction stage, one of two multibody integrators—the generalized-α method or the backward differentiation formula (BDF)—performs temporal integration of the scaled, index-reduced DAEs, outputting generalized coordinates and velocities in a single step, with no conversion to ODEs and no linearization of the original multibody system. Error-state covariance propagation employs a simplified state transition matrix to keep the computational burden manageable. In the correction stage, the posterior estimate of the error state is computed from lidar measurements, and the estimated state is projected back onto the constraint manifold through constraint-violation correction, ensuring that both position and velocity estimates consistently satisfy the system constraints at every step.

The measurement function is designed with operational flexibility in mind. In single-cable mode, the filter ingests the three-dimensional position and velocity of a target point on one arresting cable. In four-cable mode, it measures the average of the target points on all four cables—effectively the position and velocity of the capture container center, the quantity that ultimately matters for steering the net onto the descending rocket. This dual-mode design allows the same filtering machinery to serve both local cable monitoring and global capture-point control.

Numerical simulations put the method through its paces, systematically examining the effects of integrator type, time step, augmented Lagrangian parameters, and the number of cable elements. The results are telling. In a baseline error-model simulation, the estimate deviates from the true value by 1.1 meters—beyond the allowable capture tolerance of 1 meter—whereas both error-state EKF implementations significantly outperform it. With a small time step of 1 millisecond, the displacement root-mean-square errors are 0.0066 meters for the generalized-α integrator and 0.026 meters for BDF, while the velocity RMSEs come in at 0.2951 and 0.2927 meters per second, respectively. The generalized-α method proves superior in both displacement accuracy and numerical stability: it remains stable as the time step increases, whereas BDF diverges beyond 1 millisecond. Most importantly for practical deployment, when the time step grows to 5 milliseconds, the ratio of CPU time to real time for the generalized-α method drops below 1, demonstrating genuine real-time operational capability.

Scaling up to the complete four-cable-net system, the team compared two estimation strategies. A four-channel decentralized estimation based on the single-cable model achieves accuracy comparable to centralized estimation built on the full-system model, while significantly improving computational efficiency. This finding has clear engineering significance: instead of paying the price of a monolithic high-dimensional model, operators can run parallel single-cable estimators and recover essentially the same capture-relevant accuracy at a fraction of the computational cost.

The implications extend beyond a single recovery architecture. By showing that a Kalman filter can operate directly on multibody DAEs—using stabilized integrators, index reduction, and constraint-violation correction to sidestep the classical DAE-to-ODE detour—the study offers an effective solution for state estimation in flexible cable-net recovery systems and technical support for the broader theoretical development of state estimation in complex flexible multibody systems. As reusable launch vehicles grow heavier and launch cadences accelerate, the ability to sense, in real time, the invisible dynamics of a flexible capture net may well determine whether cable-net recovery matures from an appealing concept into the workhorse of routine, low-cost spaceflight.

Subject of Research: Error-state Kalman filtering with multibody integrators for state estimation of flexible cable-net recovery systems for reusable launch vehicles

Article Title: An error-state kalman filter with multibody integrator for state estimation of flexible cable-net systems

Article References: An error-state kalman filter with multibody integrator for state estimation of flexible cable-net systems. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: Kalman filtering, error-state estimation, multibody dynamics, differential-algebraic equations, flexible cable-net, reusable launch vehicles, generalized-alpha method, backward differentiation formula, Absolute Nodal Coordinate Formulation, constraint violation correction, real-time estimation, lidar measurements

The research team, led by Hu Haiyan of the university’s College of Aerospace Engineering, focused on planar deployable modules with high packaging ratios, meaning they compress into a small volume for launch and expand dramatically once in space. That property is essential for efficient use of scarce fairing space, but it comes at a cost. After deployment, residual strain in the structural members, manufacturing and assembly tolerances, and dynamic coupling effects between moving parts all conspire to distort the module away from its nominal geometry. Existing research in this field has largely concentrated on optimizing the accuracy of the final deployed configuration, treating deformation as a problem to be solved before operations begin. What has been missing, the authors argue, is a cooperative strategy that controls deformation and motion simultaneously during the assembly process itself, together with a dynamic model rigorous enough to capture the interplay between rigid bodies and continuously lengthening flexible components.

The hardware at the center of the study is a deployable assembly module built from three hexagonal mechanisms connected by thin-walled carbon-fiber shells. Each hexagonal unit integrates motors, magnets, electromagnets, and thrusters, giving the module a full complement of functions: deployment actuation, in-flight configuration adjustment, docking guidance and latching, and planar motion control. Once deployed, the module takes on a triangular configuration whose three side lengths correspond to the deployed lengths of the three carbon-fiber shell segments. Because each segment is driven independently by its own motor, the module can deliberately change its baseline geometry, retracting one side while extending another to reshape its triangle on command. This variable-baseline capability is the key innovation that allows the module to correct its own shape errors rather than simply living with them.

Modeling such a system is far from straightforward. During deployment and shape correction, the carbon-fiber shells undergo large overall motion and large deformation while their effective length changes continuously over time. Conventional finite elements with fixed length simply cannot describe this behavior accurately. The team therefore turned to the Arbitrary Lagrangian–Eulerian Absolute Nodal Coordinate Formulation, abbreviated ALE-ANCF, to build a dynamic model of the variable-length thin-shell elements. This formulation captures the additional inertial effects induced by material flowing into and out of a shell segment as its deployed length changes, as well as the geometric nonlinearities that arise during large deformation. For the rigid hexagonal mechanisms, the researchers used the natural coordinate method, and the rigid–flexible coupling was enforced through constraint equations linking the two descriptions.

The dynamic model goes beyond structure alone. It incorporates magnetic force models, which govern the attractive guidance used during docking, and thruster thrust models, which describe how the module’s nozzles generate forces and torques. When the module finally docks with a target structure, the change in topology produces a velocity discontinuity, and the framework handles this event through impulse constraints. The result is a complete set of rigid–flexible coupled multibody dynamic equations that provides an accurate mechanical foundation for the control algorithms built on top of it. This level of modeling fidelity matters because control decisions made on an inaccurate model translate directly into positioning errors at the moment of docking, when tolerances are measured in millimeters or less.

For control, the team designed a strategy based on proportional–derivative, or PD, control, split into two cooperating layers. The deformation controller’s job is to drive the module’s current triangular configuration to a target configuration. The system computes the three current side lengths from the center positions of the hexagonal mechanisms and compares them with target side lengths measured from the assembly structure it intends to join. The PD controller then adjusts the rotational speeds of the three motors to retract or extend the corresponding carbon-fiber segments, with the geometric relationship between motor rotation angles and side-length changes providing the precise conversion between commanded actuation and physical shape change. In effect, the module measures its own distortion and actively machines it away.

Motion control, by contrast, handles the module’s journey through space. It uses a six-nozzle cooperative scheme in which the controller computes the desired control force and torque from the deviations between the module’s current planar position and attitude and their target values. That force and torque demand is then allocated among the six nozzles through pulse-width modulation, allowing the module to adjust its position and attitude in a synchronized fashion rather than correcting one at the expense of the other. The two control layers can operate independently or in an integrated temporal sequence: the module first performs deformation control to eliminate geometric errors, then initiates motion control to approach the assembly structure, and finally completes docking and latching under magnetic guidance from the electromagnets and magnets built into its hexagonal units.

The validation effort combined numerical simulation with ground-based air-floating experiments, and the numbers are striking. In simulation, across three sets of different target side-length conditions, the deformation controller reduced geometric errors to below 0.013 millimeters in every case, an accuracy improvement exceeding 99.98 percent, while out-of-plane displacement and attitude variation remained extremely small. Assembly simulations that integrated deformation and motion control showed the module achieving successful docking whether the target assembly structure was fixed in place or freely floating, demonstrating robustness to the realistic scenario of a compliant, drifting target. The air-floating platform experiments, which emulate the frictionless planar dynamics of a spacecraft on a cushion of air, told a slightly less perfect but still compelling story: geometric errors were reduced to below 1.5 millimeters, with accuracy improvements exceeding 98 percent across three experimental trials.

The gap between simulated and experimental performance is itself informative. The error curves from simulation and experiment were generally consistent, with the deviations attributed mainly to friction disturbances from the air-floating platform and fluctuations in the measurement system, rather than to any flaw in the control strategy itself. The decisive experiment came when the full integrated sequence was tested. When the module executed deformation control before initiating its approach, it docked successfully with the simulated assembly structure. When deformation control was deliberately omitted, the accumulated geometric deviation grew large enough that docking failed outright. That side-by-side comparison, captured in the experimental photographs, makes the central point of the paper impossible to miss: for deployable modules, shape correction is not an optional refinement but a prerequisite for autonomous assembly.

The broader significance of this work lies in its integration of two threads that have often been pursued separately. By pairing a rigid–flexible coupled dynamic model capable of handling variable-length flexible components with a cooperative deformation-and-motion control architecture, the researchers have provided a complete pipeline from mechanical understanding to executable control law. The demonstrated error tolerances, from sub-0.013-millimeter precision in simulation to sub-1.5-millimeter performance on a physical testbed, suggest that the approach can bridge the gap between theoretical modeling and engineering practice. As space agencies and commercial operators move toward on-orbit construction of structures too large for any single launch, techniques of this kind, in which modules measure, reshape, and reposition themselves before committing to contact, are likely to become foundational. The study offers both theoretical support and a practical engineering reference for the on-orbit construction campaigns that may define the next era of space infrastructure.

Subject of Research: Dynamic modeling and cooperative deformation and motion control of planar deployable modules for on-orbit autonomous assembly of large space structures

Article Title: Dynamic analysis and experiments on the assembly control of planar deployable modules

Article References: Dynamic analysis and experiments on the assembly control of planar deployable modules. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: on-orbit assembly, deployable structures, rigid-flexible coupling, ALE-ANCF, PD control, air-floating experiment, docking control, space infrastructure, variable-baseline configuration, multibody dynamics, carbon-fiber shells, autonomous docking