A responsibility of the National Highway Traffic Safety Administration (NHTSA) is to write and enforce federal regulations concerning motor vehicle safety. Under the New Car Assessment Program, the safety of a new vehicle is evaluated by analyzing data collected from a series of crash tests. One of the tests performed is the 35-mph frontal barrier impact test, whereby a vehicle is propelled such that it crashes head-on at 35 mph into a rigid barrier that spans the width of the vehicle; the video footage in Figure 1 shows a 2016 Honda CR-V undergoing this type of test.

The test dummies (a driver and typically at least one passenger), the vehicle, and the impact barrier are highly instrumented to measure, for example, acceleration and loads; the corotational components of angular velocity experienced by the dummies’ heads are sometimes measured. These acceleration and angular velocity data may then be used to numerically simulate the motion of the dummies’ heads during the crash.

## Equations of motion

Consider the crash test dummy head shown in Figure 2. We shall parameterize the rotation of the head using a 3-2-1 set of Euler angles , , and that relate the space-fixed basis to the head’s corotational basis , where the basis vectors and are oriented such that they are consistent with the Society of Automotive Engineers’ standard for dummy coordinate systems [1].

For a 3-2-1 set of Euler angles, the basis vectors and are related via the matrix-vector expression

(1)

Consequently, the corotational components of angular velocity, , may be written in terms of the Euler angles , , and and their rates of change as

(2)

Thus, provided crash test data for from the inertial measurement unit installed at the mass center of the dummy’s head, we may determine the head orientation over the duration of the crash by numerically integrating the system (2) of first-order differential equations. It is convenient to express these equations in the matrix-vector form for numerical integration in MATLAB, where the state vector :

(3)

Now, suppose we locate the head’s center of mass using a set of Cartesian coordinates, , and therefore the acceleration of the head’s mass center has the equivalent representations

(4)

where are the corotational components of absolute acceleration obtained from the specific force data recorded by the head’s inertial measurement unit. Having already computed the head’s orientation over time, we next utilize transformation (1) and crash test data for to evaluate the space-fixed components of absolute acceleration, , that are then numerically integrated twice to determine the displacement components of the head’s center of mass:

(5)

## Simulation and animation

The initial state of the head is such that the head is facing directly forward while traveling at 35 mph in that direction: , mph, and . For convenience, but without loss of generality, we assume the motion of the head’s mass center begins at the origin: . Using data collected during a 35-mph frontal barrier impact test of a 2016 Honda CR-V, the simulated absolute motion of the driver dummy’s head is animated in Figure 3.

Visually, the simulated head motion well approximates the observed behavior captured in the high-speed footage shown in Figure 4.

We may also approximate the motion of the head relative to the vehicle using acceleration data collected at, say, the rear left of the vehicle’s floor pan. The simulated relative motion of the dummy’s head is animated in Figure 5.

Based on the dummy’s response captured in the onboard high-speed footage in Figure 6, the simulated relative head motion is a good approximation of the head’s actual motion relative to the vehicle.

## Downloads

The animation in Figure 3 was generated using the following MATLAB code:

## References

- SAE International Surface Vehicle Information Report,
*Sign Convention for Vehicle Crash Testing*, SAE standard J1733, rev. Nov. 2007. The original Dec. 1994 version is freely available to the public and can be found here.