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     

Errors induced by off-axis measurement of the elastic properties of bone

Misalignment between the axes of measurement and the material symmetry axes of bone causes error in anisotropic elastic property measurements. Measurements of Poisson's ratio were strongly affected by misalignment errors. The mean errors in the measured Young's moduli were 9.5 and 1.3 percent for cancellous and cortical bone, respectively, at a misalignment angle of 10 degrees. Mean errors of 1.1 and 5.0 percent in the measured shear moduli for cancellous and cortical bone, respectively, were found at a misalignment angle of 10 degrees. Although, cancellous bone tissue was assumed to have orthotropic elastic symmetry, the possibility of the greater symmetry of transverse isotropy was investigated. When the nine orthotropic elastic constants were forced to approximate the five transverse isotropic elastic constants, errors of over 60 percent were introduced. Therefore, it was concluded that cancellous bone is truly orthotropic and not transversely isotropic. A similar but less strong result for cortical bone tissue was obtained.     

Direct mechanics assessment of elastic symmetries and properties of trabecular bone architecture

A method is presented to find orthotropic elastic symmetries and constants directly from the elastic coefficients in the overall stiffness matrix of trabecular bone test specimens. Contrary to earlier developed techniques, this method does not require pure orthotropic behavior or additional fabric measurements. The method uses high-resolution computer reconstructions of trabecular bone specimens as input for large-scale FE-analyses to determine all the 21 elastic coefficients in the overall stiffness matrix of the specimen, using a direct mechanics approach. An optimization procedure is then used to find the coordinate transformation that yields the best orthotropic representation of this matrix. The method is illustrated here relative to two trabecular bone specimens. The techniques developed here can be used to obtain a complete characterization of the mechanical properties of trabecular architecture. With the development of in vivo reconstruction techniques, even in vivo measurements will be possible.     

The bone tissue of the canine mandible is elastically isotropic

This paper reports experimental measurements which show that canine mandibular bone tissue is elastically isotropic. Earlier work has established that human, canine and bovine cortical bone tissue of the femur, tibia and skull are elastically anisotropic and therefore the reported isotropy of mandibular tissue was unexpected. The isotropic elastic moduli of the canine mandible are represented by a Young's modulus of 7.5 GPa and a Poisson's ratio of 0.4. Earlier work gave the three orthotropic Young's moduli of the cortical one of the canine femur as 12.8 GPa, 15.6 GPa and 20.1 GPa. The experimental technique employed is elastic wave propagation at ultrasonic frequencies.     

Deviatoric and hydrostatic mode interaction in hard and soft tissue

It has been established that many hard and soft tissues have anisotropic material symmetry. It is noted here that the deviatoric and hydrostatic modes interact with each other in a general anisotropic elastic material. In the special case of isotropic, linear elastic, materials these modes are non-interactive. As a consequence of the interaction of these modes encountered in anisotropic materials, the decomposition into hydrostatic and deviatoric modes, and deviatoric mode concepts such as the von Mises effective stress are not appropriate for anisotropic materials in general. The implications of this observation for the presentation of computationally generated stress contours for hard and soft tissues are discussed. It is also pointed out that the mode coupling and mode interaction raise the question of whether anisotropic living tissues respond directly to stress or to some other physical quantity such as strain or strain energy, in view of the recent hypothesis concerning the proliferation and ossification of cartilage.     

Couple-stress moduli of a trabecular bone idealized as a 3D periodic cellular network

Trabecular bone is modeled as a cellular material with an idealized periodic structure made of open cubic cells, which is effectively orthotropic. We evaluate apparent couple-stress moduli of such a periodic material; apparent moduli refer to the moduli obtained using a domain smaller than a Representative Volume Element and they depend on boundary conditions. We conduct this analysis computationally (using ANSYS) by subjecting a unit cell of this periodic cellular material to either displacement or traction boundary conditions. Cell walls, representing bone tissue, and void space, representing bone marrow, are both modeled and they are assumed to be linear elastic. The applied loadings include a uniaxial extension (or uniaxial stress), a hydrostatic deformation (or hydrostatic stress) and a shear deformation (or shear stress) to evaluate the first stiffness (or compliance) tensor, and an applied curvature (or bending moment), a uniaxial twist (or torsion), and a triaxial twist (or triaxial torsion) to evaluate the second couple-stress stiffness (or compliance) tensor. Apparent couple-stress moduli are computed by equating the total strain energy stored in the unit cell with the energy of an equivalent homogeneous orthotropic couple-stress material for each applied loading. The moduli computed using displacement boundary conditions give upper bound, while those obtained using traction boundary conditions give lower bound on effective couple-stress moduli. These bounds are very wide due to a large mismatch in elastic moduli of bone tissue and bone marrow. These results are in agreement with our studies on composite materials with very stiff or very compliant inclusions.     

Errors in the orientation of the principal stress axes if bone tissue is modeled as isotropic

The error in the prediction of the orientation of the principal axes of stress in bone tissue is determined in the case when the tissue is modeled as elastically isotropic rather than as orthotropic, the probable symmetry of bone tissue. Results are two-dimensional and assume the same underlying strain state for both the orthotropic and isotropic cases. The maximum error is 45 degrees, and the typical error is generally significant.     

Tensile behavior of cortical bone: dependence of organic matrix material properties on bone mineral content

A porous composite model is developed to analyze the tensile mechanical properties of cortical bone. The effects of microporosity (volksman's canals, osteocyte lacunae) on the mechanical properties of bone tissue are taken into account. A simple shear lag theory, wherein tensile loads are transferred between overlapped mineral platelets by shearing of the organic matrix, is used to model the reinforcement provided by mineral platelets. It is assumed that the organic matrix is elastic in tension and elastic-perfectly plastic in shear until it fails. When organic matrix shear stresses at the ends of mineral platelets reach their yield values, the stress-strain curve of bone tissue starts to deviate from linear behavior. This is referred as the microscopic yield point. At the point where the stress-strain behavior of bone shows a sharp curvature, the organic phase reaches its shear yield stress value over the entire platelet. This is referred as the macroscopic yield point. It is assumed that after macroscopic yield, mineral platelets cannot contribute to the load bearing capacity of bone and that the mechanical behavior of cortical bone tissue is determined by the organic phase only. Bone fails when the principal stress of the organic matrix is reached. By assuming that mechanical properties of the organic matrix are dependent on bone mineral content below the macroscopic yield point, the model is used to predict the entire tensile mechanical behavior of cortical bone for different mineral contents. It is found that decreased shear yield stresses and organic matrix elastic moduli are required to explain the mechanical behavior of bones with lowered mineral contents. Under these conditions, the predicted values (elastic modulus, 0.002 yield stress and strain, and ultimate stress and strain) are within 15% of experimental data.     

Mechanical strength of bone allografts subjected to chemical sterilization and other terminal processing methods

Infectious disease transmission through the use of human donor allografts can be a catastrophic complication in an otherwise straightforward surgical procedure. The use of bone allograft in reconstructive orthopedic surgeries is increasing, yet severe complications, including death, can result if the transplanted tissues transmit a communicable disease to the tissue recipient. The BioCleanse((R)) tissue sterilization process is a fully automated, low-temperature chemical sterilization process that renders allograft tissue sterile. The purpose of this study was to evaluate the effect of a chemical tissue sterilization process on the mechanical strength of cortical bone allografts prior to implantation. Cylindrical cortical bone specimens were harvested from seven human cadaver donors and treated either by: chemical sterilization alone; chemical sterilization and terminal sterilization by gamma irradiation; chemical sterilization, lyophilization, terminal sterilization by STERRAD and rehydration; or untreated. The specimens were tested to failure in axial compression, diametral compression, shear, or bending. There were no significant differences in ultimate stress, strain, or fracture energy between the chemically sterilized and control groups in any of the testing modes.     

Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue

The ability to determine trabecular bone tissue elastic and failure properties has biological and clinical importance. To date, trabecular tissue yield strains remain unknown due to experimental difficulties, and elastic moduli studies have reported controversial results. We hypothesized that the elastic and tensile and compressive yield properties of trabecular tissue are similar to those of cortical tissue. Effective tissue modulus and yield strains were calibrated for cadaveric human femoral neck specimens taken from 11 donors, using a combination of apparent-level mechanical testing and specimen-specific, high-resolution, nonlinear finite element modeling. The trabecular tissue properties were then compared to measured elastic modulus and tensile yield strain of human femoral diaphyseal cortical bone specimens obtained from a similar cohort of 34 donors. Cortical tissue properties were obtained by statistically eliminating the effects of vascular porosity. Results indicated that mean elastic modulus was 10% lower (p<0.05) for the trabecular tissue (18.0+/-2.8 GPa) than for the cortical tissue (19.9+/-1.8 GPa), and the 0.2% offset tensile yield strain was 15% lower for the trabecular tissue (0.62+/-0.04% vs. 0.73+/-0.05%, p<0.001). The tensile-compressive yield strength asymmetry for the trabecular tissue, 0.62 on average, was similar to values reported in the literature for cortical bone. We conclude that while the elastic modulus and yield strains for trabecular tissue are just slightly lower than those of cortical tissue, because of the cumulative effect of these differences, tissue strength is about 25% greater for cortical bone.     

A Bimodular Theory for Finite Deformations: Comparison of Orthotropic Second-order and Exponential Stress Constitutive Equations for Articular Cartilage

Cartilaginous tissues, such as articular cartilage and the annulus fibrosus, exhibit orthotropic behavior with highly asymmetric tensile–compressive responses. Due to this complex behavior, it is difficult to develop accurate stress constitutive equations that are valid for finite deformations. Therefore, we have developed a bimodular theory for finite deformations of elastic materials that allows the mechanical properties of the tissue to differ in tension and compression. In this paper, we derive an orthotropic stress constitutive equation that is second-order in terms of the Biot strain tensor as an alternative to traditional exponential type equations. Several reduced forms of the bimodular second-order equation, with six to nine parameters, and a bimodular exponential equation, with seven parameters, were fit to an experimental dataset that captures the highly asymmetric and orthotropic mechanical response of cartilage. The results suggest that the bimodular second-order models may be appealing for some applications with cartilaginous tissues.     

Femoral head apparent density distribution predicted from bone stresses

A new theory relating bone morphology to applied stress is used to predict the apparent density distribution in the femoral head and neck. Cancellous bone is modeled as a self-optimizing material and cortical bone as a saturated (maximum possible bone density) response to stress in the bone tissue. Three different approaches are implemented relating bone apparent density to: (1) the von Mises stress, (2) the strain energy density in the mineralized tissue and (3) a defined closed effective stress (spherical stress). An iterative nonlinear three-dimensional finite element model is used to predict the apparent density distribution in the femoral head and neck for each of the three approaches. It is shown that the von Mises stress (an open effective stress) cannot accurately predict bone apparent density. It is shown that strain energy density and the defined closed effective stress can predict apparent density and that they give predictions consistent with the observed density pattern in the femoral head and neck.     

Determination of orthotropic bone elastic constants using FEA and modal analysis

Finite element models have been widely employed in an effort to quantify the stress and strain distribution around implanted prostheses and to explore the influence of these distributions on their long-term stability. In order to provide meaningful predictions, such models must contain an appropriate reflection of mechanical properties. Detailed geometrical and density information is now readily available from CT scanning. However, despite the use of phantoms, a method of determining mechanical properties (or elastic constants) from bone density has yet to be made available in a usable form.In this study, a cadaveric bone was CT scanned and its natural frequencies were measured using modal analysis. Using the geometry obtained from the CT scan data, a finite element mesh was created with the distribution of density established by matching the mass of the FE bone model with the mass of the cadaveric bone. The maximum values of the orthotropic elastic constants were then established by matching the predictions from FE modal analyses to the experimental natural frequencies, giving a maximum error of 7.8% over 4 modes of vibration. Finally, the elastic constants of the bone derived from the analyses were compared with those measured using ultrasound techniques. This produced a difference of <1% for both the maximum density and axial Young's Modulus. This study has thereby produced an orthotropic finite element model of a human femur. More importantly, however, is the implication that it is possible to create a valid FE model by simply comparing the FE results with the measured resonant frequency of the CT scanned bone.     

Effect of gamma irradiation on mechanical properties of human cortical bone: influence of different processing methods

The secondary sterilisation by irradiation reduces the risk of infectious disease transmission with tissue allografts. Achieving sterility of bone tissue grafts compromises its biomechanical properties. There are several factors, including dose and temperature of irradiation, as well as processing conditions, that may influence mechanical properties of a bone graft. The purpose of this study was to evaluate the effect of gamma irradiation with doses of 25 or 35?kGy, performed on dry ice or at ambient temperature, on mechanical properties of non-defatted or defatted compact bone grafts. Left and right femurs from six male cadaveric donors aged from 46 to 54?years, were transversely cut into slices of 10?mm height, parallel to the longitudinal axis of the bone. Compact bone rings were assigned to the eight experimental groups according to the different processing method (defatted or non-defatted), as well as gamma irradiation dose (25 or 35?kGy) and temperature conditions of irradiation (ambient temperature or dry ice). Axial compression testing was performed with a material testing machine. Results obtained for elastic and plastic regions of stress-strain curves examined by univariate analysis are described. Based on multivariate analysis it was found that defatting of bone rings had no significant effect on any mechanical parameter studied, whereas irradiation with both doses decreased significantly the ultimate strain and its derivative toughness. The elastic limit and resilience were significantly increased by irradiation with the dose 25?kGy, but not 35?kGy, when the time of irradiation was longer. Additionally, irradiation at ambient temperature decreased maximum load, elastic limit, resilience, and ultimate stress. As strain in the elastic region was not affected, decreased elastic limit resulted in lower resilience. The opposite phenomenon was observed in the plastic region, where in spite of the lower ultimate stress, the toughness was increased due to the increase in the ultimate strain. The results of our study suggest that there may be an association between mechanical properties of bone tissue grafts and the damage process of collagen structure during gamma irradiation. This collagen damage in cortical bone allografts containing water does not depends on the temperature of irradiation or defatting during processing if dose of gamma irradiation does not exceed 35?kGy.     

A model for the nonlinear elastic response of large arteries

The nonlinear elastic response of large arteries subjected to finite deformations due to action of biaxial principal stresses, is described by simple constitutive equations. Generalized measures of strain and stress are introduced to account for material nonlinearity. This also ensures the existence of a strain energy density function. The orthotropic elastic response is described via quasi-linear relations between strains and stresses. One nonlinear parameter which defines the measures of strain and stress, and three elastic moduli are assumed to be constants. The lateral strain parameters (equivalent to Poisson's ratios in infinitesimal deformations) are deformation dependent. This dependence is defined by empirical relations developed via the incompressibility condition, and by the introduction of a fifth material parameter. The resulting constitutive model compares well with biaxial experimental data of canine carotid arteries.     

Apparent Young's modulus of human radius using inverse finite-element method

The ability to assess the elastic and failure properties of cortical bone at the radial diaphysis has a clinical importance. A new generation of quantitative ultrasound (QUS) devices and peripheral quantitative computed tomography (p-QCT) has been developed to assess non-invasively bone material and structural properties at the distal radius. This anatomical site is characterized by a thin cortical thickness that complicates traditional mechanical testing methods on specimens. Until now, mechanical properties of cortical bone at distal radius (e.g., elastic modulus, yield stress and strain) remain rarely studied probably due to experimental difficulties. The present study introduces an inverse finite-element method strategy to measure the elastic modulus and yield properties of human cortical specimens of the radial diaphysis. Twenty millimeter-thick portions of diaphysis were cut from 40 human radii (ages 45-90) for biomechanical test. Subsequently the same portion was modeled in order to obtain a specimen-specific three dimensional finite-element model (3D-FEM). Longitudinal elastic constants at the apparent level and stress characterizations were performed by coupling mechanical parameters with isotropic linear-elastic simulations. The results indicated that the mean apparent Young's modulus for radial cortical bone was 16 GPa (SD 1.8) and the yield stress was 153 MPa (SD 33). Breaking load was 12,946 N (SD 3644), cortical thickness 2.9 mm (SD 0.6), structural effective strain at the yield (epsilon(y)=0.0097) and failure (epsilon(u)=0.0154) load were also calculated. The 3D-FEM strategy described here may help to investigate bone mechanical properties when some difficulties arise from machining mechanical sample.     

An orthotropic viscoelastic material model for passive myocardium: theory and algorithmic treatment

This contribution presents a novel constitutive model in order to simulate an orthotropic rate-dependent behaviour of the passive myocardium at finite strains. The motivation for the consideration of orthotropic viscous effects in a constitutive level lies in the disagreement between theoretical predictions and experimentally observed results. In view of experimental observations, the material is deemed as nearly incompressible, hyperelastic, orthotropic and viscous. The viscoelastic response is formulated by means of a rheological model consisting of a spring coupled with a Maxwell element in parallel. In this context, the isochoric free energy function is decomposed into elastic equilibrium and viscous non-equilibrium parts. The baseline elastic response is modelled by the orthotropic model of Holzapfel and Ogden [Holzapfel GA, Ogden RW. 2009. Constitutive modelling of passive myocardium: a structurally based framework for material characterization. Philos Trans Roy Soc A Math Phys Eng Sci. 367:3445–3475]. The essential aspect of the proposed model is the account of distinct relaxation mechanisms for each orientation direction. To this end, the non-equilibrium response of the free energy function is constructed in the logarithmic strain space and additively decomposed into three anisotropic parts, denoting fibre, sheet and normal directions each accompanied by a distinct dissipation potential governing the evolution of viscous strains associated with each orientation direction. The evolution equations governing the viscous flow have an energy-activated nonlinear form. The energy storage in the Maxwell branches has a quadratic form leading to a linear stress–strain response in the logarithmic strain space. On the numerical side, the algorithmic aspects suitable for the implicit finite element method are discussed in a Lagrangian setting. The model shows excellent agreement compared to experimental data obtained from the literature. Furthermore, the finite element simulations of a heart cycle carried out with the proposed model show significant deviations in the strain field relative to the elastic solution.     

Fatigue of immature baboon cortical bone

Strain-controlled uniaxial fatigue and monotonic tensile tests were conducted on turned femoral cortical bone specimens obtained from baboons at various ages of maturity. Fatigue loading produced a progressive loss in stiffness and an increase in hysteresis prior to failure, indicating that immature primate cortical bone responds to repeated loading in a fashion similar to that previously observed for adult human cortical bone. Bone fatigue resistance under this strain controlled testing decreased during maturation. Maturation was also associated with an increase in bone dry density, ash fraction and elastic modulus. The higher elastic modulus of more mature bone meant that these specimens were subjected to higher stress levels during testing than more immature bone specimens. Anatomical regions along the femoral shaft exhibited differences in strength and fatigue resistance.     

Load transfer along the bone–dental implant interface

In this paper the variation of normal and shear stresses along a path defined on the bone–dental implant interface is investigated. In particular, the effects of implant diameter, collar length and slope, body length, and the effects of four different types of external threads on the interfacial stress distribution are studied. The geometry of the bone is digitized from a CT scan of a mandibular incisor and the surrounding bone. The bone and the implant are assumed to be perfectly bonded. The finite element method with 2D plane strain assumption is used to compute interfacial stresses. Highest continuous interfacial stresses are encountered in the region where the implant collar engages the cortical region, and near the apex of the implant in the subcortical region. Stress concentrations in the interfacial stresses occur near the geometric discontinuities on the implant contour, and jumps in stress values occur where the elastic modulus of the bone transitions between the cortical and trabecular bone values. Among the six contour parameters, the slope and the length of the implant collar, and the implant diameter influence the interfacial stress levels the most, and the effects of changing these parameters are significantly noticed only in the cortical bone (alveolar ridge) area. External threads cause significant stress concentrations in interfacial stresses in otherwise smoothly varying regions. This work shows that the presence of external threads could cause significant variations in both normal and shear stresses along the bone–implant interface, but not reduction in shear stress as previously thought.     

The importance of the elastic and plastic components of strain in tensile and compressive fatigue of human cortical bone in relation to orthopaedic biomechanics

The longevity, success, or failure of an orthopaedic implant is dependent on its osseointegration especially within the initial six months of the initial surgery. The development of strains plays a crucial role in both bone modelling and remodelling. For remodelling, in particular, strains of substantial values are required to activate the osteoblastic and osteoclastic activity for the osseointegration of the implant. Bone, however, is subject to "damage" when strain levels exceed a certain threshold level. Damage is manifested in the form of microcracks; it is linked to increased elastic strain amplitudes and is accompanied by the development of "plastic" (irrecoverable, residual) strains. Such strains increase the likelihood for the implant to subside or loosen. The present study examines the rates (per cycle) by which these two components of strain (elastic and "plastic") develop during fatigue cycling in two loading modes, tension and compression. The results of this study show that these strain rates depend on the applied stress in both loading modes. It also shows that elastic and plastic strain rates can be linked to each other through simple power law relationships so that one can calculate or predict the latter from the former and vice versa. We anticipate that such basic bone biomechanics data would be of great benefit to both clinicians and bioengineers working in the field of FEA modelling applications and orthopaedic implant surgery.