A sloped highway median provides a run-off area for errant vehicles where they can be slowed down and stopped without the danger of being re-directed into the path of other vehicles as may occur with edge barriers. According to Gabler et al.  , cross-median accidents are the most dangerous type of highway crashes so they are usually prevented by placement of a containment barrier at the bottom of the slope to prevent the errant vehicle from crossing the median and running into the opposite carriageway. However, possible redirection of the vehicle colliding with the barrier, back to the same carriageway might cause accidents as serious as the fatal cross median crash.
The research on Highway Median Safety in recent years, with wide application of simulation software, has mainly focused on two themes: one is investigation of double-sloped median terrain’s influence on vehicle behaviour  ; the other is design of barrier, placed in the median, to improve the performance of the median  . In 2008, a systematic method, VDA, Vehicle Dynamics Analyses  was proposed by the National Crash Analysis Center (NCAC), which plotted the trajectory of a vehicle driving into the median in terms of vertical and lateral position. Based on VDA, further research used both Finite Element Analysis (FEA) software and Multi-body Dynamics software to simulate a vehicle driving into the median, from which the guideline of optimal median barrier placement was drafted  . A trade-off was found by Stine et al. between designing a median to prevent vehicular rollover and designing it for preventing median crossovers, after a total run of 111,132 simulations in CarSimTM  . Reliability of applying CarSimTM in such research was validated by Uzunsoy et al.  by comparing the result of simulation with that of physical off-road tests.
The results of these researches indicate that in the study of errant vehicle behaviour in a double sloped median, the profile and geometric dimension of the median play decisive roles, and using simulation software of either FEA or Multi-body Dynamics is effective. The focus of this paper is to determine the most suitable median slope to prevent the vehicle from being re-directed outside the median after collision with the containment barrier. A time-saving method is proposed to simulate the kinematic change of the errant vehicle during collision with the containment barrier while CarSimTM is used to simulate the vehicle’s trajectory before and after the collision.
CarSimTM provides a reliable simulation environment for the dynamics of an errant vehicle as its non-linear Multi-Body Dynamics vehicle model has been shown to reflect the vehicle dynamics property of real vehicles running on the sloped terrain surface of a road median  . However, it cannot simulate the collision with the containment barrier. Therefore additional software has been adopted in simulating the collision phase of the errant vehicle behaviour.
2.1. Vehicle Simulation Structure
Various types of mathematical vehicle collision model have been developed by researchers for two purposes, vehicle collision accident reconstruction, and vehicle post-collision active safety control strategy development     . Brach developed a 3 Degree of Freedom (DOF) collision model based on planar collision mechanics, in which the collision force does not need to be calculated  . A similar collision model, based on momentum expression of Newton’s 2nd laws, has been employed in the software PC-crash which specialises in vehicle collision simulation or reconstruction   . In contrast, another vehicle collision software EDSMAC 4, is based on a different collision model which calculates the collision force by using crash stiffness and acceleration expression of Newton’s 2nd law    . To develop active safety strategy after light impacts, Zhou et al. developed a 4 DOF collision model including roll motion and tyre forces which are not considered by Brach   . The successful development of these mathematical collision models, especially the momentum-based model, demonstrates that vehicle kinematic behaviour right after the collision could be calculated by using these models.
In this study, the kinematic change of the vehicle during the collision with the rigid containment barrier is approximated by using the collision model based on planar impact mechanics. Thus a hybrid simulation platform is used consisting of a multi-body dynamics simulation for the pre- and post-collision phases (using CarSimTM), and planar impact mechanics (using MATLABTM). The general structure of the simulation, shown in Figure 1, consists of three stages, Simulation of Vehicle Motion before Collision, Computation of Collision Model for Vehicle and Rigid Barrier, and Simulation of Vehicle Motion after Collision.
2.2. Collision Model for Vehicle and Rigid Barrier
The collision model based on planar impact mechanics developed by R.M. Brach  was chosen to calculate the post-collision velocities of vehicles since CarSimTM is incapable of simulating the collision of the vehicle with rigid barrier. Such momentum-based collision model has already been applied in the commercial software PC-crash  which specializes in the simulation of vehicle-barrier
Figure 1. General structure of simulation.
and vehicle-vehicle collision  . An alternative would have been Finite Element Analysis (FEA) but it is too time-consuming  for this project which requires at least 36 runs of simulation. An alternative vehicle collision model for the software EDSMAC 4, which computes collision force in order to obtain acceleration and velocity by using Newton’s 2nd law  , was disregarded because the collision force is not of interest.
The following assumptions have been established for the application of the collision model:
1) The plane on which the vehicle is running (motion plane) just before the collision and the motion plane right after the collision are assumed to be the same.
2) The common contact surface (Figure 2) is parallel to the x coordinate, and the right front wheel centre is on that surface.
3) The level of force acting over the common contact surface is significantly higher than other forces, such as aerodynamic force, tire-roadway friction, etc. So impulses of all forces other than contact forces are neglected  .
4) The duration of the contact force impulse is very short, which implies that during contact, accelerations are high to the extent that velocities change suddenly and displacements (changes in position and orientation) are negligible  . Therefore, the final position and orientation in the pre-collision simulation could be assumed to be the initial position and orientation in the post-collision simulation.
5) Though the deformation during the collision is inevitable, the value of mass and moments of inertia about the Z coordinate, and the position of the sprung mass centre relative to the vehicle are assumed to be unchanged.
6) Coefficient of restitution e is assumed to be in the range of 0.1 - 0.3  .
7) Larger absolute values of pre-collision yaw-angle (collision angle), produce more severe inelastic deformation of the vehicle, which leads to a smaller value of e  .
Figure 2. Free-body diagram of vehicle-barrier collision.
8) Impulse ratio μ can be approximated by an equivalent coefficient of friction. Since it cannot be determined experimentally or estimated by means such as analytical modelling of the mechanism of tangential force generation, it is taken as the coefficient of dynamic friction between the common materials of the barrier and car body for approximation   .
Based on Assumption 1, planar impact mechanics including six pre-collision vehicle kinematic states ( ) is applied to the collision model. These represent the pre-collision linear and rotational velocities of the vehicle and barrier in the earth-fixed coordinates.
The first three variables are known from the pre-collision simulation result (see Stage 1 in Figure 1), and the last three variables are assumed to be 0 for the barrier. Also according to that simulation result, other kinematic variables such as vehicle pitch rate and roll rate could be assumed to be negligible. Correspondingly, eight post-collision vehicle motion variables ( ) are to be calculated by using Equations (1)-(11). These variables represent the post-collision linear and rotational velocities of the vehicle and barrier in the earth-fixed coordinates as well as the tangential impulse (Pt) and normal impulse (Pn).
Newton’s laws in the form of impulse and momentum are applied to the bodies as shown in Figure 2. The equations of impulse and momentum for vehicle and barrier in x coordinate and y coordinate are:
According to Newton’s second law, changes in angular momentum are equal to the moments of the impulses on each body. By taking moments of momentum about vehicle mass centre, two more equations about vehicle yaw rates before and after the collision are obtained:
Coefficient of restitution e is defined by Equation (7)  , and the higher the value of e, the more elastic will be the collision  .
Equation (8) is obtained by applying Equation (7) to the vehicle and barrier in Figure 2. The relative normal velocities, and (before and after collision) could be obtained by using the rigid-body kinematic relationship between the contact point and mass centre, as shown in Equation (9) and (10).
Collecting and manipulating Equations (1)-(10), and transforming them into matrix form leads to Equation (11). The values of known variables in Equation (11) are obtained from Vehicle Assembly Screen in CarSimTM: From Assumption 2, equals to the summation of sprung mass and unsprung mass; could be read directly; and are calculated by Pythagorean theorem given the front track width and the horizontal distance between the wheel centre and mass centre; (kg) and (kg∙m2 ) are assigned values of 107; shares the same value as , so do the and , which treats the barrier as a vehicle of the same size as, but far heavier than, the bullet vehicle  . The values of , and are taken as zero, while , and are known from the pre-collision simulation result. equals 0 as shown in Figure 2. With the values of known variables substituted, post-collision vehicle motion variables could be computed from Equation (11).
Impulse ratio μ defined by Equation (12) measures how much tangential impulse is generated during the collision   .
An improper value of μ would produce unrealistic results of post-collision vehicle motion variables contrary to physical laws. Therefore, two criteria are imposed to validate the computation result.
The first criterion is to calculate the kinetic energy of the vehicle before and after the collision by using Equations (13) and (14). If the kinetic energy after the collision becomes higher than that before the collision, then the result of computation by using Equation (6) will be invalid.
The second criterion is to check that the impulse ratio assumed is less than the critical impulse ratio when the vehicle and barrier reach a common tangential velocity―that is Equation (15). The vehicle is assumed to continue to slide over the contact surface throughout the contact duration, so the ratio of μ to the critical impulse must be lower than 1  .
Assembling Equations (1)-(13) into matrix form as shown in Equation (16) and using the definition of impulse ratio, it is possible to calculate . If the ratio of to is larger than 1, then the post-collision motion vector computed by Equation (11) will be invalid.
The values of and are substituted into Equations (17) and (18) and resolved into the longitudinal speed and lateral speed of the vehicle (vehicle-fixed coordinates) after the collision, based on Assumption 4. The values of, and calculated by Equations (11), (17) and (18) and validated by Equations (13), (14) and (16) are input to the simulation of vehicle motion after collision (post-collision vehicle motion).
Equations (1)-(18) are implemented as MATLAB™ M-files for the computation of the collision model.
2.3. Simulation Procedure
CarSimTM is applied to simulate the vehicle motion before collision (pre-collision simulation): the vehicle is running at specified encroachment speed and encroachment angle before passing the median border; when the vehicle just passes the median border, all the control actions exerted on the vehicle are released, resulting in a control-free vehicle whose motion status is only subject to environmental factors such as gravity and terrain; the simulation is terminated instantly the vehicle contacts the rigid barrier. These conditions are implemented in CarSimTM as VS commands. At this point the coordinates of the vehicle in the earth-fixed coordinates (x, y, z) are recovered as well as the vehicle’s dynamic state (velocities, orientation etc.) in the vehicle-fixed coordinates (X, Y, Z). The vehicle-fixed coordinate system is a moving coordinate system which is fixed to the centre of gravity of the vehicle whereas the earth-fixed coordinate system is a coordinate system whose origin is permanently fixed to the earth (see Figure 2). The vehicle’s pre-collision kinematic state is converted into the earth fixed coordinates (V1x, V1y, ω1) and used as input into the collision model to compute the post-collision kinematic state. The resulting post-collision kinematic state (V1x, V1y, Ω1) is then transformed into the vehicle-fixed coordinates (V1X, V1Y, Ω1). These, along with the pre-collision coordinates constitute the initial state of the vehicle for the post-collision simulation in CarSimTM. Again, all controls are released for the post-collision simulation. The post-collision simulation needs to end when the vehicle just runs out of the median or contacts the barrier again. These conditions are again implemented in CarSimTM as VS commands.
3. Results and Discussion
3.1 Simulation Parameters
In this study, the effect of median slope on the effectiveness of the containment barrier has been investigated. Median slopes of 1:4, 1:5, and 1:6  were investigated. In order to manifest the effect of only the median slope, other geometric parameters, e.g., median width and median depth were kept constant. The median width was assigned a value of 75 feet (22.86 m), an average value of median widths applied in rural areas of some states in the USA  . A median depth of 6.25 feet (1.905 m) was assigned, the very value for the median profile of 1:6 median slope when its bottom width equals zero. A general view of the road layout with median and containment barrier is shown in Figure 3.
The encroachment angles used were 5˚, 10˚, and 15˚ and the encroachment speed was chosen to be 37.55 mph (about 60 km/h). These values are representative of the values of typical vehicle encroachment parameters. Encroachment angles smaller than 15˚ cover 75% of vehicle encroachment accident cases  . Encroachment angles larger than 15˚ are more inclined to cause severe collisions, which indicate a higher energy loss  and less possibility that the vehicle runs back to the carriageway. The value of the chosen encroachment speed, 60 km/h, is the mean value of the maximum and minimum of the encroachment
Figure 3. General view of median profile (Median slope 1:4).
speed interval 35.1 mph to 40 mph, which accounts for the largest proportion of the relevant traffic accidents (35%) among other intervals  . Other values of encroachment speed are not considered, as the vehicle at a speed lower than 60 km/h would not travel back to the carriageway if it is contained in the median at an encroachment speed of 60 km/h; and an encroachment speed much higher than 60 km/h would result in an severely inelastic collision or fly over the median barrier  , such phenomena are beyond the scope of this study.
The vehicle encroachment simulation tests were carried out for four classes of vehicles, taken from the CarSimTM database, namely: C-Class Hatchback 2012, D-Class Sedan, E-Class Sedan and Pickup. These represent a wide range of existing vehicle product, such as Audi A3, Ford Focus, Chevrolet Silverado, BMW 3 Series, BMW 5-Series, etc  . All the default configurations of the vehicle models in CarSimTM were kept.
From the result of pre-collision simulation, the collision angle is relatively small, ranging from 15˚ to 25˚. Correspondingly, the vehicle velocity, normal to the contact surface (Figure 2), ranges from 4.31 km/h to 7.04 km/h. For such range of normal speed during the collision, the value of e is between 0.25 and 0.3  , falling in the range of e specified in Assumption 6. Therefore for approximation, and based on Assumption 7, the following discrete function is applied:
As described in Assumption 8, the value of μ could be related to the coefficient of friction when no experimental value is available. For approximation, the value of 0.3 is assumed in this study based on some sliding friction coefficient data  .
3.2. Intermediate Results
The intermediate result refers to the kinematic variables (see Figure 1) of the pre-collision simulation and computation of collision model, some of which are input to the next stage of the method so that the final result of the whole simulation could be calculated stage by stage. Examples of intermediate results derived from the raw data obtained after the computation of the collision model are shown in Table 1 for two of the vehicles:
The first two derived variables are metrics for validating the vehicle-barrier collision. The results show that in all cases, the imposed criteria for validating
Table 1. Data of derived kinematic variables for vehicle (a) Class C HATCHBACK; and (b) Pick up, full size.
the simulation results are satisfied i.e. in terms of the loss in kinetic energy and ratio of the assumed impulse ratio to the critical impulse ratio i.e. ().
The last three variables compare the vehicle kinematic status post-collision with that of the initial condition for the pre-collision simulation. The important observations from these results are:
• As the encroachment angle increases (collision angle increases too), the rate of kinematic energy loss increases, accounted for by the fact that a larger collision angle results in an increased inelastic deformation.
• As the vehicle drives downward the sloped terrain, the yaw angle can increase considerably (up to 290% in case of the 5˚ encroachment angle). However the sloped terrain has little effect on the longitudinal speed.
• The most notable effect is on the roll angle where there is a change of up to 900% for the 1:6 median slope which is quite substantial, whereas for the flat bottomed 1:4 slope, the change is quite modest. The indication is that the flat bottom enables the suspension to level before the collision, resulting in less roll on impact.
3.3. Final Results
The CarSimTM simulation of the post-collision trajectory of the vehicle allows an assessment to be made of the susceptibility of the vehicle to be re-directed outside the median into the path of other vehicles. Figure 4 shows the trajectories of the Class C Hatchback after collision for the three median slopes (1:4, 1:5 and 1:6) at three encroachment angles (5˚, 10˚ and 15˚). The 1:4 median slope results in, by far, the largest rebound, while there is very little difference between the 1:5 and 1:6 median slope results. Nevertheless, in all cases, the vehicle remains well within the median slope boundary.
The post-collision trajectories for all four vehicles tested on the three median slopes at the same three encroachment angles described earlier are shown in Figure 5. A summary of the results is presented in Table 2 which clearly indicates that the 1:6 median slope out-performs the other two median slopes in terms of susceptibility of the vehicle to be re-directed outside the road median. The vehicle motion envelope for the 1:6 median is shown in Figure 5. This indicates that in most instances the vehicle will not deviate far from the median centre line also reducing the susceptibility to roll over  .
Figure 4. Overlapped vehicle trajectory after collision―Class C Hatchback on three median profiles.
Figure 5. Overlapped trajectory of vehicle after collision―all vehicles on three median profiles.
Table 2. Average maximum displacement of vehicle after the collision for three median profiles.
This study investigates the use of a sloped road median with a containment barrier to prevent an errant vehicle from being redirected into the path of other vehicles. A time-saving simulation method has been developed for assessing the effectiveness of different median slopes in containing such errant vehicles within the sloped median. A simple collision model based on planar impact mechanics has been implemented to simulate the collision of the vehicle with the containment barrier, greatly reducing the computational cost associated with the collision phase of the errant vehicle’s trajectory in comparison with other methods such as Finite Element analysis. This model is implemented in MATLABTM and works in conjunction with CarSimTM to simulate the trajectory of the errant vehicle. The total time for simulating one errant vehicle sloped median encroachment including pre- and post-collision simulation in CarSimTM plus collision in MatlabTM is less than 30 seconds.
By comparing vehicle x-y trajectories for median profiles with each other for different encroachment angles, the 1:6 median slope was shown to be the most effective one in containing the errant vehicle. Future study is anticipated to include 4 DOF vehicle collision models, which could predict the post-collision roll motion, in the simulation method, and to validate the simulation method by comparing the simulation results with the result for the same vehicle situation generated by professional vehicular collision software, EDSMAC 4 or PC-crash.
 New Jersey Department of Transportation Bureau of Research and US Department of Transportation Federal Highway Administration. Technical Report, February 2005.
 Stine, J.S., Hamblin, B.C., Brennan, S.N. and Donnell, E.T. (2010) Analyzing the Influence of Median Cross-Section Design on Highway Safety Using Vehicle Dynamics Simulations. Accident Analysis and Prevention, 42, 1769-1777.
 Day, T.D. (1999) An Overview of the EDSMAC4 Collision Simulation Model. SAE Technical Paper 1999-01-0102. SAE Publication SP-1407. Accident Reconstruction: Technology and Animation IX. Proceedings of the 1999 SAE International Congress & Exposition, Detroit, MI, 1-4 March 1999.
 Kuzel, M., Werner, S. and Richards, D. (2006) The Effect of Stiffness Coefficients on Output Variables in EDSMAC4 Simulations. SAE Technical Paper 2006-01-1396.
 Liu, W., et al. (2009) Springback Control of Sheet Metal Forming Based on the Response-Surface Method and Multi-Objective Genetic Algorithm. Materials Science and Engineering: A, 499, 325-328.