Udwadia–Kalaba equation
This article has multiple issues. Please help improve it or discuss these issues on the talk page.

Part of a series on 
Classical mechanics 

In theoretical physics, the Udwadia–Kalaba equation is a method for deriving the equations of motion of a constrained mechanical system.^{[1]} The equation was first described by Firdaus E. Udwadia and Robert E. Kalaba in 1992.^{[2]} The approach is based on Gauss's principle of least constraint. The Udwadia–Kalaba equation applies to both holonomic constraints and nonholonomic constraints, as long as they are linear with respect to the accelerations. The equation generalizes to constraint forces that do not obey D'Alembert's principle.^{[3]}^{[4]}^{[5]}
Background
The Udwadia–Kalaba equation was developed in 1992 and describes the motion of a constrained mechanical system that is subjected to equality constraints.^{[2]}
This differs from the Lagrangian formalism, which uses the Lagrange multipliers to describe the motion of constrained mechanical systems, and other similar approaches such as the GibbsAppell approach. The physical interpretation of the equation has applications in areas beyond theoretical physics, such as the control of highly nonlinear general dynamical systems.^{[6]}
The central problem of constrained motion
In the study of the dynamics of mechanical systems, the configuration of a given system S is, in general, completely described by n generalized coordinates so that its generalized coordinate nvector is given by
where T denotes matrix transpose. Using Newtonian or Lagrangian dynamics, the unconstrained equations of motion of the system S under study can be derived as a matrix equation (see matrix multiplication):
where the dots represent derivatives with respect to time:
It is assumed that the initial conditions q(0) and are known. We call the system S unconstrained because may be arbitrarily assigned.
The nvector Q denotes the total generalized force acted on the system by some external influence; it can be expressed as the sum of all the conservative forces as well as nonconservative forces.
The nbyn matrix M is symmetric, and it can be positive definite or semipositive definite . Typically, it is assumed that M is positive definite; however, it is not uncommon to derive the unconstrained equations of motion of the system S such that M is only semipositive definite; i.e., the mass matrix may be singular (it has no inverse matrix).^{[7]}^{[8]}
Constraints
We now assume that the unconstrained system S is subjected to a set of m consistent equality constraints given by
where A is a known mbyn matrix of rank r and b is a known mvector. We note that this set of constraint equations encompass a very general variety of holonomic and nonholonomic equality constraints. For example, holonomic constraints of the form
can be differentiated twice with respect to time while nonholonomic constraints of the form
can be differentiated once with respect to time to obtain the mbyn matrix A and the mvector b. In short, constraints may be specified that are
 nonlinear functions of displacement and velocity,
 explicitly dependent on time, and
 functionally dependent.
As a consequence of subjecting these constraints to the unconstrained system S, an additional force is conceptualized to arise, namely, the force of constraint. Therefore, the constrained system S_{c} becomes
where Q_{c}—the constraint force—is the additional force needed to satisfy the imposed constraints. The central problem of constrained motion is now stated as follows:
 given the unconstrained equations of motion of the system S,
 given the generalized displacement q(t) and the generalized velocity of the constrained system S_{c} at time t, and
 given the constraints in the form as stated above,
find the equations of motion for the constrained system—the acceleration—at time t, which is in accordance with the agreed upon principles of analytical dynamics.
Equation of motion
The solution to this central problem is given by the Udwadia–Kalaba equation. When the matrix M is positive definite, the equation of motion of the constrained system S_{c}, at each instant of time, is^{[2]}^{[9]}
where the '+' symbol denotes the pseudoinverse of the matrix . The force of constraint is thus given explicitly as
and since the matrix M is positive definite the generalized acceleration of the constrained system S_{c} is determined explicitly by
In the case that the matrix M is semipositive definite , the above equation cannot be used directly because M may be singular. Furthermore, the generalized accelerations may not be unique unless the (n + m)byn matrix
has full rank (rank = n).^{[7]}^{[8]} But since the observed accelerations of mechanical systems in nature are always unique, this rank condition is a necessary and sufficient condition for obtaining the uniquely defined generalized accelerations of the constrained system S_{c} at each instant of time. Thus, when has full rank, the equations of motion of the constrained system S_{c} at each instant of time are uniquely determined by (1) creating the auxiliary unconstrained system^{[8]}
and by (2) applying the fundamental equation of constrained motion to this auxiliary unconstrained system so that the auxiliary constrained equations of motion are explicitly given by^{[8]}
Moreover, when the matrix has full rank, the matrix is always positive definite. This yields, explicitly, the generalized accelerations of the constrained system S_{c} as
This equation is valid when the matrix M is either positive definite or positive semidefinite. Additionally, the force of constraint that causes the constrained system S_{c}—a system that may have a singular mass matrix M—to satisfy the imposed constraints is explicitly given by
Nonideal constraints
At any time during the motion we may consider perturbing the system by a virtual displacement δr consistent with the constraints of the system. The displacement is allowed to be either reversible or irreversible. If the displacement is irreversible, then it performs virtual work. We may write the virtual work of the displacement as
The vector describes the nonideality of the virtual work and may be related, for example, to friction or drag forces (such forces have velocity dependence). This is a generalized D'Alembert's principle, where the usual form of the principle has vanishing virtual work with .
The Udwadia–Kalaba equation is modified by an additional nonideal constraint term to
Examples
Inverse Kepler problem
The method can solve the inverse Kepler problem of determining the force law that corresponds to the orbits that are conic sections.^{[10]} We take there to be no external forces (not even gravity) and instead constrain the particle motion to follow orbits of the form
where , is the eccentricity, and is the semilatus rectum. Differentiating twice with respect to time and rearranging slightly gives a constraint
We assume the body has a simple, constant mass. We also assume that angular momentum about the focus is conserved as
with time derivative
We can combine these two constraints into the matrix equation
The constraint matrix has inverse
The force of constraint is therefore the expected, central inverse square law
Inclined plane with friction
Consider a small block of constant mass on an inclined plane at an angle above horizontal. The constraint that the block lie on the plane can be written as
After taking two time derivatives, we can put this into a standard constraint matrix equation form
The constraint matrix has pseudoinverse
We allow there to be sliding friction between the block and the inclined plane. We parameterize this force by a standard coefficient of friction multiplied by the normal force
Whereas the force of gravity is reversible, the force of friction is not. Therefore, the virtual work associated with a virtual displacement will depend on C. We may summarize the three forces (external, ideal constraint, and nonideal constraint) as follows:
Combining the above, we find that the equations of motion are
This is like a constant downward acceleration due to gravity with a slight modification. If the block is moving up the inclined plane, then the friction increases the downward acceleration. If the block is moving down the inclined plane, then the friction reduces the downward acceleration.
References
 ^ Udwadia, F. E.; Kalaba, R. E. (1996). Analytical Dynamics: A New Approach. Cambridge: Cambridge University Press. ISBN 0521048338.
 ^ ^{a} ^{b} ^{c} Udwadia, F. E.; Kalaba, R. E. (1992). "A new perspective on constrained motion" (PDF). Proceedings of the Royal Society of London, Series A. 439 (1906): 407–410. Bibcode:1992RSPSA.439..407U. doi:10.1098/rspa.1992.0158.
 ^ Udwadia, F. E.; Kalaba, R. E. (2002). "On the Foundations of Analytical Dynamics" (PDF). International Journal of Nonlinear Mechanics. 37 (6): 1079–1090. Bibcode:2002IJNLM..37.1079U. CiteSeerX 10.1.1.174.5726. doi:10.1016/S00207462(01)000336.
 ^ Calverley, B. (2001). "Constrained or Unconstrained, That Is the Equation". USC News.
 ^ Udwadia, F.; Kalaba, R. (2002). "What is the General Form of the Explicit Equations of Motion for Constrained Mechanical Systems?" (PDF). Journal of Applied Mechanics. 69 (3): 335–339. Bibcode:2002JAM....69..335U. CiteSeerX 10.1.1.174.6353. doi:10.1115/1.1459071.
 ^ Zhao, Xiao; Chen, YeHwa; Zhao, Han; Dong, FangFang (2018). "Udwadia–Kalaba equation for constrained mechanical systems:Formulation and applications" (pdf). Chinese Journal of Mechanical Engineering. 31 (1): 106–120. doi:10.1186/s100330180310x.
 ^ ^{a} ^{b} Udwadia, F.E.; Phohomsiri, P. (2006). "Explicit equations of motion for constrained mechanical systems with singular mass matrices and applications to multibody dynamics" (PDF). Proceedings of the Royal Society of London, Series A. 462 (2071): 2097–2117. Bibcode:2006RSPSA.462.2097U. doi:10.1098/rspa.2006.1662.
 ^ ^{a} ^{b} ^{c} ^{d} Udwadia, F.E.; Schutte, A.D. (2010). "Equations of motion for general constrained systems in Lagrangian mechanics" (PDF). Acta Mechanica. 213 (1): 111–129. doi:10.1007/s0070700902722.
 ^ Udwadia, F.E.; Kalaba, R.E. (1993). "On motion" (PDF). Journal of the Franklin Institute. 330 (3): 571–577. doi:10.1016/00160032(93)90099G.
 ^ Zhang, Bingzhan; Zhen, Shengchao; Zhao, Han; Huang, Kang; Deng, Bin; Chen, YeHwa (2015). "A novel study on Kepler's law and inverse square law of gravitation". Eur. J. Phys. 36 (3): 035018. Bibcode:2015EJPh...36c5018Z. doi:10.1088/01430807/36/3/035018.