Poroviscoelastic finite element model including continuous fiber distribution for the simulation of nanoindentation tests on articular cartilage

Effect of osmolarity and displacement rate on cartilage microfracture clusters failure into two regimes

Articular cartilage is a poroviscoelastic (PVE) material with remarkable resistance to fracture and fatigue failure. Cartilage failure mechanisms and material properties that govern failure are incompletely understood. Because cartilage is partially comprised of negatively charged glycosaminoglycans, altering solvent osmolarity can influence PVE relaxations. Therefore, this study aims to use osmolarity as a tool to provide additional data to interpret the role of PVE relaxations and identify cartilage failure regimes. Cartilage fracture was induced using a 100 μm radius spheroconical indenter at controlled displacement rates under three different osmolarity solvents. Secondarily, contact pressure (CP) and strain energy density (SED) were estimated to cluster data into two failure regimes with an expectation maximization algorithm. Critical displacement, critical load, critical time, and critical work to fracture increased with increasing osmolarity at a slow displacement rate whereas no significant effect was observed at a fast displacement rate. Clustering provided two distinct failure regimes, with regime (I) at lower normalized thickness (contact radius divided by sample thickness), and regime (II) at higher normalized thickness. Varied CP and SED in regime (I) suggest that failure in the regime is strain-governed. Constant CP and SED in regime (II) suggests that failure in the regime is dominantly governed by stress. These regimes can be interpreted as ductile versus brittle, or using a pressurized fragmentation interpretation. These findings demonstrated fundamental failure properties and postulate failure regimes for articular cartilage.

A backward pre-stressing algorithm for efficient finite element implementation of in vivo material and geometrical parameters into fibril-reinforced mixture models of articular cartilage

Classical continuum mechanics has been widely used for implementation of the material models of articular cartilage (AC) mainly with the aid of the finite element (FE) method, which, in many cases, considers the stress-free configuration as the initial configuration. On the contrary, the AC experimental tests typically begin with the pre-stressed state of both material and geometrical properties. Indeed, imposing the initial pre-stress onto AC models with the in vivo values as the initial state would result in nonphysiologically expansion of the FE mesh due to the soft nature of AC. This change in the model configuration can also affect the material behavior kinematically in the mixture models of cartilage due to the intrinsic compressibility of the tissue. Although several different fixed-point backward algorithms, as the most straightforward pre-stressing methods, have already been developed to incorporate these initial conditions into FE models iteratively, such methods focused merely on the geometrical parameters, and they omitted the material variations of the anisotropic mixture models of AC. To address this issue, we propose an efficient algorithm generalizing the backward schemes to restore stress-free conditions by optimizing both the involving variables, and we hypothesize that it can affect the results considerably. To this end, a comparative simulation was implemented on an advanced and validated multiphasic model by the new and conventional algorithms. The results are in support of the hypothesis, as in our illustrative general AC model, the material parameters experienced a maximum error of 16% comparing to the initial in vivo data when the older algorithm was employed, and it led to a maximum variation of 44% in the recorded stresses comparing to the results of the new method. We conclude that our methodology enhanced the model fidelity, and it is applicable in most of the existing FE solvers for future mixture studies with accurate stress distributions.

Computational simulation of the multiphasic degeneration of the bone-cartilage unit during osteoarthritis via indentation and unconfined compression tests

It has been experimentally proposed that the discrete regions of articular cartilage, along with different subchondral bone tissues, known as the bone-cartilage unit, are biomechanically altered during osteoarthritis degeneration. However, a computational framework capturing all of the dominant changes in the multiphasic parameters has not yet been developed. This article proposes a new finite element model of the bone-cartilage unit by combining several validated, nonlinear, depth-dependent, fibril-reinforced, and swelling models, which can computationally simulate the variations in the dominant parameters during osteoarthritis degeneration by indentation and unconfined compression tests. The mentioned dominant parameters include the proteoglycan depletion, collagen fibrillar softening, permeability, and fluid fraction increase for approximately non-advanced osteoarthritis. The results depict the importance of subchondral bone tissues in fluid distribution within the bone-cartilage units by decreasing the fluid permeation and pressure (up to a maximum of 100 kPa) during osteoarthritis, supporting the notion that subchondral bones might play a role in the pathogenesis of osteoarthritis. Furthermore, the osteoarthritis composition-based studies shed light on the significant biomechanical role of the calcified cartilage, which experienced a maximum change of 70 kPa in stress, together with relative load contributions of articular cartilage constituents during osteoarthritis, in which the osmotic pressure bore around 70% of the loads after degeneration. To conclude, the new insights provided by the results reveal the significance of the multiphasic osteoarthritis simulation and demonstrate the functionality of the proposed bone-cartilage unit model.

Effect of the anisotropic permeability in the frequency dependent properties of the superficial layer of articular cartilage

Articular cartilage is a tissue of fundamental importance for the mechanics of joints, since it provides a smooth and lubricated surface for the proper transfer of loads. From a mechanical point of view, this tissue is an anisotropic poroviscoelastic material: its characteristics at the macroscopic level depend on the complex microscopic architecture. With the ability to probe the local microscopic features, dynamic nanoindentation test is a powerful tool to investigate cartilage mechanics. In this work we focus on a length scale where the time dependent behaviour is regulated by poroelasticity more than viscoelasticity and we aim to understand the effect of the anisotropic permeability on the mechanics of the superficial layer of the articular cartilage. In a previous work, a finite element model for the dynamic nanoindentation test has been presented. In this work, we improve the model by considering the presence of an anisotropic permeability tensor that depends on the collagen fibers distribution. Our sensitivity analysis highlights that the permeability decreases with increasing indentation, thus making the tissue stiffer than the case of isotropic permeability, when solicited at the same frequency. With this improved model, a revised identification of the mechanical and physical parameters for articular cartilage is provided. To this purpose the model was used to simulate experimental data from tests performed on bovine tissue, giving a better estimation of the anisotropy in the elastic properties. A relation between the identified macroscopic anisotropic permeability properties and the microscopic rearrangement of the fiber/matrix structure during indentation is also provided.

Poroelastic Modeling of Highly Hydrated Collagen Hydrogels: Experimental Results vs. Numerical Simulation With Custom and Commercial Finite Element Solvers

This study presents a comparison between the performances of two Finite Element (FE) solvers for the modeling of the poroelastic behavior of highly hydrated collagen hydrogels. Characterization of collagen hydrogels has been a widespread challenge since this is one of the most used natural biomaterials for Tissue Engineering (TE) applications. V-Biomech® is a free custom FE solver oriented to soft tissue modeling, while Abaqus® is a general-purpose commercial FE package which is widely used for biomechanics computational modeling. Poroelastic simulations with both solvers were compared to two experimental protocols performed by Busby et al. (2013) and Chandran and Barocas (2004), also using different implementations of the frequently used Neo-Hookean hyperelastic model. The average differences between solvers outputs were under 5% throughout the different tests and hydrogel properties. Thus, differences were small enough to be considered negligible and within the variability found experimentally from one sample to another. This work demonstrates that constitutive modeling of soft tissues, such as collagen hydrogels can be achieved with either V-Biomech or Abaqus standard options (without user-subroutine), which is important for the biomechanics and biomaterials research community. V-Biomech has shown its potential for the validation of biomechanical characterization of soft tissues, while Abaqus' versatility is useful for the modeling and analysis of TE applications where other complex phenomena may also need to be captured.

Biological connective tissues exhibit viscoelastic and poroelastic behavior at different frequency regimes: Application to tendon and skin biophysics

In this study, a poroviscoelastic finite element model (FEM) was developed and used in conjunction with an AFM-based wide-bandwidth nanorheology system to predict the frequency-dependent mechanical behavior of tendon and dermis subjected to compression via nanoindentation. The aim was to distinguish between loading rates that are dominated by either poroelasticity, viscoelasticity, or the superposition of these processes. Using spherical probe tips having different radii, the force and tip displacement were measured and the magnitude, E^∗, and phase angle, ϕ, of the dynamic complex modulus were evaluated for mouse supraspinatus tendon and mouse dermis. The peak frequencies of the phase angle were associated with the characteristic time constants of poroelastic and viscoelastic material behavior. The developed FE model could predict the separate poroelastic and viscoelastic responses of these soft tissues over a 4 decade frequency range, showing good agreement with experimental results. We observed that poroelasticity was the dominant energy dissipation mechanism for mouse dermis and supraspinatus tendon at higher indentation frequencies (10² to 10⁴ Hz) whereas viscoelasticity was typically dominant at lower frequencies (<10² Hz). These findings show the underlying mechanical behavior of biological connective tissues and give insight into the role played by these different energy dissipation mechanisms in governing the function of these tissues at nanoscale.

STATEMENT OF SIGNIFICANCE

Soft biological tissues exhibit complex, load- and time-dependent mechanical behavior. Evaluating their mechanical behavior requires sophisticated experimental tools and numerical models that can capture the fundamental mechanisms governing tissue function. Using an Atomic-force-microscopy-based rheology system and finite element models, the roles of the two most dominant time-dependent mechanisms (poroelasticity and viscoelasticity) that govern the dynamic loading behavior of mouse skin and tendon have been investigated. FE models were able to predict and quantify the contribution of each mechanism to the overall dynamic response and confirming the presence of these two distinct mechanisms in the mechanical response. Overall, these results provide novel insight into the viscoelastic and poroelastic properties of mouse skin and tendon and promote better understanding of the underlying origins of each mechanism.

Indentation mapping revealed poroelastic, but not viscoelastic, properties spanning native zonal articular cartilage

Osteoarthrosis is a debilitating disease affecting millions, yet engineering materials for cartilage regeneration has proven difficult because of the complex microstructure of this tissue. Articular cartilage, like many biological tissues, produces a time-dependent response to mechanical load that is critical to cell's physiological function in part due to solid and fluid phase interactions and property variations across multiple length scales. Recreating the time-dependent strain and fluid flow may be critical for successfully engineering replacement tissues but thus far has largely been neglected. Here, microindentation is used to accomplish three objectives: (1) quantify a material's time-dependent mechanical response, (2) map material properties at a cellular relevant length scale throughout zonal articular cartilage and (3) elucidate the underlying viscoelastic, poroelastic, and nonlinear poroelastic causes of deformation in articular cartilage. Untreated and trypsin-treated cartilage was sectioned perpendicular to the articular surface and indentation was used to evaluate properties throughout zonal cartilage on the cut surface. The experimental results demonstrated that within all cartilage zones, the mechanical response was well represented by a model assuming nonlinear biphasic behavior and did not follow conventional viscoelastic or linear poroelastic models. Additionally, 10% (w/w) agarose was tested and, as anticipated, behaved as a linear poroelastic material. The approach outlined here provides a method, applicable to many tissues and biomaterials, which reveals and quantifies the underlying causes of time-dependent deformation, elucidates key aspects of material structure and function, and that can be used to provide important inputs for computational models and targets for tissue engineering.

STATEMENT OF SIGNIFICANCE

Elucidating the time-dependent mechanical behavior of cartilage, and other biological materials, is critical to adequately recapitulate native mechanosensory cues for cells. We used microindentation to map the time-dependent properties of untreated and trypsin treated cartilage throughout each cartilage zone. Unlike conventional approaches that combine viscoelastic and poroelastic behaviors into a single framework, we deconvoluted the mechanical response into separate contributions to time-dependent behavior. Poroelastic effects in all cartilage zones dominated the time-dependent behavior of articular cartilage, and a model that incorporates tension-compression nonlinearity best represented cartilage mechanical behavior. These results can be used to assess the success of regeneration and repair approaches, as design targets for tissue engineering, and for development of accurate computational models.

On mechanical mechanism of damage evolution in articular cartilage

Superficial lesions of cartilage are the direct indication of osteoarthritis. To investigate the mechanical mechanism of cartilage with micro-defect under external loading, a new plain strain numerical model with micro-defect was proposed and damage evolution progression in cartilage over time has been simulated, the parameter were studied including load style, velocity of load and degree of damage. The new model consists of the hierarchical structure of cartilage and depth-dependent arched fibers. The numerical results have shown that not only damage of the cartilage altered the distribution of the stress but also matrix and fiber had distinct roles in affecting cartilage damage, and damage in either matrix or fiber could promote each other. It has been found that the superficial cracks in cartilage spread preferentially along the tangent direction of the fibers. It is the arched distribution form of fibers that affects the crack spread of cartilage, which has been verified by experiment. During the process of damage evolution, its extension direction and velocity varied constantly with the damage degree. The rolling load could cause larger stress and strain than sliding load. Strain values of the matrix initially increased and then decreased gradually with the increase of velocity, and velocity had a greater effect on matrix than fibers. Damage increased steadily before reaching 50%, sharply within 50 to 85%, and smoothly and slowly after 85%. The finding of the paper may help to understand the mechanical mechanism why the cracks in cartilage spread preferentially along the tangent direction of the fibers.

Combined numerical and experimental biomechanical characterization of soft collagen hydrogel substrate

This work presents a combined experimental-numerical framework for the biomechanical characterization of highly hydrated collagen hydrogels, namely with 0.20, 0.30 and 0.40% (by weight) of collagen concentration. Collagen is the most abundant protein in the extracellular matrix of animals and humans. Its intrinsic biocompatibility makes collagen a promising substrate for embedding cells within a highly hydrated environment mimicking natural soft tissues. Cell behaviour is greatly influenced by the mechanical properties of the surrounding matrix, but the biomechanical characterization of collagen hydrogels has been challenging up to now, since they present non-linear poro-viscoelastic properties. Combining the stiffness outcomes from rheological experiments with relevant literature data on collagen permeability, poroelastic finite element (FE) models were developed. Comparison between experimental confined compression tests available in the literature and analogous FE stress relaxation curves showed a close agreement throughout the tests. This framework allowed establishing that the dynamic shear modulus of the collagen hydrogels is between 0.0097 ± 0.018 kPa for the 0.20% concentration and 0.0601 ± 0.044 kPa for the 0.40% concentration. The Poisson's ratio values for such conditions lie within the range of 0.495-0.485 for 0.20% and 0.480-0.470 for 0.40%, respectively, showing that rheology is sensitive enough to detect these small changes in collagen concentration and thus allowing to link rheology results with the confined compression tests. In conclusion, this integrated approach allows for accurate constitutive modelling of collagen hydrogels. This framework sets the grounds for the characterization of related hydrogels and to the use of this collagen parameterization in more complex multiscale models.

Determining Tension-Compression Nonlinear Mechanical Properties of Articular Cartilage from Indentation Testing

The indentation test is widely used to determine the in situ biomechanical properties of articular cartilage. The mechanical parameters estimated from the test depend on the constitutive model adopted to analyze the data. Similar to most connective tissues, the solid matrix of cartilage displays different mechanical properties under tension and compression, termed tension-compression nonlinearity (TCN). In this study, cartilage was modeled as a porous elastic material with either a conewise linear elastic matrix with cubic symmetry or a solid matrix reinforced by a continuous fiber distribution. Both models are commonly used to describe the TCN of cartilage. The roles of each mechanical property in determining the indentation response of cartilage were identified by finite element simulation. Under constant loading, the equilibrium deformation of cartilage is mainly dependent on the compressive modulus, while the initial transient creep behavior is largely regulated by the tensile stiffness. More importantly, altering the permeability does not change the shape of the indentation creep curves, but introduces a parallel shift along the horizontal direction on a logarithmic time scale. Based on these findings, a highly efficient curve-fitting algorithm was designed, which can uniquely determine the three major mechanical properties of cartilage (compressive modulus, tensile modulus, and permeability) from a single indentation test. The new technique was tested on adult bovine knee cartilage and compared with results from the classic biphasic linear elastic curve-fitting program.

A quantitative interpretation of the response of articular cartilage to atomic force microscopy-based dynamic nanoindentation tests

In this paper, a quantitative interpretation for atomic force microscopy-based dynamic nanoindentation (AFM-DN) tests on the superficial layers of bovine articular cartilage (AC) is provided. The relevant constitutive parameters of the tissue are estimated by fitting experimental results with a finite element model in the frequency domain. Such model comprises a poroelastic stress-strain relationship for a fibril reinforced tissue constitution, assuming a continuous distribution of the collagen network orientations. The identification procedure was first validated using a simplified transversely isotropic constitutive relationship; then, the experimental data were manually fitted by using the continuous distribution fibril model. Tissue permeability is derived from the maximum value of the phase shift between the input harmonic loading and the harmonic tissue response. Tissue parameters related to the stiffness are obtained from the frequency response of the experimental storage modulus and phase shift. With this procedure, an axial to transverse stiffness ratio (anisotropy ratio) of about 0.15 is estimated.