---
title: Rigid-Link Gaussian Kinematics
url: https://www.emergentmind.com/topics/rigid-link-gaussian-kinematics-rlgk
type: topic
---

# Rigid-Link Gaussian Kinematics

Rigid-Link Gaussian Kinematics (RLGK) provides a unified variational formulation for the dynamics and kinematics of systems comprising interconnected rigid bodies. Instead of isolating Lagrange, Newton–Euler, or Kirchhoff frameworks, RLGK bases all on Gauss’s principle of least constraint. The resulting equations of motion naturally encompass all classical rigid-body dynamics formulations through judicious selection of velocity coordinates or “quasicoordinates.” RLGK also seamlessly connects with modern representation theory on $\mathrm{SE}(3)$ and projectivized dual quaternion models for mechanism constraint varieties, linking algebraic and geometric approaches to robot mechanism synthesis and analysis [1607.00156, 1607.08119, 1503.04566].

## 1. Gauss’s Principle and the RLGK Variational Functional

Gauss’s principle of least constraint identifies the physically realized accelerations among all kinematically admissible candidates as those minimizing the quadratic deviation from the free-body accelerations. For a chain of $n$ rigid bodies with unconstrained accelerations $\{a_i^*\}$ (subject only to external forces) and constraint-compatible accelerations $\{a_i\}$, the Gauss functional is:
$$
Z = \tfrac{1}{2} \sum_{i=1}^n m_i \left\| a_i - a_i^* \right\|^2
$$
The requirement $\delta Z = 0$ under all virtual variations respecting constraints yields the actual system acceleration [1607.00156].

When applied to an $n$-link rigid chain, each $a_i$ is determined from generalized coordinates $q$, their velocities $\dot q$, and accelerations $\ddot q$ via:
\[
a_i = A_i(q, \dot q) \ddot q + b_i(q, \dot q)
\]
where $A_i$ is the twist-to-Cartesian Jacobian for body $i$ and $b_i$ groups the Christoffel-type velocity-dependent (Coriolis and centripetal) terms. The unconstrained acceleration is
\[
a_i^* = M_i^{-1} f_i^{\rm ext}
\]
with $M_i$ the spatial inertia. Substitution into $Z$ yields
\[
Z(q, \dot q, \ddot q) = \frac{1}{2} \sum_{i=1}^n \left[ A_i \ddot q + b_i - a_i^* \right]^\top M_i \left[ A_i \ddot q + b_i - a_i^* \right]
\]
The stationarity condition $\delta Z = 0$ imposes a well-defined variational principle for the governed motion [1607.00156].

## 2. Unified Equations of Motion and Quasicoordinate Choice

The RLGK principle leads to the common form of rigid-body dynamics:
\[
M(q) \ddot q + C(q, \dot q) \dot q = \tau(q, \dot q)
\]
where the terms are constructed as:
- Generalized mass matrix:
  \[
  M(q) = \sum_{i=1}^n A_i(q)^\top M_i A_i(q)
  \]
- Velocity-dependent (Coriolis/centrifugal) terms:
  \[
  C(q, \dot q) \dot q = \sum_{i=1}^n A_i(q)^\top M_i b_i(q, \dot q)
  \]
- Generalized force/torque vector:
  \[
  \tau(q, \dot q) = \sum_{i=1}^n A_i(q)^\top f_i^{\rm ext}
  \]
The essence is that $\{\dot q\}$ can be chosen as arbitrary linear velocity coordinates—joint rates (Lagrange form), twist vectors (Newton–Euler form), or screw-velocities (Kirchhoff form), among others. Each classical rigid-body formalism is recovered by the appropriate selection of quasicoordinates and the corresponding Jacobians $A_i$, maintaining mathematical equivalence at the level of $\delta Z = 0$ [1607.00156].

## 3. Group-Theoretical Structure and Trivializations on $\mathrm{SE}(3)$

The configuration space for a single rigid body is the Lie group $\mathrm{SE}(3)$, coordinatized as
\[
g = \begin{pmatrix} R & p \\ 0 & 1 \end{pmatrix}
\]
with $R \in \mathrm{SO}(3)$ and $p \in \mathbb{R}^3$. Tangent velocities reside in $\mathfrak{se}(3)$ as twists. There are two natural trivializations:
- **Left-trivialization (body frame):** Spatial inertia $M_i$ remains constant for each link, and the Coriolis term $\mathrm{ad}^\ast_V(M_iV)$ becomes sparse.
- **Right-trivialization (world frame):** Yields simplified articulated Jacobians for serial chains.

By choosing trivialization, one can often achieve desirable properties: either a constant mass matrix (useful for symbolic computation) or a skew-symmetric Coriolis matrix (beneficial for energy analysis). The compatibility of RLGK with group-theoretic approaches is central for efficient symbolic and numeric implementation in mechanism dynamics [1607.00156].

## 4. Constraint Varieties in the Projectivized Dual Quaternion Model

RLGK bridges with the dual quaternion formalism, particularly for kinematic chains (dyads, linkages). In this setting, rigid displacements correspond to points in $\mathbb{P}^7$ (projective dual quaternion space), constrained to the Study quadric $S:N(q)\in\mathbb{R}$. Dyads are represented as the intersection of 3-spaces $U$ with $S$, yielding regular ruled quadrics $Q=U\cap S$. The type of dyad (RR, PR, RP, or cylindrical) is distinguished by the incidence structure of “null” lines (generators) in $Q\cap N$ (with $N$ the null cone) and the associated fiber-projectivity invariant $\varphi$ [1607.08119].

A summary of constraint incidence types is:

| Dyad Type | Ruled Quadric $Q$ | Null Incidence                      |
|-----------|-------------------|--------------------------------------|
| RR        | $Q$ contains      | Null quadrilateral—4 lines in $Q\cap N$ |
| PR/RP     | $Q$ contains      | 2 complex conjugate null lines $+\;$ 1 real null line $e_1$ |
| Cylinder  | $Q$ with fibers   | All null lines collapse: $\varphi(Q) = e_1$ |

Special coordinate choices and projective invariants discriminate among dyad classes and encode commutation and geometric singularities [1607.08119].

## 5. RLGK in Two-Link and Dyadic Mechanism Synthesis

For concrete illustration, consider a planar 2-link arm (two revolute joints), with joint angles $(q^1,q^2)$, link parameters $(l_1,l_2,m_1,m_2,I_1,I_2)$:
- Link center positions:
  \[
  \begin{aligned}
  p_1(q) &= \tfrac12 l_1 (\cos q^1, \sin q^1)^\top \\
  p_2(q) &= l_1 (\cos q^1, \sin q^1)^\top + \tfrac12 l_2 (\cos(q^1+q^2), \sin(q^1+q^2))^\top
  \end{aligned}
  \]
- Mass matrix (block form):
  \[
  M(q) = \begin{pmatrix} a(q) & b(q) \\ b(q) & c(q) \end{pmatrix}
  \]
  with explicit expressions for $a(q)$, $b(q)$, $c(q)$ (see source for details).
- Coriolis terms:
  \[
  C(q,\dot q)\dot q = \begin{pmatrix} -m_2 l_1 l_2 \sin q^2\,\dot q^2(2\dot q^1+\dot q^2) \\ m_2 l_1 l_2 \sin q^2 (\dot q^1)^2 \end{pmatrix}
  \]
- RLGK equation instantiated:
  \[
  \begin{aligned}
  a(q)\,\ddot q^1 + b(q)\,\ddot q^2 - m_2 l_1 l_2\sin q^2\,(2\dot q^1+\dot q^2)\dot q^2 &= \tau_1 \\
  b(q)\,\ddot q^1 + c(q)\,\ddot q^2 + m_2 l_1 l_2\sin q^2\,(\dot q^1)^2 &= \tau_2
  \end{aligned}
  \]

Mechanism constraint varieties for 2R dyads (two revolute axes) comprise regular ruled quadrics in $\mathbb{P}^3$ with null quadrilaterals (four null lines at $t_1=\pm i$, $t_2=\pm i$) [1503.04566]. Parameterizations and ruling-family structure correlate with feasible joint motions. For exact linkage synthesis (e.g., interpolating four poses), the twisted-cubic theory on the quadric gives two admissible classes of 5R linkages, each characterized by a commutation relation among factors in a cubic motion polynomial [1503.04566].

## 6. Synthesis, Invariants, and Gaussian Complexes

RLGK’s correspondence with Gaussian complexes in projective kinematics associates 1- and 2-joint dyads with quadrics stratified by null incidence. The fiber-projectivity invariant $\varphi$ is crucial: for PR/RP dyads, $\varphi(Q)=e_1$ identifies the cylinder case (commuting joints), while for generic PR/RP, $\varphi(Q) \neq e_1$, distinguishing these from the RR case (null quadrilateral) [1607.08119]. This stratification refines the traditional Gaussian classification, with each quadrics’ null content encoding the dyad’s kinematic mobility and singularity structure.

Admissible transformations (frame changes) correspond to projectivities fixing the Study quadric $S$ and null cone $N$. Such transformations carry RR, PR, RP, and cylinder dyad spaces into one another, preserving the nuanced distinction of ruling families and null incidence, which is essential for mechanism classification and synthesis [1607.08119].

## 7. Connections and Implications for Mechanism Science

The RLGK formalism provides a synthesis of geometric, analytic, and algebraic perspectives on rigid-body dynamics and kinematics. Its compatibility with coordinate-invariant, group-theoretic, and projective methods enables:
- Unified derivation of classical and modern equations of motion.
- Seamless translation between Newton–Euler, Lagrange, and Kirchhoff kinetic representations.
- Direct integration with dual quaternion and Plücker coordinate representations for linkage synthesis, classification, and pose interpolation tasks.
- Rigorous algebraic classification of mechanism constraint varieties via Gaussian complexes and null incidence.
- Provision of explicit geometric invariants for mechanism stratification (e.g., fiber-projectivity, null quadrilaterals).

These features make RLGK a central framework for both theoretical understanding and computational treatment of spatial mechanisms and their synthesis, as evidenced in the cited literature [1607.00156, 1607.08119, 1503.04566].

Source: https://www.emergentmind.com/topics/rigid-link-gaussian-kinematics-rlgk