2. Materials and Methods
2.1. Spatial Inertia Tensor and the Eigen-Screw Problem
The spatial inertia tensor [
8]
integrates both translational and rotational dynamics of a rigid body and is defined as:
where
m is the mass,
is the position vector from the body frame to the center of mass,
is the
rotational inertia tensor about the center of mass,
is the skew-symmetric matrix of
, and
is the
identity matrix.
The associated eigenvalue problem [
9] seeks screw vectors [
4]
, combining angular velocity
and linear velocity
, that satisfy:
The pitch
h of the screw is defined when
as:
This scalar represents the ratio of translation along the axis to rotation about it, encapsulating the geometric essence of the screw motion.
Furthermore, in the case of a finite twist, the pitch
can be derived from geometric components:
where
defines the screw axis and
is the raw moment vector [
4]. The component of
along
isolates the pitch:
Although multiple coordinate choices for are possible, the screw axis remains uniquely defined. If is not colinear with , the extracted pitch alters the final moment vector , thereby changing the screw. Therefore, for the screw to represent a unique line in Plücker coordinates, must be proportional to . This constraint ensures that the moment vector defines a consistent axis independent of parametric representation, as required in screw theory and spatial kinematics.
2.2. Harmonic Screws and One-Body Oscillation
According to Ball’s formulation of screw theory [
3], any wrench acting on a rigid body constrained to twist about a given screw
can be represented as a component of a resultant wrench projected onto that screw. The principle of *reciprocity* states that for equilibrium to be maintained, the applied wrench must be reciprocal to the screw of allowable motion.
Formally, consider a system of wrenches acting on a body free to twist only about a given screw
. The total virtual work done by the wrenches during an infinitesimal twist about
must be zero:
which simplifies to:
Here, represents the intensity of the ith wrench, and denotes the corresponding screw coordinate projection onto . This condition ensures that the net effect of the applied wrenches does not disturb the constrained motion.
From this principle, Ball further notes that a given wrench can always be equivalently replaced by another wrench acting along a different screw, provided the body is constrained to twist only about . This leads to the concept of a harmonic screw, whereby the applied and reactive wrenches achieve a dynamic balance across a screw system, allowing the body to undergo pure oscillation.
This harmonic relationship is foundational for modeling single-body oscillations—such as those seen in sports movements—where the effective action results from rhythmic alternation of forces applied on reciprocal screw systems [
10]. The harmonic screw thus describes a self-consistent twist-wrench pairing that resonates with the mechanical affordances of the body and its constraints.
As Ball observed, this principle bears strong analogy to the classical condition for equilibrium of a particle constrained to a line under the action of multiple forces. If
P and
Q are two such forces acting at angles
l and
m to the direction of allowable motion, then:
This analogy emphasizes that in both particle and screw systems, projected force components must cancel in the direction of allowable motion. The harmonic screw generalizes this idea to spatial rigid-body systems.
In the present study, we exploit this harmonic screw formulation to identify oscillatory modes in biomechanical motion—especially where the body pivots or twists rhythmically about a constrained axis. These modes provide insight into invariant biomechanical structures underlying efficient movement.
2.3. Instrumentation and Computational Analysis
This study leveraged a previously validated biomechanical dataset [
8] that recorded full-body golf swing dynamics of two female participants representing contrasting levels of expertise (
Table 1).
The selection of this dataset was motivated by its comprehensive multi-modal instrumentation, which included high-speed motion capture synchronized with ground reaction force (GRF) data, providing a robust platform for advanced mechanical modeling.
The original data acquisition was conducted using a 12-camera Qualisys optoelectronic motion capture system (model: Oqus-300, Qualisys AB, Gothenburg, Sweden) operating at a sampling frequency of 300 Hz. Twenty-four retroreflective markers, in conjunction with four rigid-body clusters, were affixed to key anatomical landmarks following ISAK anthropometric protocols. This configuration enabled accurate 3D tracking of major body segments (e.g., head, thorax, pelvis, upper and lower limbs) and the golf club. Simultaneously, kinetic data were acquired via a Kistler force platform, aligned with the lead foot to define the global reference frame based on the initial center of pressure (COP).
Figure 2 illustrates the marker arrangement at the address phase, highlighting the detailed setup used for wrist and club tracking. These features were essential for subsequent analysis of grip torque dynamics and transmission of force impulses through the kinetic chain.
To focus the analysis on the most critical biomechanical phase of the swing, only the downswing portion—from the initiation of forward acceleration to just before club–ball contact—was examined. Using segment-based local coordinate systems, we computed the instantaneous screw axes (ISAs) of the club relative to the trunk, as well as the corresponding angular velocities and linear displacements.
The ISA trajectories [
11] were interpreted as functional representations of the club’s inertial coordination with upper-body rotation. Using a screw-theoretic framework, the twist motion of the club was decomposed into time-varying pitch and orientation parameters, allowing us to characterize the dynamic coupling between translational and rotational components. This formulation enabled identification of control strategies used by each participant and provided a biomechanically grounded measure of skill-dependent swing organization.
This methodological choice to reuse an existing, high-quality dataset ensured internal consistency while reducing inter-session variability. It also preserved the ecological fidelity of the motor behavior, making it possible to examine how perception–action couplings unfold in real-world performance contexts. This approach is consistent with the theoretical framework adopted in our prior research on symmetry and motor optimization in skilled action.
2.4. Computation of Instantaneous Screw Axis and Principal Directions of Inertia
To identify the instantaneous screw axis (ISA) and spatial principal direction of inertia during the downswing, we used a screw-theoretic transformation framework based on rigid-body marker data. The motion segment (e.g., lower limb or club) was modeled as a rigid body defined by six non-collinear reflective markers tracked over time.
Let
and
denote the 3D coordinates of the
k-th marker at two consecutive frames. The relative motion between time frames was computed using singular value decomposition (SVD) applied to centered marker configurations. Let
and
represent the positions at
and
t, respectively. After removing the centroid offset, we computed:
where
and
denote mean marker positions. The rotation matrix
was extracted via:
where
ensures proper rotation. The translational shift
was derived from the mean displacement of marker centroids.
To incorporate spatial inertia, we mapped the original inertia tensor
into the moving frame via a screw transformation
constructed as:
The eigensystem of the
rotational inertia submatrix
of the transformed inertia tensor was then computed to obtain the principal directions of rotation:
To extract the eigen-screw, the cross-coupling term
was projected onto each eigenvector
, giving rise to the screw moment:
Among these, the direction corresponding to the minimum eigenvalue of the inverse inertia matrix was selected as the dominant principal screw, representing the minimal-energy configuration for rotational-translational coupling.
This computation was repeated for each frame following downswing initiation (). The dominant screw axes were normalized and written to an output array, representing the time-evolving primary ISA during the downswing. These were exported as the matrix `Sprime.xlsx` for subsequent analysis and visualization.
The screw-based principal directions provide a biomechanically meaningful axis of movement coordination, invariant to translation and frame shifts, and consistent with screw theory’s geometric formalism.