The Traffic Accident Reconstruction Origin -Article-


[Home] [ARnews] [Contents] [Classified] [Advertisers] [Approach Angles] [E-mail Directory] [Feedback] [Organizations] [Reference Library]

The Least Squares Method Applied to Accident Reconstruction

By Bruno Schmidt, PhD

Frequently, when we perform experiments, we can change one of the variables (often called the independent variable) in the experiment, and measure the resulting change (if any) of one of the other variables (called the dependent variable) . For instance, we may:

In both of these cases, we expect a linear relationship between the dependent and the independent variables. A hypothetical Crush vs Impact Speed linear relationship is illustrated below:

 

 

Note that this is not the usual way of displaying crush and speed. The Campbell model typically shows speed on the y-axis and crush on the x-axis. The usual convention is to display the independent variable on the x-axis and the dependent variable on the y-axis. Since for this example, we are talking about changing the impact speed to observe its effect upon crush, the graph does follow the usual convention. In the Campbell discussion of crush, however, he assumes that we have a vehicle whose crush we first measure (in effect, an independent variable) and from which, knowing the crush coefficients, we determine the corresponding barrier speed (the independent variable).

 

Now, if we make measurements of resulting crush from several different controlled speeds, the variations in measurements and results that are inherent in any real-world experiment will produce a "scatter" of results about the theoretical straight line. For example, assume that we have the following results from a series of three rigid barrier crash tests:

 

 

If we plot these values on a graph, we obtain the following result:

 

 

Upon examination of the data points on the graph, we find that they do not follow a perfect line-yet the Campbell model predicts they should. If we try to find the best theoretical fit to the experimental data points, we must ask ourselves, what is the "best" straight line to draw? Two possibilities are shown in the following graph, one for a threshold for crush energy at 7 mph (the solid black line), and one at 4 mph (the dashed red line). The solid black line has a steeper slope-that is, once the threshold is reached, a change in one mph for the impact speed causes more crush than for the dashed red line:

 

 

Constructing a Line


Mathematically, a straight line drawn on an x-y graph can be expressed as

where

m is the slope (rise over run) of the line and

b is called the intercept (the value of y when x is zero).

Fitting a Line to Experimental Data


The approach to use in finding the "best" straight line consistent with the experimental data points which do not lie on a single straight line is to vary the values of m and b until the straight line is the closest to all the data points. How do we define "closest"? In the least-squares method, it is defined as the line where, if we take the values of x that correspond to the measured values of the independent variable (speed in this example), calculate the corresponding values of y from the equation (crush in this example), take the difference between those calculated values of y and the experimentally measured dependent variable values, square those differences and add them up, and then adjust the values of m and b until that sum of squares is a minimum. In effect, then, we will have found the squares of the errors and we select a straight line to make the cumulative value of those squares as small as possible. For instance, in the example that we have been looking at, the m and b values of the solid black line are 0.80 inches/mph and -5.60 inches, while the m and b values of the dashed red line are 0.70 inches/mph and -2.80 inches. If we calculate the sum of the squares of the errors for the first line, we have the following:

 

Independent variable speed

Dependent variable experimental crush

Calculated crush from y = mx + b

(m=0.8, b=-5.6)

Error

(inches)

Error Squared

(inches squared)

10 mph

2.1 inches

2.4 inches

0.3

0.09

25 mph

15.9 inches

14.4 inches

-1.5

2.25

40 mph

24.3 inches

26.4 inches

2.3

5.29

 

 

 

Sum of the squares

7.63

 

If we do the same set of calculations for the second line, we have this result:

 

Independent variable speed

Dependent variable experimental crush

Calculated crush from y = mx + b

(m=0.7, b=-2.8)

Error

(inches)

Error Squared

(inches squared)

10 mph

2.1 inches

4.2 inches

2.1

4.41

25 mph

15.9 inches

14.7 inches

-1.2

1.44

40 mph

24.3 inches

25.2 inches

0.9

0.81

 

 

 

Sum of the squares

6.66

 

This means that the dashed red line (with m = 0.70 inches/mph and b = -2.8 inches) is a "better fit" to the data than the solid black line (6.66 is less than 7.63). The technique that can be used to find the best fit is to start with a value of m and b and then systematically vary one of those two values until the sum of the squares is a minimum, then systematically vary the other parameter until the sum of the squares is a better minimum, then return to the first parameter and vary it to further improve the minimum, and continue back and forth until no better minimum can be obtained. Fortunately, modern spreadsheets will do that automatically for us. The result of this, with the corresponding equation for the best straight line, is shown in the following graph:

 

 

In Microsoft Excel, this is obtained in the Chart menu, where you select Add Trendline. One of the possibilities is a straight line trendline, although a variety of other shapes of lines are possible in Excel. In all cases, however, a least squares fit is used.

Notice that none of the data points lie exactly on the line. Nonetheless, this line still comes closest to all the data points. This is the nature of a least squares fit.

The least squares result for the final spreadsheet trendline is displayed below:

 

Independent variable speed

Dependent variable experimental crush

Calculated crush from y = mx + b

(m=0.74, b=-4.4)

Error

(inches)

Error Squared

(inches squared)

10 mph

2.1 inches

3.0 inches

0.9

0.81

25 mph

15.9 inches

14.1 inches

-1.8

3.24

40 mph

24.3 inches

25.2 inches

0.9

0.81

 

 

 

Sum of the squares

4.86

The best previous Error Squared value obtained was 6.66. Clearly the Error Squared of 4.86 calculated using the Excel trendline is smaller and represents a better fit. From this we decide the line determined by a value of m = .74 and b = -4.4 is the best fit for the data.

Finally, we put the m and b values back into the context of our experiment and see that our barrier collision experiments have produced the model:

which means that entering the value of V in mph into the equation gives the theoretical crush in inches.

Note that this equation can be solved algebraically for the situation where the crush is just zero:

This means that at any speed above 5.9 mph, the crush is greater than zero, so the equation predicts that the threshold for damage is 5.9 mph.

Conclusion

Many natural phenomena follow a linear relationship between variables. Hopefully, this discussion of least squares techniques will be helpful to reconstructionists, who wish to relate theoretical formulas to experimental data, in using a systematic, quantitative approach in fitting data to calculations.


Dr. Bruno Schmidt received his Bachelor's degree in mathematics and physics from Cornell College in 1964. His Doctorate degree in physics followed from Iowa State University in 1969. He was a Physics faculty member at Southwest Missouri State University (1969-1984) and head of the Computer Science Department at Southwest Missouri State University (1984-1992). Currently he is Vice President for Academic Affairs at Southwest Missouri State University.

Dr Schmidt has been practicing in the field of accident reconstruction since 1989. He also teaches the Applied Physics for Accident Reconstruction for Texas A&M Engineering Extension Service. He maintains active memberships in S.A.E, Texas Association of Accident Reconstruction Specialists (TAARS) and National Association of Professional Accident Reconstruction Specialists (NAPARS).

He can be reached at SchmidtBF@aol.com

Copyright ©



[Home] [ARnews] [Contents] [Classified] [Advertisers] [Approach Angles] [E-mail Directory] [Feedback] [Organizations] [Reference Library]