Comparing Predictive Accuracy and Computational Costs for Viscoelastic Modeling of Spinal Cord Tissues.

The constitutive equation used to characterize and model spinal tissues can significantly influence the conclusions from experimental and computational studies. Therefore, researchers must make critical judgments regarding the balance of computational efficiency and predictive accuracy necessary for their purposes. The objective of this study is to quantitatively compare the fitting and prediction accuracy of linear viscoelastic (LV), quasi-linear viscoelastic (QLV), and (fully) nonlinear viscoelastic (NLV) modeling of spinal-cord-pia-arachnoid-construct (SCPC), isolated cord parenchyma, and isolated pia-arachnoid-complex (PAC) mechanics in order to better inform these judgements. Experimental data collected during dynamic cyclic testing of each tissue condition were used to fit each viscoelastic formulation. These fitted models were then used to predict independent experimental data from stress-relaxation testing. Relative fitting accuracy was found not to directly reflect relative predictive accuracy, emphasizing the need for material model validation through predictions of independent data. For the SCPC and isolated cord, the NLV formulation best predicted the mechanical response to arbitrary loading conditions, but required significantly greater computational run time. The mechanical response of the PAC under arbitrary loading conditions was best predicted by the QLV formulation.

Viscoelasticity of spinal cord and meningeal tissues.

Compared to the outer dura mater, the mechanical behavior of spinal pia and arachnoid meningeal layers has received very little attention in the literature. This is despite experimental evidence of their importance with respect to the overall spinal cord stiffness and recovery following compression. Accordingly, inclusion of the mechanical contribution of the pia and arachnoid maters would improve the predictive accuracy of finite element models of the spine, especially in the distribution of stresses and strain through the cord's cross-section. However, to-date, only linearly elastic moduli for what has been previously identified as spinal pia mater is available in the literature. This study is the first to quantitatively compare the viscoelastic behavior of isolated spinal pia-arachnoid-complex, neural tissue of the spinal cord parenchyma, and intact construct of the two. The results show that while it only makes up 5.5% of the overall cross-sectional area, the thin membranes of the innermost meninges significantly affect both the elastic and viscous response of the intact construct. Without the contribution of the pia and arachnoid maters, the spinal cord has very little inherent stiffness and experiences significant relaxation when strained. The ability of the fitted non-linear viscoelastic material models of each condition to predict independent data within experimental variability supports their implementation into future finite element computational studies of the spine. STATEMENT OF SIGNIFICANCE: The neural tissue of the spinal cord is surrounded by three fibrous layers called meninges which are important in the behavior of the overall spinal-cord-meningeal construct. While the mechanical properties of the outermost layer have been reported, the pia mater and arachnoid mater have received considerably less attention. This study is the first to directly compare the behavior of the isolated neural tissue of the cord, the isolated pia-arachnoid complex, and the construct of these individual components. The results show that, despite being very thin, the inner meninges significantly affect the elastic and time-dependent response of the spinal cord, which may have important implications for studies of spinal cord injury.

The development and validation of a numerical integration method for non-linear viscoelastic modeling.

Compelling evidence that many biological soft tissues display both strain- and time-dependent behavior has led to the development of fully non-linear viscoelastic modeling techniques to represent the tissue's mechanical response under dynamic conditions. Since the current stress state of a viscoelastic material is dependent on all previous loading events, numerical analyses are complicated by the requirement of computing and storing the stress at each step throughout the load history. This requirement quickly becomes computationally expensive, and in some cases intractable, for finite element models. Therefore, we have developed a strain-dependent numerical integration approach for capturing non-linear viscoelasticity that enables calculation of the current stress from a strain-dependent history state variable stored from the preceding time step only, which improves both fitting efficiency and computational tractability. This methodology was validated based on its ability to recover non-linear viscoelastic coefficients from simulated stress-relaxation (six strain levels) and dynamic cyclic (three frequencies) experimental stress-strain data. The model successfully fit each data set with average errors in recovered coefficients of 0.3% for stress-relaxation fits and 0.1% for cyclic. The results support the use of the presented methodology to develop linear or non-linear viscoelastic models from stress-relaxation or cyclic experimental data of biological soft tissues.

Comparison of in vivo and ex vivo viscoelastic behavior of the spinal cord.

Despite efforts to simulate the in vivo environment, post-mortem degradation and lack of blood perfusion complicate the use of ex vivo derived material models in computational studies of spinal cord injury. In order to quantify the mechanical changes that manifest ex vivo, the viscoelastic behavior of in vivo and ex vivo porcine spinal cord samples were compared. Stress-relaxation data from each condition were fit to a non-linear viscoelastic model using a novel characterization technique called the direct fit method. To validate the presented material models, the parameters obtained for each condition were used to predict the respective dynamic cyclic response. Both ex vivo and in vivo samples displayed non-linear viscoelastic behavior with a significant increase in relaxation with applied strain. However, at all three strain magnitudes compared, ex vivo samples experienced a higher stress and greater relaxation than in vivo samples. Significant differences between model parameters also showed distinct relaxation behaviors, especially in non-linear relaxation modulus components associated with the short-term response (0.1-1 s). The results of this study underscore the necessity of utilizing material models developed from in vivo experimental data for studies of spinal cord injury, where the time-dependent properties are critical. The ability of each material model to accurately predict the dynamic cyclic response validates the presented methodology and supports the use of the in vivo model in future high-resolution finite element modeling efforts. STATEMENT OF SIGNIFICANCE: Neural tissues (such as the brain and spinal cord) display time-dependent, or viscoelastic, mechanical behavior making it difficult to model how they respond to various loading conditions, including injury. Methods that aim to characterize the behavior of the spinal cord almost exclusively use ex vivo cadaveric or animal samples, despite evidence that time after death affects the behavior compared to that in a living animal (in vivo response). Therefore, this study directly compared the mechanical response of ex vivo and in vivo samples to quantify these differences for the first time. This will allow researchers to draw more accurate conclusions about spinal cord injuries based on ex vivo data (which are easier to obtain) and emphasizes the importance of future in vivo experimental animal work.

Associations between in-vivo glenohumeral joint motion and morphology.

Joint morphology has a significant influence on joint motion and may contribute to the development of rotator cuff pathology, but the relationships between glenohumeral joint (GHJ) morphology and in-vivo GHJ motion are not well understood. The objectives of this study were to assess measures of joint morphology and their relationship with in-vivo joint motion in two populations: shoulders with intact rotator cuffs (n=48) and shoulders with rotator cuff pathology (n=36, including 5 symptomatic tears, 9 asymptomatic tears and 22 repaired tears). GHJ morphology was measured from CT-based three-dimensional models of the humerus and scapula. In-vivo GHJ motion was measured during shoulder abduction using biplane x-ray imaging. Associations between GHJ morphology and motion were assessed with univariate and best subsets regression. The only morphological difference identified between the populations was the critical shoulder angle (intact: 34.5 ± 4.7°, pathologic: 36.9 ± 5.0°, p=0.03), which is consistent with previous research. In intact shoulders, the superior/inferior (S/I) position of the humerus on the glenoid during shoulder abduction was significantly associated with the glenoid's S/I radius of curvature (p<0.01), conformity index (p<0.01), and stability angle (p<0.01). Furthermore, the S/I position of the humerus on the glenoid was negatively associated with the critical shoulder angle (p=0.04), which contradicts previous research. No significant associations between GHJ morphology and GHJ motion were detected in shoulders with rotator cuff tears. It is unknown if rotator cuff pathology compromises the relationships between GHJ morphology and motion, or if the absence of this relationship is a pre-existing condition that increases the likelihood of pathology.

Effects of footwear on three-dimensional tibiotalar and subtalar joint motion during running.

Running is a popular form of recreation, but injuries are common and may be associated with abnormal joint motion. The objective of this study was to determine the effect of three footwear conditions - barefoot (BF), an ultraflexible training shoe (FREE), and a motion control shoe (MC) - on 3D foot and ankle motion. Dynamic, biplane radiographic images were acquired from 12 runners during overground running. 3D rotations of the tibiotalar and subtalar joints were quantified in terms of plantarflexion/dorsiflexion (PF/DF), inversion/eversion (IN/EV) and internal/external rotation (IR/ER). Across the early stance phase (defined as footstrike to heel-off), BF running demonstrated greater tibiotalar joint range of motion for PF/DF (28.2 ± 8.3°) and IR/ER (7.0 ± 1.4°) than the shod conditions (FREE: PF/DF=15.1 ± 5.9°, IR/ER=4.8 ± 2.1°; MC: PF/DF=15.0 ± 6.2°, IR/ER=4.3 ± 0.7°). Also at the tibiotalar joint, BF running resulted in a position significantly more plantarflexed (BF: 2.0 ± 12.5°, FREE: 15.7 ± 12.2°, MC: 16.5 ± 9.3°) and internally rotated (BF: 12.9 ± 4.5°, FREE: 10.7 ± 4.3°, MC: 10.6 ± 3.9°) at footstrike compared to both shod conditions. No differences were detected between the shod conditions at any point in the early stance phase at the tibiotalar joint. The MC condition demonstrated significant differences compared to FREE at several points throughout the early stance phase at the subtalar joint, with the greatest differences seen at 30% in PF/DF (MC -1.4 ± 8.8°: FREE: -0.5 ± 9.0°), IN/EV (MC -8.1 ± 5.7°: FREE -6.3 ± 5.5°) and IR/ER (MC -9.5 ± 5.3°: FREE: -8.7 ± 5.2°). These findings indicate that footwear has subtle effects on joint motion mainly between BF and shod conditions at the tibiotalar joint and between shod conditions at the subtalar joint.