Participants
Data were collected from fifteen healthy male baseball pitchers. Mean age was 24.5 years (SD 7.5, range 17-44), body height 191 cm (SD 5, range 183-199), and body mass 79.4 kg (SD 9.2, range 62.7-102.5). Of the fifteen tested pitchers, 11 were right-handed. Most participants were pitching at a recreational level, with two participants playing at the highest level in the Netherlands. None of the participants had experienced any musculoskeletal injuries in the past six months nor had they received elbow surgery in the past. The study protocol followed the guidelines stated in the Declaration of Helsinki 13 and was approved by the Ethics Committee of the Delft University of Technology (HREC). Participants were informed of the procedure before the start of the measurements. Informed consent was obtained before involvement in the study.
Procedure
The measurements were performed at the indoor human movement laboratory of the Department of Human Movement Sciences at the Vrije Universiteit Amsterdam, The Netherlands. Fourteen reflective markers were placed on anatomical bony landmarks of the participants with double-sided tape (see table 1). Electromyography (EMG) electrodes were placed on the skin of the throwing arm and an accelerometer was attached to the sternum below the incisura jugularis. The participants wore their own shoes, athletic shorts, and baseball glove, but no shirt. Prior to performing fastball pitches, participants had to perform maximum voluntary contractions (MVC, see table 1 supplementary materials). Participants gradually built-up muscle force and held this for 3 seconds. Each MVC was repeated three times. After performing their regular warm-up, the participants were instructed to pitch fastballs at full effort. Ten fastball pitches were performed within a block of pitches, with two minutes rest between each block. Before the start and between 10 blocks of pitches, the participants were asked about their self-perceived fatigue with the following text: “Place a vertical line on the visual analog scale shown below in which way you are overall fatigued”. The visual analog scale (VAS) was ranging from totally not fatigued (0%) to extremely fatigued as possible (100%). Participants were instructed to stop when having thrown 110 fastballs or when their VAS score reached 80%. The minimum required number of pitches was 60. To investigate the effect of fatigue on variability, all the pitches were measured and included in the analyses. Participants pitched from a pitching mound (height 0.55m) towards a strike zone (height 0.71m; width 0.43m), at 18.66 m.
Data acquisition
Kinematics and ball speed
Marker positions were recorded using an OptiTrack motion capture system with twelve cameras sampling at 120 Hz (OptiTrack Flex 13, OptiTrack™, Corvalis, United States). The OptiTrack system was calibrated to define camera position and orientation and to construct a convenient global coordinate system. The ball speed was measured behind the strike zone using a stalker pro radar gun (Stalker Radar, Plano, TX, USA).
Electromyography
Muscle activity of three elbow skeletal muscles of the throwing arm was measured using bipolar surface electromyography (sEMG). The flexor pronator mass (FPM), biceps brachii (BIC), and triceps brachii (TRI) muscles were measured (table 1). The electrode locations were based on the SENIAM guidelines (Hermens et al. 1999). The reference electrode was placed on the clavicle of the non-throwing arm. Disposable bipolar electrodes (Ag-AgCl; 1 cm2 recording area; Blue Sensor Electrodes N-00-S, Ambu Inc., USA) were attached in the direction of the muscle fibers with 2 cm distance between the centers of the electrodes. Before the electrodes were attached, the skin was shaved and cleaned using alcohol. The electrode cables were fixated to the skin to avoid cable movement artifacts in the signal and to minimize the risk of loosening of the electrodes from the skin during pitch movement. The cables were connected to a BioPlux research device (Plux biosignals, Lisboa, Portugal), with 16-bits analog channels, a gain of 506, and an analog 25-500Hz band-pass filter. All consecutive fastball pitches of a participant were recorded in one EMG dataset at a sampling frequency of 2000Hz and locally stored on the BioPlux research device.
Data analysis
All data analyses were performed in Python (version 3.7, Python Software Foundation, https://www.python.org/).
Kinematics and inverse dynamics
The following bony landmarks on the throwing arm were used to construct an anatomical local coordinate system for the hand, forearm, and upper arm according to the ISB recommendations 14: third proximal interphalangeal, ulna processes styloid, radius processes styloid, lateral humeral epicondyle, medial humeral epicondyle, and the acromion. Positions of the centers of mass and the moments of inertia were estimated according to Zatsiorsky (2002)15 and De Leva et al. (1996)16. The elbow joint angles were decomposed in the rotation order of ‘flexion/extension’ – ‘ab/adduction’ (floating angle)- ‘pronation-supination’ according to Grood and Sunday (1983). Maximal external shoulder rotation (MER) was obtained from the shoulder joint angles decomposed according to the y-x-y Euler decomposition (‘plane of elevation’-‘negative elevation’-‘axial rotation’) 17.
The net joint forces and moments were calculated in the global coordinate system, using a top-down inverse dynamics analysis based on the Newton-Euler equation of motions. Subsequently, the elbow joint torque was expressed in the anatomical coordinate system of the elbow; positioned in the middle of the medial and lateral humeral epicondyles. The kinetics of the segments were calculated with the segment data and scaling factors of De Leva et al. (1996) and Zatsiorsky et al. (1990). A 2nd order polynomial function was fitted using three measured data points to obtain the exact magnitude of the peak value of the external valgus torque, which occurred around the moment of MER. The inverse dynamical model can be found here: https://github.com/ThomasBTHL/BTHL_public.
Electromyography
EMG signals were first separated into the ten-pitch series. Subsequently, these were cut into single pitches. The linear envelope was obtained by rectifying the EMG and applying a fourth-order bi-directional lowpass Butterworth filter of 20Hz. EMG data were normalized to the maximum values observed in the MVC data. To quantify the indirect effect of the biceps and triceps muscles, a co-contraction index (CCI) was calculated for the biceps and triceps muscle pair at each sample () according to Rudolph et al. (2000)18, see equation (1).
An area under the curve (AUC) was calculated over a window of 150ms for the normalized EMG data and the CCI. This window started at MER at 0ms and ranged back to -150ms. The 150ms was chosen because foot contact to ball release is ~150ms19. Therefore, when considering a 50ms electromyography delay 20, the AUC can be seen as indicative for the produced force of the elbow muscles from foot contact to ball release, in which the moment of the peak valgus torque occurs.
Synchronization
The BioPlux device, containing the EMG signals and accelerometer data, did not contain the MER event. Therefore, it was synchronized with the OptiTrack system. The z-direction of the accelerometer, pointing forwards relative to the thorax, was synchronized with the forward acceleration of the trunk coordinate system. For each pitch, the data were synced on the peak linear accelerations and stored in a Python pickle.
Moving window approach
In addition to ball speed, four outcome variables were analyzed in relation to repetitive pitching: the elbow valgus torque magnitude, valgus torque variability, the FPM AUC, and the biceps-triceps CCI. A moving window of ten pitches was applied to all variables and moved over a single subsequent pitch. The mean of the ten values within each window was used to quantify the course of ball speed and the valgus torque magnitude, and the standard deviation of the ten values within each window was used to quantify the course of the within-individual valgus torque variability over the individual sessions of 60-110 pitches (see figure 1). For the EMG outcome variables, the means of the ten values of the moving 10-pitch windows were quantified as the FPM activity and the biceps-triceps CCI.
Statistical analysis
To statistically explore the relationship between ball speed and the four outcome variables and repetitive pitching five linear mixed models (LMM) were examined. The LMM deals with missing data and data of different samples 21, which is an advantage as some participants indicated to be fatigued >80% after already 60 balls, whereas others did not reach this level of fatigue after 110 balls. The fixed factor was the number of subsequent moving windows for each individual series of pitches and the random factor was the participant. Given the multilevel structure of the data (level 1 pitch window number nested in level 2 participants), it was considered necessary to build three models: (1) a basic model with a random intercept across participants, (2) a model with pitch window number as predictor and random intercept across participants and (3) a model with pitch window number as predictor, a random effect of pitch window number over participants (random slope) and random intercepts. To select the best-fitted model, the models were compared using a chi-square likelihood ratio test with a significant level of 0.05. If the models were significantly different, the model with the smallest AIC value was used. The maximum likelihood was used as the estimation method. The nlme package for R was used to perform the LMM analysis 22. All statistical analyses were performed in R (version 4.2.0) 23 and Rstudio (version 2022.2.0.443) 24.