Fabric dependence of wave propagation in anisotropic porous media

Current diagnosis of bone loss and osteoporosis is based on the measurement of the bone mineral density (BMD) or the apparent mass density. Unfortunately, in most clinical ultrasound densitometers: 1) measurements are often performed in a single anatomical direction, 2) only the first wave arriving to the ultrasound probe is characterized, and 3) the analysis of bone status is based on empirical relationships between measurable quantities such as speed of sound (SOS) and broadband ultrasound attenuation (BUA) and the density of the porous medium. However, the existence of a second wave in cancellous bone has been reported, which is an unequivocal signature of poroelastic media, as predicted by Biot’s poroelastic wave propagation theory. In this paper, the governing equations for wave motion in the linear theory of anisotropic poroelastic materials are developed and extended to include the dependence of the constitutive relations upon fabric—a quantitative stereological measure of the degree of structural anisotropy in the pore architecture of a porous medium. This fabric-dependent anisotropic poroelastic approach is a theoretical framework to describe the microarchitectural-dependent relationship between measurable wave properties and the elastic constants of trabecular bone, and thus represents an alternative for bone quality assessment beyond BMD alone.     

Interstitial fluid flow in the osteon with spatial gradients of mechanical properties: a finite element study

Bone remodelling is the process that maintains bone structure and strength through adaptation of bone tissue mechanical properties to applied loads. Bone can be modelled as a porous deformable material whose pores are filled with cells, organic material and interstitial fluid. Fluid flow is believed to play a role in the mechanotransduction of signals for bone remodelling. In this work, an osteon, the elementary unit of cortical bone, is idealized as a hollow cylinder made of a deformable porous matrix saturated with an interstitial fluid. We use Biot’s poroelasticity theory to model the mechanical behaviour of bone tissue taking into account transverse isotropic mechanical properties. A finite element poroelastic model is developed in the COMSOL Multiphysics software. Elasticity equations and Darcy’s law are implemented in this software; they are coupled through the introduction of an interaction term to obtain poroelasticity equations. Using numerical simulations, the investigation of the effect of spatial gradients of permeability or Poisson’s ratio is performed. Results are discussed for their implication on fluid flow in osteons: (i) a permeability gradient affects more the fluid pressure than the velocity profile; (ii) focusing on the fluid flow, the key element of loading is the strain rate; (iii) a Poisson’s ratio gradient affects both fluid pressure and fluid velocity. The influence of textural and mechanical properties of bone on mechanotransduction signals for bone remodelling is also discussed.     

The estimated elastic constants for a single bone osteonal lamella

Micromechanical estimates of the elastic constants for a single bone osteonal lamella and its substructures are reported. These estimates of elastic constants are accomplished at three distinct and organized hierarchical levels, that of a mineralized collagen fibril, a collagen fiber, and a single lamella. The smallest collagen structure is the collagen fibril whose diameter is the order of 20 nm. The next structural level is the collagen fiber with a diameter of the order of 80 nm. A lamella is a laminate structure, composed of multiple collagen fibers with embedded minerals and consists of several laminates. The thickness of one laminate in the lamella is approximately 130 nm. All collagen fibers in a laminate in the lamella are oriented in one direction. However, the laminates rotate relative to the adjacent laminates. In this work, all collagen fibers in a lamella are assumed to be aligned in the longitudinal direction. This kind of bone with all collagen fibers aligned in one direction is called a parallel fibered bone. The effective elastic constants for a parallel fibered bone are estimated by assuming periodic substructures. These results provide a database for estimating the anisotropic poroelastic constants of an osteon and also provide a database for building mathematical or computational models in bone micromechanics, such as bone damage mechanics and bone poroelasticity.     

Impact of the porous microstructure on the overall elastic properties of the osteonal cortical bone

Mechanical properties of bones are largely determined by their microstructure. The latter comprises a large number of diverse pores. The present paper analyzes a connection between structure of the porous space of the osteonal cortical bone and bone's overall anisotropic elastic moduli. The analysis is based on recent developments in the theory of porous materials that predict the anisotropic effective moduli of porous solids in terms of pores' shapes, orientations and densities. Bone's microstructure is modeled using available micrographs. The calculated anisotropic elastic constants for porous cortical bone are, mostly, in agreement with available experimental data. The influence of each of the pore types on the overall moduli is examined. The results of the analysis can also be used to estimate the extent of mineralization (hydroxyapatite content) if the overall porosity and the effective moduli are known and, vice versa, to estimate porosity from the measured moduli and the extent of mineralization.     

On the infusion of a therapeutic agent into a solid tumor modeled as a poroelastic medium

The direct infusion of an agent into a solid tumor, modeled as a spherical poroelastic material with anisotropic dependence of the tumor hydraulic conductivity upon the tissue deformation, is treated both by solving the coupled fluid/elastic equations, and by expressing the solution as an asymptotic expansion in terms of a small parameter, ratio between the driving pressure force in the fluid system, and the elastic properties of the solid. Results at order one match almost perfectly the solutions of the full system over a large range of infusion pressures. Comparison with experimental results is acceptable after the hydraulic conductivity of the medium is properly calibrated. Given the uncertain estimates of some model constants, the order zero solution of the expansion, for which fluid and porous matrix are decoupled, yields acceptable values and trends for all the physical fields of interest, rendering the coupled analysis (in the limit of small displacements) of little use. When the deformation of the tissue becomes large nonlinear elasticity theory must be resorted to.     

Structural models for human spinal motion segments based on a poroelastic view of the intervertebral disk

Analytical and finite element models (FEMs) were used to quantify poroelastic material properties for a human intervertebral disk. An axisymmetric FEM based on a poroelastic view of disk constituents was developed for a representative human spinal motion segment (SMS). Creep and steady-state response predicted by FEMs agreed with experimental observations, i.e., long-time creep occurs with flow in the SMS, whereas for rapid steady-state loading an "undrained," nearly incompressible response is evident. A relatively low value was determined for discal permeability. Transient and long-term creep FE analyses included the study of deformation, pore fluid flow, stress, and pore fluid pressure. Relative fluid motion associated with transient creep is related to nuclear nutrition and the overall mechanical response in the normal disk. Degeneration of the disk may be associated with an increase in permeability.     

Collagen Network of Articular Cartilage Modulates Fluid Flow and Mechanical Stresses in Chondrocyte

The extracellular matrix of articular cartilage modulates the mechanical signals sensed by the chondrocytes. In the present study, a finite element model (FEM) of the chondrocyte and its microenvironment was reconstructed using the information from fourier transform infrared imaging spectroscopy. This environment consisted of pericellular, territorial (mainly proteoglycans), and inter-territorial (mainly collagen) matrices. The chondrocyte, pericellular, and territorial matrix were assumedto be mechanically isotropic and poroelastic, whereas the inter-territorial matrix, due to its high collagen content, was assumed to be transversely isotropic and poroelastic. Under instantaneous strain-controlled compression, the FEM indicated that the fluid pressure within the chondrocyte increased nonlinearly as a function of the in-plane Young’s modulus of the collagen network. Under instantaneous force-controlled compression, the chondrocyte experienced the highest fluid pressure when the in-plane Young’s modulus of the collagen network was ~4 MPa. Based on the present results, the mechanical characteristics of the collagen network of articular cartilage can modify fluid flow and stresses in chondrocytes. Therefore, the integrity of the collagen network may be an important determinant in cell stimulation and in the control of the matrix maintenance.     

Estimates of the peak pressures in bone pore water.

The maximum pore fluid pressures due to uniaxial compression are determined for both the vascular porosity (Haversian and Volkmann's canals) and the lacunar-canalicular porosity of live cortical bone. It is estimated that the peak pore water pressure will be 19 percent of the applied axial stress in the vascular porosity and 12 percent of the applied axial stress in the lacunar-canalicular porosity for an impulsive step loading. However, the estimated relaxation time for the vascular porosity (1.36 microseconds) is three orders of magnitude faster than that estimated for the lacunar-canalicular porosity (4.9 ms). Thus, under physiological loading, which has a stress rise time generally larger than 1 ms, pressures higher than the vascular pressure cannot be sustained in the vascular porosity due to the swift pressure relaxation in this porosity (unless the fluid drainage through the boundary is obstructed). The model also predicts a slight hydraulic stiffening of the bulk modulus due to longer draining time of the lacunar-canalicular porosity. The undrained bulk modulus is 6 percent higher than the drained bulk modulus in this case.     

Estimation of the effective transversely isotropic elastic constants of a material from known values of the material's orthotropic elastic constants

A method is illustrated for determining the effective transversely isotropic (or isotropic) elastic constants from measured orthotropic elastic constants. This method consists of constructing upper and lower bounds on the effective transversely isotropic (or isotropic) elastic constants using the known orthotropic values. This method is illustrated using three sets of elastic constants for bone. Fortunately, the upper and lower bounds are very close. Thus very good approximations for the effective transversely isotropic (or isotropic) elastic constants for cortical and cancellous bone are obtained from previously published data on the orthotropic elastic constants for those tissue types. This work is undertaken to build a greater database for the transversely isotropic elastic constants of bone with the intention of employing them in a transversely isotropic model of bone poroelasticity. An interesting aspect of the present result is that the Voigt and Reuss bounds are very tight for these anisotropic materials. This is not always the case for these bounds. Received: 14 November 2001 / Accepted: 25 February 2002     

Effects of the torso boundary and internal conductivity interfaces in electrocardiography: An evaluation of the ‘Infinite medium’ approximation

The analytic, eccentric spheres model of the torso was used to examine the validity of approximating the ‘infinite medium’ potential by integrating ‘finite medium potentials’ over the torso surface. Although idealized, the analytic model is sophisticated enough for all important torso conductivity and geometry parameters to be preserved in the formulation. The model generates both ‘finite medium’ potentials (for which the torso is surrounded by air) and also ‘infinite medium’ potentials (for which the outermost layer of the torso extends outward to infinity). The finite medium torso potentials were integrated over the torso surface in accordance with the approximation used by many investigators in an effort to make the surface distribution more representative of the primary cardiac sources. The resulting potential distribution was compared with the true infinite medium potential, in which the effects of internal inhomogeneities (secondary sources) were taken into account. The difference between the two representations was found to be significant, and caution should be used when interpreting such data.     

Elastic anisotropy of bone lamellae as a function of fibril orientation pattern

In this study, the homogenized anisotropic elastic properties of single bone lamellae are computed using a finite element unit cell method. The resulting stiffness tensor is utilized to calculate indentation moduli for multiple indentation directions in the lamella plane which are then related to nanoindentation experiments. The model accounts for different fibril orientation patterns in the lamellae—the twisted and orthogonal plywood pattern, a 5-sublayer pattern and an X-ray diffraction-based pattern. Three-dimensional sectional views of each pattern facilitate the comparison to transmission electron (TEM) images of real lamella cuts. The model results indicate, that the 5-sublayer- and the X-ray diffraction-based patterns cause the lamellae to have a stiffness maximum between 0° and 45° to the osteon axis. Their in-plane stiffness characteristics are qualitatively matching the experimental findings that report a higher stiffness in the osteon axis than in the circumferential direction. In contrast, lamellae owning the orthogonal or twisted plywood fibril orientation patterns have no preferred stiffness alignment. This work shows that the variety of fibril orientation patterns leads to qualitative and quantitative differences in the lamella elastic mechanical behavior. The study is a step toward a deeper understanding of the structure—mechanical function relationship of bone lamellae.     

A dynamical study of the mechanical stimuli and tissue differentiation within a CaP scaffold based on micro-CT finite element models

The control of the mechanical stimuli transmitted to the cells is critical for the design of functional scaffolds for tissue engineering. The objective of this study was to investigate the dynamics of the mechanical stimuli transmitted to the cells during tissue differentiation in an irregular morphology scaffold under compressive load and perfusion flow. A calcium phosphate-based glass porous scaffold was used. The solid phase and the fluid flow within the pores were modeled as linear elastic solid material and Newtonian fluid, respectively. In the fluid model, different levels of viscosity were used to simulate tissue differentiation. Compressive strain of 0.5% and fluid flow with constant inlet velocity of 10 μm/s or constant inlet pressure of 3 Pa were applied. Octahedral shear strain and fluid shear stress were used as mechano-regulatory stimuli. For constant inlet velocity, stimuli equivalent to bone were predicted in 80% of pore volume for the case of low tissue viscosity. For the cases of high viscosity, fluctuations between stimuli equivalent to tissue formation and cell death were predicted due to the increase in the fluid shear stress when tissue started to fill pores. When constant pressure was applied, stimuli equivalent to bone were predicted in 62% of pore volume when low tissue viscosity was used and 42% when high tissue viscosity was used. This study predicted critical variations of fluid shear stress when cells differentiated. If these variations are not controlled in vitro, they can impede the formation of new matured tissue.     

Estimating the dielectric constant of the channel protein and pore

When modelling biological ion channels using Brownian dynamics (BD) or Poisson–Nernst–Planck theory, the force encountered by permeant ions is calculated by solving Poisson’s equation. Two free parameters needed to solve this equation are the dielectric constant of water in the pore and the dielectric constant of the protein forming the channel. Although these values can in theory be deduced by various methods, they do not give a reliable answer when applied to channel-like geometries that contain charged particles. To determine the appropriate values of the dielectric constants, here we solve the inverse problem. Given the structure of the MthK channel, we attempt to determine the values of the protein and pore dielectric constants that minimize the discrepancies between the experimentally-determined current–voltage curve and the curve obtained from BD simulations. Two different methods have been applied to determine these values. First, we use all possible pairs of the pore dielectric constant of water, ranging from 20 to 80 in steps of 10, and the protein dielectric constant of 2–10 in steps of 2, and compare the simulated results with the experimental values. We find that the best agreement is obtained with experiment when a protein dielectric constant of 2 and a pore water dielectric constant of 60 is used. Second, we employ a learning-based stochastic optimization algorithm to pick out the optimum combination of the two dielectric constants. From the algorithm we obtain an optimum value of 2 for the protein dielectric constant and 64 for the pore dielectric constant.     

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.     

Linear poroelastic cancellous bone anisotropy: trabecular solid elastic and fluid transport properties

The mechanical performance of cancellous bone is characterized using experiments which apply linear poroelasticity theory. It is hypothesized that the anisotropic organization of the solid and pore volumes of cancellous bone can be physically characterized separately (no deformable boundary interactive effects) within the same bone sample. Due to its spongy construction, the in vivo mechanical function of cancellous or trabecular bone is dependent upon fluid and solid materials which may interact in a hydraulic, convective fashion during functional loading. This project provides insight into the organization of the tissue, ie., the trabecular connectivity, by defining the separate nature of this biphasic performance. Previous fluid flow experiments [Kohles et al., 2001, Journal of Biomechanics, 34(11), pp. 1197-1202] describe the pore space via orthotropic permeability. Ultrasonic wave propagation through the trabecular network is used to describe the solid component via orthotropic elastic moduli and material stiffness coefficients. The linear poroelastic nature of the tissue is further described by relating transport (fluid flow) and elasticity (trabecular load transmission) during regression analysis. In addition, an empirical relationship between permeability and porosity is applied to the collected data. Mean parameters in the superior-inferior (SI) orientation of cubic samples (n=20) harvested from a single bovine distal femur were the largest (p<0.05) in comparison to medial-lateral (ML) and anterior-posterior (AP) orientations: Apparent elastic modulus (2,139 MPa), permeability (4.65x10(-10) m2), and material stiffness coefficient (13.6 GPa). A negative correlation between permeability as a predictor of structural elastic modulus supported a parametric relationship in the ML (R2=0.4793), AP (R2=0.3018), and SI (R2=0.6445) directions (p<0.05).     

Influence of boundary conditions on computed apparent elastic properties of cancellous bone

High-resolution finite element models of trabecular bone can be used to study trabecular structure–function relationships, elasticity, multiaxial strength, and tissue remodelling in more detail than experiments. Beside effects of the model size, scan/analysis resolution, segmentation process, etc., the type of the applied boundary conditions (BCs) have a strong influence on the predicted elastic properties. Appropriate BCs have to be applied on hexahedral digital finite element models in order to obtain effective elastic properties. Homogeneous displacement BCs as proposed by Van Rietbergen et al. (J Biomech 29(12):1653–1657, 1996) lead to “apparent” rather than to “effective” elastic properties. This study provides some answers concerning such differences by comparing various BC types (uniform displacement, mixed BCs, periodic BCs), different volume element definitions (original and mirrored models), and several bone volume fractions (BVTV ranging from 6.5 to 37.6%). First, the mixed BCs formulated by Hazanov (Arch Appl Mech 68(6):385–394, 1998) are theoretically extended to shear loading of a porous media. Second, six human bone samples are analyzed, their orthotropic Young’s moduli, shear moduli, and Poisson’s ratios computed and compared. It is found that the proposed mixed BCs give exactly the same effective elastic properties as periodic BCs if a periodic and orthotropic micro-structured material is used and thus denoted as “periodicity compatible” mixed uniform BCs (PMUBCs). As bone samples were shown to be nearly orthotropic for volume element side lengths ≥5 mm the proposed mixed BCs turn out to be the best choice because they give again essentially the same overall elastic properties as periodic BCs. For bone samples of smaller dimensions ( < 5 mm) with a strong anisotropy (beyond orthotropy) uniform displacement BCs remain applicable but they can significantly overestimate the effective stiffness. In Memoriam, Prof. Christian Huet.     

Assessment of composition and anisotropic elastic properties of secondary osteon lamellae

Measurement of the elastic properties of single osteon lamellae is still one of the most demanding tasks in bone mechanics to be solved. By means of site-matched Raman microspectroscopy, acoustic microscopy and nanoindentation the structure, chemical composition and anisotropic elasticity of individual lamellae in secondary osteons were investigated. Acoustic impedance images (911-MHz) and two-dimensional Raman spectra were acquired in sections of human femoral bone. The samples were prepared with orientations at various observation angles theta relative to the femoral long axis. Nanoindentations provided local estimations of the elastic modulus and landmarks necessary for spatial fusion of the acoustic and spectral Raman images. Phosphate nu(1) (961 cm(-1)) and amide I (1665 cm(-1)) band images representing spatial distributions of mineral and collagen were fused with the acoustic images. Acoustic impedance was correlated with the indentation elastic modulus E(IT) (R(2)=0.61). Both parameters are sensitive to elastic tissue anisotropy. The lowest values were obtained in the direction perpendicular to the femoral long axis. Acoustic images exhibit a characteristic bimodal lamellar pattern of alternating high and low impedance values. Since this undulation was not associated with a variation of the phosphate nu(1)-band intensity in the Raman images, it was attributed to variations of the lamellar orientation. After threshold segmentation and conversion to elastic modulus the orientation and transverse isotropic elastic constants were derived for individual ensembles of apparent thin and thick lamellae. Our results suggest that this model represents the effective anisotropic properties of an asymmetric twisted plywood structure made of transverse isotropic fibrils. This is the first report that proves experimentally the ability of acoustic microscopy to map tissue elasticity in two dimensions with micrometer resolution. The combination with Raman microspectroscopy provides a unique way to study bone and mineral metabolism and the relation with mechanical function at the ultrastructural tissue level.     

FGF2 and dexamethasone increase the production of hyaluronan in two-dimensional culture of elastic cartilage-derived cells: in vitro analyses and in vivo cartilage formation

Elastic cartilage-derived cells cultured two-dimensionally with FGF2 and corticosteroid produce gel-type masses that become mature cartilage when injected into a subcutaneous pocket. This unique method has previously been clinically applied for treatments of nasal augmentation. However, the components of the gel-type mass and the mechanism of its synthesis remain unknown. Here, we have investigated the components of the gel-type mass produced by elastic cartilage-derived cells, and whether this gel-type mass can be produced by using other cell sources or other media. Human elastic cartilage-derived cells from auricular cartilage, hyaline cartilage-derived cells from articular cartilage, and mesenchymal stem cells from synovium were cultured in three media: “redifferentiation medium” containing FGF2 and dexamethasone; “chondrogenic medium” containing bone morphogenetic protein-2, transforming growth factor-β3, and dexamethasone specific for in vitro chondrogenesis of mesenchymal stem cells; control medium. The elastic cartilage-derived cells cultured in redifferentiation medium produced a gelatinous matrix positive for Alcian blue. During culture, the amount of chondroitin 4-sulfate, chondroitin 6-sulfate, and especially hyaluronan increased. However, the expression of RNAs for most chondrogenic genes did not increase. We also reproduced cartilage tissue formation by the injection of elastic cartilage-derived cells with the gelatinous mass into the subcutaneous space of the nude mouse. The synthesis of gelatinous matrix in vitro and the formation of cartilage tissue in vivo could be obtained only for the combination of elastic cartilage-derived cells with redifferentiation medium. This study was supported in part by grants from the “Japan Society for the Promotion of Science (19591752)” and “Center of Excellence Program for Frontier Research on Molecular Destruction and Reconstruction of Tooth and Bone in Tokyo Medical and Dental University” to Takeshi Muneta, and the “Japan Society for the Promotion of Science (18591657)” to Ichiro Sekiya.