Project Report Sample
Project Report Sample
Project Report
by
17 February 2014
CONTENTS
Contents
1 Statement of the Problem
2 Mathematical Modeling
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
3
3
3
4
5
5
4 Conclusion
Bibliography
Appendices
2 Mathematical Modeling
The Newtons second law of motion forms the basis of mathematical modeling for this problem (?). Its mathematical expression is normally written
as
F = ma
(1)
where F is the net force acting on the body (N or kg m/s2 ), m is the mass of
the object (kg) and a is its acceleration (m/s2 ).
The second law can be recast in the format of
a=
F
m
(2)
2 MATHEMATICAL MODELING
The Newtons second law of motion can also be used to determine the
terminal velocity of a free-falling body near the earths surface. In our problem the falling body is the bungee jumper, Figure 1.
A mathematical model can be derived by expressing the acceleration
as the time rate of change of the velocity (dv/dt), and substituting it into
Eq. (2) to yield
dv
F
=
dt
m
(3)
where v (m/s) is velocity. Thus, the rate of change of the velocity is equal
to the net force acting on the body normalized to its mass. If the net force
is positive, the object will accelerate. If it is negative, the object will decelerate. If the net force is zero, the objects velocity will remain the same at a
constant level.
We can now express the net force in terms of measurable variables and
parameters. For a body falling within the vacinity of the earth, the net force
is composed of two opposing forces: the downward pull of gravity FD and
the upward force of air resistance FU , see Figure 1:
F = FD + FU
(4)
If force in the downward direction is assigned positive sign, the second law
can be used to formulate the force due to gravity as
FD = mg
(5)
(6)
where cd is the drag coefficient (kg/m). Thus, the greater the fall velocity,
the greater the upward for due to air resistance. The parameter cd account
for properties of the falling object, such as shape or surface roughness, that
affect air resistance. For this case, cd might be a function of the type of
clothing or the orientation used by the jumper during free fall.
The net force is the difference between the downward and upward force.
Therefore, Eqs. (3) through (6) can be combined to yield
c
dv
= g d v2
dt
m
(7)
Eq. (7) is the mathematical model for this problem which relates the acceleration of a falling object to the forces acting on it. It is a differential equation
because it is written in terms of the differential rate of change (dv/dt) of the
variable that we are interested in predicting.
e x e x
e x + e x
(9)
Analytical solution to the problem may be obtained by inserting parameters into Eq. (8) to yield
!
r
r
9.81 68.1
9.81 0.25
v( t) =
tanh
t = 51.6928 tanh (0.18977t)
0.25
68.1
which can be used to compute the results in the table below:
Time t (s)
0
2
4
6
8
10
12
Velocity v (m/s)
0
18.7292
33.1118
42.0762
46.9575
49.4214
50.6175
51.6938
(10)
where v and t are differences in velocity and time computed over finite
intervals, v(ti ) is velocity at an initial time ti , and v(ti+1 ) is velocity at some
later time ti+1 . Note that dv/dt v/t is approximate because t is
finite.
Remember from calculus,
v
dv
= lim
dt
t0 t
Eq. (10) represents the reverse process and it can be substituted into Eq. (7)
to give
v ( t i +1 ) v ( t i )
c
= g d v ( t i )2
t i +1 t i
m
This equation can be rerranged to yield
i
h
c
(11)
v ( t i +1 ) = v ( t i ) + g d v ( t i )2 ( t i +1 t i )
m
The term in brackets of the above equation is the RHS of the differential
equation itself Eq. (3). That is, it provides a means to compute the rate of
change or slope of v. Thus, the equation can be rewritten as
dvi
t
(12)
dt
where the nomenclature vi designates velocity at time ti and t = ti+1 ti .
Eq. (12) can be solved using the Eulers method.
v i +1 = v i +
3.2.1
Algorithm
Flowchart
The flowchart for the algorithm developed to solve the free-falling bungee
jumper above is shown in Figure 2.
3.2.3
Program Code
The Matlab script, see Palm (2011) and Chapra & Canale (2009), to compute
the velocity of the free-falling bungee jumper as given below:
% bungee.m
% Compute the free-fall velocity of an object
%
% Input:
% g = acceleration due to gravity
% m = mass
% cd = drag coefficient
%
g = 9.81;
m = 68.1;
cd = 0.25;
% Compute velocity after 12 s
disp( )
fprintf(Velocity after 12 s is %8.4f m/s\n, ...
4 CONCLUSION
sqrt(g*m/cd) * tanh(sqrt(g*cd/m)*12))
% Plot velocity againts time
t = [0:1:30];
v = sqrt(g*m/cd) * tanh(sqrt(g*cd/m)*t);
plot(t,v,t,v,o)
title(Plot of v versus t);
xlabel(Values of t)
ylabel(Values of v)
grid
was saved into the file bungee.m and run at the Matlab prompt thus:
>> bungee
The result was displayed as
Velocity after 12 s is
50.6175 m/s
together with the plot of velocity versus time. This plot of the results computed over a period of 30 s is shown in Figure 3.
Plot of v versus t
60
50
Values of v
40
30
20
10
10
15
Values of t
20
25
30
4 Conclusion
According to the model, Eq. (7), the jumper accelerates rapidly, Figure 3.
A velocity of 50.6175 m/s is attained after 12 s. Note also that after a
4 CONCLUSION
Bibliography
1. C HAPRA , S. C. & C ANALE , R. P. (2009): Numerical Methods for Engineers,
6ed, ISBN 0-39-095080-7, McGraw-Hill
2. M AYO , W. & C WIAKALA , M. (1994): Schaums Outline of Programming With
Fortran 77, ISBN 978-0070411555, McGraw-Hill
3. PALM , W. J. (2011): Introduction to MATLAB for Engineers, 3ed, ISBN 978-007-353487-9, McGraw-Hill
4. R AO , S. S. (2002): Applied Numerical Methods for Engineers and Scientists,
ISBN 0-13-089480-X, Prentice Hall
APPENDICES