The mechanical properties of neonatal rat spinal cord in vitro,and comparisons with adult

A number of studies have investigated the mechanical properties of adult spinal cord under tension, however it is not known whether age has an effect on these properties. This is of interest to those aiming to understand the clinical differences between adults and children with spinal cord injury (e.g. severity and recovery), and those developing experimental or computational models for paediatric spinal cord injury. Entire spinal cords were freshly harvested from neonatal rats (14 days) and tested in vitro under uniaxial tension at a range of strain rates (0.2, 0.02, 0.002/s) to a range of strains (2%, 3.5%, 5%), with relaxation responses being recorded for 15 min. These mechanical properties were compared to previously reported data from similar experiments on adult rat spinal cords, and the peak stress and the stress after 15 min of relaxation were found to be significantly higher for spinal cords from adults than neonates (p<0.001). A non-linear viscoelastic model was developed and was observed to adequately predict the mechanical behaviour of this tissue. The model developed in this study may be of use in computational models of paediatric spinal cord. The significant differences between adult and neonatal spinal cord properties may explain the higher initial severity of spinal cord injury in children and may have implications for the development of experimental animal models for paediatric spinal cord injury, specifically for those aiming to match the injury severity with adult experimental models.     

Compression behavior of porcine spinal cord white matter

Spinal cord injury often results from a compressive load; however, the compression behavior of spinal cord white matter has not been clearly established. Quantifying the compression behavior is important for advancing our understanding of spinal cord injury mechanics and facilitating the use of finite element models to study injury. The objective of this study was to characterize the unconfined compression behavior of isolated white matter segments and determine the constitutive model which best captured the stress-strain behavior. Spinal cord white matter samples were harvested immediately following sacrifice from juvenile Yorkshire pigs (n=104). The samples were compressed to 40% strain at four strain rates (0.005, 0.05, 0.5, and 5.0/s) and allowed to relax for 60s. The effects of preload, peak strain, sample aspect ratio, and time post mortem on peak stress, and constitutive model parameters were also examined. Strain rate had a significant effect on peak stress (p<0.001). A first-order Ogden model best captured the loading response of spinal cord white matter (R(2)=0.99) and a viscoelastic material model combining a first-order Ogden model with a 3-term Prony series effectively captured the effect of strain rate and the relaxation response. This study showed spinal cord white matter to be less stiff than previously estimated by inverse finite element methods, which will have a significant effect on finite element model predictions of the magnitude and distribution of stresses and strains in the spinal cord. This study is the first to quantify the unconfined compression response of spinal cord white matter.     

Nonlinear viscoelastic characterization of bovine trabecular bone

The time-independent elastic properties of trabecular bone have been extensively investigated, and several stiffness–density relations have been proposed. Although it is recognized that trabecular bone exhibits time-dependent mechanical behaviour, a property of viscoelastic materials, the characterization of this behaviour has received limited attention. The objective of the present study was to investigate the time-dependent behaviour of bovine trabecular bone through a series of compressive creep–recovery experiments and to identify its nonlinear constitutive viscoelastic material parameters. Uniaxial compressive creep and recovery experiments at multiple loads were performed on cylindrical bovine trabecular bone samples (\(n = 19\)). Creep response was found to be significant and always comprised of recoverable and irrecoverable strains, even at low stress/strain levels. This response was also found to vary nonlinearly with applied stress. A systematic methodology was developed to separate recoverable (nonlinear viscoelastic) and irrecoverable (permanent) strains from the total experimental strain response. We found that Schapery’s nonlinear viscoelastic constitutive model describes the viscoelastic response of the trabecular bone, and parameters associated with this model were estimated from the multiple load creep–recovery (MLCR) experiments. Nonlinear viscoelastic recovery compliance was found to have a decreasing and then increasing trend with increasing stress level, indicating possible stiffening and softening behaviour of trabecular bone due to creep. The obtained parameters from MLCR tests, expressed as second-order polynomial functions of stress, showed a similar trend for all the samples, and also demonstrate stiffening–softening behaviour with increasing stress.     

Linear and quasi-linear viscoelastic characterization of ankle ligaments

The objective of this study was to produce linear and nonlinear viscoelastic models of eight major ligaments in the human ankle/foot complex for use in computer models of the lower extremity. The ligaments included in this study were the anterior talofibular (ATaF), anterior tibiofibular (ATiF), anterior tibiotalar (ATT), calcaneofibular (CF), posterior talofibular (PTaF), posterior tibiofibular (PTiF), posterior tibiotalar (PTT), and tibiocalcaneal (TiC) ligaments. Step relaxation and ramp tests were performed. Back-extrapolation was used to correct for vibration effects and the error introduced by the finite rise time in step relaxation tests. Ligament behavior was found to be nonlinear viscoelastic, but could be adequately modeled up to 15 percent strain using Fung's quasilinear viscoelastic (QLV) model. Failure properties and the effects of preconditioning were also examined.     

On the viscoelastic character of liver tissue: experiments and modelling of the linear behaviour

This paper presents the results of recent investigations into the viscoelastic character of liver tissue, especially the linear response. The linear domain of liver tissue has been determined by a strain sweep oscillation test. The fluid-like viscoelastic behaviour of liver tissue was confirmed by stress relaxation experiments. To demonstrate the consistency and quality of the experimental data, the relaxation spectrum was used to reconstruct the relaxation modulus, as well as the storage and loss moduli. The generalised Maxwell model was successfully employed to model the mechanical behaviour in the linear viscoelastic domain.     

A constituent-based model for the nonlinear viscoelastic behavior of ligaments

This paper presents a constitutive model for predicting the nonlinear viscoelastic behavior of soft biological tissues and in particular of ligaments. The constitutive law is a generalization of the well-known quasi-linear viscoelastic theory (QLV) in which the elastic response of the tissue and the time-dependent properties are independently modeled and combined into a convolution time integral. The elastic behavior, based on the definition of anisotropic strain energy function, is extended to the time-dependent regime by means of a suitably developed time discretization scheme. The time-dependent constitutive law is based on the postulate that a constituent-based relaxation behavior may be defined through two different stress relaxation functions: one for the isotropic matrix and one for the reinforcing (collagen) fibers. The constitutive parameters of the viscoelastic model have been estimated by curve fitting the stress relaxation experiments conducted on medial collateral ligaments (MCLs) taken from the literature, whereas the predictive capability of the model was assessed by simulating experimental tests different from those used for the parameter estimation. In particular, creep tests at different maximum stresses have been successfully simulated. The proposed nonlinear viscoelastic model is able to predict the time-dependent response of ligaments described in experimental works (Bonifasi-Lista et al., 2005, J. Orthopaed. Res., 23, pp. 67-76; Hingorani et al., 2004, Ann. Biomed. Eng., 32, pp. 306-312; Provenzano et al., 2001, Ann. Biomed. Eng., 29, pp. 908-214; Weiss et al., 2002, J. Biomech., 35, pp. 943-950). In particular, the nonlinear viscoelastic response which implies different relaxation rates for different applied strains, as well as different creep rates for different applied stresses and direction-dependent relaxation behavior, can be described.     

The role of viscoelasticity of collagen fibers in articular cartilage: theory and numerical formulation

The relative importance of fluid-dependent and fluid-independent transient mechanical behavior in articular cartilage was examined for tensile and unconfined compression testing using a fibril reinforced model. The collagen matrix of articular cartilage was modeled as viscoelastic using a quasi-linear viscoelastic formulation with strain-dependent elastic modulus, while the proteoglycan matrix was considered as linearly elastic. The collagen viscoelastic properties were obtained by fitting experimental data from a tensile test. These properties were used to investigate unconfined compression testing, and the sensitivity of the properties was also explored. It was predicted that the stress relaxation observed in tensile tests was not caused by fluid pressurization at the macroscopic level. A multi-step tensile stress relaxation test could be approximated using a hereditary integral in which the elastic fibrillar modulus was taken to be a linear function of the fibrillar strain. Applying the same formulation to the radial fibers in unconfined compression, stress relaxation could not be simulated if fluid pressurization were absent. Collagen viscoelasticity was found to slightly weaken fluid pressurization in unconfined compression, and this effect was relatively more significant at moderate strain rates. Therefore, collagen viscoelasticity appears to play an import role in articular cartilage in tensile testing, while fluid pressurization dominates the transient mechanical behavior in compression. Collagen viscoelasticity plays a minor role in the mechanical response of cartilage in unconfined compression if significant fluid flow is present.     

Biphasic poroviscoelastic simulation of the unconfined compression of articular cartilage: II--Effect of variable strain rates

This study investigated the abilities of the linear biphasic poroviscoelastic (BPVE) model and the linear biphasic poroelastic (BPE) model to simulate the effect of variable ramp strain rates on the unconfined compression stress relaxation response of articular cartilage. Curve fitting of experimental data showed that the BPVE model was able to successfully account for the ramp strain rate-dependent viscoelastic behavior of articular cartilage under unconfined compression, while the BPE model was able to account for the complete viscoelastic response at a slow strain rate, but only the long-term viscoelastic response at faster strain rates. We concluded that the short-term viscoelastic behavior of articular cartilage, when subjected to a fast ramp strain rate, is primarily governed by a fluid flow-independent (intrinsic) viscoelastic mechanism, whereas the long-term viscoelastic behavior is governed by a fluid flow-dependent (biphasic) viscoelastic mechanism. Furthermore, a linear viscoelastic representation of the solid stress was found to be a valid model assumption for the simulation of ramp strain rate-dependent relaxation behaviors of articular cartilage within the range of ramp strain rates investigated.     

Prediction of internal pressure profile of compression bandages using stress relaxation parameters

The efficacy of compression therapy using compression bandages is highly dependent on the level of compression applied and the sustenance of the pressure during the course of treatment. This study attempts to predict the pressure profile generated by compression bandages using constitutive equations describing relaxation behavior of viscoelastic materials. It is observed that this pressure profile is highly correlated with the stress relaxation behavior of the bandage. To model the pressure profile, the stress relaxation behavior of compression bandages was studied and modeled using three mechanical models: the Maxwell model, the standard linear solid model and the two-component Maxwell model with a nonlinear spring. It was observed that the models with more component values explained the experimental relaxation curves better. The parameters used for modelling relaxation behavior were used to describe the pressure profile, which is significantly dependent on the longitudinal stress relaxation behavior of the bandage, using the modified Laplace's law equation. This approach thus helps in evaluating the bandage performance with time during compression therapy as novel wound care management.     

Viscoelastic properties of passive skeletal muscle in compression—Cyclic behaviour

Skeletal muscle relaxation behaviour in compression has been previously reported, but the anisotropic behaviour at higher loading rates remains poorly understood. In this paper, uniaxial unconfined cyclic compression tests were performed on fresh porcine muscle samples at various fibre orientations to determine muscle viscoelastic behaviour. Mean compression level of 25% was applied and cycles of 2% and 10% amplitude were performed at 0.2–80 Hz. Under cycles of low frequency and amplitude, linear viscoelastic cyclic relaxation was observed. Fibre/cross-fibre results were qualitatively similar, but the cross-fibre direction was stiffer (ratio of 1.2). In higher amplitude tests nonlinear viscoelastic behaviour with a frequency dependent increase in the stress cycles amplitude was found (factor of 4.1 from 0.2 to 80 Hz).The predictive capability of an anisotropic quasi-linear viscoelastic model previously fitted to stress-relaxation data from similar tissue samples was investigated. Good qualitative results were obtained for low amplitude cycles but differences were observed in the stress cycle amplitudes (errors of 7.5% and 31.8%, respectively, in the fibre/cross-fibre directions). At higher amplitudes significant qualitative and quantitative differences were evident. A nonlinear model formulation was therefore developed which provided a good fit and predictions to high amplitude low frequency cyclic tests performed in the fibre/cross-fibre directions. However, this model gave a poorer fit to high frequency cyclic tests and to relaxation tests. Neither model adequately predicts the stiffness increase observed at frequencies above 5 Hz.Together with data previously presented, the experimental data presented here provide a unique dataset for validation of future constitutive models for skeletal muscle in compression.     

The biomechanical effects of limb lengthening and botulinum toxin type A on rabbit tendon

Numerous studies have examined the effects of distraction osteogenesis (DO) on bone, but relatively fewer have explored muscle adaptation, and even less have addressed the concomitant alterations that occur in the tendon. The purpose herein was to characterize the biomechanical properties of normal and elongated rabbit (N=20) tendons with and without prophylactic botulinum toxin type A (BTX-A) treatment. Elastic and viscoelastic properties of Achilles and Tibialis anterior (TA) tendons were evaluated through pull to failure and stress relaxation tests.All TA tendons displayed nonlinear viscoelastic responses that were strain dependent. A power law formulation was used to model tendon viscoelastic responses and tendon elastic responses were fit with a microstructural model. Distraction-elongated tendons displayed increases in compliance and stress relaxation rates over undistracted tendons; BTX-A administration offset this result. The elastic moduli of distraction-lengthened TA tendons were diminished (p=0.010) when distraction was combined with gastrocnemius (GA) BTX-A administration, elastic moduli were further decreased (p=0.004) and distraction following TA BTX-A administration resulted in TA tendons with moduli not different from contralateral control (p>0.05). Compared to contralateral control, distraction and GA BTX-A administration displayed shortened toe regions, (p=0.031 and 0.038, respectively), while tendons receiving BTX-A in the TA had no differences in the toe region (p>0.05). Ultimate tensile stress was unaltered by DO, but stress at the transition from the toe to the linear region of the stress–stretch curve was diminished in all distraction-elongated TA tendons (p<0.05). The data suggest that prophylactic BTX-A treatment to the TA protects some tendon biomechanical properties.     

Stress relaxation behaviors of articular cartilages in porcine temporomandibular joint

In this study, we tested the compressive stress relaxation behaviors of the mandibular condylar and temporal cartilages in the porcine temporomandibular joint (TMJ). The aim was to determine the quantitative and qualitative similarities and differences of compressive stress relaxation behaviors between the two cartilages. Ten porcine TMJs were used; the articular surface was divided into 5 regions: anterior, central, posterior, lateral and medial. Compressive relaxation test was carried out at a strain level of 5% in each region of the two cartilages. The stress relaxation was monitored over a period of 5 min. In all the regions of the two cartilages, the time-dependent stress relaxation curves showed a marked drop in stress within the initial 10 s, which can be fitted by a standard linear viscoelastic model. The instantaneous moduli in the temporal cartilage were dominantly larger than those in the condylar cartilage, while the condylar cartilage had slightly larger relaxation moduli than the temporal cartilage except for the medial region. The both cartilages showed the regional differences in the compressive stress relaxation behavior, and in the temporal cartilage the lateral and medial regions revealed the largest values for the instantaneous and relaxation moduli. The present results demonstrate that the viscoelastic properties of compressive stress relaxation in both cartilages are region-specific, which might have an important implication for stress distribution and transmission along with the TMJ disc.     

The effect of auxin on stress relaxation in isolated Avena coleoptiles

In order to characterize further the mechanical properties of coleoptile cell walls, stress relaxation measurements were made on methanol-boiled sections of Avena coleoptiles. Relaxation was measured both in mechanically conditioned specimens and in specimens which had not been previously extended. In both cases the relaxation was proportional to log time. Mechanical conditioning increased the relaxation modules and decreased the relative rate of relaxation. In contrast, pretreatment of the live coleoptiles with indoleacetic acid reduced the relaxation modulus and the absolute rate of relaxation but did not affect the relative rate of relaxation. Essentially similar pictures of the mechanical properties of coleoptile walls are obtained from stress relaxation and creep tests; the wall behaves as a nonlinear viscoelastic material.     

A nonlinear model of passive muscle viscosity

The material properties of passive skeletal muscle are critical to proper function and are frequently a target for therapeutic and interventional strategies. Investigations into the passive viscoelasticity of muscle have primarily focused on characterizing the elastic behavior, largely neglecting the viscous component. However, viscosity is a sizeable contributor to muscle stress and extensibility during passive stretch and thus there is a need for characterization of the viscous as well as the elastic components of muscle viscoelasticity. Single mouse muscle fibers were subjected to incremental stress relaxation tests to characterize the dependence of passive muscle stress on time, strain and strain rate. A model was then developed to describe fiber viscoelasticity incorporating the observed nonlinearities. The results of this model were compared with two commonly used linear viscoelastic models in their ability to represent fiber stress relaxation and strain rate sensitivity. The viscous component of mouse muscle fiber stress was not linear as is typically assumed, but rather a more complex function of time, strain and strain rate. The model developed here, which incorporates these nonlinearities, was better able to represent the stress relaxation behavior of fibers under the conditions tested than commonly used models with linear viscosity. It presents a new tool to investigate the changes in muscle viscous stresses with age, injury and disuse.     

The role of mucus gel viscosity, spinnability, and adhesive properties in clearance by simulated cough

We investigated the role of the viscoelastic and adhesive properties of mucus gel simulants on the clearance of mucus by simulated cough. Mucus-like gels with widely varying viscoelastic properties were prepared from polysaccharides crosslinked with sodium borate. Cough was simulated by opening a solenoid valve connecting a model trachea to a pressurized tank. The clearance of gels lining the model trachea was quantified by observing marker particle transport. Viscosity elastic modulus, relaxation time and yield stress were measured with a steady-shear viscoelastometer. Spinnability (thread formation) was determined with a filancemeter. Adhesivity (surface tension) was measured by the platinum ring technique. The viscoelastic and adhesive properties of the mucus gel simulants spanned the ranges observed for bronchial secretions from patients with COPD. The relationship between simulated cough clearance and the viscoelastic and adhesive properties of the gels was analyzed by stepwise linear regression of the non-zero data matrix. The major independent variable relating to clearance was viscosity. Secondary, but highly significant dependences, were also found for spinnability and adhesivity. Elastic modulus, relaxation time and yield stress had no independent effect on cough clearance over the investigated range. The results indicate that, in the absence of airway surface liquid, cough-type clearance relates primarily with mucus gel viscosity. For a given viscosity, clearance is also impaired by spinnability, i.e. the capacity of the mucus to form threads. At constant viscosity and spinnability, clearance is further impaired by an increase in the adhesivity of the mucus. The negative dependence of each of these physical factors can be rationalized in terms of their inhibitory effect on wave formation in the mucus lining layer during high velocity airflow interaction.     

Non-linear viscoelastic behavior of abdominal aortic aneurysm thrombus

The objective of this work was to determine the linear and non-linear viscoelastic behavior of abdominal aortic aneurysm thrombus and to study the changes in mechanical properties throughout the thickness of the thrombus. Samples are gathered from thrombi of seven patients. Linear viscoelastic data from oscillatory shear experiments show that the change of properties throughout the thrombus is different for each thrombus. Furthermore the variations found within one thrombus are of the same order of magnitude as the variation between patients. To study the non-linear regime, stress relaxation experiments are performed. To describe the phenomena observed experimentally, a non-linear multimode model is presented. The parameters for this model are obtained by fitting this model successfully to the experiments. The model cannot only describe the average stress response for all thrombus samples but also the highest and lowest stress responses. To determine the influence on the wall stress of the behavior observed the model proposed needs to implemented in the finite element wall stress analysis.     

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.     

Nonlinear viscoelasticity of cheese

Cheese viscoelasticity is commonly measured using steady uniaxial compression, steady uniaxial extension, shear and compressive creep, and stress relaxation in shear and compression. Viscoelastic properties for many cheeses have also been studied using small amplitude oscillatory shear (SAOS). However, there is little on the measurement of nonlinear viscoelastic properties. The large deformation test usually conducted on cheese to study nonlinear viscoelasticity is uniaxial compression, but this test hardly departs from linear behavior. In this work, the nonlinear viscoelasticity of four cheese varieties was studied using large amplitude oscillatory shear (LAOS). A sliding plate rheometer incorporating a shear stress transducer was used. The data were evaluated using spectral analysis, and results are presented mainly in the form of shear stress versus shear rate loops.     

A comparison of cartilage stress-relaxation models in unconfined compression: QLV and stretched exponential in combination with fluid flow

Cartilage exhibits nonlinear viscoelastic behaviour. Various models have been proposed to explain cartilage stress relaxation, but it is unclear whether explicit modelling of fluid flow in unconfined compression is needed. This study compared Fung's quasi-linear viscoelastic (QLV) model with a stretched-exponential model of cartilage stress relaxation and examined each of these models both alone and in combination with a fluid-flow model in unconfined compression. Cartilage explants were harvested from bovine calf patellofemoral joints and equilibrated in tissue culture for 5 days before stress-relaxation testing in unconfined compression at 5% nominal strain. The stretched exponential models fit as well as the QLV models. Furthermore, the average stretched exponential relaxation time determined by this model lies within the range of experimentally measured relaxation times for extracted proteoglycan aggregates, consistent with the hypothesis that the stretched exponential model represents polymeric mechanisms of cartilage viscoelasticity.     

Contribution of the extracellular matrix to the viscoelastic behavior of the urinary bladder wall

We previously reported that when the stress relaxation response of urinary bladder wall (UBW) tissue was analyzed using a single continuous reduced relaxation function (RRF), we observed non-uniformly distributed, time-dependent residuals (Ann Biomed Eng 32(10):1409-1419, 2004). We concluded that the single relaxation spectrum was inadequate and that a new viscoelastic model for bladder wall was necessary. In the present study, we report a new approach composed of independent RRFs for smooth muscle and the extracellular matrix components (ECM), connected through a stress-dependent recruitment function. In order to determine the RRF for the ECM component, biaxial stress relaxation experiments were first performed on decellularized extracellular matrix network of the bladder obtained from normal and spinal cord injured rats. While it was assumed that smooth muscle followed a single spectrum RRF, modeling the UBW ECM required a dual-Gaussian spectrum. Experimental results revealed that the ECM stress relaxation response was insensitive to the initial stress level. Thus, the average ECM RRF parameters were determined by fitting the average stress relaxation data. The resulting stress relaxation behavior of whole bladder tissue was modeled by combining the ECM RRF with the RRF for the smooth muscle component using an exponential recruitment function representing the recruitment of collagen fibers at higher stress levels. In summary, the present study demonstrated, for the first time, that stress relaxation response of bladder tissue can be better modeled when divided into the contributions of the extracellular matrix and smooth muscle components. This modeling approach is suitable for prediction of mechanical behaviors of the urinary bladder and other organs that exhibit rapid tissue remodeling (i.e., smooth muscle hypertrophy and altered ECM synthesis) under various pathological conditions.