Introduction

The Russian bar, a captivating circus discipline born in the mid-20th century, showcases the synchronized performance of three artists: two porters and one flyer. In a classical Russian bar act, the porters support a flexible composite bar on their shoulders while the flyer performs acrobatic maneuvers by jumping on the bar. The flyer, akin to a trampoline athlete, executes aerial maneuvers such as candle jumps and saltos, described as straight jumps with backward or forward rotations. The porters accumulate energy by flexing their hips and lower extremity joints, and then extending them to propel the flyer during the propulsion phase. During the flight, the porters ensure the flyer’s safety by keeping the bar below the flyer’s hips and feet. In the landing phase, the porters absorb the flyer’s impact on the bar by flexing their hips, ankles and knees to optimize energy absorption.

One of the major concerns about the practice of Russian bar is the potential development of chronic low back pain among the porters.1–4 This might result from high compressive forces, or asymmetrical forces on the porters’ shoulders and back. Prior research has indicated that asymmetric lifting carries a higher risk of developing low back pain compared to symmetric lifting.5–7 Circus companies are committed to addressing artists’ injuries and implementing preventive measures. However, research investigating the spinal loads imposed on artists’ bodies is lacking. Conducting comprehensive studies to understand the dynamic loading conditions experienced by porters is essential to defining effective prevention strategies for Russian bar artists.

Human musculoskeletal models serve as valuable tools for assessing joint forces applied to the body8 using motion capture and by acquiring external forces.9 These models have been used in sport activities including dance,10 triple-jump11 and golf.12 Many biomechanical simulations have been created to assess the load on the lumbar spine during a variety of lifting or throwing activities.5,13,14 However, no biomechanical model has been specifically tailored for Russian bar circus artists. The present study aims to quantify the load exerted on the lumbar spine of Russian bar porters using a musculoskeletal model. We hypothesized the following: 1) Russian bar generates asymmetric loads on the lumbar spine, 2) saltos generate higher loads on porter’s lumbar spine than candle and consecutive jumps and 3) landing phases generate higher lumbar spine loads than propulsion and stationary phases. By undertaking this investigation, our goal is to enhance our understanding of the forces experienced by Russian bar porters and contribute to the development of injury prevention strategies for circus companies.

Methods

Experimental data collection

Ethics and participants

Written consent was obtained from three healthy Russian bar artists (two male porters, designated as porter 1 and 2, aged 31 and 41 years old, respectively, and a 34-year-old female flyer) following approval of the ethics committees of École nationale du cirque (CER 2122–15C), École de technologie supérieure (H20220512) and Cégep de Lanaudiere (2022-09-08-01). Mass was measured for each participant (porter 1 = 74.4 kg, porter 2 = 81.7 kg, flyer = 57.8 kg). All were professional experienced circus artists having practiced Russian bar together for five years.

Experimental protocol

Video recordings were conducted during all activities. An eighteen-camera motion analysis system (OptiTrack, Natural Point, Oregon, US) tracked the movements of the porters, the position of the flyer and the deflection of the bar in 3D space. Ninety-two reflective markers were placed on one porter’s body (including twelve rigid marker sets called rigid bodies) following the instructions of Beaucage-Gauvreau,9 while eighteen markers were attached along the length of the Russian bar (Figure 1b). As noted in Figure 1a, due to limited camera availability, only one porter could be recorded at a time. Additionally, six markers were placed on the flyer’s hips and ankles to follow the flyer’s height and position. Eight electromyographic (EMG) sensors (Delsys Inc., Massachusetts, US) were used to record muscle activity on specific muscles: rectus abdominis (RA), external oblique (EO), thoracic erector spinae (TES) and lumbar erector spinae (LES).

The figure shows four parts: 1) Three artists performing Russian bar. Two porters have the bar on one shoulder and the flyer prepares to jump on the bar. 2) An artist doing a T-pose. He is wearing OptiTrak sensors on the body. 3) Four illustrations showing a typical sequence of a candle jump. 4) Four illustrations showing a typical sequence of a salto jump.
Figure 1.

a) Experimental setup with the Russian bar trio. b) Porter equipped with markers and rigid bodies on the head, arms, forearms, thighs and shins. c) Illustration of a candle jump. d) Illustration of a salto jump.

Experimental data were collected for three distinct jump activities performed by the flyer:

  1. candle jump: a straight jump (Figure 1c);

  2. salto jump: a rotation in the sagittal plane (Figure 1d);

  3. consecutive jumps: a series of successive candle jumps.

Figure 2 illustrates typical phases of these three activities according to the bar deflection. The graphs are derived from the deflection measurement taken during the data collection. All activities start with an oscillation phase characterized by small-amplitude up and down motions where the flyer’s feet remain in contact with the bar. Once sufficient speed and balance are reached, the flyer initiates the jump from the bar during the propulsion phase. Subsequently, the flyer lands on the bar during the landing phase. It’s important to note that a salto jump is always preceded by a single candle jump (Figure 1d), and that consecutive jumps consist of successive candle jumps of consistent heights. These consecutive candle jumps form a so-called stationary phase, where there is not a proper propulsion phase as the height of the flyer is gradually increased using the bar’s flexibility almost exclusively. The landing phase is identical to that of a candle jump and, therefore, will not be analyzed in this study.

Three graphs show the deflection of the Russian bar during a candle jump, salto jump and consecutive jumps, respectively. Jump phases are shown in the different graphs. The candle jump and salto jump exhibit three phases as oscillation, propulsion and reception phase. Propulsion and reception phases are defined as peak in the Russian bar deflection. The consecutive jump graphs show two phases, the oscillation and stationary phase expressed by consecutive peaks of the same magnitude.
Figure 2.

Bar deflection over time for three jump activities: a) candle jump, b) salto jump and c) consecutive jumps. Pictograms visually depict the activity at various phases and time points, with shaded areas representing the phases analysed in this study.

Two force plates (Bertec portable model 6090–06, Bertec Corporation, Ohio, US) measured the ground reaction forces (GRFs) generated by each of the porter’s feet during candle jumps lower than two meters in height. The force plates were removed for larger jumps as the porters constantly adjust their position to ensure the flyer’s safety. The force plate data were ultimately used for model validation.

Methods for determination of shoulder forces

Due to challenges in directly measuring external forces on the porters’ shoulders, a two-step indirect approach was employed. The first step involved characterizing the bar mechanical behavior under a static load (Figure 3a), while the second step examined the influence of dynamic loading on its mechanical response (Figure 3b).

Three parts: 1) Experimental setup with a Russian bar bending above two rigid wooden supports. Five arrows explain the different objects shown on the picture which are: OptiTrack camera, markers, weight, rigid support and load cell. 2) Two pictures show a flyer who jumps on a bar during a dynamic test. The bar rests on two small wooden support. Two individuals ensure the extremity of bar doesn’t fall off the support while the flyer jumps. c) Graph showing linear relation between force and deflections.
Figure 3.

a) Experimental setup for estimating the static force applied on the porters’ shoulder as a function of bar deflection. b) Experimental setup and sequence for estimating the dynamic force applied on the porters’ shoulder as a function of bar deflection. c) Graphic depicting the linear relation between the shoulder’s force and bar deflection in both static and dynamic conditions. This graphic also illustrates that the bar exhibits slight flexion without any applied force due to the weight of the bar and loss of structural integrity from use over time (0.03 m deflection without any force applied).

A linear relationship between static bar deflection and shoulder force was observed (Figure 4c), leading to F = K∂ + C (Equation 1) where F (N) denotes the shoulder force, K (N/m) denotes the bending stiffness of the bar, ∂ (m) represents the bar deflection and C (N) is a constant representing the mathematical error of the linear regression. (This value is negligible, as it is below 2% of the maximum force in each trial.) Using this experimental protocol, a bending stiffness of 1,703.3 N/m was determined (Figure 3c).

Jump activities can also be accurately described by observing the force applied to the porter’s shoulder, as depicted in Figure 4. Jump phases (i.e., propulsion, landing and stationary phases) were further divided into loading and unloading steps. Between unloading and loading steps, the porters are in a standby phase.

A line graph shows the force applied (in Newtons) over several cycles of an acrobatic movement. The Y-axis represents force in Newtons (N), ranging from 0 to 2000 along the jump cycle. The graph is divided into three phases representing the oscillation and standby phases, the loading phase where force increases and the unloading phase where force decrease. Above the graph, illustrations depict an acrobat performing a sequence of movements.
Figure 4.

Shoulder force during a typical candle jump activity showing loading (red line) and unloading steps (green line) for a propulsion and a landing phase.

Musculoskeletal modeling

OpenSim (version 4.4)15 was used to assess the lumbar spine loading of the Russian bar porters. OpenSim provides sophisticated computational resources and is the most-used open source platform in the biomechanics field. Accordingly, it was selected for the current research study. The full-body model of the porter is an adapted version of a comprehensive lifting model originally developed and validated for analyzing spinal loads during lifting tasks.13 The modifications made to the original model from Beaucage-Gauvreau are presented in the Appendix. The use of a full-body model was mandatory for our study as the porters interact with external objects and use a wide range of motion. This model was validated for lifting tasks involving forward bending and is simpler than other full-body models, thus reducing computational time. The trunk exhibits three degrees-of-freedom: flexion/extension, lateral bending and axial rotation. The trunk’s musculature is represented by eight major muscle groups: erector spinae (ES), rectus abdominis (RA), external obliques (EO), internal obliques (IO), multifidus (MF), quadratus lumborum (QL), psoas major (PS) and latissimus dorsi (LD). In a normal anatomical position, the angle convention is zero degrees and increases with joint flexion.

External forces

Three external forces were applied to the model: the force exerted by the bar on the shoulder and the two ground reaction forces. Shoulder forces were applied to the torso on a point located fifteen centimeters medial to the right acromion in the frontal plane. To ensure consistency and simplicity, the force was applied in the sagittal plane, perpendicular to the bar. The direction was obtained by calculating the normal vector of the bar at the location of the shoulder (Figure 5a and 5b). Focusing on the sagittal plane captured primary biomechanical effects while minimizing complexity. Force intensity was computed using the bar deflexion measured during jump activities and the stiffness coefficient of 1,703.3 N/m.

Three parts: 1) Isometric view of the OpenSim model showing an arrow starting from his right shoulder and pointing downward. 2) Lateral view of the OpenSim model showing an arrow starting from his right shoulder and pointing downward. 3) OpenSim model with the lumbar local body axis such that the X-axis points forward, the Z-axis points to the right and the Y-axis points upward.
Figure 5.

a) Isometric and b) lateral views of the model during a loading phase. The pink arrow represents the load applied by the bar on the right shoulder of the porter. c) Lumbar column of the musculoskeletal model (L1 in blue, L2 in yellow, L3 in red, L4 in green and L5 in purple) and the L4 local axis system.

Computation of ground reaction forces with the model was accomplished through the implementation of inverse dynamics-based simulations using an actuator method, as described in previous studies.16,17 This computational approach enabled the accurate estimation and quantification of the forces exerted on the ground during Russian bar, providing valuable insights into the interaction between the performers and their surrounding environment. This indirect method of estimating GRFs was validated using measurements acquired from the two force plates during a static pose and candle jumps of less than two meters in height.

Data processing

Data processing was done with MATLAB (2022a) and OpenSim. Statistical analysis was performed using Minitab (21.1). Kinematic and kinetic data were low-pass filtered using a fourth-order, zero-lag Butterworth filter with a cut-off frequency 6 Hz. EMG signals were band-pass filtered (10–400 Hz), rectified and filtered at 5 Hz, normalized to the peak activation (measured during performed maximal isometric contractions) and then normalized to the peak activation of each trial. Muscular activation of the model was filtered and normalized to maximum peak activation.

Normality of data was assessed using the Shapiro-Wilk test and showed a non-normal distribution. Statistical analysis was performed using non-parametric paired and unpaired Kruskal-Wallis tests for comparing L4-L5 joint reaction forces obtained from the different jump activities (candle, salto and consecutive jumps) and jump phases (propulsion, landing and stationary phases). Statistical tests were conducted with 0.05 significance value. L4-L5 was selected for the analysis due to its recognized role in low back pain and for comparison with existing literature.18

Results

Overall, 47 jumps including seven saltos, eleven candle jumps and 29 consecutive jumps were recorded. Propulsion and reception phases were extracted from all recordings of candle jumps and saltos (EMG, kinematics and joint reaction forces), while stationary phases were extracted from consecutive jumps. For the kinematic analysis, propulsion and reception phases were further subdivided into loading and unloading steps. Compressive force is the force along the Y-axis expressed in local body frame. The shearing force is defined as the Euclidean norm of the shear force along the X-axis and the shear force along the Z-axis.

Jump analysis

Kinematic analysis

The propulsion and landing phases displayed distinctive articular kinematic patterns, as illustrated in Figure 6, indicating that both porters used different strategies for propelling and receiving the flyer. During the propulsion phases, the porters exhibited hip flexion during the loading step and hip extension during the unloading step. Conversely, in the landing phases, hip flexion steadily increased throughout the cycle. Propulsion phases showed left and right hip flexions angles of 81.4° ± 4.2° and 75.2° ± 14.3°, respectively. Stationary phases mirrored propulsion phases though with reduced motion at the left (43.0° ± 18.9°) and right (40.9° ± 16.1°) hips. This suggests that the porter predominantly leveraged the bar’s flexibility for upward propulsion during consecutive jumps. In the landing phase, the range of hip flexion was slightly higher than in the propulsion phase (left: 93.3° ± 4.9°; right: 91.1° ± 4.9°). Small changes in lumbar lordosis angle were recorded (<10°) as the porters kept the lumbar column “straight”.

The figure consists of three sections, each representing a phase of movement (propulsion phase, stationary phase and landing phase). Propulsion phase: Both hips show a steady increase in flexion angle from 0° to 90° during the loading phase (red lines) and decrease back to 0° during the unloading phase (green lines). Standard deviation is shown with dashed lines. Stationary phase: Hip flexion angles increase from 0° to around 50° during loading and decrease back to 0° during unloading for both hips, with a smaller range of motion compared to the propulsion phase. Landing phase: Flexion angles rise gradually from 0° to 60° during loading and further increase during unloading for both hips, indicating a different movement pattern for the landing phase with continuous flexion.
Figure 6.

Hip joint angle during the three phases of jumps for the two porters: propulsion (eighteen jumps); stationary (29 jumps); landing (eighteen jumps).

Joint reaction forces

Table 1 presents joint reaction forces at L4-L5, expressed in the local coordinate system of the L4 vertebra (Figure 5). No significant differences were found in compressive forces between candle and salto jumps during propulsion (p=0.143) or landing (p=0.364) phases. However, paired Kruskal-Wallis test applied to both candle jumps and saltos indicate that the propulsion phases generated significantly higher forces than the landing phase (p<0.001). Additionally, the stationary phases exhibited significantly lower forces than the propulsion (p<0.001) and landing phases (p=0.038).

Table 1.

Mean value of the peak internal forces and moments in the L4-L5 joint and hip flexion across all jump phases.

Jump

Phase

Hip flexion max (deg)

Joint reaction forces

Joint reaction moments

Compression (N)

Shear (N)

Mx (Nm)

My (Nm)

Mz (Nm)

Candle (11 jumps)

Propulsion

90.8

10,775.1

654.5

75.9

38.7

76.8

SD

8.7

2,530.9

275.2

38.9

22.4

26.3

Landing

90.8

7,807.4

557.9

38.1

11.5

63.3

SD

8.7

1,142.5

84.2

38.1

31.1

12.0

Salto (7 jumps)

Propulsion

91.6

12,571.7

961.6

121.5

60.2

89.8

SD

7.9

1,714.5

419.5

55.9

30.0

29.8

Landing

92.7

7,309.7

696.2

40.7

7.4

69.0

SD

7.9

645.1

86.7

11.4

23.9

5.3

Consecutive (29 jumps)

Stationary

37.5

6,353.6

434.9

22.8

2.4

42.9

SD

21.3

1,414.1

51.2

10.8

15.7

8.4

Salto jumps during propulsion induced significantly higher lumbar shear forces than candle jumps (p=0.04). Salto jumps showed higher forces during the propulsion than the reception phases (p=0.004). A paired test showed that there were no significant differences in shear forces between the propulsion and landing phases (p=0.163). Shear forces for both propulsion and landing phases were significantly lower during stationary jumps than in candle jumps (p<0.015).

Across all phases, the bending moment (Mz) in the sagittal plane exhibited the highest magnitude, closely followed by the moment in the frontal plane (Mx). The moment Mz was significantly higher in the propulsion phases than in other phases (p=0.029 vs landing, p<0.001 vs stationary). The moment in the frontal plane Mx was significantly higher in the propulsion phases than in other phases (p=0.002 vs landing, p=0.014 vs stationary). Mx was also significantly higher during the propulsion phases of salto jumps than during the propulsion phases of candle jumps (p=0.04). No other significant differences were found for Mx and Mz.

Discussion

This study utilized a complete musculoskeletal human model to estimate the load applied on the lumbar spine of Russian bar porters. As anticipated, the largest moments were exhibited in the sagittal plane (bending moments Mz) for all jump activities, reflecting the need for the porter to flex their hips to control the bar’s reaction. Consistent with our first hypothesis, the moment in the frontal plane (lateral bending) was similar in magnitude to the bending moment, suggesting an asymmetric load on the porter’s back. Our second hypothesis was not supported. Salto jumps exhibited significantly higher compressive forces and moments than successive jumps. These jumps also showed higher shear forces than candle and successive jumps. However, there were no significant differences in lumbar compressive forces and moments between salto and candle jumps. Our third hypothesis could not be confirmed with our results, as the propulsion phase did not generate greater forces than landing phases. Ultimately, jump height and hip flexion significantly contributed to the L4-L5 compressive forces, particularly in the case of successive jumps (R2=0.74 and 0.73, respectively).

Our findings indicate that a Russian bar porter endures high compressive lumbar forces, reaching up to 12,571 ± 1,715 N, but relatively low moments, with a maximum value of 121.5 ± 55.9 Nm. This can be explained by the observation that the porter maintains their low back in a straight position (low lumbar spine range of motion < 10°) limiting the lever arm on the vertebrae and aligning their low back with the external force applied on their shoulders in the sagittal plane. In comparison, Eltoukhy et al. validated a model for assessing weightlifting exercises, showing values of 6,224 ± 1,753 N (snatch), 7,963 ± 2,784 N (deadlift) and 8,701 ± 3,263 N (clean), with bending moments in the sagittal plane of 734 ± 331 Nm (deadlift) to 1,731 ± 1,410 Nm (clean).19 Schmid et al. reported compressive values up to three times body weight for lifting tasks, whereas our study identified forces reaching up to ten times body weight.20 While our calculated forces may appear remarkably high, previous studies have seldom delved into such high efforts. Comparable loading conditions were found in strongman activities,21 where participants sustained compressive forces surpassing 12,000 N and moments exceeding 300 Nm during the super yoke walk and keg walk (with a 155 kg keg on the right shoulder), respectively.

The model effectively estimates lumbar spine loading in Russian bar exercises. However, prudence is advised when interpreting our quantitative results, given the reliance on assumptions and the constraints posed by our small sample size. Specifically, there are certain limitations in the estimation of external loads, as both of our external forces are indirectly estimated. While the bar characterization reliably estimated the vertically-applied loads in the sagittal plane, its accuracy diminished when assessing loads applied in the frontal plane. The interaction between the bar and the porter’s body is intricately complex, as the bar is “blocked” between the neck and the arm/forearm. This interaction is heavily influenced by the porter’s musculature and specific morphology. While applying the external force through a point near the acromion may not be the most precise representation, the agreement between calculated and experimentally measured ground reaction forces supports this simplification. Another interesting point that is not addressed in this study is the coordination between all three artists, which might influence lumbar biomechanical responses. A new study should specifically focus on the coordination between the flyer’s jump and the porters’ movements.

This model will assist coaches and artists in evaluating new bar prototypes and techniques. It helps quantify mechanical stress on artists, thereby aiding in injury prevention by identifying risk factors. Furthermore, in light of our findings, artists and coaches should consider various factors. Given the direct correlation between shoulder force and the mechanical characteristics of the bar, opting for a more flexible bar is advisable. This modification would primarily alleviate shoulder loading but might limit artistic possibilities. Additionally, a bar that is too flexible could pose a risk, potentially touching the ground during landing and endangering the flyer. Consequently, an optimal range of bar flexibility values must be found. As the forces involved in this discipline are exceedingly intense, special precautions are recommended, particularly for young and novice artists. Coaches should implement suitable exercises and training programs to mitigate exposure to high forces or ensure that artists are physically prepared to withstand these forces. Repetitive force on the spinal column can potentially lead to a decline in low back function.22 As illustrated by Kazemi,23 joint reaction forces increase with muscular fatigue, highlighting the need to limit training duration to prevent such effects. Enhanced co-contraction has the potential to reduce intervertebral forces.24

Conclusion

In conclusion, the proposed model has provided valuable insights into the loads exerted on the lumbar spine of Russian bar porters, making a substantial contribution to the understanding of the lumbar spine biomechanics associated with this captivating acrobatic discipline. Future research should focus on refining the modeling of external forces and muscle actions to improve accuracy. Exploring various loading scenarios, including more acrobatic maneuvers at different heights, would be intriguing and could further enhance our comprehension of the discipline while supporting the well-being of circus artists.

Acknowledgements

The authors thank all those who contributed to this paper. We would like to thank all the institutions that help the development of the project, including ETS and ENC for the technical help. We would also like to thank MITACS for the funding.

Funding

This work was supported by the MITACS Acceleration program under Grant IT31369.

References

1 Long AS, Ambegaonkar JP, Fahringer PM. Injury Reporting Rates and Injury Concealment Patterns Differ Between High-school Cirque Performers and Basketball Players. Med Probl Perform Art. 2011 Dec 1; 26(4):200–205. https://doi.org/10.21091/mppa.2011.4032https://doi.org/10.21091/mppa.2011.4032

2 Munro D. Injury Patterns and Rates Amongst Students at the National Institute of Circus Arts: An Observational Study. Med Probl Perform Art. 2014 Dec 1; 29(4): 235–240. https://doi.org/10.21091/mppa.2014.4046https://doi.org/10.21091/mppa.2014.4046

3 Stubbe JH, Richardson A, van Rijn RM. Prospective cohort study on injuries and health problems among circus arts students. BMJ Open Sport Exerc Med. BMJ Specialist Journals; 2018 Jun 1; 4(1):e000327. https://doi.org/10.1136/bmjsem-2017-000327https://doi.org/10.1136/bmjsem-2017-000327

4 Wolfenden HE, Angioi M. Musculoskeletal Injury Profile of Circus Artists: A Systematic Review of the Literature. Med Probl Perform Art. 2017 Mar; 32(1):51–59. https://doi.org/10.21091/mppa.2017.1008https://doi.org/10.21091/mppa.2017.1008

5 Kim HK, Zhang Y. Estimation of lumbar spinal loading and trunk muscle forces during asymmetric lifting tasks: application of whole-body musculoskeletal modelling in OpenSim. Ergonomics. 2017 Apr 3; 60(4):563–576. https://doi.org/10.1080/00140139.2016.1191679https://doi.org/10.1080/00140139.2016.1191679

6 Marras WS, Lavender SA, Leurgans SE, Fathallah FA, Ferguson SA, Allread WG, Rajulu SL. Biomechanical risk factors for occupationally related low back disorders. Ergonomics. 1995 Feb; 38(2):377–410. https://doi.org/10.1080/00140139508925111https://doi.org/10.1080/00140139508925111

7 Marras WS, Ferguson SA, Burr D, Davis KG, Gupta P. Spine loading in patients with low back pain during asymmetric lifting exertions. Spine J. 2004 Jan 2; 4(1):64–75. https://doi.org/10.1016/S1529-9430(03)00424-8https://doi.org/10.1016/S1529-9430(03)00424-8

8 Reinbolt JA, Seth A, Delp SL. Simulation of human movement: applications using OpenSim. Procedia IUTAM. 2011 Jan 1; 2:186–198. https://doi.org/10.1016/j.piutam.2011.04.019https://doi.org/10.1016/j.piutam.2011.04.019

9 Beaucage-Gauvreau E, Brandon SC, Robertson WS, Fraser R, Freeman BJ, Graham RB, Thewlis D, Jones CF. A braced arm-to-thigh (BATT) lifting technique reduces lumbar spine loads in healthy and low back pain participants. J Biomech. 2020 Feb 13; 100:109584. https://doi.org/10.1016/j.jbiomech.2019.109584https://doi.org/10.1016/j.jbiomech.2019.109584

10 Alderson J, Hopper L, Elliott B, Ackland T. Risk Factors for Lower Back Injury in Male Dancers Performing Ballet Lifts. J Dance Med Sci. 2009 Sep 1; 13:83–9. https://doi.org/10.1177/1089313X0901300303https://doi.org/10.1177/1089313X0901300303

11 Allen SJ, King MA, Yeadon MR (Fred). Optimisation of phase ratio in the triple jump using computer simulation. Hum Mov Sci. 2016 Apr 1; 46:167–176. https://doi.org/10.1016/j.humov.2015.12.012https://doi.org/10.1016/j.humov.2015.12.012

12 Mahadas S, Mahadas K, Hung GK. Biomechanics of the golf swing using OpenSim. Comput Biol Med. 2019 Feb 1; 105:39–45. https://doi.org/10.1016/j.compbiomed.2018.12.002https://doi.org/10.1016/j.compbiomed.2018.12.002

13 Beaucage-Gauvreau E, Robertson WS, Brandon SC, Fraser R, Freeman BJ, Graham RB, Thewlis D, Jones CF. Validation of an OpenSim full-body model with detailed lumbar spine for estimating lower lumbar spine loads during symmetric and asymmetric lifting tasks. Comput Methods Biomech Biomed Engin. 2019 Apr 4; 22(5): 451–464. https://doi.org/10.1080/10255842.2018.1564819https://doi.org/10.1080/10255842.2018.1564819

14 Zaman R, Xiang Y, Cruz J, Yang J. Three-dimensional asymmetric maximum weight lifting prediction considering dynamic joint strength. Proc Inst Mech Eng H: J Eng Med. 2021 Apr 1; 235(4):437–446. https://doi.org/10.1177/0954411920987035https://doi.org/10.1177/0954411920987035

15 Delp SL, Anderson FC, Arnold AS, Loan P, Habib A, John CT, Guendelman E, Thelen DG. OpenSim: Open-Source Software to Create and Analyze Dynamic Simulations of Movement. IEEE Trans Biomed Eng. 2007 Nov; 54(11):1940–1950. https://doi.org/10.1109/TBME.2007.901024https://doi.org/10.1109/TBME.2007.901024

16 Fluit R, Andersen MS, Kolk S, Verdonschot N, Koopman HF. Prediction of ground reaction forces and moments during various activities of daily living. J Biomech. 2014 Jul 18; 47(10):2321–2329. https://doi.org/10.1016/j.jbiomech.2014.04.030https://doi.org/10.1016/j.jbiomech.2014.04.030

17 Skals S, Jung MK, Damsgaard M, Andersen MS. Prediction of ground reaction forces and moments during sports-related movements. Multibody Syst Dyn. 2017 Mar; 39(3):175–195. https://doi.org/10.1007/s11044-016-9537-4https://doi.org/10.1007/s11044-016-9537-4

18 Arjmand N, Shirazi-Adl A, Bazrgari B. Wrapping of trunk thoracic extensor muscles influences muscle forces and spinal loads in lifting tasks. Clin Biomech. 2006 Aug 1; 21(7):668–675. https://doi.org/10.1016/j.clinbiomech.2006.03.006https://doi.org/10.1016/j.clinbiomech.2006.03.006

19 Eltoukhy M, Travascio F, Asfour S, Elmasry S, Heredia-Vargas H, Signorile J. Examination of a lumbar spine biomechanical model for assessing axial compression, shear, and bending moment using selected Olympic lifts. J Orthop. 2016 Sep 1; 13(3):210–219. https://doi.org/10.1016/j.jor.2015.04.002https://doi.org/10.1016/j.jor.2015.04.002

20 Schmid S, Connolly L, Moschini G, Meier ML, Senteler M. Skin marker-based subject-specific spinal alignment modeling: A feasibility study. J Biomech. 2022 May 1; 137:111102. https://doi.org/10.1016/j.jbiomech.2022.111102https://doi.org/10.1016/j.jbiomech.2022.111102

21 McGill SM, McDermott A, Fenwick CM. Comparison of Different Strongman Events: Trunk Muscle Activation and Lumbar Spine Motion, Load, and Stiffness. J Strength Cond Res. 2009 Jul; 23(4):1148. https://doi.org/10.1519/JSC.0b013e318198f8f7https://doi.org/10.1519/JSC.0b013e318198f8f7

22 Marras WS, Ferguson SA, Lavender SA, Splittstoesser RE, Yang G. Cumulative Spine Loading and Clinically Meaningful Declines in Low-Back Function. Hum Factors J Hum Factors Ergon Soc. 2014 Feb; 56(1):29–43. https://doi.org/10.1177/0018720813496814https://doi.org/10.1177/0018720813496814

23 Kazemi Z, Mazloumi A, Arjmand N, Keihani A, Karimi Z, Ghasemi MS, Kordi R. A Comprehensive Evaluation of Spine Kinematics, Kinetics, and Trunk Muscle Activities During Fatigue-Induced Repetitive Lifting. Hum Factors. 2022 Sep 1; 64(6):997–1012. https://doi.org/10.1177/0018720820983621https://doi.org/10.1177/0018720820983621

24 Stokes IA, Gardner-Morse MG, Henry SM. Intra-abdominal pressure and abdominal wall muscular function: spinal unloading mechanism. Clin Biomech. 2010 Nov; 25(9):859–866. https://doi.org/10.1016/j.clinbiomech.2010.06.018https://doi.org/10.1016/j.clinbiomech.2010.06.018

25 Actis JA, Honegger JD, Gates DH, Petrella AJ, Nolasco LA, Silverman AK. Validation of lumbar spine loading from a musculoskeletal model including the lower limbs and lumbar spine. J Biomech. 2018 Feb 8; 68:107–114. https://doi.org/10.1016/j.jbiomech.2017.12.001https://doi.org/10.1016/j.jbiomech.2017.12.001

26 Favier CD, Finnegan ME, Quest RA, Honeyfield L, McGregor AH, Phillips AT. An open-source musculoskeletal model of the lumbar spine and lower limbs: a validation for movements of the lumbar spine. Comput Methods Biomech Biomed Engin. 2021 Sep 10; 24(12):1310–1325. https://doi.org/10.1080/10255842.2021.1886284https://doi.org/10.1080/10255842.2021.1886284

27 Molinaro DD, King AS, Young AJ. Biomechanical analysis of common solid waste collection throwing techniques using OpenSim and an EMG-assisted solver. J Biomech. 2020 May 7; 104:109704. https://doi.org/10.1016/j.jbiomech.2020.109704https://doi.org/10.1016/j.jbiomech.2020.109704

Appendix

Inverse kinematic (IK), inverse dynamic (ID), static optimization (SO) and joint reaction analysis (JRA)

Inverse kinematic (IK) was conducted using the OpenSim graphical user interface. Inverse dynamic (ID), static optimization (SO) and joint reaction analysis (JRA) were carried out using available MATLAB scripts provided by OpenSim. MATLAB scripts enabled the results analysis and graph generation.

Bar characterization

A dedicated experimental test bench comprising two rigid supports with a three-axis load cell (Mini85, ATI industrial automation, Michigan, US) screwed onto one of the supports was built (Figure 3a). The Russian bar was positioned on top of the rigid supports and the load cell. Loads of 40, 80, 120, 160 and 200 kg were successively applied at the middle of the Russian bar using straps and gym weights while the OptiTrak motion capture system (previously described) was used to track the position of the eighteen optical markers placed on the bar. This setup allowed for the measurement of the force acting on the support as a function of the bar deflection. The force at the porter’s shoulder was estimated by calculating the component of the force perpendicular to the bar’s shape.

In a second step, a dynamic test showed no significative impact of loading speed on structural deformation occurring during a static loading. Consequently, by monitoring bar deflection during jump activities using the eighteen optical markers, it was possible to calculate the force exerted on the porter’s shoulder using Equation 1.

Model modification and validation

One major modification regarding maximal isometric force helped to estimate spinal loads for Russian bar porters. It was observed that the model’s overall strength was inadequate for the intended purpose, as some calculations did not converge during the initial round of testing. Subsequent tests indicated that increasing the maximum isometric force enhanced the convergence of the calculations. Factors of two, three, five and ten were evaluated. According to our tests, the results did not vary significantly across the different factors applied. Consequently, a conservative factor of ten was applied to all maximum isometric forces in the model.

The model was validated using two different methods: a method using muscular activity comparison and a method using GRF comparison. In the first method, muscle activity recorded via EMG was compared phase-by-phase, across all jumps, to model-predicted muscle activations. Each muscle activity recorded with EMG was normalized to its peak maximal voluntary contraction (MVC), while the model’s muscle activation was the sum of the muscle fascicle activations normalized to their peak activation across each jump activity. The EMG validation relied on two criteria: the peak cross-correlation between the two signals signal using Matlab xcorr function13 and the agreement criteria proposed by Actis et al.25 The agreement criteria involved calculating the duration expressed as a percentage of each jump cycle when both EMG and muscle activity predicted by the model are either above or below a threshold of 0.5.25 The combined use of these criteria greatly enhanced the model’s validation credibility.

In the second method, ground reaction forces computed by the actuator-based method were directly compared to the force platform acquisitions during the propulsion and landing phases of candle jumps of less than two meters in height. Validation was assessed using two criteria: the root mean square error (RMSE) and the peak cross-correlation between the GRF total forces.

EMG validation

Figure 7 and Figure 8 display the normalized muscular activation values of eight muscle groups from the model estimations and experimental measurements for the propulsion, reception and stationary phases, averaged across all jump activities. Peak cross-correlation values ranged between 0.64 and 0.89, and agreement values ranged between 52% and 86% (Table 2).

The figure compares experimental electromyographic (EMG) signals (in red) with model predictions (in black) for several muscles (right rectus abdominis (RRA), left rectus abdominis (LRA), right external obliques (REO) and left external obliques (LEO)) across three movement phases. The graphs depict variations in muscle activation throughout the cycle, accompanied by a sequence of motion images at the bottom, illustrating body posture at key moments during the jump cycle.
Figure 7.

Mean muscular activation values recorded by EMG sensors (red) and estimated by the model (black) for four muscles groups: right rectus abdominis (RRA), left rectus abdominis (LRA), right external obliques (REO) and left external obliques (LEO). The dashed lines represent one half standard deviation from model values. The horizontal dotted lines represent the one-half threshold value for agreement calculation.

The figure compares experimental electromyographic (EMG) signals (in red) with model predictions (in black) for several muscles (right lumbar erector spinae (RLES), left lumbar erector spinae (LLES), right thoracic erector spinae (RTES) and left thoracic erector spinae (LTES)) across three movement phases. The graphs depict variations in muscle activation throughout the cycle, accompanied by a sequence of motion images at the bottom, illustrating body posture at key moments during the jump cycle.
Figure 8.

Mean muscular activation recorded by EMG sensors (red) and estimated by the model (black) for four muscles groups: right lumbar erector Spinae (RLES), left lumbar erector spinae (LLES), right thoracic erector spinae (RTES) and left thoracic erector spinae (LTES). The dashed lines represent one half standard deviation from model values. The horizontal dotted lines represent the one-half threshold value for agreement calculation.

Table 2.

Average peak cross-correlation and agreement values between muscle activation from model estimations and experimental measurements across all jump activities.

Jump

Validation

Propulsion

Landing

Stationary

Rectus Abdominis

Right

Correlation

0.73

0.77

0.74

Agreement

68%

67%

63%

Rectus Abdominis

Left

Correlation

0.76

0.68

0.77

Agreement

64%

66%

64%

External Obliques

Right

Correlation

0.82

0.83

0.83

Agreement

54%

59%

61%

External Obliques

Left

Correlation

0.75

0.77

0.81

Agreement

59%

55%

63%

Lumbar Erector Spinae

Right

Correlation

0.70

0.67

0.66

Agreement

58%

66%

63%

Lumbar Erector Spinae

Left

Correlation

0.77

0.75

0.78

Agreement

72%

70%

71%

Thoracic Erector Spinae

Right

Correlation

0.78

0.72

0.70

Agreement

64%

62%

61%

Thoracic Erector Spinae

Left

Correlation

0.89

0.85

0.64

Agreement

75%

52%

86%

The correlations between EMG and OpenSim muscular activation were predominantly high for most muscle groups, with the lowest peak cross-correlation at 0.64 and the highest reaching 0.86. In the original model validation by Beaucage-Gauvreau, peak cross-correlations were reported up to 0.93, and values above 0.82 were observed for all back muscle groups.13 The agreement values, falling within the range of 52 and 89%, were deemed acceptable. As emphasized by Actis et al.,25 agreement can be influenced by normalization and processing methods. Furthermore, the use of EMG sensors presented challenges due to noise generated by the impact of the Russian bar and the physical strain on the porters. This resulted in some sensors being displaced during the trials and noisy signal. In a broader context, our findings seem to be low compared to studies specifically dedicated to validating lumbar models,25–27 underscoring the need for an enhanced model meticulously tailored for circus artists given the highly specific demands of their activities.

GRF validation

Figure 9 illustrates the comparison between the estimated GRFs and the forces measured experimentally for both propulsion and landing phases of candle jumps (less than two meters in height). The RMSE for the propulsion phase was 191 N, with a mean difference of 13.6%. In the landing phases, there were more discrepancies, with an RMSE of 245 N and a mean difference of 19.9%. Despite differences between measured and estimated values, most results fall within one standard deviation, as supported by a cross-correlation of 0.99. Figure 10 shows that differences in the propulsion phases could be associated with a poor distribution of forces between feet from the actuator-based method. For the landing phases, there were fewer differences measured between feet, and the model seems to overestimate forces for both feet.

Two graphs showing a comparison between total estimated GRFs and total measured GRFs accompanied by a sequence of motion images at the bottom, illustrating the body posture at key moments during jump cycle.
Figure 9.

Comparison between estimated GRFs (red lines) and GRFs measured experimentally (green lines) on both force plates during propulsion and landing phases (two porters and eighteen jumps).

Two graphs showing a comparison between estimated GRF and measured GRF for the two different feet accompanied by a sequence of motion images at the bottom, illustrating body posture at key moments during jump cycle.
Figure 10.

Comparison between estimated GRF (solid lines) and measured GRF (dashed lines) for the right foot (red lines) and left foot (green lines) during propulsion and landing phases.