Viscoelastic properties of gerbil tympanic membrane at very low frequencies

The tympanic membrane transfers sound waves in the ear canal to mechanical vibrations in the middle ear and cochlea. Good estimates of the mechanical properties of the tympanic membrane are important to obtain realistic models. Up till now, only limited resources about tympanic membrane viscoelastic properties are available in the literature. This study aimed to quantify the viscoelastic properties of gerbil tympanic membrane. Step indentations were applied with a custom indenter on four fresh, intact tympanic membranes and the resulting force relaxation was measured. The reduced relaxation functions were then fitted with two viscoelastic model representations: a 5-parameter Maxwell model and a model with a continuous relaxation spectrum. The average relaxation function is described by an initial rapid decrease of 6.5% with characteristic time 0.77 s, followed by a long term decrease with characteristic time 46 s that gradually tends stable till a total relaxation of 15%. The relaxation curves in the time domain were transformed to complex moduli in the frequency domain. It was found that these transformations yield information on strain-rate dependence only from quasi-static to the very lowest acoustic frequencies. Finally, relaxation and hysteresis were simulated in a finite element model with viscoelastic material properties.     

A method for measuring linearly viscoelastic properties of human tympanic membrane using nanoindentation

A viscoelastic nanoindentation technique was developed to measure both in-plane and through-thickness viscoelastic properties of human tympanic membrane (TM). For measurement of in-plane Young's relaxation modulus, the TM sample was clamped on a circular hole and a nanoindenter tip was used to apply a concentrated force at the center of the TM sample. In this setup, the resistance to nanoindentation displacement can be considered due primarily to the in-plane stiffness. The load-displacement curve obtained was used along with finite element analysis to determine the in-plane viscoelastic properties of TM. For measurements of Young's relaxation modulus in the through-thickness (out-of-plane) direction, the TM sample was placed on a relatively rigid solid substrate and nanoindentation was made on the sample surface. In this latter setup, the resistance to nanoindentation displacement arises primarily due to out-of-plane stiffness. The load-displacement curve obtained in this manner was used to determine the out-of-plane relaxation modulus using the method appropriate for viscoelastic materials. From our sample tests, we obtained the steady-state values for in-plane moduli as approximately 17.4 MPa and approximately 19.0 MPa for posterior and anterior portions of TM samples, respectively, and the value for through-thickness modulus as approximately 6.0 MPa for both posterior and anterior TM samples. Using this technique, the local out-of-plane viscoelastic modulus can be determined for different locations over the entire TM, and the in-plane properties can be determined for different quadrants of the TM.     

Mechano-acoustic determination of Young's modulus of articular cartilage

The compressive stiffness of an elastic material is traditionally characterized by its Young's modulus. Young's modulus of articular cartilage can be directly measured using unconfined compression geometry by assuming the cartilage to be homogeneous and isotropic. In isotropic materials, Young's modulus can also be determined acoustically by the measurement of sound speed and density of the material. In the present study, acoustic and mechanical techniques, feasible for in vivo measurements, were investigated to quantify the static and dynamic compressive stiffness of bovine articular cartilage in situ. Ultrasound reflection from the cartilage surface, as well as the dynamic modulus were determined with the recently developed ultrasound indentation instrument and compared with the reference mechanical and ultrasound speed measurements in unconfined compression (n=72). In addition, the applicability of manual creep measurements with the ultrasound indentation instrument was evaluated both experimentally and numerically. Our experimental results indicated that the sound speed could predict 47% and 53% of the variation in the Young's modulus and dynamic modulus of cartilage, respectively. The dynamic modulus, as determined manually with the ultrasound indentation instrument, showed significant linear correlations with the reference Young's modulus (r(2)=0.445, p<0.01, n=70) and dynamic modulus (r(2)=0.779, p<0.01, n=70) of the cartilage. Numerical analyses indicated that the creep measurements, conducted manually with the ultrasound indentation instrument, were sensitive to changes in Young's modulus and permeability of the tissue, and were significantly influenced by the tissue thickness. We conclude that acoustic parameters, i.e. ultrasound speed and reflection, are indicative to the intrinsic mechanical properties of the articular cartilage. Ultrasound indentation instrument, when further developed, provides an applicable tool for the in vivo detection of cartilage mechano-acoustic properties. These techniques could promote the diagnostics of osteoarthrosis.     

Dynamic stiffness and crossbridge action in muscle

Small sinusoidal vibrations at 300 HZ were applied to frog sartorius muscle to measure the dynamic stiffness (Young's modulus) throughout the course of tetanus. For a peak-to-peak amplitude of 0.4% the dynamic Young's modulus increased from 1.5 X 10(5) Nm-2 in the resting state to 2 X 10(7) Nm-2 in tetanus. After correction for the external connective tissue, the dynamic Young's modulus of the muscle was almost directly proportional to the tension throughout the development of tetanus. The ratio of dynamic Young's modulus to tensile stress thus remained constant (with a value at 300 Hz of approximately 100), consistently with Huxley and Simmon's identification of the crossbridges as the source of both tension and stiffness. For a single crossbridge the ratio of stiffness to tension was 8.2 X 10(7) m-1 at 300 Hz; it is deduced from literature data that the limiting value at high frequencies is about 1.6 X 10(8) m-1. This ratio is interpreted on Harrington's (1971) model to show that crossbridge action can be explained by a helix-coil transition of about 80 out of the 260 residues in each S-2 myosin strand. It is also shown that a helix-coil model can account for the observed rapid relaxation of muscle without invoking any complex behaviour of the crossbridge head.     

Dynamic stiffness and crossbridge action in muscle

Small sinusoidal vibrations at 300 Hz were applied to frog sartorius muscle to measure the dynamic stiffness (Young's modulus) throughout the course of tetanus. For a peak-to-peak amplitude of 0.4% the dynamic Young's modulus increased from 1.5×10⁵ Nm^–2 in the resting state to 2×10⁷ Nm^–2 in tetanus. After correction for the external connective tissue, the dynamic Young's modulus of the muscle was almost directly proportional to the tension throughout the development of tetanus. The ratio of dynamic Young's modulus to tensile stress thus remained constant (with a value at 300 Hz of approximately 100), consistently with Huxley and Simmons' identification of the crossbridges as the source of both tension and stiffness.For a single crossbridge the ratio of stiffness to tension was 8.2×10⁷ m^–1 at 300 Hz; it is deduced from literature data that the limiting value at high frequencies is about 1.6×10⁸ m^–1. This ratio is interpreted on Harrington's (1971) model to show that crossbridge action can be explained by a helix-coil transition of about 80 out of the 260 residues in each S-2 myosin strand. It is also shown that a helix-coil model can account for the observed rapid relaxation of muscle without invoking any complex behaviour of the crossbridge head.     

Local Viscoelastic Properties of Live Cells Investigated Using Dynamic and Quasi-Static Atomic Force Microscopy Methods

The measurement of viscoelasticity of cells in physiological environments with high spatio-temporal resolution is a key goal in cell mechanobiology. Traditionally only the elastic properties have been measured from quasi-static force-distance curves using the atomic force microscope (AFM). Recently, dynamic AFM-based methods have been proposed to map the local in vitro viscoelastic properties of living cells with nanoscale resolution. However, the differences in viscoelastic properties estimated from such dynamic and traditional quasi-static techniques are poorly understood. In this work we quantitatively reconstruct the local force and dissipation gradients (viscoelasticity) on live fibroblast cells in buffer solutions using Lorentz force excited cantilevers and present a careful comparison between mechanical properties (local stiffness and damping) extracted using dynamic and quasi-static force spectroscopy methods. The results highlight the dependence of measured viscoelastic properties on both the frequency at which the chosen technique operates as well as the interactions with subcellular components beyond certain indentation depth, both of which are responsible for differences between the viscoelasticity property maps acquired using the dynamic AFM method against the quasi-static measurements.     

Local Viscoelastic Properties of Live Cells Investigated Using Dynamic and Quasi-Static Atomic Force Microscopy Methods

The measurement of viscoelasticity of cells in physiological environments with high spatio-temporal resolution is a key goal in cell mechanobiology. Traditionally only the elastic properties have been measured from quasi-static force-distance curves using the atomic force microscope (AFM). Recently, dynamic AFM-based methods have been proposed to map the local in vitro viscoelastic properties of living cells with nanoscale resolution. However, the differences in viscoelastic properties estimated from such dynamic and traditional quasi-static techniques are poorly understood. In this work we quantitatively reconstruct the local force and dissipation gradients (viscoelasticity) on live fibroblast cells in buffer solutions using Lorentz force excited cantilevers and present a careful comparison between mechanical properties (local stiffness and damping) extracted using dynamic and quasi-static force spectroscopy methods. The results highlight the dependence of measured viscoelastic properties on both the frequency at which the chosen technique operates as well as the interactions with subcellular components beyond certain indentation depth, both of which are responsible for differences between the viscoelasticity property maps acquired using the dynamic AFM method against the quasi-static measurements.     

Dynamic response of immature bovine articular cartilage in tension and compression, and nonlinear viscoelastic modeling of the tensile response

Very limited information is currently available on the constitutive modeling of the tensile response of articular cartilage and its dynamic modulus at various loading frequencies. The objectives of this study were to (1) formulate and experimentally validate a constitutive model for the intrinsic viscoelasticity of cartilage in tension, (2) confirm the hypothesis that energy dissipation in tension is less than in compression at various loading frequencies, and (3) test the hypothesis that the dynamic modulus of cartilage in unconfined compression is dependent upon the dynamic tensile modulus. Experiment 1: Immature bovine articular cartilage samples were tested in tensile stress relaxation and cyclical loading. A proposed reduced relaxation function was fitted to the stress-relaxation response and the resulting material coefficients were used to predict the response to cyclical loading. Adjoining tissue samples were tested in unconfined compression stress relaxation and cyclical loading. Experiment 2: Tensile stress relaxation experiments were performed at varying strains to explore the strain-dependence of the viscoelastic response. The proposed relaxation function successfully fit the experimental tensile stress-relaxation response, with R2 = 0.970+/-0.019 at 1% strain and R2 = 0.992+/-0.007 at 2% strain. The predicted cyclical response agreed well with experimental measurements, particularly for the dynamic modulus at various frequencies. The relaxation function, measured from 2% to 10% strain, was found to be strain dependent, indicating that cartilage is nonlinearly viscoelastic in tension. Under dynamic loading, the tensile modulus at 10 Hz was approximately 2.3 times the value of the equilibrium modulus. In contrast, the dynamic stiffening ratio in unconfined compression was approximately 24. The energy dissipation in tension was found to be significantly smaller than in compression (dynamic phase angle of 16.7+/-7.4 deg versus 53.5+/-12.8 deg at 10(-3) Hz). A very strong linear correlation was observed between the dynamic tensile and dynamic compressive moduli at various frequencies (R2 = 0.908+/-0.100). The tensile response of cartilage is nonlinearly viscoelastic, with the relaxation response varying with strain. A proposed constitutive relation for the tensile response was successfully validated. The frequency response of the tensile modulus of cartilage was reported for the first time. Results emphasize that fluid-flow dependent viscoelasticity dominates the compressive response of cartilage, whereas intrinsic solid matrix viscoelasticity dominates the tensile response. Yet the dynamic compressive modulus of cartilage is critically dependent upon elevated values of the dynamic tensile modulus.     

Time-dependent nanomechanics of cartilage

In this study, atomic force microscopy-based dynamic oscillatory and force-relaxation indentation was employed to quantify the time-dependent nanomechanics of native (untreated) and proteoglycan (PG)-depleted cartilage disks, including indentation modulus E_ind, force-relaxation time constant τ, magnitude of dynamic complex modulus |E^∗|, phase angle δ between force and indentation depth, storage modulus E′, and loss modulus E″. At ∼2 nm dynamic deformation amplitude, |E^∗| increased significantly with frequency from 0.22 ± 0.02 MPa (1 Hz) to 0.77 ± 0.10 MPa (316 Hz), accompanied by an increase in δ (energy dissipation). At this length scale, the energy dissipation mechanisms were deconvoluted: the dynamic frequency dependence was primarily governed by the fluid-flow-induced poroelasticity, whereas the long-time force relaxation reflected flow-independent viscoelasticity. After PG depletion, the change in the frequency response of |E^∗| and δ was consistent with an increase in cartilage local hydraulic permeability. Although untreated disks showed only slight dynamic amplitude-dependent behavior, PG-depleted disks showed great amplitude-enhanced energy dissipation, possibly due to additional viscoelastic mechanisms. Hence, in addition to functioning as a primary determinant of cartilage compressive stiffness and hydraulic permeability, the presence of aggrecan minimized the amplitude dependence of |E^∗| at nanometer-scale deformation.     

The viscoelastic standard nonlinear solid model: predicting the response of the lumbar intervertebral disk to low-frequency vibrations

Due to the mathematical complexity of current musculoskeletal spine models, there is a need for computationally efficient models of the intervertebral disk (IVD). The aim of this study is to develop a mathematical model that will adequately describe the motion of the IVD under axial cyclic loading as well as maintain computational efficiency for use in future musculoskeletal spine models. Several studies have successfully modeled the creep characteristics of the IVD using the three-parameter viscoelastic standard linear solid (SLS) model. However, when the SLS model is subjected to cyclic loading, it underestimates the load relaxation, the cyclic modulus, and the hysteresis of the human lumbar IVD. A viscoelastic standard nonlinear solid (SNS) model was used to predict the response of the human lumbar IVD subjected to low-frequency vibration. Nonlinear behavior of the SNS model was simulated by a strain-dependent elastic modulus on the SLS model. Parameters of the SNS model were estimated from experimental load deformation and stress-relaxation curves obtained from the literature. The SNS model was able to predict the cyclic modulus of the IVD at frequencies of 0.01 Hz, 0.1 Hz, and 1 Hz. Furthermore, the SNS model was able to quantitatively predict the load relaxation at a frequency of 0.01 Hz. However, model performance was unsatisfactory when predicting load relaxation and hysteresis at higher frequencies (0.1 Hz and 1 Hz). The SLS model of the lumbar IVD may require strain-dependent elastic and viscous behavior to represent the dynamic response to compressive strain.     

Quasi-Static and Dynamic Nanoindentation of Some Selected Biomaterials

This study details an investigation of the viscoelastic behavior of some biomaterials （nacre, cattle horn and beetle cuticle） at lamellar length scales using quasi-static and dynamic nanoindentation techniques in the materials＇ Transverse Direction （TD） and Longitudinal Direction （LD）. Our results show that nacre exhibits high fracture toughness moving towards a larger cam- paniForm as the stress frequency varies from 10 Hz to 200 Hz. Elytra cuticle exhibits the least fracture toughness presenting little energy dissipation in TD. It was initially speculated that the fracture toughness of the subject materials would be directly related to energy-dissipating mechanisms （mechanical hysteresis）, but not the maximum value of the loss tangent tan＆ However, it was found that the materials＇ elastic modulus and hardness are similar in both the TD and LD when assessed using the quasi-static nanoindentation method, but not dynamic nanoindentation. It is believed that the reported results can be useful in the design of new crack arrest and damping materials based on biological counterparts.     

Estimation of the Young's modulus of articular cartilage using an arthroscopic indentation instrument and ultrasonic measurement of tissue thickness

We evaluated whether the use of cartilage thickness measurement would improve the ability of the arthroscopic indentation technique to estimate the intrinsic stiffness of articular cartilage. First, cartilage thickness and ultrasound reflection from the surface of bovine humeral head were registered in situ using a high-frequency ultrasound probe. Subsequently, cartilage was indented in situ at the sites of the ultrasound measurements using arthroscopic instruments with plane-ended and spherical-ended indenters. Finally, full-thickness cartilage disks (n=30) were extracted from the indented sites (thickness=799-1654microm) and the equilibrium Young's modulus was determined with a material testing device in unconfined compression geometry. We applied analytical and numerical indentation models for the theoretical correction of experimental indentation measurements. An aspect-ratio (the ratio of indenter radius to cartilage thickness) correction improved the correlation of the indenter force with the equilibrium Young's modulus from r(2)=0.488 to r(2)=0.642-0.648 (n=30) for the plane-ended indenter (diameter=1.000mm, height=0.300mm) and from r(2)=0.654 to r(2)=0.684-0.692 (n=30) for the spherical-ended indenter (diameter=0.500mm, height=0.100mm), depending on the indentation model used for the correction. The linear correlation between the ultrasound reflection and the Young's modulus was r(2)=0.400 (n=30). These results suggest that with large indenters, knowledge of the cartilage thickness improves the reliability of the indentation measurements, especially in pathological situations where cartilage thickness may be significantly lower than normal. Ultrasound measurements also provide diagnostically important information about cartilage thickness as well as knowledge of the integrity of the superficial zone of cartilage.     

Quasi-linear viscoelastic properties of costal cartilage using atomic force microscopy

Costal cartilage (CC) is one of the load-bearing tissues of the rib cage. Literature on material characterisation of the CC is limited. Atomic force microscopy (AFM) has been extremely successful in characterising the elastic properties of soft biomaterials such as articular cartilage and hydrogels, which are often the material of choice for cartilage models. But AFM data on CC are absent in the literature. In this study, AFM indentations using spherical beaded tips were performed on human CC to isolate the mechanical properties. A novel method was developed for modelling the relaxation indentation experiments based on Fung's quasi-linear viscoelasticity and a continuous relaxation spectrum. This particular model has been popular for uniaxial compression test data analysis. Using the model, the mean Young's modulus of CC was found to be about 2.17, 4.11 and 5.49?MPa for three specimens. A large variation of modulus was observed over the tissue. Also, the modulus values decreased with distance from the costochondral junction.     

Wideband MRE and static mechanical indentation of human liver specimen: Sensitivity of viscoelastic constants to the alteration of tissue structure in hepatic fibrosis

Despite the success of elastography in grading hepatic fibrosis by stiffness related noninvasive markers the relationship between viscoelastic constants in the liver and tissue structure remains unclear. We therefore studied the mechanical properties of 16 human liver specimens with different degrees of fibrosis, inflammation and steatosis by wideband magnetic resonance elastography (MRE) and static indentation experiments providing the specimens? static Young?s modulus (E), dynamic storage modulus (G′) and dynamic loss modulus (G″). A frequency-independent shear modulus μ and a powerlaw exponent α were obtained by fitting G′ and G″ using the two-parameter sprinpot model. The mechanical parameters were compared to the specimens? histology derived parameters such as degree of Fibrosis (F), inflammation score and fat score, amount of hydroxyproline (HYP) used for quantification of collagen, blood markers and presurgery in vivo function tests.     

Viscoelastic properties of single attached cells under compression

The viscoelastic properties of single, attached C2C12 myoblasts were measured using a recently developed cell loading device. The device allows global compression of an attached cell, while simultaneously measuring the associated forces. The viscoelastic properties were examined by performing a series of dynamic experiments over two frequency decades (0.1-10 Hz) and at a range of axial strains (approximately 10-40%). Confocal laser scanning microscopy was used to visualize the cell during these experiments. To analyze the experimentally obtained force-deformation curves, a nonlinear viscoelastic model was developed. The nonlinear viscoelastic model was able to describe the complete series of dynamic experiments using only a single set of parameters, yielding an elastic modulus of 2120 +/- 900 Pa for the elastic spring, an elastic modulus of 1960 +/- 1350 for the nonlinear spring, and a relaxation time constant of 0.3 +/- 0.12 s. To our knowledge, it is the first time that the global viscoelastic properties of attached cells have been quantified over such a wide range of strains. Furthermore, the experiments were performed under optimal environmental conditions and the results are, therefore, believed to reflect the viscoelastic mechanical behavior of cells, such as would be present in vivo.     

A comparison of healthy human and swine articular cartilage dynamic indentation mechanics

Articular cartilage is a multicomponent, poroviscoelastic tissue with nonlinear mechanical properties vital to its function. A consequent goal of repair or replacement of injured cartilage is to achieve mechanical properties in the repair tissue similar to healthy native cartilage. Since fresh healthy human articular cartilage (HC) is not readily available, we tested whether swine cartilage (SC) could serve as a suitable substitute for mechanical comparisons. To a first approximation, cartilage tissue and surgical substitutes can be evaluated mechanically as viscoelastic materials. Stiffness measurements (dynamic modulus, loss angle) are vital to function and are also a non-destructive means of evaluation. Since viscoelastic material stiffness is strongly strain rate dependent, stiffness was tested under different loading conditions related to function. Stiffness of healthy HC and SC specimens was determined and compared using two non-destructive, mm-scale indentation test modes: fast impact and slow sinusoidal deformation. Deformation resistance (dynamic modulus) and energy handling (loss angle) were determined. For equivalent anatomic locations, there was no difference in dynamic modulus. However, the HC loss angle was ~35% lower in fast impact and ~12% higher in slow sinusoidal mode. Differences seem attributable to age (young SC, older HC) but also to species anatomy and biology. Test mode-related differences in human-swine loss angle support use of multiple function-related test modes. Keeping loss angle differences in mind, swine specimens could serve as a standard of comparison for mechanical evaluation of e.g. engineered cartilage or synthetic repair materials.     

Determination of Poisson's ratio of articular cartilage by indentation using different-sized indenters

Articular cartilage is often characterized as an isotropic elastic material with no interstitial fluid flow during instantaneous and equilibrium conditions, and indentation testing commonly used to deduce material properties of Young's modulus and Poisson's ratio. Since only one elastic parameter can be deduced from a single indentation test, some other test method is often used to allow separate measurement of both parameters. In this study, a new method is introduced by which the two material parameters can be obtained using indentation tests alone, without requiring a secondary different type of test. This feature makes the method more suitable for testing small samples in situ. The method takes advantages of the finite layer effect. By indenting the sample twice with different-sized indenters, a nonlinear equation with the Poisson's ratio as the only unknown can be formed and Poisson's ratio obtained by solving the nonlinear equation. The method was validated by comparing the predicted Poisson's ratio for urethane rubber with the manufacturer's supplied value, and comparing the predicted Young's modulus for urethane rubber and an elastic foam material with modulii measured by unconfined compression. Anisotropic and nonhomogeneous finite-element (FE) models of the indentation were developed to aid in data interpretation. Applying the method to bovine patellar cartilage, the tissue Young's modulus was found to be 1.79 +/- 0.59 MPa in instantaneous response and 0.45 +/- 0.26 MPa in equilibrium, and the Poisson's ratio 0.503 +/- 0.028 and 0.463 +/- 0.073 in instantaneous and equilibrium, respectively. The equilibrium Poisson's ratio obtained in our work was substantially higher than those derived from biphasic indentation theory and those optically measured in an unconfined compression test. The finite element model results and examination of viscoelastic-biphasic models suggest this could be due to viscoelastic, inhomogeneity, and anisotropy effects.     

Systematic error in mechanical measures of damage during four-point bending fatigue of cortical bone

Accumulation of fatigue microdamage in cortical bone specimens is commonly measured by a modulus or stiffness degradation after normalizing tissue heterogeneity by the initial modulus or stiffness of each specimen measured during a preloading step. In the first experiment, the initial specimen modulus defined using linear elastic beam theory (LEBT) was shown to be nonlinearly dependent on the preload level, which subsequently caused systematic error in the amount and rate of damage accumulation measured by the LEBT modulus degradation. Therefore, the secant modulus is recommended for measurements of the initial specimen modulus during preloading. In the second experiment, different measures of mechanical degradation were directly compared and shown to result in widely varying estimates of damage accumulation during fatigue. After loading to 400,000 cycles, the normalized LEBT modulus decreased by 26% and the creep strain ratio decreased by 58%, but the normalized secant modulus experienced no degradation and histology revealed no significant differences in microcrack density. The LEBT modulus was shown to include the combined effect of both elastic (recovered) and creep (accumulated) strain. Therefore, at minimum, both the secant modulus and creep should be measured throughout a test to most accurately indicate damage accumulation and account for different damage mechanisms. Histology revealed indentation of tissue adjacent to roller supports, with significant sub-surface damage beneath large indentations, accounting for 22% of the creep strain on average. The indentation of roller supports resulted in inflated measures of the LEBT modulus degradation and creep. The results of this study suggest that investigations of fatigue microdamage in cortical bone should avoid the use of four-point bending unless no other option is possible.     

Reduced nucleus pulposus glycosaminoglycan content alters intervertebral disc dynamic viscoelastic mechanics

The intervertebral disc functions over a range of dynamic loading regimes including axial loads applied across a spectrum of frequencies at varying compressive loads. Biochemical changes occurring in early degeneration, including reduced nucleus pulposus glycosaminoglycan content, may alter disc mechanical behavior and thus may contribute to the progression of degeneration. The objective of this study was to determine disc dynamic viscoelastic properties under several equilibrium loads and loading frequencies, and further, to determine how reduced nucleus glycosaminoglycan content alters dynamic mechanics. We hypothesized that (1) dynamic stiffness would be elevated with increasing equilibrium load and increasing frequency, (2) the disc would behave more elastically at higher frequencies, and finally, (3) dynamic stiffness would be reduced at low equilibrium loads under all frequencies due to nucleus glycosaminoglycan loss. We mechanically tested control and chondroitinase ABC injected rat lumbar motion segments at several equilibrium loads using oscillatory loading at frequencies ranging from 0.05 to 5 Hz. The rat lumbar disc behaved non-linearly with higher dynamic stiffness at elevated compressive loads irrespective of frequency. Phase angle was not affected by equilibrium load, although it decreased as frequency was increased. Reduced glycosaminoglycan decreased dynamic stiffness at low loads but not at high equilibrium loads and led to increased phase angle at all loads and frequencies. The findings of this study demonstrate the effect of equilibrium load and loading frequencies on dynamic disc mechanics and indicate possible mechanical mechanisms through which disc degeneration can progress.     

Calibration of hyperelastic material properties of the human lumbar intervertebral disc under fast dynamic compressive loads

Under fast dynamic loading conditions (e.g. high-energy impact), the load rate dependency of the intervertebral disc (IVD) material properties may play a crucial role in the biomechanics of spinal trauma. However, most finite element models (FEM) of dynamic spinal trauma uses material properties derived from quasi-static experiments, thus neglecting this load rate dependency. The aim of this study was to identify hyperelastic material properties that ensure a more biofidelic simulation of the IVD under a fast dynamic compressive load. A hyperelastic material law based on a first-order Mooney-Rivlin formulation was implemented in a detailed FEM of a L2-L3 functional spinal unit (FSU) to represent the mechanical behavior of the IVD. Bony structures were modeled using an elasto-plastic Johnson-Cook material law that simulates bone fracture while ligaments were governed by a viscoelastic material law. To mimic experimental studies performed in fast dynamic compression, a compressive loading velocity of 1 m/s was applied to the superior half of L2, while the inferior half of L3 was fixed. An exploratory technique was used to simulate dynamic compression of the FSU using 34 sets of hyperelastic material constants randomly selected using an optimal Latin hypercube algorithm and a set of material constants derived from quasi-static experiments. Selection or rejection of the sets of material constants was based on compressive stiffness and failure parameters criteria measured experimentally. The two simulations performed with calibrated hyperelastic constants resulted in nonlinear load-displacement curves with compressive stiffness (7335 and 7079 N/mm), load (12,488 and 12,473 N), displacement (1.95 and 2.09 mm) and energy at failure (13.5 and 14.7 J) in agreement with experimental results (6551 ± 2017 N/mm, 12,411 ± 829 N, 2.1 ± 0.2 mm and 13.0 ± 1.5 J respectively). The fracture pattern and location also agreed with experimental results. The simulation performed with constants derived from quasi-static experiments showed a failure energy (13.2 J) and a fracture pattern and location in agreement with experimental results, but a compressive stiffness (1580 N/mm), a failure load (5976 N) and a displacement to failure (4.8 mm) outside the experimental corridors. The proposed method offers an innovative way to calibrate the hyperelastic material properties of the IVD and to offer a more realistic simulation of the FSU in fast dynamic compression.