Bayu Jayawardhana

Arjan van der Schaft, Shodhan Rao, Bayu Jayawardhana

Motivated by recent progress on the interplay between graph theory, dynamics, and systems theory, we revisit the analysis of chemical reaction networks described by mass action kinetics. For reaction networks possessing a thermodynamic equilibrium we derive a compact formulation exhibiting at the same time the structure of the complex graph and the stoichiometry of the network, and which admits a direct thermodynamical interpretation. This formulation allows us to easily characterize the set of equilibria and their stability properties. Furthermore, we develop a framework for interconnection of chemical reaction networks. Finally we discuss how the established framework leads to a new approach for model reduction.

Shodhan Rao, Arjan van der Schaft, Bayu Jayawardhana

In this paper we derive a compact mathematical formulation describing the dynamics of chemical reaction networks that are complex-balanced and are governed by mass action kinetics. The formulation is based on the graph of (substrate and product) complexes and the stoichiometric information of these complexes, and crucially uses a balanced weighted Laplacian matrix. It is shown that this formulation leads to elegant methods for characterizing the space of all equilibria for complex-balanced networks and for deriving stability properties of such networks. We propose a method for model reduction of complex-balanced networks, which is similar to the Kron reduction method for electrical networks and involves the computation of Schur complements of the balanced weighted Laplacian matrix.

Claudio De Persis, Bayu Jayawardhana

In this paper we investigate a passivity approach to collective coordination and synchronization problems in the presence of quantized measurements and show that coordination tasks can be achieved in a practical sense for a large class of passive systems.

Rodlfo Reyes-Báez, Arjan van der Schaft, Bayu Jayawardhana

Based on recent advances in contraction methods in systems and control, in this paper we present the virtual differential passivity based control (v-dPBC) technique. This is a constructive design method that combines the concept of virtual systems and of differential passivity. We apply the method to the tracking control problem of flexible joints robots (FJRs) which are formulated in the port-Hamiltonian (pH) framework. Simulations on a single flexible joint link are presented for showing the performance of a controller obtained with this approach.

Rodolfo Reyes-Báez, Alejandro Donaire, Arjan van der Schaft et al.

In this work we propose a family of trajectory tracking controllers for marine craft in the port-Hamiltonian (pH) framework using virtual differential passivity based control (v-dPBC). Two pH models of marine craft are considered, one in a body frame and another in an inertial frame. The structure and workless forces of pH models are exploited to design two virtual control systems which are related to the original marine craft's pH models. These virtual systems are rendered differentially passive with an imposed steady-state trajectory, both by means of a control scheme. Finally, the original marine craft pH models in closed-loop with above controllers solve the trajectory tracking problem. The performance of the closedloop system is evaluated on numerical simulations.

Shodhan Rao, Arjan van der Schaft, Karen van Eunen et al.

In this paper we propose a model-order reduction method for chemical reaction networks governed by general enzyme kinetics, including the mass-action and Michaelis-Menten kinetics. The model-order reduction method is based on the Kron reduction of the weighted Laplacian matrix which describes the graph structure of complexes in the chemical reaction network. We apply our method to a yeast glycolysis model, where the simulation result shows that the transient behaviour of a number of key metabolites of the reduced-order model is in good agreement with those of the full-order model.

Hector Garcia de Marina, Bayu Jayawardhana, Ming Cao

This paper presents a distributed algorithm for controlling the deployment of a team of mobile agents in formations whose shapes can be characterized by a broad class of polygons, including regular ones, where each agent occupies a corner of the polygon. The algorithm shares the appealing properties of the popular distance-based rigid formation control, but with the additional advantage of requiring the control of fewer pairs of neighboring agents. Furthermore, the scale of the polygon can be controlled by only one pair of neighboring agents. We also exploit the exponential stability of the controlled formation in order to steer the formation as a whole with translations and rotations in a prescribed way. We provide both theoretical analysis and illustrative simulations.

Rodolfo Reyes-Báez, Arjan van der Schaft, Bayu Jayawardhana

In this paper, we propose a novel trajectory tracking controller for fully-actuated mechanical port-Hamiltonian (pH) systems, which is based on recent advances in contraction-based control theory. Our proposed controller renders a desired sliding manifold (where the reference trajectory lies) attractive by making the corresponding error system partially contracting. Finally, we present numerical simulation results where a SCARA robot is commanded by our proposed tracking control law.

Ruiyue Ouyang, Vincent Andrieu, Bayu Jayawardhana

In this paper we investigate the dissipativity property of a certain class of Duhem hysteresis operator, which has clockwise (CW) input-output (I/O) behavior. In particular, we provide sufficient conditions on the Duhem operator such that it is CW and propose an explicit construction of the corresponding storage function satisfying dissipation inequality of CW systems. The result is used to analyze the stability of a second order system with hysteretic friction which is described by a Dahl model.

Muhammad Zakiyullah Romdlony, Bayu Jayawardhana

In this paper, we propose a new robustness notion that is applicable for certifying systems' safety with respect to external disturbance signals. The proposed input-to-state safety (ISSf) notion allows us to certify systems' safety in the presence of the disturbances which is analogous to the notion of input-to-state stability (ISS) for analyzing systems' stability.

Rodolfo Reyes-Baez, Arjan van der Schaft, Bayu Jayawardhana

In this paper we present three distributed control laws for the coordination of networked Euler-Lagrange (EL) systems. We first reformulate the passivity-based control design method in \cite{Arcak} by considering that each edge is associated with an \emph{artificial spring system} instead of the usual diffusive coupling among the communicating agents. With this configuration, the networked EL system possesses a "symmetric" feedback structure which together with the strict passivity of both agents' and edges' dynamics lead to a strictly passive network dynamics. Subsequently we present the networked version of two different passivity-based tracking controllers %local controllers that are particular cases of our method and the one in \cite{Arcak}. Numerical simulation is presented to show the performance of the proposed methods.

Jurrien Keulen, Bayu Jayawardhana, Arjan van der Schaft

This paper presents a port-Hamiltonian formulation of hysteretic energy storage elements. First, we revisit the passivity property of backlash-driven storage elements by presenting a family of storage functions associated to the dissipativity property of such elements. We explicitly derive the corresponding available storage and required supply functions `a la Willems [1], and show the interlacing property of the aforementioned family of storage functions sandwiched between the available storage and required supply functions. Second, using the proposed family of storage functions, we present a port-Hamiltonian formulation of hysteretic inductors as prototypical storage elements in port-Hamiltonian systems. In particular, we show how a Hamiltonian function can be chosen from the family of storage functions and how the hysteretic elements can be expressed as port-Hamiltonian system with feedthrough term, where the feedthrough term represents energy dissipation. Correspondingly, we illustrate its applicability in describing an RLC circuit (in parallel and in series) containing a hysteretic inductor element.

Ruiyue Ouyang, Bayu Jayawardhana

In this paper, we investigate the stability of a feedback interconnection between a linear system and a Duhem hysteresis operator, where the linear system and the Duhem hysteresis operator satisfy either the counter-clockwise (CCW) or clockwise (CW) input-output dynamics. More precisely, we present sufficient conditions for the stability of the interconnected system that depend on the CW or CCW properties of the linear system and the Duhem operator. Based on these results we introduce a control design methodology for stabilizing a linear plant with a hysteretic actuator or sensor without requiring precise information on the hysteresis operator.

Yanxin Si, Bayu Jayawardhana, J. Nathan Kutz et al.

Experimental Modal Analysis (EMA) has been widely used to identify structural dynamic properties, including natural frequencies, damping ratios, and mode shapes, for structural integrity assessment. The Poly-reference Least Squares Complex Frequency (pLSCF) method is one of the most widely adopted approaches for EMA because of its strong ability to separate closely spaced modes and its robustness to measurement noise. However, pLSCF-based EMA is generally limited to low-dimensional cases with a small number of measurement points, as its computational cost increases rapidly for high-dimensional or continuous structural measurements, particularly with increasing model order. To overcome this limitation, this paper develops a high-dimensional EMA framework based on Dynamic Mode Decomposition with control (DMDc), a powerful data-driven technique originally developed in fluid dynamics, for modal identification under high-dimensional measurement scenarios. Specifically, the relationship between pLSCF and time-delay DMDc is clarified through the discrete state-space representation of the auto-regressive with exogenous inputs (ARX) model for linear systems. By showing that both methods describe the same physical dynamics of the structure, this study provides a physics-based rationale for applying time-delay DMDc to EMA. The capability and advantages of time-delay DMDc for modal parameter identification in both low- and high-dimensional measurements are validated through numerical simulations of a 6-DOF system and experiments on a cantilever beam using a digital camera. The results demonstrate that time-delay DMDc enables robust and reliable modal parameter identification, effectively addressing high-dimensional EMA problems that are difficult for conventional pLSCF and highlighting its potential for real-world structural dynamics applications.

Moh Kamalul Wafi, Ahmad Ataka, Yul Y. Nazaruddin et al.

This paper studies distributed trajectory tracking for networks of nonholonomic mobile robots under adversarial information exchange. An exact global input--output feedback linearization scheme is developed to regulate planar position outputs, yielding linear error dynamics without prescribing internal state trajectories. To mitigate corrupted neighbor information, a resilient desired-signal construction is proposed that combines local redundancy with trusted in-neighbor signals, without requiring adversary detection or isolation. When sufficient redundancy is available, the method suppresses adversarial influence and recovers nominal tracking performance. If redundancy conditions are violated, adversarial effects enter as bounded disturbances and the tracking error remains ultimately bounded. Simulation results on star, cyclic, and path topologies validate the analysis and demonstrate the superior resilience of cyclic networks due to distributed information propagation.

Santiago Sanchez-Escalonilla, Samuele Zoboli, Bayu Jayawardhana

In this paper, we restructure the Neural Interconnection and Damping Assignment - Passivity Based Control (Neural IDA-PBC) design methodology, and we formally analyze its closed-loop properties. Neural IDA-PBC redefines the IDA-PBC design approach as an optimization problem by building on the framework of Physics Informed Neural Networks (PINNs). However, the closed-loop stability and robustness properties under Neural IDA-PBC remain unexplored. To address the issue, we study the behavior of classical IDA-PBC under approximations. Our theoretical analysis allows deriving conditions for practical and asymptotic stability of the desired equilibrium point. Moreover, it extends the Neural IDA-PBC applicability to port-Hamiltonian systems where the matching conditions cannot be solved exactly. Our renewed optimization-based design introduces three significant aspects: i) it involves a novel optimization objective including stability and robustness constraints issued from our theoretical analysis; ii) it employs separate Neural Networks (NNs), which can be structured to reduce the search space to relevant functions; iii) it does not require knowledge about the port-Hamiltonian formulation of the system's model. Our methodology is validated with simulations on three standard benchmarks: a double pendulum, a nonlinear mass-spring-damper and a cartpole. Notably, classical IDA-PBC designs cannot be analytically derived for the latter.

Santiago Sanchez-Escalonilla, Rodolfo Reyes-Baez, Bayu Jayawardhana

In this work we exploit the universal approximation property of Neural Networks (NNs) to design interconnection and damping assignment (IDA) passivity-based control (PBC) schemes for fully-actuated mechanical systems in the port-Hamiltonian (pH) framework. To that end, we transform the IDA-PBC method into a supervised learning problem that solves the partial differential matching equations, and fulfills equilibrium assignment and Lyapunov stability conditions. A main consequence of this, is that the output of the learning algorithm has a clear control-theoretic interpretation in terms of passivity and Lyapunov stability. The proposed control design methodology is validated for mechanical systems of one and two degrees-of-freedom via numerical simulations.

Haiwen Wu, Bayu Jayawardhana, Hector Garcia de Marina et al.

We present three classes of distributed formation controllers for achieving and maintaining the 2D/3D formation shape of manipulator end-effectors to cope with different scenarios due to availability of modeling parameters. We firstly present a distributed formation controller for manipulators whose system parameters are perfectly known. The formation control objective is achieved by assigning virtual springs between end-effectors and by adding damping terms at joints, which provides a clear physical interpretation of the proposed solution. Subsequently, we extend it to the case where manipulator kinematic and system parameters are not exactly known. An extra integrator and an adaptive estimator are introduced for gravitational compensation and stabilization, respectively. Simulation results with planar manipulators and with seven degree-of-freedom humanoid manipulator arms are presented to illustrate the effectiveness of the proposed approach.

Tinghua Li, Bayu Jayawardhana, Amar Kamat et al.

We present the design and experimental validation of source seeking control algorithms for a unicycle mobile robot that is equipped with novel 3D-printed flexible graphene-based piezoresistive airflow sensors. Based solely on a local gradient measurement from the airflow sensors, we propose and analyze a projected gradient ascent algorithm to solve the source seeking problem. In the case of partial sensor failure, we propose a combination of Extremum-Seeking Control with our projected gradient ascent algorithm. For both control laws, we prove the asymptotic convergence of the robot to the source. Numerical simulations were performed to validate the algorithms and experimental validations are presented to demonstrate the efficacy of the proposed methods.

Haiwen Wu, Bayu Jayawardhana, Hector Garcia de Marina et al.

This paper addresses the problem of end-effector formation control for manipulators that are subjected to external disturbances: input disturbance torques and disturbance forces at each end-effector. The disturbances are assumed to be non-vanishing and are superposition of finite number of sinusoidal and step signals. The formation control objective is achieved by assigning virtual springs between end-effectors, by adding damping terms at joints, and by incorporating internal model-based dynamic compensators to counteract the effect of the disturbances; all of which presents a clear physical interpretation of the proposed approach. Simulation results are presented to illustrate the effectiveness of the proposed approach.

Nelson P. K. Chan, Bayu Jayawardhana, Hector Garcia de Marina

This paper focuses on securing a triangular shape (up to translation) for a team of three mobile robots that uses heterogeneous sensing mechanism. Based on the available local information, each robot employs the popular gradient-based control law to attain the assigned individual task(s). In the current work, robots are assigned either distance and signed area task(s) or bearing task(s). We provide a sufficient condition on the gain ratio $R_{\text{Ad}}$ between the signed area and the distance control term such that the desired formation shape, an isosceles triangle, is reached from all feasible starting positions. Numerical simulations are provided to support the theoretical analyses.

Nelson P. K. Chan, Bayu Jayawardhana, Hector Garcia de Marina

In this letter, we investigate the formation control problem of mobile robots moving in the plane where, instead of assuming robots to be simple points, each robot is assumed to have the form of a disk with equal radius. Based on interior angle measurements of the neighboring robots' disk, which can be obtained from low-cost vision sensors, we propose a gradient-based distributed control law and show the exponential convergence property of the associated error system. By construction, the proposed control law has the appealing property of ensuring collision avoidance between neighboring robots. We also present simulation results for {a team} of four circular mobile robots forming a rectangular shape.

Rodolfo Reyes-Báez, Arjan van der Schaft, Bayu Jayawardhana et al.

In this work we present a constructive method to design a family of virtual contraction based controllers that solve the standard trajectory tracking problem of flexible-joint robots (FJRs) in the port-Hamiltonian (pH) framework. The proposed design method, called virtual contraction based control (v-CBC), combines the concepts of virtual control systems and contraction analysis. It is shown that under potential energy matching conditions, the closed-loop virtual system is contractive and exponential convergence to a predefined trajectory is guaranteed. Moreover, the closed-loop virtual system exhibits properties such as structure preservation, differential passivity and the existence of (incrementally) passive maps.

Hector Garcia de Marina, Johan Siemonsma, Bayu Jayawardhana et al.

In this paper, we present the design and implementation of a robust motion formation distributed control algorithm for a team of mobile robots. The primary task for the team is to form a geometric shape, which can be freely translated and rotated at the same time. This approach makes the robots to behave as a cohesive whole, which can be useful in tasks such as collaborative transportation. The robustness of the algorithm relies on the fact that each robot employs only local measurements from a laser sensor which does not need to be off-line calibrated. Furthermore, robots do not need to exchange any information with each other. Being free of sensor calibration and not requiring a communication channel helps the scaling of the overall system to a large number of robots. In addition, since the robots do not need any off-board localization system, but require only relative positions with respect to their neighbors, it can be aimed to have a full autonomous team that operates in environments where such localization systems are not available. The computational cost of the algorithm is inexpensive and the resources from a standard microcontroller will suffice. This fact makes the usage of our approach appealing as a support for other more demanding algorithms, e.g., processing images from onboard cameras. We validate the performance of the algorithm with a team of four mobile robots equipped with low-cost commercially available laser scanners.

Hector Garcia de Marina, Bayu Jayawardhana, Ming Cao

Recently it has been reported that biased range-measurements among neighboring agents in the gradient distance-based formation control can lead to predictable collective motion. In this paper we take advantage of this effect and by introducing distributed parameters to the prescribed inter-distances we are able to manipulate the steady-state motion of the formation. This manipulation is in the form of inducing simultaneously the combination of constant translational and angular velocities and a controlled scaling of the rigid formation. While the computation of the distributed parameters for the translational and angular velocities is based on the well-known graph rigidity theory, the parameters responsible for the scaling are based on some recent findings in bearing rigidity theory. We carry out the stability analysis of the modified gradient system and simulations in order to validate the main result.

Hector Garcia de Marina, Ming Cao, Bayu Jayawardhana

Despite the great success of using gradient-based controllers to stabilize rigid formations of autonomous agents in the past years, surprising yet intriguing undesirable collective motions have been reported recently when inconsistent measurements are used in the agents' local controllers. To make the existing gradient control robust against such measurement inconsistency, we exploit local estimators following the well known internal model principle for robust output regulation control. The new estimator-based gradient control is still distributed in nature and can be constructed systematically even when the number of agents in a rigid formation grows. We prove rigorously that the proposed control is able to guarantee exponential convergence and then demonstrate through robotic experiments and computer simulations that the reported inconsistency-induced orbits of collective movements are effectively eliminated.

Hector Garcia de Marina, Bayu Jayawardhana, Ming Cao

Recently it has been reported that range-measurement inconsistency, or equivalently mismatches in prescribed inter-agent distances, may prevent the popular gradient controllers from guiding rigid formations of mobile agents to converge to their desired shape, and even worse from standing still at any location. In this paper, instead of treating mismatches as the source of ill performance, we take them as design parameters and show that by introducing such a pair of parameters per distance constraint, distributed controller achieving simultaneously both formation and motion control can be designed that not only encompasses the popular gradient control, but more importantly allows us to achieve constant collective translation, rotation or their combination while guaranteeing asymptotically no distortion in the formation shape occurs. Such motion control results are then applied to (a) the alignment of formations orientations and (b) enclosing and tracking a moving target. Besides rigorous mathematical proof, experiments using mobile robots are demonstrated to show the satisfying performances of the proposed formation-motion distributed controller.

Hector Garcia de Marina, Bayu Jayawardhana, Ming Cao

This paper presents the analysis on the influence of distance mismatches on the standard gradient-based rigid formation control for second-order agents. It is shown that, similar to the first-order case as recently discussed in the literature, these mismatches introduce two undesired group behaviors: a distorted final shape and a steady-state motion of the group formation. We show that such undesired behaviors can be eliminated by combining the standard formation control law with distributed estimators. Finally, we show how the mismatches can be effectively employed as design parameters in order to control a combined translational and rotational motion of the formation.