The effect of stress concentration on bone remodeling: theoretical predictions

Theoretical predictions of internal bone remodeling around an elliptical hole are studied. The internal remodeling theory due to Cowin and Hegedus is employed. The bone is modeled as an initially homogeneous adaptive elastic plate with an elliptical hole under a superposed steady compressive load. It is shown that there will exist a final inhomogeneous remodeling distribution around the hole that will disappear away from the hole. The remodeling is such that the compressive stress concentration causes the bone structure to evolve to one of greater density and stiffer elastic coefficients. The speed of remodeling around the hole and its variation with respect to distance is investigated and discussed. It is shown that the rate of bone reinforcement in the area of compressive stress concentration is much higher than the rate of bone resorption in the area of existing tensile stress. Special cases of a circular hole and vertical and horizontal slots are studied and discussed.     

Understanding site-specific residual strain and architecture in bovine cortical bone

Living bone is considered as adaptive material to the mechanical functions, which continually undergoes change in its histological arrangement with respect to external prolonged loading. Such remodeling phenomena within bone depend on the degree of stimuli caused by the mechanical loading being experienced, and therefore, are specific to the sites. In the attempts of understanding strain adaptive phenomena within bones, different theoretical models have been proposed. Also, the existing literatures mostly follow the measurement of surface strains using strain gauges to experimentally quantify the strains experienced in the functional environment. In this work, we propose a novel idea of understanding site-specific functional adaptation to the prolonged load in bone on the basis of inherited residual strains and structural organization. We quantified the residual strains and amount of apatite crystals distribution, i.e., the degree of orientation, using X-ray diffraction procedures. The sites of naturally existing hole in bone, called foramen, are considered from bovine femur and metacarpal samples. Significant values of residual strains are found to exist in the specimens. Trends of residual strains noted in the specimens are mostly consistent with the degree of orientation of the crystallites. These features explain the response behavior of bone to the mechanical loading history near the foramen sites. Preferential orientation of crystals mapped around a femoral foramen specimen showed furnished tailored arrangement of the crystals around the hole. Effect of external loading at the femoral foramen site is also explained by the tensile loading experiment.     

Predicting bone remodeling around tissue- and bone-level dental implants used in reduced bone width

The objective of this study was to predict time-dependent bone remodeling around tissue- and bone-level dental implants used in patients with reduced bone width. The remodeling of bone around titanium tissue-level, and titanium and titanium–zirconium alloy bone-level implants was studied under 100 N oblique load for one month by implementing the Stanford theory into three-dimensional finite element models. Maximum principal stress, minimum principal stress, and strain energy density in peri-implant bone and displacement in x- and y- axes of the implant were evaluated. Maximum and minimum principal stresses around tissue-level implant were higher than bone-level implants and both bone-level implants experienced comparable stresses. Total strain energy density in bone around titanium implants slightly decreased during the first two weeks of loading followed by a recovery, and the titanium–zirconium implant showed minor changes in the axial plane. Total strain energy density changes in the loading and contralateral sides were higher in tissue-level implant than other implants in the cortical bone at the horizontal plane. The displacement values of the implants were almost constant over time. Tissue-level implants were associated with higher stresses than bone-level implants. The time-dependent biomechanical outcome of titanium–zirconium alloy bone-level implant was comparable to the titanium implant.     

Reasons why dynamic compression plates are inferior to locking plates in osteoporotic bone: a finite element explanation

While locking plate fixation is becoming increasingly popular for complex and osteoporotic fractures, for many indications compression plating remains the standard choice. This study compares the mechanical behaviour of the more recent locking compression plate (LCP) device, with the traditional dynamic compression plates (DCPs) in bone of varying quality using finite element modelling. The bone properties considered include orthotropy, inhomogeneity, cortical thinning and periosteal apposition associated with osteoporosis. The effect of preloads induced by compression plating was included in the models. Two different fracture scenarios were modelled: one with complete reduction and one with a fracture gap. The results show that the preload arising in DCPs results in large principal strains in the bone all around the perimeter of the screw hole, whereas for LCPs large principal strains occur primarily on the side of the screw proximal to the load. The strains within the bone produced by the two screw types are similar in healthy bone with a reduced fracture gap; however, the DCP produces much larger strains in osteoporotic bone. In the presence of a fracture gap, the DCP results in a considerably larger region with high tensile strains and a slightly smaller region with high compressive strains. These findings provide a biomechanical basis for the reported improved performance of locking plates in poorer bone quality.     

Non-interacting modes for stress, strain and energy in anisotropic hard tissue.

The six non-interacting modes for stress, strain and energy in an orthotropic elastic model of human femoral cortical bone tissue are discussed and illustrated. The stress and strain modes are illustrated using the representation of the stress and strain fields around a circular hole in a flat plate of cortical bone subjected to a uniaxial field of tension as the example. The six modes play a role in the stress analysis of orthotropic elastic materials similar to the roles played by the hydrostatic and deviatoric non-interacting stress, strain and energy modes in isotropic elasticity. The biomechanical significance of the six non-interacting modes for stress, strain and energy in hard tissue is both practical and suggestive. The modes suggest a practical scheme for the representation of stress and strain fields in hard tissue. The existence of the modes suggests physical insights, for example, possible failure mechanisms or adaptation strategies possessed by the hard tissues.     

Numerical analysis of an osseointegrated prosthesis fixation with reduced bone failure risk and periprosthetic bone loss

Currently available implants for direct attachment of prosthesis to the skeletal system after transfemoral amputation (OPRA system, Integrum AB, Sweden and ISP Endo/Exo prosthesis, ESKA Implants AG, Germany) show many advantages over the conventional socket fixation. However, restraining biomechanical issues such as considerable bone loss around the stem and peri-prosthetic bone fractures are present. To overcome these limiting issues a new concept of the direct intramedullary fixation was developed. We hypothesize that the new design will reduce the peri-prosthetic bone failure risk and adverse bone remodeling by restoring the natural load transfer in the femur. Generic CT-based finite element models of an intact femur and amputated bones implanted with 3 analyzed implants were created and loaded with a normal walking and a forward fall load. The strain adaptive bone remodeling theory was used to predict long-term bone changes around the implants and the periprosthetic bone failure risk was evaluated by the von Mises stress criterion. The results show that the new design provides close to physiological distribution of stresses in the bone and lower bone failure risk for the normal walking as compared to the OPRA and the ISP implants. The bone remodeling simulations did not reveal any overall bone loss around the new design, as opposed to the OPRA and the ISP implants, which induce considerable bone loss in the distal end of the femur. This positive outcome shows that the presented concept has a potential to considerably improve safety of the rehabilitation with the direct fixation implants.     

Periosteal and endosteal control of bone remodeling under torsional loading.

The shape changes that occur in the mid-diaphysis of a long bone due to adaptive remodeling induced by increasing or decreasing the axial and/or torsional loading of the bone are investigated using a simple model. In this model the mid-diaphysis of a long bone is represented as a hollow thick-walled right-circular cylinder, and different optimal strategies for bone remodeling are considered. It is shown that if such a thick-walled right-circular cylinder capable of surface remodeling is subjected to an axial compressive load and a twisting torque, then the remodeling patterns depend on whether the periosteal surface or the endosteal surface controls the limits of the remodeling process. It is shown that the effect of increasing the torque is always opposite to the effect of increasing the compressive load. Thus, similar remodeling patterns are obtained by increasing one type of loading and decreasing the other. Aside from the restriction of idealized cylindrical geometry, the only assumptions made are that the bone tissue is linearly elastic and that there exists a finite range of remodeling equilibrium stresses. Only those loading situations which maintain the bone in remodeling equilibrium are considered in this work. It follows that the results presented are independent of the specific type of rule governing the temporal evolution of the bone shape, since any such rule applies only in situations where there is active remodeling and, hence, no remodeling equilibrium.     

Multi-factorial analysis of variables influencing the bone loss of an implant placed in the maxilla: Prediction using FEA and SED bone remodeling algorithm

The aim of this study was to investigate the interactions of implant position, implant–abutment connection and loading condition influencing bone loss of an implant placed in the maxilla using finite element (FE) analysis and mathematical bone remodeling theory. The maxilla section contours were acquired using CT images to construct FE models containing RS (internal retaining-screw) and the TIS (taper integrated screwed-in) implants placed in SC (along the axis of occlusal force) and RA (along the axis of residual ridge) positions. The adaptive strain energy density (SED) algorithm was combined with FE approach to study the preliminary bone remodeling around implant systems under different load conditions. The simulated results showed that the implant position obviously influenced the bone loss. An implant placed in the RA position resulted in substantially increased bone loss. Implant receiving a lateral load slightly increased bone loss compared with an axial load. The implant type did not significantly influence bone loss. It was found that buccal site suffered the most bone loss around the implant, followed by distal, lingual and mesial sites. The implant position primarily influenced bone loss and it was found most obviously at the buccal site. Implant placed along the axial load direction of a proposed prosthesis could obtain less bone loss around the implant. Attaining proper occlusal adjustments to reduce the lateral occlusal force is recommended in implant–bone–prosthesis system. Abutments of internal engagement with or without taper-fit did not affect the bone loss in the surrounding bone.     

Trabecular bone remodeling around smooth and porous implants in an equine patellar model

The objective of this investigation was to examine the stress-morphology relationships for trabecular bone around implants with different surface characteristics. Stainless steel spheres with either a polished surface or a sintered-bead porous coating were implanted unilaterally into equine patellae and maintained for a 6 month period. Stereological methods were used to quantify the trabecular bone morphology and finite element analyses were performed to predict the trabecular bone stresses. In general, the remodeling response around the smooth implants was greater than that around those porous implants that exhibited bone ingrowth. In accordance with these differences, the finite element models predicted greater changes in the stresses adjacent to the smooth implants due to the nonlinear boundary conditions. However, it did not appear that the trajectorial theory, in its simplest form, was applicable to the remodeling induced by the implants. A linear relationship between the change in bone areal density and the change in von Mises effective stress provides support for the hypothesis that the architecture of trabecular bone corresponds to an optimal structure. The results also demonstrated that, under certain circumstances, small changes in the stress state may result in large changes in the principal material orientation.     

Finite element modelling of implant designs and cortical bone thickness on stress distribution in maxillary type IV bone

The aims of this study were to examine the effect of implant neck design and cortical bone thickness using 3D finite element analysis and to analyse the stability of clinical evidence based on micromotion and principal stress. Four commercial dental implants for a type IV bone and maxillary segments were created. Various parameters were considered, including the osseointegration condition, loading direction and cortical bone thickness. Micromotion and principal stresses were used to evaluate the failure of osseointegration and bone overloading, respectively. It was found that the maximum stress of the peri-implant bone decreased as cortical bone thickness increased. The micromotion level in full osseointegration is less than that in non-osseointegration and it also decreases as cortical bone thickness increases. The cortical bone thickness should be measured before surgery to help select a proper implant. In the early stage of implantation, the horizontal loading component induces stress concentration in bone around the implant neck more easily than does the vertical loading component, and this may result in crestal bone loss.     

Finite element evaluation of the AIA shear specimen for bone

An elastic-plastic finite element analysis is performed on the AIA shear specimen to evaluate its effectiveness to yield ultimate shear strength values. The effect of geometry, material properties, and yield criteria are discussed in the light of applications to human femoral cortical bone. Specimen dimensions are noted as follows: W, width, D, hole diameter and H, distance between holes. As the H/D ratio increases the stress distribution tends more toward pure shear at the same time the overshoot in the shear distribution increases. An H/D ratio equal to 1.2-1.5 is optimal. The H/W parameter does not affect the overshoot noticeably but it does slightly affect the purity of shear. The material parameters do affect the performance of the shear specimen. However, the effect of the material parameters are far more pronounced in the anisotropic case than it is in the isotropic case. In the isotropic case, the Young modulus does not affect the overshoot. The increase in Poisson's ratio does slightly decrease the overshoot. For the anisotropic case, the increase in the ratio of shear modulus to Young modulus in the transverse direction (G/E2) results in an increase in the overshoot (in the shear distribution). The increase in the ratio of the Young modulus in the transverse direction to that of the axial direction (E2/E1) also results in an increase in the overshoot. Creating a notch at the top of the hole is shown to have the effect of decreasing the overshoot. Its effect on the purity of the shear is rather slight. It is found that plasticity is initiated at the sides of the two holes where the tensile normal stresses are maximum. The plastic region first expands around the perimeter of the hole then radially outward; and finally, it expands into the significant region. If the W/H parameter is less than 5, a sizable portion of the width of the specimen around the hole can go plastic with the significant region still being in the elastic state. Such a situation can cause tearing of the specimen across the width. A W/H ratio of 6 or more can prevent that danger. It is also found that the onset of plasticity brings about higher overshoot and higher purity of shear. The notched shear specimen performs better in actual tests and is more reliable in producing shear failures. The shear strength results obtained from AIA shear tests tend to confirm those shear strength results obtained from torsion tests.     

A three-dimensional finite element stress analysis for tunnel placement and buttons in anterior cruciate ligament reconstructions

This communication reports the results of a three-dimensional finite element (FE) model of stresses in a surgically altered femur and tibia. The model incorporated a novel approach in implementing orthotropic and inhomogeneous bone properties and non-uniform distributed loading. Cortical, cancellous, and subchondral bone of the femur and the tibia were modeled. Mechanical properties for the cortical and cancellous bone were mapped from published data characterizing the anisotropy and inhomogeneity of the bone properties. Mesh adequacy was determined using stress convergence and strain energy error convergence. Qualitatively, the results of the study compare well with experimental principal compressive strains from the literature. With respect to tunnel placement in anterior cruciate ligament reconstruction, the model predicted stress-shielding at the postero-lateral region of the tunnel wall, and increased stress at the postero-medial region of the tunnel wall. The stresses in the cancellous bone beneath the tunnel were, in general, lower than those above the tunnel. Prolonged stress shielding leads to bone resorption of the posterior tunnel wall leading to tunnel enlargement, and possible compromise of the ACL reconstruction. The stresses on the femoral cortex produced from a button-type fixation were noticeable for low levels of loading; the stress levels were very similar in models incorporating bone properties of patients aged 45 and 65. Repeated compression of the femoral cortex at these stress levels may cause microdamage to the cortex eventually resulting in fatigue failure.     

The open section effect in a long bone with a longitudinal defect - a theoretical modeling study

A longitudinal defect dramatically alters the stress distribution within a long bone. The altered stress distribution can influence the structural properties of the bone and the stimulus for repair and remodeling of the defect and the surrounding bone. For applied torsion, the defect interrupts the normal shear flow around the bone. Reversal of the shear flow along the inner cortex of the bone is the primary characteristic of the "open-section" effect. Stress concentration effects also produce large stresses at the defect corners. A finite element model of a femur mid-diaphysis with a rectangular defect in the posterior cortex was developed to quantify the femur stress distribution and torsional stiffness for defect widths ranging from one-tenth of the femur outer diameter (0.1 OD) to 0.3 OD, and defect lengths ranging from 0.5 to 5 OD. Defects with a length of 1 OD or shorter had little influence on the femur torsional stiffness or the femur shear-stress distribution. The torsional stiffness decreased most dramatically as the defect length increased from 2 to 3 OD, but began to approach an asymptote near 5 OD. Shear flow reversal peaked at the center of the defect for defects longer than 1 OD, and the magnitude of the reversal began to approach an asymptote near 5 OD. For each defect, the largest stresses within the bone, developed at the defect corners. The results indicate that the open-section effect decreases the torsional stiffness and stress concentration effects decrease the torsional strength of a long bone with a longitudinal defect.     

Stress relaxation behaviour of trabecular bone specimens

The present study defines several conditions under which stress relaxation tests can be performed and investigates the viscoelastic behaviour of trabecular bone in compression through a series of stress relaxation tests at three strain levels and in three loading directions of each cubic specimen. A visoelastic model is proposed to characterize the behaviour of trabecular bone and a spectrum of relaxation times is determined. Trabecular bone from the femoral head is non-linearly viscoelastic and displays anisotropic behaviour, which cannot be more symmetric elastically than orthotropic.     

Elastic constants of composites formed from PMMA bone cement and anisotropic bovine tibial cancellous bone

An ultrasonic pulse-transit time technique is used to determine the nine orthotropic engineering constants of 32 cement-cancellous bone composites as a function of volume fractions of bone ranging from 0.0 to 0.4. The composites are manufactured using well-aligned bovine cancellous bone from the proximal end of the tibia and low viscosity bone cement. Selected composites are also subjected to mechanical compression tests to compare with the ultrasonic results. There is excellent correlation between the dynamic or ultrasonically determined moduli and the static or mechanically determined moduli; the dynamic moduli are approximately twice the static moduli and this difference is thought to be due to the effect of strain rate. An orthotropic model is assumed requiring nine independent elastic constants to be determined. The dynamic Young's modulus in the direction of major trabecular alignment, E1, increases linearly from 4.9 to 10.4 GPa as bone volume fraction increases from 0 to 0.4; dynamic E2 and E3 values increase from 4.9 to 7 GPa as bone volume fractions increase from 0 to 0.4, with E2 being slightly higher than E3. The dynamic shear modulus, G12, increases from 1.8 to 3.0 GPa, and G31 and G23 increase slightly from 1.8 to 2.2 GPa as bone volume fractions increase from 0 to 0.4. The Poisson's ratios are more sensitive than the Young's moduli and shear moduli to experimental error in the velocity measurements. The mechanically tested modulus (static modulus) in the direction of major trabecular alignment, E1, increases with volume fraction of bone from 2.4 to 4.4 GPa as the bone volume fraction increases from 0 to 0.25; static E2 and E3 values are either equal to or lower than that of pure PMMA.     

Prediction of microdamage formation using a mineral-collagen composite model of bone

Age-related changes in bone quality are mainly manifested in the reduced toughness. Since the post-yield deformation of bone is realized through microdamage formation (e.g., microcracking and diffuse damage), it is necessary to understand the mechanism of microdamage formation in bone in order to elucidate underlying mechanisms of age-related bone fractures. In this study, a two-dimensional shear lag model was developed to predict stress concentration fields around an initial crack in a mineral-collagen composite. In this model, non-linear elasticity was assumed for the collagen phase, and linear elasticity for the mineral. Based on the pattern of the stress concentration fields, the condition for microdamage formation was discussed. The results of our analyses indicate that: (1) an initial crack formed in mineral phase may cause stress concentration in the adjacent mineral layers; (2) the pattern of stress concentration fields depends not only on the spatial but also mechanical properties of the collagen and mineral phases; (3) the pattern of the stress concentration fields could determine either coalescence or scattering of nano cracks around the initial crack.     

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.     

Predictions of bone remodeling around dental implant systems

This study presents the implementation of a mathematical bone remodeling algorithm to bone adaptation in the premolar area of the mandible around various dental implant systems, and thus sheds a new perspective to the complex interactions in dental implant mechanics. A two-dimensional, plane strain model of the bone was built from a CT-scan. The effect of implant contour on internal bone remodeling was investigated by considering four dental implant systems with contours similar to commercially available ones and another four with cylindrical and conical cross-sections. The remodeling algorithm predicts non-homogeneous density/elastic modulus distribution; and, implant contour has some effect on how this is distributed. Bone density is predicted to increase on the tips of the threads of the implants, but to decrease inside the grooves. Threadless implants favor to develop a softer bone around their periphery, compared to implant systems that have threads. The overall contour (dimensions and the shape) of an implant affect the bone density redistribution, but the differences between different implant systems are relatively small.     

On stiffness of scaffolds for bone tissue engineering—a numerical study

Tissue scaffolds are typically designed and fabricated to match native bone properties. However, it is unclear if this would lead to the best tissue ingrowth outcome within the scaffold as neo-tissue keeps changing the stiffness of entire construct. This paper presents a numerical method to address this issue for design optimization and assessment of tissue scaffolds. The elasticity tensors of two different types of bones are weighted by different multipliers before being used as the targets in scaffold design. A cost function regarding the difference between the effective elasticity tensor, calculated by the homogenization technique, and the target tensor, is minimized by using topology optimization procedure. It is found that different stiffnesses can lead to different remodeling results. The comparison confirms that bone remodeling is at its best when the scaffold elastic tensor matches or is slightly higher than the elastic properties of the host bone.     

A density distribution algorithm for bone incorporating local orthotropy,modal analysis and theories of cellular solids

An algorithm for bone remodeling is presented which allows for both a redistribution of density and a continuous change of principal material directions for the orthotropic material properties of bone. It employs a modal analysis to add density for growth and a local effective strain based analysis to redistribute density. General re-distribution functions are presented. The model utilizes theories of cellular solids to relate density and strength. The code predicts the same general density distributions and local orthotropy as observed in reality.