Finite element study of human lumbar disc nucleus replacements

Currently, there are a number of nucleus replacements under development. The important concern is how well these implants duplicate the mechanical function of the native nucleus. This finite element model study aimed to investigate the influence of different nucleus replacements on the mechanical response of the disc. Models included partial, full, over-sized, partially saturated, elastic and poroelastic solid replacements. Over-sized nucleus replacements up to 25% yielded results that were comparable to those in the intact state. Differences were much greater in cases with under-sized nucleus replacements. The effect was most pronounced for the 75% under-sized replacement that resembled the condition with a full nucleotomy. Nucleus implants with elastic properties substantially altered load transmission when 10% under-sized and over-sized replacements were considered. Compared to intact, the under-sized implants should be avoided when using biphasic materials with properties similar to the native nucleus, whereas for elastic replacements both under- and over-sized implants should not be used.     

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.     

Finite element analysis of Stryker Xia pedicle screw in artificial bone samples with and without supplemental cement augmentation

A validated, using in vitro biomechanical testing, finite element model was used to evaluate the affects of (1) cement augmentation and (2) an intact posterior cortex in osteoporotic bone. The presence of augmentation and/or a posterior cortical cortex increased the stabilization of the pedicle screw 2–5 fold. Placement of cement influenced failure load and toggle; with distal placement having the largest increase in failure load and decrease in cephalad–caudad toggle. The presence of posterior cortex caused a decrease in the amount of toggle, a proximal shift of the center of rotation and an increase in the maximum failure force.     

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 analysis of the temporomandibular joint during lateral excursions of the mandible

One of the most significant characteristics of the temporomandibular joint (TMJ) is that it is in fact composed of two joints. Several finite element simulations of the TMJ have been developed but none of them analysed the different responses of its two sides during nonsymmetrical movement. In this paper, a lateral excursion of the mandible was introduced and the biomechanical behaviour of both sides was studied. A three-dimensional finite element model of the joint comprising the bone components, both articular discs, and the temporomandibular ligaments was used. A fibre-reinforced porohyperelastic model was introduced to simulate the behaviour of the articular discs, taking into account the orientation of the fibres in each zone of these cartilage components. The mandible movement during its lateral excursion was introduced as the loading condition in the analysis. As a consequence of the movement asymmetry, the discs were subjected to different load distributions. It was observed that the maximal shear stresses were located in the lateral zone of both discs and that the lateral attachment of the ipsilateral condyle-disc complex suffered a large distortion, due to the compression of this disc against the inferior surface of the temporal bone. These results may be related with possible consequences of a common disorder called bruxism. Although it would be necessary to perform an exhaustive analysis of this disorder, including the contact forces between the teeth during grinding, it could be suggested that a continuous lateral movement of the jaw may lead to perforations of both discs in their lateral part and may damage the lateral attachments of the disc to the condyle.     

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.     

Constraints and fixation for implanted joint replacements

Observations made on both cemented and uncemented joint prostheses after about two years of use show a layer of fibrous tissue next to the bone. This fibrous layer smooths over surface features up to l mm in size, and must be assumed to be weak in tension and shear. Prosthesis-bone interfaces should therefore be designed to transmit all forces acting on the prosthesis as compressive stresses. The forces acting across prostheses are controlled by the constraints exerted by the articulating surfaces; if the useful ligaments present are allowed to do their job of transmitting tension, the articulating surfaces can be designed to transmit only such forces and moments as can safely be transmitted across the prosthesis-bone interfaces.     

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 bone remodelling with piezoelectric effects using an open-source framework

Bone tissue exhibits piezoelectric properties and thus is capable of transforming mechanical stress into electrical potential. Piezoelectricity has been shown to play a vital role in bone adaptation and remodelling processes. Therefore, to better understand the interplay between mechanical and electrical stimulation during these processes, strain-adaptive bone remodelling models without and with considering the piezoelectric effect were simulated using the Python-based open-source software framework. To discretise numerical attributes, the finite element method (FEM) was used for the spatial variables and an explicit Euler scheme for the temporal derivatives. The predicted bone apparent density distributions were qualitatively and quantitatively evaluated against the radiographic scan of a human proximal femur and the bone apparent density calculated using a bone mineral density (BMD) calibration phantom, respectively. Additionally, the effect of the initial bone density on the resulting predicted density distribution was investigated globally and locally. The simulation results showed that the electrically stimulated bone surface enhanced bone deposition and these are in good agreement with previous findings from the literature. Moreover, mechanical stimuli due to daily physical activities could be supported by therapeutic electrical stimulation to reduce bone loss in case of physical impairment or osteoporosis. The bone remodelling algorithm implemented using an open-source software framework facilitates easy accessibility and reproducibility of finite element analysis made.

     

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 analysis of trabecular bone structure: a comparison of image-based meshing techniques

In this study, we investigate if finite element (FE) analyses of human trabecular bone architecture based on 168 μm images can provide relevant information about the bone mechanical characteristics. Three human trabecular bone samples, one taken from the femoral head, one from the iliac crest, and one from the lumbar spine, were imaged with micro-computed tomography (micro-CT) using a 28 μm resolution. After reconstruction the resolution was coarsened to 168 μm. First, all reconstructions were thresholded and directly converted to FE-models built of hexahedral elements. For the coarser resolutions of two samples, this resulted in a loss of trabecular connections and a subsequent loss of stiffness. To reduce this effect, a tetrahedral element meshing based on the marching cubes algorithm, as well as a modified hexahedron meshing, which thresholds the image such that load carrying bone mass is preserved, were employed. For each sample elastic moduli and tissue Von Mises stresses of the three different 168 μm models were compared to those from the hexahedron 28 μm model. For one sample the hexahedron meshing at 168 μm produced excellent results. For the other two samples the results obtained from the hexahedral models at 168 μm resolution were poor. Considerably better results were attained for these samples when using the mass-compensated or tetrahedron meshing techniques. We conclude that the accuracy of the FE-models at 168 μm strongly depends on the bone morphology, in particular its trabecular thickness. A substantial loss of trabecular connections during the hexahedron meshing process indicates that poor FE results will be obtained. In this case the tetrahedron or mass-compensated hexahedron meshing techniques can reduce the loss of connections and produce better results than the plain hexahedron meshing techniques.     

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 displacement analysis of a lung

There is considerable physiological evidence that the regional variation in pleural pressure and expansion of the lung is largely determined by gravity. In this paper a method is given based on the technique of finite elements which determines theoretically the mechanical behaviour of a lung-shaped body loaded by its own weight. The results of this theoretical analysis have been compared with actual measurements of alveolar size and pleural pressures in animal lungs.     

Finite element analysis of moment-rotation relationships for human cervical spine

A comprehensive, geometrically accurate, nonlinear C0-C7 FE model of head and cervical spine based on the actual geometry of a human cadaver specimen was developed. The motions of each cervical vertebral level under pure moment loading of 1.0 Nm applied incrementally on the skull to simulate the movements of the head and cervical spine under flexion, tension, axial rotation and lateral bending with the inferior surface of the C7 vertebral body fully constrained were analysed. The predicted range of motion (ROM) for each motion segment were computed and compared with published experimental data. The model predicted the nonlinear moment-rotation relationship of human cervical spine. Under the same loading magnitude, the model predicted the largest rotation in extension, followed by flexion and axial rotation, and least ROM in lateral bending. The upper cervical spines are more flexible than the lower cervical levels. The motions of the two uppermost motion segments account for half (or even higher) of the whole cervical spine motion under rotational loadings. The differences in the ROMs among the lower cervical spines (C3-C7) were relatively small. The FE predicted segmental motions effectively reflect the behavior of human cervical spine and were in agreement with the experimental data. The C0-C7 FE model offers potentials for biomedical and injury studies.     

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.