Computational Modeling of the Spine

Published on 02/04/2015 by admin

Filed under Neurosurgery

Last modified 22/04/2025

Print this page

rate 1 star rate 2 star rate 3 star rate 4 star rate 5 star
Your rating: none, Average: 3 (1 votes)

This article have been viewed 2630 times

Chapter 22 Computational Modeling of the Spine

An anatomic experimental approach has been the gold standard for the biomechanical evaluation of spinal structures. In vitro biomechanical investigations commonly use cadavers, animal carcasses, or synthetic models. Among these options, cadaveric investigations produce the most reliable and realistic results. However, despite its high clinical relevance, fresh-frozen cadaveric spinal specimens may be very expensive (>$1000 in the United States) or difficult to acquire. Animal models such as calves, sheep, and pigs, on the other hand, are relatively inexpensive (~$100). Animal models offer two main advantages. First, bovine, ovine, and porcine spine models exhibit gross anatomic structures (e.g., facets, processes, ligaments) similar to those of the human spine, which permits the convenient application of the instrumentation used in surgery. Secondly, between-subject variability is minimal with animal models, which reduces the statistical sample size demand. Synthetic models usually involve polyurethane foam surrogates, which are available in any density and geometry, ranging from blocks to any anatomic shape. They are quite inexpensive (<$100). In addition, synthetic models can be produced homogeneously, presenting little or no interspecimen or intraspecimen variation. Although synthetic models can represent the human anatomy better than animal models, and have less specimen-related variability, they are considered the least clinically relevant models because they cannot match the material properties of the natural bone.

Some investigators, facing the shortcomings of all of these approaches, turn to computational techniques. With the recent advances in computer and imaging technologies, finite element models of spinal segments or vertebrae can be developed and validated rapidly and inexpensively. Finite element analysis, a mathematical method performed by fictively dividing structures into simpler substructures such as triangles or cubes, can determine the stress or strain response of complex structures under imposed conditions by performing calculations on the simpler substructures. Once the anatomic model is introduced into the finite element analysis software, kinetic and kinematic analyses, as well as thermal, flow, and time-related problems, can be solved.

Computational techniques, especially finite element analysis, have provided substantial help for researchers in the area of musculoskeletal biomechanics for decades. The reduced costs and time (in most cases) and the ability to simulate sophisticated loading conditions and vary certain parameters (such as treatment strategies) in a perfectly controlled analysis environment are the major advantages of finite element analysis. Finite element modeling will become much more pervasive with developments in the hardware and software technology that provide the ability to build complex and comprehensive models that truly reconstruct real-life structures and conditions of the spine.

In this chapter, common finite element modeling procedures (i.e., obtaining three-dimensional [3D] models of spinal segments, determining structural and material properties, and assigning contact, loads, and boundary conditions) and some biomechanical applications are discussed.

Finite Element Models

Even though the finite element method has been applied to structural engineering problems since the 1950s, it was first applied to the biomechanics of the spine in the 1970s.1,2 With improvements in numerical techniques and computer technology, its use in the field has gained momentum and it has become an integral part of spine biomechanics that is viewed as a complementary and insightful analysis (Fig. 22-1).

Modeling of Spinal Elements

Finite element modeling of spine segments or vertebrae, similar to other engineering solutions, involves the transfer of the physical sample into the 3D reconstructed form in the computer environment, assignment of structural and material properties, definition of contact properties, and application of loads and boundary conditions. Although the finite element modeling procedure is essentially the same in most of the studies, many different techniques and properties are used in these steps.

Three-Dimensional Reconstruction and Structural Modeling

In the finite element method, geometric accuracy of the model has a significant effect on problem solving and can cause marked differences in the results.

Several methods have been commonly used to model spinal units, including touch probe digitizer,3 laser scanner,4 CT,5 and MRI.6 Even though precise, highly accurate geometric models can be obtained using the touch probe digitizer and laser scanner methods, noninvasive techniques such as CT and MRI often are preferred because of their ability to retrieve data on internal features of the specimens—and even from live subjects—with ease. Some techniques also combine some of the methods mentioned earlier for the reconstruction of spinal units.7

In thorough finite element models, vertebral bodies usually are modeled in several sections, as cortical shell, cancellous core, dorsal bony elements, bony end plates, and cartilage layers, due to their unique material properties.5,7,8

In almost all studies, intervertebral discs are formed by three significant parts: the nucleus pulposus, the anulus fibrosus, and the cartilaginous end plates. The anulus is further divided into subcomponents such as ground substance and anulus fibers. Anulus ground substance is modeled as several layers surrounding the nucleus. The anulus fibers are angled 30 degrees and 150 degrees with respect to the transverse plane in adjacent layers.9,10 Some studies have used more realistic approaches in the design of these fibers by considering the increase in fiber angle from ventral to dorsal, and also fiber angle variation between anulus layers5,11 motivated by the histologic findings.12,13

After a solid model of the spine is built, each of its components is carefully divided into simpler substructures called “elements” for mathematical calculation purposes. Generally, hard tissues such as cortical shell, trabecular core, and dorsal elements are composed of tetrahedral (pyramid) and pentahedron (wedge) elements.14,15 With soft tissues, such as the anulus ground substance or cartilage layers, it is necessary to use hexahedral (rectangular prism) elements due to low-elasticity modulus and high Poisson ratio properties. In most studies, bilinear link or spring elements (cables), which are only capable of sustaining tension, are employed for ligaments and anulus fibers.5, 911

Material Definition

Material definition is one of the most sensitive and intrinsic parts of finite element modeling of the spine. Due to variability within and among test specimens, different researchers have defined many material properties for the different components the spine.

Vertebral Bodies, End Plates, and Facet Cartilages

Even though bone is a poroelastic, anisotropic, and viscoelastic structure, in reality,16 vertebral cortical shell and cancellous bone usually are considered as linear elastic isotropic or transversely isotropic. Dorsal bony elements most commonly are modeled as linear elastic and isotropic (Tables 22-1 to 22-3). Researchers have either neglected the end plate due to its low elasticity module or modeled it as a composite structure (i.e., bony and cartilaginous layers). End plates most commonly were modeled as linear elastic and isotropic bone with a cartilaginous layer (Table 22-4). Cartilage layers also are added to the facets in most modeling paradigms (Table 22-5).

Articular Facet Joints

The facet joint plays an important role in load bearing and restriction of motion in segmental movements. It functions in tandem with the intervertebral disc in load transfer between adjacent vertebral bodies.17,18 The articulating structure must be modeled with realistic attributes and proper procedures in computational modeling of spinal segments because it has a significant effect on the quality and quantity of the motion. The articulation characteristics and relative motion depend on many factors, including cartilage layers, gap distance, the condition of the intervertebral disc, loading type, and geometric features of the articulating surface.

Kumaresan et al.19 compared different techniques for modeling of cervical facet joints. They reconstructed the joint capsule using four different methods: slide-line, contact plane, hyperelastic, and fluid models. The slide-line and contact plane models lacked the synovial fluid and synovial membrane. In the slide-line model, interaction was defined by using slide-line elements (gap elements), whereas in the contact plane model it was achieved by defining a contact plane between the cartilage surfaces. A friction coefficient of 0.01 was assigned for both models. The other models, hyperelastic and fluid, did not have elements for contact modeling; instead, they included synovial fluid between the cartilage layers modeled using noncompressible hyperelastic solid elements and hydrostatic noncompressible fluid elements, respectively. They concluded that fluid modeling of the facet joint matched both actual facet joint anatomy and its function better than the other three models. Similarly, Wheeldon et al.20 used full anatomic parts, such as facet joint cartilage, facet joint synovial fluid, and facet joint membrane, for modeling of the articular facet joints.

Shirazi-Adl et al.17,18 viewed the articulation process as a general moving contact problem and assigned 1 mm for the initial gap between the articulating surfaces. Guan et al.7,21 used 3D surface-to-surface contact to simulate the articulation phenomenon and assigned the initial gap between the surfaces according to CT images. Three-dimensional eight-noded compression-only gap elements were used to simulate the articulation between the facets in a study by Goel et al.22 Based on the CT images, they assigned the value of 0.45 mm for the initial gap between the surfaces.

The facet joints also are modeled with 3D gap contact elements with cartilage layers modeled using a parameter of “softened contact,” which exponentially adjusts force transfer between facet surfaces, according to the size of the gap.9,2327 An initial gap of 0.5 mm was assigned, and at full closure, the joint was assumed to have the same stiffness as the surrounding bone material.

Schmidt et al.28,29 modeled the facet joints as nonlinear, with frictionless contact. They assigned the value of 0.6 mm for the initial gap between the cartilage layers. In other studies,30,31 the facet joints were simulated using frictionless surface-to-surface contact elements in combination with the penalty algorithm with a normal contact stiffness of 200 N/mm, and the initial gap between the cartilage layers assumed to be 0.4 mm.

Some researchers viewed the facet joints as a 3D contact problem with friction.8,3234 To allow random motion of the surfaces, such as separation, sliding, and rotation, they defined a finite sliding interaction. For the friction characteristic, they chose a classic isotropic Coulomb friction model and assigned a relatively high friction coefficient of 0.1 as a worst-case scenario.

Lu et al.35 modeled the facet joints using sliding contact elements and assumed the contact pressure to be between 0 and 5000 megapascals (MPa), depending on the gap size. Pressure has been considered to be 0 MPa at a gap size of 0.5 mm and 5000 MPa when the gap is closed completely.

Loads and Boundary Conditions

Application of loads and boundary conditions varies according to the aim of the study. Generally, loading types used in the computational analysis of the spine are simple static or incrementally altering quasistatic loads. In general, the inferior end plate of the lowest vertebra is fixed rigidly, and loads (or displacements) are applied to the superior end plate of the uppermost vertebral body.

In most cases it is easy to simulate tension or compression loads. It is more complicated to apply moments, add muscle forces and upper body weight effects, or realize coupled motions.

Shirazi-adl et al.17,18 and Polikeit et al.8,34 applied a linearly varying distribution of axial loads to the top of the upper vertebra to simulate pure bending moments, with no net axial force. Kumaresan et al.36 and Wheeldon et al.20 succeeded in producing pure moments by using a force couple acting on the rigid plate attached to the superior vertebra.

Sharma et al.37 simulated pure flexion-extension moments using a pair of concentrated axial loads on top of the upper vertebra. They also studied the effect of anteroposterior shear during sagittal motion by a horizontal force applied at the center of the upper vertebral body. Two equal and opposite forces at the nodes at the edge of the sagittal or medial periphery of the upper vertebral body were applied to produce a force couple yielding a pure moment.3842

Researchers have used the follower load technique, suggested by Patwardhan et al.,43 extensively for the last decade to simulate upper-body weight and the stabilizing effects of muscle forces seen in vivo without inducing any intervertebral rotations.9,2528,30,44,45

Biomechanical Applications

Perhaps the greatest advantage of the finite element method for researchers in the area of musculoskeletal biomechanics is its ability to predict the consequences of surgical applications; age-related changes such as disc degeneration, osteophytes, and osteoporosis; or effects of daily activities on a spinal segment. It also is a great tool for the development and analysis of fusion instrumentation or motion preservation devices.

Disc Degeneration and Vertebral Osteoporosis Simulation

Disc degeneration is one of the most widely studied questions using finite element analysis of the spine because of the complex natural history of the disease and its high impact on clinical practice. Since the 1980s, a large body of literature on the modeling of the normal and degenerated disc has been published. The difficulty in proper simulation of the disc and disc degeneration arises from the hierarchical structural organization of the disc, the coexistence of multiple physical states (i.e., void, solid, fluid) presenting a nonhomogeneous material property distribution, and the influence of the degeneration cascade on these variables at various levels. The “weakest link” in modeling of the healthy disc and disc degeneration is validation of the model, because a significant challenge is presented by experimental evaluation of the material properties of the disc components.

In one of the earliest studies, Shirazi-Adl et al.45 modeled the degenerated disc by simply omitting a nucleus section. In another study, Shirazi-Adl et al.18 modeled the healthy disc as noncompressible fluid and suggested the loss of noncompressibility as a means to simulate the degenerated disc.

Goel et al.10 modeled the aging disc with a great deal of complexity by realistically simulating the radial and circumferential fissures in the anulus fibrosus that appear with age and discectomy. They simulated the radial crack in two ways: “out-in” injury started from the outermost layer and went through the innermost anulus layer with step-wise removal of the anulus elements, whereas “in-out” injury started from the innermost layer and went through the last layer of anulus simulated by replacing material of the anulus’s elements with nucleus material. They simulated the circumferential injury by removing the fourth anulus layer along the circumference, starting from posterior and extending through lateral.

Rohlmann et al.24 investigated the mechanical effects of three different grades of intervertebral disc degeneration—mild, moderate, and severe—compared with the healthy disc. They modeled the degeneration by reducing the disc height to 20%, 40%, and 60% of the healthy disc for mildly, moderately, and severely degenerated discs, respectively. In the finite element model, the laxity in all disc fibers and surrounding ligaments, seen on in vivo conditions after reduction of the disc height, was simulated by offsetting their nonlinear stiffness curves. They assigned the compressibility values for the nucleus pulposus of 0.0005, 0.0503, 0.0995, and 0.15, respectively, for healthy disc and mildly, moderately, and severely degenerated intervertebral discs. They assumed no effect on the material properties of the anulus fibrosus from degeneration of the disc.

Hussain et al.41,42 modeled two different degenerated discs, moderately and severely degenerated, to study the effects of disc degeneration on the facet loads. They simulated the degeneration according to the signs of Thompson disc grades 2 and 3 (moderate) and 4 and 5 (severe). They decreased the disc height, increased the disc area, and reduced the nucleus portion to simulate the degenerated disc. They also reduced the Poisson’s ratio in severely degenerated discs and increased the Young’s modulus of anulus and nucleus in both moderately and severely degenerated discs.

Galbusera et al.11 meticulously modeled the degeneration of the disc based on five different degenerative characteristics: (1) loss of water content in nucleus and anulus, (2) calcification and thickness reduction in end plates, (3) loss of disc height, (4) formation of osteophytes, and (5) diffuse sclerosis. They obtained ten different spinal segments by scaling the size of the finite element model to represent different stages of each degeneration level (mild, moderate, and severe). They changed the initial void ratio to model the degeneration in the nucleus and anulus, and they formed three different types of anular lesions: a rim lesion, a radial tear, and a circumferential tear. They modeled end plate cartilage degeneration by reducing its height and changing its material properties. Degeneration of the disc also was modeled with reduction of its height. When simulating the osteophytes, they used solid elements with the same material attributes of bone and anulus between them. They simulated sclerosis by interchanging the material properties of cancellous bone with the cortical bone in certain areas. They divided the lower half of L4 and the upper half of L5 into four regions and randomly assigned the sclerosis deficiency to these sections.

Polikeit et al.46 investigated the effects of osteoporosis and disc generation on the load transfer characteristics of an L2-3 segment. Initially, they modeled an osteoporotic vertebra by decreasing the elasticity modulus of trabecular bone by 66% and other bony parts by 33% compared with healthy model. Next, they simulated the progress of osteoporosis by modeling the cancellous bone transversely isotropic (i.e., the elastic modulus in the horizontal direction was one third of that in the vertical direction, using the same material properties as for the isotropic case) and increasing the anisotropy gradually. They postulated five steps of degeneration based on the isotropic osteoporotic vertebra model. At the first step, they changed the material definition of the nucleus pulposus from noncompressible to deformable. In the second step, the nucleus pulposus was hardened by assigning the same material properties as the ground substance of the anulus, and subsequently the Young modulus of the whole disc was doubled. In the fourth step, the innermost anulus layer was removed, other layers were thickened, and increased space between the fibers resulted. In the last step, they further reduced the elasticity modulus of the last remaining fiber layers. They also modeled a worst-case scenario by using transversely isotropic osteoporotic vertebral bodies in conjunction with the degenerated intervertebral disc.

Modeling and Analysis of the Instrumented Spine

A strong statistical power can be obtained in experimental approaches to spinal biomechanics by using a paired design analysis in which the same specimen is used in all experimental and treatment groups. This makes it possible to limit the variables and provide controlled analysis conditions. However, paired analysis design may not always be an option, especially in destructive studies, and the measurement of some parameters may not be possible without introducing confounding factors. Finite element analysis is an effective tool for observing changes in outcome parameters in response to varying input while controlling all other factors.

Analysis of Fusion

Guan et al.21 studied the range of motion, disc pressure, and facet joint loading following ventral fusion at L4-5 or L5-S1 under flexion-extension loading. They simulated fusion by replacing the whole disc with either a bone graft or a porous tantalum. They assumed an elasticity modulus of 100 MPa and 3.3 gigapascals (GPa), respectively, and a Poisson’s ratio of 0.2 and 0.1, respectively, for bone graft and porous tantalum. They found that caudal fusion, compared to rostral, and tantalum, compared to bone graft, caused less range of motion, and that rostral fusion induced more facet joint loads and less adjacent disc stress than caudal fusion.

Park et al.33 investigated the biomechanical effects of the number of ventral fixation rods, with or without dorsal fixation, on the spinal stability of thoracolumbar burst fracture constructs. They prepared ten separate models, composed of combinations of the following parameters: ventral rod number (0 to 2), with or without midcolumn compression and with or without dorsal fixation. They used different corpectomy procedures for the compression. They placed a cage in the center of L1. To simulate the surgical effects of the posterior fixation, the supraspinal ligament, interspinal ligament, and ligamentum flavum were removed. Pedicle screws and a ventral screw also were placed to simulate dorsal and ventral fixation. They concluded that two-rod ventral fixation ensured enough spinal stability and one-rod ventral fixation with posterior fixation hindered excessive motion.

Kim et al.47 focused on the biomechanical alterations after removal of the dorsal stabilization device in a fused lumbar segment. They modeled the fusion mass between the transverse processes of L3 and L4 as bilateral rectangular columns. To model the decompressed condition, the supraspinous ligament, the interspinous ligament, a portion of the spinous process, the distal region of the L3 lamina, and the ligamentum flavum were removed. Fixation was achieved with posterior pedicle screws. They inserted the screws into the pedicles at an inward angle of 10 degrees with respect to the horizontal plane and used a rod 5.5 mm in diameter. After intact testing, the fusion state with or without posterior support was studied. They concluded that removal of pedicle screws could lower the stiffness value of the fused segment, which would, however, decrease the stress levels at the adjacent discs.

Polikeit et al.8 investigated the effects of intervertebral cage applications on the load transfer and stress levels of adjacent tissues. They used a box-shaped cage model with titanium and carbon fiber for the material properties. They simulated the surgical effects by removing the anterior longitudinal ligament, nucleus pulposus, and ventral portion of the anulus. The gap elements with a high friction coefficient of 0.8 were used to define contact between cage surfaces and end plates. They found that intervertebral cage application increased the stress and changed the load transfer at the adjacent segments.

In an optimization study, Zhong et al.32 sought to devise a new intervertebral cage design that would allow more bone graft volume. They used a dorsal approach to implant the titanium cages and then carried out laminectomy and partial discectomy procedures on the model. They removed the dorsal elements, supraspinous ligament, interspinous ligament, ligamentum flavum, and part of the disc to simulate the procedure. The posterior support was modeled with pedicle screws and rods. The cages were placed between the vertebral bodies within two-thirds of the area of the disc. Cages could bear only compression loads. In the optimization step, they prioritized the gain of volume for bone ingrowth without sacrificing the biomechanical characteristics. They obtained a new cage with a volume of 1603 mm3 in comparison with the initial cage model with a volume of 2058 mm3. The new design was satisfying in terms of range of motion and stress levels at the adjacent disc, but had higher stress levels than the original cages.

Analysis of Arthroplasty

Goel et al.23 used a hybrid approach, a combination of finite element method and biomechanical testing, to investigate the mechanical effects of the Charité artificial disc (DePuy Spine, Inc., Raynham, MA) on the adjacent segments. They assumed that the surfaces of the artificial disc were smooth and flat. Furthermore, they also simulated the surgical procedure by resecting the anterior longitudinal ligament and ventral anulus and excising the nucleus pulposus. They concluded that arthroplasty increases mobility in flexion and extension by 19% and 44%, respectively, with minimal adjacent disc pressure changes. They also found that shear stress levels were higher at the superior bone interface with the artificial disc than at the inferior interface.

Rohlmann et al.48 investigated the biomechanical effects of artificial disc implant position and height, as well as gradual removal of the disc and contribution of the anterior longitudinal ligament to the biomechanical behavior of the lumbar spine. They used three different disc heights—10, 12, and 14 mm—placed in various locations in the anterior-posterior plane. Removal of the disc was modeled in steps: the ventral anulus, the nucleus, and, finally, excision of the whole disc. To study the contribution of the anterior longitudinal ligament, two models were used, with and without the ligament, simulating a ligament repair after implantation. They concluded that implant position had a great effect on intersegmental rotation for standing and flexion. The increase in disc height increased intersegmental rotation, and when lateral portions of anulus were preserved, intersegmental rotations were found similar to the intact model.

Schmidt et al.30 studied two hypotheses: (1) that increasing the artificial disc implanted segment number increased the mobility of the spine and facet joint loads, and (2) that incorrect positioning of the implanted disc caused instability. They used the SB Charité 3 Disc in the finite element model. They used chrome-cobalt alloy for the metallic end plates, with an elasticity modulus of 300 GPa and a Poisson’s ratio of 0.27. They chose ultra-high-molecular-weight polyethylene for the core with a Young’s modulus of 2 GPa and a Poisson ratio of 0.3. The anterior longitudinal ligament, ventral anulus, and nucleus were removed. They obtained a perfect osteointegration between the metallic end plates and the vertebral bodies by using a perfect bond contact. Eleven cases were studied, involving combinations of two, three, or four levels of implantation. The effect of misalignment of the discs also was studied. The results showed that the increase in motion ranged from 51% for two implants up to 91% for four implants. The increase in facet joint loads was 24% for two implants and up to 38% for four implants for the centrally positioned, extension situation. They found that misalignments in the position of the artificial disc induced high facet joint loads, unfavorable kinematics, and lift-off phenomena.

Analysis of Dynamic Systems

Zander et al.49 studied a hybrid construct with a dynamic implant placed adjacent to a fused segment with a spinal fixator. A pedicle screw–based dynamic implant was placed at the L3-4 segment. An L2-3 fusion initially was modeled with a bone graft and a rigid dorsal fixator. The location of the fusion was tested by moving the fusion and fixators to the L4-5 level. They concluded that the dynamic implant decreased intersegmental rotation of the implanted level in walking, flexion, and extension, and reduced facet joint loads in axial rotation.

Rohlman et al.9 compared the mechanical performance of dorsal dynamic implants and dorsal rigid fixators. They inserted a pair of dorsal dynamic implants with fictional material attributes to the L3-4 segment. To compare the dynamic and rigid devices, the dynamic system was substituted with the rigid system. They also divided the dynamic implant investigation by studying the construct both with and without distraction to the height of a healthy disc. Furthermore, they studied the effect of implant stiffness on intersegmental rotation by changing the longitudinal rod stiffness in several steps. They found that a dynamic implant decreased the intersegmental rotation and facet joint loads at the implanted level, and that it reduced intradiscal pressure in a healthy disc in extension and standing. They also found that intradiscal pressures were reduced only for the rigid fixator after distraction.

In a recent probabilistic study, Rohlmann et al.27 investigated the effects of ten variables inherent in a pedicle screw–based motion preservation system on the intradiscal pressure, facet loads, and intervertebral rotations. They used different lengths of implants in the study with titanium screws and polycarbonate-urethane rods. The parameters that varied were elastic modulus, diameter of the rods, angular rigidity of the screw head, interval between the rod and the pedicle screw entry point, orientation of the screws, whether cranial or caudal levels were bridged, and several kinds of defects. They observed that motion preservation systems reduced levels of intersegmental rotation, intradiscal pressure, and facet joint loads in most cases. They concluded that the loading conditions, Young’s modulus and the diameter of the rod, the distraction of the segment, and the torsional stiffness of the connection between the screws and rod had a greater effect on the output parameters.

In a parametric study, Schmidt et al.29 sought the optimal axial and bending stiffness values for a dorsal dynamic implant while validating their findings with an in vitro study. They used a simplified model of the stabilization system, including rod and screws, to investigate the range of axial and bending stiffness of the rod, which represented the nonfixed and fixed situations. They found that lower stiffness values of the rod provided acceptable stiffening to the segment and higher stiffness values of the rod had negligible effects.

Analysis of Cement Augmentation

Polikeit et al.34 investigated the influence of cement augmentation on load transfer by using an L2-3 finite element model. After adapting the material properties of the osteoporotic vertebral body, they considered polymethylmethacrylate (PMMA) as the bone cement material. In the first case, they simulated the effects of cement augmentation of the whole trabecular bone by changing the cancellous bone material properties to that of PMMA in L2 and L3, subsequently. They modeled unipedicular and bipedicular cement augmentations with one or two vertically placed cement barrels filling up one sixth of the vertebral body volume according to observations on radiographs of treated patients. The results of the study showed that augmentation increased the intradiscal pressure and deflection of the end plate into the adjacent vertebra. They found that augmentation changed the stress levels and distribution in the adjacent vertebrae.

Baroud et al.50 investigated the effects of vertebroplasty on load transfer in the segment. They modeled the material properties of bone cement composition after vertebroplasty as linear isotropic elastic, 46 times stronger and 12 times stiffer than the osteoporotic bone. To simulate the cement-infiltrated vertebra, they modeled the whole volume of L5 as augmented cement. They found that augmentation reduced bulging of the end plate of the augmented vertebra, which caused stiffening of the intervertebral joint.

Rohmann et al.51 studied vertebroplasty and kyphoplasty for compression fractures and to investigate the importance of the center of gravity shift compensation of the upper torso on muscle forces and intradiscal pressure. They prepared an intact model, a vertebroplasty model, and a kyphoplasty model. They simulated the vertebroplasty model by reducing the ventral height of L3 by 35% compared with an intact disc, forming a wedge fracture. They simulated the kyphoplasty model by reducing the ventral height of L3 by 90% of the intact disc, assuming a nearly full fracture. In addition, they simulated the upper torso center of gravity shift compensation with three different scenarios. They concluded that the shift of the upper body center of gravity played a greater role in vertebral body fractures at the adjacent segments compared with kyphoplasty or vertebroplasty.

Summary

Although experimental techniques traditionally have commanded a higher position in the hierarchical order of scientific approaches to investigating real-life questions in spine biomechanics, it is not always possible to use them, due to limitations of technology or specimen morphology. Finite element analysis has become an invaluable and relatively inexpensive tool for researchers to investigate questions that are either costly to answer using an experimental model or require measurements that cannot be performed without destroying the specimens. Computational studies, therefore, have contributed significantly to the literature to help in the understanding of complicated biomechanical principles such as facet loads, the relationship between osteoporosis and disc degeneration, as well as the effects on intervertebral load transfer, and implant design.

Despite the advantages of finite element analysis, one must be aware that information obtained using this approach does have limitations. Finite element analysis is founded on four major pillars: modeling, materials, interactions, and boundary conditions. Thanks to the progress in imaging technology, the modeling capabilities have been improved significantly, which has made possible the modeling of the details of bone, even at microscopic levels. The greatest limiting factor in finite element analysis is the material property definitions. Material properties of the elements of the spine (e.g., bone, ligaments, cartilage) are approximated from the previously existing information in the literature. This is associated with a totally separate set of limitations. Interactions between the facet surfaces or the implant-bone interface are based on assumptions such as gaps, friction, or fused surfaces. Loading usually is defined as pure moment, with or without the presence of axial compression. Therefore, the quality of the information obtained in finite element analysis depends on how realistic these assumptions were. The results of the finite element analysis should be treated as predictions rather than absolute facts. If the hypotheses were reasonably appropriate to the question studied, one can gain the most from the results of the computational analyses by observing such characteristics as relative changes in study groups, characteristics of the parameter behaviors during the course of loading, or qualitative analysis of stress or strain distributions.

Nevertheless, computational approaches will continue to complement the experimental science of spine biomechanics and become an increasingly useful tool for the study of complex problems, including the design of implants. Educational environments that provide surgeons with, at minimum, the foundations of finite element analysis are critical. There is no doubt that patient-specific finite element studies can help surgeons better prepare for surgical procedures. The field will most certainly advance by the development of faster modeling capabilities and techniques that can determine patient-specific material definitions.

References

1. Yoganandan N., Kumaresan S., Voo L., Pintar F.A. Finite element applications in human cervical spine modeling. Spine (Phila Pa 1976). 1996;21(15):1824-1834.

2. Moroney S.P., Schultz A.B., Miller J.A., Andersson G.B. Load-displacement properties of lower cervical spine motion segments. J Biomech. 1988;21(9):769-779.

3. Lee K.K., Teo E.C., Qiu T.X. Finite element modeling of L2-L3 using digitizer. Int J Comput Applications Technol, special issue on Biomedical Engineering and IT. 2003;20:1-9.

4. Heuer F., Schmidt H., Claes L., Wilke H.J. A new laser scanning technique for imaging intervertebral disc displacement and its application to modeling nucleotomy. Clin Biomech (Bristol, Avon). 2008;23(3):260-269.

5. Schmidt H., Heuer F., Drumm J., et al. Application of a calibration method provides more realistic results for a finite element model of a lumbar spinal segment. Clin Biomech (Bristol, Avon). 2007;22(4):377-384.

6. Pfirrmann C.W., Metzdorf A., Zanetti M., et al. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine (Phila Pa 1976). 2001;26(17):1873-1878.

7. Guan Y., Yoganandan N., Zhang J., et al. Validation of a clinical finite element model of the human lumbosacral spine. Med Biol Eng Comput. 2006;44(8):633-641.

8. Polikeit A., Ferguson S.J., Nolte L.P., Orr T.E. Factors influencing stresses in the lumbar spine after the insertion of intervertebral cages: finite element analysis. Eur Spine J. 2003;12(4):413-420.

9. Rohlmann A., Burra N.K., Zander T., Bergmann G. Comparison of the effects of bilateral posterior dynamic and rigid fixation devices on the loads in the lumbar spine: a finite element analysis. Eur Spine J. 2007;16(8):1223-1231.

10. Goel V.K., Monroe B.T., Gilbertson L.G., Brinckmann P. Interlaminar shear stresses and laminae separation in a disc. Finite element analysis of the L3-L4 motion segment subjected to axial compressive loads. Spine (Phila Pa 1976). 1995;20(6):689-698.

11. Galbusera F., Schmidt H., Neidlinger-Wilke C., et al. The mechanical response of the lumbar spine to different combinations of disc degenerative changes investigated using randomized poroelastic finite element models. Eur Spine J. 2011;20(4):563-571.

12. Eberlein R., Holzapfel G.A., Schulze-Bauer CAJ. An anisotropic constitutive model for annulus tissue and enhanced finite element analysis of intact lumbar disc bodies. Comput Methods Biomech Biomed Engin. 2001;4:209-230.

13. Cassidy J.J., Hiltner A., Baer E. Hierarchical structure of the intervertebral disc. Connect Tissue Res. 1989;23(1):75-88.

14. Rundell S.A., Guerin H.L., Auerbach J.D., Kurtz S.M. Effect of nucleus replacement device properties on lumbar spine mechanics. Spine (Phila Pa 1976). 2009;34(19):2022-2032.

15. Charriere E., Sirey F., Zysset P.K. A finite element model of the L5-S1 functional spinal unit: development and comparison with biomechanical tests in vitro. Comput Methods Biomech Biomed Engin. 2003;6(4):249-261.

16. Cowin S.C. Bone poroelasticity. J Biomech. 1999;32(3):217-238.

17. Shirazi-Adl A., Ahmed A.M., Shrivastava S.C. Mechanical response of a lumbar motion segment in axial torque alone and combined with compression. Spine (Phila Pa 1976). 1986;11(9):914-927.

18. Shirazi-Adl A., Ahmed A.M., Shrivastava S.C. A finite element study of a lumbar motion segment subjected to pure sagittal plane moments. J Biomech. 1986;19(4):331-350.

19. Kumaresan S., Yoganandan N., Pintar F.A. Finite element modeling approaches of human cervical spine facet joint capsule. J Biomech. 1998;31(4):371-376.

20. Wheeldon J.A., Stemper B.D., Yoganandan N., Pintar F.A. Validation of a finite element model of the young normal lower cervical spine. Ann Biomed Eng. 2008;36(9):1458-1469.

21. Guan Y., Yoganandan N., Maiman D.J., Pintar F.A. Internal and external responses of anterior lumbar/lumbosacral fusion: nonlinear finite element analysis. J Spinal Disord Tech. 2008;21(4):299-304.

22. Goel V.K., Kong W., Han J.S., et al. A combined finite element and optimization investigation of lumbar spine mechanics with and without muscles. Spine (Phila Pa 1976). 1993;18(11):1531-1541.

23. Goel V.K., Grauer J.N., Patel T., et al. Effects of charite artificial disc on the implanted and adjacent spinal segments mechanics using a hybrid testing protocol. Spine (Phila Pa 1976). 2005;30(24):2755-2764.

24. Rohlmann A., Zander T., Schmidt H., et al. Analysis of the influence of disc degeneration on the mechanical behaviour of a lumbar motion segment using the finite element method. J Biomech. 2006;39(13):2484-2490.

25. Dreischarf M., Rohlmann A., Bergmann G., Zander T. Optimised loads for the simulation of axial rotation in the lumbar spine. J Biomech. 2011;44(12):2323-2327.

26. Rohlmann A., Zander T., Rao M., Bergmann G. Realistic loading conditions for upper body bending. J Biomech. 2009;42(7):884-890.

27. Rohlmann A., Nabil Boustani H., Bergmann G., Zander T. Effect of a pedicle-screw-based motion preservation system on lumbar spine biomechanics: a probabilistic finite element study with subsequent sensitivity analysis. J Biomech. 2010;43(15):2963-2969.

28. Schmidt H., Heuer F., Wilke H.J. Interaction between finite helical axes and facet joint forces under combined loading. Spine (Phila Pa 1976). 2008;33(25):2741-2748.

29. Schmidt H., Heuer F., Wilke H.J. Which axial and bending stiffnesses of posterior implants are required to design a flexible lumbar stabilization system? J Biomech. 2009;42(1):48-54.

30. Schmidt H, Galbusera F, Rohlmann A, et al: Effect of multilevel lumbar disc arthroplasty on spine kinematics and facet joint loads in flexion and extension: a finite element analysis. Eur Spine J April 2, 2010 [Epub ahead of print].

31. Schmidt H., Midderhoff S., Adkins K., Wilke H.J. The effect of different design concepts in lumbar total disc arthroplasty on the range of motion, facet joint forces and instantaneous center of rotation of a L4-5 segment. Eur Spine J. 2009;18(11):1695-1705.

32. Zhong Z.C., Wei S.H., Wang J.P., et al. Finite element analysis of the lumbar spine with a new cage using a topology optimization method. Med Eng Phys. 2006;28(1):90-98.

33. Park W.M., Park Y.S., Kim K., Kim Y.H. Biomechanical comparison of instrumentation techniques in treatment of thoracolumbar burst fractures: a finite element analysis. J Orthop Sci. 2009;14(4):443-449.

34. Polikeit A., Nolte L.P., Ferguson S.J. The effect of cement augmentation on the load transfer in an osteoporotic functional spinal unit: finite-element analysis. Spine (Phila Pa 1976). 2003;28(10):991-996.

35. Lu Y.M., Hutton W.C., Gharpuray V.M. Do bending, twisting, and diurnal fluid changes in the disc affect the propensity to prolapse? A viscoelastic finite element model. Spine (Phila Pa 1976). 1996;21(22):2570-2579.

36. Kumaresan S., Yoganandan N., Pintar F.A., et al. Biomechanical study of pediatric human cervical spine: a finite element approach. J Biomech Eng. 2000;122(1):60-71.

37. Sharma M., Langrana N.A., Rodriguez J. Role of ligaments and facets in lumbar spinal stability. Spine (Phila Pa 1976). 1995;20(8):887-900.

38. Ruberte L.M., Natarajan R.N., Andersson G.B. Influence of single-level lumbar degenerative disc disease on the behavior of the adjacent segments–a finite element model study. J Biomech. 2009;42(3):341-348.

39. Hussain M., Nassr A., Natarajan R.N., et al. Biomechanical effects of anterior, posterior, and combined anterior-posterior instrumentation techniques on the stability of a multilevel cervical corpectomy construct: a finite element model analysis. Spine J. 2011;11(4):324-330.

40. Hussain M., Natarajan R.N., Chaudhary G., et al. Simulation of inhomogeneous rather than homogeneous poroelastic tissue material properties within disc annulus and nucleus better predicts cervical spine response: a C3-T1 finite element model analysis under compression and moment loadings. Spine (Phila Pa 1976). 2011;36(4):E245-E255.

41. Hussain M., Natarajan R.N., An H.S., Andersson G.B. Reduction in segmental flexibility because of disc degeneration is accompanied by higher changes in facet loads than changes in disc pressure: a poroelastic C5-C6 finite element investigation. Spine J. 2010;10(12):1069-1077.

42. Hussain M., Natarajan R.N., An H.S., Andersson G.B. Patterns of height changes in anterior and posterior cervical disc regions affects the contact loading at posterior facets during moderate and severe disc degeneration: a poroelastic C5-C6 finite element model study. Spine (Phila Pa 1976). 2010;35(18):E873-E881.

43. Patwardhan A.G., Havey R.M., Meade K.P., et al. A follower load increases the load-carrying capacity of the lumbar spine in compression. Spine (Phila Pa 1976). 1999;24(10):1003-1009.

44. Rohlmann A., Zander T., Bergmann G. Comparison of the biomechanical effects of posterior and anterior spine-stabilizing implants. Eur Spine J. 2005;14(5):445-453.

45. Shirazi-Adl S.A., Shrivastava S.C., Ahmed A.M. Stress analysis of the lumbar disc-body unit in compression. A three-dimensional nonlinear finite element study. Spine (Phila Pa 1976). 1984;9(2):120-134.

46. Polikeit A., Nolte L.P., Ferguson S.J. Simulated influence of osteoporosis and disc degeneration on the load transfer in a lumbar functional spinal unit. J Biomech. 2004;37(7):1061-1069.

47. Kim H.J., Chun H.J., Moon S.H., et al. Analysis of biomechanical changes after removal of instrumentation in lumbar arthrodesis by finite element analysis. Med Biol Eng Comput. 2010;48(7):703-709.

48. Rohlmann A., Zander T., Bergmann G. Effect of total disc replacement with ProDisc on intersegmental rotation of the lumbar spine. Spine (Phila Pa 1976). 2005;30(7):738-743.

49. Zander T., Rohlmann A., Burra N.K., Bergmann G. Effect of a posterior dynamic implant adjacent to a rigid spinal fixator. Clin Biomech (Bristol, Avon). 2006;21(8):767-774.

50. Baroud G., Nemes J., Heini P., Steffen T. Load shift of the intervertebral disc after a vertebroplasty: a finite-element study. Eur Spine J. 2003;12(4):421-426.

51. Rohlmann A., Zander T., Bergmann G. Spinal loads after osteoporotic vertebral fractures treated by vertebroplasty or kyphoplasty. Eur Spine J. 2006;15(8):1255-1264.

52. Brolin K., Halldin P. Development of a finite element model of the upper cervical spine and a parameter study of ligament characteristics. Spine (Phila Pa 1976). 2004;29(4):376-385.

53. Cheung J.T., Zhang M., Chow D.H. Biomechanical responses of the intervertebral joints to static and vibrational loading: a finite element study. Clin Biomech (Bristol, Avon). 2003;18(9):790-799.

54. Lavaste F., Skalli W., Robin, et al. Three-dimensional geometrical and mechanical modelling of the lumbar spine. J Biomech. 1992;25(10):1153-1164.

55. Wang J.L., Parnianpour M., Shirazi-Adl A., Engin A.E. Viscoelastic finite-element analysis of a lumbar motion segment in combined compression and sagittal flexion. Effect of loading rate. Spine (Phila Pa 1976). 2000;25(3):310-318.

56. Denoziere G., Ku D.N. Biomechanical comparison between fusion of two vertebrae and implantation of an artificial intervertebral disc. J Biomech. 2006;39(4):766-775.

57. Williams J.R., Natarajan R.N., Andersson G.B. Inclusion of regional poroelastic material properties better predicts biomechanical behavior of lumbar discs subjected to dynamic loading. J Biomech. 2007;40(9):1981-1987.

58. Anasetti F., Galbusera F., Aziz H.N., et al. Spine stability after implantation of an interspinous device: an in vitro and finite element biomechanical study. J Neurosurg Spine. 2010;13(5):568-575.

59. Kumaresan S., Yoganandan N., Pintar F.A. Finite element analysis of the cervical spine: a material property sensitivity study. Clin Biomech (Bristol, Avon). 1999;14(1):41-53.

60. Argoubi M., Shirazi-Adl A. Poroelastic creep response analysis of a lumbar motion segment in compression. J Biomech. 1996;29(10):1331-1339.

61. Rohlmann A., Bauer L., Zander T., et al. Determination of trunk muscle forces for flexion and extension by using a validated finite element model of the lumbar spine and measured in vivo data. J Biomech. 2006;39(6):981-989.

62. Jebaseelan D.D., Jebaraj C., Yoganandan N., Rajasekaran S. Validation efforts and flexibilities of an eight-year-old human juvenile lumbar spine using a three-dimensional finite element model. Med Biol Eng Comput. 2010;48(12):1223-1231.

63. Ng H.W., Teo E.C. Nonlinear finite-element analysis of the lower cervical spine (C4-C6) under axial loading. J Spinal Disord. 2001;14(3):201-210.

64. Little J.P., de Visser H., Pearcy M.J., Adam C.J. Are coupled rotations in the lumbar spine largely due to the osseo-ligamentous anatomy? A modeling study. Comput Methods Biomech Biomed Engin. 2008;11(1):95-103.

65. Lee S.H., Im Y.J., Kim K.T., et al. Comparison of cervical spine biomechanics after fixed- and mobile-core artificial disc replacement: a finite element analysis. Spine (Phila Pa 1976). 2011;36(9):700-708.

66. Yoganandan N., Kumaresan S.C., Voo L., et al. Finite element modeling of the C4-C6 cervical spine unit. Med Eng Phys. 1996;18(7):569-574.

67. Schmidt H., Shirazi-Adl A., Galbusera F., Wilke H.J. Response analysis of the lumbar spine during regular daily activities—a finite element analysis. J Biomech. 2010;43(10):1849-1856.

68. Galbusera F., Schmidt H., Noailly J., et al. Comparison of four methods to simulate swelling in poroelastic finite element models of intervertebral discs. J Mech Behav Biomed Mater. 2011;4(7):1234-1241.

69. Maurel N., Lavaste F., Skalli W. A three-dimensional parameterized finite element model of the lower cervical spine. Study of the influence of the posterior articular facets. J Biomech. 1997;30(9):921-931.

70. Schmidt H., Heuer F., Simon U., et al. Application of a new calibration method for a three-dimensional finite element model of a human lumbar annulus fibrosus. Clin Biomech (Bristol, Avon). 2006;21(4):337-344.

71. Noailly J., Lacroix D., Planell J.A. Finite element study of a novel intervertebral disc substitute. Spine (Phila Pa 1976). 2005;30(20):2257-2264.

72. Noailly J., Wilke H.J., Planell J.A., Lacroix D. How does the geometry affect the internal biomechanics of a lumbar spine bi-segment finite element model? Consequences on the validation process. J Biomech. 2007;40(11):2414-2425.

73. Lu Y.M., Hutton W.C., Gharpuray V.M. Can variations in intervertebral disc height affect the mechanical function of the disc? Spine (Phila Pa 1976). 1996;21(19):2208-2216. discussion 2217

74. Hussain M., Natarajan R.N., Chaudhary G., et al. Relative contributions of strain-dependent permeability and fixed charged density of proteoglycans in predicting cervical disc biomechanics: a poroelastic C5-C6 finite element model study. Med Eng Phys. 2011;33(4):438-445.

75. Goel V.K., Clausen J.D. Prediction of load sharing among spinal components of a C5-C6 motion segment using the finite element approach. Spine (Phila Pa 1976). 1998;23(6):684-691.

76. Malandrino A., Planell J.A., Lacroix D. Statistical factorial analysis on the poroelastic material properties sensitivity of the lumbar intervertebral disc under compression, flexion and axial rotation. J Biomech. 2009;42(16):2780-2788.

77. Ferguson S.J., Ito K., Nolte L.P. Fluid flow and convective transport of solutes within the intervertebral disc. J Biomech. 2004;37(2):213-221.

78. Schmidt H., Heuer F., Claes L., Wilke H.J. The relation between the instantaneous center of rotation and facet joint forces: a finite element analysis. Clin Biomech (Bristol, Avon). 2008;23(3):270-278.

79. Meijer G.J., Homminga J., Hekman E.E., et al. The effect of three-dimensional geometrical changes during adolescent growth on the biomechanics of a spinal motion segment. J Biomech. 2010;43(8):1590-1597.

80. Fagan M.J., Julian S., Siddall D.J., Mohsen A.M. Patient-specific spine models. Part 1: Finite element analysis of the lumbar intervertebral disc—a material sensitivity study. Proc Inst Mech Eng H. 2002;216(5):299-314.