Integration of the Equations of Motion

The system is described by a single degree of freedom, , the distance along the track.
However, the free body diagrams give the second derivative of , the acceleration. This
requires both the position and velocity to be calculated as state variables, since the integration
routine to be used will solve for the derivative of the state variables. Thus, the acceleration will
be calculated from the integration routine, and used with the equations produced from the free
body diagram to calculate the forces on the car. The state variables are intoduced below.

A Runge-Kutta-Fehlburg integration routine will be used to solve for the state variables. The required derivatives of the state variables are shown here.

The value of is calculated by means of the free body diagram presented in section 3.1. The component of the forces and accelerations in the direction are equated.

Recall that and are perpendicular to . Evaluating the dot products results in the following.

Equation 3.6 is now solved for .

The result shows that must also be calculated. Further analysis of the free body diagram is made. First the components of the forces and accelerations in the direction are found.

Evaluation of the dot product results in the following.

The components of forces and accelerations in the direction are now determined.

The dot product is then evaluated.

And finally, the components in the direction are found.

The dot product is evaluated to produce equation 3.10.

Equations 3.6, 3.8, 3.9, and 3.10 can be combined to determine . First, equations 3.8 and 3.9 are combined to eliminate
. The result is then combined with equation 3.10 to eliminate . Finally, equation 3.6 was used to eliminate . The resulting seat force equation is shown here.

Many quantities in the equations for and are geometric variables, and depend on which element the car is traveling in. These include ,
, and
. The independent variable, , is used to calculate the distance into the current element, as shown in equation 3.12. Let represent the distance into the element, and let represent the total length of that element.

The values of the state variables come from the integration routine. The values of and used in the calculation of the current and are those from the previous iteration. The time step used should be small, as a smaller time step will result in more accurate results. The limiting factor in this will typically be the computer on which the program is run. Many data values are created, and a smaller time step requires more calculations to be performed and more memory to store the results.

A discussion of the Runge-Kutta-Fehlburg method can be found in many numerical methods texts. C code is also available to perform the routine[13].