There are three distinct parts to this exercise, all of which are simple extensions of the previous exercise (hpvloop), and they should be solved and debugged separately:
1. A data file named Cheetah.data contains
the basic private data of the Cheetah
hpv, which broke the world speed record for a human powered vehicle
in 1992 - 68.7mph! We wish to write a function getbike which
will get the data from the file, and use the Constructor of class
HPvehicle to assign this data to the private variables of our
hpv Object. The format of the file Cheetah.data follows:
Notice that the file has three data items followed
by six lines of explanatory information on the data. The function
getbike should both get the data into the respective variables
and then read and display all the remaining lines of information.
Refer to the week7 directory for details.
2. Extend the root finding technique to enable you to determine the steady state velocity attained if you apply a specified power to the pedals. Essentially this requires reformulating the power equation that we developed in a previous exercise in terms of a difference function as follows:
Thus the required velocity is the root of the difference function Δ(V). This difference function should be interpreted as follows: "The value of V that makes the difference function Δ(V) = 0 is the value of the steady state velocity resulting from a steady power P applied at the pedals." Note that function find_root has now been renamed find_vel, and has been augmented to accept the required power power as an argument. It represents the value of the steady power applied to the pedals, and is read from the keyboard in the main function module to be used in the new function dif_power, as shown. Note also that you need not rewrite the function find_power since the new function dif_power will invoke function find_power without change, as shown in the following structure diagram of hpvfile:
3. The third extension of program 5 (hpvloop) will let the program create a datafile of (velocity, power) values in your home directory which will subsequently be used for off-line plotting. Thus you will modify the main function to enable creating, writing to, and closing a data file using the various file commands that we learn about in lab7. This second part of the exercise will give you experience in transferring files from condor to your local computer using the SFTP (Secure File Transfer Protocol) program, editing local textfiles with the Wordpad program, and plotting data using the MATLAB application program.
The following is a complete typical program execution output:
get bike data from the "Cheetah.data" file
data read: 95.000, 0.046, 0.005
=============================
Cheetah data, Colorado, 1992.
total mass = 95.0 kg
CdA = 0.046 sq.m
Cr = 0.005
=============================
hpv initialised as follows:
mass (hpv + rider): 95.000[kg]
cda (coeff.drag*area): 0.046[sq.m]
cr (coeff.rolling resist): 0.005
local conditions initialised as follws:
gravity acceleration: 9.807[m/s/s]
air density: 1.180[kg/cu.m]
wind velocity: 0.000[mph]
slope: 0.000
-----------------------------------------------------
enter wind velocity [mph] - positive for headwind
-1
value entered is -1.000[mph]
enter slope (height/distance) - positive for uphill
-0.02
value entered is -0.020
new local conditions:
wind velocity = -1.000[mph]
slope = -0.020
evaluate power vs velocity
enter low and high velocity bounds [mph]
0 70
values entered are [0.000,70.000] [mph]
enter the number of points to evaluate
20
number of points is 20
velocity (mph) power (W)
0.000 -0.000
3.684 -22.820
7.368 -45.054
11.053 -65.988
14.737 -84.906
18.421 -101.094
22.105 -113.837
25.789 -122.420
29.474 -126.128
33.158 -124.247
36.842 -116.060
40.526 -100.854
44.211 -77.913
47.895 -46.523
51.579 -5.968
55.263 44.466
58.947 105.494
62.632 177.831
66.316 262.192
70.000 359.292
20 rows of data written to "hpv.out" file
find velocity as a function of applied power
enter applied power [watts]
60
enter required accuracy in velocity [mph]
0.1
required accuracy is 0.100[mph]
for applied power of 60.000[watts]
steady state velocity is 56.294[mph]
|
Notice the various items of data entry to this program shown in Italics script. Your teacher has used his own input data for this example, being a one mph tailwind and a downhill slope of 0.02 (or 2%). Even though he was only able to produce 60 watts of power (typical of a laid back professor), he was still able to reach a steady state velocity of 56mph! Use your own data for the execution run, as specified in the first programming exercise (hpv)
Notice the line in which it is stated that 20 lines of data have been written to the hpv.out file. This file will be saved in your home directory, and contains the data that we wish to use for plotting. Each column of data in the datafile is separated by a horizontal tab symbol, so that it can be read into a spreadsheet or graphing program.
A printout of the datafile hpv.out is shown below.
0.000 -0.000
3.684 -22.820
7.368 -45.054
11.053 -65.988
14.737 -84.906
18.421 -101.094
22.105 -113.837
25.789 -122.420
29.474 -126.128
33.158 -124.247
36.842 -116.060
40.526 -100.854
44.211 -77.913
47.895 -46.523
51.579 -5.968
55.263 44.466
58.947 105.494
62.632 177.831
66.316 262.192
70.000 359.292
|
Once the datafile is available in your home directory then you should use the SFTP program to "fetch" the file and save it on the local computer disk. At this stage you should be able to examine the file, and any out-of-range data lines can be edited out. Finally you should read the file into the MATLAB application and plot power vs velocity, as shown below. The final graph can be printed on the Laserprinter in the lab. Notice on the graph that the horizontal line representing the actual applied power is included, allowing one to read off the graph the value of velocity attained at that power.
As before, your program will be separated into
two files, the program file hpvfile.cpp and the header file hpvfile.h. Both
of these files should be in your home directory by 10:00 am of
the due date. Furthermore, the printout of the graph should be
submitted in the box outside Stocker Room 265 by 10:00 am of the
due date.
