Applied Bionics and Biomechanics
Volume 2018 (2018), Article ID 8626102, 8 pages
https://doi.org/10.1155/2018/8626102

Research Article

A Numerical Investigation of Risk Factors Affecting Lumbar Spine Injuries Using a Detailed Lumbar Model

Jiajia Zheng,¹ Liang Tang,¹ and Jingwen Hu²

¹School of Technology, Beijing Forestry University, Beijing 100083, China
²University of Michigan Transportation Research Institute, Ann Arbor, MI 48109, USA

Correspondence should be addressed to Liang Tang; tang-l04@mails.tsinghua.edu.cn and Jingwen Hu; jwhu@umich.edu

Received 9 December 2017; Accepted 5 March 2018; Published 17 April 2018

Academic Editor: Kiros Karamanidis

Copyright © 2018 Jiajia Zheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Recent field data showed that lumbar spine fractures occurred more frequently in late model vehicles than early ones in frontal crashes. However, the lumbar spine designs of the current crash test dummies are not accurate in human anatomy and have not been validated against any human/cadaver impact responses. In addition, the lumbar spines of finite element (FE) human models, including GHBMC and THUMS, have never been validated previously against cadaver tests. Therefore, this study developed a detailed FE lumbar spine model and validated it against cadaveric tests. To investigate the mechanism of lumbar spine injury in frontal crashes, effects of changing the coefficient of friction (COF), impact velocity, cushion thickness and stiffness, and cushion angle on the risk of lumbar spine injuries were analyzed based on a Taguchi array of design of experiments. The results showed that impact velocity is the most important factor in determining the risk of lumbar spine fracture (). After controlling the impact velocity, increases in the cushion thickness can effectively reduce the risk of lumbar spine fracture ().

1. Introduction

Safety designs in newer cars generally provided better protections to occupants than those in older cars. However, several recent field data analyses have shown that lumbar spine fractures occurred more frequently in late model vehicles than the early ones in frontal crashes [1, 2].

This increasing trend of lumbar spine fractures in frontal crashes is extremely concerning, because none of the current crash test programs, including FMVSS 208, US-NCAP, and IIHS, considered lumbar spine injury in their safety evaluation process. This is partially due to the fact that none of the current crash test dummies can accurately estimate lumbar spine fracture risks. Current crash test dummies were designed to focus on estimating injury risks to the head, neck, chest, and lower extremities, and their lumbar spine designs were not based on real human anatomy and were not validated against any human/cadaver impact responses. As a result, although lumbar spine loadings could be measured in a crash test, they may not necessarily reflect the real injury mechanism or injury risk in frontal crashes. In addition, lumbar spines of the available FE human models, such as GHBMC (Global Human Body Model Consortium) model and THUMS (Total Human Model for Safety) model, have never been validated previously against cadaver tests. More recently, Arun et al. [3] applied the stiffness values obtained from cadaver tests directly to the lumbar spine of the GHBMC simplified model. However, it is a rigid body-based lumbar spine model. It is lacking of detailed anatomical structures of the lumbar spine and cannot estimate strain and stress in the vertebrae.

Factors affecting the lumbar spine fracture have been discussed extensively in the literature. Results indicated that older occupants (65+ years) were five times more likely to sustain spinal injury compared to younger occupants [4], women had more lumbar spine fractures than men [5], and lower bone quality is associated with an increased number of lumbar spine fractures within the CIREN cases analyzed [6]. Lumbar spine posture was also found to be an important factor affecting lumbar spine injuries. It was reported that slouched posture may cause an increase in stress on the lumbar vertebrae [7] and more reclined postures are associated with a higher lumbar vertebrae fracture risk by reconstructing real-world motor vehicle crashes [8]. Crash pulse magnitude and shape are also crucial for determining lumbar spine injury risk. For example, cadaver tests under axial loading showed that triangular pulses produced pelvic fractures with no lumbar spine fractures, while a sigmoid-shaped pulse produced no pelvic fractures but did produce an L1 burst fracture [9]. In addition, more severe crash pulses may lead to higher lumbar spine injury risks [10].

The objectives of the present study were to develop a detailed FE lumbar spine model, validate it against cadaveric tests, and use it to investigate the effects of changing the coefficient of friction (COF), impact velocity, cushion thickness and stiffness, and cushion angle on the risk of lumbar spine injuries. This study could provide better understanding of how to design countermeasures to reduce occupant lumbar spine injuries in the new generation of vehicle models.

2. Lumbar Spine Model Development and Validation

2.1. Model Development

In this study, the GHBMC model was selected as the baseline model to be modified due to its accurate geometrical representation. In the original GHBMC lumbar spine model (Figure 1(a)), vertebrae (T12-L5) were represented by a rigid material, intervertebral discs were modeled by shell elements (Figure 2(a)), and no ligaments were modeled. In the modified model (Figure 1(b)), the intervertebral disc was separated into three parts including nucleus, annulus, and collagenous fibers (Figure 2(b)). A 0.5 mm layer of shell element was added to the outside of the original vertebrae, representing the cortical bone. Detailed ligaments modeled by beam elements were added between the vertebrae. Nucleus, annulus, and cancellous bones were modeled by solid elements. Collagenous fibers and ligaments were modeled by beam elements. The facet joint articulations were simulated as the tiebreak contact between adjacent endplates. Note that muscles were ignored in the current study.

(a) Original GHBMC model

(b) Modified model

Figure 1: Comparison between original and modified lumbar spine models.

(a) Disc in GHBMC (shell elements)

(b) Disc in current model (three parts)

Figure 2: Comparison of the discs between GHBMC and modified lumbar models.

The modified lumbar model (Figure 1(b)) was composed of 64,338 elements and 12,295 nodes. The cancellous bones and cortical bones were modeled by MAT_SIMPLIFIED_JOHNSON_COOK in LS_DYNA [11, 12] with 1.3% and 0.94% effective plastic failure strains, respectively. The endplates of the discs were modeled by an elastic material, and the endplates of vertebrae were modeled by MAT_SIMPLIFIED_JOHNSON_COOK. The nucleus and annulus were modeled by MAT_MOONEY-RIVLIN_RUBBER [11, 12]. Accordingly, the shear modulus constant A of the nucleus and annulus were 0.64 MPa and 0.24 MPa, and the constant B of nucleus and annulus were −0.16 MPa and −0.06 MPa, respectively [12]. A linear elastic material was used for ligaments. Selections of stiffness of ligaments (Table 1) were based on a previous study [13]. The collagenous fibers were represented by a nonlinear load-displacement curve obtained from the literature [14]. Because external lamellae are stiffer than internal lamellae, the fibers in different layers were weighted according to a previous research [15]. To tune the model responses, the material properties of the bones and ligaments were adjusted slightly to match the test results during model validations. The material properties of all the parts in the modified lumbar model are listed in Table 2.

Table 1: Stiffness of ligaments used in the current study (unit: N/mm).

Ligaments
T12-L1
L1-L2
L2-L3
L3-L4
L4-L5
L5-S1

ALL
32.9
32.4
20.8
39.5
40.5
32.9

PLL
10
17.1
36.6
10.6
25.8
10

LF
24.2
23
25.1
34.5
27.2
24.2

CL
31.7
42.5
33.9
32.3
30.6
31.7

ISL
12.1
10
9.6
18.1
8.7
12.1

SSL
15.1
23
24.8
34.8
18
15.1

ITL
15.1
23
24.8
34.8
18
15.1

Table 2: Material properties used in the current study.

Segments
MAT
R0 (t/mm³)
E (MPa)/K (N/mm)
Poisson ratio
A (MPa)
B (MPa)
N
C
PSFAIL
SIGMAX (MPa)

Cortical bone
98
1.83E − 09
11,740
0.3
106.7
100
0.1
1
9.40E − 03
150

Cancellous bone
98
1.80E − 10
259
0.25
1.71
20
1
1
1.30E − 02
1.985

Endplate vertebra
98
1.06E − 09
9450
0.3
5.67
100
1
3
1.89E − 02
7.1

Endplate disc
1
1.20E − 09
200
0.3

Disc annulus
27
1.20E − 09
0.49

0.24
−0.06

Disc nucleus
27
1.00E − 09
0.495

0.64
−0.16

Fiber
71
1.20E − 09
Curves

Ligaments
74
1.00E − 09

In LS_DYNA material keywords, the strain energy density of MAT_MOONEY-RIVLIN_RUBBER is defined as a function of constant A, constant B, and Poisson ratio [16].

2.2. Validation against Nonfailure Tests

Cadaveric tests from Demetropoulos et al. [17] were used to validate the modified lumbar model. It is important to adjust the initial posture of the lumbar spine to be consistent with the experimental data before validation [18]. Therefore, a presimulation was conducted to adjust the spine curvature to match the specimen pretest condition, as shown in Figure 3. Impact simulations were performed at different loading modes corresponding to the testing conditions. As a result, six simulations under compression, anterior shear, posterior shear, lateral shear, extension, and lateral bending were performed with the modified lumbar spine model. Each simulation was run at a displacement rate of 100 mm/sec. The maximum displacements were set to be the same as those in the tests. During these simulations, T12 was fixed to the fixture on the top, and L5 was attached to the fixture at the bottom. In the compression and shear conditions, load was applied from the lower fixture with the upper fixture rigidly constrained. In the bending condition, displacement was applied from the lower fixture, and a bending moment was applied to the superior end of the lumbar spine by a cable. Force, moment, and angular displacement were measured at the upper fixture.

(a) Initial lumbar spine posture in experiment

(b) Initial lumbar spine posture in simulation

Figure 3: Initial posture of the lumbar spine of cadaver experiment and simulation.

Figure 4 shows load-displacement curves in each loading condition for the modified lumbar spine model. In compression, posterior shear, lateral shear, extension and lateral bending, simulation results shown in red line fell within test corridors reasonably well.

(a) Compression

(b) Anterior shear

(c) Posterior shear

(d) Lateral shear

(e) Extension

(f) Lateral bending

Figure 4: Response comparisons between tests and simulations (red: simulation results, grey: test corridors).

2.3. Validation against Failure Tests

Cadaver tests with tissue failure conducted by Duma et al. [19] were used to validate the modified lumbar model. The validation simulations were set up under the same configurations as those in the dynamic compression tests (Figure 5). The first simulation used a lumbar spine motion segment (Figure 5(a)), and the second simulation used the entire lumbar spine (Figure 5(b)). During the simulations, the superior end of the lumbar spine was attached to the upper plate, and the inferior end was attached to the lower plate. A prescribed motion was applied on the upper plate while the lower plate was fixed. All the failure simulations were loaded at 1.0 m/s. Force and moment were calculated in these simulations under different loadings.

(a) Lumbar spine motion segment

(b) Whole lumbar spine

Figure 5: Test configurations for dynamic compression with failure.

As shown in Figure 6, the model-predicted force-displacement curve for the motion segment configuration fell within the experimental corridor. Good correlations between the simulation and the test for the entire lumbar spine compression were also achieved, as shown in Figure 7. In addition, the locations of fractures in the simulation were consistent with those in the test (underneath the T12 and L3), as shown in Figures 7(c) and 7(d).

(a) Force-displacement curve

(b) Fracture location in simulation (underneath the upper vertebra)

Figure 6: Lumbar spine motion segment test and simulation results under compression.

(a) Force-displacement curve

(b) Moment-displacement curve

(c) Fracture location in simulation (T12)

(d) Fracture location in simulation (L3)

Figure 7: Test and simulation results of whole lumbar spines.

3. Design of Experiment Analysis

A DoE analysis based on Taguchi Array was used to study the effects of multiple factors on the risk of lumbar spine injuries. In this study, 5 factors were investigated, including cushion stiffness, cushion thickness, cushion angle, the coefficient of friction, and impact velocity.

The setup of the simulations in the DoE analysis is shown in Figure 8 [20]. The upper mass was attached to the upper platform to simulate the mass of the torso, head-neck, and upper extremities. The interaction between the upper platform and the upper fixture was in a form of a laterally oriented cylinder. The lumbar spine was fixed at cranial and caudal ends to the upper fixture and lower fixture, respectively [20]. A cushion foam, of which the thickness was defined as h, was attached to the fixed plate as shown in Figure 8, and its density was varied for the DoE analysis. The fixed plate was constrained at a cushion angle of . The initial impact velocity was applied to the upper mass, upper platform, and impact cylinder, through the upper fixture, lumbar spine specimen, and lower fixture, and finally reached the lower platform, cushion foam, and the fixed plate (Figure 8).

Figure 8: Setup used in Taguchi design of experiment analysis.

The initial setup conditions of the impact velocity, cushion angle, thickness, and density were 0.1 m/s, 0 degree, 20 mm, and 20 kg/m³, respectively. The coefficient of friction (COF) between the lower platform and the seat cushion was ranged from 0 to 0.9. Four different levels were assigned for each factor (Table 3). Because most lumbar injuries were reported as vertebral fractures, the maximum principal strain in the bony structures was selected to evaluate the risk of lumbar fracture [21]. These setups resulted in a total of 16 simulations based on the Taguchi Array, as shown in Table 4. One-way analysis of variance (ANOVA) and analysis of covariance (ANCONA) were performed using SPSS 20.0.

Table 3: Magnitude of different levels for the five factors selected for DOE analysis.

Level
Velocity (m/s)
Cushion angle (degrees)
Cushion thickness (mm)
Cushion density (kg/m³)
Coefficient of friction

1
0.1
0
20
20
0

2
0.3
10
60
30
0.3

3
0.5
20
100
40
0.6

4
0.7
30
140
50
0.9

Table 4: Maximum principal strain of 16 simulations in DOE analysis.

Number
Velocity (mm/ms)
Cushion angle (°)
Coefficient of friction
Cushion thickness (mm)
Cushion density (kg/m³)
Maximum principal strain (E-3)

1
0.5
10
0.9
140
20
4.08

2
0.5
20
0.6
20
40
5.17

3
0.7
20
0.3
60
20
4.6

4
0.1
20
0.9
100
30
2.97

5
0.5
30
0
60
30
2.72

6
0.5
0
0.3
100
50
3.43

7
0.3
20
0
140
50
3.22

8
0.1
30
0.3
140
40
2.51

9
0.7
0
0.6
140
30
5.14

10
0.1
0
0
20
20
3.82

11
0.1
10
0.6
60
50
3.99

12
0.7
30
0.9
20
50
6.08

13
0.3
10
0.3
20
30
4.21

14
0.3
0
0.9
60
40
3.29

15
0.3
30
0.6
100
20
2.36

16
0.7
10
0
100
40
5.7

The average values of the maximum principal strains of different factor categories are shown in Figure 9. Impact velocity () had the most significant influence on the risk of lumbar injuries. With impact velocity varying between 0.1 and 0.5 m/s, the maximum principal strain fluctuated between 0.0032 and 0.0038, but it increased rapidly to 0.0053 when impact velocity rose to 0.7 m/s. Except the impact velocity, no other factor was statistically significant ().

Figure 9: The average maximum principal strain with respect to all 5 factors simulated using the Taguchi array.

Because the impact velocity dominated the lumbar fracture risk in the Taguchi array, a one-way ANOVA was then conducted for the remaining 4 factors by controlling the impact velocity as a constant. The average values of the maximum principal strains of different factor categories are shown in Figure 10. The cushion thickness () became significant for the lumbar injuries. An increase in the cushion thickness will lead to more energy absorption and in turn lower lumbar spine fracture risk. Even though other factors were not significant, some general trends have also been noted. For example, with an increase in the COF, the lumbar spine fracture risk generally increased. For the cushion angle, the highest injury risk occurred when the lumbar spine orientation was perpendicular to the cushion. These trends are widely consistent to the findings from other studies on lumbar spine injuries [22] and cervical spine injuries [23].

Figure 10: The average maximum principal strain with respect to the left 4 factors without impact velocity (impact velocity was 0.7 m/s).

4. Conclusions

A detailed modified FE lumbar spine model based on the GHBMC model was developed and validated against available cadaveric tests which appeared in the literature. Risk factors affecting lumbar spine injuries were investigated using the modified model. Results of the DoE analysis demonstrated that the impact velocity is a significant factor influencing the lumbar injuries. After controlling the impact velocity, cushion thickness is another significant factor influencing lumbar injuries. An increase of cushion thickness or decrease of cushion stiffness will reduce the lumbar spine fracture risk. Additionally, minimizing the COF between the padding and the lumbar spine can reduce the lumbar spine injury risk.

One of the limitations of this study is that human factors such as BMI (body mass index), sex, and age were not considered. Older occupants have higher risk of lumbar spine fractures, and their lumbar spine curvatures may be different to the younger adults [24]. In addition, muscles were not included in the current lumbar model. Future studies may investigate the effects of human factors and muscle activations on the lumbar spine injury risks using parametric human models [25, 26].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the financial support of the National Science Foundation of China under Grant no. 51605032, the National Science Foundation of Beijing under Grant no. 3174052, and the Fundamental Research Funds for the Central Universities no. 2017ZY31.

References

F. A. Pintar, W. Curry, D. Halloway, W. Malik, and N. Yoganandan, “Thoracolumbar spine fractures in frontal crashes: an update of field data,” Traffic Injury Prevention, vol. 15, pp. S254–S256, 2014.

M. R. Zonfrillo, C. M. Locey, S. R. Scarfone, and K. B. Arbogast, “Motor vehicle crash–related injury causation scenarios for spinal injuries in restrained children and adolescents,” Traffic Injury Prevention, vol. 15, Supplement 1, pp. S49–S55, 2014.

M. W. J. Arun, P. Hadagali, K. Driesslein, W. Curry, N. Yoganandan, and F. A. Pintar, “Biomechanics of lumbar motion-segments in dynamic compression,” Stapp Car Crash Journal, vol. 61, pp. 1–25, 2017.

L. E. Bilston, E. C. Clarke, and J. Brown, “Spinal injury in car crashes: crash factors and the effects of occupant age,” Injury Prevention, vol. 17, no. 4, pp. 228–232, 2011.

F. A. Pintar, N. Yoganandan, D. J. Maiman, M. Scarboro, and R. W. Rudd, “Thoracolumbar spine fractures in frontal impact crashes,” Association for the Advancement of Automotive Medicine, vol. 56, pp. 277–283, 2012.

M. Saffarzadeh, R. C. Hightower, J. W. Talton, A. N. Miller, J. D. Stitzel, and A. A. Weaver, “Multicenter analysis of ciren occupant lumbar bone mineral density and correlation with age and fracture incidence,” Traffic Injury Prevention, vol. 17, Supplement 1, pp. 34–41, 2016.

C. Du, Z. Mo, S. Tian et al., “Biomechanical investigation of thoracolumbar spine in different postures during ejection using a combined finite element and multi-body approach,” International Journal for Numerical Methods in Biomedical Engineering, vol. 30, no. 11, pp. 1121–1131, 2014.

D. A. Jones, J. P. Gaewsky, M. E. Kelley, A. A. Weaver, A. N. Miller, and J. D. Stitzel, “Lumbar vertebrae fracture injury risk in finite element reconstruction of ciren and nass frontal motor vehicle crashes,” Traffic Injury Prevention, vol. 17, Supplement 1, pp. 109–115, 2016.

F. A. Pintar, N. Yoganandan, and J. Moore, “Spine and pelvis fractures in axial Z-acceleration,” in Proceedings of 40th International Workshop on Human Subjects for Biomechanical Research, Savannah, GA, USA, 2012.

T. Adolph, M. Wisch, A. Eggers et al., “Analyses of thoracic and lumbar spine injuries in frontal impacts,” in IRCOBI Conference, Gothenburg, Sweden, 2013.

M. El-Rich, P. J. Arnoux, E. Wagnac, C. Brunet, and C. E. Aubin, “Finite element investigation of the loading rate effect on the spinal load-sharing changes under impact conditions,” Journal of Biomechanics, vol. 42, no. 9, pp. 1252–1262, 2009.

E. Wagnac, P. J. Arnoux, A. Garo, and C. E. Aubin, “Finite element analysis of the influence of loading rate on a model of the full lumbar spine under dynamic loading conditions,” Medical & Biological Engineering & Computing, vol. 50, no. 9, pp. 903–915, 2012.

F. A. Pintar, N. Yoganandan, T. Myers, A. Elhagediab, and A. Sances Jr., “Biomechanical properties of human lumbar spine ligaments,” Journal of Biomechanics, vol. 25, no. 11, pp. 1351–1356, 1992.

A. Shirazi-Adl, “Biomechanics of the lumbar spine in sagittal/lateral moments,” Spine, vol. 19, no. 21, pp. 2407–2414, 1994.

D. Brickley-parsons and M. J. Glimcher, “Is the chemistry of collagen in intervertebral discs an expression of Wolff’s law? A study of the human lumbar spine,” Spine, vol. 9, no. 2, pp. 148–163, 1984.

B. Gladman, LS-DYNA Keywords User’s Manual Version, Livermore Software Technology Corporation, Livermore, CA, USA, 2007.

C. K. Demetropoulos, K. H. Yang, M. J. Grimm, T. B. Khalil, and A. I. King, “Mechanical properties of the cadaveric and hybrid iii lumbar spines,” Sae Technical Papers, vol. 28, no. 3, pp. 151–203, 1998.

J. D. M. Drake and J. P. Callaghan, “Do flexion/extension postures affect the in vivo passive lumbar spine response to applied axial twist moments?” Clinical Biomechanics, vol. 23, no. 5, pp. 510–519, 2008.

S. M. Duma, A. R. Kemper, D. M. Mcneely, P. G. Brolinson, and F. Matsuoka, “Biomechanical response of the lumbar spine in dynamic compression,” Biomedical Sciences Instrumentation, vol. 42, pp. 476–481, 2006.

B. D. Stemper, S. G. Storvik, N. Yoganandan et al., “A new pmhs model for lumbar spine injuries during vertical acceleration,” Journal of Biomechanical Engineering, vol. 133, no. 8, article 081002, 2011.

E. M. Porumb and F. Analoui, “A numerical investigation of factors affecting cervical spine injuries during rollover crashes,” Spine, vol. 33, no. 23, pp. 2529–2535, 2008.

R. Moradi, H. K. Beheshti, and H. M. Lankarani, “Lumbar load attenuation for rotorcraft occupants using a design methodology for the seat impact energy-absorbing system,” Open Engineering, vol. 2, no. 4, pp. 562–577, 2012.

J. Hu, K. H. Yang, C. C. Chou, and A. I. King, “A numerical investigation of factors affecting cervical spine injuries during rollover crashes,” Spine, vol. 33, no. 23, pp. 2529–2535, 2008.

C. S. Parenteau, P. Zhang, S. Holcombe, and S. Wang, “Characterization of vertebral angle and torso depth by gender and age groups with a focus on occupant safety,” Traffic Injury Prevention, vol. 15, no. 1, pp. 66–72, 2014.

E. Hwang, J. Hu, C. Chen et al., “Development, evaluation, and sensitivity analysis of parametric finite element whole-body human models in side impacts,” Stapp Car Crash Journal, vol. 60, pp. 473–508, 2016.

K. Zhang, L. Cao, A. Fanta et al., “An automated method to morph finite element whole-body human models with a wide range of stature and body shape for both men and women,” Journal of Biomechanics, vol. 60, pp. 253–260, 2017.