Finite Element Analysis of Interbody Cages in a Human Lumbar Spine

An innovative surgical procedure is vertebral stabilization by interbody cages. It is currently being used to separate and stabilize vertebral bodies and to promote bony fusion of the vertebrae onto or through the cages. This surgery, at some spine levels, can be performed through a laparoscope as an outpatient procedure with low morbidity. Because the procedure is new, little structural information is available on the interbody cages. The objective of this study was to evaluate the human lumbar spine stabilized by interbody cages biomechanically. The finite element method was used to compare cage designs by considering stresses in the cage and in the bone as well as relative displacements between the cage and the adjacent bone at the interface. The biomechanical evaluation considered different bone densities and considered axial, torsional, and bending loads on the lumbar spine. Stress analysis predicts local regions of stress concentration that could be damaging to cancellous bone and will likely require a remodeling response for local damage. This study predicts relative micromotion that could cause the bone resorption and fibrous tissue formation on the contact surfaces of the cage. The geometric constraints caused by the use of two cages will reduce the relative motion and therefore be more likely to allow bone ingrowth at the posterocentral contact region. Finite element analysis suggests that cages are a promising method for separation and stabilization of the vertebral bodies.     

Finite Element Analysis of Interbody Cages in a Human Lumbar Spine

Abstract

An innovative surgical procedure is vertebral stabilization by interbody cages. It is currently being used to separate and stabilize vertebral bodies and to promote bony fusion of the vertebrae onto or through the cages. This surgery, at some spine levels, can be performed through a laparoscope as an outpatient procedure with low morbidity. Because the procedure is new, little structural information is available on the interbody cages. The objective of this study was to evaluate the human lumbar spine stabilized by interbody cages biomechanically. The finite element method was used to compare cage designs by considering stresses in the cage and in the bone as well as relative displacements between the cage and the adjacent bone at the interface. The biomechanical evaluation considered different bone densities and considered axial, torsional, and bending loads on the lumbar spine. Stress analysis predicts local regions of stress concentration that could be damaging to cancellous bone and will likely require a remodeling response for local damage. This study predicts relative micromotion that could cause the bone resorption and fibrous tissue formation on the contact surfaces of the cage. The geometric constraints caused by the use of two cages will reduce the relative motion and therefore be more likely to allow bone ingrowth at the posterocentral contact region. Finite element analysis suggests that cages are a promising method for separation and stabilization of the vertebral bodies.     

Therapeutic electric stimulation does not affect immune status in healthy individuals – a preliminary report

Background

For the treatment of low back pain, the following three scenarios of posterior lumbar interbody fusion (PLIF) were usually used, i.e., PLIF procedure with autogenous iliac bone (PAIB model), PLIF with cages made of PEEK (PCP model) or titanium (Ti) (PCT model) materiel. But the benefits or adverse effects among the three surgical scenarios were still not fully understood.

Method

Finite element analysis (FEA), as an efficient tool for the analysis of lumbar diseases, was used to establish a three-dimensional nonlinear L1-S1 FE model (intact model) with the ligaments of solid elements. Then it was modified to simulate the three scenarios of PLIF. 10?Nm moments with 400?N preload were applied to the upper L1 vertebral body under the loading conditions of extension, flexion, lateral bending and torsion, respectively.

Results

Different mechanical parameters were calculated to evaluate the differences among the three surgical models. The lowest stresses on the bone grafts and the greatest stresses on endplate were found in the PCT model. The PCP model obtained considerable stresses on the bone grafts and less stresses on ligaments. But the changes of stresses on the adjacent discs and endplate were minimal in the PAIB model.

Conclusions

The PCT model was inferior to the other two models. Both the PCP and PAIB models had their own relative merits. The findings provide theoretical basis for the choice of a suitable surgical scenario for different patients.     

A comparison of the peri-implant bone stress generated by the preload with screw tightening between ‘bonded’ and ‘contact’ model

A number of finite element analyses (FEAs) for the dental implant were performed without regard for preload and with all interfaces ‘fixed-bonded’. The purpose of this study was comparing the stress distributions between the conventional FEA model with all contacting interfaces ‘fixed-bonded’ (bonded model) and the model with the interfaces of the components in ‘contact’ with friction simulated as a preloaded implant (contact model). We further verified the accuracy of the result of the FEA using model experiment. In the contact model, the stress was more widely distributed than in the bonded model. From the model study, the preload induced by screw tightening generated strain at the peri-implant bone, even before the application of external force. As a result, the bonded model could not reproduce the mechanical phenomena, whereas the contact model is considered to be appropriate for analysing mechanical problems.     

Possible role of stresses in inducing spinal stenosis--a long term complication following disk excision

A three-dimensional finite element model of an intact ligamentous lumbar motion segment (L3-4) was used to predict stresses in the pars interarticularis regions of the modeled vertebral bodies. The changes in stresses following disk excision, as compared to the intact model, also were computed. The predicted results show an increase in stresses in the posterior bony elements following disk excision. In some patients over a long period of time this increase in stresses, in association with other clinical factors, may lead to bony hypertrophy of the structures that surround the nerve roots. Ultimately, over a long period of time the increase in pressure on the entrapped nerve root may induce recurrent pain and other complications reported in the literature.     

A finite element study of the effect of diametral interface gaps on the contact areas and pressures in uncemented cylindrical femoral total hip components

Uncemented femoral total hip components rely entirely on contact with the prepared femur for their initial fixation. The contact areas and stresses between a straight tubular bone and a metal cylindrical prosthesis 12.5 cm long and 13 mm in diameter were calculated in a finite element model which includes uniform diametral gaps varying from 20 to 500 microns, using transverse loads from 100 to 2000 N. Frictionless three-dimensional contact elements were used between the bone and the prosthesis. Contact stresses were high and irregular in all cases, and the contact areas were small. Two regions of contact were apparent for lower loads and larger gaps. A third region of contact occurred near the distal tip of the implant at higher loads. This region of contact markedly increased the contact stresses at the distal tip of the prosthesis. A 20 microns overlap between bone and implant was modelled to assess a slight interference fit. The contact stress distribution in this case was markedly different from the stress distribution with a 20 microns diametral gap. The data collectively indicates that gaps of less than 20 microns between bone and implant can substantially change contact stress distributions.     

The effect of the interface on the bone stresses beneath tibial components

It was proposed that the stresses in the layer of bone immediately beneath a tibial component are an important determinant of fixation durability. Using finite element analysis, (ANSYS), the stresses were determined as a function of the amount of bone resection, the localization or completeness of implant-bone contact, and the interface material. The model was of two-dimensional sagittal slices consisting of quadrilateral elements (1 mm) with a range of seventeen material properties determined by CT scans. Typical prosthesis designs shifted the center of pressure more centrally rather than posteriorly, and thus increased anterior bone stresses. Resection up to 10 mm could actually decrease bone stresses due to an increase in bone surface area as long as complete coverage was obtained. A cement interface or direct metal on bone produced identical stresses. However a 1 mm complian: interface significantly reduced stresses in regions of high elastic modulus gradient. For rigid interfaces, the contact can be irregular, which leads to areas of over and understressing of bone. These conclusions have implications related to implant design.     

Hemiarthroplasty of hip joint: An experimental validation using porcine acetabulum

Biphasic properties of articular cartilage allow it to be an excellent bearing material and have been studied through several simplified experiments as well as finite element modelling. However, three-dimensional biphasic finite element (FE) models of the whole joint are rare. The current study was carried out to experimentally validate FE methodology for modelling hemiarthroplasty. Material properties such as equilibrium elastic modulus and permeability of porcine acetabular cartilage were initially derived by curve-fitting an experimental deformation curve with that obtained using FE. These properties were then used in the hemiarthroplasty hip joint modelling. Each porcine acetabular cup was loaded with 400N using a 34mm diameter CoCr femoral head. A specimen-specific FE model of each acetabular cup was created using μCT and a series of software processes. Each model was analysed under conditions similar to those tested experimentally. Contact stresses and contact areas predicted by the model, immediately after loading, were then compared with the corresponding experimentally measured values. Very high peak contact stresses (maximum experimental: 14.09MPa) were recorded. A maximum difference of 12.42% was found in peak contact stresses. The corresponding error for contact area was 20.69%. Due to a fairly good agreement in predicted and measured values of contact stresses and contact areas, the integrated methodology developed in this study can be used as a basis for future work. In addition, FE predicted total fluid load support was around 80% immediately after loading. This was lower than that observed in conforming contact problems involving biphasic cartilage and was due to a smaller local contact area and variable clearance making fluid exudation easier.     

Responses to the change in the environment in pairs of male rats genetically selected for activity level

Laboratory Wistar strain rats were genetically selected for high (+A) and low (-A) activity level. In thirteen pairs of adult males of the 23rd filial generation reactions to changes in the external environment were studied. The animals were housed in breeding cages four each. Two parallel studies were conducted: in pairs simultaneously placed into a novel environment (NOV), empty cages of the same dimensions as the home cage (HC), in the second, behaviour of the second pair that remained in the HC, after removal of two cage-mates, was tested. Once a minute, for a period of one hour, the type of activity was recorded and noted whether it was an element effected in contact with the partner or without any contact. The animals +A and -A differed in the frequency of various types of activity and immobility, in the ratio between behavioural manifestations shown in or without contact as well as in the response to the type of modified environment. To changes in the situation, whether removed cage-mates from the HC or placed into NOV +A animals reacted with a high wave of environment exploration which gradually habituated. -A rats equally responded with exploration but on a lower level. In +rats we recorded more frequently exploration without contact with the partner in HC and NOV in comparison with -A, more frequent grooming, less immobility in contact and with no contact. Between +A partners there was a greater number of contacts in NOV than in HC whereas in the -A group the incidence of contact did not differ between HC and NOV. ANOVA revealed the influence of factors of genetics and environment and interaction in several behavioural categories. The simple and in time economical method demonstrated the possibility of use for the detection of differences between +A and -A lines even at relatively small changes in the external stimulatory situation.     

Prediction of biomechanical parameters in the lumbar spine during static sagittal plane lifting.

A combined approach involving optimization and the finite element technique was used to predict biomechanical parameters in the lumbar spine during static lifting in the sagittal plane. Forces in muscle fascicles of the lumbar region were first predicted using an optimization-based force model including the entire lumbar spine. These muscle forces as well as the distributed upper body weight and the lifted load were then applied to a three-dimensional finite element model of the thoracolumbar spine and rib cage to predict deformation, the intradiskal pressure, strains, stresses, and load transfer paths in the spine. The predicted intradiskal pressures in the L3-4 disk at the most deviated from the in vivo measurements by 8.2 percent for the four lifting cases analyzed. The lumbosacral joint flexed, while the other lumbar joints extended for all of the four lifting cases studied (rotation of a joint is the relative rotation between its two vertebral bodies). High stresses were predicted in the posterolateral regions of the endplates and at the junctions of the pedicles and vertebral bodies. High interlaminar shear stresses were found in the posterolateral regions of the lumbar disks. While the facet joints of the upper two lumbar segments did not transmit any load, the facet joints of the lower two lumbar segments experienced significant loads. The ligaments of all lumbar motion segments except the lumbosacral junction provided only marginal moments. The limitations of the current model and possible improvements are discussed.     

Finite element estimates of interface stress in the trans-tibial prosthesis using gap elements are different from those using automated contact

When compared with automated contact methods of finite element (FE) analyses, gap elements have certain inherent disadvantages in simulating large slip of compliant materials on stiff surfaces. However, automated contact has found limited use in the biomechanical literature. A non-linear, three-dimensional, geometrically accurate, FE analysis of the trans-tibial limb-socket prosthetic system was used to compare an automated contact interface model with a gap element model, and to evaluate the sensitivity of automated contact to interfacial coefficient of friction (COF). Peak normal stresses and resultant shear stresses were higher in the gap element model than in the automated contact model, while the maximum axial slip was less. Under proximally directed load, compared with automated contact, gap elements predicted larger areas of stress concentration that were located more distally. Gap elements did not predict any relative slip at the distal end, and also transmitted a larger proportion of axial load as shear stress. Both models demonstrated non -linear sensitivity to COF, with larger variation at lower magnitudes of COF. By imposing physical connections between interface surfaces, gap elements distort the interface stress distributions under large slip. Automated contact methods offer an attractive alternative in applications such as prosthetic FE modeling, where the initial position of the limb in the socket is not known, where local geometric features have high design significance, and where large slip occurs under load.     

Stress analysis of compression plate fixation and its effects on long bone remodeling

A three-dimensional finite element model is generated for an intact plexiglass tube with an attached six-hole stainless steel compression plate. The results for a wide range of loads, including cyclic external loads and static tensile preloads in the plate and screws, are examined as specifically related to plate-induced osteopenia. The model demonstrates that disuse osteopenia, resulting from a reduction in magnitude of cyclic axial stress, should be limited to the central region between the inner screws. Also, the addition of a static preload negates any reduced axial stress levels in this region, thus raising questions on the relative importance of static and cyclic stresses for the internal remodeling of bone.     

How plate positioning impacts the biomechanics of the open wedge tibial osteotomy; A finite element analysis

A numerical model of the medial open wedge tibial osteotomy based on the finite element method was developed. Two plate positions were tested numerically. In a configuration, (a), the plate was fixed in a medial position and (b) in an anteromedial position. The simulation took into account soft tissues preload, muscular tonus and maximal gait load.

The maximal stresses observed in the four structural elements (bone, plate, wedge, screws) of an osteotomy with plate in medial position were substantially higher (1.13–2.8 times more) than those observed in osteotomy with an anteromedial plate configuration. An important increase (1.71 times more) of the relative micromotions between the wedge and the bone was also observed. In order to avoid formation of fibrous tissue at the bone wedge interface, the osteotomy should be loaded under 18.8% (～50 kg) of the normal gait load until the osteotomy interfaces union is achieved.     

How plate positioning impacts the biomechanics of the open wedge tibial osteotomy; a finite element analysis

A numerical model of the medial open wedge tibial osteotomy based on the finite element method was developed. Two plate positions were tested numerically. In a configuration, (a), the plate was fixed in a medial position and (b) in an anteromedial position. The simulation took into account soft tissues preload, muscular tonus and maximal gait load.The maximal stresses observed in the four structural elements (bone, plate, wedge, screws) of an osteotomy with plate in medial position were substantially higher (1.13-2.8 times more) than those observed in osteotomy with an anteromedial plate configuration. An important increase (1.71 times more) of the relative micromotions between the wedge and the bone was also observed. In order to avoid formation of fibrous tissue at the bone wedge interface, the osteotomy should be loaded under 18.8% (approximately 50 kg) of the normal gait load until the osteotomy interfaces union is achieved.     

Finite-element modelling of femoral shaft fracture fixation techniques post total hip arthroplasty.

The presence of a femoral prosthesis superior to a shaft fracture severely complicates fixation and treatment. This study uses two-dimensional, multithickness, plane stress finite-element models of a femur with prosthesis to investigate the stresses developed with the application of three popular fixation techniques: revision to a long stem prosthesis, lateral plating with a cortical bone allograft strut and cerclage wires, and custom plate application with proximal Parham band fixation with distal cortical screws (Ogden plate). The plate and bone contact as well as the fracture site contact were modelled by using orthotropic elements with custom-fit moduli so that only the normal stress to the interface was significant. A thermal analogy was used to model the cerclage and Parham band preloads so that representative preloads in the proximal fixation of the two types of plate treatments could be modelled. A parametric study was performed with the long-prosthesis model to show variations in stem lengths of one, two and three femoral diameters distal to the fracture site. The Ogden plate model showed a transfer of tensile stress near the proximal-most band, with the highest tensile stress being at the fracture site with evidence of stress shielding of the proximal lateral cortex. The cortical bone strut model showed a transfer of tensile stress to the bone strut but showed less shielding of the proximal cortex. The cerclage wires at the base of the bone strut showed the highest changes in load with the distalmost wire increasing to almost four times its original preload.(ABSTRACT TRUNCATED AT 250 WORDS)     

A finite element evaluation of the moment arm hypothesis for altered vertebral shear failure force

The mechanism of vertebral shear failure is likely a bending moment generated about the pars interarticularis by facet contact, and the moment arm length (MAL) between the centroid of facet contact and the location of pars interarticularis failure has been hypothesised to be an influential modulator of shear failure force. To quantitatively evaluate this hypothesis, anterior shear of C3 over C4 was simulated in a finite element model of the porcine C3–C4 vertebral joint with each combination of five compressive force magnitudes (0–60% of estimated compressive failure force) and three postures (flexed, neutral and extended). Bilateral locations of peak stress within C3's pars interarticularis were identified along with the centroids of contact force on the inferior facets. These measurements were used to calculate the MAL of facet contact force. Changes in MAL were also related to shear failure forces measured from similar in vitro tests. Flexed and extended vertebral postures respectively increased and decreased the MAL by 6.6% and 4.8%. The MAL decreased by only 2.6% from the smallest to the largest compressive force. Furthermore, altered MAL explained 70% of the variance in measured shear failure force from comparable in vitro testing with larger MALs being associated with lower shear failure forces. Our results confirmed that the MAL is indeed a significant modulator of vertebral shear failure force. Considering spine flexion is necessary when assessing low-back shear injury potential because of the association between altered facet articulation and lower vertebral shear failure tolerance.     

Dynamic characterization of human breast cancer cells using a piezoresistive microcantilever

In this paper, frequency response (dynamic compression and recovery) is suggested as a new physical marker to differentiate between breast cancer cells (MCF7) and normal cells (MCF10A). A single cell is placed on the laminated piezoelectric actuator and a piezoresistive microcantilever is placed on the upper surface of the cell at a specified preload displacement (or an equivalent force). The piezoelectric actuator excites the single cell in a sinusoidal fashion and its dynamic deformation is then evaluated from the displacement converted by measuring the voltage output through a piezoresistor in the microcantilever. The microcantilever has a flat contact surface with no sharp tip, making it possible to measure the overall properties of the cell rather than the local properties. These results indicate that the MCF7 cells are more deformable in quasi-static conditions compared with MCF10A cells, consistent with known characteristics. Under conditions of high frequency of over 50 Hz at a 1 μm preload displacement, 1 Hz at a 2 μm preload displacement, and all frequency ranges tested at a 3 μm preload displacement, MCF7 cells showed smaller deformation than MCF10A cells. MCF7 cells have higher absorption than MCF10A cells such that MCF7 cells appear to have higher deformability according to increasing frequency. Moreover, larger preload and higher frequencies are shown to enhance the differences in cell deformability between the MCF7 cells and MCF10A cells, which can be used as a physical marker for differentiating between MCF10A cells and MCF7 cells, even for high-speed screening devices.     

Effect of meniscus replacement fixation technique on restoration of knee contact mechanics and stability

The menisci are important biomechanical components of the knee. We developed and validated a finite element model of meniscal replacement to assess the effect of surgical fixation technique on contact behavior and knee stability. The geometry of femoral and tibial articular cartilage and menisci was segmented from magnetic resonance images of a normal cadaver knee using MIMICS (Materialise, Leuven, Belgium). A finite element mesh was generated using HyperWorks (Altair Inc, Santa Ana, CA). A finite element solver (Abaqus v6.9, Simulia, Providence, RI) was used to compute contact area and stresses under axial loading and to assess stability (reaction force generated during anteroposterior translation of the femur). The natural and surgical attachments of the meniscal horns and peripheral rim were simulated using springs. After total meniscectomy, femoral contact area decreased by 26% with a concomitant increase in average contact stresses (36%) and peak contact stresses (33%). Replacing the meniscus without suturing the horns did little to restore femoral contact area. Suturing the horns increased contact area and reduced peak contact stresses. Increasing suture stiffness correlated with increased meniscal contact stresses as a greater proportion of tibiofemoral load was transferred to the meniscus. A small incremental benefit was seen of simulated bone plug fixation over the suture construct with the highest stiffness (50 N/mm). Suturing the rim did little to change contact conditions. The nominal anteroposterior stiffness reduced by 3.1 N/mm after meniscectomy. In contrast to contact area and stress, stiffness of the horn fixation sutures had a smaller effect on anteroposterior stability. On the other hand suturing the rim of the meniscus affected anteroposterior stability to a much larger degree. This model emphasizes the importance of the meniscus in knee biomechanics. Appropriate meniscal replacement fixation techniques are likely to be critical to the clinical success of meniscal replacement. While contact conditions are mainly sensitive to meniscus horn fixation, the stability of the knee under anteroposterior shear loads appeared to be more sensitive to meniscal rim fixation. This model may also be useful in predicting the effect of biomaterial mechanical properties and meniscal replacement shape on knee contact conditions.     

Development of a finite element musculoskeletal model with the ability to predict contractions of three-dimensional muscles

Representation of realistic muscle geometries is needed for systematic biomechanical simulation of musculoskeletal systems. Most of the previous musculoskeletal models are based on multibody dynamics simulation with muscles simplified as one-dimensional (1D) line-segments without accounting for the large muscle attachment areas, spatial fibre alignment within muscles and contact and wrapping between muscles and surrounding tissues. In previous musculoskeletal models with three-dimensional (3D) muscles, contractions of muscles were among the inputs rather than calculated, which hampers the predictive capability of these models. To address these issues, a finite element musculoskeletal model with the ability to predict contractions of 3D muscles was developed. Muscles with realistic 3D geometry, spatial muscle fibre alignment and muscle-muscle and muscle-bone interactions were accounted for. Active contractile stresses of the 3D muscles were determined through an efficient optimization approach based on the measured kinematics of the lower extremity and ground force during gait. This model also provided stresses and strains of muscles and contact mechanics of the muscle-muscle and muscle-bone interactions. The total contact force of the knee predicted by the model corresponded well to the in vivo measurement. Contact and wrapping between muscles and surrounding tissues were evident, demonstrating the need to consider 3D contact models of muscles. This modelling framework serves as the methodological basis for developing musculoskeletal modelling systems in finite element method incorporating 3D deformable contact models of muscles, joints, ligaments and bones.     

The effect of direct and indirect force transmission on peri-implant bone stress – a contact finite element analysis

In almost all finite element (FE) studies in dentistry, virtual forces are applied directly to dentures. The purpose of this study was to develop a FE model with non-linear contact simulation using an antagonist as force transmitter and to compare this with a similar model that uses direct force transmission. Furthermore, five contact situations were created in order to examine their influence on the peri-implant bone stresses, which are relevant to the survival rate of implants. It was found that the peri-implant bone stresses were strongly influenced by the kind of force transmission and contact number.