A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

The human spinal segment is an inherently complex structure, a combination of flexible and semi-rigid articulating elements stabilised by seven principal ligaments. An understanding of how mechanical loading is shared among these passive elements of the segment is required to estimate tissue failure stresses. A 3D rigid body model of the complete lumbar spine has been developed to facilitate the prediction of load sharing across the passive elements. In contrast to previous multibody models, this model includes a non-linear, six degrees of freedom intervertebral disc, facet bony articulations and all spinal ligaments. Predictions of segmental kinematics and facet joint forces, in response to pure moment loading (flexion–extension), were compared to published in vitro data. On inclusion of detailed representation of the disc and facets, the multibody model fully captures the non-linear flexibility response of the spinal segment, i.e. coupled motions and a mobile instantaneous centre of rotation. Predicted facet joint forces corresponded well with reported values. For the loading case considered, the model predicted that the ligaments are the main stabilising elements within the physiological motion range; however, the disc resists a greater proportion of the applied load as the spine is fully flexed. In extension, the facets and capsular ligaments provide the principal resistance. Overall patterns of load distribution to the spinal ligaments are in agreement with previous predictions; however, the current model highlights the important role of the intraspinous ligament in flexion and the potentially high risk of failure. Several important refinements to the multibody modelling of the passive elements of the spine have been described, and such an enhanced passive model can be easily integrated into a full musculoskeletal model for the prediction of spinal loading for a variety of daily activities.     

Biomechanical implications of lumbar spinal ligament transection

Many lumbar spine surgeries either intentionally or inadvertently damage or transect spinal ligaments. The purpose of this work was to quantify the previously unknown biomechanical consequences of isolated spinal ligament transection on the remaining spinal ligaments (stress transfer), vertebrae (bone remodelling stimulus) and intervertebral discs (disc pressure) of the lumbar spine. A finite element model of the full lumbar spine was developed and validated against experimental data and tested in the primary modes of spinal motion in the intact condition. Once a ligament was removed, stress increased in the remaining spinal ligaments and changes occurred in vertebral strain energy, but disc pressure remained similar. All major biomechanical changes occurred at the same spinal level as the transected ligament, with minor changes at adjacent levels. This work demonstrates that iatrogenic damage to spinal ligaments disturbs the load sharing within the spinal ligament network and may induce significant clinically relevant changes in the spinal motion segment.     

A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

The human spinal segment is an inherently complex structure, a combination of flexible and semi-rigid articulating elements stabilised by seven principal ligaments. An understanding of how mechanical loading is shared among these passive elements of the segment is required to estimate tissue failure stresses. A 3D rigid body model of the complete lumbar spine has been developed to facilitate the prediction of load sharing across the passive elements. In contrast to previous multibody models, this model includes a non-linear, six degrees of freedom intervertebral disc, facet bony articulations and all spinal ligaments. Predictions of segmental kinematics and facet joint forces, in response to pure moment loading (flexion-extension), were compared to published in vitro data. On inclusion of detailed representation of the disc and facets, the multibody model fully captures the non-linear flexibility response of the spinal segment, i.e. coupled motions and a mobile instantaneous centre of rotation. Predicted facet joint forces corresponded well with reported values. For the loading case considered, the model predicted that the ligaments are the main stabilising elements within the physiological motion range; however, the disc resists a greater proportion of the applied load as the spine is fully flexed. In extension, the facets and capsular ligaments provide the principal resistance. Overall patterns of load distribution to the spinal ligaments are in agreement with previous predictions; however, the current model highlights the important role of the intraspinous ligament in flexion and the potentially high risk of failure. Several important refinements to the multibody modelling of the passive elements of the spine have been described, and such an enhanced passive model can be easily integrated into a full musculoskeletal model for the prediction of spinal loading for a variety of daily activities.     

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.     

Analysis of simulated single ligament transection on the mechanical behaviour of a lumbar functional spinal unit.

The influence of the different lumbar spinal ligaments on intersegmental rotation is not fully understood. In order to explore this effect, a finite element model of the functional spinal unit L3/L4 was loaded with pure moments in the three main anatomic planes. The two extremes--minimum and maximum--ligament stiffness values reported in the literature were applied. After virtual transection of each of the spinal ligaments in turn, the intersegmental rotation and forces in the remaining ligaments were calculated. On flexion, the highest force was found for the posterior longitudinal ligament; on extension and lateral bending for the anterior longitudinal ligament; and on axial rotation for the facet capsular ligament. The strongest influence on intersegmental rotation is exerted by the interspinous ligament on flexion, by the anterior longitudinal ligament on extension and lateral bending, and by the facet capsular ligaments on axial rotation. Ligament stiffness has a strong influence on intersegmental rotation and forces in the ligaments, so that finite element models of spinal segments must be validated by experimental data. This study should help to elucidate the role of the various ligaments.     

Spinal facet joint biomechanics and mechanotransduction in normal, injury and degenerative conditions

The facet joint is a crucial anatomic region of the spine owing to its biomechanical role in facilitating articulation of the vertebrae of the spinal column. It is a diarthrodial joint with opposing articular cartilage surfaces that provide a low friction environment and a ligamentous capsule that encloses the joint space. Together with the disc, the bilateral facet joints transfer loads and guide and constrain motions in the spine due to their geometry and mechanical function. Although a great deal of research has focused on defining the biomechanics of the spine and the form and function of the disc, the facet joint has only recently become the focus of experimental, computational and clinical studies. This mechanical behavior ensures the normal health and function of the spine during physiologic loading but can also lead to its dysfunction when the tissues of the facet joint are altered either by injury, degeneration or as a result of surgical modification of the spine. The anatomical, biomechanical and physiological characteristics of the facet joints in the cervical and lumbar spines have become the focus of increased attention recently with the advent of surgical procedures of the spine, such as disc repair and replacement, which may impact facet responses. Accordingly, this review summarizes the relevant anatomy and biomechanics of the facet joint and the individual tissues that comprise it. In order to better understand the physiological implications of tissue loading in all conditions, a review of mechanotransduction pathways in the cartilage, ligament and bone is also presented ranging from the tissue-level scale to cellular modifications. With this context, experimental studies are summarized as they relate to the most common modifications that alter the biomechanics and health of the spine-injury and degeneration. In addition, many computational and finite element models have been developed that enable more-detailed and specific investigations of the facet joint and its tissues than are provided by experimental approaches and also that expand their utility for the field of biomechanics. These are also reviewed to provide a more complete summary of the current knowledge of facet joint mechanics. Overall, the goal of this review is to present a comprehensive review of the breadth and depth of knowledge regarding the mechanical and adaptive responses of the facet joint and its tissues across a variety of relevant size scales.     

Load-sharing in the lumbosacral spine in neutral standing & flexed postures – A combined finite element and inverse static study

Understanding load-sharing in the spine during in-vivo conditions is critical for better spinal implant design and testing. Previous studies of load-sharing that considered actual spinal geometry applied compressive follower load, with or without moment, to simulate muscle forces. Other studies used musculoskeletal models, which include muscle forces, but model the discs by simple beams or spherical joints and ignore the articular facet joints.This study investigated load-sharing in neutral standing and flexed postures using a detailed Finite Element (FE) model of the ligamentous lumbosacral spine, where muscle forces, gravity loads and intra-abdominal pressure, as predicted by a musculoskeletal model of the upper body, are input into the FE model. Flexion was simulated by applying vertebral rotations following spine rhythm measured in a previous in-vivo study, to the musculoskeletal model. The FE model predicted intradiscal pressure (IDP), strains in the annular fibers, contact forces in the facet joints, and forces in the ligaments. The disc forces and moments were determined using equilibrium equations, which considered the applied loads, including muscle forces and IDP, as well as forces in the ligaments and facet joints predicted by the FE model. Load-sharing was calculated as the portion of the total spinal load carried along the spine by each individual spinal structure. The results revealed that spinal loads which increased substantially from the upright to the flexed posture were mainly supported by the discs in the upright posture, whereas the ligaments’ contribution in resisting shear, compression, and moment was more significant in the flexed posture.     

Diurnal variations in intervertebral disc height affect spine flexibility,intradiscal pressure and contact compressive forces in the facet joints

Diurnal changes of intervertebral disc height are caused by high compressive loading during the day, which expulses fluid from the disc, and by osmotic pressure, which imbibes fluid into the disc at low loading. The aim of the present study was to determine the magnitude of diurnal changes in spine flexibility, intradiscal pressures and contact forces in the facet joints. A validated osseoligamentous finite element model of the lumbar spine was used to determine these quantities for morning and evening situations. Disc height varied by 10% for these two situations. Spine flexibility and facet joint forces were markedly higher in the evening than in the morning. Intradiscal pressures were higher in the morning than in the evening. The different spine flexibilities in the morning and evening should be taken into account during kinematical measurements. Predicted facet joint forces may be used for the designing and pre-clinical testing of artificial facet joint replacements.     

Three-dimensional static modeling of the lumbar spine

This paper presents three-dimensional static modeling of the human lumbar spine to be used in the formation of anatomically-correct movement patterns for a fully cable-actuated robotic lumbar spine which can mimic in vivo human lumbar spine movements to provide better hands-on training for medical students. The mathematical model incorporates five lumbar vertebrae between the first lumbar vertebra and the sacrum, with dimensions of an average adult human spine. The vertebrae are connected to each other by elastic elements, torsional springs and a spherical joint located at the inferoposterior corner in the mid-sagittal plane of the vertebral body. Elastic elements represent the ligaments that surround the facet joints and the torsional springs represent the collective effect of intervertebral disc which plays a major role in balancing torsional load during upper body motion and the remaining ligaments that support the spinal column. The elastic elements and torsional springs are considered to be nonlinear. The nonlinear stiffness constants for six motion types were solved using a multiobjective optimization technique. The quantitative comparison between the angles of rotations predicted by the proposed model and in the experimental data confirmed that the model yields angles of rotation close to the experimental data. The main contribution is that the new model can be used for all motions while the experimental data was only obtained at discrete measurement points.     

Stepwise reduction of functional spinal structures increase range of motion and change lordosis angle

Many investigators have performed studies on specific defect situations or determined the contribution on isolated structures. Investigating the contribution of functional structures requires obtaining the kinematic response directly on spinal segments. The purpose of this study was to quantify the function of anatomical components on lumbar segments for different loading magnitudes. Eight spinal segments (L4-5) with a median age of 52 years (ranging from 38 to 59 years) and a low degree of disc degeneration were utilized for the in vitro testing. Specimens were mounted in a custom-built spine tester and loaded with pure moments (1-10 N m) to move within three anatomical planes at a loading rate of 1.0 degrees /s. Anatomy was successively reduced by: ligaments, facet capsules, joints and nucleus. Data were evaluated for range of motion, neutral zone and lordosis angle. Transection of posterior ligaments predominantly increased specimen flexion for all bending moments applied. Supraspinous ligament also indicated to resist in extension slightly, whereas the facet capsules did not. Facet joints contributed to axial rotation, but not in lateral bending. The anterior longitudinal ligament was found to slightly resist in axial rotation, but strongly in extension. Nucleotomy caused largest increase of all movements. The unloaded posture of the specimens changed after ligament dissection, indicating ligament pretension. The region of lumbar spine is interesting for finite element (FE) simulation due to the high evidence of disc degeneration and injuries. This study may help to understand the function of specific anatomical structures and assists in FE model calibration. We suggest to start a calibration procedure for such models with the smallest functional structure (annulus) and to cumulatively add further structures.     

Investigation of optimal follower load path generated by trunk muscle coordination

It has been reported that the center of rotation of each vertebral body is located posterior to the vertebral body center. Moreover, it has been suggested that an optimized follower load (FL) acts posterior to the vertebral body center. However, the optimal position of the FL with respect to typical biomechanical characteristics regarding spinal stabilization, such as joint compressive force, shear force, joint moment, and muscle stress, has not been studied. A variation in the center of rotation of each vertebra was formulated in a three-dimensional finite element model of the lumbar spine with 117 pairs of trunk muscles. Then, the optimal translation of the FL path connecting the centers of rotations was estimated by solving the optimization problem that was to simultaneously minimize the compressive forces, the shear forces, and the joint moments or to minimize the cubic muscle stresses. An upright neutral standing position and a standing position with 200N in both hands were considered. The FL path moved posterior, regardless of the optimization criteria and loading conditions. The FL path moved 5.0 and 7.8mm posterior in upright standing and 4.1mm and 7.0mm posterior in standing with 200N in hands for each optimization scheme. In addition, it was presented that the optimal FL path may have advantages in comparison to the body center FL path. The present techniques may be important in understanding the spine stabilization function of the trunk muscles.     

Effects of eight different ligament property datasets on biomechanics of a lumbar L4-L5 finite element model

Ligaments assist trunk muscles in balancing external moments and providing spinal stability. In absence of the personalized material properties for ligaments, finite element (FE) models use dispersed data from the literature. This study aims to investigate the relative effects of eight different ligament property datasets on FE model responses. Eight L4-L5 models distinct only in ligament properties were constructed and loaded under moment (15 N m) alone or combined with a compressive follower load (FL). Range of motions (RoM) of the disc-alone model matched well in vitro data. Ligament properties significantly affected only sagittal RoMs (∼3.0–7.1° in flexion and ∼3.8–5.8° in extension at 10 N m). Sequential removal of ligaments shifted sagittal RoMs in and out of the corresponding in vitro ranges. When moment was combined with FL, center of rotation matched in vivo data for all models (3.8 ± 0.9 mm and 4.3 ± 1.8 mm posterior to the disc center in flexion and extension, respectively). Under 15 N m sagittal moments, ligament strains were often smaller or within the in vitro range in flexion whereas some posterior ligament forces approached their failure forces in some models. Ligament forces varied substantially within the models and affected the moment-sharing and internal forces on the disc and facet joints. Intradiscal pressure (IDP) had the greatest variation between models in extension. None of the datasets yielded results in agreement with all reported measurements. Results emphasized the important role of ligaments especially under larger moments and the need for their accurate representation in search for valid spinal models.     

Numerical analysis of the influence of nucleus pulposus removal on the biomechanical behavior of a lumbar motion segment

Nucleus replacement was deemed to have therapeutic potential for patients with intervertebral disc herniation. However, whether a patient would benefit from nucleus replacement is technically unclear. This study aimed to investigate the influence of nucleus pulposus (NP) removal on the biomechanical behavior of a lumbar motion segment and to further explore a computational method of biomechanical characteristics of NP removal, which can evaluate the mechanical stability of pulposus replacement. We, respectively, reconstructed three types of models for a mildly herniated disc and three types of models for a severely herniated disc based on a L4–L5 segment finite element model with computed tomography image data from a healthy adult. First, the NP was removed from the herniated disc models, and the biomechanical behavior of NP removal was simulated. Second, the NP cavities were filled with an experimental material (Poisson's ratio = 0.3; elastic modulus = 3 MPa), and the biomechanical behavior of pulposus replacement was simulated. The simulations were carried out under the five loadings of axial compression, flexion, lateral bending, extension, and axial rotation. The changes of the four biomechanical characteristics, i.e. the rotation degree, the maximum stress in the annulus fibrosus (AF), joint facet contact forces, and the maximum disc deformation, were computed for all models. Experimental results showed that the rotation range, the maximum AF stress, and joint facet contact forces increased, and the maximum disc deformation decreased after NP removal, while they changed in the opposite way after the nucleus cavities were filled with the experimental material.     

Geometric and mechanical properties of human cervical spine ligaments

This study characterized the geometry and mechanical properties of the cervical ligaments from C2-T1 levels. The lengths and cross-sectional areas of the anterior longitudinal ligament, posterior longitudinal ligament, joint capsules, ligamentum flavum, and interspinous ligament were determined from eight human cadavers using cryomicrotomy images. The geometry was defined based on spinal anatomy and its potential use in complex mathematical models. The biomechanical force-deflection, stiffness, energy, stress, and strain data were obtained from 25 cadavers using in situ axial tensile tests. Data were grouped into middle (C2-C5) and lower (C5-T1) cervical levels. Both the geometric length and area of cross section, and the biomechanical properties including the stiffness, stress, strain, energy, and Young's modulus, were presented for each of the five ligaments. In both groups, joint capsules and ligamentum flavum exhibited the highest cross-sectional area (p < 0.005), while the longitudinal ligaments had the highest length measurements. Although not reaching statistical significance, for all ligaments, cross-sectional areas were higher in the C5-T1 than in the C2-C5 group; and lengths were higher in the C2-C5 than in the C5-T1 group with the exception of the flavum (Table 1 in the main text). Force-deflection characteristics (plots) are provided for all ligaments in both groups. Failure strains were higher for the ligaments of the posterior (interspinous ligament, joint capsules, and ligamentum flavum) than the anterior complex (anterior and posterior longitudinal ligaments) in both groups. In contrast, the failure stress and Young's modulus were higher for the anterior and posterior longitudinal ligaments compared to the ligaments of the posterior complex in the two groups. However, similar tendencies in the structural responses (stiffness, energy) were not found in both groups. Researchers attempting to incorporate these data into stress-analysis models can choose the specific parameter(s) based on the complexity of the model used to study the biomechanical behavior of the human cervical spine.     

Cruciate coupling and screw-home mechanism in passive knee joint during extension--flexion

The screw-home mechanism and coupling between forces in cruciate ligaments during passive knee joint flexion were investigated for various boundary conditions, flexion axis alignments and posterior cruciate ligaments (PCL)/anterior cruciate ligament (ACL) conditions. A developed non-linear 3D finite element model was used to perform detailed elasto-static response analyses of the human tibiofemoral joint as a function of flexion angle varying from 10 degrees hyper-extension to 90 degrees flexion. The tibia rotated internally as the femur flexed and externally as the femur extended. The re-alignment of the flexion axis by +/-5 degrees rotation about the axial (distal-proximal) axis, transection of the ACL and changes in cruciate ligament initial strains substantially influenced the 'screw-home' motion. On the other hand, restraint on this coupled rotation diminished ACL forces in flexion. A remarkable coupling was predicted between ACL and PCL forces in flexion; forces in both cruciate ligaments increased as the initial strain or pretension in one of them increased whereas they both diminished as one of them was cut or became slack. This has important consequences in joint functional biomechanics following a ligament injury or replacement surgery and, hence, in the proper management of joint disorders.     

Analysis of the interaction between vertebral lateral deviation and axial rotation in scoliosis

There is a lack of clear biomechanical analyses to explain the interaction of the lateral and axial deformity of the spine in idiopathic scoliosis. A finite element model which represented an isolated ligamentous spine with realistic elastic properties and idealized geometry was used to analyse this interaction. Three variations of this model were used to investigate two different hypotheses about the etiology of scoliosis and to define the forces required to produce a scoliosis deformity. The first hypothesis is that coupling within a motion segment produces the interaction between lateral deviation and axial rotation. The second hypothesis is that posterior tethering by soft tissues in the growing spine produces the observed interaction. Modeling of both hypotheses failed to produce the clinically observed pattern of interaction. Therefore, to find which biomechanical forces were required to produce an idealized scoliosis, prescribed displacements were applied to the model. Production of a double curve scoliosis of 10 degrees Cobb angles required lateral forces on the order of 20 N acting 40 mm anterior to the vertebral body centers. There do not appear to be any anatomic structures capable of producing such forces. Therefore, it seems unlikely that scoliosis deformity can be explained in terms of forces acting on the spine, and understanding of its origins may come from examination of other mechanisms such as asymmetric thoracic growth, or asymmetric vertebral development.     

Effects of lumbo-pelvic rhythm on trunk muscle forces and disc loads during forward flexion: A combined musculoskeletal and finite element simulation study

Previous in-vivo studies suggest that the ratio of total lumbar rotation over pelvic rotation (lumbo-pelvic rhythm) during trunk sagittal movement is essential to evaluate spinal loads and discriminate between low back pain and asymptomatic population. Similarly, there is also evidence that the lumbo-pelvic rhythm is key for evaluation of realistic muscle and joint reaction forces and moments predicted by various computational musculoskeletal models. This study investigated the effects of three lumbo-pelvic rhythms defined based on in-vivo measurements on the spinal response during moderate forward flexion (60°) using a combined approach of musculoskeletal modeling of the upper body and finite element model of the lumbosacral spine. The muscle forces and joint loads predicted by the musculoskeletal model, together with the gravitational forces, were applied to the finite element model to compute the disc force and moment, intradiscal pressure, annular fibers strain, and load-sharing. The results revealed that a rhythm with high pelvic rotation and low lumbar flexion involves more global muscles and increases the role of the disc in resisting spinal loads, while its counterpart, with low pelvic rotation, recruits more local muscles and engages the ligaments to lower the disc loads. On the other hand, a normal rhythm that has balanced pelvic and lumbar rotations yields almost equal disc and ligament load-sharing and results in more balanced synergy between global and local muscles. The lumbo-pelvic rhythm has less effect on the intradiscal pressure and annular fibers strain. This work demonstrated that the spinal response during forward flexion is highly dependent on the lumbo-pelvic rhythm. It is therefore, essential to adapt this parameter instead of using the default values in musculoskeletal models for accurate prediction of muscle forces and joint reaction forces and moments. The findings provided by this work are expected to improve knowledge of spinal response during forward flexion, and are clinically relevant towards low back pain treatment and disc injury prevention.     

Effect of an interspinous implant on loads in the lumbar spine.

Interspinous process implants are increasingly used to treat canal stenoses. Little information exists about the effects of implant height and stiffness on the biomechanical behavior of the lumbar spine. Therefore, a three-dimensional nonlinear finite element model of the osseoligamentous lumbar spine (L1 to L5) was created with a slightly degenerated disc at L3/L4. An interspinous implant was inserted at that segment. Implants with different heights and stiffnesses were studied. The model was loaded with the upper body weight and muscle forces to simulate walking and 25 degrees extension. Implant forces are influenced strongly by the height and negligibly by the elastic modulus of the implant. Intersegmental rotation at implant level is markedly reduced, while intradiscal pressure is slightly increased. Implant size and stiffness have only a minor effect on intradiscal pressure. The maximum von Mises stress in the vertebral arch is strongly increased by the implant.     

Investigation of impact loading rate effects on the ligamentous cervical spinal load-partitioning using finite element model of functional spinal unit C2–C3

The cervical spine functions as a complex mechanism that responds to sudden loading in a unique manner, due to intricate structural features and kinematics. The spinal load-sharing under pure compression and sagittal flexion/extension at two different impact rates were compared using a bio-fidelic finite element (FE) model of the ligamentous cervical functional spinal unit (FSU) C2–C3. This model was developed using a comprehensive and realistic geometry of spinal components and material laws that include strain rate dependency, bone fracture, and ligament failure. The range of motion, contact pressure in facet joints, failure forces in ligaments were compared to experimental findings. The model demonstrated that resistance of spinal components to impact load is dependent on loading rate and direction. For the loads applied, stress increased with loading rate in all spinal components, and was concentrated in the outer intervertebral disc (IVD), regions of ligaments to bone attachment, and in the cancellous bone of the facet joints. The highest stress in ligaments was found in capsular ligament (CL) in all cases. Intradiscal pressure (IDP) in the nucleus was affected by loading rate change. It increased under compression/flexion but decreased under extension. Contact pressure in the facet joints showed less variation under compression, but increased significantly under flexion/extension particularly under extension. Cancellous bone of the facet joints region was the only component fractured and fracture occurred under extension at both rates. The cervical ligaments were the primary load-bearing component followed by the IVD, endplates and cancellous bone; however, the latter was the most vulnerable to extension as it fractured at low energy impact.     

The effect of three-dimensional geometrical changes during adolescent growth on the biomechanics of a spinal motion segment

During adolescent growth, vertebrae and intervertebral discs undergo various geometrical changes. Although such changes in geometry are well known, their effects on spinal stiffness remains poorly understood. However, this understanding is essential in the treatment of spinal abnormalities during growth, such as scoliosis.A finite element model of an L3–L4 motion segment was developed, validated and applied to study the quantitative effects of changing geometry during adolescent growth on spinal stiffness in flexion, extension, lateral bending and axial rotation. Height, width and depth of the vertebrae and intervertebral disc were varied, as were the width of the transverse processes, the length of the spinous process, the size of the nucleus, facet joint areas and ligament size. These variations were based on average growth data for girls, as reported in literature.Overall, adolescent growth increases the stiffness with 36% (lateral bending and extension) to 44% (flexion). Two thirds of this increase occurs between 10 and 14 years of age and the last third between 14 years of age and maturity.Although the height is the largest geometrical change during adolescent growth, its effect on the biomechanics is small. The depth increase of the disc and vertebrae significantly affects the stiffness in all directions, while the width increase mainly affects the lateral bending stiffness. Hence, when analysing the biomechanics of the growing adolescent spine (for instance in scoliosis research), the inclusion of depth and width changes, in addition to the usually implemented height change, is essential.