This example shows how to work with units in physics calculations. Calculate the terminal velocity of a falling paratrooper in both SI and imperial units. Solve the motion of the paratrooper, taking into account the gravitational force and the drag force.
Imagine a paratrooper jumping out of an airplane. Assume there are only two forces acting on the paratrooper: the gravitational force and an opposing drag force from the parachute. The drag force is proportional to the velocity squared of the paratrooper.
The net force acting on the paratrooper can be expressed as
is the mass of the paratrooper
g is the gravitational acceleration
is the velocity of the paratrooper
is the drag constant
Define the differential equation describing the equation of motion.
syms g m c_d syms v(t) eq = m*diff(v(t),t) + m*g == c_d*v(t)^2
Assume that the parachute opens immediately at so that the equation
eq is valid for all values of . Solve the differential equation analytically using
dsolve with the initial condition . The solution represents the velocity of the paratrooper as a function of time.
velocity = simplify(dsolve(eq, v(0) == 0))
Find the SI unit of the drag constant .
The SI unit of force is the Newton . In terms of the base units, the Newton is . Since these are equivalent, they have a unit conversion factor of 1.
u = symunit; unitConversionFactor(u.N, u.kg*u.m/u.s^2)
The drag force must have the same unit in Newton as the gravitational force . Using dimensional analysis, solve for the unit of .
syms drag_units_SI drag_units_SI = simplify(solve(drag_units_SI * (u.m / u.s)^2 == u.N))
Describe the motion of the paratrooper by defining the following values.
Mass of the paratrooper
Substitute these values into the velocity equation and simplify the result.
vel_SI = subs(velocity,[g,m,c_d],[9.81*u.m/u.s^2, 70*u.kg, 40*drag_units_SI])
vel_SI = simplify(vel_SI)
Compute a numerical approximation of the velocity to 3 significant digits.
digits(3) vel_SI = vpa(vel_SI)
The paratrooper approaches a constant velocity when the gravitational force is balanced by the drag force. This is called the terminal velocity and it occurs when the drag force from the parachute cancels out the gravitational force (there is no further acceleration). Find the terminal velocity by taking the limit of .
vel_term_SI = limit(vel_SI, t, Inf)
Finally, convert the velocity function from SI units to imperial units.
vel_Imperial = rewrite(vel_SI,u.ft)
Convert the terminal velocity.
vel_term_Imperial = rewrite(vel_term_SI,u.ft)
To plot the velocity as a function of time, express the time
t in seconds and replace
T s, where
T is a dimensionless symbolic variable.
syms T vel_SI = subs(vel_SI, t, T*u.s)
vel_Imperial = rewrite(vel_SI, u.ft)
Separate the expression from the units by using
separateUnits. Plot the expression using
fplot. Convert the units to strings for use as plot labels by using
[data_SI, units_SI] = separateUnits(vel_SI); [data_Imperial, units_Imperial] = separateUnits(vel_Imperial);
The velocity of the paratrooper approaches steady state when . Show how the velocity approaches terminal velocity by plotting the velocity over the range .
subplot(1,2,1) fplot(data_SI,[0 2]) title('Velocity in SI Units') xlabel('Time in s') ylabel(['Velocity in ' symunit2str(units_SI)]) subplot(1,2,2) fplot(data_Imperial,[0 2]) title('Velocity in Imperial Units') xlabel('Time in s') ylabel(['Velocity in ' symunit2str(units_Imperial)])