Inside the Racket: Modelling Frame Vibration with MATLAB and ANSYS

Inside the Racket: Modelling Frame Vibration with MATLAB and ANSYS

πŸ’» The first open parametric framework for tennis racket frame vibration analysis.

At a glance

Tools: MATLAB Β· ANSYS Β· Finite Element Analysis

Skills Demonstrated: Parametric Modelling, Quantitative Data Analysis, Validation & Testing, Python/MATLAB scripting, Structured Problem Solving, Technical Communication

Result: Validated to <0.5% against published experimental data across three independent benchmarks β€” deformation, frequency, and stress

The Result

Four findings with direct design implications emerged from the 13-configuration study. Each one maps a specific geometric parameter to a measurable structural outcome.

πŸ“ 383.6 Hz simulated. 382 Hz experimental. Error: <0.5%.

Validated against Russell (2014), Banwell (2014), and Kotze (2000) β€” three independent published benchmarks, three aligned results.

The Problem

Most published FEA studies of tennis rackets use a single fixed geometry β€” one specific racket, modelled by hand. That answers "how does this racket vibrate?" but not the more useful question: how does racket shape drive vibrational behaviour?

You can't compare what you can't vary. You can't optimise what you can't compare.

When I searched the literature, no open parametric geometry framework existed. That gap defined the project.

Why This Matters for Racket Pulse

Racket Pulse's premise β€” that a sensor can read racket speed, spin, and impact point from frame vibration β€” only holds if vibration signatures are measurable, repeatable, and distinguishable across different impact conditions.

This dissertation provided that proof.

The validation against Russell, Banwell, and Kotze quantified what Racket Pulse needed to assume: the physics is real, the signatures are reproducible, and a sensor at the throat would see a clean, structured signal.

β†’ Read the Racket Pulse case study

The Parametric Framework

Six input parameters. One complete racket geometry, ready for FEA.

Parameter
Baseline
Range tested
Unit
Head size
100
85 – 115
inΒ²
Aspect ratio
1.35
1.20 – 1.55
β€”
Total length
27
26 – 28
in
Beam width
23
18 – 28
mm
Beam height
26
20 – 32
mm
Wall thickness
2.5
1.5 – 3.5
mm

Three approaches were considered: manual CAD in SpaceClaim, importing manufacturer drawings, or parametric MATLAB scripting. Manual CAD was rejected because it locks the study to a single geometry. Manufacturer drawings were unavailable. Parametric scripting won: without it, only one racket could have been studied.

The script went through three iterations. The first produced symmetrical elliptical curves with no control over throat shape. After reviewing Grant et al. (2022, 2024), independent variables were exposed for head width, aspect ratio, beam dimensions and shaft length, at which point the geometries began matching real racket proportions.

A Quick Look β€” Interactive Geometry Explorer

Input

image

Adjust any of the six parameters and the frame outline updates in real time using the same parametric equations behind the MATLAB script. The baseline configuration matches the racket used to seed the FEA study.

Output

image
image
image

The Pipeline

  1. Generate geometry in MATLAB from six input parameters β†’ export as DXF
  2. Import into ANSYS SpaceClaim β†’ convert 2D outline to 3D surface body using beam height and wall thickness
  3. Apply orthotropic graphite/epoxy material properties (Taraborrelli et al., 2019)
  4. Apply clamped-handle boundary condition simulating a player's grip
  5. Run mesh convergence study β€” early runs showed >5% frequency variation between mesh refinements; converged element size adopted and held constant across the full study
  6. Solve modal analysis β†’ extract first six natural frequencies and mode shapes
  7. Extract regional stress data at throat, hoop, shaft, and handle
⚠️ The mesh convergence step was unplanned. Early results varied by more than 5% between mesh refinements. Recognising that and responding rather than accepting the default output was one of the more important engineering calls in the project.

Key Findings

Validation

image

Three independent validation points all aligned:

Metric
Simulated
Published
Source
Error
Mode V torsional frequency
383.6 Hz
382 Hz
Russell (2014)
<0.5%
Mode III bending frequency
185.4 Hz
100–200 Hz range
Banwell (2014)
Within range
Throat junction peak stress
41,180 Pa
41,180 Pa
Kotze (2000)
0%

Throat Stress is Localised and Asymmetric

image

The fundamental bending mode loaded one side of the throat substantially more than the other two struts. This asymmetry does not appear in any of the literature reviewed at project start. It emerged from the simulation and prompted a new discussion section on throat geometry as a critical design variable β€” one not in the original scope.

Parameter Sensitivities (13-configuration OFAT study)

Parameter
Finding
Frame length
+4% length β†’ 137% increase in Z-deformation. Consistent with classical LΒ³ scaling.
Head size
85 inΒ² β†’ 1,084% spike in Z-deformation. Torsional stiffness effectively collapses. Diminishing returns above baseline.
Beam height
32 mm β†’ 44,703 Pa peak throat stress. Primary lever for managing stress concentration.
Frame length (frequency)
Extended frames push torsional response toward the 80–200 Hz forearm loading band Hennig (2007) linked to overuse injury.

What I'd Do Next

  • String-bed modelling: using LINK180 line elements to capture string-frame coupling, currently abstracted out
  • Dynamic transient FEA: to capture inertial effects during the 4–6 ms ball-impact window
  • Experimental validation: via Digital Image Correlation (DIC) and accelerometry on physical rackets
  • Optimisation coupling: wrap the framework in a loop to minimise throat stress for a target frame mass, turning the analysis tool into a design tool

Four Things I Learned

  1. Parametric thinking opens the door to comparative analysis. Without the flexible geometry model, only one racket could have been studied and the conclusions would have been correspondingly narrower.
  2. Trust validation, not assumption. The unplanned mesh convergence study was the difference between a credible result and a flawed one.
  3. The most interesting results are the ones you didn't predict. The asymmetric throat stress wasn't in any source literature. Looking carefully at unexpected outputs rather than treating them as noise is a discipline I'll carry into every analytical role.
  4. A small original contribution is worth more than a polished replication. Reproducing a single known result would have been easier. Building the parametric framework produced something nobody else had.

Links

← Back to portfolio