The effect of matrix tension-compression nonlinearity and fixed negative charges on chondrocyte responses in cartilage

Thorough analyses of the mechano-electrochemical interaction between articular cartilage matrix and the chondrocytes are crucial to understanding of the signal transduction mechanisms that modulate the cell metabolic activities and biosynthesis. Attempts have been made to model the chondrocytes embedded in the collagen-proteoglycan extracellular matrix to determine the distribution of local stress-strain field, fluid pressure and the time-dependent deformation of the cell. To date, these models still have not taken into account a remarkable characteristic of the cartilage extracellular matrix given rise from organization of the collagen fiber architecture, now known as the tension-compression nonlinearity (TCN) of the tissue, as well as the effect of negative charges attached to the proteoglycan molecules, and the cell cytoskeleton that interacts with mobile ions in the interstitial fluid to create osmotic and electro-kinetic events in and around the cells. In this study, we proposed a triphasic, multi-scale, finite element model incorporating the Conewise Linear Elasticity that can describe the various known coupled mechanical, electrical and chemical events, while at the same time representing the TCN of the extracellular matrix. The model was employed to perform a detailed analysis of the chondrocytes' deformational and volume responses, and to quantitatively describe the mechano-electrochemical environment of these cells. Such a model describes contributions of the known detailed micro-structural and composition of articular cartilage. Expectedly, results from model simulations showed substantial effects of the matrix TCN on the cell deformational and volume change response. A low compressive Poisson's ratio of the cartilage matrix exhibiting TCN resulted in dramatic recoiling behavior of the tissue under unconfined compression and induced significant volume change in the cell. The fixed charge density of the chondrocyte and the pericellular matrix were also found to play an important role in both the time-dependent and equilibrium deformation of the cell. The pericellular matrix tended to create a uniform osmolarity around the cell and overall amplified the cell volume change. It is concluded that the proposed model can be a useful tool that allows detailed analysis of the mechano-electrochemical interactions between the chondrocytes and its surrounding extracellular matrix, which leads to more quantitative insights in the cell mechano-transduction.     

A Conewise Linear Elasticity mixture model for the analysis of tension-compression nonlinearity in articular cartilage

A biphasic mixture model is developed that can account for the observed tension-compression nonlinearity of cartilage by employing the continuum-based Conewise Linear Elasticity (CLE) model of Curnier et al. (J. Elasticity, 37, 1-38, 1995) to describe the solid phase of the mixture. In this first investigation, the orthotropic octantwise linear elasticity model was reduced to the more specialized case of cubic symmetry, to reduce the number of elastic constants from twelve to four. Confined and unconfined compression stress-relaxation, and torsional shear testing were performed on each of nine bovine humeral head articular cartilage cylindrical plugs from 6 month old calves. Using the CLE model with cubic symmetry, the aggregate modulus in compression and axial permeability were obtained from confined compression (H-A = 0.64 +/- 0.22 MPa, k2 = 3.62 +/- 0.97 x 10(-16) m4/N.s, r2 = 0.95 +/- 0.03), the tensile modulus, compressive Poisson ratio, and radial permeability were obtained from unconfined compression (E+Y = 12.75 +/- 1.56 MPa, v- = 0.03 +/- 0.01, kr = 6.06 +/- 2.10 x 10(-16) m4/N.s, r2 = 0.99 +/- 0.00), and the shear modulus was obtained from torsional shear (mu = 0.17 +/- 0.06 MPa). The model was also employed to predict the interstitial fluid pressure successfully at the center of the cartilage plug in unconfined compression (r2 = 0.98 +/- 0.01). The results of this study demonstrate that the integration of the CLE model with the biphasic mixture theory can provide a model of cartilage that can successfully curve-fit three distinct testing configurations while producing material parameters consistent with previous reports in the literature.     

Experimental verification of the roles of intrinsic matrix viscoelasticity and tension-compression nonlinearity in the biphasic response of cartilage

A biphasic-CLE-QLV model proposed in our recent study [2001, J. Biomech. Eng., 123, pp. 410-417] extended the biphasic theory of Mow et al. [1980, J. Biomech. Eng., 102, pp. 73-84] to include both tension-compression nonlinearity and intrinsic viscoelasticity of the cartilage solid matrix by incorporating it with the conewise linear elasticity (CLE) model [1995, J. Elasticity, 37, pp. 1-38] and the quasi-linear viscoelasticity (QLV) model [Biomechanics: Its foundations and objectives, Prentice Hall, Englewood Cliffs, 1972]. This model demonstrates that a simultaneous prediction of compression and tension experiments of articular cartilage, under stress-relaxation and dynamic loading, can be achieved when properly taking into account both flow-dependent and flow-independent viscoelastic effects, as well as tension-compression nonlinearity. The objective of this study is to directly test this biphasic-CLE-QLV model against experimental data from unconfined compression stress-relaxation tests at slow and fast strain rates as well as dynamic loading. Twelve full-thickness cartilage cylindrical plugs were harvested from six bovine glenohumeral joints and multiple confined and unconfined compression stress-relaxation tests were performed on each specimen. The material properties of specimens were determined by curve-fitting the experimental results from the confined and unconfined compression stress relaxation tests. The findings of this study demonstrate that the biphasic-CLE-QLV model is able to describe the strain-rate-dependent mechanical behaviors of articular cartilage in unconfined compression as attested by good agreements between experimental and theoretical curvefits (r2 = 0.966 +/- 0.032 for testing at slow strain rate; r2 = 0.998 +/- 0.002 for testing at fast strain rate) and predictions of the dynamic response (r2 = 0.91 +/- 0.06). This experimental study also provides supporting evidence for the hypothesis that both tension-compression nonlinearity and intrinsic viscoelasticity of the solid matrix of cartilage are necessary for modeling the transient and equilibrium responses of this tissue in tension and compression. Furthermore, the biphasic-CLE-QLV model can produce better predictions of the dynamic modulus of cartilage in unconfined dynamic compression than the biphasic-CLE and biphasic poroviscoelastic models, indicating that intrinsic viscoelasticity and tension-compression nonlinearity of articular cartilage may play important roles in the load-support mechanism of cartilage under physiologic loading.     

Tensile properties and fiber alignment of human supraspinatus tendon in the transverse direction demonstrate inhomogeneity,nonlinearity, and regional isotropy

A recent study (Lake et al., 2009); reported the properties of human supraspinatus tendon (SST) tested along the predominant fiber direction. The SST was found to have a relatively disperse distribution of collagen fibers, which may represent an adaptation to multiaxial loads imposed by the complex loading environment of the rotator cuff. However, the multiaxial mechanical properties of human SST remain unknown. The objective of this study, therefore, was to evaluate the mechanical properties, fiber alignment, change in alignment with applied load, and structure–function relationships of SST in transverse testing. Samples from six SST locations were tested in uniaxial tension with samples oriented transverse to the tendon long-axis. Polarized light imaging was used to quantify collagen fiber alignment and change in alignment under applied load. The mechanical properties of samples taken near the tendon–bone insertion were much greater on the bursal surface compared to the joint surface (e.g., bursal moduli 15–30 times greater than joint; p<0.001). In fact, the transverse moduli values of the bursal samples were very similar to values obtained from samples tested along the tendon long-axis (Lake et al., 2009). This key and unexpected finding suggests planar mechanical isotropy for bursal surface samples near the insertion, which may be due to complex in vivo loading. Organizationally, fiber distributions became less aligned along the tendon long-axis in the toe-region of the stress–strain response. Alignment changes occurred to a slightly lesser degree in the linear-region, suggesting that movement of collagen fibers may play a role in mechanical nonlinearity. Transverse mechanical properties were significantly correlated with fiber alignment (e.g., for linear-region modulus r_s=0.74, p<0.0001), demonstrating strong structure–function relationships. These results greatly enhance current understanding of the properties of human SST and provide clinicians and scientists with vital information in attempting to treat or replace this complex tissue.     

On the anisotropy and inhomogeneity of permeability in articular cartilage

Articular cartilage is known to be anisotropic and inhomogeneous because of its microstructure. In particular, its elastic properties are influenced by the arrangement of the collagen fibres, which are orthogonal to the bone-cartilage interface in the deep zone, randomly oriented in the middle zone, and parallel to the surface in the superficial zone. In past studies, cartilage permeability has been related directly to the orientation of the glycosaminoglycan chains attached to the proteoglycans which constitute the tissue matrix. These studies predicted permeability to be isotropic in the undeformed configuration, and anisotropic under compression. They neglected tissue anisotropy caused by the collagen network. However, magnetic resonance studies suggest that fluid flow is "directed" by collagen fibres in biological tissues. Therefore, the aim of this study was to express the permeability of cartilage accounting for the microstructural anisotropy and inhomogeneity caused by the collagen fibres. Permeability is predicted to be anisotropic and inhomogeneous, independent of the state of strain, which is consistent with the morphology of the tissue. Looking at the local anisotropy of permeability, we may infer that the arrangement of the collagen fibre network plays an important role in directing fluid flow to optimise tissue functioning.     

Effects of tension-compression nonlinearity on solute transport in charged hydrated fibrous tissues under dynamic unconfined compression

Cartilage is a charged hydrated fibrous tissue exhibiting a high degree of tension-compression nonlinearity (i.e., tissue anisotropy). The effect of tension-compression nonlinearity on solute transport has not been investigated in cartilaginous tissue under dynamic loading conditions. In this study, a new model was developed based on the mechano-electrochemical mixture model [Yao and Gu, 2007, J. Biomech. Model Mechanobiol., 6, pp. 63-72, Lai et al., 1991, J. Biomech. Eng., 113, pp. 245-258], and conewise linear elasticity model [Soltz and Ateshian, 2000, J. Biomech. Eng., 122, pp. 576-586; Curnier et al., 1995, J. Elasticity, 37, pp. 1-38]. The solute desorption in cartilage under unconfined dynamic compression was investigated numerically using this new model. Analyses and results demonstrated that a high degree of tissue tension-compression nonlinearity could enhance the transport of large solutes considerably in the cartilage sample under dynamic unconfined compression, whereas it had little effect on the transport of small solutes (at 5% dynamic strain level). The loading-induced convection is an important mechanism for enhancing the transport of large solutes in the cartilage sample with tension-compression nonlinearity. The dynamic compression also promoted diffusion of large solutes in both tissues with and without tension-compression nonlinearity. These findings provide a new insight into the mechanisms of solute transport in hydrated, fibrous soft tissues.     

Anisotropy of age-related toughness loss in human cortical bone: a finite element study

A mechanistic understanding of the role of bone quality on fracture processes is essential for determining the underlying causes of age-related changes in the mechanical response of the human bone. In this study, a previously developed cohesive finite element model was used to investigate the effects of age-related changes and the orientation of crack growth on the toughening behavior of human cortical bone. The change in the anisotropy of toughening mechanisms with age was also studied. Finite element method (FEM) simulations showed that the initiation toughness decreased by 3% and 8%/decade for transverse and longitudinal crack growth, respectively. In contrast, fracture resistance curve slope for transverse and longitudinal crack growth decreased by 2% and 3%/decade, respectively. Initiation fracture toughness values were higher for the transverse than for the longitudinal for a given age. On the other hand, propagation fracture toughness values were higher for longitudinal than for transverse crack growth for a given age. With respect to age, the toughness ratio for crack initiation decreased by 6%/decade, but that for propagation showed almost no change (less than 1%). In light of these findings, an analytical model evaluating the crack arresting feature of cement lines, is proposed to explain the factors that determine crack penetration into osteons or its deflection by cement lines.     

A finite deformation theory for cartilage and other soft hydrated connective tissues--I. Equilibrium results

The determination of valid stress-strain relations for articular cartilage under finite deformation conditions is a prerequisite for constructing models for synovial joint lubrication. Under physiological conditions of high strain rates and/or high stresses in the joint, large strains occur in cartilage. A finite deformation theory valid for describing cartilage, as well as other soft hydrated connective tissues under large loads, has been developed. This theory is based on the choice of a specific Helmholtz energy function which satisfies the generalized Coleman-Noll (GCN0) condition and the Baker-Ericksen (B-E) inequalities established in finite elasticity theory. In addition, the finite deformation biphasic theory includes the effects of strain-dependent porosity and permeability. These nonlinear effects are essential for properly describing the biomechanical behavior of articular cartilage, even when strain rates are low and strains are infinitesimal. The finite deformation theory describes the large strain behavior of cartilage observed in one-dimensional confined compression experiments at equilibrium, and it reduces to the linear biphasic theory under infinitesimal strain and slow strain rate conditions. Using this theory, we have determined the material coefficients of both human and bovine articular cartilages under large strain conditions at equilibrium. The theory compares very well with experimental results.     

Anisotropy and inter-condyle heterogeneity of cartilage under large-strain shear

We were the first to examine the mechanical responses of skeletally mature bovine femoral cartilage under large-strain simple shear (up to ±20%) using a multiaxial shear testing device. Since shear loading is critical in both tissue failure and chondrocyte responses, we aimed to probe (1) anisotropy with respect to the split-line direction (principal alignment of the collagen fibers near the articulating surface), (2) heterogeneity between femoral condyles, and (3) the influence of local cartilage thickness. We harvested a total of 48 cuboid cartilage specimens from four bovine knee joints. With each specimen we applied shear strains both parallel and perpendicular to the local split-line direction at a rate of 75 μm/min and calculated the peak-to-peak shear stresses, shear strain–energy dissipation densities, and peak effective shear moduli. The Wilcoxon signed rank test revealed that the medial condyle was anisotropic in some mechanical measures at applied shear strains above 5%, while the lateral condyle was mechanically isotropic at all applied shear strains. The Kruskal–Wallis test revealed no significant differences in the median mechanical behavior of the lateral and medial condyles. Spearman׳s rank correlations revealed statistically significant negative monotonic correlations among thickness and most of our mechanical measures for both lateral and medial condyles at most applied strains and directions of applied shear. These results suggest that large-strain analyses account for nonlinear, anisotropic and location-dependent effects not fully realized at small strains. Our findings may inspire new experiments and models that consider anisotropy and heterogeneity of cartilage in ways previously ignored.     

Solute transport in cartilage undergoing cyclic deformation

There are no blood vessels in cartilage to transport nutrients and growth factors to chondrocytes dispersed throughout the cartilage matrix. Insulin-like growth factor-I (IGF-I) is a large molecule with an important role in cartilage growth and metabolism, however, it first must reach the chondrocytes to exert its effect. While diffusion of IGF-I through cartilage is possible, it has been speculated that cyclic loading can enhance the rate of solute transport within cartilage. To better understand this process, here a one-dimensional axisymmetric mathematical model is developed to examine the transport of solutes through a cylindrical plug of cartilage undergoing cyclic axial deformation in the range of 10(-3) -1 Hz. This study has revealed the role of timescales in interpreting transport results in cartilage. It is shown that dynamic strains can either enhance or inhibit IGF-I transport at small timescales (< 20 min after onset of loading), depending on loading frequency. However, on longer timescales it is found that dynamic loading has negligible effect on IGF-I transport. Most importantly, in all cases examined the steady state IGF-I concentration did not exceed the fixed boundary value, in contrast to the predictions of Mauk et al. (2003).     

Anisotropy and inhomogeneity of the trabecular structure can describe the mechanical strength of osteoarthritic cancellous bone

Osteoarthritic cancellous bone was studied to investigate the development of this pathology, and the functional changes it induces in the bone. In order to predict how the morphological alterations of the tissue induced by the pathology can change the mechanical properties of the structure, two different strategies have been used in the literature: (1) emphasising the influence of structural anisotropy; (2) stressing the highly inhomogeneous characteristics of cancellous bone. The aim of the present study was to verify the theory that mechanical strength of osteoarthritic cancellous bone depends both on tissue anisotropy and inhomogeneity.Twenty-five specimens were extracted from osteoarthritic femoral heads, along selected directions, and analysed by means of a microtomograph. The same specimens were mechanically tested in compression to determine the mechanical strength. The most representative structural parameters, confirmed by a stepwise analysis, were used to define four models to describe the measured mechanical strength. The models were applied neglecting (global analysis) or considering (local analysis) tissue inhomogeneities to verify whether the correlation with ultimate stress could be improved.The coefficient of determination increased from 0.53, considering only bone volume fraction, up to 0.88, combining it with off-axis angle and normalised eigenvalue. A further improvement was found performing a local analysis (R²=0.90), which corresponded to a decrease of 17% in the residual error.The proposed approach of considering both tissue anisotropy and inhomogeneity improved the accuracy in predicting the mechanical behaviour of cancellous bone tissue and should be suitable for more general loading conditions.     

The deformation behavior and viscoelastic properties of chondrocytes in articular cartilage

Chondrocytes in articular cartilage utilize mechanical signals in conjunction with other environmental factors to regulate their metabolic activity. However, the sequence of biomechanical and biochemical events involved in the process of mechanical signal transduction has not been fully deciphered. A fundamental step in determining the role of various factors in regulating chondrocyte activity is to characterize accurately the biophysical environment within the tissue under physiological conditions of mechanical loading. Microscopic imaging studies have revealed that chondrocytes as well as their nuclei undergo shape and volume changes in a coordinated manner with deformation of the tissue matrix. Through micromechanical experiments, it has been shown that the chondrocyte behaves as a viscoelastic solid material with a mechanical stiffness that is several orders of magnitude lower than that of the cartilage extracellular matrix. These properties seem to be due to the structure of the chondrocyte cytoskeleton, and in part, the viscoelastic properties of the cell nucleus. The mechanical properties of the pericellular matrix that immediately surrounds the chondrocyte significantly differ from those of the chondrocyte and the extracellular matrix, suggesting that the pericellular matrix plays an important role in defining the mechanical environment of the chondrocyte. These experimentally measured values for chondrocyte and cartilage mechanical properties have been used in combination with theoretical constitutive modeling of the chondrocyte within articular cartilage to predict the non-uniform and time-varying stress-strain and fluid flow environment of the cell. The ultimate goal of these studies has been to elucidate the sequence of biomechanical and biochemical events through which mechanical stress influences chondrocyte activity in both health and in disease.     

Methodology and sensitivity studies for finite element modeling of the inferior glenohumeral ligament complex

The objectives of this research were to develop a methodology for three-dimensional finite element (FE) modeling of the inferior glenohumeral ligament complex (IGHL complex) as a continuous structure, to determine optimal mesh density for FE simulations, to examine strains and forces in the IGHL complex in clinically relevant joint positions, and to perform sensitivity studies to assess the effects of assumed material properties. A simple translation test in the anterior direction was performed on a cadaveric shoulder, with the humerus oriented at 60 degrees of glenohumeral abduction and 0 degrees of flexion/extension, at 0 degrees , 30 degrees and 60 degrees of humeral external rotation. The geometries of the relevant structures were extracted from volumetric CT data to create a FE model. Experimentally measured kinematics were applied to the FE model to simulate the simple translation test. First principal strains, insertion site forces and contact forces were analyzed. At maximum anterior humeral translation, strains in the IGHL complex were highly inhomogeneous for all external rotation angles. The motion of the humerus with respect to the glenoid during the simple translation test produced a tangential load at the proximal and distal edges of the IGHL complex. This loading was primarily in the plane of the inferior glenohumeral ligament complex, producing an in-plane shear-loading pattern. There was a significant increase in strain with increasing angle of external rotation. The largest insertion site forces occurred at the axillary pouch insertion to the humerus (36.7N at 60 degrees of external rotation) and the highest contact forces were between the anterior band of the IGHL complex and the humeral cartilage (7.3N at 60 degrees of external rotation). Strain predictions were highly sensitive to changes in the ratio of bulk to shear modulus of the IGHL complex, while predictions were moderately sensitive to changes in elastic modulus of the IGHL complex. Changes to the material properties of the humeral cartilage had little effect on predicted strains. The methodologies developed in this research and the results of the mesh convergence and sensitivity studies provide a basis for the subject-specific modeling of the mechanics of the IGHL complex.     

Uncalcified cartilage resorption in human fetal cartilage canals

In the human fetus, epiphyses appear as a solid avascular cartilaginous mass until the eleventh week of development. Around the third fetal month of development, vascular canals coming from the perichondrium are recognized in the mineralized epiphyseal cartilage. Whether cartilage canals develop by passive inclusion or active chondrolysis is still a matter of controversy. We studied the relationships between the intracanalar cells and the surrounding matrix on human fetal epiphyses embedded in glycol methacrylate. At the blind end of canals both stellate fibroblast-like cells and vacuolated macrophages are observed. These cellular foci show all characteristics of active chondrolysis (loss of metachromasia, lacunae containing cells intimately associated with matrix, and presence of granular debris). Similar resorptive foci have been observed in the pannus of rheumatoid joints and in the embryonic chick growth plate composed of uncalcified cartilage. A cellular cooperation (fibroblast/macrophage) is necessary for uncalcified cartilage breakdown. In the human fetus, monocytes/macrophages have been recognized in the peripheral blood as early as the twelfth week of gestation. Our observations support the view that chondrolysis due to both fibroblasts (of mesenchymal origin) and macrophages is the basic mechanism for cartilage canal development.     

The micromechanical environment of intervertebral disc cells determined by a finite deformation,anisotropic, and biphasic finite element model

Cellular response to mechanical loading varies between the anatomic zones of the intervertebral disc. This difference may be related to differences in the structure and mechanics of both cells and extracellular matrix, which are expected to cause differences in the physical stimuli (such as pressure, stress, and strain) in the cellular micromechanical environment. In this study, a finite element model was developed that was capable of describing the cell micromechanical environment in the intervertebral disc. The model was capable of describing a number of important mechanical phenomena: flow-dependent viscoelasticity using the biphasic theory for soft tissues; finite deformation effects using a hyperelastic constitutive law for the solid phase; and material anisotropy by including a fiber-reinforced continuum law in the hyperelastic strain energy function. To construct accurate finite element meshes, the in situ geometry of IVD cells were measured experimentally using laser scanning confocal microscopy and three-dimensional reconstruction techniques. The model predicted that the cellular micromechanical environment varies dramatically between the anatomic zones, with larger cellular strains predicted in the anisotropic anulus fibrosus and transition zone compared to the isotropic nucleus pulposus. These results suggest that deformation related stimuli may dominate for anulus fibrosus and transition zone cells, while hydrostatic pressurization may dominate in the nucleus pulposus. Furthermore, the model predicted that micromechanical environment is strongly influenced by cell geometry, suggesting that the geometry of IVD cells in situ may be an adaptation to reduce cellular strains during tissue loading.     

Expression of NGF, Trka and p75 in human cartilage

Nerve growth factor (NGF) exerts its action through two types of receptor: high-affinity tyrosine kinase A receptor (trkA) and low-affinity p75 receptor. NGF has a neurotrophic role in central and peripheral nervous system development, but there is also clear evidence of its involvement in the developing skeleton. The aim of the present immunohistochemical study was to investigate the expression and distribution of NGF, trkA, and p75 in normal cartilaginous tissues from adult subjects: articular and meniscal cartilage of the knee, cartilage from the epiglottis, and intervertebral disc tissue. Detection of NGF mRNA was also performed by in situ hybridization. Immunoreaction for NGF and the two receptors in articular chondrocytes, chondrocyte-like cells of meniscus and annulus fibrosus, and chondrocytes of the epiglottis demonstrated that they are all expressed in hyaline, fibrous and elastic cartilaginous tissues, suggesting that they could be involved in cartilage physio-pathology.     

Shear deformation kinematics during cartilage articulation: effect of lubrication, degeneration, and stress relaxation

During joint articulation, the biomechanical behavior of cartilage not only facilitates load-bearing and low-friction, but also provides regulatory cues to chondrocytes. Elucidation of cartilage kinematics under combined compression and shearing conditions clarifies these cues in health and disease. The objectives of this study were to elucidate the effects of lubricant, tissue degeneration, and stress relaxation duration on cartilage shear kinematics during articulation. Human osteochondral cores with normal and mildly degenerate surface structures were isolated. Paired blocks from each core were apposed, compressed, allowed to stress relax for 5 or 60 min, and shear tested with a micro-scale video microscopy system using phosphate-buffered saline (PBS) or synovial fluid as lubricant. During applied lateral motion, local and overall shear strain (Exz) of articular cartilage were determined. The applied lateral displacement at which Exz reached 50% of the peak (Deltax(1/2)) was also determined. Quantitatively, surface Exz increased at the onset of lateral motion and peaked just as surfaces detached and slid. With continued lateral motion, surface Exz was maintained. After short stress relaxation, effects of lubrication on Exz and Deltax(1/2) were not apparent. With prolonged stress relaxation, Exz and Deltax(1/2) near the articular surface increased markedly when PBS was used as lubricant. Similar patterns were observed for overall Exz and Deltax(1/2). With degeneration, surface Exz was consistently higher for all cases after the onset of lateral motion. Thus, cartilage shear kinematics is markedly affected by lubricant, cartilage degeneration, and loading duration. Changes in these factors may be involved in the pathogenesis of osteoarthritis.     

A nonlinear finite element model of cartilage growth

The long range objective of this work is to develop a cartilage growth finite element model (CGFEM), based on the theories of growing mixtures that has the capability to depict the evolution of the anisotropic and inhomogeneous mechanical properties, residual stresses, and nonhomogeneities that are attained by native adult cartilage. The CGFEM developed here simulates isotropic in vitro growth of cartilage with and without mechanical stimulation. To accomplish this analysis a commercial finite element code (ABAQUS) is combined with an external program (MATLAB) to solve an incremental equilibrium boundary value problem representing one increment of growth. This procedure is repeated for as many increments as needed to simulate the desired growth protocol. A case study is presented utilizing a growth law dependent on the magnitude of the diffusive fluid velocity to simulate an in vitro dynamic confined compression loading protocol run for 2 weeks. The results include changes in tissue size and shape, nonhomogeneities that develop in the tissue, as well as the variation that occurs in the tissue constitutive behavior from growth.     

Prestressing in finite deformation abdominal aortic aneurysm simulation

In abdominal aortic aneurysm (AAA) simulation the patient-specific geometry of the object of interest is very often reconstructed from in vivo medical imaging such as CT scans. Such geometries represent a deformed configuration stressed by typical in vivo conditions. However, commonly, such structures are considered stress-free in simulation. In this contribution we sketch and compare two methods to introduce a physically meaningful stress/strain state to the obtained geometry for simulations in the finite strain regime and demonstrate the necessity of such prestressing techniques. One method is based on an inverse design analysis to calculate a stress-free reference configuration. The other method developed here is based on a modified updated Lagrangian formulation. Formulation of both methods is provided. Applicability and accurateness of both approaches are compared and evaluated utilizing fully three-dimensional patient-specific AAA structures in the finite strain regime.