1
|
Hadagali P, Cronin D. Enhancing the Biofidelity of an Upper Cervical Spine Finite Element Model within the Physiologic Range of Motion and Its Effect On the Full Ligamentous Neck Model Response. J Biomech Eng 2022; 145:1143325. [PMID: 35864785 DOI: 10.1115/1.4055037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Indexed: 11/08/2022]
Abstract
Contemporary finite element neck models are developed in a neutral posture; however, evaluation of injury risk for out-of-position impacts requires neck model repositioning to non-neutral postures, with much of the motion occurring in the upper cervical spine (UCS). Current neck models demonstrate a limitation in predicting the intervertebral motions within the UCS within the range of motion (ROM), while recent studies have highlighted the importance of including the tissue strains resulting from repositioning FE neck models to predict injury risk. In the current study, the ligamentous cervical spine from a contemporary neck model (GHBMC M50 v4.5) was evaluated in flexion, extension and axial rotation by applying moments from 0 to 1.5 Nm in 0.5 Nm increments, as reported in experimental studies and corresponding to the physiologic loading of the UCS. Enhancements to the UCS model were identified, including the C0-C1 joint-space and alar ligament orientation. Following geometric enhancements, an analysis was undertaken to determine the UCS ligament laxities, using a sensitivity study followed by an optimization study. The ligament laxities were optimized to UCS-level experimental data from the literature. The mean percent difference between UCS model response and experimental data improved from 55% to 23% with enhancements. The enhanced UCS model was integrated with a ligamentous cervical spine (LS) model and assessed with independent experimental data. The mean percent difference between the LS model and the experimental data improved from 46% to 35% with the integration of the enhanced UCS model.
Collapse
Affiliation(s)
- Prasannaah Hadagali
- Department of Mechanical and Mechatronics Engineering, University of Waterloo, 200 University Ave. West, Waterloo, Ontario, Canada, N2L 3G1
| | - Duane Cronin
- Department of Mechanical and Mechatronics Engineering, University of Waterloo, 200 University Ave. West, Waterloo, Ontario, Canada, N2L 3G1
| |
Collapse
|
2
|
Liang Z, Mo F, Zheng Z, Li Y, Tian Y, Jiang X, Liu T. Quantitative cervical spine injury responses in whiplash loading with a numerical method of natural neural reflex consideration. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2022; 219:106761. [PMID: 35344767 DOI: 10.1016/j.cmpb.2022.106761] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/12/2021] [Revised: 03/07/2022] [Accepted: 03/15/2022] [Indexed: 06/14/2023]
Abstract
BACKGROUND AND OBJECTIVE Neural reflex is hypothesized as a regulating step in spine stabilizing system. However, neural reflex control is still in its infancy to consider in the previous finite element analysis of head-neck system for various applications. The purpose of this study is to investigate the influences of neural reflex control on neck biomechanical responses, then provide a new way to achieve an accurate biomechanical analysis for head-neck system with a finite element model. METHODS A new FE head-neck model with detailed active muscles and spinal cord modeling was established and globally validated at multi-levels. Then, it was coupled with our previously developed neuromuscular head-neck model to analyze the effects of vestibular and proprioceptive reflexes on biomechanical responses of head-neck system in a typical spinal injury loading condition (whiplash). The obtained effects were further analyzed by comparing a review of epidemiologic data on cervical spine injury situations. RESULT The results showed that the active model (AM) with neural reflex control obviously presented both rational head-neck kinematics and tissue injury risk referring to the previous experimental and epidemiologic studies, when compared with the passive model (PM) without it. Tissue load concentration locations as well as stress/strain levels were both changed due to the muscle activation forces caused by neural reflex control during the whole loading process. For the bony structures, the AM showed a peak stress level accounting for only about 25% of the PM. For the discs, the stress concentrated location was transferred from C2-C6 in the PM to C4-C6 in the AM. For the spinal cord, the strain concentrated locations were transferred from C1 segment to around C4 segment when the effects of neural reflex control were implemented, while the gray matter and white matter peak strains were reduced to 1/3 and 1/2 of the PM, respectively. All these were well correlated with epidemiological studies on clinical cervical spine injuries. CONCLUSION In summary, the present work demonstrated necessity of considering neural reflex in FE analysis of a head-neck system as well as our model biofidelity. Overall results also verified the previous hypothesis and further quantitatively indicated that the muscle activation caused by neural reflex is providing a protection for the neck in impact loading by decreasing the strain level and changing the possible injury to lower spinal cord level to reduce injury severity.
Collapse
Affiliation(s)
- Ziyang Liang
- College of Mechanical and Vehicle Engineering, Hunan University, Changsha, Hunan 410082, China; Department of Orthopedics, The Second Xiangya Hospital of Central South University, Changsha, Hunan 410011, China
| | - Fuhao Mo
- College of Mechanical and Vehicle Engineering, Hunan University, Changsha, Hunan 410082, China.
| | - Zhefen Zheng
- College of Mechanical and Vehicle Engineering, Hunan University, Changsha, Hunan 410082, China
| | - Yuandong Li
- College of Mechanical and Vehicle Engineering, Hunan University, Changsha, Hunan 410082, China
| | - Ye Tian
- College of Mechanical and Vehicle Engineering, Hunan University, Changsha, Hunan 410082, China
| | - Xiaobing Jiang
- Department of Spine Surgery, Guangzhou University of Chinese medicine, Guangzhou, Guangdong 510405, China
| | - Tang Liu
- Department of Orthopedics, The Second Xiangya Hospital of Central South University, Changsha, Hunan 410011, China
| |
Collapse
|
3
|
Zheng Z, Wang Y, Wang T, Wu Y, Li Y. A Systematic Review and Meta-Analysis of the Facet Joint Orientation and Its Effect on the Lumbar. JOURNAL OF HEALTHCARE ENGINEERING 2022; 2022:2486745. [PMID: 35242295 PMCID: PMC8888088 DOI: 10.1155/2022/2486745] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Revised: 12/21/2021] [Accepted: 01/10/2022] [Indexed: 11/18/2022]
Abstract
Facet tropism is recognized as the difference in the positioning of the facet joints in association with each other in the sagittal plane. This guides to an imbalanced biomechanical force over the facet joints and the intervertebral disc during rotation and other physiological activities. A systematic review and meta-analysis of Web of Science, EMBASE, PubMed, Cochrane Library, SCOPUS, and CINHAL from 2004 to 2021 to recognize the related research studies was performed. The data for meta-analysis were obtained from multiple studies to get the combined effect of the facet tropism on the lumbar disc herniation (LDH) and the degenerative lumbar spondylolisthesis (LDS). 117 articles were incorporated in the systematic review, where 41 studies were selected for meta-analysis, out of which 7 studies were found eligible as per the inclusion criteria. When degenerative lumbar spondylolisthesis was compared with the normal group, 95% CI was observed at 1.94 (1.59, 2.28). There was a comparison of disc herniation with the normal group in L4/L5, with a 95% CI of 0.60 (0.05, 1.14). The L5/S1 disc herniation was compared with the normal group and was found to be 0.21 (-0.48, 0.90). Therefore, it was observed that facet tropism is related to lumbar disc herniation and degenerative lumbar spondylolisthesis. Our meta-analysis demonstrated a unique link between the facet tropism and the lumbar disk degeneration along with degenerative lumbar spondylolisthesis.
Collapse
Affiliation(s)
- Zhirui Zheng
- The Second Affiliated Hospital of Harbin Medical University Orthopedic Surgery Three Ward, Harbin, China
| | - Youqiang Wang
- The Second Affiliated Hospital of Harbin Medical University Orthopedic Surgery Three Ward, Harbin, China
| | - Tong Wang
- The Second Affiliated Hospital of Harbin Medical University Orthopedic Surgery Three Ward, Harbin, China
| | - Yue Wu
- The Second Affiliated Hospital of Harbin Medical University Orthopedic Surgery Three Ward, Harbin, China
| | - Yuhui Li
- The Second Affiliated Hospital of Harbin Medical University Orthopedic Surgery Three Ward, Harbin, China
| |
Collapse
|
4
|
A biomechanical investigation of lumbar interbody fusion techniques. J Mech Behav Biomed Mater 2021; 125:104961. [PMID: 34781226 DOI: 10.1016/j.jmbbm.2021.104961] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2020] [Revised: 11/01/2021] [Accepted: 11/06/2021] [Indexed: 11/24/2022]
Abstract
The anterior, posterior, transforaminal, and circumferential lumbar interbody fusions (ALIF, PLIF, TLIF, CLIF/360) are used to treat spondylolisthesis, trauma, and degenerative pathologies. This study aims to investigate the biomechanical effects of the lumbar interbody fusion techniques on the spine. A validated T12-sacrum lumbar spine finite-element model was used to simulate surgical fusion of L4-L5 segment using ALIF, PLIF with one and two cages, TLIF with unilateral and bilateral fixation, and CLIF/360. The models were simulated under pure-moment and combined (moment and compression) loadings to investigate the effect of different lumbar interbody fusion techniques on range of motion, forces transferred through the vertebral bodies, disc pressures, and endplate stresses. The range of motion of the lumbar spine was decreased the most for fusions with bilateral posterior instrumentations (TLIF, PLIF, and CLIF/360). The increase in forces transmitted through the vertebrae and increase in disc pressures were directly proportional to the range of motion. The discs superior to fusion were under higher pressure, which was attributed to adjacent segment degeneration in the superior discs. The increase in endplate stresses was directly proportional to the cross-sectional area and was greater in caudal endplates at the fusion level, which was attributed to cage subsidence. The response of the models was in line with overall clinical observations from the patients and can be further used for future studies, which aim to investigate the effect of geometrical and material variations in the spine. The model results will assist surgeons in making informed decisions when selecting fusion procedures based on biomechanical effects.
Collapse
|
5
|
Yoganandan N, Purushothaman Y, Humm J. Normalized vertebral-level specific range of motion corridors for female spines in rear impact. TRAFFIC INJURY PREVENTION 2021; 22:S137-S140. [PMID: 34699297 DOI: 10.1080/15389588.2021.1982593] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
OBJECTIVE It is well known that the biomechanical responses of female and male spines are different in rear impacts. Female-specific finite element models are being developed as improvements over generic models. Such advancements need female-specific segmental responses for validation. The objectives of the study were to develop vertebral level-specific range of motion corridors from female human cadaver head-neck complexes exposed to rear impact loading. METHODS Previously conducted experiments from five human cadaver head-neck complexes were used in this analysis-based study. Briefly, the female head-neck complexes were isolated at the second thoracic vertebral level from the whole body such that the skin and the surrounding tissues of the osteoligamentous complex were intact. The distal end was fixed to the platform of a min-sled testing device. The anterior angulation of T1 was at 25 degrees with respect to the horizontal axis to simulate the normal driver posture. The occipital condyles were directly superior to the T1 body, and the Frankfort plane was horizontal. Rear impact loading were applied at a velocity of 2.6 m/s. The range of motion was defined as the inter-segmental angle at each level of the subaxial spinal column, and it was obtained by tracking the motion of the retroreflective targets that were secured on vertebral bodies and lateral masses of C2 through C7 vertebrae. Data were normalized with respect to the fifth percentile female total body mass, and corridors were developed using the equal stress equal velocity approach and expressed as mean ± 1 standard deviation corridors for each segment. RESULTS The segmental motions of the subaxial cervical spinal column were such that the upper regions responded with flexion while the lower regions responded with extension during the initial accelerative loading phase of the impact, resulting in a non-physiological curvature. During the later phase, all segments were in extension. individual corridors are presented as temporal responses in the body of the manuscript. A comparison of the mean temporal responses at each segment are presented to depict the angulation motion differences within the spinal column. CONCLUSIONS The present corridors are unique to the female spines. Because female spines have significantly (p < 0.05) different biomechanical responses when compared to male spines, local anatomical differences exist between male and female spines, and field data and clinical studies show female bias to whiplash associated disorders under the rear impact of loading, the present set of corridors serve as a fundamental dataset for the validation of female-specific finite element models. Current computational models can also use these corridors for improved validation to add confidence in their outputs.
Collapse
Affiliation(s)
- Narayan Yoganandan
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, Wisconsin
- Department of Veterans Affairs, Zablocki VA Medical Center, Neuroscience Research, Milwaukee, Wisconsin
| | | | - John Humm
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, Wisconsin
| |
Collapse
|
6
|
Humm JR, Yoganandan N. Sagittal plane moment responses of the THOR-05F anthropomorphic test device. TRAFFIC INJURY PREVENTION 2021; 22:S155-S158. [PMID: 34686084 DOI: 10.1080/15389588.2021.1982602] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
OBJECTIVE Anthropomorphic test devices (ATD) are used in crashworthiness studies to advance safety in automotive, military, aviation, and other environments. The Test Device for Human Occupant Restraint (THOR) is an advancement over the widely used Hybrid III ATD. The female version THOR-05F is different from the male as it is not a scaled-down version of the male, and it is based on the recognition that the cervical spines (necks) of females have a different response than males. The objective of this study is to evaluate its response at dynamic rates of loading and compare it with previous postmortem human surrogate (PMHS) responses under sagittal plane bending. METHODS The head/neck assembly was separated from the thorax, and a lower neck plate was attached to the head/neck assembly to mount the preparation to the frame of an electro-hydraulic testing device. A custom upper neck interface plate was attached to a novel angular displacement test device that converted the linear motion of the vertical electrohydraulic piston to moment loading at the occipital condyle joint. The neck was preconditioned by applying a sinusoidal 10-degree flexion-extension cycle for 90 s and then three repeat dynamic tests at a target rate of 90 Nm/s. Flexion and extension tests were performed with and without the front and rear neck cables of the THOR-05F neck. Targets were fixed to the upper neck adapter plate, occipital condyle joint, mid-spine aluminum puck, and lower neck adapter plate. The targets' three-dimensional positions were measured using a seven-camera optical motion capture system. Upper neck load cell and occipital condyle potentiometer data were sampled at 20 kHz, and loading rates were determined by calculating the sagittal moment slope between 15% and 85% of the signal. RESULTS The mean occipital condyle angle versus sagittal moment response from the 12 tests (three tests each with and without cables and under flexion and extension) are given in the body of the manuscript. With and without cables, the loading rates for flexion tests were 89.3 ± 0.5 Nm/s and 86.3 ± 0.4 Nm/s, and for extension tests they were 90.8 ± 1.2 Nm/s and 88.0 ± 1.5 Nm/s. The average peak sagittal moments were 34.2 ± 0.3 Nm and 30.3 ± 0.2 Nm for flexion and 50.6 ± 0.3 Nm and 47.0 ± 0.3 Nm for extension tests. The mean peak occipital condyle angles were 23.5 ± 0.2 deg and 25.3 ± 0.1 deg for flexion and 22.7 ± 0.2 deg and 25.8 ± 0.1 deg for extension. CONCLUSION Using the angular motion as a basis and comparing it with the previously conducted PMHS tests, the THOR-05F neck has approximately twice the stiffness of the human under sagittal plane bending.
Collapse
Affiliation(s)
- John R Humm
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, Wisconsin
| | - Narayan Yoganandan
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, Wisconsin
- Department of Orthopaedic Surgery, Medical College of Wisconsin, Milwaukee, Wisconsin
- Department of Veterans Affairs, Neuroscience Research, Zablocki VA Medical Center, Milwaukee, Wisconsin
| |
Collapse
|
7
|
Sato F, Miyazaki Y, Morikawa S, Ferreiro Perez A, Schick S, Brolin K, Svensson M. The Effect of Seat Back Inclination on Spinal Alignment in Automotive Seating Postures. Front Bioeng Biotechnol 2021; 9:684043. [PMID: 34409020 PMCID: PMC8365515 DOI: 10.3389/fbioe.2021.684043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2021] [Accepted: 06/28/2021] [Indexed: 11/25/2022] Open
Abstract
Experimental studies have demonstrated a relationship between spinal injury severity and vertebral kinematics, influenced by the initial spinal alignment of automotive occupants. Spinal alignment has been considered one of the possible causes of gender differences in the risk of sustaining spinal injuries. To predict vertebral kinematics and investigate spinal injury mechanisms, including gender-related mechanisms, under different seat back inclinations, it is needed to investigate the effect of the seat back inclination on initial spinal alignment in automotive seating postures for both men and women. The purpose of this study was to investigate the effect of the seat back inclination on spinal alignments, comparing spinal alignments of automotive seating postures in the 20° and 25° seat back angle and standing and supine postures. The spinal columns of 11 female and 12 male volunteers in automotive seating, standing, and supine postures were scanned in an upright open magnetic resonance imaging system. Patterns of their spinal alignments were analyzed using Multidimensional Scaling presented in a distribution map. Spinal segmental angles (cervical curvature, T1 slope, total thoracic kyphosis, upper thoracic kyphosis, lower thoracic kyphosis, lumbar lordosis, and sacral slope) were also measured using the imaging data. In the maximum individual variances in spinal alignment, a relationship between the cervical and thoracic spinal alignment was found in multidimensional scaling analyses. Subjects with a more lordotic cervical spine had a pronounced kyphotic thoracic spine, whereas subjects with a straighter to kyphotic cervical spine had a less kyphotic thoracic spine. When categorizing spinal alignments into two groups based on the spinal segmental angle of cervical curvature, spinal alignments with a lordotic cervical spine showed significantly greater absolute average values of T1 slope, total thoracic kyphosis, and lower thoracic kyphosis for both the 20° and 25° seat back angles. For automotive seating postures, the gender difference in spinal alignment was almost straight cervical and less-kyphotic thoracic spine for the female subjects and lordotic cervical and more pronounced kyphotic thoracic spine for the male subjects. The most prominent influence of seatback inclination appeared in Total thoracic kyphosis, with increased angles for 25° seat back, 8.0° greater in spinal alignments with a lordotic cervical spine, 3.2° greater in spinal alignments with a kyphotic cervical spine. The difference in total thoracic kyphosis between the two seatback angles and between the seating posture with the 20° seat back angle and the standing posture was greater for spinal alignments with a lordotic cervical spine than for spinal alignments with a kyphotic cervical spine. The female subjects in this study had a tendency toward the kyphotic cervical spine. Some of the differences between average gender-specific spinal alignments may be explained by the findings observed in the differences between spinal alignments with a lordotic and kyphotic cervical spine.
Collapse
Affiliation(s)
- Fusako Sato
- Safety Research Division, Japan Automobile Research Institute, Tsukuba, Japan.,Department of Mechanics and Maritime Sciences, Chalmers University of Technology, Gothenburg, Sweden
| | - Yusuke Miyazaki
- Department of Systems and Control Engineering, Tokyo Institute of Technology, Tokyo, Japan
| | | | | | - Sylvia Schick
- Department of Forensic Epidemiology, Institute of Legal Medicine, Ludwig Maximilian University of Munich, Munich, Germany
| | - Karin Brolin
- Lightness by Design Aktiebolag (AB), Stockholm, Sweden
| | - Mats Svensson
- Department of Mechanics and Maritime Sciences, Chalmers University of Technology, Gothenburg, Sweden
| |
Collapse
|
8
|
Oakley PA, Betz JW, Harrison DE, Siskin LA, Hirsh DW. Radiophobia Overreaction: College of Chiropractors of British Columbia Revoke Full X-Ray Rights Based on Flawed Study and Radiation Fear-Mongering. Dose Response 2021; 19:15593258211033142. [PMID: 34421439 PMCID: PMC8375354 DOI: 10.1177/15593258211033142] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2021] [Revised: 06/17/2021] [Accepted: 06/22/2021] [Indexed: 01/11/2023] Open
Abstract
Fears over radiation have created irrational pressures to dissuade radiography use within chiropractic. Recently, the regulatory body for chiropractors practicing in British Columbia, Canada, the College of Chiropractors of British Columbia (CCBC), contracted Pierre Côté to review the clinical use of X-rays within the chiropractic profession. A "rapid review" was performed and published quickly and included only 9 papers, the most recent dating from 2005; they concluded, "Given the inherent risks of radiation, we recommend that chiropractors do not use radiographs for the routine and repeat evaluation of the structure and function of the spine." The CCBC then launched an immediate review of the use of X-rays by chiropractors in their jurisdiction. Member and public opinion were gathered but not presented to their members. On February 4, 2021, the College announced amendments to their Professional Conduct Handbook that revoked X-ray rights for routine/repeat assessment and management of patients with spine disorders. Here, we highlight current and historical evidence that substantiates that X-rays are not a public health threat. We also point out critical and insurmountable flaws in the single paper used to support irrational and unscientific policy that discriminates against chiropractors who practice certain forms of evidence-based X-ray-guided methods.
Collapse
Affiliation(s)
| | | | | | | | | | - International Chiropractors Association Rapid Response Research Review Subcommittee
- Private Practice, Newmarket, ON, Canada
- Private Practice, Boise, ID, USA
- CBP NonProfit, Inc, Eagle, ID, USA
- Private Practice, Green Brook, NJ, USA
- Private Practice, Laurel, MD, USA
| |
Collapse
|
9
|
Corrales MA, Cronin DS. Importance of the cervical capsular joint cartilage geometry on head and facet joint kinematics assessed in a Finite element neck model. J Biomech 2021; 123:110528. [PMID: 34082236 DOI: 10.1016/j.jbiomech.2021.110528] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Revised: 05/06/2021] [Accepted: 05/09/2021] [Indexed: 11/29/2022]
Abstract
Finite element human neck models (NMs) aim to predict neck response and injury at the tissue level; however, contemporary models are most often assessed using global response such as head kinematics. Additionally, many NMs are developed from subject-specific imaging with limited soft tissue resolution in small structures such as the facet joints in the neck. Such details may be critical to enable NMs to predict tissue-level response. In the present study, the capsular joint cartilage (CJC) geometry in a contemporary NM was enhanced (M50-CJC) based on literature data. The M50-CJC was validated at the segment and full neck levels and assessed using relative facet joint kinematics (FJK), capsular ligament (CL) and intervertebral disc (IVD) strains, a relative vertebral rotation assessment (IV-NIC) and head kinematics in frontal and rear impact. The validation ratings at the segment level increased from 0.60 to 0.64, with improvements for modes of deformation associated with the facet joints, while no difference was noted at the head kinematic level. The improved CJC led to increased FJK rotation (188%) and IVD strain (152.2%,) attributed to the reduced facet joint gap. Further enhancements of the capsular joint representation or a link between the FJK and CL injury risk are recommended. Enhancements at the tissue level demonstrated a large effect on the IVD strain, but were not apparent in global metrics such as head kinematics. This study demonstrated that a biofidelic and detailed geometrical representation of the CJC contributes significantly to the predicted joint response, which is critical to investigate neck injury risk at the tissue level.
Collapse
Affiliation(s)
- M A Corrales
- Department of Mechanical and Mechatronics Engineering, University of Waterloo, 200 University Avenue West, Waterloo, Canada
| | - D S Cronin
- Department of Mechanical and Mechatronics Engineering, University of Waterloo, 200 University Avenue West, Waterloo, Canada.
| |
Collapse
|
10
|
Yoganandan N, Purushothaman Y, Choi H, Baisden J, Rajasekaran D, Banerjee A, Jebaseelan D, Kurpad S. Biomechanical Study of Cervical Disc Arthroplasty Devices Using Finite Element Modeling. JOURNAL OF ENGINEERING AND SCIENCE IN MEDICAL DIAGNOSTICS AND THERAPY 2021; 4:021004. [PMID: 35832636 PMCID: PMC8597568 DOI: 10.1115/1.4049907] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Revised: 01/11/2021] [Indexed: 06/15/2023]
Abstract
Many artificial discs for have been introduced to overcome the disadvantages of conventional anterior discectomy and fusion. The purpose of this study was to evaluate the performance of different U.S. Food and Drug Administration (FDA)-approved cervical disc arthroplasty (CDA) on the range of motion (ROM), intradiscal pressure, and facet force variables under physiological loading. A validated three-dimensional finite element model of the human intact cervical spine (C2-T1) was used. The intact spine was modified to simulate CDAs at C5-C6. Hybrid loading with a follower load of 75 N and moments under flexion, extension, and lateral bending of 2 N·m each were applied to intact and CDA spines. From this work, it was found that at the index level, all CDAs except the Bryan disc increased ROM, and at the adjacent levels, motion decreased in all modes. The largest increase occurred under the lateral bending mode. The Bryan disc had compensatory motion increases at the adjacent levels. Intradiscal pressure reduced at the adjacent levels with Mobi-C and Secure-C. Facet force increased at the index level in all CDAs, with the highest force with the Mobi-C. The force generally decreased at the adjacent levels, except for the Bryan disc and Prestige LP in lateral bending. This study demonstrates the influence of different CDA designs on the anterior and posterior loading patterns at the index and adjacent levels with head supported mass type loadings. The study validates key clinical observations: CDA procedure is contraindicated in cases of facet arthroplasty and may be protective against adjacent segment degeneration.
Collapse
Affiliation(s)
- Narayan Yoganandan
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI 53226; VA Medical Center, Medical College of Wisconsin, Milwaukee, WI 53226; Department of Neurosurgery, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226
| | - Yuvaraj Purushothaman
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI 53226
| | - Hoon Choi
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI 53226
| | - Jamie Baisden
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI 53226
| | - Deepak Rajasekaran
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI 53226
| | - Anjishnu Banerjee
- Division of Biostatistics, Medical College of Wisconsin, Milwaukee, WI 53226
| | - Davidson Jebaseelan
- School of Mechanical Engineering, Vellore Institute of Technology, Chennai Campus, Chennai 600036, India
| | - Shekar Kurpad
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI 53226; VA Medical Center, Medical College of Wisconsin, Milwaukee, WI 53226
| |
Collapse
|
11
|
Yoganandan N, Baisden JL, John J, Saravana Kumar G, Banerjee A, Choi H. Vertebral Level-dependent Kinematics of Female and Male Necks Under G+x Loading. Mil Med 2021; 186:619-624. [PMID: 33499461 DOI: 10.1093/milmed/usaa312] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Revised: 08/11/2020] [Accepted: 08/28/2020] [Indexed: 11/12/2022] Open
Abstract
INTRODUCTION Size-matched volunteer studies report gender-dependent variations in spine morphology, and head mass and inertia properties. The objective of this study was to determine the influence of these properties on upper and lower cervical spine temporal kinematics during G+x loading. METHODS Parametrized three-dimensional head-neck finite element models were used, and impacts were applied at 1.8 and 2.6 m/s at the distal end. Details are given in the article. Contributions of population-based variations in morphological and mass-related variables on temporal kinematics were evaluated using sensitivity analysis. Influence of variations on time to maximum nonphysiological curve formation, and flexion of upper and extension of the lower spines were analyzed for male-like and female-like spines. RESULTS Upper and lower spines responded with initial flexion and extension, resulting in a nonphysiological curve. Time to maximum nonphysiological curve and range of motions (ROMs) of the cervical column ranged from 45 to 66 ms, and 30 to 42 deg. Vertebral depth and location of the head center of gravity (cg) along anteroposterior axis were most influential variables for the upper spine flexion. Location of head cg along anteroposterior axis had the greatest influence on the time of the curve. Both anteroposterior and vertical locations of head cg, disc height, vertebral depth, head mass, and size were influential for the lower spine extension kinematics. CONCLUSIONS Models with lesser vertebral depth, that is, female-like spines, experienced greater range of motions and pronounced nonphysiological curves. This results in greater distraction/stretch of the posterior upper spine complex, a phenomenon attributed to suboccipital headaches. Forward location of head cg along anteroposterior axis had the greatest influence on upper and lower spine motions and time of formation of the curve. Any increased anteroposterior location of cg attributable to head supported mass may induce greater risk of injuries/neck pain in women during G+x loading.
Collapse
Affiliation(s)
- Narayan Yoganandan
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, USA
| | - Jamie L Baisden
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, USA
| | - Jobin John
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, USA
| | - Gurunathan Saravana Kumar
- Department of Engineering Design, Indian Institute of Technology Madras, Sardar Patel Road, Chennai, 600036, India
| | - Anjishnu Banerjee
- Division of Biostatistics, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, USA
| | - Hoon Choi
- Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, USA
| |
Collapse
|
12
|
Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis. Biomech Model Mechanobiol 2020; 20:389-401. [PMID: 33221991 PMCID: PMC7979651 DOI: 10.1007/s10237-020-01403-7] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2020] [Accepted: 11/04/2020] [Indexed: 12/13/2022]
Abstract
There is an increased interest in studying the biomechanics of the facet joints. For in silico studies, it is therefore important to understand the level of reliability of models for outputs of interest related to the facet joints. In this work, a systematic review of finite element models of multi-level spinal section with facet joints output of interest was performed. The review focused on the methodology used to model the facet joints and its associated validation. From the 110 papers analysed, 18 presented some validation of the facet joints outputs. Validation was done by comparing outputs to literature data, either computational or experimental values; with the major drawback that, when comparing to computational values, the baseline data was rarely validated. Analysis of the modelling methodology showed that there seems to be a compromise made between accuracy of the geometry and nonlinearity of the cartilage behaviour in compression. Most models either used a soft contact representation of the cartilage layer at the joint or included a cartilage layer which was linear elastic. Most concerning, soft contact models usually did not contain much information on the pressure-overclosure law. This review shows that to increase the reliability of in silico model of the spine for facet joints outputs, more needs to be done regarding the description of the methods used to model the facet joints, and the validation for specific outputs of interest needs to be more thorough, with recommendation to systematically share input and output data of validation studies.
Collapse
|
13
|
Umale S, Yoganandan N, Kurpad SN. Development and validation of osteoligamentous lumbar spine under complex loading conditions: A step towards patient-specific modeling. J Mech Behav Biomed Mater 2020; 110:103898. [PMID: 32957203 DOI: 10.1016/j.jmbbm.2020.103898] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2020] [Revised: 05/13/2020] [Accepted: 05/30/2020] [Indexed: 01/04/2023]
Abstract
Finite-element models are used to investigate the biomechanics of normal, diseased and surgically fused spines. Generally, nominal spine geometries are used to understand the biomechanics, which has created a need for a technique that develops patient-specific lumbar spine geometries. In the current study, a lumbar spine (T12-Sacrum) was developed using a technique that facilitates geometrical morphing, which assists in incorporating patient-specific morphologies into the model. The model evaluations can be used to propose a biomechanically suitable lumbar spine fusion procedure for patients. This study focuses on the validation of the base model under pure-moment, pure-compression and combined-compression-and-moment loadings. Experimental data from the literature were used to validate the response of the model. The L1-L2, L2-L3, L3-L4, L4-L5 and L5-sacrum segments demonstrated a range of motion of 4.5, 4.0, 5.4, 5.0 and 8.9° in flexion; 3.0, 2.5, 3.6, 3.1 and 5.2° in extension; 6.2, 5.8, 6.4, 5.0 and 6.1° in right and left lateral bending; and 2.9, 3.0, 2.9, 1.9 and 2.5° in right and left axial rotation, all under 10 Nm pure-moment loading. The L1-L2, L2-L3, L3-L4, L4-L5 and L5-sacrum discs demonstrated compressions of 1.1, 1.4, 1.6, 1.4 and 0.9 mm under 1200 N follower- or pure-compression loading. With the combined loading of 280 N follower and 7.5 Nm moment, the L1-L5 model demonstrated 11.7, 7.2, 18.3 and 10.4 degrees of range of motion in flexion, extension, bending and rotation, respectively. The model results were in good agreement with corridors from six different experimental studies and can be used for future clinical studies.
Collapse
Affiliation(s)
- Sagar Umale
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, USA.
| | - Narayan Yoganandan
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, USA; Clement J. Zablocki VA Medical Center, Milwaukee, WI, USA; Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, USA
| | - Shekar N Kurpad
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, USA; Clement J. Zablocki VA Medical Center, Milwaukee, WI, USA; Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, USA
| |
Collapse
|
14
|
Hindman BJ, Dexter F, Gadomski BC, Bucx MJ. Sex-Specific Intubation Biomechanics: Intubation Forces Are Greater in Male Than in Female Patients, Independent of Body Weight. Cureus 2020; 12:e8749. [PMID: 32714687 PMCID: PMC7377029 DOI: 10.7759/cureus.8749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
Background Studies of head, neck, and cervical spine morphology and tissue material properties indicate that cervical spine biomechanics differ between adult males and females. These differences result in sex-specific cervical spine kinematics and injury patterns in response to standardized loading conditions. Because direct laryngoscopy and endotracheal intubation require the application of a load to the cervical spine, intubation biomechanics should be sex-specific. The aim of this study was to determine if intubation forces during direct laryngoscopy differ between male and female patients and, if so, is the difference independent of body weight. Methods We pooled original data from three previously published adult clinical intubation studies that used methodologically reliable intubation force measurements (measured total laryngoscope force applied to the tongue, and force values were insensitive to or accounted for other laryngoscope blade forces). All patients had undergone direct laryngoscopy and orotracheal intubation with a Macintosh 3 blade under general anesthesia. Patient data included sex, age, height, weight, and maximum intubation force. Least squares multivariable linear regression was performed between the dependent variable (maximum intubation force) and two independent variables (patient sex and patient weight). A third term was added for the interaction between patient sex and weight. Results Among all patients (males n=42, females n=59), the median intubation force was 42.2 N (25th to 75th percentiles: 31.5 to 57.4 N). While controlling for patient body weight, intubation force differed between the sexes; P=0.011, with greater intubation force in male patients. While controlling for patient sex, there was a positive association between patient body weight and intubation force; P=0.009. In addition, there was a significant interaction between patient sex and weight; P=0.002, with intubation force in male patients having greater dependence on body weight. The difference in intubation force between male and female patients who had the same body weight exceeded 5 N when body weight exceeded 75 kg, and intubation force differences between male and female patients increased as patient body weight increased. Additional analyses using robust regression and using body mass index instead of weight provided comparable results. Conclusion In adult patients, the biomechanics of direct laryngoscopy and intubation are sex-specific. Our findings support the need to control for patient sex and weight in future clinical and laboratory studies of the human cervical spine and head and neck biomechanics.
Collapse
|
15
|
John JD, Saravana Kumar G, Yoganandan N. Rear-Impact Neck Whiplash: Role of Head Inertial Properties and Spine Morphological Variations on Segmental Rotations. J Biomech Eng 2019; 141:2733244. [DOI: 10.1115/1.4043666] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2018] [Indexed: 12/17/2022]
Abstract
Whiplash injuries continue to be a concern in low-speed rear impact. This study was designed to investigate the role of variations in spine morphology and head inertia properties on cervical spine segmental rotation in rear-impact whiplash loading. Vertebral morphology is rarely considered as an input parameter in spine finite element (FE) models. A methodology toward considering morphological variations as input parameters and identifying the influential variations is presented in this paper. A cervical spine FE model, with its morphology parametrized using mesh morphing, was used to study the influence of disk height, anteroposterior vertebral depth, and segmental size, as well as variations in head mass, moment of inertia, and center of mass locations. The influence of these variations on the characteristic S-curve formation in whiplash response was evaluated using the peak C2–C3 flexion marking the maximum S-curve formation and time taken for the formation of maximum S-curve. The peak C2–C3 flexion in the S-curve formation was most influenced by disk height and vertebral depth, followed by anteroposterior head center of mass location. The time to maximum S-curve was most influenced by the anteroposterior location of head center of mass. The influence of gender-dependent variations, such as the vertebral depth, suggests that they contribute to the greater segmental rotations observed in females resulting in different S-curve formation from men. These results suggest that both spine morphology and head inertia properties should be considered to describe rear-impact responses.
Collapse
Affiliation(s)
- Jobin D. John
- Department of Engineering Design, Indian Institute of Technology Madras, Chennai 600036, India e-mail:
| | - Gurunathan Saravana Kumar
- Mem. ASME Department of Engineering Design, Indian Institute of Technology, Madras Chennai 600036, India e-mail:
| | - Narayan Yoganandan
- Fellow ASME Center for NeuroTrauma Research, Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI 53226 e-mail:
| |
Collapse
|
16
|
DeVogel N, Yoganandan N, Banerjee A, Pintar FA. Hierarchical process using Brier Score Metrics for lower leg injury risk curves in vertical impact. BMJ Mil Health 2019; 166:318-323. [PMID: 30709924 DOI: 10.1136/jramc-2018-001124] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2018] [Revised: 01/02/2019] [Accepted: 01/03/2019] [Indexed: 11/04/2022]
Abstract
INTRODUCTION Parametric survival models are used to develop injury risk curves (IRCs) from impact tests using postmortem human surrogates (PMHS). Through the consideration of different output variables, input parameters and censoring, different IRCs could be created. The purpose of this study was to demonstrate the feasibility of the Brier Score Metric (BSM) to determine the optimal IRCs and derive them from lower leg impact tests. METHODS Two series of tests of axial impacts to PMHS foot-ankle complex were used in the study. The first series used the metrics of force, time and rate, and covariates of age, posture, stature, device and presence of a boot. Also demonstrated were different censoring schemes: right and exact/uncensored (RC-UC) or right and uncensored/left (RC-UC-LC). The second series involved only one metric, force, and covariates age, sex and weight. It contained interval censored (IC) data demonstrating different censoring schemes: RC-IC-UC, RC-IC-LC and RC-IC-UC-LC. RESULTS For each test set combination, optimal IRCs were chosen based on metric-covariate combination that had the lowest BSM value. These optimal IRCs are shown along with 95% CIs and other measures of interval quality. Forces were greater for UC than LC data sets, at the same risk levels (10% used in North Atlantic Treaty Organisation (NATO)). All data and IRCs are presented. CONCLUSIONS This study demonstrates a novel approach to examining which metrics and covariates create the best parametric survival analysis-based IRCs to describe human tolerance, the first step in describing lower leg injury criteria under axial loading to the plantar surface of the foot.
Collapse
Affiliation(s)
- Nicholas DeVogel
- Division of Biostatistics, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - N Yoganandan
- Department of Neurosurgery, Center for Neuro-Trauma Research, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - A Banerjee
- Division of Biostatistics, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - F A Pintar
- Joint Department of Biomedical Engineering, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| |
Collapse
|
17
|
Cervical spine morphology and ligament property variations: A finite element study of their influence on sagittal bending characteristics. J Biomech 2019; 85:18-26. [PMID: 30704760 DOI: 10.1016/j.jbiomech.2018.12.044] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2018] [Revised: 12/07/2018] [Accepted: 12/26/2018] [Indexed: 12/30/2022]
Abstract
Cervical spine finite element models reported in biomechanical literature usually represent a static morphology. Not considering morphology as a model parameter limits the predictive capabilities for applications in personalized medicine, a growing trend in modern clinical practice. The objective of the study was to investigate the influence of variations in spinal morphology on the flexion-extension responses, utilizing mesh-morphing-based parametrization and metamodel-based sensitivity analysis. A C5-C6 segment was used as the baseline model. Variations of intervertebral disc height, facet joint slope, facet joint articular processes height, vertebral body anterior-posterior depth, and segment size were parametrized. In addition, material property variations of ligaments were considered for sensitivity analysis. The influence of these variations on vertebral rotation and forces in the ligaments were analyzed. The disc height, segmental size, and body depth were found to be the most influential (in the cited order) morphology variations; while among the ligament material property variations, capsular ligament and ligamentum flavum influenced vertebral rotation the most. Changes in disc height influenced forces in the posterior ligaments, indicating that changes in the anterior load-bearing column of the spine could have consequences on the posterior column. A method to identify influential morphology variations is presented in this work, which will help automation efforts in modeling to focus on variations that matter. This study underscores the importance of incorporating influential morphology parameters, easily obtained through computed tomography/magnetic resonance images, to better predict subject-specific biomechanical responses for applications in personalized medicine.
Collapse
|
18
|
Humm JR, Yoganandan N, Driesslein KG, Pintar FA. Three-dimensional kinematic corridors of the head, spine, and pelvis for small female driver seat occupants in near- and far-side oblique frontal impacts. TRAFFIC INJURY PREVENTION 2018; 19:S64-S69. [PMID: 30517033 DOI: 10.1080/15389588.2018.1498973] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2018] [Revised: 07/05/2018] [Accepted: 07/05/2018] [Indexed: 06/09/2023]
Abstract
OBJECTIVES Analyses of recent automotive accident data indicate an increased risk of injury for small female occupants compared to males in similar accidents. Females have been shown to be more susceptible to spinal injuries than males. To protect this more vulnerable population, advanced anthropomorphic test devices (ATDs) and computer human body models are being developed and require biofidelity curves for validation. The aim of this study is to generate female-specific 3D kinematic corridors in near- and far-side oblique frontal impacts for the head, spine, and pelvis. METHODS Eight specimens were procured and prescreened for mass, stature, and quantitative computed tomography bone mineral density and preexisting injuries to minimize biologic variability. Sets of 4 noncolinear retroreflective targets were placed on the back of the head; dorsal spine at T1, T8, and L2; and posterior sacrum. Instrumented computed tomography scans were obtained to measure the orientation and position of the markers relative to anatomic fiducials. The specimens were placed on a buck representative of a generic automotive driver's seat environment designed to minimize lower-extremity and pelvic motion. The buck was oriented such that the buck centerline was seated 30° from the impact vector in either a near- or far-side oblique frontal configuration. Preposition of the occupant was specified to the 50th percentile male H-point location, thigh and tibial angles, and torso angle. Impact was delivered via a servo-acceleration sled to the base of the buck with a 30 km/h 9 g trapezoidal pulse. Occupants were restrained by a standard 3-point belt that had a custom load-limiter device set to 2 kN at the D-ring side of the shoulder belt. Target motion was recorded at 1 kHz using a 3D optical motion capture system. Anatomic motion of the head, spine, and pelvis was calculated relative to the seat, and the average response was determined from 4 near-side and 4 far-side tests. The borders of the corridor were determined by calculating a standard deviational ellipse in the x, y, and z planes at each time step. RESULTS Plots of the biofidelity corridors for near- and far-side tests are shown in planes parallel to the seat from the lateral, rear, and overhead directions. Averaged peak excursions in the fore/aft and lateral directions are compared for the near- and far-side corridors. Near-side female and male tests are similarly compared. CONCLUSIONS In general, average peak excursions were greater in the far-side configuration than in the near-side configuration. Peak excursion results compared well with similar tests conducted on male postmortem human subjects (PMHS). The kinematic corridors developed in the current study serve as a set of biofidelity corridors for the development of current and future physical and computational surrogates.
Collapse
Affiliation(s)
- John R Humm
- a Department of Neurosurgery , Medical College of Wisconsin , Milwaukee , Wisconsin
| | - Narayan Yoganandan
- a Department of Neurosurgery , Medical College of Wisconsin , Milwaukee , Wisconsin
- b Department of Orthopaedic Surgery , Medical College of Wisconsin , Milwaukee , Wisconsin
| | - Klaus G Driesslein
- a Department of Neurosurgery , Medical College of Wisconsin , Milwaukee , Wisconsin
| | - Frank A Pintar
- a Department of Neurosurgery , Medical College of Wisconsin , Milwaukee , Wisconsin
| |
Collapse
|