Vertical Jump Simulation

Simulates vertical jump using Runge-Kutta method

vj_simulate(
  mass = 75,
  weight = mass * gravity_const,
  push_off_distance = 0.4,
  gravity_const = 9.81,
  time_step = 0.001,
  save_trace = TRUE,
  fgen_func = fgen_get_output,
  iter = FALSE,
  max_iter = 1000,
  ...
)

Arguments

mass	Numeric value. Default is 75kg
weight	Numeric value. Default is mass * gravity_const
push_off_distance	Numeric value. Default is 0.4m
gravity_const	Numeric value. Default is 9.81
time_step	Numeric value. Time step used in simulation. Default is 0.001
save_trace	Logical. Default is TRUE
fgen_func	Function used to represent Force Generator. Default is fgen_get_output.
iter	Logical. Default is FALSE
max_iter	Numeric value. Default value 1000
...	Forwarded to fgen_func

Value

List object with summary data frame and trace data frame

Examples

vertical_jump <- vj_simulate(
  mass = 85,
  push_off_distance = 0.4,
  time_step = 0.001
)

round(t(vertical_jump$summary), 3)
#>                                           [,1]
#> mass                                    85.000
#> weight                                 833.850
#> push_off_distance                        0.400
#> gravity_const                            9.810
#> take_off_distance                        0.398
#> take_off_time                            0.353
#> take_off_ground_reaction_force          57.803
#> take_off_propulsive_force             -776.047
#> take_off_velocity                        2.366
#> height                                   0.285
#> mean_GRF_over_distance                1431.465
#> mean_GRF_over_time                    1403.633
#> mean_velocity                            1.128
#> work_done                              570.006
#> impulse                                201.134
#> mean_power                            1614.747
#> mean_RFD                             13244.599
#> mean_RPD                             10156.443
#> peak_GRF                              1946.956
#> peak_GRF_time                            0.147
#> peak_GRF_distance                        0.028
#> peak_velocity                            2.416
#> peak_velocity_time                       0.341
#> peak_velocity_distance                   0.369
#> peak_power                            2945.368
#> peak_power_distance                      0.252
#> peak_power_time                          0.290
#> peak_RFD                             16527.988
#> peak_RFD_distance                        0.005
#> peak_RFD_time                            0.097
#> GRF_at_100ms                          1515.783
#> GRF_at_200ms                          1760.750
#> peak_RPD                             27244.399
#> peak_RPD_distance                        0.019
#> peak_RPD_time                            0.133
#> diff_peak_to_take_off_velocity           0.049
#> ratio_peak_to_take_off_velocity          1.021
#> diff_mean_to_take_off_velocity          -1.238
#> ratio_mean_to_take_off_velocity          0.477
#> diff_mean_to_peak_velocity              -1.288
#> ratio_mean_to_peak_velocity              0.467
#> diff_mean_to_peak_power              -1330.621
#> ratio_mean_to_peak_power                 0.548
#> diff_mean_GRF_over_distance_to_time     27.832
#> ratio_mean_GRF_over_distance_to_time     1.020
#> mean_velocity_as_TOV_half                1.183

plot(
  x = vertical_jump$trace$distance,
  y = vertical_jump$trace$velocity,
  type = "l"
)

plot(
  x = vertical_jump$trace$time,
  y = vertical_jump$trace$ground_reaction_force,
  type = "l"
)

plot(
  x = vertical_jump$trace$distance,
  y = vertical_jump$trace$ground_reaction_force,
  type = "l"
)