Confocal arthroscopy-based patient-specific constitutive models of cartilaginous tissues—II: prediction of reaction force history of meniscal cartilage specimens

The theoretical framework developed in a companion paper (Part I) is used to derive estimates of mechanical response of two meniscal cartilage specimens. The previously developed framework consisted of a constitutive model capable of incorporating confocal image-derived tissue microstructural data. In the present paper (Part II) fibre and matrix constitutive parameters are first estimated from mechanical testing of a batch of specimens similar to, but independent from those under consideration. Image analysis techniques which allow estimation of tissue microstructural parameters form confocal images are presented. The constitutive model and image-derived structural parameters are then used to predict the reaction force history of the two meniscal specimens subjected to partially confined compression. The predictions are made on the basis of the specimens' individual structural condition as assessed by confocal microscopy and involve no tuning of material parameters. Although the model does not reproduce all features of the experimental curves, as an unfitted estimate of mechanical response the prediction is quite accurate. In light of the obtained results it is judged that more general non-invasive estimation of tissue mechanical properties is possible using the developed framework.     

Confocal arthroscopy-based patient-specific constitutive models of cartilaginous tissues--I: development of a microstructural model

Current development of a laser scanning confocal arthroscope within our school will enable 3D microscopic imaging of joint tissues in vivo. Such an instrument could be useful, for example, in assessing the microstructural condition of the living tissues without physical biopsy. It is envisaged also that linked to a suitable microstructural constitutive formulation, such imaging could allow non-invasive patient-specific estimation of tissue mechanical performance. Such a procedure could have applications in surgical planning and simulation, and assessment of engineered tissue replacements, where tissue biopsy is unacceptable. In this first of two papers the development of a suitable constitutive framework for generating such estimates is reported. A microstructure-based constitutive formulation for cartilaginous tissues is presented. The model extends existing fibre composite-type models and accounts for strain-rate sensitivity of the tissue mechanical response through incorporation of a viscoelastic fibre phase. Importantly, the model is constructed so as to allow direct incorporation of structural data from confocal images. A finite element implementation of the formulation suitable for incorporation within commercial codes is also presented.     

Inferring spatial variations of microstructural properties from macroscopic mechanical response

Disease alters tissue microstructure, which in turn affects the macroscopic mechanical properties of tissue. In elasticity imaging, the macroscopic response is measured and is used to infer the spatial distribution of the elastic constitutive parameters. When an empirical constitutive model is used, these parameters cannot be linked to the microstructure. However, when the constitutive model is derived from a microstructural representation of the material, it allows for the possibility of inferring the local averages of the spatial distribution of the microstructural parameters. This idea forms the basis of this study. In particular, we first derive a constitutive model by homogenizing the mechanical response of a network of elastic, tortuous fibers. Thereafter, we use this model in an inverse problem to determine the spatial distribution of the microstructural parameters. We solve the inverse problem as a constrained minimization problem and develop efficient methods for solving it. We apply these methods to displacement fields obtained by deforming gelatin–agar co-gels and determine the spatial distribution of agar concentration and fiber tortuosity, thereby demonstrating that it is possible to image local averages of microstructural parameters from macroscopic measurements of deformation.     

Modeling the effect of collagen fibril alignment on ligament mechanical behavior

Ligament mechanical behavior is primarily regulated by fibrous networks of type I collagen. Although these fibrous networks are typically highly aligned, healthy and injured ligament can also exhibit disorganized collagen architecture. The objective of this study was to determine whether variations in the collagen fibril network between neighboring ligaments can predict observed differences in mechanical behavior. Ligament specimens from two regions of bovine fetlock joints, which either exhibited highly aligned or disorganized collagen fibril networks, were mechanically tested in uniaxial tension. Confocal microscopy and FiberFit software were used to quantify the collagen fibril dispersion and mean fibril orientation in the mechanically tested specimens. These two structural parameters served as inputs into an established hyperelastic constitutive model that accounts for a continuous distribution of planar fibril orientations. The ability of the model to predict differences in the mechanical behavior between neighboring ligaments was tested by (1) curve fitting the model parameters to the stress response of the ligament with highly aligned fibrils and then (2) using this model to predict the stress response of the ligament with disorganized fibrils by only changing the parameter values for fibril dispersion and mean fibril orientation. This study found that when using parameter values for fibril dispersion and mean fibril orientation based on confocal imaging data, the model strongly predicted the average stress response of ligaments with disorganized fibrils (\(R^{2}=0.97\)); however, the model only successfully predicted the individual stress response of ligaments with disorganized fibrils in half the specimens tested. Model predictions became worse when parameters for fibril dispersion and mean fibril orientation were not based on confocal imaging data. These findings emphasize the importance of collagen fibril alignment in ligament mechanics and help advance a mechanistic understanding of fibrillar networks in healthy and injured ligament.     

A constitutive model for vascular tissue that integrates fibril, fiber and continuum levels with application to the isotropic and passive properties of the infrarenal aorta

A fundamental understanding of the mechanical properties of the extracellular matrix (ECM) is critically important to quantify the amount of macroscopic stress and/or strain transmitted to the cellular level of vascular tissue. Structural constitutive models integrate histological and mechanical information, and hence, allocate stress and strain to the different microstructural components of the vascular wall. The present work proposes a novel multi-scale structural constitutive model of passive vascular tissue, where collagen fibers are assembled by proteoglycan (PG) cross-linked collagen fibrils and reinforce an otherwise isotropic matrix material. Multiplicative kinematics account for the straightening and stretching of collagen fibrils, and an orientation density function captures the spatial organization of collagen fibers in the tissue. Mechanical and structural assumptions at the collagen fibril level define a piece-wise analytical stress-stretch response of collagen fibers, which in turn is integrated over the unit sphere to constitute the tissue's macroscopic mechanical properties. The proposed model displays the salient macroscopic features of vascular tissue, and employs the material and structural parameters of clear physical meaning. Likewise, the constitutive concept renders a highly efficient multi-scale structural approach that allows for the numerical analysis at the organ level. Model parameters were estimated from isotropic mean-population data of the normal and aneurysmatic aortic wall and used to predict in-vivo stress states of patient-specific vascular geometries, thought to demonstrate the robustness of the particular Finite Element (FE) implementation. The collagen fibril level of the multi-scale constitutive formulation provided an interface to integrate vascular wall biology and to account for collagen turnover.     

Experimental-numerical analysis of minipig's multi-rooted teeth

The paper pertains to the analysis of the biomechanical behaviour of the periodontal ligament (PDL) by using a combined experimental and numerical approach. Experimental analysis provides information about a two-rooted pig premolar tooth in its socket with regard to morphological configuration and deformational response. The numerical analysis developed for the present investigation adopts a specific anisotropic hyperelastic formulation, accounting for tissue structural arrangement. The parameters to be adopted for the PDL constitutive model are evaluated with reference to data deducted from experimental in vitro tests on different specimens taken from literature. According to morphometric data relieved, solid models are provided as basis for the development of numerical models that adopt the constitutive formulation proposed. A reciprocal validation of experimental and numerical data allows for the evaluation of reliability of results obtained. The work is intended as preliminary investigation to study the correlation between mechanical status of PDL and induction to cellular activity in orthodontic treatments.     

Constitutive modeling of mouse carotid arteries using experimentally measured microstructural parameters

Changes in the local mechanical environment and tissue mechanical properties affect the biological activity of cells and play a key role in a variety of diseases, such as cancer, arthritis, nephropathy, and cardiovascular disease. Constitutive relations have long been used to predict the local mechanical environment within biological tissues and to investigate the relationship between biological responses and mechanical stimuli. Recent constitutive relations for soft tissues consider both material and structural properties by incorporating parameters that describe microstructural organization, such as fiber distributions, fiber angles, fiber crimping, and constituent volume fractions. The recently developed technique of imaging the microstructure of a single artery as it undergoes multiple deformations provides quantitative structural data that can reduce the number of estimated parameters by using parameters that are truly experimentally intractable. Here, we employed nonlinear multiphoton microscopy to quantify collagen fiber organization in mouse carotid arteries and incorporated measured fiber distribution data into structurally motivated constitutive relations. Microscopy results demonstrate that collagen fibers deform in an affine manner over physiologically relevant deformations. The incorporation of measured fiber angle distributions into constitutive relations improves the model's predictive accuracy and does not significantly reduce the goodness of fit. The use of measured structural parameters rather than estimated structural parameters promises to improve the predictive capabilities of the local mechanical environment, and to extend the utility of intravital imaging methods for estimating the mechanical behavior of tissues using in situ structural information.     

Mechanical properties and microstructure of intraluminal thrombus from abdominal aortic aneurysm.

Accurate estimation of the wall stress distribution in an abdominal aortic aneurysm (AAA) may prove clinically useful by predicting when a particular aneurysm will rupture. Appropriate constitutive models for both the wall and the intraluminal thrombus (ILT) found in most AAA are necessary for this task. The purpose of this work was to determine the mechanical properties of ILT within AAA and to derive a more suitable constitutive model for this material. Uniaxial tensile testing was carried out on 50 specimens, including 14 longitudinally oriented and 14 circumferentially oriented specimens from the luminal region of the ILT, and 11 longitudinally oriented and 11 circumferentially oriented specimens from the medial region. A two-parameter, large-strain, hyperelastic constitutive model was developed and used to fit the uniaxial tensile testing data for determination of the material parameters. Maximum stiffness and strength were also determined from the data for each specimen. Scanning electron microscopy (SEM) was conducted to study the regional microstructural difference. Our results indicate that the microstructure of ILT differs between the luminal, medial, and abluminal regions, with the luminal region stronger and stiffer than the medial region. In all cases, the constitutive model fit the experimental data very well (R2>0.98). No significant difference was found for either of the two material parameters between longitudinal and circumferential directions, but a significant difference in material parameters, stiffness, and strength between the laminal and medial regions was determined (p<0.01). Therefore, our results suggest that ILT is an inhomogeneous and possibly isotropic material. The two-parameter, hyperelastic, isotropic, incompressible material model derived here for ILT can be easily incorporated into finite element models for simulation of wall stress distribution in AAA.     

Identifying a minimal rheological configuration: a tool for effective and efficient constitutive modeling of soft tissues

We describe a modeling methodology intended as a preliminary step in the identification of appropriate constitutive frameworks for the time-dependent response of biological tissues. The modeling approach comprises a customizable rheological network of viscous and elastic elements governed by user-defined 1D constitutive relationships. The model parameters are identified by iterative nonlinear optimization, minimizing the error between experimental and model-predicted structural (load-displacement) tissue response under a specific mode of deformation. We demonstrate the use of this methodology by determining the minimal rheological arrangement, constitutive relationships, and model parameters for the structural response of various soft tissues, including ex vivo perfused porcine liver in indentation, ex vivo porcine brain cortical tissue in indentation, and ex vivo human cervical tissue in unconfined compression. Our results indicate that the identified rheological configurations provide good agreement with experimental data, including multiple constant strain rate load/unload tests and stress relaxation tests. Our experience suggests that the described modeling framework is an efficient tool for exploring a wide array of constitutive relationships and rheological arrangements, which can subsequently serve as a basis for 3D constitutive model development and finite-element implementations. The proposed approach can also be employed as a self-contained tool to obtain simplified 1D phenomenological models of the structural response of biological tissue to single-axis manipulations for applications in haptic technologies.     

An approach to the mechanical constitutive modelling of arterial tissue based on homogenization and optimization

This paper is concerned with characterizing the quasistatic mechanical behaviour of arterial tissue undergoing finite deformation through hyperelastic constitutive functions. Commonly the parameters of constitutive functions are established by a process of optimization based on experimental data. Instead we construct a finite element model of a representative volume element of the material and compute its homogenized response to a range of deformations. These data are then used to provide objective functions for optimizing the parameters of two analytical models from the literature.     

Experimentally validated microstructural 3D constitutive model of coronary arterial media

Accurate modeling of arterial response to physiological or pathological loads may shed light on the processes leading to initiation and progression of a number of vascular diseases and may serve as a tool for prediction and diagnosis. In this study, a microstructure based hyperelastic constitutive model is developed for passive media of porcine coronary arteries. The most general model contains 12 independent parameters representing the three-dimensional inner fibrous structure of the media and includes the effects of residual stresses and osmotic swelling. Parameter estimation and model validation were based on mechanical data of porcine left anterior descending (LAD) media under radial inflation, axial extension, and twist tests. The results show that a reduced four parameter model is sufficient to reliably predict the passive mechanical properties. These parameters represent the stiffness and the helical orientation of each lamellae fiber and the stiffness of the interlamellar struts interconnecting these lamellae. Other structural features, such as orientational distribution of helical fibers and anisotropy of the interlamellar network, as well as possible transmural distribution of structural features, were found to have little effect on the global media mechanical response. It is shown that the model provides good predictions of the LAD media twist response based on parameters estimated from only biaxial tests of inflation and extension. In addition, good predictive capabilities are demonstrated for the model behavior at high axial stretch ratio based on data of law stretches.     

A microstructural model for the anisotropic drained stiffness of articular cartilage

A constitutive model for articular cartilage is developed to study directional load sharing within the soft biological tissue. Cartilage is idealized as a composite structure whose static mechanical response is dominated by distortion of a sparse fibrous network and by changes in fixed charge density. These histological features of living cartilage are represented in a microstructural analog of the tissue, linking the directionality of mechanical stiffness to the orientation of microstructure. The discretized 'model tissue' is used to define a stiffness tensor relating drained stress and strain over a regime of large deformation. The primary goal of this work was to develop a methodology permitting more complete treatment of anisotropy in the stiffness of cartilage. The results demonstrate that simple oriented microscopic behaviors can combine to produce complicated larger scale response. For the illustrative example of a homogeneous specimen subjected to confined compression, the model predicts a nonlinear anisotropic drained response, with inherent uncertainty at cellular size scales.     

Incorporation of experimentally-derived fiber orientation into a structural constitutive model for planar collagenous tissues

Structural constitutive models integrate information on tissue composition and structure, avoiding ambiguities in material characterization. However, critical structural information (such as fiber orientation) must be modeled using assumed statistical distributions, with the distribution parameters estimated from fits to the mechanical test data. Thus, full realization of structural approaches continues to be limited without direct quantitative structural information for direct implementation or to validate model predictions. In the present study, fiber orientation information obtained using small angle light scattering (SALS) was directly incorporated into a structural constitutive model based on work by Lanir (J. Biomech., v. 16, pp. 1-12, 1983). Demonstration of the model was performed using existing biaxial mechanical and fiber orientation data for native bovine pericardium (Sacks and Chuong, ABME, v.26, pp. 892-902, 1998). The structural constitutive model accurately predicted the complete measured biaxial mechanical response. An important aspect of this approach is that only a single equibiaxial test to determine the effective fiber stress-strain response and the SALS-derived fiber orientation distribution were required to determine the complete planar biaxial mechanical response. Changes in collagen fiber crimp under equibiaxial strain suggest that, at the meso-scale, fiber deformations follow the global tissue strains. This result supports the assumption of affine strain to estimate the fiber strains. However, future evaluations will have to be performed for tissue subjected to a wider range of strain to more fully validate the current approach.     

Quantifying nonlinear anisotropic elastic material properties of biological tissue by use of membrane inflation

Determination of material parameters for soft tissue frequently involves regression of material parameters for nonlinear, anisotropic constitutive models against experimental data from heterogeneous tests. Here, parameter estimation based on membrane inflation is considered. A four parameter nonlinear, anisotropic hyperelastic strain energy function was used to model the material, in which the parameters are cast in terms of key response features. The experiment was simulated using finite element (FE) analysis in order to predict the experimental measurements of pressure versus profile strain. Material parameter regression was automated using inverse FE analysis; parameter values were updated by use of both local and global techniques, and the ability of these techniques to efficiently converge to a best case was examined. This approach provides a framework in which additional experimental data, including surface strain measurements or local structural information, may be incorporated in order to quantify heterogeneous nonlinear material properties.     

A computational model for understanding the micro-mechanics of collagen fiber network in the tunica adventitia

Abdominal aortic aneurysm is a prevalent cardiovascular disease with high mortality rates. The mechanical response of the arterial wall relies on the organizational and structural behavior of its microstructural components, and thus, a detailed understanding of the microscopic mechanical response of the arterial wall layers at loads ranging up to rupture is necessary to improve diagnostic techniques and possibly treatments. Following the common notion that adventitia is the ultimate barrier at loads close to rupture, in the present study, a finite element model of adventitial collagen network was developed to study the mechanical state at the fiber level under uniaxial loading. Image stacks of the rabbit carotid adventitial tissue at rest and under uniaxial tension obtained using multi-photon microscopy were used in this study, as well as the force–displacement curves obtained from previously published experiments. Morphological parameters like fiber orientation distribution, waviness, and volume fraction were extracted for one sample from the confocal image stacks. An inverse random sampling approach combined with a random walk algorithm was employed to reconstruct the collagen network for numerical simulation. The model was then verified using experimental stress–stretch curves. The model shows the remarkable capacity of collagen fibers to uncrimp and reorient in the loading direction. These results further show that at high stretches, collagen network behaves in a highly non-affine manner, which was quantified for each sample. A comprehensive parameter study to understand the relationship between structural parameters and their influence on mechanical behavior is presented. Through this study, the model was used to conclude important structure–function relationships that control the mechanical response. Our results also show that at loads close to rupture, the probability of failure occurring at the fiber level is up to 2%. Uncertainties in usually employed rupture risk indicators and the stochastic nature of the event of rupture combined with limited knowledge on the microscopic determinants motivate the development of such an analysis. Moreover, this study will advance the study of coupling microscopic mechanisms to rupture of the artery as a whole.

     

Effects of fabrication on the mechanics,microstructure and micromechanical environment of small intestinal submucosa scaffolds for vascular tissue engineering

In small intestinal submucosa scaffolds for functional tissue engineering, the impact of scaffold fabrication parameters on success rate may be related to the mechanotransductory properties of the final microstructural organization of collagen fibers. We hypothesized that two fabrication parameters, 1) preservation (P) or removal (R) of a dense collagen layer present in SIS and 2) SIS in a final dehydrated (D) or hydrated (H) state, have an effect on scaffold void area, microstructural anisotropy (fiber alignment) and mechanical anisotropy (global mechanical compliance). We further integrated our experimental measurements in a constitutive model to explore final effects on the micromechanical environment inside the scaffold volume. Our results indicated that PH scaffolds might exhibit recurrent and large force fluctuations between layers (up to 195 pN), while fluctuations in RH scaffolds might be larger (up to 256 pN) but not as recurrent. In contrast, both PD and RD groups were estimated to produce scarcer and smaller fluctuations (not larger than 50 pN). We concluded that the hydration parameter strongly affects the micromechanics of SIS and that an adequate choice of fabrication parameters, assisted by the herein developed method, might leverage the use of SIS for functional tissue engineering applications, where forces at the cellular level are of concern in the guidance of new tissue formation.     

Biaxial tensile testing and constitutive modeling of human supraspinatus tendon

The heterogeneous composition and mechanical properties of the supraspinatus tendon offer an opportunity for studying the structure-function relationships of fibrous musculoskeletal connective tissues. Previous uniaxial testing has demonstrated a correlation between the collagen fiber angle distribution and tendon mechanics in response to tensile loading both parallel and transverse to the tendon longitudinal axis. However, the planar mechanics of the supraspinatus tendon may be more appropriately characterized through biaxial tensile testing, which avoids the limitation of nonphysiologic traction-free boundary conditions present during uniaxial testing. Combined with a structural constitutive model, biaxial testing can help identify the specific structural mechanisms underlying the tendon's two-dimensional mechanical behavior. Therefore, the objective of this study was to evaluate the contribution of collagen fiber organization to the planar tensile mechanics of the human supraspinatus tendon by fitting biaxial tensile data with a structural constitutive model that incorporates a sample-specific angular distribution of nonlinear fibers. Regional samples were tested under several biaxial boundary conditions while simultaneously measuring the collagen fiber orientations via polarized light imaging. The histograms of fiber angles were fit with a von Mises probability distribution and input into a hyperelastic constitutive model incorporating the contributions of the uncrimped fibers. Samples with a wide fiber angle distribution produced greater transverse stresses than more highly aligned samples. The structural model fit the longitudinal stresses well (median R(2) ≥ 0.96) and was validated by successfully predicting the stress response to a mechanical protocol not used for parameter estimation. The transverse stresses were fit less well with greater errors observed for less aligned samples. Sensitivity analyses and relatively affine fiber kinematics suggest that these errors are not due to inaccuracies in measuring the collagen fiber organization. More likely, additional strain energy terms representing fiber-fiber interactions are necessary to provide a closer approximation of the transverse stresses. Nevertheless, this approach demonstrated that the longitudinal tensile mechanics of the supraspinatus tendon are primarily dependent on the moduli, crimp, and angular distribution of its collagen fibers. These results add to the existing knowledge of structure-function relationships in fibrous musculoskeletal tissue, which is valuable for understanding the etiology of degenerative disease, developing effective tissue engineering design strategies, and predicting outcomes of tissue repair.