The Traffic Accident Reconstruction Origin -ARnews-


Re: Rotating Vehicles Drag Factor

Brian McHenry (mchenry@interpath.com)
Sat, 15 Mar 1997 11:54:59 -0500 (EST)

To approximate the drag factor for a rotating vehicle. I caution against attempting to use the arbitrary % drag factors mentioned in this thread: "50% to 80% of the standard drag factor" or "70-75% of the full drag".
Our recent paper "CRASH-97 - Refinement of the Trajectory Solution Procedure" SAE Paper 97-0949
includes a discussion of approximations for drag factors and the use of a "Runge-Kutta form of differential eqaution for the numerical integration" (also mentioned in this thread).
As we point out in the paper, with the astounding power and performance of common PC computers (Pentiums, etc.) the use of a simulation to approximate the spinout of a vehicle to rest is a relatively easy task.
Any user of the CRASH3/EDCRASH or other CRASH based program has the trajectory solution procedure from SMAC included.
The CRASH3 trajectory simulation procedure is basically the SMAC trajectory solution procedure.

The following is from an Appendix of SAE 97-0949 which contains a discussion of rotating vhicle drag factors, etc.
(The paper is (c) SAE, please visit SAE at
http://www.sae.org
to obtain a copy of the full paper).

APPENDIX 1: DISCUSSION OF SPIN2 from SAE 97-0949

The SPIN2 procedure of the original CRASH program uses as a starting point the relationships developed by Marquard [38]. The Marquard procedure takes into account the fact that the linear and angular (i.e., yaw rotation) displacements of a four-wheeled vehicle subsequent to a collision each occur under conditions of intermittent deceleration when the wheels are free to rotate.
By approximating the linear and angular deceleration rates of a vehicle with either
(1) all wheels freely rotating or
(2) all wheels locked during different phases of spinout motion,
Marquard developed approximate relationships between the total linear and angular displacements during the travel from separation to rest and the corresponding linear and angular velocities of a vehicle at separation from its collision partner, for the two cited cases of rotational resistance.
In the CRASH program [1], the SPIN2 routine was developed to extend the relatively simple Marquard relationships to include the cases of partial braking and/or damage-locked individual wheels.
Evaluations of the resulting, modified relationships by means of trial applications to spinout trajectories generated with SMAC [20] revealed several shortcomings of the initial SPIN2 relationships.

First, a residual linear velocity frequently exists at the end of the rotational (i.e., yawing) motion.
Next, the general shapes of plots of linear and angular velocity vs. time changed substantially as functions of the ratio of linear and angular velocity at separation from the collision.
Finally, the transitions between the different deceleration rates of linear and angular
motions were found to occur gradually rather than abruptly. Slope changes in the plots of linear and angular velocity vs. time were found to generally occur in the form of rounded "corners" in the curves.
To improve the accuracy of approximations of separation velocities, provisions for the introduction of a residual linear velocity at the end of the rotational motion and the development of empirical coefficients, in the form of polynomial functions of the ratio of linear to angular velocity at separation, were incorporated in the SPIN2 analytical relationships of the CRASH program.
Since the separation velocity ratio is initially unknown, a solution procedure was developed whereby several trial values of the ratio, based on an approximate equation, were used to test multiple solutions.
The cited analytical developments, reported in [2], involved only limited efforts which were aimed primarily at demonstrating the feasibility of the CRASH concept. Polynomial functions to generate empirical coefficients were developed, on the basis of 18 single-vehicle SMAC runs with relatively high linear and angular velocities for starting (i.e., separation) conditions. In the more common, real-life accident case, a relatively small rotation (i.e., yawing) velocity may exist at separation. In such a case the initial direction of the velocity vector with respect to the longitudinal axis of the vehicle will obviously affect the sequence and the duration of the linear and angular deceleration rates
of the vehicle.
In consideration of known shortcomings of the SPIN2 aspect of the CRASH program, a subcontract to refine SPIN2 was undertaken in 1979 [24]. A representative sample of actual accident cases was selected from the NCSS [39] files for use in the study. A total of 50 cases were selected and then reconstructed with the SMAC computer program. For each of the SMAC reconstructions, separation information was used to formulate a basis for a refinement of the SPIN2 empirical coefficients. A careful examination of the time-history plots of linear and angular velocities for all of the cases in the sample revealed a significant number
of cases in which the SMAC-predicted behaviour deviated from the analytical assumptions upon which the SPIN2 routine is based. Attempts were undertaken within the research project to discriminate characteristics of separation conditions. Unfortunately, only partial success was achieved in the attempts to accommodate deviations by means of the use of logic and discriminators. As a result, a realistic appraisal of residual scatter in the empirical fits led to the conclusion in [24]:
"To achieve a general improvement in the reliability and accuracy of approximations of the angular and linear velocities at separation, a step-by-step time history form of trajectory solution should be implemented."
Subseqent work which has been performed on investigation and refinement of the SPIN empirical coefficients [30, 40,41] and the corresponding modifications to CRASH is subject to the effects of ‘scatter’. Any proposed refinements of the SPIN empirical coefficients and any reconstruction
techniques which are based on the refinements of the SPIN empirical approach will ultimately fail in some applications to individual case reconstructions due to the possibililty that the particular case being investigated may be characteristic of a "scatter" point.
The research cited in this paper strongly supports the conclusion from 1981 that implementation of a trajectory solution procedure should utilize an iterative time-history simulation.

REFERENCES
1. McHenry, R.R., "The CRASH Program - A simplified Collision Reconstruction Program" Proceedings of the Motor Vehicle Collision Investigation Symposium, Calspan, 1975
2. McHenry, R.R. Lynch, J.P., "User’s Manual for the Crash Computer Program" Calspan Report No. ZQ-5708-V-3, Contract No. DOT-HS-5-01124, Jan 1976
24. McHenry, R.R., McHenry, B.G., "National Crash Severity Study - Quality Contol, Task V: Analysis to Refine Spinout Aspects of CRASH", Calspan Field Services, Inc. ZP-6003-V-4; DOT-HS-6-01442, January 1981
30. Fonda, A.G., "Nonconservation of Momentum During Impact", SAE Paper 95-0355
39. Kahane, C.J, et al, "The National Crash Severity Study", Sixth International Technical Conference on Experimental Safety Vehicles (1976) 495-516
40. Fonda, A.G., Metz, L.D., "Post-Impact Spin, 1968-1993", SAE Paper 93-0653
41. Fonda, A.G., "Energy and Major Diversion in Accident Reconstruction", SAE Paper 96-0888

Brian McHenry
mchenry@interpath.com


NOTE: You are reading in an archived session of ARnews. It is possible that this topic is still being discussed. To see if this topic is still active, or of there were any more recent posts on this topic, check later archives of ARnews.

If there is no current post, and you would like to add to this topic, link to the Current ARnews Discussion and begin a new thread. Be sure that if you are starting a new post that the thread title does not contain the abbreviation RE: Placing RE: at the beginning of a new post will confuse Hypermail and prevent others from answering your post in the future.

For example, to continue this discussion look for a thread titled

Rotating Vehicles Drag Factor

If this thread does not exist in the current archive, you can begin another one by using that title.


here

Speed to Slide to Stop - 2 surfaces

: This option performs the same function as its parent, only 2 surfaces can now be computed. Two surfaces were chosen since the chances of a vehicle cross more than that are fairly slim. The calculations for more than two surfaces is identical to the calculation for two.

Critical Curve

: Over the past couple of years, this calculation has fallen into some disrepute. There are some questions whether the calculation is over- or underestimating the speed.

s = 3.87 SQR(rf)
where


Free Fall

: The main consideration when using this option is to ensure that you get your signs right.

s = 2.74d/SQR(m (d+h))
where

Combined Speeds

: This option takes two speeds and combines them to provide a final speed. Caution should be exercised when it comes to combining speeds which are, in and of themselves, final speeds. A good example is a speed derived from a fall. That's the final speed, regardless of what happens after the fall.

Velocity - Known Distance, Known Time

: This is one of the basic formulae used in physics. It determines the velocity to travel a known distance during a known time.

v = d/t
where

Velocity Gained/Lost

: When the acceleration/deceleration rate is known and the time is known, this calculation will yield the velocity gained or lost.

Velocity At Any Time

: This option will provide the velocity at any time during an acceleration/deceleration when the initial velocity, acceleration/deceleration rate, and time are known.

Velocity At Any Distance

: This option will provide the velocity at any distance during an acceleration/deceleration when the acceleration/deceleration rate and distance are known.

iCRASH - Damage Only

The iCRASH subprogram computes speed changes experienced during a vehicle-to-vehicle or a vehicle-to-fixed object collision. It makes use of the locations and extent of structural crush and is based on the same energy calculations used in the CRASH III program. Users should be aware of CRASH III techniques and limitations before selecting this option. There are two options under this menu selection. A brief discussion of each follows:

New: Tells iCAR to reset all the numerical fields to 0 and clears all buffers. A fresh entry screen will then be presented.

Rerun: Retains all data which was entered during the most recent New or Rerun compilation and allows the user to make changes to any or all data items.

Showit: This option first prompts the user to determine if a printout is required, then generates a screen image showing vehicle damage, PDOF, and D. At present, the screen print utility is limited to use by Hewlett-Packard LaserJet products. A driver for dot-matrix printers is in the works.

Additional Notes: Whenever New is selected, the user will be given the option to create an ASCII file which will record the results from any number of runs. Also, during each run, the user is given the option to alter the vehicle stiffness constants. The decision to offer this option was based on the fact that the default constants, based primarily on the CRASH III model, have some age on them now. More up-to-date information is becoming available. If the user has this information, it should be used. When a larger body of this information becomes available to me, I will alter the internal constants to reflect this change.

Speed Estimates from Damage

Sometimes it is useful to compute a "rule of thumb" estimate of a vehicle's speed based on damage. There is one main menu and one follow-on menu as shown below.

Main Menu

Follow on Menu

Calculate Stiffness Values

Calculate stiffness values based on crash test data. The generated values would be used in the CRASH program.

Time/Distance Calculations

The following selections are used to determine where a particular vehicle was relative to another before a collision and to answer questions as to whether maneuvers by either party could have had a positive effect.

Utilities

Radius of Curvature

: Uses the standard radius of curvature formula. The formula assumes that the curve being measured is a regular (constant) curve. A quick way to check this is to make two measurements in addition to the middle ordinate measurement. The two measurements must occur at like points (i.e., identical distances along the chord). If the measurements are the same, it's a regular curve; otherwise, it is an irregular curve and this calculation is invalid.

R = (cý / 8m) + (m / 2) where R = radius c = chord m = middle ordinate

Convert FPS to MPH

: Feet-per-second to miles-per-hour. mph = fps x 0.6818

Convert MPH to FPS

: Miles-per-hour to feet-per-second. fps = mph x 1.467

Conversion Formulae

: This option serves as a reference source. There are many conversions. Most of these are simple and not worth valuable computing time when a simple calculator will do the job. This option just provides a listing of input to get a desired output. No more, no less. Compute Drag Factor: For a vehicle with all wheels locked, the drag factor is the same as the coefficient of friction. The following three selections all determine acceleration and/or deceleration factors. They vary only what it known prior to the calculation.

Compute AD Factor

: Speed and Distance Known adf = sý / 30d
where

Compute AD Factor

: Velocity and Time Known adf = v / 32.2t
where

Compute AD Factor

: Distance and Time Known
adf = d / 16.1t^2
where

Compute Acceleration Factor

: Vi, Vf, and Time Known: Computes an acceleration factor when the initial velocity, final velocity, and time are known.
a = (vf - vo) / (32.2t)
where

Compute Deceleration Factor

: Vi, Vf, and Time Known: Computes a deceleration factor when the initial velocity, final velocity, and time are known.
a = (vo - vf) / (32.2t)
where

Compute AD Rate

: AD Factor Known: Once the factor is known, the rate is computed by multiplying the factor time the gravity (32.2 feet per second per second).

DOS Operations

iCAR provides the user with two important DOS-related operations. Both are important and necessary. Being able to escape into the DOS shell without leaving a program provides significant flexibility and a result recording operation which runs in the background allows the user to concentrate on the difficult task at hand without having to worry about the details of file manipulation.

Exit to DOS

: Temporarily suspends iCAR, clears the screen, and displays the DOS prompt, from which you can run other programs or DOS commands. You must remember, however, that iCAR is still resident, so your computer will not have as much memory as it would normally. To return from the shell, simply type EXIT at the DOS prompt.

Create an ASCII Session Record

: Creates an ASCII text file using a name of your choosing, activates an internal boolean variable, and will faithfully record calculations results for the duration of the session. Please remember that the file name chosen must follow normal DOS protocol, and that any identical file name will be overwritten. Choose a file name that means something to you (i.e., case0914.jp).

Appendix A - Vehicle Parameters

The appropriateness of the set of eight frontal stiffness coefficients used by the CRASH3 program were examined in an SAE paper entitled "A Comparison Between NHTSA Crash Test Data and CRASH3 Frontal Stiffness Coefficients". The authors (Messrs. Strother, Woolley, and James) generated a new set of stiffness categories which are shown in the following table.

Vehicle Type A(lbf/in) B(lbf/in^2) A(lbf/in) B(lbf/in^2)
Stiffness CatagoryCRASH3
Users Guide
CRASH 3
Users Guide
Strother,
et al
Strother,
et al
1) Subcompact302.047.0 237.958.9
2) Compact259.043.0240.060.0
3) Intermediate317.056.0247.558.95
4) Full-size356.034.0236.751.5
5&6 Largest325.037.0247.257.9
7) Vans383.0126.0349.799.8
MPV's383.0126.0350.9100.5
8) Pickups 480.050.0425.6 72.5
9) Front Wheel Drive373.038.0240.458.2

Appendix B - Vehicle Size Categories

          SIZE           WHEELBASE
          ----           -------------
          1               80.9 -  94.8
          2               94.8 - 101.6
          3              101.6 - 110.4
          4              110.4 - 117.5
          5              117.5 - 123.2
          6              123.2 - 150.0
          7              109.0 - 130.0 VANS
          8              PICKUPS [Select 1 to 6 based on wheelbase]
          9              JEEPS [Select 1 to 6 based on wheelbase]
          11             IMMOVABLE BARRIER

Appendix C - Vehicle Measurement Terms

James Perry is a Traffic Safety Investigator with Dynamic Science in Anaheim, California.
He holds an MS in Information Systems from Nova Southeastern University and a BGS in Psychology/Police Sciences from the Univerity of Nebraska-Omaha.
He can be reached at:
299 West Cerritos Avenue, Anaheim, CA 92805
jperry1@ix.netcom.com


Copyright ©



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