Energy and Power in (Sports) Biomechanics Department Center of for Sport Sensory-Motor and Exercise Interaction Science SPORTSCI 306 Technique Anvendt Assessment Biomekanik Uwe Uwe Kersting Kersting MiniModule Lecture 04 06-2007 2008 1 Objectives Review and learn basic energy formulae Determine total body energy by centre of mass calculations The use of power equations in performance assessment Apply energy estimations to an example of changes to sport equipment Learn about basic physics principles of rotational movements 2 Contents 1. Calculation of total body energy 2. Estimates of rotational energy 3. Estimates of energy consumption with varied running footwear 4. Energy calculations in high jump 5. Rotational movements 6. Power calculations in jump tests 3 1
Energy in a purely physics sense Mechanical energy: E kin = 0.5 m v 2 E pot = m g h E rot = 0.5 I ω 2 E el = 0.5 k x 2 Total body energy: n E tot = Σ E kin i + E pot i + E rot I : n = # of segments i = 1 4 CoM velocity? E kin tot = Σ m i v i 2 + potential Energy h E pot = Σ m i g h i Question: Is the total Energy equal to CoM based calculations? 5 COM calculation from segment positions x 0 = (Σ m i * x i ) M y 0 = (Σ m i * y i ) M 6 2
Rotational Energy contribution? E rot tot = Σ I i ω i 2 Example Leg: m leg = 0.0465 BM V leg sprint = 20 m/s I leg = 0.064 kg m 2 ω leg sprint = 500 /s 7 Mechanics Hook s law: Force is proportionate to deformation deformation: force: stiffness: Hook s: s [m] F [N] k [N/m] F s = - k * s F ext k s 8 Example: Running with and without shoes Oxygen consumption: + 4 5 % when running with shoes (6 7 min in a marathon!) Total work 10 7 J Shoe mass: 100 g Lifting height 0.2 m Maximal speed of foot during swing: 10 m/s Step length: 2 m => 20 ksteps Additional work due to gravity Additional work to accelerate the shoe 9 3
Add. Work [%] 10 5 Additional Energy needed 400g 300g 200g 100g 5 10 15 Max foot speed [m/s] 10 Application to high jump In high jump the initial energy has to be transformed into jump energy. 11 There is no linear relationship between approach velocity and vertical take-off velocity. (Dapena et al., 1990) Transformation is coupled with a decrease in total energy during the last ground contact. (Brüggemann & Arampatzis, 1991) 12 4
Purpose To examine the approach and take-off strategies of high jumpers at the world class level. To determine how to estimate the optimal take-off behavior from given initial characteristics. 13 Methods 14 Data Aquisition and Analysis four stationary PAL video cameras (50 Hz), two for each side, synchronised calibration with 2 x 2 x 3 m cube digitisation with Peak Motus (19 points) calculation of body angles, CM position, CM velocity and total energy by fast information program using DLT 15 5
Digitised Frames 16 Calculations V H = H 1 + 2 sin 2 α 2 g E p2 E k2 sin H = + 2 α g 2 g (1) (2) 17 where: E k 2 = E kin 2 m ; E p 2 = E pot 2 m T in = α E dec (3) 18 6
The effective height (H) can be expressed as a function of: initial energy, energy loss and transformation index. E p2 (E T1 - Edec - E p2) sin H = + 2 (Τin Edec) g g (4) 19 RESULTS 20 42 40 38 Cluster Analysis Group1 (n=16) Group2 (n=10) 36 E T2 [J/kg] 34 32 30 E T1 28 42 J/kg 26 40 J/kg 24 38 J/kg 32 J/kg 34 J/kg 22 36 J/kg 2 3 4 5 6 7 8 9 10 11 12 13 14 E dec [J/kg] 21 7
Parameters Group1 (n=16) Group2 (n=10) Initial energy [J/kg] 35.15 (1.38) 39.04 (2.06) * Energy decrease [J/kg] 4.26 (0.99) 9.23 (2.06) * Transformation index [ /J/kg] 12.05 (3.34) 5.47 (1.01) * Final energy [J/kg] 30.88 (1.14) 29.81 (0.74) Take-off angle [ ] 48.83 (1.99) 48.67 (2.94) 22 18 16 14 46 50 α [ ] Group1 (n=16) Group2 (n=10) T In [ /J/kg] 12 10 8 6 4 2 3 4 5 6 7 8 9 10 11 12 13 14 E dec [J/kg] 23 Two groups were identified that showed both varying initial conditions and jumping strategies. 24 8
Group1 showed a lower horizontal velocity and therefore, a lower initial energy prior to the last ground contact. Goup2 had a markedly increased initial energy and demonstrated a lower transformation index. 25 18 16 Formula: T in = exp(-c * E dec ) 14 T In [ /J/kg] 12 10 8 R 2 = 0.98 6 4 2 2 4 6 8 10 12 14 E dec [J/kg] 26 Height [m] 2,8 2,7 2,6 2,5 2,4 2,3 2,2 2,1 2,0 1,9 1,8 1,7 With that: 1,6 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 E dec [J/kg] Group1 (n=16) Group2 (n=10) 27 9
Parameter [m/s] Group1 (n=16) Group2 (n=10) Horizontal TD velocity 7.11 (0.29) 7.83 (0.43) * Decrease in h. velocity 3.22 (0.30) 4.05 (0.56) * Horizontal TO velocity 3.89 (0.27) 3.77 (0.22) Vertical TO velocity 4.44 (0.11) 4.29 (0.26) 28 Discussion 29 Both of the identified strategies resulted in comparable effective heights. The faster approaching athletes could not benefit from their higher energy produced during the run-up. 30 10
Since body orientation at touch-down shows no significant differences a possible explanation could be an altered muscle stiffness. The influence of changes in muscle stiffness has to be further investigated. 31 For both groups the optimum energy loss that leads to the maximum vertical velocity after the take-off can be estimated. 32 3,0 2,8 Sotomayor E T1 =41.07 J/kg Partyka E T1 =36.53 J/kg Forsyth E T1 =35.82 J/kg Matusevitch E T1 =36.65 J/kg 2,6 Jump height [m] 2,4 2,2 2,0 1,8 1,6 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 E dec [J/kg] 33 11
Training concepts may be developed to increase performance at both the world class level and in sub-elite jumpers. 34 35 Summary Relationships between kinematic parameters and jump performance are usually not very strong. Therefore, predictions are almost impossible. Energy approach provides much simpler analysis and can be used on a daily basis (?). Muscle stiffness/joint stiffness has to be included. 36 12
Initial knee angle, incoming speed and jump height Results from simulation study (Alexander, 1980) 37 Some basic rules on rotation 1. Gathering angular momentum 2. Regulating rotations 3. Aerials Summersaults with twist 4. Assessment of jump performances - simple tests - using a force platform 38 Aerials Off-centre force during take-off creates torque. The amount of angular momentum in the airborne phase cannot be changed. F Board F Com Resulting in: Momentum: M = m v (Impulse) F Rot Angular Mom.: A = I ω 39 13
How does a snowboarder/skier produce angular momentum? Virtually no friction to generate F rot! From TO angular momentum is preserved; CoM follows parabolic path (i.e., projectile motion) 40 Regulation of rotation Athletes have to change body configuration to alter I (Moment of inertia); Angular Momentum unchanged when airborne: A = I ω reducing I increases ω. Used for slowing down or accelerating rotations. How is forward rotation generated in ski jumps? 41 Aerials with twist If body position changes body rotates about a different local axis. Angular velocity is a vector along the rotational axis (right hand rule!) ω ω ω ω 42 14
Result (Yeadon 2000, etc.) a) in athlete s coordinate system b) in laboratory coordinate system 43 Jump Performance Common tests: - Jump and reach - Assumptions: - Jump belt test: Assumptions: 44 Determining Jump Height Common tests - Determine flight time: (contact mats or foot switches, video??) Assumptions: signal V flight time time 45 15
Impulse-Momentum Relationship Momentum (of a moving object): Mom = m * v(t) at any point in time t m v Impulse: I = F * t; actually: I = F * dt 46 In other terms BW impulse : BWI = t contact * BW Total Impulse: TI = F(t) dt = S Fi * Dt Resulting momentum: Mom = TI - BWI 47 Realisation in excel To calculate from measured force signal: Acceleration of CoM Velocity of CoM Displacement of CoM (all in vertical (z) direction; however applies for horizontal forces as well)..\labs\jumps_template.xls 48 16
F = m * a What else to get from GRF v(t) = s(t) = v0 + a dt s0 + v dt s(t) = s0 + v0 (t-t0) + a dt 2 49 Example: Result Counter Movement Jump Msubj [kg] 81.00 Caine hbox [m] 0.00 tcont [s] 0.70 Fmax [N] 2474.45 vtd [m/s] 0.00 vto [m/s] 3.37 maxp [W/kg] 74.00 hj [m] 0.57 neg Work[J/kg] -2.60 pos Work[J/kg] 7.11 gain [J/kg] 9.71 Force [N] 3000 2500 2000 1500 1000 500 0 0 0.2 0.4 0.6 0.8 vcm [m/s] 4 2 0 0 0.2 0.4 0.6 0.8-2 W [J/kg] 10 5 0 0 0.2 0.4 0.6 0.8-5 P [W/kg] 100 50 0 0 0.2 0.4 0.6 0.8-50 -100 50 17