Grand Prix Racing - | The Science of Fast Pinewood Cars |
Mathematical modeling is a great way to capture important behaviors of physical systems. Much work goes into researching and understanding how to construct the model. But once that has been done, the concepts that are encoded in its algebra and equations of motion can be used to make measureable predictions. Those predictions, in turn, can be verified by experiment and used to either make better physical systems or change the model to reflect new understanding.
What your car does on the track is called behavior. The behavior of a model describes what happens under what circumstances. Sometimes, behavior is more than a description. It can be measured. Not all behavior is easily measured. Just think about drift. But much is. That behavior which is measureable, we want to model.
What kind of behavior is measureable? Certainly, speed and finish time are measureable in a real race with the right physical instruments. But so are all the forces on the car and its energy uses. But these measurements would be difficult to make on a real car. They and many others can be measured virtually through modeling.
Examination of a Grand Prix car's physical behavior on a Grand Prix track yielded three primary models and a general energy equation that can be used with arbitrary track shapes in computer simulations.
Though most of the physical aspects of the race were examined to produce these models, some effects were found too small to make a difference. For example, the rotational inertia of a car body though examined, was not included in the closed dynamic parameters model. On the other hand, some aspects can not be modeled in a simple closed system but require iterative solutions on a case-by-race basis. Race dynamics on a curved ramp/flat transition are one example.
|
Many kinds of measurements are needed to model a Grand Prix race. These parameters describe the track, your car and the physics between them. The following tables summarize a "simple" behavior model for a Grand Prix race. This model was constructed from two parts. A ramp part and a flat part as shown below. The transistion between the two is an angular join (that is, not a smooth transition) so its effect on the behavior of the race can be minimized. Equations of motion for general angular joins between ramps and catenary transitions are developed in this manual to study transition effects by computer.
In order to present a physically correct model below, the energy lost at the abrupt transition has been included as a cosO factor in the equation for the initial speed on the flat, v0. If you want to use this summary to calculate dynamics closely simmulating a smooth transition track, set that cosO term to the value one (1.0) - only in the v0 equation. The energy loss at the abrupt angular join can be quite high, but on a smooth transition the energy loss (due to ill positioned weight and centrifugal force) is much less.
Below, the most relevant parameters, equations and relationships are presented.
Though air density and gravitational acceleration are not strickly constants, they are treated as such in this model.
SYMBOL | UNITS | VALUE | DESCRIPTION |
---|---|---|---|
p | ozs2/in4 | 0.000001877865 | Density of air, 1 atm, 60 F |
g | in/s2 | 386.088 | Acceleration due to gravity at sea level |
The first four track parameters should be changed to be those of your Grand Prix track. The rest are defined in terms of these. Measurements from the Official AWANA Grand Prix track were used as default values.
SYMBOL | UNITS | DEFAULT | DESCRIPTION |
---|---|---|---|
O | degrees | -21.1781 | Ramp angle, measured from the base; negative for down hill |
H | inches | 39.9 | Ramp height at starting line |
Lr | inches | 110.5 | Length of ramp surface |
Lf | inches | 249.5 | Length of flat surface |
h(y) = H + ysinO | in | 0 - 39.9 | Height of ramp in terms of distance along the ramp |
0 <= y <= Lr | in | 0 - 110.5 | Range of arc parameter, distance along the ramp |
h(y) = 0 | in | 0 | Height of flat part of track |
0 <= y <= Lf | in | 0 - 249.5 | Range of arc parameter along the flat (it starts over at 0) |
Some of the following car parameters are easily measured. Others like the moment of inertia of a wheel are not. This manual contains means by which they can be determined. Default values reflect those of AWANA Grand Prix cars.
SYMBOL | UNITS | DEFAULT | DESCRIPTION |
---|---|---|---|
a | scalar | 0.5 | Aerodynamic drag coefficient |
A | in2 | 1.5 | Frontal area |
m | ozs2/in | 0.01295 | Total mass (5 oz) |
b | ozs2/in | 0.01209 | Mass of car body (total - wheels) |
mf | ozs2/in | 0.000216 | Mass of a front wheel (0.0833 oz) |
mr | ozs2/in | 0.000216 | Mass of a rear wheel (0.0833 oz) |
If | ozins2 | 0.0000486 | Moment of inertia of a front wheel |
Ir | ozins2 | 0.0000486 | Moment of inertia of a rear wheel |
nf | scalar | 0 - 1.520 (AWANA kit axles < 0.2) | Drag coefficient of a front axle |
nr | scalar | 0 - 1.520 (AWANA kit axles < 0.2) | Drag coefficient of a rear axle |
uf | in | 0 - 0.02 (AWANA kit wheels on sanded track < 0.15) | Coefficient of front rolling (tread) friction |
ur | in | 0 - 0.02 (AWANA kit wheels on sanded track < 0.15) | Coefficient of rear rolling (tread) friction |
Rf | in | 0.59 | Outer radius of front wheels |
Rr | in | 0.59 | Outer radius of rear wheels |
rf | in | 0.049 | Inner radius of front wheel bore |
rr | in | 0.049 | Inner radius of rear wheel bore |
CMy | in | 0 - 5.8 | Center of mass distance above rear axle normal to the car's base |
CMx | in | -0.59 - 2.41 | Center of mass distance in front of rear axle parallel to the car's base |
N | in | 0 - 5.2 (1.5 kit) | Length of car nose (front to front axle) |
B | in | 0 - 5.8 (3.875 kit) | Distance between the front and rear axles |
* | rad | -pi/2 - pi/2 (0 kit) | Angle between track and a plane across the front and rear axles |
Ff = bg(CMx-CMytan(O+*))/B | oz | 0 - 4.8334 | Weight supported by a front axle (O=0 on flat) |
Fr = bg(B-CMx+CMytan(O+*))/B | oz | 0 - 4.8334 | Weight supported by a rear axle (O=0 on flat) |
k = apA/2 | ozs2/in2 | 0 - 0.000008 | Air mass displacement per inch of track |
A few of the ramp parameters like M, D and D1 are defined to make the expressions clearer and shorter. All are defined in terms of the track and car parameters above. Their values change with those parameters, so they have no defaults. Each expression has a meaningful physical interpretation.
Work of aerodynamic drag
/ y Wa = | apAv(y)2dy/2 / 0
EXPRESSION | UNITS | DESCRIPTION |
---|---|---|
M = m + 2If/Rf2 + 2Ir/Rr2 | ozs2/in | Inertial Mass |
Dr = (nfFfrf/Rf + 2uf(Ff/2+mfg) + nrFrrr/Rr + 2ur(Fr/2+mrg))cosO | oz | Total axle and tread drag using Ff and Fr computed on the ramp |
E = -mgLsinO = -mg(L-y)sinO + Mv2/2 + Dry + Wa | ozin | Ramp Energy Relation |
F = M a = -mgsinO - Dr - kv2 | oz | Equation of motion along a ramp |
vt = \[-(mgsinO+Dr)/k] | in/s | Terminal velocity on ramp |
Lm = M/k | in | Distance needed to move air mass equal to car's inertial mass |
T = Lm/vt | s | Time needed to move air mass equal to car's inertia at terminal velocity on the ramp |
t(y) = Tarccosh(exp(y/Lm)) | s | Time it takes the car to travel a distance, y, from the starting gate |
v(t) = vt tanh(t/T) | in/s | Car's speed at time, t, after the start |
v0 = vt \[1-exp(-2L1/Lm)] | in/s | Speed at end of ramp (no energy loss due to abrupt join) |
a(t) = (mgsinO+Dr) (tanh(t/T)2 - 1)/M | in/s2 | Car's acceleration in terms of time, t |
y(t) = Lm ln(cosh(t/T)) | in | Distance along ramp at time, t |
Note: arccosh(X) = ln(X+\[XX-1])
Note: tanh(X) = (e^X - e^-X)/(e^X + e^-X)
Some of the symbols used in the flat track expressions are ramp parameters. All are in terms of track and car parameters. Note the expressions that measure coasting behaviors.
EXPRESSION | UNITS | DESCRIPTION |
---|---|---|
v0 = vt \[1-exp(-2L1/Lm)]cosO | in/s | Speed at start of flat accounting for energy lost at the join |
Df = nfFfrf/Rf + 2uf(Ff/2+mfg) + nrFrrr/Rr + 2ur(Fr/2+mrg) | oz | Total axle and tread drag using Ff and Fr computed on the flat (in particular, O = 0) |
E = Mv02/2 = Mv2/2 + Dfy + Wa | ozin | Energy Relation |
F = M a = -Df - kv2 | oz | Equation of motion for coasting along flat |
vi = \[Df/k] | in/s | Speed of displaced air with energy matched to work done on car by drag |
Tf = Lm/vi | s | Time needed to displace an air mass equal to the car's inertial mass at the speed, vi |
yc = -Lm ln(cos(tc/Tf)) | in | Coasting distance on flat |
tc = Tfarctan(v0/vi) | s | Coasting time on flat |
t(y) = tc - Tf arccos(exp((y-yc)/Lm)) | s | Time to coast a distance, y, on flat |
v(t) = vi tan((tc-t)/Tf) | in/s | Speed at time, t, on flat |
a(t) = -Df(tan((tc-t)/Tf)2 + 1)/M | in/s2 | Car's acceleration in terms of time, t |
y(t) = yc + Lm ln(cos((tc-t)/Tf)) | in | Distance along flat at time, t |
For race results, the ramp and flat equations must be used properly. The car's trajectory must be used on the ramp and flat. This means using trajectory lengths instead of the track lengths, but the same angle for the ramp. Parameters from the track, car, ramp and flat are all used.
EXPRESSION | UNITS | DESCRIPTION |
---|---|---|
L1 = Lr + (N+B-CMx)cos* + CMysin* + (Rr + CMycos* + CMxsin*)tan(O/2) | in | Length of trajectory on ramp |
L2 = Lf - (N+B-CMx)cos* - CMysin* + (Rr + CMycos* + CMxsin*)tan(O/2) | in | Length of trajectory on flat |
t = T arccosh(exp(L1/Lm)) + tc - Tf arccos(exp((L2-yc)/Lm)) | s | Time at finish line |
v = vi \[exp(-2(L2-yc)/Lm)-1] | in/s | Speed at finish line |
a = -D(tan((tc-t)/Tf)2 + 1)/M | in/s2 | Deceleration at finish line |
vmax = v0 | in/s | Maximum speed (bottom of ramp) |
amax = -(mgsinO+Dr)/M | in/s2 | Maximum acceleration (starting line) |
With these expressions, a virtual race between two or more cars can be run. Each car can be placed using its time. Intervals between cars are estimated from time difference and the slower speed. That is,
interval = time_difference*slower_speed
Here the time_difference is that between the fastest car and one of the slower ones in the virtual race (slower time minus fastest time). The slower_speed is that of the slower car. So the interval is how far behind the slower car was when the fastest car crossed the finish line. This is a good estimate upto a couple of feet, unless one of the cars is reaching its coasting distance limit and almost stopping on the finish line.
By breaking up the simplified Grand Prix track into two parts, we get a close look at the differences in accelerations, speeds and times and at the differences in the math it took to get there. The equations are very similar with hyperbolic versions of the functions occuring on the ramp part and trigonometric ones on the flat part. Also, because of the initial velocity on the flat track, we got the coasting distance and time with no extra effort.
Notice that if the friction was so high that the coasting distance had been less than the length of the track, our equations would have "broken down" and the "real" car would have stopped on the track. The time function for the finish line would have broken down. Its arccos function has exp((L2-yc)/LM) as an arguement. L2-yc would have been greater than zero making the exp function greater than one. The arccos function is not defined for numbers greater than one. Therefore, the time of the race could not be determined with that expression; the race time is infinite! There are other situations like this where an expression is not valid. For reasonable values of the parameters, this is not a problem. But knowing what "reasonable" is may be.
Is this track model realistic? Yes, it is actually quite good. When used without the correction for energy loss at the abrupt join to simulate a smooth transition (which most tracks have) it matches kinematic computer results and actual race times pretty closely. Since the many factors and behaviors of drift are not included in this model, one expects the calculated times to be less than actual race times. The times might also be less because a smooth transition generates centrifugal force and affects the distribution of this higher axle force to the front and rear wheels differently depending on the balance point. The difference in times can provide clues about the importance of drift and the slightly increased transition friction in the race.
In addition, some factors were not included in this model because they were found to make very little difference and would have required more terms and complexity to calculate. They are examined in the content of this manaul. All of these model expressions and the two other models mentioned above are derived in this manual also.
There is help for evaluating these model equations.
[Pit Area] | [Title Page] |
Grand Prix Racing - | The Science of Fast Pinewood Cars |
Copyright © 1997, 2004 by Michael Lastufka, All rights reserved worldwide. |