Affine Control Interpretation of the Golf Swing — Part 4
Modal Approximation, Wrenches, and Pendulum Dynamics
These entries remain visible until evidence-backed adjudication changes their governed status.
- Open / High: Critique: Input-Dependent Boundary Conditions (The Grip Paradox) (
crit-input-dependent-boundary-conditions)
Part IV is the "engine room" of the theory. It contains the raw mathematics that prove the ideas in Parts I-III are solid.
Why This Exists
In science, you can't just claim a theory works—you have to prove it mathematically. This section derives the equations of motion from scratch to show that our "Drift vs. Input" separation isn't just a nice idea, but a mathematical necessity.
Toy Models
We use simplified "toy models" (like a pendulum) to illustrate the math.
Appendices
- Established (textbook): Modal analysis, wrench decomposition in screw theory, double-pendulum dynamics, block-matrix partitioning of equations of motion.
- Novel application: Applying these classical tools to derive explicit symbolic solutions for the drift/input decomposition in golf swing biomechanics.
- Notation note: This article uses the \((r, f)\) subscript convention for block matrix partitions. See the notation cross-reference section for mapping to other notation systems.
Recall the control-affine form \(\dot{x} = f(x) + G(x)u\) established in Part 1, where \(f\) encodes drift and \(G\) encodes control directions (Isidori 1995). The preceding parts of this manuscript have established the conceptual framework (Part I), the diagnostic tools (Part II), and the classification taxonomy (Part III) for the affine control interpretation of the golf swing.
This Part IV collects the detailed mathematical proofs that substantiate those assertions (Goldstein et al. 2002; Spong et al. 2005; Featherstone 2008). We explicitly derive the equations of motion for the coupled rigid–flexible system, demonstrate the preservation of the affine structure under modal truncation, and provide illustrative low-dimensional examples to build intuition.
Appendix a: Mathematical Derivations
The derivations in this appendix treat the grip boundary condition as a fixed mechanical constraint: the hand is clamped with constant stiffness and damping (the “Effective Plant” impedance). Physiologically, grip stiffness modulates with muscle activation, which is part of the input \(u\). If grip impedance were state-dependent or input-dependent, the boundary conditions defining the modal shapes would depend on \(u\), making the inertia matrix \(M(q,\eta,u)\) input-dependent and violating the affine structure. To preserve the validity of the control-affine form \(\dot{x} = f(x) + G(x)u\), we assume constant boundary conditions. This is justified during the ballistic phases of the swing (takeaway through impact), where grip adjustments are minimal; it may fail during the address or follow-through. All results in this Part IV assume this constant impedance context.
This appendix provides detailed derivations supporting the unified control-affine formulation. The goal is to show explicitly how the coupled rigid–flexible dynamics produce the affine decomposition
\[\dot{x} = f(x) + G(x)u\]
and to document all assumptions required for the formulation.
Coupled Rigid–Flexible Equations of Motion
We partition the generalized coordinates into
\[q \in \mathbb{R}^{n}, \qquad \eta \in \mathbb{R}^{m},\]
where \(q\) describes the rigid-body degrees of freedom and \(\eta\) contains the modal coordinates of the flexible shaft.
Let
\[v = \begin{bmatrix} \dot{q} \\ \dot{\eta} \end{bmatrix}\]
denote the generalized velocity.
The equations of motion for the coupled system may be written in compact form (Murray et al. 1994; Spong et al. 2005):
\[M(q,\eta)\,\dot{v} + C(q,\dot{q},\eta,\dot{\eta})\,v + G(q) + F_{s}(\eta,\dot{\eta}) = \begin{bmatrix} \tau \\[0.2em] 0 \end{bmatrix},\]
where:
- \(M(q,\eta)\) is the full inertia matrix of the rigid–flexible system,
- \(C(q,\dot{q},\eta,\dot{\eta})\) contains Coriolis and centrifugal terms,
- \(g(q)\) contains gravitational torques acting on the rigid-body DOFs,
- \(F_{s}(\eta,\dot{\eta}) = \begin{bmatrix} 0 \\ K_s \eta + C_s \dot{\eta} \end{bmatrix}\) contains shaft restoring and damping forces,
- \(\tau = B(q)\,u\) is the generalized torque vector induced by the input \(u\).
The flexible DOFs have no direct torque inputs.
Block Structure of the Inertia and Coriolis Matrices
The inertia matrix has the block form
\[M(q,\eta) = \begin{bmatrix} M_{rr}(q,\eta) & M_{rf}(q,\eta) \\ M_{fr}(q,\eta) & M_{ff}(q,\eta) \end{bmatrix},\]
where:
- \(M_{rr}\) is the rigid-body inertia matrix,
- \(M_{ff}\) is the modal inertia matrix for the flexible shaft,
- \(M_{rf}\) and \(M_{fr}\) represent rigid–flexible coupling.
This article uses the \((r, f)\) subscript convention (r = rigid, f = flexible). The companion article Affine Control Interpretation of the Golf Swing (affine-nature-golf-swing.qmd) uses the equivalent \((q, \eta)\) convention (q = joint coordinates, \(\eta\) = flexible modes): \(M_{rr} \equiv M_{qq}\), \(M_{ff} \equiv M_{\eta\eta}\), \(M_{rf} \equiv M_{q\eta}\), \(M_{fr} \equiv M_{\eta q}\).
Similarly, the Coriolis matrix has the block structure:
\[C(q,\dot{q},\eta,\dot{\eta}) = \begin{bmatrix} C_{rr} & C_{rf} \\ C_{fr} & C_{ff} \end{bmatrix},\]
with the usual property that
\[\dot{M} - 2C \quad \text{is skew-symmetric}.\]
These blocks arise from standard manipulator equations extended to include flexible-body contributions.
Solving for Accelerations
Writing \(\dot{v} = [\ddot{q}^{T}, \ddot{\eta}^{T}]^{T}\), the equations of motion become:
\[M \begin{bmatrix} \ddot{q} \\ \ddot{\eta} \end{bmatrix} = - C \begin{bmatrix} \dot{q} \\ \dot{\eta} \end{bmatrix} - G(q) - F_s(\eta,\dot{\eta}) + \begin{bmatrix} B(q)\,u \\ 0 \end{bmatrix}.\]
Assuming \(M\) is invertible (true for all admissible configurations of a well-defined multibody system), we solve for the generalized accelerations:
\[\begin{bmatrix} \ddot{q} \\ \ddot{\eta} \end{bmatrix} = - M^{-1} \left( C\,v + G + F_{s} \right) + M^{-1} \begin{bmatrix} B(q)\,u \\ 0 \end{bmatrix}.\]
Drift Acceleration
We define the drift acceleration as
\[a_{\mathrm{drift}}(x) = -M^{-1}(q,\eta) \left( C(q,\dot{q},\eta,\dot{\eta})\,v + G(q) + F_{s}(\eta,\dot{\eta}) \right).\]
All passive mechanical terms are explicitly included:
- inertia coupling,
- Coriolis and centripetal terms,
- gravity,
- elastic and damping forces from shaft deformation.
Importantly,
\[a_{\mathrm{drift}}(x) \quad \text{is independent of } u.\]
This constitutes the drift invariance property formalized in the main text.
Input Acceleration Mapping
The input-driven acceleration is the part proportional to \(u\). We can expand the inverse inertia matrix \(M^{-1}\) in block form. Let \(H = M^{-1}\). Its structure reflects the full coupling: \[H(q,\eta) = \begin{bmatrix} H_{qq} & H_{q\eta} \\ H_{\eta q} & H_{\eta\eta} \end{bmatrix}.\] The acceleration equation becomes: \[\begin{bmatrix} \ddot{q} \\ \ddot{\eta} \end{bmatrix}_{\text{input}} = \begin{bmatrix} H_{qq} & H_{q\eta} \\ H_{\eta q} & H_{\eta\eta} \end{bmatrix} \begin{bmatrix} B(q)u \\ 0 \end{bmatrix} = \begin{bmatrix} H_{qq} B(q) u \\ H_{\eta q} B(q) u \end{bmatrix}.\] Thus, the input torque \(u\) drives the rigid body accelerations \(\ddot{q}\) through the effective inertia \(H_{qq}\), and simultaneously drives the flexible mode accelerations \(\ddot{\eta}\) through the inertial coupling term \(H_{\eta q}\).
The input matrix is:
\[A_{\mathrm{input}}(x) = M^{-1}(q,\eta) \begin{bmatrix} B(q) \\ 0 \end{bmatrix} = \begin{bmatrix} H_{qq} B(q) \\ H_{\eta q} B(q) \end{bmatrix}.\]
Linearity of the torque input is preserved even in the presence of:
- flexible-body modes,
- rigid–flexible coupling,
- configuration-dependent \(B(q)\),
- velocity-dependent \(C\) terms.
First-Order State-Space (Affine) Form
Define the full state
\[x = \begin{bmatrix} q \\ \dot{q} \\ \eta \\ \dot{\eta} \end{bmatrix}.\]
The first-order dynamics follow directly from the definitions:
\[\dot{x} = \begin{bmatrix} \dot{q} \\ \ddot{q} \\ \dot{\eta} \\ \ddot{\eta} \end{bmatrix} = f(x) + G(x)u,\]
with
\[f(x) = \begin{bmatrix} \dot{q} \\[0.2em] \left[a_{\mathrm{drift}}(x)\right]_{1:n} \\[0.2em] \dot{\eta} \\[0.2em] \left[a_{\mathrm{drift}}(x)\right]_{n+1:n+m} \end{bmatrix},\]
\[G(x) = \begin{bmatrix} 0 \\[0.2em] \left[A_{\mathrm{input}}(x)\right]_{1:n,:} \\[0.2em] 0 \\[0.2em] \left[A_{\mathrm{input}}(x)\right]_{n+1:n+m,:} \end{bmatrix}.\]
This establishes the nonlinear control-affine structure used throughout the manuscript.
Appendix B: Modal Approximation for the Flexible Shaft
This appendix details the finite-dimensional modal approximation used to model shaft deformation in the unified rigid–flexible dynamics of the golfer–club system. The goal is to show how a continuous beam model reduces to the generalized modal coordinates \(\eta \in \mathbb{R}^{m}\) used in the main text, and why this representation preserves the linearity of the torque input in the final control-affine form.
Continuous Beam Model
The golf shaft is modeled as a uniform Euler–Bernoulli beam (Meirovitch 2010; Shabana 2020) of length \(L\), with transverse displacement field
\[w(s,t), \qquad s \in [0,L],\]
where \(s\) is the arc-length coordinate from the butt (grip) end. For small deflections, the beam equation is
\[\rho A \,\frac{\partial^2 w}{\partial t^2} + c_b \frac{\partial w}{\partial t} + EI \,\frac{\partial^4 w}{\partial s^4} = f_{\mathrm{base}}(s,t),\]
where:
- \(\rho A\) is the linear mass density,
- \(c_b\) is distributed damping,
- \(E I\) is the bending rigidity,
- \(f_{\mathrm{base}}(s,t)\) contains forces transmitted from the handle motion.
Boundary conditions depend on grip modeling. We assume:
\[w(0,t) = 0, \qquad w_s(0,t) = 0,\]
for a clamped handle, and free-tip conditions at \(s=L\):
\[EI\,w_{ss}(L,t) = 0, \qquad EI\,w_{sss}(L,t) = 0.\]
These constraints define the mode shapes used in the modal expansion.
Theoretically, the “clamped” boundary condition (\(w(0,t)=0, w_s(0,t)=0\)) implies an infinitely stiff grip. In reality, grip impedance varies with muscle activation, which is part of the input \(u\). If the boundary conditions depended on \(u\), the mode shapes \(\phi_i\) and consequently the mass matrix \(M\) would become input-dependent (\(M(u)\)), breaking the affine structure. To preserve the validity of the decomposition \(\dot{x} = f(x) + G(x)u\), we formally adopt the Constant Impedance Assumption: the grip is treated as a fixed mechanical constraint defining the “Effective Plant” for the swing. The input \(u\) is the torque applied to this plant, not a parameter that reshapes it.
Modal Expansion
The transverse displacement is approximated using a finite set of vibration modes (Meirovitch 2010):
\[w(s,t) \approx \sum_{i=1}^{m} \phi_i(s)\,\eta_i(t),\]
where:
- \(\phi_i(s)\) are the eigenfunctions of the beam operator with the boundary conditions above,
- \(\eta_i(t)\) are the modal amplitudes,
- \(m\) is the number of modes retained.
The mode shapes satisfy the eigenvalue problem:
\[EI \,\phi_i'''' = \omega_i^2 \rho A \,\phi_i,\]
with orthogonality relations
\[\int_0^L \rho A \,\phi_i(s)\,\phi_j(s)\,ds = \delta_{ij}.\]
Substituting the modal expansion into the beam equation and projecting onto each \(\phi_i\) yields the modal equations of motion.
In practice, how many modes \(m\) are required to capture shaft dynamics in the golf swing? For typical golf club shafts (length 0.9 m, hollow steel or composite), empirical studies show that the first 2–3 bending modes capture more than 95% of the shaft’s kinetic and strain energy during the swing. Higher modes have natural frequencies (20+ Hz) that are damped by grip compliance and rarely excited by the joint-driven excitation of the swing (which typically has energy content up to 10–15 Hz). For theoretical completeness, we retain the full modal series; in practice, truncating at \(m = 2\) or \(m = 3\) yields excellent fidelity while reducing computational cost.
Modal Equations of Motion
The projected dynamics have the form:
\[\ddot{\eta}_i + 2\,\zeta_i \omega_i\,\dot{\eta}_i + \omega_i^2 \eta_i = Q_i(q,\dot{q}),\]
for \(i = 1,\dots,m\), where:
- \(\omega_i\) is the natural frequency of mode \(i\),
- \(\zeta_i\) is the modal damping ratio (derived from \(c_b\)),
- \(Q_i(q,\dot{q})\) are generalized modal forces induced by handle motion.
Collecting terms,
\[\ddot{\eta} + C_s \dot{\eta} + K_s \eta = Q(q,\dot{q}),\]
where:
\[K_s = \mathrm{diag}(\omega_1^2,\dots,\omega_m^2),\qquad C_s = \mathrm{diag}(2\zeta_1\omega_1,\dots,2\zeta_m\omega_m).\]
These matrices match those used in the main text, where the shaft forces enter the drift term.
Preservation of Linear Torque Input
From the full equations of motion:
\[M(q,\eta)\,\dot{v} + C(q,\dot{q},\eta,\dot{\eta})\,v + G(q) + F_s(\eta,\dot{\eta}) = \begin{bmatrix} B(q)u \\ 0 \end{bmatrix},\]
the modal forces enter entirely through the drift term:
\[F_s(\eta,\dot{\eta}) = \begin{bmatrix} 0 \\ K_s \eta + C_s \dot{\eta} \end{bmatrix}.\]
Because \(u\) appears only in the block \(B(q)u\) associated with rigid-body torques:
\[u \mapsto \ddot{\eta} \quad \text{is filtered through} \quad \ddot{q}.\]
Hence: The modal representation does not create any nonlinear torque dependence.
All flexible-body dynamics remain embedded in the drift vector field \(f(x)\).
Appendix C: Pendulum Examples
This appendix presents simplified pendulum-based examples that illustrate the control-affine structure and the superposition of drift and input forces in explicit form. We treat:
- a planar (2D) pendulum with and without a flexible shaft.
The same affine structure extends to the spatial (3D) case — the Coriolis and gyroscopic terms become more complex, but the torque input still appears linearly — however the detailed 3D derivation is deferred to future work.
In each case, we show that the equations of motion can be written as
\[\dot{x} = f(x) + G(x)u,\]
with all flexible-shaft contributions confined to the drift term, preserving linearity in the torque input \(u\).
Planar Rigid Pendulum With Torque Input
Consider a simple planar pendulum of length \(L\) and mass \(m\), pivoted at the origin and moving in the vertical plane. Let \(\theta\) be the angle from the downward vertical, positive counterclockwise. A torque input \(u\) is applied at the pivot.
The kinetic and potential energy are:
\[T = \frac{1}{2} m L^2 \dot{\theta}^2, \qquad V = m g L (1 - \cos\theta),\]
and the Lagrangian is \(\mathcal{L} = T - V\).
The Euler–Lagrange equation (Goldstein et al. 2002) with generalized torque \(u\) is:
\[m L^2 \ddot{\theta} + m g L \sin\theta = u.\]
Defining the state \(x = [\theta, \dot{\theta}]^T\), we obtain the first-order system:
\[\dot{x} = \begin{bmatrix} \dot{\theta} \\ -\frac{g}{L} \sin\theta \end{bmatrix} + \begin{bmatrix} 0 \\ \frac{1}{m L^2} \end{bmatrix} u.\]
Thus,
\[f(x) = \begin{bmatrix} \dot{\theta} \\ -\frac{g}{L} \sin\theta \end{bmatrix}, \qquad G(x) = \begin{bmatrix} 0 \\ \frac{1}{m L^2} \end{bmatrix},\]
and the system is clearly control-affine.
Planar Pendulum With Flexible Shaft
We now attach a single bending mode to the pendulum to emulate a flexible shaft aligned with the pendulum rod. Let \(\eta\) be the modal coordinate representing transverse deflection of the shaft relative to the rigid rod.
The resulting equations of motion can be written schematically as:
\[M(\theta,\eta) \begin{bmatrix} \ddot{\theta} \\[0.2em] \ddot{\eta} \end{bmatrix} + C(\theta,\dot{\theta},\eta,\dot{\eta}) \begin{bmatrix} \dot{\theta} \\[0.2em] \dot{\eta} \end{bmatrix} + G(\theta,\eta) + \begin{bmatrix} 0 \\[0.2em] k_\eta \eta + c_\eta \dot{\eta} \end{bmatrix} = \begin{bmatrix} u \\[0.2em] 0 \end{bmatrix},\]
where \(k_\eta\) and \(c_\eta\) are the modal stiffness and damping, and all coupling terms have been absorbed into the matrices \(M,C,G\).
Solving for accelerations:
\[\begin{bmatrix} \ddot{\theta} \\[0.2em] \ddot{\eta} \end{bmatrix} = -M^{-1} \left( C \begin{bmatrix} \dot{\theta} \\ \dot{\eta} \end{bmatrix} + G + \begin{bmatrix} 0 \\ k_\eta \eta + c_\eta \dot{\eta} \end{bmatrix} \right) + M^{-1} \begin{bmatrix} 1 \\ 0 \end{bmatrix} u.\]
Defining
\[a_{\mathrm{drift}}(x) = -M^{-1} \left( C \begin{bmatrix} \dot{\theta} \\ \dot{\eta} \end{bmatrix} + G + \begin{bmatrix} 0 \\ k_\eta \eta + c_\eta \dot{\eta} \end{bmatrix} \right), \qquad A_{\mathrm{input}}(x) = M^{-1} \begin{bmatrix} 1 \\ 0 \end{bmatrix},\]
we obtain
\[\begin{bmatrix} \ddot{\theta} \\ \ddot{\eta} \end{bmatrix} = a_{\mathrm{drift}}(x) + A_{\mathrm{input}}(x)\,u.\]
Superposition still holds. Even though the dynamics are now richer and nonlinear in \((\theta,\eta)\), the torque input still appears linearly as \(A_{\mathrm{input}}(x)\,u\). The drift term contains all flexible-shaft effects; any two solutions corresponding to inputs \(u_1,u_2\) differ by the same linear map \(A_{\mathrm{input}}(x)\) applied to the difference in inputs.
To make the block structure reproducible, fix representative parameters: rod length \(L = 1.1\) m, rod/hand mass \(m_h = 0.40\) kg concentrated at the wrist, clubhead mass \(m_c = 0.200\) kg at the tip, and a single first-bending mode with modal mass \(m_\eta = 0.10\) kg, modal stiffness \(k_\eta = 2500\) N/m (first-mode frequency \(f_1 = \tfrac{1}{2\pi}\sqrt{k_\eta/m_\eta} \approx 25.2\) Hz — typical for a steel iron shaft), modal damping \(c_\eta = 0.8\) N·s/m, and shape function \(\Phi(L) = 1.0\) at the clubhead.
At the midsection configuration (\(\theta = -\tfrac{\pi}{4}\), \(\eta = 0\)) the inertia blocks evaluate to \[ M_{qq} = m_h L^2 + m_c L^2 + \text{mode contribution} \approx 0.726\ \mathrm{kg\,m^2}, \qquad M_{q\eta} = m_c L \cdot \Phi(L) \approx 0.220\ \mathrm{kg\,m}, \] \[ M_{\eta q} = M_{q\eta}, \qquad M_{\eta\eta} = m_\eta + m_c \Phi(L)^2 \approx 0.300\ \mathrm{kg}. \] The coupling ratio is then \[ \Gamma = - M_{\eta\eta}^{-1} M_{\eta q} \approx -\frac{0.220}{0.300} \approx -0.733\ \mathrm{m^{-1}} \] (the dimension is \(\mathrm{m^{-1}}\) here because \(\eta\) is a linear deflection while \(\theta\) is angular; in a model where \(\eta\) were itself angular, \(\Gamma\) would be dimensionless). The Schur complement / articulated-body inertia is \[ \Delta \;=\; M_{qq} - M_{q\eta} M_{\eta\eta}^{-1} M_{\eta q} \;\approx\; 0.726 - \frac{(0.220)^2}{0.300} \;\approx\; 0.565\ \mathrm{kg\,m^2}, \] so the flexible mode reduces the effective hand-swing inertia by roughly \(0.161/0.726 \approx 22\%\) — the “inertial relief” of the free flexible mode, exactly as claimed in Part 1 (Control-affine state space form). Verifying the sign of \(\Gamma\): apply unit torque \(u = 1\) N·m; hand acceleration is \(\ddot\theta = H_{qq} u = \Delta^{-1} \approx 1.77\ \mathrm{rad/s^2}\), and the modal acceleration is \(\ddot\eta_{\text{input}} = \Gamma H_{qq} u \approx -1.30\ \mathrm{m/s^2}\) — the shaft deflects backward (opposite the hand’s motion) as expected for inertial lag. A positive-\(\Gamma\) convention would predict forward-running deflection, which is physically wrong.
These numbers are reproducible with a 20-line Python script that builds \(M\), inverts it block-wise, and prints \(\Gamma, \Delta, H_{qq}\) — we leave the script as an exercise; the numerical values above were computed from the parameters listed at the top of this box.
Summary: Superposition and Affine Structure
Across both cases—planar rigid pendulum and planar flexible pendulum—the same pattern holds (and extends by the same reasoning to the spatial cases):
- Flexible dynamics (shaft modes) modify the drift term \(f(x)\) by adding state-dependent inertia, Coriolis, and restoring forces.
- The torque input \(u\) enters linearly through an input matrix \(G(x)\) that depends only on configuration and flexible coordinates, not on \(u\).
- For any fixed state \(x\), the acceleration (and hence \(\dot{x}\)) is an affine function of \(u\); the passive part and active part superpose exactly.
These examples provide explicit, low-dimensional realizations of the general control-affine structure used in the main text and highlight why the drift–input decomposition, ZTCF family (see canonical definitions), and ZVCF constructions remain valid and interpretable in the presence of shaft flexibility.
Transition to Numerical Validation (Part V)
The derivations in this appendix confirm the algebraic consistency of the AffineDrift framework. We have shown that the control-affine form \(\dot{x} = f(x) + G(x)u\) is not an approximation but a direct consequence of Lagrangian mechanics, preserved even under modal truncation of the flexible shaft.
However, algebraic exactness in continuous time does not guarantee numerical stability in discrete-time simulation. The golf swing involves stiff differential equations (high-frequency shaft modes) and rapid gyroscopic coupling. It is not immediately obvious that the delicate subtraction operations required for the ZTCF (\(u=0\)) and ZVCF (\(v=0\)) will remain robust under integration error.
Part V addresses this practical challenge (Hairer and Wanner 1996). It documents the Simulink forward-dynamics implementation used to validate these concepts, serving as a numerical proof-of-concept that bridges the gap between the derivations of Part IV and the empirical reality of the experimental applications.