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.     

Anisotropic mechanical properties of magnetically aligned fibrin gels measured by magnetic resonance elastography

The anisotropic mechanical properties of magnetically aligned fibrin gels were measured by magnetic resonance elastography (MRE) and by a standard mechanical test: unconfined compression. Soft anisotropic biomaterials are notoriously difficult to characterize, especially in vivo. MRE is well-suited for efficient, non-invasive, and non-destructive assessment of shear modulus. Direction-dependent differences in shear modulus were found to be statistically significant for gels polymerized at magnetic fields of 11.7 and 4.7 T compared to control gels. Mechanical anisotropy was greater in the gels polymerized at the higher magnetic field. These observations were consistent with results from unconfined compression tests. Analysis of confocal microscopy images of gels showed measurable alignment of fibrils in gels polymerized at 11.7 T. This study provides direct, quantitative measurements of the anisotropy in mechanical properties that accompanies fibril alignment in fibrin gels.     

The role of fibril reinforcement in the mechanical behavior of cartilage

Collagen fibril reinforcement was incorporated into a nonlinear poroelastic model for articular cartilage in unconfined compression. It was found that the radial fibrils play a predominant role in the transient mechanical behavior but a less important role in the equilibrium response of cartilage. The radial fibrils are in tension and can be highly stressed during compression, in contrast to low compressive stresses in all directions for the proteoglycan matrix after a small initial compression. The strain dependent fibril stiffening produces strong nonlinear transient response; the fibrils provide extra stiffness to balance a rising fluid pressure and to restrain stress increase in the proteoglycans. The fibril reinforcement, induced by the fluid pressure and flow, also accounts for a complex pattern of strain-magnitude and strain-rate dependence of cartilage stiffness.     

Fibril reinforced poroelastic model predicts specifically mechanical behavior of normal,proteoglycan depleted and collagen degraded articular cartilage

Degradation of collagen network and proteoglycan (PG) macromolecules are signs of articular cartilage degeneration. These changes impair cartilage mechanical function. Effects of collagen degradation and PG depletion on the time-dependent mechanical behavior of cartilage are different. In this study, numerical analyses, which take the compression-tension nonlinearity of the tissue into account, were carried out using a fibril reinforced poroelastic finite element model. The study aimed at improving our understanding of the stress-relaxation behavior of normal and degenerated cartilage in unconfined compression. PG and collagen degradations were simulated by decreasing the Young's modulus of the drained porous (nonfibrillar) matrix and the fibril network, respectively. Numerical analyses were compared to results from experimental tests with chondroitinase ABC (PG depletion) or collagenase (collagen degradation) digested samples. Fibril reinforced poroelastic model predicted the experimental behavior of cartilage after chondroitinase ABC digestion by a major decrease of the drained porous matrix modulus (-64+/-28%) and a minor decrease of the fibril network modulus (-11+/-9%). After collagenase digestion, in contrast, the numerical analyses predicted the experimental behavior of cartilage by a major decrease of the fibril network modulus (-69+/-5%) and a decrease of the drained porous matrix modulus (-44+/-18%). The reduction of the drained porous matrix modulus after collagenase digestion was consistent with the microscopically observed secondary PG loss from the tissue. The present results indicate that the fibril reinforced poroelastic model is able to predict specifically characteristic alterations in the stress-relaxation behavior of cartilage after enzymatic modifications of the tissue. We conclude that the compression-tension nonlinearity of the tissue is needed to capture realistically the mechanical behavior of normal and degenerated articular cartilage.     

Whey Protein/Arabic Gum Gels Formed by Chemical or Physical Gelation Process

Whey proteins (WP) gelation process with addition of Arabic gum (AG) was studied. Two different driving processes were employed to induce gelation: (1) heating of 12% whey protein isolate (WPI) solutions (w/w) or (2) acidification of previous thermal denatured WPI solutions (5% w/w) with glucono-δ-lactone (GDL). Protein concentrations were different because they were minimal to form gel in these two processes, but denaturation conditions were the same (90 °C/30 min). Water-holding capacity and mechanical properties of the gels were evaluated. The BST equation was used to evaluate the nonlinear part of the stress–strain data. Cold-set gels were weaker than heat-set gels at the pH range near the isoelectric point (pI) of the main whey proteins, but heated gels were more deformable (did not exhibit rupture point) and showed greater elasticity modulus. However, gels formed by heating far from the pI (pH 6.7 or 3.5) showed more fragile structure, indicating that, in these mixed gels, there are prevailing biopolymers interactions. Cold-set and heat-set gels at pH near or below the WP pI showed strain-weakening behavior, but heated gels at neutral pH showed strong strain-hardening behavior. Such results suggest that differences in stress–strain curve at the nonlinear part of the data could be correlated to structure particularities obtained from different gelation processes.     

Collagen fiber orientation and geometry effects on the mechanical properties of secondary osteons.

The effects of collagen fiber orientation and osteon geometry on the mechanical properties of secondary osteons under axial compression/tension and combined loadings (compression, bending and torsion) were investigated using a composite-beam finite-element model. Three cross-sectional shapes of secondary osteons were studied to show the effect of geometry. The results of stiffness are presented using the tension and compression properties for each lamella. The model shows that the mechanical properties of osteons are enhanced in bending and torsion when collagen fibers are oriented within 30 degrees of the loading axis. Osteons with alternating lamellar orientation are not well adapted to resist torsional moments, but alternate collagen fiber orientation has virtually no effect on the bending stiffness of osteons. Fiber orientation affects the mechanical properties less significantly when osteons are non-circular. Collagen fiber orientation and osteon geometry interact to determine the mechanical behavior of the osteon, and may act in a compensatory manner in the adaptive process.     

Rheological properties and mechanical stability of new gel-entrapment systems applied in bioreactors

The mechanical stability of gels applied for entrapment and retention of biocatalysts in bioreactors is of crucial importance for successful scale-up applications. Gel abrasion in agitated reactors will depend on liquid shear, bubble shear, and wall shear, as well as collisions between the gel particles. As a simplified standardized model system, abrasion of gel beads was studied in 1-m-high bubble columns with controlled aeration, and quantified by measuring the loss of gel material into solution. Gel beads were also taken out to measure stress-strain response during controlled compression. More general rheological properties of different gels were studied by applying a variety of regimes of controlled compression of standardized gel cylinders: Gel strength was measured by recording the fracture properties and the Young's modulus. Viscoelastic properties were revealed by recording creep during compression as well as recovery after compression. Oscillation tests up to 1000 cyclic compressions were applied to compare the fatigue of different gels. Results obtained for Ca-alginate gels, gels of chemically modified polyvinyl alcohol with stilbazolium groups (PVA-SbQ) as well as mixed gels of Ca-alginate and PVA-SbQ are compared with previously published data for kappa-carrageenan, agar, and polyethylene glycol (PEG) gels. It is concluded that material fatigue rather than mechanical properties such as stiffness or fracture stress should be considered when selecting a suitable gel material on the basis of abrasion resistance. The very soft and superelastic PVA-SbQ gel showed no significant fatigue in mechanical tests and no abrasion was detected in the standardized model system used. Ca-alginate gels, however, showed severe irreversible changes due to fatigue at oscillating loads and creep at constant load. Due to their similarities with kappa-carrageenan gels in mechanical tests, it is likely that Ca-alginate would also be sensitive to abrasion. Mixed gels of Ca-alginate and PVA-SbQ represent a complex system with intermediate properties, showing significant fatigue and creep, but elastic properties from the PVA-SbQ gel make it less sensitive than the pure Ca-alginate gel.     

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.     

A fibril reinforced nonhomogeneous poroelastic model for articular cartilage: inhomogeneous response in unconfined compression

The depth dependence of material properties of articular cartilage, known as the zonal differences, is incorporated into a nonlinear fibril-reinforced poroelastic model developed previously in order to explore the significance of material heterogeneity in the mechanical behavior of cartilage. The material variations proposed are based on extensive observations. The collagen fibrils are modeled as a distinct constituent which reinforces the other two constituents representing proteoglycans and water. The Young's modulus and Poisson's ratio of the drained nonfibrillar matrix are so determined that the aggregate compressive modulus for confined geometry fits the experimental data. Three nonlinear factors are considered, i.e. the effect of finite deformation, the dependence of permeability on dilatation and the fibril stiffening with its tensile strain. Solutions are extracted using a finite element procedure to simulate unconfined compression tests. The features of the model are then demonstrated with an emphasis on the results obtainable only with a nonhomogeneous model, showing reasonable agreement with experiments. The model suggests mechanical behaviors significantly different from those revealed by homogeneous models: not only the depth variations of the strains which are expected by qualitative analyses, but also, for instance, the relaxation-time dependence of the axial strain which is normally not expected in a relaxation test. Therefore, such a nonhomogeneous model is necessary for better understanding of the mechanical behavior of cartilage.     

In Vitro Fracture Testing of Submicron Diameter Collagen Fibril Specimens

Mechanical testing of collagenous tissues at different length scales will provide improved understanding of the mechanical behavior of structures such as skin, tendon, and bone, and also guide the development of multiscale mechanical models. Using a microelectromechanical-systems (MEMS) platform, stress-strain response curves up to failure of type I collagen fibril specimens isolated from the dermis of sea cucumbers were obtained in vitro. A majority of the fibril specimens showed brittle fracture. Some displayed linear behavior up to failure, while others displayed some nonlinearity. The fibril specimens showed an elastic modulus of 470 ± 410 MPa, a fracture strength of 230 ± 160 MPa, and a fracture strain of 80% ± 44%. The fibril specimens displayed significantly lower elastic modulus in vitro than previously measured in air. Fracture strength/strain obtained in vitro and in air are both significantly larger than those obtained in vacuo, indicating that the difference arises from the lack of intrafibrillar water molecules produced by vacuum drying. Furthermore, fracture strength/strain of fibril specimens were different from those reported for collagenous tissues of higher hierarchical levels, indicating the importance of obtaining these properties at the fibrillar level for multiscale modeling.     

Effects of decorin proteoglycan on fibrillogenesis,ultrastructure, and mechanics of type I collagen gels

The proteoglycan decorin is known to affect both the fibrillogenesis and the resulting ultrastructure of in vitro polymerized collagen gels. However, little is known about its effects on mechanical properties. In this study, 3D collagen gels were polymerized into tensile test specimens in the presence of decorin proteoglycan, decorin core protein, or dermatan sulfate (DS). Collagen fibrillogenesis, ultrastructure, and mechanical properties were then quantified using a turbidity assay, 2 forms of microscopy (SEM and confocal), and tensile testing. The presence of decorin proteoglycan or core protein decreased the rate and ultimate turbidity during fibrillogenesis and decreased the number of fibril aggregates (fibers) compared to control gels. The addition of decorin and core protein increased the linear modulus by a factor of 2 compared to controls, while the addition of DS reduced the linear modulus by a factor of 3. Adding decorin after fibrillogenesis had no effect, suggesting that decorin must be present during fibrillogenesis to increase the mechanical properties of the resulting gels. These results show that the inclusion of decorin proteoglycan during fibrillogenesis of type I collagen increases the modulus and tensile strength of resulting collagen gels. The increase in mechanical properties when polymerization occurs in the presence of the decorin proteoglycan is due to a reduction in the aggregation of fibrils into larger order structures such as fibers and fiber bundles.     

Increasing strain and strain rate strengthen transient stiffness but weaken the response to subsequent compression for articular cartilage in unconfined compression

Strain amplitude and strain rate dependent nonlinear behavior and load-induced mechanical property alterations of full-thickness bovine articular cartilage attached to bone were investigated in unconfined compression. A sequence of test compressions of finite deformation (ranging from 0.9% to 34.5% nominal strain) was performed at strain rates ranging from approximately 0.053%/s to 5.8%/s. Peak and equilibrium loads were analyzed to determine strain amplitude and strain rate dependence of linear versus nonlinear responses. The test protocol was designed to reveal changes in mechanical properties due to these finite deformations by interspersing small-amplitude witness ramps of approximately 1.1% deformation and approximately 0.44%/s strain rate between the test ramps ("witness" meaning to assess any mechanical property changes). We found that peak loads displayed high nonlinearity, stiffening with both increasing compression amplitude and more so with increasing strain rate. The response to witness ramps suggested that mechanical weakening occurred when compression amplitude reached 1.9-2.9% strain and beyond, and that weakening was much more significant at higher strain rate. These findings delineate regimes of linear versus nonlinear behavior of cartilage, and indicate the types of loads which can cause mechanical property alterations. Biological implications of this study are that strain amplitude and strain rate dependent stiffening may be essential to bear physiological loads and to protect cells and matrix from mechanical damage. Structural changes reflected by mechanical weakening at small compression could also initiate remodeling or disease processes.     

The asymmetry of transient response in compression versus release for cartilage in unconfined compression

Observations in compression tests of articular cartilage have revealed unequal load increments for compression and release of the same amplitude applied to a disk with an identical previously imposed compression (in equilibrium). The mechanism of this asymmetric transient response is investigated here using a nonlinear fibril-reinforced model. It is found that the asymmetry is predominantly produced by the fibril stiffening with its tensile strain. In addition, allowing the hydraulic permeability to decrease significantly with compressive dilatation of cartilage increases the transient fibril strain, resulting in a stronger asymmetry. Large deformation also enhances the asymmetry as a consequence of stronger fibril stiffening.     

Nonlinear mechanical property of tracheal cartilage: a theoretical and experimental study

BACKGROUND: Despite being the stiffest airway of the bronchial tree, the trachea undergoes significant deformation due to intrathoracic pressure during breathing. The mechanical properties of the trachea affect the flow in the airway and may contribute to the biological function of the lung. METHOD: A Fung-type strain energy density function was used to investigate the nonlinear mechanical behavior of tracheal cartilage. A bending test on pig tracheal cartilage was performed and a mathematical model for analyzing the deformation of tracheal cartilage was developed. The constants included in the strain energy density function were determined by fitting the experimental data. RESULT: The experimental data show that tracheal cartilage is a nonlinear material displaying higher strength in compression than in tension. When the compression forces varied from -0.02 to -0.03N and from -0.03 to -0.04N, the deformation ratios were 11.03+/-2.18% and 7.27+/-1.59%, respectively. Both were much smaller than the deformation ratios (20.01+/-4.49%) under tension forces of 0.02 to 0.01N. The Fung-type strain energy density function can capture this nonlinear behavior very well, whilst the linear stress-strain relation cannot. It underestimates the stability of trachea by exaggerating the displacement in compression. This study may improve our understanding of the nonlinear behavior of tracheal cartilage and it may be useful for the future study on tracheal collapse behavior under physiological and pathological conditions.     

Mechanical properties of human patellar tendon at the hierarchical levels of tendon and fibril

Tendons are strong hierarchical structures, but how tensile forces are transmitted between different levels remains incompletely understood. Collagen fibrils are thought to be primary determinants of whole tendon properties, and therefore we hypothesized that the whole human patellar tendon and its distinct collagen fibrils would display similar mechanical properties. Human patellar tendons (n = 5) were mechanically tested in vivo by ultrasonography. Biopsies were obtained from each tendon, and individual collagen fibrils were dissected and tested mechanically by atomic force microscopy. The Young's modulus was 2.0 ± 0.5 GPa, and the toe region reached 3.3 ± 1.9% strain in whole patellar tendons. Based on dry cross-sectional area, the Young's modulus of isolated collagen fibrils was 2.8 ± 0.3 GPa, and the toe region reached 0.86 ± 0.08% strain. The measured fibril modulus was insufficient to account for the modulus of the tendon in vivo when fibril content in the tendon was accounted for. Thus, our original hypothesis was not supported, although the in vitro fibril modulus corresponded well with reported in vitro tendon values. This correspondence together with the fibril modulus not being greater than that of tendon supports that fibrillar rather than interfibrillar properties govern the subfailure tendon response, making the fibrillar level a meaningful target of intervention. The lower modulus found in vitro suggests a possible adverse effect of removing the tissue from its natural environment. In addition to the primary work comparing the two hierarchical levels, we also verified the existence of viscoelastic behavior in isolated human collagen fibrils.     

Photo-cross-linking of type I collagen gels in the presence of smooth muscle cells: mechanical properties,cell viability,and function

The effectiveness of photomediated cross-linking of type I collagen gels in the presence of rat aortic smooth muscle cells (RASMC) as a method to enhance gel mechanical properties while retaining native collagen triple helical structure and maintaining high cell viability was investigated. Collagen was chemically modified to incorporate an acrylate moiety. Collagen methacrylamide was cast into gels in the presence of a photoinitiator along with RASMC. The gels were cross-linked using visible light irradiation. Neither acrylate modification nor the cross-linking reaction altered collagen triple helical content. The cross-linking reaction, however, moved the denaturation temperature beyond the physiologic range. A twelve-fold increase in shear modulus was observed after cross-linking. Cell viability in the range of 70% (n = 4, p > 0.05) was observed in the photo-cross-linked gels. Moreover the cells were able to contract the cross-linked gel in a manner commensurate with that observed for natural type I collagen. Methacrylate-mediated photo-cross-linking is a facile route to improve mechanical properties of collagen gels in the presence of cells while maintaining high cell viability. This enhances the potential for type I collagen gels to be used as scaffolds for tissue engineering.     

Viscoelastic behaviour and failure of bovine cancellous bone under constant strain rate

A programme has been established to characterize the long-term behaviour of cancellous bone. Fresh bovine cancellous specimens of dimensions 10 x 10 x 10 mm3 and 10 x 40 x 3.6 mm3 were manufactured and used within the testing programme. Results published in the literature indicate that the long-term behaviour of cancellous bone is well described by a power law, which is a very similar response of typical polymers. So far, dynamic mechanical tests (DMA) in three-point bending, under frequencies between 0.01 and 100 Hz at room temperature, confirmed the published results in a qualitative way. Nevertheless, the measured dimensionless damping, tan delta, was slightly higher than the values reported in the literature for the compact bone. The relaxation curves were obtained from dynamic tests and confirmed that bone relaxation modulus can be described by a power law function of time. Tests under constant compression strain rate were performed at four different strain rates: 0.15/s, 0.015/s, 0.0015/s and 0.00015/s and strain rate dependent behaviour was observed. An average elastic bending modulus of 300 MPa was obtained.     

Fitted hyperelastic parameters for Human brain tissue from reported tension,compression, and shear tests

The mechanical properties of human brain tissue are the subject of interest because of their use in understanding brain trauma and in developing therapeutic treatments and procedures. To represent the behavior of the tissue, we have developed hyperelastic mechanical models whose parameters are fitted in accordance with experimental test results. However, most studies available in the literature have fitted parameters with data of a single type of loading, such as tension, compression, or shear. Recently, Jin et al. (Journal of Biomechanics 46:2795−2801, 2013) reported data from ex vivo tests of human brain tissue under tension, compression, and shear loading using four strain rates and four different brain regions. However, they do not report parameters of energy functions that can be readily used in finite element simulations. To represent the tissue behavior for the quasi-static loading conditions, we aimed to determine the best fit of the hyperelastic parameters of the hyperfoam, Ogden, and polynomial strain energy functions available in ABAQUS for the low strain rate data, while simultaneously considering all three loading modes. We used an optimization process conducted in MATLAB, calling iteratively three finite element models developed in ABAQUS that represent the three loadings. Results showed a relatively good fit to experimental data in all loading modes using two terms in the energy functions. Values for the shear modulus obtained in this analysis (897−1653 Pa) are in the range of those presented in other studies. These energy-function parameters can be used in brain tissue simulations using finite element models.