Utilization of Finite Element Analysis Techniques for Adolescent Idiopathic Scoliosis Surgical Planning

Adolescent Idiopathic Scoliosis, a three-dimensional deformity of the thoracolumbar spine, affects approximately 1-3% of patients ages 10-18. Surgical correction and treatment of the spinal column is a costly and high-risk task that is consistently complicated by factors such as patient-specific spinal deformities, curve flexibility, and surgeon experience. The following dissertation utilizes finite element analysis to develop a cost-effective, building-block approach by which surgical procedures and kinematic evaluations may be investigated. All studies conducted are based off a volumetric, thoracolumbar finite element (FE) model developed from computer-aided design (CAD) anatomy whose components are kinematically validated with in-vitro data. Spinal ligament stiffness properties derived from the literature are compared for kinematic assessment of a thoracic functional spinal unit (FSU) and benchmarked with available in-vitro kinematic data. Once ligament stiffness properties were selected, load sharing among soft tissues (e.g., ligaments and intervertebral disc) within the same FSU is then assessed during individual steps of a posterior correction procedure commonly used on scoliosis patients. Finally, the entire thoracolumbar spine is utilized to mechanically induce a mild scoliosis profile through an iterative preload and growth procedure described by the Hueter-Volkmann law. The mild scoliosis model is then kinematically compared with an asymptomatic counterpart. The thoracic deformation exhibited in the mild scoliosis model compared well with available CT datasets. Key findings of the studies confirm the importance of appropriately assigning spinal ligament properties with traditional toe and linear stiffness regimes to properly characterize thoracic spine FE models. Stiffness properties assigned within spinal FE models may also alter how intact ligaments and intervertebral discs respond to external loads during posterior correction procedures involving serial ligament removal, and thus can affect any desired post-surgical outcomes. Lastly, the thoracolumbar spine containing mild scoliosis experiences up to a 37% reduction in global range of motion compared to an asymptomatic spine, while also exhibiting larger decreases in segmental axial rotations at apical deformity levels. Future studies will address kinematic behavior of a severe scoliosis deformity and set the stage for column-based osseoligamentous load sharing assessments during surgical procedures

Vertical Drop Testing and Simulation of Anthropomorphic Test Devices

A series of 14 vertical impact tests were conducted using Hybrid III 50th Percentile and Hybrid II 50th Percentile Anthropomorphic Test Devices (ATDs) at NASA Langley Research Center. The purpose of conducting these tests was threefold: to compare and contrast the impact responses of Hybrid II and Hybrid III ATDs under two different loading conditions, to compare the impact responses of the Hybrid III configured with a nominal curved lumbar spine to that of a Hybrid III configured with a straight lumbar spine, and to generate data for comparison with predicted responses from two commercially available ATD finite element models. The two loading conditions examined were a high magnitude, short duration acceleration pulse, and a low magnitude, long duration acceleration pulse, each created by using different paper honeycomb blocks as pulse shape generators in the drop tower. The test results show that the Hybrid III results differ from the Hybrid II results more for the high magnitude, short duration pulse case. The comparison of the lumbar loads for each ATD configuration show drastic differences in the loads seen in the spine. The analytical results show major differences between the responses of the two finite element models. A detailed discussion of possible sources of the discrepancies between the two analytical models is also provided

Evaluation of Material Models within LS-DYNA(Registered TradeMark) for a Kevlar/Epoxy Composite Honeycomb

The performance of material models to simulate a novel composite honeycomb Deployable Energy Absorber (DEA) was evaluated using the nonlinear explicit dynamic finite element code LS-DYNA(Registered TradeMark). Prototypes of the DEA concept were manufactured using a Kevlar/Epoxy composite material in which the fibers are oriented at +/-45 degrees with respect to the loading axis. The development of the DEA has included laboratory tests at subcomponent and component levels such as three-point bend testing of single hexagonal cells, dynamic crush testing of single multi-cell components, and impact testing of a full-scale fuselage section fitted with a system of DEA components onto multi-terrain environments. Due to the thin nature of the cell walls, the DEA was modeled using shell elements. In an attempt to simulate the dynamic response of the DEA, it was first represented using *MAT_LAMINATED_COMPOSITE_FABRIC, or *MAT_58, in LS-DYNA. Values for each parameter within the material model were generated such that an in-plane isotropic configuration for the DEA material was assumed. Analytical predictions showed that the load-deflection behavior of a single-cell during three-point bending was within the range of test data, but predicted the DEA crush response to be very stiff. In addition, a *MAT_PIECEWISE_LINEAR_PLASTICITY, or *MAT_24, material model in LS-DYNA was developed, which represented the Kevlar/Epoxy composite as an isotropic elastic-plastic material with input from +/-45 degrees tensile coupon data. The predicted crush response matched that of the test and localized folding patterns of the DEA were captured under compression, but the model failed to predict the single-cell three-point bending response

Calibration of Airframe and Occupant Models for Two Full-Scale Rotorcraft Crash Tests

Two full-scale crash tests of an MD-500 helicopter were conducted in 2009 and 2010 at NASA Langley's Landing and Impact Research Facility in support of NASA s Subsonic Rotary Wing Crashworthiness Project. The first crash test was conducted to evaluate the performance of an externally mounted composite deployable energy absorber under combined impact conditions. In the second crash test, the energy absorber was removed to establish baseline loads that are regarded as severe but survivable. Accelerations and kinematic data collected from the crash tests were compared to a system integrated finite element model of the test article. Results from 19 accelerometers placed throughout the airframe were compared to finite element model responses. The model developed for the purposes of predicting acceleration responses from the first crash test was inadequate when evaluating more severe conditions seen in the second crash test. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used to calibrate model results for the second full-scale crash test. This combination of heuristic and quantitative methods was used to identify modeling deficiencies, evaluate parameter importance, and propose required model changes. It is shown that the multi-dimensional calibration techniques presented here are particularly effective in identifying model adequacy. Acceleration results for the calibrated model were compared to test results and the original model results. There was a noticeable improvement in the pilot and co-pilot region, a slight improvement in the occupant model response, and an over-stiffening effect in the passenger region. This approach should be adopted early on, in combination with the building-block approaches that are customarily used, for model development and test planning guidance. Complete crash simulations with validated finite element models can be used to satisfy crash certification requirements, thereby reducing overall development costs

Simulating the Response of a Composite Honeycomb Energy Absorber

NASA has sponsored research to evaluate an externally deployable composite honeycomb designed to attenuate loads in the event of a helicopter crash. The concept, designated the Deployable Energy Absorber (DEA), is an expandable Kevlar(Registered TradeMark) honeycomb. The DEA has a flexible hinge that allows the honeycomb to be stowed collapsed until needed during an emergency. Evaluation of the DEA began with material characterization of the Kevlar(Registered TradeMark)-129 fabric/epoxy, and ended with a full-scale crash test of a retrofitted MD-500 helicopter. During each evaluation phase, finite element models of the test articles were developed and simulations were performed using the dynamic finite element code, LS-DYNA(Registered TradeMark). The paper will focus on simulations of two full-scale impact tests involving the DEA, a mass-simulator and a full-scale crash of an instrumented MD-500 helicopter. Isotropic (MAT24) and composite (MAT58) material models, which were assigned to DEA shell elements, were compared. Based on simulations results, the MAT58 model showed better agreement with test

Development and Calibration of a System-Integrated Rotorcraft Finite Element Model for Impact Scenarios

Two full-scale crash tests of an MD-500 helicopter were conducted in 2009 and 2010 at NASA Langley's Landing and Impact Research Facility in support of NASA s Subsonic Rotary Wing Crashworthiness Project. The first crash test was conducted to evaluate the performance of an externally mounted composite deployable energy absorber (DEA) under combined impact conditions. In the second crash test, the energy absorber was removed to establish baseline loads that are regarded as severe but survivable. The presence of this energy absorbing device reduced the peak impact acceleration levels by a factor of three. Accelerations and kinematic data collected from the crash tests were compared to a system-integrated finite element model of the test article developed in parallel with the test program. In preparation for the full-scale crash test, a series of sub-scale and MD-500 mass simulator tests were conducted to evaluate the impact performances of various components and subsystems, including new crush tubes and the DEA blocks. Parameters defined for the system-integrated finite element model were determined from these tests. Results from 19 accelerometers placed throughout the airframe were compared to finite element model responses. The model developed for the purposes of predicting acceleration responses from the first crash test was inadequate when evaluating more severe conditions seen in the second crash test. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used to calibrate model results for the full-scale crash test without the DEA. This combination of heuristic and quantitative methods identified modeling deficiencies, evaluated parameter importance, and proposed required model changes. The multidimensional calibration techniques presented here are particularly effective in identifying model adequacy. Acceleration results for the calibrated model were compared to test results and the original model results. There was a noticeable improvement in the pilot and copilot region, a slight improvement in the occupant model response, and an over-stiffening effect in the passenger region. One lesson learned was that this approach should be adopted early on, in combination with the building-block approaches that are customarily used, for model development and pretest predictions. Complete crash simulations with validated finite element models can be used to satisfy crash certification requirements, potentially reducing overall development costs