The prediction of stress-relaxation of ligaments and tendons using the quasi-linear viscoelastic model

Recent studies have questioned the ability of the quasi-linear viscoelastic (QLV) model to predict stresses and strains in response to loading conditions other than those used to fit the model. The objective of this study was to evaluate the ability of several models in the literature to predict the elastic stress response of ligament and tendon at strain levels higher than the levels used to fit the model. The constitutive models were then used to evaluate the ability of the QLV model to predict the overall stress response during stress relaxation. The models expressing stress as an exponential function of strain significantly overestimated stress when used at higher strain levels. The polynomial formulation of the Mooney–Rivlin model more accurately predicted the stress–strain behavior of ligament and tendon. The results demonstrate that the ability of the QLV model to accurately predict the stress-relaxation response is dependent in part on the accuracy of the function used to model the elastic response of the soft tissue.     

An approach to determine pressure profile generated by compression bandage using quasi-linear viscoelastic model

Understanding the stress relaxation behavior of the compression bandage could be very useful in determining the behavior of the interface pressure exerted by the bandage on a limb during the course of the compression treatment. There has been no comprehensive study in the literature to investigate the pressure profile (interface pressure with time) generated by a compression bandage when applied at different levels of strain. The present study attempts to describe the pressure profile, with the use of a quasi-linear viscoelastic model, generated by a compression bandage during compression therapy. The quasi-linear viscoelastic (QLV) theory proposed by Fung (Fung, 1972, "Stress Strain History Relations of Soft Tissues in Simple Elongation," Biomechanics: Its Foundations and Objectives, Y. C. Fung, N. Perrone, and M. Anliker, eds., Prentice-Hall, Englewood Cliffs, NJ, pp. 181-207). was used to model the nonlinear time- and history-dependent relaxation behavior of the bandage using the ramp strain approach. The regression analysis was done to find the correlation between the pressure profile and the relaxation behavior of the bandage. The parameters of the QLV model, describing the relaxation behavior of the bandage, were used to determine the pressure profile generated by the bandage at different levels of strain. The relaxation behaviors of the bandage at different levels of strain were well described by the QLV model parameters. A high correlation coefficient (nearly 0.98) shows a good correlation of the pressure profile with the stress relaxation behavior of the bandage.The prediction of the pressure profile using the QLV model parameters were in agreement with the experimental data. The pressure profile generated by a compression bandage could be predicted using the QLV model describing the nonlinear relaxation behavior of the bandage. This new application of the QLV theory helps in evaluating the bandage performance during compression therapy as scientific wound care management.     

The effect of overshooting the target strain on estimating viscoelastic properties from stress relaxation experiments

BACKGROUND: Tendon's mechanical behaviors have frequently been quantified using the quasi-linear viscoelastic (QLV) model. The QLV parameters are typically estimated by fitting the model to a single-step stress relaxation experiment. Unfortunately, overshoot of the target strain occurs to some degree in most experiments. This has never been formally investigated even though failing to measure, minimize, or compensate for overshoot may cause large errors in the estimation of parameters. Therefore, the objective of this study was to investigate the effect of overshoot on the estimation of QLV parameters. METHOD OF APPROACH: A simulated experiment was first performed to quantify the effect of different amounts of overshoot on the estimated QLV parameters. Experimental data from tendon was then used to determine if the errors associated with overshoot could be reduced when a direct fit is used (i.e., the actual strain history was used in the curve fit). RESULTS: We found that both the elastic and viscous QLV parameters were incorrectly estimated if overshoot was not properly accounted for in the fit. Furthermore, the errors associated with overshoot were partially reduced when overshoot was accounted for using a direct fit. CONCLUSIONS: A slow ramp rate is recommended to limit the amount of overshoot and a direct fit is recommended to limit the errors associated with overshoot, although other approaches such as adjusting the control system to limit overshoot could also be utilized.     

Viscous elements have little impact on measured passive length-tension properties of human gastrocnemius muscle-tendon units in vivo

Several studies have measured the elastic properties of a single human muscle-tendon unit in vivo. However the viscoelastic behavior of single human muscles has not been characterized. In this study, we adapted QLV theory to model the viscoelastic behavior of human gastrocnemius muscle-tendon units in vivo. We also determined the influence of viscoelasticity on passive length-tension properties of human gastrocnemius muscle-tendon units. Eight subjects participated in the experiment, which consisted of two parts. First, the stress relaxation response of human gastrocnemius muscle-tendon units was determined at a range of knee and ankle angles. Subsequently, passive ankle torque and ankle angle were collected during cyclic dorsiflexion and plantarflexion at a range of knee angles. Viscous parameters were determined by fitting the stress relaxation experiment data with a two-term exponential function, and elastic parameters were estimated by fitting the QLV model and viscous parameters to the cyclic experiment data. The model fitted the experimental data well at slow speeds (RMSE: 1.7 ± 0.5N) and at fast speeds (RMSE: 1.9 ± 0.2N). Muscle-tendon units demonstrated a large amount of stress relaxation. Nonetheless, viscoelastic passive length-tension curves estimated with the QLV model were similar to elastic passive length-tension curves obtained using a model that ignored viscosity. There was little difference in the elastic passive length-tension curves at different loading rates. We conclude that (a) the QLV model can be used to quantify viscoelastic behaviors of relaxed human gastrocnemius muscle-tendon units in vivo, and (b) over the range of velocities we examined, the velocity of loading has little effect on the passive length-tension properties of human gastrocnemius muscle-tendon units.     

An improved method to analyze the stress relaxation of ligaments following a finite ramp time based on the quasi-linear viscoelastic theory

The quasi-linear viscoelastic (QLV) theory proposed by Fung (1972) has been frequently used to model the nonlinear time- and history-dependent viscoelastic behavior of many soft tissues. It is common to use five constants to describe the instantaneous elastic response (constants A and B) and reduced relaxation function (constants C, tau 1, and tau 2) on experiments with finite ramp times followed by stress relaxation to equilibrium. However, a limitation is that the theory is based on a step change in strain which is not possible to perform experimentally. Accounting for this limitation may result in regression algorithms that converge poorly and yield nonunique solutions with highly variable constants, especially for long ramp times (Kwan et al. 1993). The goal of the present study was to introduce an improved approach to obtain the constants for QLV theory that converges to a unique solution with minimal variability. Six goat femur-medial collateral ligament-tibia complexes were subjected to a uniaxial tension test (ramp time of 18.4 s) followed by one hour of stress relaxation. The convoluted QLV constitutive equation was simultaneously curve-fit to the ramping and relaxation portions of the data (r2 > 0.99). Confidence intervals of the constants were generated from a bootstrapping analysis and revealed that constants were distributed within 1% of their median values. For validation, the determined constants were used to predict peak stresses from a separate cyclic stress relaxation test with averaged errors across all specimens measuring less than 6.3 +/- 6.0% of the experimental values. For comparison, an analysis that assumed an instantaneous ramp time was also performed and the constants obtained for the two approaches were compared. Significant differences were observed for constants B, C, tau 1, and tau 2, with tau 1 differing by an order of magnitude. By taking into account the ramping phase of the experiment, the approach allows for viscoelastic properties to be determined independent of the strain rate applied. Thus, the results obtained from different laboratories and from different tissues may be compared.     

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.     

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.     

Mechanical responses of tendons to repeated extensions and wait periods

Tendon specimens were repeatedly extended to peak strains of either 2, 3, 4, or 6 percent. During the three 1800 s (30 min.) periods of cyclic extension, the peak loads relaxed with decreases in hysteresis and increases in slack strain. During the two 1800 s wait periods of no extension, the specimens recovered with increases in peak load and hysteresis and decreases in slack strain. However, the recovery during the wait periods was eradicated in the first few subsequent extensions and the relaxation continued as if there were no 1800 s wait periods. Stress-strain responses were well fit with power relations.     

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.     

The shear mechanical properties of diabetic and non-diabetic plantar soft tissue

Changes in the plantar soft tissue shear properties may contribute to ulceration in diabetic patients, however, little is known about these shear parameters. This study examines the elastic and viscoelastic shear behavior of both diabetic and non-diabetic plantar tissue. Previously compression tested plantar tissue specimens (n=54) at six relevant plantar locations (hallux, first, third, and fifth metatarsal heads, lateral midfoot, and calcaneus) from four cadaveric diabetic feet and five non-diabetic feet were utilized. Per in vivo data (i.e., combined deformation patterns of compression followed by shear), an initial static compressive strain (36-38%) was applied to the tissue followed by target shear strains of 50% and 85% of initial thickness. Triangle waves were used to quantify elastic parameters at both strain levels and a stress relaxation test (0.25 s ramp and 300 s hold) was used to quantify the viscoelastic parameters at the upper strain level. Several differences were found between test groups including a 52-62% increase in peak shear stress, a 63% increase in toe shear modulus, a 47% increase in final shear modulus, and a 67% increase in middle slope magnitude (sharper drop in relaxation) in the diabetic tissue. Beyond a 54% greater peak compressive stress in the third metatarsal compared to the lateral midfoot, there were no differences in shear properties between plantar locations. Notably, this study demonstrates that plantar soft tissue with diabetes is stiffer than healthy tissue, thereby compromising its ability to dissipate shear stresses borne by the foot that may increase ulceration risk.     

Fractional order viscoelasticity of the aortic valve cusp: an alternative to quasilinear viscoelasticity

BACKGROUND: Quasilinear viscoelasticity (QLV) theory has been widely and successfully used to describe the time-dependent response of connective tissues. Difficulties remain, however, particularly in material parameter estimation and sensitivities. In this study, we introduce a new alternative: the fractional order viscoelasticity (FOV) theory, which uses a fractional order integral to describe the relaxation response. FOV implies a fractal-like tissue structure, reflecting the hierarchical arrangement of collagenous tissues. METHOD OF APPROACH: A one-dimensional (I-D) FOV reduced relaxation function was developed, replacing the QLV "box-spectrum" function with a fractional relaxation function. A direct-fit, global optimization method was used to estimate material parameters from stress relaxation tests on aortic valve tissue. RESULTS: We found that for the aortic heart valve, FOV had similar accuracy and better parameter sensitivity than QLV, particularly for the long time constant (tau2). The mean (n = 5) fractional order was 0.29, indicating that the viscoelastic response of the tissue was strongly fractal-like. RESULTS SUMMARY: mean QLV parameters were C = 0.079, tau1 = 0.004, tau2 = 76, and mean FOV parameters were beta = 0.29, tau = 0.076, and rho = 1.84. CONCLUSIONS: FOV can provide valuable new insights into tissue viscoelastic behavior Determining the fractional order can provide a new and sensitive quantitative measure for tissue comparison.     

Methods for quasi-linear viscoelastic modeling of soft tissue: application to incremental stress-relaxation experiments

The quasi-linear viscoelastic (QLV) model was applied to incremental stress-relaxation tests and an expression for the stress was derived for each step. This expression was used to compare two methods for normalizing stress data prior to estimating QLV parameters. The first and commonly used normalization method was shown to be strain-dependent. Thus, a second normalization method was proposed and shown to be strain-independent and more sensitive to QLV time constants. These analytical results agreed with representative tendon data. Therefore, this method for normalizing stress data was proposed for future studies of incremental stress-relaxation, or whenever comparing stress-relaxation at different strains.     

The relation between collagen fibril kinematics and mechanical properties in the mitral valve anterior leaflet

We have recently demonstrated that the mitral valve anterior leaflet (MVAL) exhibited minimal hysteresis, no strain rate sensitivity, stress relaxation but not creep (Grashow et al., 2006, Ann Biomed Eng., 34(2), pp. 315-325; Grashow et al., 2006, Ann Biomed. Eng., 34(10), pp. 1509-1518). However, the underlying structural basis for this unique quasi-elastic mechanical behavior is presently unknown. As collagen is the major structural component of the MVAL, we investigated the relation between collagen fibril kinematics (rotation and stretch) and tissue-level mechanical properties in the MVAL under biaxial loading using small angle X-ray scattering. A novel device was developed and utilized to perform simultaneous measurements of tissue level forces and strain under a planar biaxial loading state. Collagen fibril D-period strain (epsilonD) and the fibrillar angular distribution were measured under equibiaxial tension, creep, and stress relaxation to a peak tension of 90 N/m. Results indicated that, under equibiaxial tension, collagen fibril straining did not initiate until the end of the nonlinear region of the tissue-level stress-strain curve. At higher tissue tension levels, epsilonD increased linearly with increasing tension. Changes in the angular distribution of the collagen fibrils mainly occurred in the tissue toe region. Using epsilonD, the tangent modulus of collagen fibrils was estimated to be 95.5+/-25.5 MPa, which was approximately 27 times higher than the tissue tensile tangent modulus of 3.58+/-1.83 MPa. In creep tests performed at 90 N/m equibiaxial tension for 60 min, both tissue strain and epsilonD remained constant with no observable changes over the test length. In contrast, in stress relaxation tests performed for 90 min epsilonD was found to rapidly decrease in the first 10 min followed by a slower decay rate for the remainder of the test. Using a single exponential model, the time constant for the reduction in collagen fibril strain was 8.3 min, which was smaller than the tissue-level stress relaxation time constants of 22.0 and 16.9 min in the circumferential and radial directions, respectively. Moreover, there was no change in the fibril angular distribution under both creep and stress relaxation over the test period. Our results suggest that (1) the MVAL collagen fibrils do not exhibit intrinsic viscoelastic behavior, (2) tissue relaxation results from the removal of stress from the fibrils, possibly by a slipping mechanism modulated by noncollagenous components (e.g. proteoglycans), and (3) the lack of creep but the occurrence of stress relaxation suggests a "load-locking" behavior under maintained loading conditions. These unique mechanical characteristics are likely necessary for normal valvular function.     

Method for characterizing viscoelasticity of human gluteal tissue

Characterizing compressive transient large deformation properties of biological tissue is becoming increasingly important in impact biomechanics and rehabilitation engineering, which includes devices interfacing with the human body and virtual surgical guidance simulation. Individual mechanical in vivo behaviour, specifically of human gluteal adipose and passive skeletal muscle tissue compressed with finite strain, has, however, been sparsely characterised. Employing a combined experimental and numerical approach, a method is presented to investigate the time-dependent properties of in vivo gluteal adipose and passive skeletal muscle tissue. Specifically, displacement-controlled ramp-and-hold indentation relaxation tests were performed and documented with magnetic resonance imaging. A time domain quasi-linear viscoelasticity (QLV) formulation with Prony series valid for finite strains was used in conjunction with a hyperelastic model formulation for soft tissue constitutive model parameter identification and calibration of the relaxation test data. A finite element model of the indentation region was employed. Strong non-linear elastic but linear viscoelastic tissue material behaviour at finite strains was apparent for both adipose and passive skeletal muscle mechanical properties with orthogonal skin and transversal muscle fibre loading. Using a force-equilibrium assumption, the employed material model was well suited to fit the experimental data and derive viscoelastic model parameters by inverse finite element parameter estimation. An individual characterisation of in vivo gluteal adipose and muscle tissue could thus be established. Initial shear moduli were calculated from the long-term parameters for human gluteal skin/fat: G(∞,S/F)=1850 Pa and for cross-fibre gluteal muscle tissue: G(∞,M)=881 Pa. Instantaneous shear moduli were found at the employed ramp speed: G(0,S/F)=1920 Pa and G(0,M)=1032 Pa.     

In situ estimation of tendon material properties: Differences between muscles of the feline hindlimb

Recent experiments to characterize the short-range stiffness (SRS)–force relationship in several cat hindlimb muscles suggested that the there are differences in the tendon elastic moduli across muscles [Cui, L., Perreault, E.J., Maas, H., Sandercock, T.G., 2008. Modeling short-range stiffness of feline lower hindlimb muscles. J. Biomech. 41 (9), 1945–1952.]. Those conclusions were inferred from whole muscle experiments and a computational model of SRS. The present study sought to directly measure tendon elasticity, the material property most relevant to SRS, during physiological loading to confirm the previous modeling results. Measurements were made from the medial gastrocnemius (MG), tibialis anterior (TA) and extensor digitorum longus (EDL) muscles during loading. For the latter, the model indicated a substantially different elastic modulus than for MG and TA. For each muscle, the stress–strain relationship of the external tendon was measured in situ during the loading phase of isometric contractions conducted at optimum length. Young's moduli were assessed at equal strain levels (1%, 2% and 3%), as well as at peak strain. The stress–strain relationship was significantly different between EDL and MG/TA, but not between MG and TA. EDL had a more apparent toe region (i.e., lower Young's modulus at 1% strain), followed by a more rapid increase in the slope of the stress–strain curve (i.e., higher Young's modulus at 2% and 3% strain). Young's modulus at peak strain also was significantly higher in EDL compared to MG/TA, whereas no significant difference was found between MG and TA. These results indicate that during natural loading, tendon Young's moduli can vary considerably across muscles. This creates challenges to estimating muscle behavior in biomechanical models for which direct measures of tendon properties are not available.     

A simplified approach to quasi-linear viscoelastic modeling

The fitting of quasi-linear viscoelastic (QLV) constitutive models to material data often involves somewhat cumbersome numerical convolution. A new approach to treating quasi-linearity in 1-D is described and applied to characterize the behavior of reconstituted collagen. This approach is based on a new principle for including nonlinearity and requires considerably less computation than other comparable models for both model calibration and response prediction, especially for smoothly applied stretching. Additionally, the approach allows relaxation to adapt with the strain history. The modeling approach is demonstrated through tests on pure reconstituted collagen. Sequences of "ramp-and-hold" stretching tests were applied to rectangular collagen specimens. The relaxation force data from the "hold" was used to calibrate a new "adaptive QLV model" and several models from literature, and the force data from the "ramp" was used to check the accuracy of model predictions. Additionally, the ability of the models to predict the force response on a reloading of the specimen was assessed. The "adaptive QLV model" based on this new approach predicts collagen behavior comparably to or better than existing models, with much less computation.     

A visco-hyperelastic constitutive approach for modeling polyvinyl alcohol sponge

This study proposes the quasi-linear viscoelastic (QLV) model to characterize the time dependent mechanical behavior of poly(vinyl alcohol) (PVA) sponges. The PVA sponges have implications in many viscoelastic soft tissues, including cartilage, liver, and kidney as an implant. However, a critical barrier to the use of the PVA sponge as tissue replacement material is a lack of sufficient study on its viscoelastic mechanical properties. In this study, the nonlinear mechanical behavior of a fabricated PVA sponge is investigated experimentally and computationally using relaxation and stress failure tests as well as finite element (FE) modeling. Hyperelastic strain energy density functions, such as Yeoh and Neo-Hookean, are used to capture the mechanical behavior of PVA sponge at ramp part, and viscoelastic model is used to describe the viscose behavior at hold part. Hyperelastic material constants are obtained and their general prediction ability is verified using FE simulations of PVA tensile experiments. The results of relaxation and stress failure tests revealed that Yeoh material model can define the mechanical behavior of PVA sponge properly compared with Neo-Hookean one. FE modeling results are also affirmed the appropriateness of Yeoh model to characterize the mechanical behavior of PVA sponge. Thus, the Yeoh model can be used in future biomechanical simulations of the spongy biomaterials. These results can be utilized to understand the viscoelastic behavior of PVA sponges and has implications for tissue engineering as scaffold.     

Mechanical characterization of electrospun polycaprolactone (PCL): a potential scaffold for tissue engineering

This paper investigates the mechanical behavior of electrospun polycaprolactone (PCL) under tensile loading. PCL in bulk form degrades slowly and is biocompatible, two properties that make it a viable option for tissue engineering applications in biomedicine. Of particular interest is the use of electrospun PCL tubes as scaffolds for tissue engineered blood vessel implants. Stress relaxation and tensile tests have been conducted with specimens at room temperature (21 degrees C) and 37 degrees C. Additionally, to probe the effects of moisture on mechanical behavior, specimens were tested either dry (in air) or submerged in water. In general, the electrospun PCL was found to exhibit rate dependence, as well as some dependence on the test temperature and on whether the sample was wet or dry. Two different models were investigated to describe the experimentally observed material behavior. The models used were Fung's theory of quasilinear viscoelasticity (QLV) and the eight-chain model developed for rubber elastomers by Arruda and Boyce (1993, "A Three-Dimensional Constitutive Model for the Large Stretch Behavior of Rubber Elastic Materials," J. Mech. Phys. Solids, 41(2), pp. 389-412). The implementation and fitting results, as well as the advantages and disadvantages of each model, are presented. In general, it was found that the QLV theory provided a better fit.