## Grip Impedance and the Guitar Analogy {#sec-impact-grip-impedance}

A loosely resting finger on a guitar string can absorb energy, alter a boundary,
and suppress modes that have motion at the finger location. Touching near a node
can instead leave selected harmonics relatively intact. The useful analogy is
**mode-dependent coupling to a dissipative contact**, not a rule that a looser
golf grip always sounds worse or a tighter grip always sounds better.

A club is a distributed shaft joined to a complex head. A hand near the grip may
strongly affect low shaft modes but weakly affect a localized face or crown mode.
It can alter damping, effective inertia, stiffness, contact area, and slipping
simultaneously. A clubhead is also not an ideal string: bending, torsion, shell
motion and acoustic radiation have different spatial patterns.

For a stationary support, represent the linearized grip by a driving-point
mechanical impedance $\mathbf Z_g(\omega)$ relating the wrench into the
boundary to the velocity at its port. A matrix is needed for translations, rotations and their coupling.
For a passive boundary its Hermitian real part is nonnegative on the represented
ports. In a scalar spring–damper–mass example, using $e^{i\omega t}$,

$$
Z_g(\omega)=c_g+i\omega m_g+\frac{k_g}{i\omega}.
$$

This is an illustrative network, not the impedance of a measured hand. Identify
force and velocity at the same port, preserve phase, and distinguish active
muscle work from passive dissipation. A moving prescribed grip can do work even
if its displacement during the collision is small.

For the following frozen linear reduction, let $\boldsymbol\eta$ be small
structural displacements and $\mathbf v_g=\mathbf L_g\dot{\boldsymbol\eta}$
the grip-port velocity relative to a stationary support. Express the port's
force, moment, linear velocity and angular velocity at one declared origin and
in compatible axes. If $\mathbf w_{\rm into}=\mathbf Z_g\mathbf v_g$ is the
wrench exerted by the club on the boundary, its reaction on the structure is
$-\mathbf L_g^T\mathbf w_{\rm into}$. The transpose follows from power:
$\dot{\boldsymbol\eta}^T\mathbf L_g^T\mathbf w_{\rm into}
=\mathbf v_g^T\mathbf w_{\rm into}$. Eliminating this boundary therefore adds $i\omega\mathbf L_g^T\mathbf Z_g\mathbf L_g$ to dynamic stiffness. The
contact mobility, in velocity per force, is

$$
\begin{aligned}
\mathbf Y_c(\omega)&=i\omega\mathbf B_c^T\mathbf D(\omega)^{-1}\mathbf B_c,\\
\mathbf D(\omega)&=-\omega^2\mathbf M+i\omega(\mathbf C+\mathbf G)\\
&\quad+\mathbf K+i\omega\mathbf L_g^T\mathbf Z_g\mathbf L_g.
\end{aligned}
$$ {#eq-impact-frequency-mobility}

Do not convert this complex, frequency-dependent object into one real
"effective mass" without specifying frequency and identification convention.
Resonances, antiresonances and dissipative phase are central to the question.
The frozen-frequency description also needs comparison against the time-varying
swing when rotation, tension or contact topology changes appreciably.
The displayed inverse is only defined away from poles of that frozen operator;
use a compatible causal response and initial conditions for a transient.
Do not include the same hand mass, stiffness or damping both in the structural
matrices and again in an impedance that already condenses those dynamics.

### Why a Moving Support Needs Another Input {#sec-impact-moving-grip}

A fixed-support driving-point impedance does not by itself specify the response
to an accelerating support. Consider an illustrative one-dimensional boundary:
club-port displacement $x$, hand-mass displacement $y$, support displacement $b$,
a grip spring/damper $(k_g,c_g)$ between $x$ and $y$, and another $(k_b,c_b)$
between $y$ and $b$. These are small motions about a stated equilibrium, not
identified physiological parameters. The force exerted by the club on this
boundary is $F_g=k_g(x-y)+c_g(\dot x-\dot y)$; the hand mass obeys
$m_g\ddot y=F_g-k_b(y-b)-c_b(\dot y-\dot b)$.

With $D_g=k_g+i\omega c_g$, $D_b=k_b+i\omega c_b$, and
$D_h=D_g+D_b-m_g\omega^2\ne0$, elimination gives

$$
\begin{aligned}
\widehat F_g&=D_{cc}\widehat x+D_{cb}\widehat b,\\
D_{cc}&=D_g\left(1-\frac{D_g}{D_h}\right),\\
D_{cb}&=-\frac{D_gD_b}{D_h},\\
Z_{cc}&=\frac{D_{cc}}{i\omega},\qquad \omega\ne0.
\end{aligned}
$$ {#eq-impact-moving-boundary}

The fixed-support impedance is $Z_{cc}$, but the second coefficient is generally
not $-D_{cc}$ when $m_g\ne0$. Replacing $\dot x$ by $\dot x-\dot b$ in that
one-port law would therefore lose the support-acceleration response. The
massless two-terminal spring/damper reduction is a special case. Even
$\widehat x=\widehat b$ can require force to accelerate the intervening mass.
A moving hand/support model needs its additional input or retained states,
and the energy audit must include work at both driven ports. The force on the
club is $-F_g$; positive force alone does not determine its sign of work.
This example explains a possible pathway without assigning it a golfer-specific
size or proving an acoustic benefit.


Chiementin and colleagues identified grip-dependent modes of one club with one
right-handed participant. Their bending-mode table and discussion generally
show stronger damping under firmer grip, whereas their conclusion says damping
decreases. That inconsistency must be retained in the evidence ledger, not
silently converted into a monotone physiological law. Their sub-500 Hz modal
study does not qualify a broadband driver-sound model [@chiementin2019grip].

## Contact-Force Regularity and Acoustic Bandwidth {#sec-impact-force-regularity}

A contact model can predict ball impulse adequately while misrepresenting the
force waveform that excites audible modes. Consider the unilateral clipped
Kelvin–Voigt baseline: for compression $\delta>0$ its force is
$F_n=\max(0,k_n\delta+c_n\dot\delta)$, with zero force outside contact.
At a closing first touch with finite incoming speed $v_0>0$, its incoming
force limit is $c_n v_0$. Assigning zero force at the single instant
$\delta=0$ does not remove this right-hand limit.

This finite force jump is **not a Dirac impulse**. Momentum and accumulated
work can remain continuous across first touch. Its spectral consequences
nevertheless differ from those of a continuous-onset force. With
$\widehat F(\omega)=\int F(t)e^{-i\omega t}\,dt$, a compactly supported,
piecewise smooth force with finitely many jumps $\Delta F_j$ satisfies

$$
\begin{aligned}
i\omega\widehat F(\omega)&=\sum_j\Delta F_j e^{-i\omega t_j}\\
&\quad+\int F'_{\mathrm{regular}}(t)e^{-i\omega t}\,dt.
\end{aligned}
$$ {#eq-impact-force-jump-spectrum}

This follows by integrating the distributional derivative. If the regular
derivative is integrable, its transform tends to zero, leaving an oscillatory
$1/\omega$ leading contribution from the jumps. Cancellation can create spectral
zeros; this is not a monotone amplitude law or a pointwise lower bound. Units
agree: $\widehat F$ and $\Delta F/\omega$ both have units N s. For a continuous
force whose first derivative has the corresponding bounded-variation and
integrability properties, a second integration gives an $O(\omega^{-2})$ bound.
Other events and the analysis-window endpoints can determine the overall tail,
so smooth onset alone does not qualify the full pulse spectrum.

Equal impulse or restitution therefore does not imply equal modal excitation.
An equal-area rectangular pulse and triangular pulse provide a simple
analytical comparison: the former has force jumps, whereas the latter has
continuous force and slope jumps. No player, club material or perception model
is needed for that distinction. It motivates testing force-history bandwidth
in addition to final ball speed, spin and work.

A continuous-onset candidate is $F_n=\delta^n(k+\lambda\dot\delta)$ for
$\delta>0$ and $n>0$. Simbody documents a sphere/sphere-plane Hunt–Crossley
implementation with $n=3/2$, $\lambda=(3/2)kc$ and recoverable storage
$U=(2/5)k\delta^{5/2}$; its restitution-versus-speed parameterization is a
low-speed approximation. These geometry, parameter and speed conventions are
part of the model, not measured properties of a multilayer golf ball and
flexing clubface [@simbody2019huntcrossley].

Continuity does not guarantee nonadhesion under arbitrary external forcing:
$k+\lambda\dot\delta$ can become negative during sufficiently rapid unloading.
Carvalho and Martins report this limitation and propose extensions and an
exponential contact law. Their institutional abstract supports that caution;
it does not supply a fully inspected derivation for adopting an exact
restitution mapping here [@carvalho2019restitution]. A free, frictionless
two-body calibration must not silently become a guarantee for a rotating,
externally loaded shaft/grip system. Clipping a candidate force also requires
an explicit storage and dissipation audit.

Before comparing acoustic outcomes, resolve first touch, unloading and
separation, and establish separate convergence of impulse, peaks, waveform
amplitude and phase over a declared frequency band. Retain the force-law
identity, contact patch, time resolution, sampling, filtering and window.
Windowing can suppress a numerical artifact without validating its physical
origin. A smoother candidate law needs held-out force measurements; a more
pleasant simulated tone is not evidence that its mechanics are more accurate.

The [source-identified coupled-contact study](https://github.com/D-sorganization/Tools/blob/279926e95e8a17cc42702a53c30ea03bf4638bec/docs/development/impact-acoustics/FRICTION_ENTRY.md)
resolves first-touch and sliding-entry reference events, while explicitly
retaining unresolved production event-history and cutoff-work limits.
Its [regularity research note](https://github.com/D-sorganization/Tools/blob/279926e95e8a17cc42702a53c30ea03bf4638bec/docs/development/impact-acoustics/CONTACT_FORCE_REGULARITY.md)
connects these findings to the contact and acoustic epics. These are bounded
synthetic results and analytical requirements. A force spectrum still needs
structural and acoustic transfer functions before it becomes pressure, and
blinded listener evidence before it supports a sweetness claim.

## From Impact to Radiated Sound {#sec-impact-acoustic-radiation}

The modeling chain is contact traction $\rightarrow$ structural motion
$\rightarrow$ air pressure $\rightarrow$ microphone/ear $\rightarrow$ perception.
A force spectrum alone is not sound, and a shaft accelerometer is not a
microphone. An arbitrary audio oscillator driven by ball speed would be
sonification, not a prediction of impact acoustics.

For a mass-normalized uncoupled structural mode $j$,

$$
\begin{aligned}
\ddot\eta_j+2\zeta_j\omega_j\dot\eta_j+\omega_j^2\eta_j&=Q_j(t),\\
Q_j(t)&=\int_{S_c}\boldsymbol\phi_j(\mathbf x)^T
\mathbf t_c(\mathbf x,t)\,dS.
\end{aligned}
$$

Strike position changes modal participation through the mode shape
$\boldsymbol\phi_j$ and the contact patch. Contact duration and force shape
change the excitation spectrum. At nonzero initial modal amplitude or velocity,
the response is the sum of the homogeneous motion and forced response; their
phase can reinforce or cancel particular components. A free-head impulse model
omits these states.

For fixed coefficients,

$$
\widehat\eta_j(\omega)=
\frac{\widehat Q_j(\omega)}{\omega_j^2-\omega^2+2i\zeta_j\omega_j\omega}
$$

describes the forced part. Coupled modes require the matrix operator instead.
Continue the calculation after ball separation: much of the audible ringdown
occurs when no ball–face force remains. Changes to the boundary after separation
can change later sound while leaving ball launch unchanged.

For an exterior linear acoustic field with sound speed $c_a$, air density
$\rho_a$ and normal surface velocity $v_n$, a frequency-domain model satisfies

$$
\begin{aligned}
(\nabla^2+k_a^2)\widehat p&=0,\\
\frac{\partial\widehat p}{\partial n}&=-i\omega\rho_a\widehat v_n,\\
k_a&=\omega/c_a,
\end{aligned}
$$

with $n$ pointing from the solid into the air and an outgoing-wave condition
consistent with $e^{i\omega t}$. A boundary-element or qualified reduced
radiation model maps the entire vibrating surface to observer pressure. Thin
plates can be useful validation fixtures; the baffled-plate Rayleigh integral
must not be applied to a hollow driver without checking its geometric and
wavelength assumptions. Include cavity openings or internal air coupling where
they materially affect the response.

Mase and colleagues provide a golf-specific finite-element/boundary-element
modeling precedent. Delaye and colleagues use a titanium plate as an intermediate
acoustic validation fixture [@mase2012sound; @delaye2016plate]. The plate study's
microphone was not calibrated for pressure: it compared relative spectral
amplitudes and decay, with substantial discrepancies, rather than validating
absolute radiated pressure. Its material parameters were also adjusted against
the plate, so transfer to a new club remains a prediction to test. Earlier thesis work
also records poor correlation for a driver despite success on a simpler plate;
this is an explicit warning against treating a plausible simulated tone as
validation [@sharpe2010sound].

A useful measured transfer model is

$$
\begin{aligned}
\widehat p_{\rm mic}(\omega)&=H_{\rm mic}(\omega)H_{\rm room}(\omega)\\
&\quad\times\sum_j R_j(\omega,\mathbf x_{\rm mic})\widehat\eta_j(\omega)\\
&\quad+\widehat p_{\rm other}(\omega).
\end{aligned}
$$

Here $R_j$ includes radiation and propagation, with units consistent with modal
normalization. The other-source term can contain ball vibration, air motion,
turf/tee contact and background noise. Source movement changes distance,
directivity and arrival time; evaluate Doppler and moving-source corrections
where the required accuracy warrants them. One microphone cannot generally
separate all these sources without additional measurements or assumptions.

Report pressure in pascals and calibrated sound levels, analysis window, sample
rate, filtering and microphone geometry. Peak pressure, sound exposure, spectral
centroid, band energy, decay times and modal frequencies describe different
features. Psychoacoustic loudness and sharpness require their specified
algorithms; an FFT amplitude or spectral centroid is not automatically either
quantity. Sound-pressure level uses the declared RMS window and reference
pressure, conventionally 20 micropascals in air.

## Sound Quality Is a Perception Result {#sec-impact-sweetness}

Roberts and colleagues found correlations between impact sound and perceived
feel in elite golfers. A companion study examined vibration and feel. Together
they justify measuring auditory and tactile responses separately; they do not
establish that one player's body mechanically adds mass to the strike
[@roberts2005sound; @roberts2005vibration].

For the reported player-to-player sound difference, first record repeated shots
with the same head, shaft, grip and ball model. Measure actual delivery and face
position instead of assuming that similar-looking swings are matched. Record
grip forces, shaft strain and head/shaft vibration where feasible; their absence
limits mechanism attribution.

Use two complementary listening tests: calibrated original-level playback to
measure the complete audible difference, and **level-matched** playback to
investigate timbre beyond loudness. Randomize identity and order, mask player
skill and ball-flight outcomes, balance headphone/speaker presentation, and
test repeatability. Crossed or swapped audio/visual cues can expose expectation
effects. Audio-only, vibration-only and combined presentations answer different
questions. Define the population and ratings before data collection.

A **matched-state** experiment fixes measured rigid delivery and then varies a
declared elastic or boundary state. A matched-input experiment applies the same
control program and allows delivery to change. Neither design is achieved by
matching clubhead speed alone. Statistical matching cannot remove confounding
from unmeasured contact location, equipment differences or elastic state.

If observers distinguish players after all measured launch variables are
controlled, that is evidence of a residual auditory difference, not yet proof
of grip damping. Establish the mechanism by a reproducible physical
intervention that changes the predicted modes or transfer response and produces
the corresponding measured acoustic change.