1. Introduction
Lawn Tennis (tennis) is widely considered to have been ‘invented’ in 1874, when Major Wingfield patented the game and started marketing sets for playing on suitable grass [
1]. Early tennis rackets were similar to those used in Real Tennis, and were typically asymmetrical or ‘lopsided’ [
2]. The ball bounces low in Real Tennis, and a racket with a lopsided head is desirable to make it easier to bring the hitting surface close to the ground. These lopsided rackets disappeared from tennis as the game developed, with symmetrical frames common by the end of the nineteenth century. Most early tennis rackets were wooden, with incremental developments until the 1960s [
2,
3,
4]. Tennis became more competitive with the introduction of the Open Era in 1968, when professionals and amateurs started competing together for prize money. In 1968, the Wimbledon Championships prize money for Gentlemen’s Singles and Ladies’ Singles winners were £2000 and £750 (~£30 k and £11 k value today), respectively, compared to over two million in 2019 [
5]. Increased competition, larger prize funds and demand from players helped drive the development of the tennis racket. With better rackets, tennis developed from a game of serve-and-volley to baseline play, with the fastest male and female servers exceeding ball speeds of 55 m/s (>200 km/h).
Engineers explored aluminium and fibre–polymer composites as alternatives to wood in the 1970s, as they experimented with new racket shapes [
2,
3,
4]. In 1976, Howard Head patented a frame with an ‘oversized head’ [
6] that made it easier to play and laid foundations for the modern racket. In response to these developments, the International Tennis Federation (ITF) began specifying rules for the racket in 1978, to protect the nature of the game and ensure that the player was the primary determinant of a match outcome. Racket size was first limited in 1981, with current limits of 737 mm for overall length, 292 mm for width, and 394 mm for hitting surface length [
7]. Most modern rackets are now made from fibre–polymer composites, which offer high specific modulus and manufacturing versatility, providing the engineer with more freedom over parameters such as the mass distribution and stiffness [
4].
There are many publications on the mechanics of tennis and the role of the racket (e.g., [
8,
9,
10,
11,
12,
13,
14,
15]), with reviews highlighting the importance of considering player-racket interactions [
16,
17]. Racket inertial properties are particularly important, as they influence the player’s ability to accelerate the racket through the stroke [
18,
19,
20,
21]. A better knowledge of how the racket has developed over time could improve our understanding of its contribution to player performance, with implications for product development and regulation [
7], injury prevention strategies [
9,
22,
23], and exercise promotion initiatives [
24], as well as spectator experience, and education purposes [
25,
26,
27]. Haake et al. [
3] published the most comprehensive work to date on the historical development of the tennis racket, by characterising 150 rackets from the 1870s to 2007 and illustrating how their dimensions, inertial properties and natural frequencies have changed. This work goes further than Haake et al. [
3], by undertaking a wider range of measurements on a larger and more diverse group of racket dating from 1874 to 2017, while also reporting on materials.
2. Materials and Methods
Data from 525 rackets from 1874 to 2017 were collected and recorded in Excel
® (Microsoft, Redmond, WA, USA) (see
Table S1 in the Supplementary Material) and analysed in MATLAB
® (Matworks, Natick, MA, USA). One hundred and eight rackets were characterised at the Wimbledon Lawn Tennis Museum, where it was possible to access old, rare and valuable samples, and one was from Manchester Metropolitan University. Measurements (and other records) of the other 416 rackets came from prior work [
28,
29]. Of the 525 rackets, 417 (79.5%) were from the museum, 91 (17%) were from a brand’s headquarters [
28], 13 (2.5%) were from the university and four (1%) were from the ITF.
Racket dating was principally carried out using a book [
2] and the museum catalogue. If a racket featured in the book [
2] and a range of less than ten years was given (135 (26%) rackets), the earliest date (year of release) was used. In 323 (62%) cases, dating was carried out using the museum catalogue, mainly (282 cases) because the racket was not in the book [
2]. For racket models that were manufactured for more than nine years (17 models (e.g., Demon, Slazenger; Tournament, Spalding; Maxply, Dunlop) and 41 (8%) rackets), dates were taken from the catalogue to give the best estimation of the sample, as the design may have changed over time. The book [
2] did not include rackets postdating 1990, and dates for 67 (13%) rackets (from the work reported in [
28]) were obtained by consulting the manufacturer and websites [
30,
31]. Some rackets may have been dated incorrectly (the catalogue stated “circa” before the date for 271 (52%) of the rackets), with errors estimated to be within five years in most instances, which was deemed acceptable given that >500 rackets spanning >140 years were characterised.
Visual inspection identified racket materials, with cross checking against the book [
2], catalogue, websites and graphics. No distinction was made between types of wood (ash, maple, mahogany, beech etc.), grades of metal nor the constituent materials of a fibre–polymer composite (epoxy, glass, graphene, carbon, aramid etc.), with rackets categorised as either, wood, steel, aluminium or composite, or a combination of two of these. Based on recommendations from previous work [
28], string materials were not considered as it is challenging to identify them from visual inspection and the original strings may have been replaced.
Figure 1a–c defines the racket measurements and
Table 1 summarises the methods, which are detailed in [
29]. Each racket was photographed, as were any distinguishing features (
Figure 1d–g) (relative to a conventional modern racket) such as holes in the frame, extended string bed/no yoke, asymmetry, adjustable string patterns and unconventional yokes/throat regions. These distinguishing features were grouped in four categories: Head features, throat features, handle features and stringing patterns. For older rackets with wooden handles, grip length was identified using any grooves or engravings, a different wood/colour or polished coverings on the handle, as illustrated in [
29]. A scale (Smart Weigh SWS1KG Elite Series) and support were used to locate the racket’s centre of mass (CoM) [
29]. The racket was held horizontal with one end on the scale and the other the support, with the CoM location obtained from the product of its length and the ratio of the scale reading to the total mass. The method was compared to the more common approach of locating the ‘balance point’ for five rackets (1932 to 1987 with CoM locations from 321 to 335 mm), whereby the racket is balanced on a narrow support under the CoM [
3,
28]: the largest difference between the two methods was 2 mm (<1%).
Racket frame depth (
Figure 1c) was estimated as the mean of the minimum and maximum depth measurements. Racket frame thickness (
Figure 1c) was estimated as half the difference between the external and internal head width
(
Figure 1b). Head length was not used to estimate frame thickness, because some rackets (e.g., Dayton Steel Racket Corporation) had no yoke with strings connected to the handle. Racket head area (h
a) was approximated as an ellipse (
Figure 1b) using Equation (1),
where
is the internal head length. To account for changes in head size when investigating the number of strings, string bed density (s
d) was calculated using Equation (2).
where
and
are the number of cross and main strings, respectively.
Transverse moment of inertia (MOI) was defined as the MOI acting about a lateral in-plane axis passing through the butt (
Figure 1a), which gives higher values than taking MOI about the lateral in-plane axis passing through the CoM. The parallel axis theorem can be applied to calculate MOI about the lateral in-plane axis at any location along the length of the frame, using the values provided in
Table S1 in the Supplementary Material (e.g., for obtaining ‘swingweight’ as defined
Figure 1a). Polar MOI, or ‘twistweight’, is the MOI acting about the longitudinal axis of the racket (
Figure 1a). Transverse and Polar moments of inertia were estimated using models from [
29], with a Polar MOI model that estimates the racket head as a circle selected for its simplicity. The Polar MOI model could not be used for the 23 (4%) lopsided (or asymmetric) frames. Measured MOIs for the 416 (79.5%) rackets from previous work [
28,
29] were plotted against date alongside the estimated values, to check the models gave the correct trends over time.
The frequency and damping ratio (
) of the first mode of each racket were obtained from modal analysis. The racket was suspended by a long string with a single axis accelerometer (TE Connectivity Model 805) strapped to the handle (antinode), and the string-bed was struck lightly on the longitudinal axis away from the node with a ball to excite vibrations. The accelerometer was connected via a signal conditioner (PCB Model 480E09) to a digital oscilloscope (PicoScope 2000 Series), and then to a computer where vibrations could be viewed live and natural frequencies (i.e., fundamental mode of the freely suspended racket) obtained using the oscilloscope software (Picosope 6). The accelerometer sampled at 3 kHz and recorded for 5 s with the trigger located at 10% of the time window, giving a frequency resolution of 0.2 Hz. The vibration signals for the 91 (17%) rackets characterised previously at the brand were obtained using a wireless accelerometer (MetaWear CPro, Mbientlab Inc) sampling at 800 Hz, as described in [
28]. A check using a racket model (Ti.S6, HEAD) located within both the brand and museum confirmed that the two accelerometers gave the same result for natural frequency (176 Hz).
The recorded signals for all rackets were post-processed in MATLAB
® to identify the natural frequency and damping ratio. The script performed a Fourier transform of the signal with the MATLAB
® function
fft, and the first peak in the power spectrum (f
p) was assigned to the natural frequency of the racket. Similar to [
28], the wired accelerometer was calibrated using a shaker excited at frequencies of 100, 150 and 200 Hz. Each calibration test was repeated three times, with the frequency from the accelerometer always matching the excitation value. Damping ratio was calculated with the half-power bandwidth method (
Figure 2), as commonly used for estimating damping in multi-DOF systems [
32]. The method first involves identifying the two half-power frequencies (f
1 and f
2), which are located either side of the peak frequency were the amplitude (A
1) is equal to the amplitude at the peak (A
p) divided by √2 (i.e.,
). The half-power damping ratio can then be obtained using Equation (3),
where
is the frequency range. The number of data points in the region of the power spectrum that corresponds to the natural frequency (± 10 Hz of f
p) was increased using a
spline function in MATLAB
®, to improve the frequency resolution (0.02 Hz) and the accuracy of the damping estimation, as shown in
Figure 2. When applying the half-power bandwidth method, for a given natural frequency, damping ratio will increase with the frequency range (wider and less pronounced peak), while for a given frequency range, damping ratio will increase as natural frequency decreases. Each racket was tested twice, with the mean natural frequency (to 1 Hz) and damping ratio reported. The uncertainty in the calculation of damping ratio was dependent on both the damping ratio and natural frequency of the racket, as outlined in Equation (4),
For example, the uncertainty in damping ratio for a wooden racket with a natural frequency of 100 Hz and damping ratio 1% would be 2.2% of the calculated value. Four hundred and ninety-four (94%) rackets had an uncertainty in damping ratio below 5%, 28 (4%) rackets had an uncertainty between 5 and 10%, and 3 (<1%) rackets had an uncertainty between 10 and 15%. The highest uncertainty in damping ratio was associated with metal rackets, because they have both low frequency and low damping.
The resistance of a racket to bending can be increased using stiffer materials or larger cross sections. Fibre–polymer composites allow for rackets with large hollow cross sections, combining a high second moment of area with low mass. Measuring the cross section of hollow rackets was not possible without cutting them or taking them to a three-dimensional imaging device, such as a computed tomography scanner, so the ratio of the frame depth to thickness was used as a proxy for the second moment of area. The natural frequency of a freely suspended racket is often used as an analogue to its stiffness [
33], although it is also affected by mass [
34].
To reduce the dimensionality of the number of measured variables, a principal component analysis (PCA) [
35] was conducted for all racket frame measurements (racket, head (int.) and grip length, racket width (int.), frame thickness and depth, frequency and damping, mass, CoM location, and transverse and polar MOI). The three principal components that accounted for the most variance in the measured variables were then further analysed. Firstly, they were examined to see if they varied between material or brand (only including the 12 brands that were represented by at least six rackets, equating to 271 (52%) samples). Material and brand were introduced as between-factors in a multivariate ANOVA, with PCA1, PCA2 and PCA3 as the dependent variables. Partial eta squared (η
p2) was also used to quantify the effect sizes of these analyses, where η
p2 > 0.01 is a small effect, η
p2 > 0.06 is a medium effect and η
p2 > 0.14 is a large effect [
36]. A bivariate Pearson’s correlation was conducted on PCA1 and PCA2 to see which of the measured variables were best correlated to them, in terms of their correlation coefficients, and the variables where r > 0.7 are presented.
In most instances, the measured parameters were plotted against date with a data point for each racket colour coded for the material categorisation. A moving average of the data was calculated in MATLAB®, using the smooth function (moving average filter) and setting a method (rloess) that assigns lower weight to any outliers, and plotted as a trend line. A moving standard deviation over 50 consecutive points was calculated using the movstd MATLAB® function, and was shown in the background of the graphs shaded green.
4. Discussion
Collected racket parameters were approximated by the main principal components (PCA1). Head sizes and natural frequencies of tennis rackets have changed most significantly during their development and correlate best with PCA1. We suggest that these main changes were likely to be facilitated by the increased stiffness of composite materials. Wood was the common material for tennis rackets until the 1960s, when composites were used as reinforcement on a wooden frame. By the 1970s, engineers were experimenting with steel, aluminium and fibre–polymer composites, as they are stronger than wood and allowed for rackets with larger heads. Oversized rackets made tennis easier to play, but metal frames had lower natural frequencies than their composite counterparts and offered limited damping, so that vibrations from ball contact away from the node (the ‘sweet spot’) lasted longer [
40]. By the 1980s, composites had surpassed wood as the common racket material, providing engineers with more freedom over frame shape, stiffness, and mass distribution.
Many racket properties, such as the frame cross section shape (
Figure 6b), mass (
Figure 7a), Polar MOI (
Figure 7c) and natural frequency (
Figure 8a) were more varied in modern designs, giving the player greater choice to suit their style, experience and preference. The analysis presented here did not reveal a lasting brand that departed radically from established principles of tennis racket design. All brands tended to follow the same fundamentals with regards to head size, stiffness and material. Differences revealed in the analysis appear to have been due to the longevity of particular brands and the era of their presence within the marketplace. Observed differences were due to restrictions in racket material for a brand prevalent in the wooden era compared to modern tennis racket brands.
This study provides the first presentation of how string bed densities, damping and external cross section dimensions of tennis rackets have changed over time. While other racket parameters have been previously documented by Haake et al. [
3], over three times more rackets were characterised here with a wider variety of measurements, analysis methods and reporting of materials. Findings agreed with Haake et al. [
3] in terms of trends in properties since the 1870s, including overall length, grip and head length, head width, mass, CoM location, natural frequency and transverse MOI (reported as swingweight (
Figure 1a) in [
3]), but this study also shows how materials have driven the development of the tennis racket. Haake et al. [
3] only includes rackets up to 2007; data presented here indicates that the rate of change for some properties, such as head size (
Figure 5a) and mass (
Figure 6a), may have slowed recently, since 2007. Haake et al. [
3] showed scatter in their results for Polar MOI with no clear trend over time, while those reported here show little change until 1970 followed by a steady increase with the move to rackets with wider heads, as expected [
29,
38]. Haake et al. [
3] estimated Polar MOI from measurements of transverse and lateral MOI (using a Babolat Racquet Diagnostic Centre) [
3], and the assumption that a racket is planar, combined with uncertainty from two measures may have influenced their results [
28]. Holding a racket by the handle can change modal frequencies and shapes and substantially increase damping [
40,
41,
42], and future work could investigate the dynamic properties of rackets from different eras in more depth, if freely suspended, hand-held and with mass added at the handle to represent the hand.
While over 500 rackets were characterized, they do not represent all of the rackets that have ever been made. The reason each racket entered the museum collection is unknown, and could be due to factors such as availability, survival, design, and significance. The museum store is climate controlled, but the extent to which time, use or manufacturing inconsistencies affected racket properties, like natural frequency and damping, is unknown. The measurement techniques are time-consuming when applied to many rackets, and future work could look to make them quicker. Image processing could be used to obtain dimensions automatically from photographs [
29], removing the need for manual measurements while improving estimates for properties like frame thickness and head area, and facilitating measurement of string spacing. Automated three-dimensional imaging devices [
43] could also be used to scan rackets, to give a better representation of shape, including wall thickness in hollow metal and composite frames, which would allow for calculation of second moments of area. More efficient protocols could facilitate testing of more rackets, to produce a better picture of how they have developed over time. Repeat testing of different examples of the same racket could help to account for any manufacturing inconsistencies and deterioration.
The data collected here shows a distinct transition in racket design between the 1960s and 1990s, with almost all rackets now being made from fibre–polymer composites (
Figure 3). It is possible that some parameters may have reached a new equilibrium today with, for instance, racket masses of around ~330 g, natural frequencies of ~143 Hz and transverse MOIs of ~0.047 kgm
2. It might be beneficial to investigate how constituent materials and fiber orientations of composite rackets could influence performance, such as ball velocity and spin, and player perceptions and preferences [
44,
45], as well as injury risk [
9,
22,
23]. While fibre–polymer composites have improved the tennis racket, their damping is limited, and manufacturing is currently labour intensive. Finite element models [
13] could be produced for rackets from different eras, to further our understanding of how their development has influenced player performance, while also serving as design tools. Such finite element models could facilitate the application of new materials and manufacturing techniques to the tennis racket. Natural fibre composites, as used in some surfboards and fishing rods, could reduce the cost and environmental impact of making rackets [
46], while auxetic composites [
47] and those additively manufactured with continuous fibres [
48,
49,
50] could bring more design freedom and support mass customisation for different players. Sensors incorporated within the frame could help brands develop a better understanding of how players use their rackets, and support mass customisation strategies. Simple measures, such as those described here, can be used to monitor changes in racket design.