67
Numerical Simulation of Skydiving Motion*
Our goal for this chapter is to understand how we created the previously shown graphs of acceleration, velocity, and position of our example skydiver, even though the net force and acceleration changes throughout. We will use a numerical simulation that ties together just about everything we have learned so far in this unit to achieve this goal. We already know that the initial velocity is zero and therefore the initial drag force is zero. With no drag force in the first moment of the jump, the diver is in free fall and the acceleration is just g in the downward direction, or -9.8 m/s/s . We can then calculate the velocity after a short time interval as:
We have made theassumption that the acceleration during this interval was constant, even though it wasn’t, but if we choose a time interval that is very small compared to the time over which the acceleration changes significantly, then our result is a good approximation. A time interval of one second will satisfy this condition in our case, so we now calculate the velocity at the end of the first one second interval:
Now that we have a velocity we can calculate the air resistance at the start of the second interval using our previously stated values for human drag coefficient, cross-sectional area, and the standard value for air density:
Now that we have a drag force due to air resistance we can use Newton's Second Law to calculate the acceleration at the start of the second interval. We have only two forces, drag and gravity and we will use our previously stated skydiver mass of 80 kg:
Now we just continue this iterative process of using acceleration and velocity values from the previous interval to calculate new velocity, drag force, and acceleration for next interval.
Using the data produced by the simulation we can graph the drag force. Showing the weight on the same graph we can see how the drag force approaches the weight.
We can also use the data to create motion graphs for the skydive and see that the acceleration gradually transitions from -9.8 m/s/s to zero as drag force increases.
We see that velocity is always negative and the speed is always increasing, but the slope becomes less steep because the acceleration is decreasing with time:
Finally we can see that the position graph eventually becomes linear as terminal velocity is reached. (Note that we have converted our initial position of 12,000 ft to the equivalent 3660 m)