Injury tolerances for oblique impact helmet testing

Similar documents
Intracranial Response in Helmet Oblique Impacts

S Kleiven. doi: /ijcr

The UCD community has made this article openly available. Please share how this access benefits you. Your story matters!

Validating Critical Limits of the Universal Brain Injury Criterion

Importance of Windscreen Modelling Approach for Head Injury Prediction. Victor S. Alvarez, Svein Kleiven

INFLUENCE OF DIRECTION AND DURATION OF IMPACTS TO THE HUMAN HEAD EVALUATED USING THE FINITE ELEMENT METHOD

ROLE OF ROTATIONAL ACCELERATION AND DECELERATION PULSES ON BRAIN STRAINS IN LATERAL IMPACT

Influence of Impact Velocity and Angle in a Detailed Reconstruction of a Bicycle Accident

THE IMPORTANCE OF ROTATIONAL KINEMATICS IN PEDESTRIAN HEAD TO WINDSHIELD IMPACTS.

Design and Finite Element Analysis of Helmeted Head Form in Impact

AN INVESTIGATION OF BRAIN INJURY RISK IN VEHICLE CRASHES (SECOND REPORT)

Injury Biomechanics Workshop at IRCOBI Asia 2018

Simulation and Test Validation of Windscreen Subject to Pedestrian Head Impact

On the Role of a Nonlinear Stress-Strain Relation in Brain Trauma

DEVELOPMENT OF AN IMPROVED DUMMY HEAD FOR USE IN HELMET CERTIFICATION TESTS

NONLINEAR STATIC SIMULATION OF AUTOMOTIVE BUMPER OF A PASSENGER CAR IN LOW SPEED IMPACT CRASH SIMULATION

Material property determination of the lining layers of a versatile helmet

Impact Analysis of Frontal Car Bumper using Long Fibre Reinforced Thermoplastics

Investigating the dynamic response of a punch to human head using finite element analysis

Computational Simulation and Visualization of Traumatic Brain Injuries

MODELLING AND ANALYSIS OF SOCCER HEADING AND PROTECTIVE HEADGEAR TO UNDERSTAND AND PREVENT MILD TRAUMATIC BRAIN INJURY

MECHANICAL PROPERTIES OF POLYTETRAFLOUROETHYLENE ELASTOMER MEMBRANE FOR DYNAMIC CELL CULTURE TESTING ABSTRACT INTRODUCTION

Development of Rotational Brain Injury Criterion (BRIC) Human Injury Research Division

ANALYTICAL PENDULUM METHOD USED TO PREDICT THE ROLLOVER BEHAVIOR OF A BODY STRUCTURE

METHODOLOGY FOR ESTIMATING THORACIC IMPACT RESPONSE IN FRONTAL CRASH TESTS

Finite Element Models for European Testing: Side Impact Barrier to WG13 Pedestrian Impactors to WG17

IMPACT TESTT ON MOTORCYCLE HELMET FOR DIFFERENT ANGLES USING FINETE ELEMENT ANALYSIS

Crash Simulation and Analysis of a Car Body Using ANSYS LS- DYNA

Development and Validation of the FAT Finite Element Model for the Side Impact Dummy EUROSID-1

MECHANICAL CHARACTERIZATION OF BRAIN TISSUE

MODELLING THE SHAKEN BABY SYNDROME

EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS

AN LS-DYNA USER DEFINED MATERIAL MODEL FOR LOOSELY WOVEN FABRIC WITH NON-ORTHOGONAL VARYING WEFT AND WARP ANGLE

Stapp Car Crash Journal, Vol. 46 (November 2002), pp. Copyright 2002 The Stapp Association

Finite-element analysis of bicycle helmet oblique impacts

Finite Element Simulation of Bar-Plate Friction Welded Joints Steel Product Subjected to Impact Loading

ME 2570 MECHANICS OF MATERIALS

PARAMETRIC STUDY OF THE STRUCTURAL CAPACITY OF REINFORCED CONCRETE CONTAINMENT SUBJECTED TO LARGE COMMERICAL AIRCRFAT IMPACT

NUMERICAL SIMULATION OF DAMAGE IN THERMOPLASTIC COMPOSITE MATERIALS

EXAMINATION OF THE IMPACT OF HELMETS ON THE LEVEL OF TRANSFERRED LOADS TO THE HEAD UNDER BALLISTIC AND BLAST LOADS

CONSTITUTIVE MODELING AND OPTIMAL DESIGN OF POLYMERIC FOAMS FOR CRASHWORTHINESS

NUMERICAL INVESTIGATION OF A THREE-DIMENSIONAL DISC-PAD MODEL WITH AND WITHOUT THERMAL EFFECTS

Journal of Biomechanics, 47 (5):

Low-level Fall Simulation: Toddler Falling From Chair

Reducing effect of softball-to-head impact by incorporating slip-surface in helmet

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

ME 176 Final Exam, Fall 1997

Physical Science and Engineering. Course Information. Course Number: ME 100

The design of a formula student front impact attenuator

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.

ANALYSIS AND SIMULATION OF AN AIRPLANE SEAT DURING VERTICAL IMPACTS

Natural Frequency Analysis of the Human Head in Normal models

Compact energy absorbing cellular structure

Crash helmet testing and design specifications

Anisotropic modeling of short fibers reinforced thermoplastics materials with LS-DYNA

Chapter 5. Vibration Analysis. Workbench - Mechanical Introduction ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.

Calibration and Experimental Validation of LS-DYNA Composite Material Models by Multi Objective Optimization Techniques

A prestudy of the potential of using finite element analysis for understanding horse accidents

International Journal of Scientific and Research Publications, Volume 5, Issue 8, August ISSN

N = Shear stress / Shear strain

Structural Integrity Assessment of a Rupture Disc Housing with Explicit FE- Simulation

IMPACT PROPERTY AND POST IMPACT VIBRATION BETWEEN TWO IDENTICAL SPHERES

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

Stress Analysis and Validation of Superstructure of 15-meter Long Bus under Normal Operation

An orthotropic damage model for crash simulation of composites

Solid-fluid interaction in force and energy transmission in shaken baby syndrome

arxiv: v1 [physics.med-ph] 19 Apr 2009

Finite Element Analyses on Dynamic Behavior of Railway Bridge Due To High Speed Train

Impact loading of ductile rectangular plates

Numerical simulation the bottom structures. grounding test by LS-DYNA

A ROLLOVER TEST OF BUS BODY SECTIONS USING ANSYS

Shock enhancement of cellular structures under impact loading: Part II Analysis

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure

ROTATIONAL BRAIN INJURY TOLERANCE CRITERIA AS A FUNCTION OF VEHICLE CRASH PARAMETERS

Quasi-static simulation approaches on rollover impact of a bus structure

DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION

The Simulation of Dropped Objects on the Offshore Structure Liping SUN 1,a, Gang MA 1,b, Chunyong NIE 2,c, Zihan WANG 1,d

Elastic Response of bones at varying strain rates

STRUCTURAL ANALYSIS OF THE LIFTING DEVICE DETECTOR SUPPORTS FOR THE LHCb VERTEX LOCATOR (VELO)

PDDC 1 st Semester Civil Engineering Department Assignments of Mechanics of Solids [ ] Introduction, Fundamentals of Statics

QUASI-ANALYTIC ACCELERATION INJURY RISK FUNCTIONS: APPLICATION TO CAR OCCUPANT RISK IN FRONTAL COLLISIONS

NUMERICAL SIMULATION OF BLAST RESISTANT STEEL PLATE STRENGTHENED WITH COMPOSITE

Constitutive models: Incremental plasticity Drücker s postulate

Two Tier projects for students in ME 160 class

A numerical parametric study on the plate inclination angle of steel plate shear walls

Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method

Modelling of bird strike on the engine fan blades using FE-SPH

Ultimate shear strength of FPSO stiffened panels after supply vessel collision

Evaluation of human thorax FE model in various impact scenarios

EXPERIMENTAL VERIFICATION OF THE USE OF DUMMIES FOR TESTS IN EMPIRICAL RESEARCH

Comparison of crash tests and simulations for various vehicle restraint systems

STUDY ON THE RESPONSE OF ELASTOMERIC BEARINGS WITH 3D NUMERICAL SIMULATIONS AND EXPERIMENTAL VALIDATION

COMPARISON OF EXPERIMENTAL RESULTS WITH FEM ONES OF RECTANGULAR CFRP TUBES FOR FRONT SIDE MEMBERS OF AUTOMOBILES

SIMPLIFIED MODELING OF THIN-WALLED TUBES WITH OCTAGONAL CROSS SECTION AXIAL CRUSHING. Authors and Correspondance: Abstract:

ANSYS Mechanical Basic Structural Nonlinearities

Some Aspects Of Dynamic Buckling of Plates Under In Plane Pulse Loading

MAAE 2202 A. Come to the PASS workshop with your mock exam complete. During the workshop you can work with other students to review your work.

Mathematical Modelling of Car Crash Test

The Determination of Vehicle Speeds from Delta-V in Two Vehicle Planar Collisions

DESIGN OF A FREE MOTION HEADFORM IMPACTOR SUPPORT STRUCTURE FOR SCHOOL BUS AND MOTOR COACH TESTING. Thesis

Transcription:

1 2 3 4 5 Injury tolerances for oblique impact helmet testing M Aare, S Kleiven and P Halldin Division of Neuronic Engineering, CTV Centre for Technology within Healthcare, Royal Institute of Technology (KTH) and Karolinska Institute, Stockholm, Sweden. 6 7 8 9 11 12 13 14 15 16 17 18 19 21 22 Abstract: The most frequently sustained severe injuries in motorcycle crashes are injuries to the head, and many of these are caused by rotational force. Rotational force is most commonly the result of oblique impacts to the head. Good testing methods for evaluating the effects of such impacts are currently lacking. There is also a need for improving our understanding of the effects of oblique impacts on the human head. Helmet standards currently in use today do not measure rotational effects in test dummy heads. However rotational force to the head results in large shear strains arising in the brain, which has been proposed as a cause of traumatic brain injuries like diffuse axonal injuries (DAI). This paper investigates a number of well-defined impacts, simulated using a detailed finite element (FE) model of the human head, an FE model of the Hybrid III dummy head and an FE model of a helmet. The same simulations were performed on both the FE human head model and the FE Hybrid III head model, both fitted with helmets. Simulations on both these heads were performed to describe the relationship between load levels in the FE Hybrid III head model and strains in the brain tissue in the FE human head model. In this study, the change in rotational velocity and the head injury criterion (HIC) value were chosen as appropriate measurements. It was concluded that both rotational and translational effects are important when predicting the strain levels in the human brain. Keywords: Injury tolerances, Oblique impacts, Helmet, Head injuries and FE (Finite Element) 23 24 25 26 27 28 29 3 31 32 33 34 35 36 37 38 39 INTRODUCTION Motorcycle riders are very exposed in traffic accidents since their level of protection is limited. One way to provide better safety for motorcycle riders is to improve the helmets. The most frequently sustained severe injuries in motorcycle crashes are namely injuries to the head [1]. Many of these injuries are caused by rotational forces [2] that are most commonly generated as a result of oblique impacts found in motorcycle crashes [3]. The most frequent impact is an impact close to the visor attachment points with an average impact thereof speed of 44 km/h and an average angle of 28 to the impacted object [3]. Most safety helmets are designed to meet the requirements established by standardised tests (e.g. ECE22.5 [4], FMVSS218 [5], and BS6658 [6]). In these tests, the helmet is typically dropped vertically onto a flat or curved rigid surface that is set up for a tangential impact Corresponding Author: M Aare, S Kleiven and P Halldin Marinens v. 3, 1364 Haninge, Sweden Tel: +46-7495362 Email: aare@kth.se to the helmet surface. During the drop tests, the translational acceleration of a head form is measured. The British Standard BS6658 includes an oblique impact, but there is no measurement of the rotational effects, which makes it difficult to correlate these effects to injuries in the human brain. The oblique impact tests in the BS 6658 ensure: that projecting visor mounts and other projections shear off easily when there is an impact with a series of parallel bars; and that the tangential force on the helmet shell, when it impacts with a rough flat surface, is no larger than that of typical shell materials used in 1985 (the year of introduction of the test). The ECE22.5 standard also includes a test that can be seen as an oblique impact test. However, this is a test for projection and surface friction, and does not include rotational measurements in the dummy head. Today s standards have resulted in helmets with a good degree of protection for radial impacts while their protection for oblique impacts remains unknown. In an oblique impact to a helmeted head, the rotational energy can be absorbed by (a) frictional energy between 4 41 42 43 44 45 46 47 48 49 5 51 52 53 54 55 56 57 58 59 6 61 62 63 Woodhead Publishing Ltd 268 1 IJCrash 4 Vol. 9 No. pp.

M Aare, S Kleiven and P Halldin 64 65 66 67 68 69 7 71 72 73 74 75 76 77 78 79 8 81 82 83 84 85 86 87 88 89 9 91 92 93 94 95 96 97 98 99 1 2 3 4 5 6 7 the helmet shell and the impacting surface, (b) internal energy in the liner and comfort foam due to shear deformation, and (c) frictional energy between the liner and the head. The standard test procedures are based on vary vague reasons. More realistic would be to test helmets for impacts similar to the most commonly observed impacts in real life motorcycle accidents. Aare and Halldin (3) [7] proposed a new method to test helmets for oblique impacts, but they did not propose any injury tolerances for such a test. The injury tolerances and criteria used when testing oblique impacts to helmets should predict injuries frequently seen in traffic accidents. The most frequently sustained brain injuries in motor vehicle accidents that result in either fatality or the need for long-term rehabilitation are subdural hematomas (SDH) and diffuse axonal injuries (DAI) [2]. The main causes of SDH are ruptured arteries or bridging veins. DAI is caused by the tearing of neuronal axons in the brain tissue. It has been proposed that there are correlations between these injuries and rotational effects to the head [2]; [8]. In these studies injury thresholds for purely angular motions was proposed. If an angular motion is combined with a translational motion, these thresholds probably have to be decreased, as shown in the study by DiMasi et al. (1995) [9] and Ueno and Melvin (1995) []. To establish tolerances or criteria for an oblique impact to the helmeted head it is possible to use a detailed finite element (FE) model of the human head [11]. This method should be seen as a complement to biological experiments on human cadavers. It has been suggested that there is a correlation between the strain in the brain and DAI [12]; [13]. If the FE model of the human head is well correlated with relevant experiments, then the strain computed in the model can be compared with accepted tissue-based injury thresholds. Bain and Meaney () [12] proposed a threshold of % strain in the brain tissue for the onset of the malfunction of the nerves in the brain, which could be seen as a first stage of DAI. Currently there is a lack of good testing methods for evaluating the effects of oblique impacts. There is also a need for improved understanding of the effects on the human head subjected to oblique impacts. In current helmet standards, no rotational effects are measured in the head form, partly because there are no accepted global injury thresholds for a combination of rotations and translations. The objective of this study was therefore to study if, and how, rotational and translational parameters influence the strain levels in the brain for three well-defined and commonly observed [3] oblique impacts to helmets. The aim of this study was to propose new injury tolerances for a specific set of oblique helmet impacts. Further it was hypothesized that it is possible to predict the strain in brain tissue using the peak change of rotational velocity and the HIC value. METHODS This is a numerical study using the non-linear and dynamic finite element (FE) code LS-DYNA [14]. In this study, FE models of (1) the human head [15], (2) a Hybrid III dummy head [16], and (3) an experimental helmet were used. Simulations where performed on both the human head and on the Hybrid III dummy head, wearing the helmet. Injury tolerances The maximal principal strain in the brain tissue was chosen as a predictor of injuries, as it has been shown to correlate with DAI [12]. A strain of % was shown to be critical to the brain tissue. As the strains in the brain tissue are proportional to the HIC value for pure translations, and also proportional to the change in rotational velocity for pure rotations of short impact durations [11], it is suggested that the output data is fitted to the following formula: ε = k 1 ω + k 2 HIC [1] where ε is the maximum strain in the brain tissue, ω is the peak resultant change in rotational velocity, HIC is the head injury criteria [17] and k 1 and k 2 are constants. The strain (ε) is taken from the FE human head model, whereas ω and HIC are taken from FE Hybrid III head model. FE model of the human head The FE model of the human head (Figure 1) developed at the Royal Institute of Technology (KTH) Stockholm, Sweden, is based on data from the Visual Human Database 8 9 1 111 112 113 114 115 116 117 118 119 1 121 122 123 124 125 126 127 128 129 13 131 132 133 134 135 136 137 138 139 14 141 142 143 144 145 (a) (b) (c) Figure 1 Finite element mesh of (a) the human head, (b) falx and tentorium including transverse and superior sagittal sinuses with bridging veins, and (c) cranium. IJCrash 4 Vol. 9 No. 2 268 Woodhead Publishing Ltd

Injury tolerances for oblique impact helmet testing 146 147 148 149 15 151 152 153 154 155 156 157 158 159 16 161 162 163 164 165 166 167 168 169 17 171 172 173 174 175 176 177 178 179 18 181 182 183 184 185 [18]. The FE model is anatomically detailed and includes the scalp, the skull, the brain, the meninges, the cerebrospinal fluid (CSF) and eleven pairs of parasagittal bridging veins [15], [19], [] and [21]. FE model of the Hybrid III dummy head The FE model of the Hybrid III dummy head is very similar to the real Hybrid III dummy, with respect to geometry, mass, inertia, and material properties [16]. The model used in this study does not include the neck. This is a simplification used in all standard helmet tests. The solid part of the head (the aluminium) was modelled as rigid, whereas the rubber skin was modelled using a viscoelastic material model [16]. The Hybrid III head model used in this study was meshed using 2311 elements (two elements through the thickness in the solid part making a total of 1288 elements and two elements through the thickness in the rubber skin making a total of 23 elements). FE model of the helmet The FE model of the helmet includes a linearly-elastic material model of the ABS thermoplastic shell, validated against quasi-static tensile tests (thickness: 4.7 mm, Young s modulus: 1.64 GPa, density: 1161 kg/m 3 and Poisson s ratio:.45). The energy absorbing liner material consists of expanded polystyrene (EPS) with a density of 4 kg/ m 3 modelled using a crushable foam material [14]. This foam material model requires a stress-strain curve. Therefore hydrostatic compressive tests were performed on expanded polystyrene (Figure 2). The data from this test were implemented in the material model. Poisson s ratio was estimated to zero and Young s modulus implemented in the material model was estimated to 8 MPa. Young s modulus implemented in this material model represents the stiffness during unloading. The stiffness during loading is taken from the implemented load-curve. The tensile stress cut-off was estimated at 1 MPa and the damping coefficient estimated at.5 since the material was assumed not to be strain-rate dependent [22]. The helmet model was not fitted with a chinstrap, as it was concluded when studying high-speed movies from Stress (MPa) 1.4 1.2 1..8.6.4.2 Hydrostatic test...1.2.3.4.5.6.7.8.9 Volumetric strain (1-V/Vinitial) Figure 2 Stress-strain curve from the hydrostatic test. experiments performed by Aare M. and Halldin P [7] that the inclusion of a chin strap made no significant difference to well-fitting helmets during short duration impacts. The contact definition between (a) the FE model of the human head and the helmet and (b) the FE model of the Hybrid III dummy head and the helmet was the surface-to-surface interference [14]. This means that when the head is larger than the space inside the helmet, an initial pressure or contact force arises. When the head is smaller than the space inside the helmet, there is no initial pressure. However, the inside surface of the helmet fitted with the Hybrid III head model was shaped to fit the geometry of the head perfectly. The inside surface of the helmet fitted with the human head model was shaped to fit the geometry of that head perfectly. This means that there was no initial gap between the heads and the helmets anywhere throughout the contact surfaces. This also means that the geometry of the inside surface if the helmets were slightly different. However, the thickness of the liner in the two different helmets was similar. FE model of the helmet and the Hybrid III head combined The complete model of the helmet and the Hybrid III dummy was validated against both radial and oblique impact tests on helmets (Figure 3). Some of these experimental data was presented by Aare and Halldin (3) [7]. Numerical simulations Three different impacts were tested, where impact 1 is to the top of the head inducing sagital plane rotation, Impact 2 is lateral inducing axial rotation, and Impact 3 is lateral inducing coronal plane rotation (Figure 3). For Impacts 1 and 3, three different impact velocities were used (5, 7 and 9 m/s), and for Impact 2 an additional impact velocity of 3 m/s was also used in the testing. The reason for simulating impact 2 at 3 m/s was that all the other impact velocities caused strains in the brain tissue larger than %. The impact angles were induced by altering the speed of the head and the speed of the impactor (Table 1). For all three impacts and impact velocities, four different impact angles between the head and the impactor were tested (,, and ). These specific scenarios were chosen to cover the range of the most commonly observed impacts in real life motorcycle accidents [3]. Data from the helmeted Hybrid III head was correlated to load levels in the human brain, by performing simulations on both the FE human head and the FE Hybrid III head. Comparisons were made between the strains in the brain tissue of the FE human head and the change of rotational velocity and the HIC-value in the FE Hybrid III head for identical impacts. These comparisons were done to find the relationship between the global kinematics of a Hybrid III head and the strains in the human brain. 186 187 188 189 19 191 192 193 194 195 196 197 198 199 1 2 3 4 5 6 7 8 9 2 211 212 213 214 215 216 217 218 219 2 221 222 223 224 225 226 227 228 229 23 231 232 233 234 235 236 237 238 239 Woodhead Publishing Ltd 268 3 IJCrash 4 Vol. 9 No.

M Aare, S Kleiven and P Halldin Impact 1 Impact 2 Impact 3 Figure 3 The three different impacts in the numerical simulations. Table 1 Vertical and horizontal velocities of the different impact speeds and angles used in the numerical simulations. Impact velocity: 3 m/s Impact velocity: 5 m/s Impact velocity: 7 m/s Impact velocity: 9 m/s Vertical Horizontal Vertical Horizontal Vertical Horizontal Vertical Horizontal Angle of impact speed of the impact speed of the impact speed of the impact speed of the impact speed (m/s) plate (m/s) speed (m/s) plate (m/s) speed (m/s) plate (m/s) speed (m/s) plate (m/s) 1.5 2.6 2.5 4.33 3.5 6.6 4.5 7.79 2.12 2.12 3.54 3.54 4.95 4.95 6.36 6.36 2.6 1.5 4.33 2.5 6.6 3.5 7.79 4.5 3.. 5.. 7.. 9.. 24 241 242 243 244 245 246 247 248 249 25 251 252 253 254 255 256 257 258 259 26 261 262 263 264 265 266 267 268 269 27 RESULTS Figure 4 shows a comparison between experiments and the numerical simulations from impact 1, velocity 7 m/s and impact angle. The results from the numerical simulations are presented in Table 2. Comparing the results from the various different impacts, it is quite clear that the rotational effects have a major influence on the strain levels in the human brain (Figures 5, 6 and 7). Each point in these figures represents one simulated impact. These points where plotted to give the reader a better picture of which impact scenarios that are critical. The maximum strain in the brain was usually found in the white brain tissue, though in different elements. The constants k 1 and k 2 in Equation 1 were computed for Impacts 1, 2 and 3 respectively, using least square regression analyses of Equation 1 (Table 3). The regression coefficient is a figure that describes how well the data correlates to the equation. A regression coefficient of one represents a perfect match and a regression coefficient of zero represents no correlation at all. Different isostrain curves can be plotted using the computed constants k 1 and k 2 in Equation 1 (Figures 5, 6 and 7). DISCUSSION This study has presented a method for computing testspecific injury tolerances using an FE model of the human head. The two regression parameters, HIC and change in rotational velocity, vary in their influence on the head response depending on the impact location and orientation of the helmeted head (Figures 5, 6 and 7). In Impact 1, the change in rotational velocity exhibits a fairly sharp cut-off point for determining safe and unsafe conditions, whereas in Impact 2 and 3, the HIC value (or at least a translational component) also needs to be considered when predicting injury. It is therefore difficult to use the one injury tolerance indicator for all types of impacts. It is suggested that one specific set of constants be used for each impact direction and speed. Regression analysis was used in this study to fit the data to equation 1. The regression coefficient is a figure that describes how well the data correlates to the equation. A regression coefficient of one represents a perfect match and a regression coefficient of zero represents no correlation at all. The regression analysis showed good correlations for Impacts 1 and 2 (Table 2). For this impact, a formula of higher order would probably fit the data better, especially for very low changes in rotational velocity. It is likely that the curves plotted in figures 5, 6 and 7 are more accurate in the central parts of the diagrams close to the load levels that are used in the simulations. The simulated impacts in this study were not only chosen to cover the most commonly observed impacts in motorcycle accidents, but also to cover rotations around all three Cartesian axes. Impact 2 was according to statistics, the most common impact in real-life accidents. The head is vulnerable to accelerations in the lateral direction [23], [24] and therefore the HIC values are more critical for Impacts 2 and 3 than for Impact 1. Even when the head was dropped vertically (), rotations were induced in all the impacts and in particular in Impacts 2 and 3. These rotations occur because the impact point is not situated directly under the centre of gravity. For rotational impulses of short duration, the change 271 272 273 274 275 276 277 278 279 28 281 282 283 284 285 286 287 288 289 29 291 292 293 294 295 296 297 298 299 3 31 32 33 34 IJCrash 4 Vol. 9 No. 4 268 Woodhead Publishing Ltd

Injury tolerances for oblique impact helmet testing Rotational velocity (rad/s) 35 3 25 15 5 Simulation Experiment Rotational acceleration (krad/s 2 ) 12 8 6 4 2 2 4 Simulation Experiment 5 15 5 15 Time (ms) (a) 6 8 Time (ms) (b) Translational acceleration (g) 18 16 14 1 8 6 4 Simulation Experiment 5 15 Time (ms) (c) Figure 4 Comparison of simulations and experiments from impact 1 (Figure 3), velocity 7 m/s and impact angle, where (a) is the change of rotational velocity, (b) is rotational acceleration, and (c) is translational acceleration. 35 36 37 38 39 3 311 312 313 314 315 316 317 318 319 3 321 322 323 324 325 in angular velocity has been shown to correspond better than all other parameters with the intracranial strains found in the FE model [11]. This is in agreement with Holbourn s hypothesis [13]. For translational impulses on the other hand, the HIC [17] and the HIP [25] have shown the best correlation with the strain levels found in the model [11]. In Table 2 all data needed to plot Figures 5, 6 and 7 are displayed. Additionally, the rotational accelerations are displayed in Table 2, since these parameters are used in some other injury criteria [8] and might therefore be interesting to some readers. The translational acceleration is also displayed in Table 2, since this parameter is correlated to the HIC-value and might therefore also be of interest. An FE model of the human head was used to compute the maximal principal strain in the brain, and thereby analyse the risk for DAI. There are also other tissue-level measurements that may be used as injury predictors, such as strain rate, the product of strain and strain rate, von- Mises stress, and strain energy. In this study, the strain has been chosen to analyse the risk for DAI, as this measurement has been experimentally verified previously [12], [26]. An FE model needs to be well correlated to relevant experimental studies on the human brain. Another important issue in modelling of the human head is the selection of appropriate material properties for various intracranial structures. The choice of shear properties for the brain tissue is difficult as the range of published values varies several orders of magnitude [21]; [27]. However, FE model used in this study has been carefully validated and shown to have a good correlation with experiments found in the literature including rotational injuries correlated to strain in the brain as well as local brain motion experiments [28], [29]. Another limitation with the test method presented here (and indeed with most helmet testing methods) is that only the head is used. A more realistic simulation would probably involve using the whole body, or at least the head, neck and torso. Involving for example the neck would change the dynamics of the impact, as the boundary conditions for the head would change. Ruan et. al. (1991) [3] showed that a single hinge coupling between the head and the neck had a limited effect on the intra-cranial 326 327 328 329 33 331 332 333 334 335 336 337 338 339 34 341 342 343 344 345 346 Woodhead Publishing Ltd 268 5 IJCrash 4 Vol. 9 No.

M Aare, S Kleiven and P Halldin Table 2 Results from the numerical simulations. FE model of the FE model of the Hybrid III head human head Peak Peak resultant resultant rotational Maximum Impact Impact Impact translational acc strain in brain area* velocity angle ω (rad/s) HIC-value acc (g) (krad/s 2 ) tissue 1 5 3 26.6 3 2 9.1.16 1 5 45 21.7 651 132 8.5.13 1 5 6 12.2 974 165 4.6.8 1 5 9 6.8 1183 191 2.6. 1 7 3 35.6 652 145 11.8.23 1 7 45 26.5 1363 188 7.8.19 1 7 6 14.6 1939 217 4.9.11 1 7 9 9.5 2279 243 3.6.12 1 9 3 45. 83 177 14.2.3 1 9 45 34.4 26 2 11.6.24 1 9 6 21.2 3178 264 9.2.14 1 9 9 12.3 3791 294 4.8.12 2 3 3.8 5 52 6.3.13 2 3 45 19.5 95 65 6.5.13 2 3 6 16.9 139 75 5.9.14 2 3 9 12.9 183 8 4.7.12 2 5 3 32.9 134 73 9.2.22 2 5 45 31.3 262 97..23 2 5 6 27.4 394 115 9.4.23 2 5 9 21.5 536 126 8.8.19 2 7 3 44.4 255 92 11.4.32 2 7 45 42.7 515 126 12.9.32 2 7 6 36.9 811 151 12.3.31 2 7 9 29. 1252 165 11.6.25 2 9 3 54.1 384 9 13.4.4 2 9 45 52.8 883 152 16.1.43 2 9 6 46.3 1658 196 15..43 2 9 9 34.3 3152 28 13.4.34 3 5 3 25.2 318 7.4.12 3 5 45 16.6 614 133 8.4.14 3 5 6 7.3 827 157 2.9.15 3 5 9 15.3 922 171 5.6.19 3 7 3 33.7 664 134 9.9.18 3 7 45 22.3 1244 173 8.7.18 3 7 6 12.2 167 195 6.8. 3 7 9. 184 2 8.6.25 3 9 3 43.6 94 162 12..24 3 9 45 3.4 2119 2 8.1.23 3 9 6 16.9 3 239 9.8.27 3 9 9 24.8 3553 285..34 * Defined in figure 3 347 348 349 35 351 352 353 pressure during impact. Beusenberg et. al. (1) [31] simulated the influence of different neck models using data from head impacts recorded from the National Football League in the US. In that study, it was shown that for oblique impacts to the front of the head, the coupling between the head and the neck is only important for rotations in the sagital plane. These rotations are however strongly dependent on the boundary conditions of the neck. Hering and Derler () [32] performed both radial and oblique helmet impact tests on both a detached Hybrid III dummy head and a complete Hybrid III dummy. It could be concluded from their study, that the influence of the neck and body on the rotational effects in the head was different for different impact locations. They also 354 355 356 357 358 359 36 IJCrash 4 Vol. 9 No. 6 268 Woodhead Publishing Ltd

Injury tolerances for oblique impact helmet testing 6 Impact 1 5 m/s Change of resultant rotational velocity (rad/s) 5 4 3 7 m/s 9 m/s Strain =.25 Strain =. Strain =.15 5 15 25 3 35 4 HIC Figure 5 Results from impact 1, where the axes represents computed levels in the FE Hybrid III head model, and the straight lines are isostrain curves representing different strain levels in the brain tissue in the FE human head model. 6 Impact 2 Change in resultant rotational velocity (rad/s) 5 4 3 3 m/s 5 m/s 7 m/s 9 m/s Strain =.25 Strain =. Strain =.15 5 15 25 3 35 4 HIC Figure 6 Results from impact 2, where the axes represents computed levels in the FE Hybrid III head model, and the straight lines are isostrain curves representing different strain levels in the brain tissue in the FE human head model. 361 362 363 364 365 366 367 concluded that because the Hybrid III neck is much stiffer than the human neck, this presented a problem. It is proposed here that the influence of the neck ought to be investigated in future studies. CONCLUSIONS When comparing the FE Hybrid III head model kinematics with strains in the FE human brain tissue during oblique impacts, it can be concluded that rotational kinematics are as important as translational kinematics and should therefore be included in future head injury criteria. In Impact 1, changes in the rotational velocity provide a critical parameter. In Impact 2 and 3, changes in rotational velocity as well as the HIC value are important indicators (reiterating here, Impact 1 is to the top of the head inducing sagital plane rotation, Impact 2 is lateral inducing axial rotation, and Impact 3 is lateral inducing coronal plane rotation). 368 369 37 371 372 373 374 375 376 Woodhead Publishing Ltd 268 7 IJCrash 4 Vol. 9 No.

M Aare, S Kleiven and P Halldin 6 Impact 3 5 5 m/s 7 m/s 9 m/s Change of resultant rotational velocity (rad/s) 4 3 Strain =.15 5 15 25 3 35 4 HIC Figure 7 Results from impact 3, where the axes represents computed levels in the FE Hybrid III head model, and the straight lines are isostrain curves representing different strain levels in the brain tissue in the FE human head model. Strain =. Strain =.25 377 378 379 38 381 382 383 384 385 386 387 388 389 39 391 392 393 394 Table 3 Constants k 1 and k 2 and the regression coefficient r 2 from the regressions analysis. This study shows that it is possible to make a good prediction of the strain in the brain tissue using the peak change of rotational velocity and the HIC value. Testspecific injury thresholds should therefore include both rotational and translational parameters, such as the change in rotational velocity and HIC. The results presented here can be helpful in preventing head injuries through their consideration in the setting of future standards for oblique impact helmet tests. ACKNOWLEDGEMENT The authors wish to thank the Länsförsäkringar Corporation Research Found and VINNOVA for financial support. REFERENCES Impact 1 Impact 2 Impact 3 k 1 6.14* 3 7.26* 3 3.92* 3 k 2 1.32* 5 3.5* 5 7.41* 5 r 2.93.97.68 1. AARE, M and VON HOLST, H. (3). Injuries from Motorcycle- and Moped crashes in Sweden from 1987 to 1999, Injury Control and Safety Promotion, vol., No. 3, pp.131 138. 2. GENNARELLI, T A. (1983). Head Injury in Man and Experimental Animals: Clinical Aspects, Acta Neurochirurgica, Suppl. 32, pp. 1 13. 1983. 3. OTTE, D, CHINN, B, DOYLE, D, MÄKITUPA, S, STURROCK, K and SCHULLER, E. (1999). Contribution to Final Report of COST 327 Project, University of Hanover. 4. [ECE22.5] United Nations, Economic commission for Europe, working party WP29 on the Construction of Vehicles, Regulation 22, Uniform provisions concerning the approval of protective helmets for drivers and passengers of motorcycles and mopeds, Geneva, Initially passed 1958, amendment 3 1988, amendment 4 1995, amendment 5 1999. 5. [FMV218] Federal motor vehicle standard No 218, U.S. Government Printing Office, 1997. 6. [BS6658] British Standards Institution. (1985). Protective helmets for vehicle users, BS 6658, London. 7. AARE, M and HALLDIN, P. (3). A New Laboratory Rig for Evaluating Helmets subject to Oblique Impacts, Traffic Injury Prevention, Vol. 4, Issue 3, pp. 24 248. 8. MARGULIES, S S and THIBAULT, L E. (1992). A proposed Tolerance Criterion for Diffuse Axonal Injuries in Man. J. of Biomechanics, 25 (89), pp. 917 923. 9. DIMASI, F P, EPPINGER, R H and BANDAK, F A. (1995). Computational Analysis of Head Impact Response Under Car Crash Loadings. Proceedings to the 39th STAPP Car Crash Conference, San Diego, SAE 952718. 9. DiMasi, F P, EPPINGER, R H and BANDAK, F A. (1995). Computational Analysis of Head Impact Response Under Car Crash Loadings. Proceedings to the 39th STAPP Car Crash Conference, San Diego, SAE 952718.. UENO, K and MELVIN, J. (1995). Finite Element Model Study of Head Impact Based on Hybrid III Head Acceleration: The Effects of Rotational and Translational 395 396 397 398 399 4 41 42 43 44 45 46 47 48 49 4 411 412 413 414 415 416 417 418 419 4 421 422 423 424 425 426 427 428 IJCrash 4 Vol. 9 No. 8 268 Woodhead Publishing Ltd

Injury tolerances for oblique impact helmet testing 429 43 431 432 433 434 435 436 437 438 439 44 441 442 443 444 445 446 447 448 449 45 451 452 453 454 455 456 457 458 459 46 461 462 463 464 465 466 467 468 469 47 Acceleration, Journal of Biomechanical Engineering, August 1995, vol. 117 pp. 319 328. 11. KLEIVEN, S and VON HOLST, H. (3). Review and evaluation of head injury criteria, Proc. RTO Specialist Meeting, the NATO, Coblenz, Germany. 12. BAIN, B C and MEANEY, D F. (). Tissue-Level Thresholds for Axonal Damage in an Experimental Model of Central Nervous System White Matter Injury, J. Biomech. Eng., 16, pp. 615 622. 13. HOLBOURN, A H S. (1943). Mechanics of head injury, Lancet 2, October 9, pp. 438 441. 14. HALLQUIST, J O. (1998). LS-DYNA3D Theoretical Manual. Livermore Software Technology Corporation. 15. KLEIVEN, S. (3). Influence of Impact Direction to the Human Head in Prediction of Subdural Hematoma. Journal of Neurotrauma (4), pp. 365 379. 16. FREDRIKSSON, L A. (1996). A Finite Element Data Base for Occupant Substitutes, Doctoral Thesis. Dissertations No. 447, Division of Solid Mechanics, Dept. of Mech. Engineering, Linköping University, Sweden. 17. National Highway Traffic Safety Administration (NHTSA), Department of Transportation (DOT), (1972). Occupant Crash Protection Head Injury Criterion S6.2 of MVSS 571.8, Docket 69-7, Notice 17. NHTSA, Washington, DC. 18. Visible Human Project. Web site: www.nih.gov 14-8-3. 19. KLEIVEN, S and VON HOLST, H. (1). Consequences of brain size following impact in prediction of subdural hematoma evaluated with numerical techniques, in: IRCOBI Conference 1, Isle of Man (UK), 161 172.. KLEIVEN, S and VON HOLST, H. (2). Consequences of head size following trauma to the human head. Journal of Biomechanics 35 (2), 153 16, 2. 21. KLEIVEN, S and HARDY, W N. (2). Correlation of an FE Model of the Human Head with Experiments on Localized Motion of the Brain Consequences for Injury Prediction. SAE Paper No. 2S-76, Society of Automotive Engineers. 46 th Stapp Car Crash Journal. 22. YETTRAM, A L, GODFREY, N P M and CHINN, B P. (1994). Materials for motorcycle crash helmets a finite element parametric study, Plastic, Rubber and Composites Processing and Applications, 22(1994) pp. 215 221 23. GENNARELLI, T A, THIBAULT, L E, ADAMS, J H. (1982). Diffuse Axonal Injury and Traumatic Coma in the Primate, Ann. Neurol. 12, pp. 564 574. 24. GENNARELLI, T A, THIBAULT, L E, TOMEI, G. (1987). Directional dependence of axonal brain injury due to centroidal and non-centroidal acceleration, SAE Paper No. 872197, Proc. 31 st Stapp Car Crash Conference, Society of Automotive Engineers, Warrendale, PA, pp. 49 53. 25. NEWMAN, J A, SHEWCHENKO, N and WELBOURNE, E. (). A Proposed New Biomechanical Head Injury Assessment Function - The Maximum Power Index, Proc. 44th Stapp Car Crash Conf., SAE Paper No. -1-SC16. 26. BAIN, B C, BILLIAR, A C, SHREIBER, D I, MCINTOSH, T K and MEANEY, D F. (1996). In Vivo Mechanical Thresholds for Traumatic Axonal Damage. AGARD AMP Specialists Meeting on Impact Head Injury, New Mexico 7 9 November 1996. 27. DONNELLY, B R. (1998). Brain tissue material properties: A comparison of results. Biomechanical research: Experimental and computational, Proc. 26 th Int. Workshop, pp. 47 57. 28. HARDY, W N, FOSTER, C D, MASON, M J, YANG, K H, KING, A I and TASHMAN, S. (1). Investigation of Head Injury Mechanisms Using Neutral Density Technology and High- Speed Biplanar X-ray. SAE Paper No. 1S-52, in: 45th Stapp Car Crash Conf., Society of Automotive Engineers, pp. 337 368. 29. KING, A I, HARDY, W N, MASON, M J and TASHMAN, S. (2). Comparison of Relative Motion Between the Brain and Skull of the Human Cadaver for Rotation in the Coronal and Sagittal Planes, 4th World Congress of Biomechanics. August 2, Calgary, Alberta, Canada. 3. RUAN, J S, KHALIL, T and KING, A I. (1991). Human Head Dynamic Response to Side Impact by Finite Element Modeling, J. Biomechanical engineering, 113(3) pp. 276 283. 31. BEUSENBERG, M, SHEWCHENKO, N, NEWMAN, J A, DE LANGE, R and CAPPON, H. (1). In: IRCOBI Conference 1, Isle of Man (UK), 295 3. 32. HERING, A M and DERLER, S. (). Motorcycle helmet drop tests using a Hybrid III dummy, in: IRCOBI Conference, Montpellier, France, September, pp. 37 321. 471 472 473 474 475 476 477 478 479 48 481 482 483 484 485 486 487 488 489 49 491 492 493 494 495 496 497 498 499 5 51 52 53 54 55 56 57 58 59 5 511 512 Woodhead Publishing Ltd 268 9 IJCrash 4 Vol. 9 No.