Finite element and experimental models of cemented hip joint reconstructions can produce similar bone and cement strains in pre-clinical tests

Finite element (FE) models could be used for pre-clinical testing of cemented hip replacement implants against the damage accumulation failure scenario. To accurately predict mechanical failure, the models should accurately predict stresses and strains. This should be the case for various implants. In the current study, two FE models of composite hip reconstructions with two different implants were validated relative to experimental bone and cement strains. The objective was an overall agreement within 10% between experimental and FE strains. Two stem types with different clinical results were analyzed: the Lubinus SPII and the Mueller Curved with loosening rates of 4% and 16% after 10 yr, respectively (Prognosis of total hip replacement. 63rd Annual Meeting of the American Academy of orthopaedic surgeons, Atlanta, USA). For both implant types, six stems were implanted in composite femurs. All specimens were subjected to bending. The Mueller Curved specimens were additionally subjected to torsion. Bone strains were recorded at 10 locations on the cortex and cement strains at three locations within the cement mantle. An FE model was built for both stem types and the experiments were simulated. Bone and cement strains were calculated at the experimental gauge locations. Most FE bone strains corresponded to the mean experimental strains within two standard deviations; most FE cement strains within one standard deviation. Linear regression between the FE and mean experimental strains produced slopes between 0.82 and 1.03, and R(2) values above 0.98. Particularly for the Mueller Curved, agreement improved considerably when FE strains were compared to the strains from the experimental specimen used to build the FE model. The objective of overall agreement within 10% was achieved, indicating that both FE models were successfully validated. This prerequisite for accurately predicting long-term failure has been satisfied.     

Finite element analysis of a three-dimensional open-celled model for trabecular bone

Based on a regular array of cubic unit cells, each containing a body-centered spherical void, we created an idealized three-dimensional model for both subchondral trabecular bone and a class of porous foams. By considering only face-to-face stacking of unit cells, the inherent symmetry was such that, except at the surface, the displacements and stresses within any one unit cell were representative of the entire porous structure. Using prescribed displacements the model was loaded in both uniaxial compressive strain and uniaxial shear strain. Based on the response to these loads, we found the tensor of elastic constants for an equivalent homogeneous elastic solid with cubic symmetry. We then compared the predicted modulus with our experimental values for bovine trabecular bone and literature values for an open-celled latex rubber foam.     

Design and analysis of robust total joint replacements: finite element model experiments with environmental variables

Computer simulation of orthopaedic devices can be prohibitively time consuming, particularly when assessing multiple design and environmental factors. Chang et al. (1999) address these computational challenges using an efficient statistical predictor to optimize a flexible hip implant, defined by a midstem reduction, subjected to multiple environmental conditions. Here, we extend this methodology by: (1) explicitly considering constraint equations in the optimization formulation, (2) showing that the optimal design for one environmental distribution is robust to alternate distributions, and (3) illustrating a sensitivity analysis technique to determine influential design and environmental factors. A thin midstem diameter with a short stabilizing distal tip minimized the bone remodeling signal while maintaining satisfactory stability. Hip joint force orientation was more influential than the effect of the controllable design variables on bone remodeling and the cancellous bone elastic modulus had the most influence on relative motion, both results indicating the importance of including uncontrollable environmental factors. The optimal search indicated that only 16 to 22 computer simulations were necessary to predict the optimal design, a significant savings over traditional search techniques.     

Finite element modeling of trabecular bone damage

This paper presents a finite element-based, computational model for analysis of structural damage to trabecular bone tissues. A modulus reduction method was formulated from elasto-plasticity theory, and was used to account for site-specific trabecular bone tissue damage. Trabecular bone tissue damage is illustrated using a large-scale, anatomically accurate, two-dimensional, microstructural finite element model of a human thoracic vertebral body. Four models with varying specifications for damage accumulation were subjected to compressive loading and unloading cycles. The numerical results and experimental validation demonstrated that the modulus reduction method reproduced the non-linear mechanical behaviour of vertebal trabecular bone. The iterative computational approach presented provides a methodology to study trabecular bone damage, and should provide researchers with a computational approach to study bone fracture and repair and to predict vertebral fragility.     

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.     

Finite element analysis of the impact of bone nanostructure on its piezoelectric response

Biomechanics and Modeling in Mechanobiology - The piezoelectric response of bone at the submicron scale is analyzed under mechanical loadings using the finite element (FE) method. A new algorithm...     

Finite element analysis of covered microstents

Currently available neuroendovascular devices are inadequate for effective treatment of many wide-necked or fusiform intracranial aneurysms and intracranial carotid-cavernous fistulae (CCF). Placing a covered microstent across the intracranial aneurysm neck and CCF rent could restore normal vessel morphology by preventing blood flow into the aneurysm lumen or CCF rent. To fabricate covered microstents, our research group has developed highly flexible ultra thin (approximately 150 microm) silicone coverings and elastomerically captured them onto commercially available metal stents without stitching. Preliminary in vivo studies were conducted by placing these covered microstents in the common carotid artery of rabbits. The feasibility of using covered stents was demonstrated. However, the cover affected the deployment pressure and the stents failed occasionally during deployment due to tearing of the cover. Appropriate modeling of covered stents will assist in designing suitable coverings, and help to reduce the failure rate of covered microstents. The purpose of this study is to use the finite element method to determine the mechanical properties of the covered microstent and investigate the effects of the covering on the mechanical behavior of the covered microstent. Variations in the mechanical properties of the covered microstent such as deployment pressure, elastic recoil and longitudinal shortening due to change in thickness and material properties of the cover have been investigated. This work is also important for custom design of covered microstents such as adding cutout holes to save adjacent perforating arteries.     

Finite element simulations of a hip joint with femoroacetabular impingement

In this study, a three-dimensional finite element (FE) model based on the specific anatomy of a patient presenting a femoroacetabular impingement of the ‘cam’-type is developed. The FE meshes of the structures of interest are obtained from arthrographic magnetic resonance images. All soft tissues are considered linear elastic and isotropic, and the bones were assumed rigid. A compression of the femur on the acetabular cavity as well as flexural movements and internal rotations are applied. Stresses and contact pressures are evaluated in this patient-specific model in order to better interpret the mechanism of aggression of the femoral and acetabular cartilages. The corresponding results are presented and discussed. The values obtained for the contact pressures are similar to those reported by other models based on idealised geometries. An FE analysis of a non-cam hip is also performed for comparison with the pathological case.     

Declining needs for total joint replacements for rheumatoid arthritis

This millennium brings new views to rheumatology. Total joint replacement surgery is needed less often as active treatment strategies combined with availability of new medications has led to more effective rheumatoid arthritis control. This was beautifully shown in a recent issue of Arthritis Research & Therapy by a Swedish study that uses data from national registers and compares incidence rates for total hip and knee arthroplasties before and after the establishment of biologic agents use for rheumatoid arthritis     

Finite element modelling of polymethylmethacrylate flow through cancellous bone

A mathematical model based on the Finite Element Method is developed to simulate the non-linear flow of acrylic bone cement through cancellous bone. The cancellous bone bed is modelled as a bed of parallel capillaries filled with equal spaced toroidal trabeculae. By manipulating the relative size of the torus and the capillary, the flow within bone of varying porosity is simulated. An apparent permeability based on the volume weighted average viscosity and Darcy's law is developed to describe the flow of the acrylic through the cancellous bone bed. The model predicts a cancellous bone permeability of 5.6 x 10(-9)-8.3 x 10(-9) m2 for linear flow. The non-linear behavior of the acrylic cement results in an increase of apparent permeability when compared to the permeability computed for linear flow. Estimates of penetration are achieved by running the model in a quasi-steady state fashion with pressure applied over a fixed time increment. Close agreement is shown between model predictions of penetration depth and experimental results available in the literature.     

Finite element prediction of fatigue damage growth in cancellous bone

Cyclic stresses applied to bones generate fatigue damage that affects the bone stiffness and its elastic modulus. This paper proposes a finite element model for the prediction of fatigue damage accumulation and failure in cancellous bone at continuum scale. The model is based on continuum damage mechanics and incorporates crack closure effects in compression. The propagation of the cracks is completely simulated throughout the damaged area. In this case, the stiffness of the broken element is reduced by 98% to ensure no stress-carrying capacities of completely damaged elements. Once a crack is initiated, the propagation direction is simulated by the propagation of the broken elements of the mesh. The proposed model suggests that damage evolves over a real physical time variable (cycles). In order to reduce the computation time, the integration of the damage growth rate is based on the cycle blocks approach. In this approach, the real number of cycles is reduced (divided) into equivalent blocks of cycles. Damage accumulation is computed over the cycle blocks and then extrapolated over the corresponding real cycles. The results show a clear difference between local tensile and compressive stresses on damage accumulation. Incorporating stiffness reduction also produces a redistribution of the peak stresses in the damaged region, which results in a delay in damage fracture.     

Finite element modeling for strain rate dependency of fracture resistance in compact bone

Crack growths in compact bones driven by various strain rate levels were studied using finite element modeling. The energy resistance curves in bovine femur cortical bones were characterized, whereas the orthotropic viscoelasticity in bone materials was accounted for to assess the effect of strain rate on the energy resistance curve. The models were also used to justify the anticipated plane strain response as a result of rather thick specimens used in experiments. Similarities were found between the experimental and model results when crack resistance ability exhibited in bones with slow loading rates, while unstable crack growth existed in bones with rapid loading rates. The critical energy release rates slightly decreased with the increase in strain rates. The hybrid experimental and computational method introduced in this study could be beneficial for application in fracture study in which standard experiments cannot be validly performed.     

Finite element analysis of ramming in Ovis canadensis

The energy produced during the ramming of bighorn sheep (Ovis canadensis) would be expected to result in undesirable stresses in their frontal skull, which in turn would cause brain injury; yet, this animal seems to suffer no ill effects. In general, horn is made of an α-keratin sheath covering a bone. Despite volumes of data on the ramming behavior of Ovis canadensis, the extent to which structural components of horn and horn-associated structure or tissue absorb the impact energy generated by the ramming event is still unknown. This study investigates the hypothesis that there is a mechanical relationship present among the ramming event, the structural constituents of the horn, and the horn-associated structure. The three-dimensional complex structure of the bighorn sheep horn was successfully constructed and modeled using a computed tomography (CT) scan and finite element (FE) method, respectively. Three different three-dimensional quasi-static models, including a horn model with trabecular bone, a horn model with compact bone that instead of trabecular bone, and a horn model with trabecular bone as well as frontal sinuses, were studied. FE simulations were used to compare distributions of principal stress in the horn and the frontal sinuses and the strain energy under quasi-static loading conditions. It was noticed that strain energy due to elastic deformation of the complex structure of horn modeled with trabecular bone and with trabecular bone and frontal sinus was different. In addition, trabecular bone in the horn distributes the stresses over a larger volume, suggesting a mechanical link between the structural constituents and the ramming event. This phenomenon was elucidated through the principal stress distribution in the structure. This study will help designers in choosing appropriate material combinations for the successful design of protective structures against a similar impact.     

Finite element scaling analysis of human craniofacial growth

The study of form change is central to traditional cephalometric research. Unfortunately, traditional cephalometric studies operate within systems of measurement that are based on registration and orientation. Measurements produced in registered systems are insufficient for the craniofacial biologist who is interested in locating morphological differences between forms. In this article we apply a registration-free method called finite element scaling analysis in a study of the form change occurring during growth of the normal human craniofacial complex. The method provides form change data that can be summarized at various morphological levels. Twenty normal male individuals are used to analyze the form change that occurs from age 4 to ages 5, 7, 8, 9, 10, 12, 13, and 15 years. The magnitude and direction of growth expressed as shape and size change specific to craniofacial landmarks are presented. Although exceptions occur, our analysis shows that localized size change is, on the average, greater than localized shape change. The relation between size and shape change during growth shows allometry (shape change increasing during growth along with size change) but at a lesser magnitude and slower rate. We conclude that although shape change occurs throughout ontogeny, the magnitude and rate of shape change in relation to size change diminishes as age increases. This analysis represents new insights into the understanding of human craniofacial growth at various levels of morphological integration.     

Finite element analysis of an additional implant for intramedullary nailing]

In the surgical treatment of fractured femurs, the fracture is bridged by a medullary nail fixed in the bone with interlocking screws. Failure of bone substance in the region of the interlocking screws is the most common complication in the treatment of osteoporotic bone. With the aim of preventing this complication, an additional implant was developed. A finite element analysis of an ideal bone/implant system was carried out to investigate the role of the additional implant. Three defined finite element models were generated, and the associated stress situations compared. The first model is a standard fixation without the additional implant. In the second model, the additional implant is integrated within the bone/implant system. The third model uses a modified form of the additional implant. The results show that both additional implants reduce the stresses occurring, both in the bone substance and at the screws. The modified form of the additional implant proved to be the most favorable version. In the case of the original additional implant, the negative effect of the sharp edges of the thread was demonstrable.     

Finite element analysis of heel pad with insoles

To design optimal insoles for reduction of pedal tissue trauma, experimental measurements and computational analyses were performed. To characterize the mechanical properties of the tissues, indentation tests were performed. Pedal tissue geometry and morphology were obtained from magnetic resonance scan of the subject's foot. Axisymmetrical finite element models of the heel of the foot were created with 1/4 of body weight load applied. The stress, strain and strain energy density (SED) fields produced in the pedal tissues were computed. The effects of various insole designs and materials on the resulting stress, strain, and SED in the soft pedal tissues were analyzed. The results showed: (a) Flat insoles made of soft material provide some reductions in the maximum stress, strain and SED produced in the pedal tissues. These maximum values were computed near the calcaneus. (b) Flat insoles, with conical/cylindrical reliefs, provided more reductions in these maximum values than without reliefs. (c) Custom insoles, contoured to match the pedal geometry provide most reductions in the maximum stress, strain and SED. Also note, the maximum stress, strain and SED computed near the calcaneus were found to be about 10 times the corresponding peak values computed on the skin surface. Based on the FEA analysis, it can be concluded that changing insole design and using different material can significantly redistribute the stress/strain inside the heel pad as well as on the skin surface.     

Finite element study of a tissue-engineered cartilage transplant in human tibiofemoral joint

Most tissue-engineered cartilage constructs are more compliant than native articular cartilage (AC) and are poorly integrated to the surrounding tissue. To investigate the effect of an implanted tissue-engineered construct (TEC) with these inferior properties on the mechanical environment of both the engineered and adjacent native tissues, a finite element study was conducted. Biphasic swelling was used to model tibial cartilage and an implanted TEC with the material properties of either native tissue or a decreased elastic modulus and fixed charged density. Creep loading was applied with a rigid impermeable indenter that represented the femur. In comparison with an intact joint, compressive strains in the transplant, surface contact stress in the adjacent native AC and load partitioning between different phases of cartilage were affected by inferior properties of TEC. Results of this study may lead to a better understanding of the complex mechanical environment of an implanted TEC.