The viscoelastic standard nonlinear solid model: predicting the response of the lumbar intervertebral disk to low-frequency vibrations

Due to the mathematical complexity of current musculoskeletal spine models, there is a need for computationally efficient models of the intervertebral disk (IVD). The aim of this study is to develop a mathematical model that will adequately describe the motion of the IVD under axial cyclic loading as well as maintain computational efficiency for use in future musculoskeletal spine models. Several studies have successfully modeled the creep characteristics of the IVD using the three-parameter viscoelastic standard linear solid (SLS) model. However, when the SLS model is subjected to cyclic loading, it underestimates the load relaxation, the cyclic modulus, and the hysteresis of the human lumbar IVD. A viscoelastic standard nonlinear solid (SNS) model was used to predict the response of the human lumbar IVD subjected to low-frequency vibration. Nonlinear behavior of the SNS model was simulated by a strain-dependent elastic modulus on the SLS model. Parameters of the SNS model were estimated from experimental load deformation and stress-relaxation curves obtained from the literature. The SNS model was able to predict the cyclic modulus of the IVD at frequencies of 0.01 Hz, 0.1 Hz, and 1 Hz. Furthermore, the SNS model was able to quantitatively predict the load relaxation at a frequency of 0.01 Hz. However, model performance was unsatisfactory when predicting load relaxation and hysteresis at higher frequencies (0.1 Hz and 1 Hz). The SLS model of the lumbar IVD may require strain-dependent elastic and viscous behavior to represent the dynamic response to compressive strain.     

A comparison between porcine,ovine, and bovine intervertebral disc anatomy and single lamella annulus fibrosus tensile properties

This project aimed to compare gross anatomical measures and biomechanical properties of single lamellae from the annulus fibrosus of ovine and porcine lumbar vertebrae, and bovine tail vertebrae. The morphology of the vertebrae of these species differ significantly both from each other and from human, yet how these differences alter biomechanical properties is unknown. Geometric parameters measured in this study included: 1) absolute and relative intervertebral (IVD) and vertebral body height and 2) absolute and relative intervertebral disc (IVD) anterior‐posterior (AP) and medial‐lateral (ML) widths. Single lamella tensile properties included toe‐region stress and stretch ratio, stiffness, and tensile strength. As expected, the bovine tail IVD revealed a more circular shape compared with both the ovine and porcine lumbar IVD. The bovine tail also had the largest IVD to vertebral body height ratio (due to having the highest absolute IVD height). Bovine tail lamellae were also found to be strongest and stiffest (in tension) while ovine lumbar lamellae were weakest and most compliant. Histological analysis revealed the greatest proportion of collagen in the bovine corroborating findings of increased strength and stiffness. The observed differences in anatomical shape, connective tissue composition, and tensile properties need to be considered when choosing an appropriate model for IVD research. J. Morphol. 277:244–251, 2016. © 2015 Wiley Periodicals, Inc.     

The biological basis of degenerative disc disease: proteomic and biomechanical analysis of the canine intervertebral disc

IntroductionIn the present study, we sought to quantify and contrast the secretome and biomechanical properties of the non-chondrodystrophic (NCD) and chondrodystrophic (CD) canine intervertebral disc (IVD) nucleus pulposus (NP).MethodsWe used iTRAQ proteomic methods to quantify the secretome of both CD and NCD NP. Differential levels of proteins detected were further verified using immunohistochemistry, Western blotting, and proteoglycan extraction in order to evaluate the integrity of the small leucine-rich proteoglycans (SLRPs) decorin and biglycan. Additionally, we used robotic biomechanical testing to evaluate the biomechanical properties of spinal motion segments from both CD and NCD canines.ResultsWe detected differential levels of decorin, biglycan, and fibronectin, as well as of other important extracellular matrix (ECM)-related proteins, such as fibromodulin and HAPLN1 in the IVD NP obtained from CD canines compared with NCD canines. The core proteins of the vital SLRPs decorin and biglycan were fragmented in CD NP but were intact in the NP of the NCD animals. CD and NCD vertebral motion segments demonstrated significant differences, with the CD segments having less stiffness and a more varied range of motion.ConclusionsThe CD NP recapitulates key elements of human degenerative disc disease. Our data suggest that at least some of the compromised biomechanical properties of the degenerative disc arise from fibrocartilaginous metaplasia of the NP secondary to fragmentation of SLRP core proteins and associated degenerative changes affecting the ECM. This study demonstrates that the degenerative changes that naturally occur within the CD NP make this animal a valuable animal model with which to study IVD degeneration and potential biological therapeutics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0733-z) contains supplementary material, which is available to authorized users.     

TNFα Transport Induced by Dynamic Loading Alters Biomechanics of Intact Intervertebral Discs

Objective

Intervertebral disc (IVD) degeneration is an important contributor to the development of back pain, and a key factor relating pain and degeneration are the presence of pro-inflammatory cytokines and IVD motion. There is surprisingly limited understanding of how mechanics and inflammation interact in the IVD. This study investigated interactions between mechanical loading and pro-inflammatory cytokines in a large animal organ culture model to address fundamental questions regarding (i.) how inflammatory mediators arise within the IVD, (ii.) how long inflammatory mediators persist, and (iii.) how inflammatory mediators influence IVD biomechanics.

Methods

Bovine caudal IVDs were cultured for 6 or 20-days under static & dynamic loading with or without exogenous TNFα in the culture medium, simulating a consequence of inflammation of the surrounding spinal tissues. TNFα transport within the IVD was assessed via immunohistochemistry. Changes in IVD structural integrity (dimensions, histology & aggrecan degradation), biomechanical behavior (Creep, Recovery & Dynamic stiffness) and pro-inflammatory cytokines in the culture medium (ELISA) were assessed.

Results

TNFα was able to penetrate intact IVDs when subjected to dynamic loading but not static loading. Once transported within the IVD, pro-inflammatory mediators persisted for 4–8 days after TNFα removal. TNFα exposure induced changes in IVD biomechanics (reduced diurnal displacements & increased dynamic stiffness).

Discussion

This study demonstrated that exposure to TNFα, as might occur from injured surrounding tissues, can penetrate healthy intact IVDs, induce expression of additional pro-inflammatory cytokines and alter IVD mechanical behavior. We conclude that exposure to pro-inflammatory cytokine may be an initiating event in the progression of IVD degeneration in addition to being a consequence of disease.     

Effect of microgravity on the biomechanical properties of lumbar and caudal intervertebral discs in mice

Prolonged exposure to microgravity has shown to have deleterious effects on the human spine, indicated by low back pain during spaceflight and increased incidence of post-spaceflight herniated nucleus pulposus. We examined the effect of microgravity on biomechanical properties of lumbar and caudal discs from mice having been on 15-day shuttle mission STS-131. Sixteen C57BL/C mice (spaceflight group, n=8; ground-based control group, n=8) were sacrificed immediately after spaceflight. Physiological disc height (PDH) was measured in situ, and compressive creep tests were performed to parameterize biomechanical properties into endplate permeability (k), nuclear swelling pressure strain dependence (D), and annular viscoelasticity (G). For caudal discs, the spaceflight group exhibited 32% lower PDH, 70% lower D and crept more compared to the control mice (p=0.03). For lumbar discs, neither PDH nor D was significantly different between murine groups. Initial modulus, osmotic pressure, k and G for lumbar and caudal discs did not appear influenced by microgravity (p>0.05). Decreases in both PDH and D suggest prolonged microgravity effectively diminished biomechanical properties of caudal discs. By contrast, differences were not noted for lumbar discs. This potentially deleterious interaction between prolonged weightlessness and differential ranges of motion along the spine may underlie the increased cervical versus lumbar disc herniation rates observed among astronauts.     

Untersuchungen zur Regeneration des Hinterendes beiAnaitides mucosa (Polychaeta,Phyllodocidae)

Caudal regeneration was investigated in decerebrateAnaitides mucosa and in brain-intact individuals. Both groups show an identical capacity to regenerate lost caudal segments. Furthermore there is no difference in males and females. Low temperature (5 °C) inhibits the regeneration of caudal segments, but it is necessary for normal oogenesis. Under conditions of high temperature (15 °C), caudal regeneration is very extensive. At the same time degeneration of most oocytes occurs.     

Neuromuscular manifestations of viscoelastic tissue degradation following high and low risk repetitive lumbar flexion

Cumulative lumbar disorder is common in individuals engaged in long term performance of repetitive and static occupational/sports activities with the spine. The triggering source and of the disorder, the tissues involved in the failure and the biomechanical, neuromuscular, and biological processes active in the initiation and development of the disorder are not known. The hypothesis is forwarded that static and repetitive (cyclic) lumbar flexion-extension and the associated repeated stretch of the various viscoelastic tissues (ligaments, fascia, facet capsule, discs, etc.) causes micro-damage in their collagen fibers followed by an acute inflammation, triggering pain and reflexive muscle spasms/hyper-excitability. Continued exposure to activities, over time, converts the acute inflammation into a chronic one, viscoelastic tissues remodeling/degeneration, modified motor control strategy and permanent disability. Changes in lumbar stability are expected during the development of the disorder. A series of experimental data from in-vivo feline is reviewed and integrated with supporting evidence from the literature to gain a valuable insight into the multi-factorial development of the disorder. Prolonged cyclic lumbar flexion-extension at high loads, high velocities, many repetitions and short in between rest periods induced transient creep/laxity in the spine, muscle spasms and reduced stability followed, several hours later, by an acute inflammation/tissue degradation, muscular hyper-excitability and increased stability. The major findings assert that viscoelastic tissues sub-failure damage is the source and inflammation is the process which governs the mechanical and neuromuscular characteristic symptoms of the disorder. A comprehensive model of the disorder is presented. The experimental data validates the hypothesis as well as provide insights into the development of potential treatment and prevention of the disorder.     

Role of endplates in contributing to compression behaviors of motion segments and intervertebral discs

The purpose of this study was to gain an improved understanding of the mechanical behavior of the intervertebral disc in the presence and absence of the vertebral endplates. Mechanical behaviors of rat caudal motion segments, vertebrae and isolated disc explants under two different permeability conditions were investigated and viscoelastic behaviors were evaluated using a stretched-exponential function to describe creep and recovery behaviors. The results demonstrated that both vertebrae and discs underwent significant deformations in the motion segment even under relatively low-loading conditions. Secondly, disruption of the collagenous network had minimal impact on equilibrium deformations of disc explants as compared to disc deformations occurring in the motion segments provided that vertebral deformations were accounted for; however, differences in endplate permeability conditions had a significant effect on viscoelastic behaviors. Creep occurred more quickly than recovery for motion segment and explant specimens. In addition, disc explants and motion segments both exhibited non-recoverable deformations under axial compression under low- and high-loading conditions. Results have important implications for interpreting the role of vertebral endplates in contributing to disc mechanical behaviors and direct application to mechanobiology studies involving external loading to rodent tail intervertebral discs.     

Effect of overload on changes in mechanical and structural properties of the annulus fibrosus of the intervertebral disc

The research focussed on analysing structural and mechanical properties in the intervertebral disc (IVD), caused by long-term cyclic loading. Spinal motion segments were divided into two groups: the control (C), and the group in which it was analysed the impact of posterior column in the load-bearing system of the spine—specimens with intact posterior column (IPC) and without posterior column (WPC). To evaluate the structural and mechanical changes, the specimens were tested with simulation of 100,000 compression-flexion load cycles after which it was performed macroscopic analysis. Mechanical properties of the annulus fibrosis (AF) from the anterior and posterior regions of the IVD were tested at the uniaxial tension test. The stiffness coefficient values were statistically 32% higher in the WPC group (110 N/mm) than in the IPC (79 N/mm). The dynamics of increase in this parameter does not correspond with the course of decrease in height loss. WPC segments revealed clear structural changes that mainly involve the posterior regions of the IVD (bulging and delamination with the effect of separation of collagen fibre bundles). Pathological changes also caused decreases in the value of stress in the AF. The greatest changes in the stress value about group C (7.43 ± 4.49 MPa) were observed in the front part of the fibrous ring, where this value was for IPC 4.49 ± 4.78 MPa and WPC 2.56 ± 1.01 MPa. The research indicates that the applied load model allows simulating damage that occurs in pathological IVD. And the posterior column’s presence affects this change’s dynamics, structural and mechanical properties of AF.

     

Flexible device for vertebral body replacement

A novel vertebral prosthesis is presented. The prosthesis was developed for surgical procedures requiring the resection of a complete vertebral body and the adjacent intervertebral discs, the design objective being to develop a flexible implant that would be robust enought to withstand the in vivo stress environment of the human spine. In theory, a flexible implant should preserve a more normal range of motion and apply less stress to surrounding tissue than a rigid implant. A prototype implant was constructed so as to combine a rigid stainless steel structure with flexible silicon rubber elements in order to form an implant with static and dynamic mechanical characteristics similar to those of the anterior spinal column. Implant flexibility characteristics were determined from ex vivo stress-strain behaviour during bending and compressive creep testing. Results from the bending tests indicated good agreement for the lateral and sagittal bending characteristics in comparison with in vitro bending tests of human lumbar motion segments. Comparison of the implant compressive creep responsé with similar in vitro tests on human lumbar intervertebral discs also demonstrated similarities in the time-dependent mechanical parameters.     

Influence of single-level lumbar degenerative disc disease on the behavior of the adjacent segments—A finite element model study

The current study investigated mechanical predictors for the development of adjacent disc degeneration. A 3-D finite element model of a lumbar spine was modified to simulate two grades of degeneration at the L4–L5 disc. Degeneration was modeled by changes in geometry and material properties. All models were subjected to follower preloads of 800 N and moment loads in the three principal directions of motion using a hybrid protocol. Degeneration caused changes in the loading and motion patterns of the segments above and below the degenerated disc. At the level (L3–L4) above the degenerated disc, the motion increased due to moderate degeneration by 21% under lateral bending, 26% under axial rotation and 28% under flexion/extension. At the level (L5-S1) below the degenerated disc, motion increased only during lateral bending by 20% due to moderate degeneration. Both the L3–L4 and L5-S1 segment showed a monotonic increase in both the maximum von Mises stress and shear stress in the annulus as degeneration progressed for all loading directions, expect extension at L3–L4. The most significant increase in stress was observed at the L5-S1 level during axial rotation with nearly a ten-fold increase in the maximum shear stress and 103% increase in the maximum von Mises stress. The L5-S1 segment also showed a progressive increase in facet contact force for all loading directions with degeneration. Nucleus pressure did not increase significantly for any loading direction at either the caudal or cephalic adjacent segment. Results suggest that single-level degeneration can increase the risk for injury at the adjacent levels.     

Is bone density associated with intervertebral disc pressure in healthy and degenerated discs?

The coupling of the intervertebral disc (IVD) and vertebra as a biomechanical unit suggests that changes in the distribution of pressure within the IVD (intradiscal pressure, IDP) as a result of disc degeneration can influence the distribution of bone density within the vertebra, and vice versa. The goal of this study was to assess the correspondence between IDP and bone density in the adjacent vertebrae, with emphasis on how this correspondence differs between healthy and degenerated IVDs. Bone density of the endplates and subchondral bone in regions adjacent to the anterior and posterior annulus fibrosus (aAF and pAF, respectively) and nucleus pulposus (NP) was measured via quantitative computed tomography (QCT) in 61 spine segments (T7-9, T9-11, T10-12; 71 ± 14 years). IDP was measured in the aAF, NP, and pAF regions in 26 of the spine segments (68 ± 16 years) while they were tested in flexed (5°) or erect postures. Disc degeneration was assessed by multiple grading schemes. No correlation was found between bone density and IDP in either posture (p > 0.104). Regional variations in IDP and, to a greater extent bone density, were found to change with advancing degeneration: both IDP (p = 0.045) and bone density (p = 0.024) decreased in the NP region relative to the aAF region. The finding of only a modest correspondence between degeneration-associated changes in IDP and bone density may arise from complexity in how IDP relates to mechanical force transmission through the endplate and from limitations of the available IVD grading schemes in estimating the mechanical behavior of the IVD.     

Preparation of Intact Bovine Tail Intervertebral Discs for Organ Culture

The intervertebral disc (IVD) is the joint of the spine connecting vertebra to vertebra. It functions to transmit loading of the spine and give flexibility to the spine. It composes of three compartments: the innermost nucleus pulposus (NP) encompassing by the annulus fibrosus (AF), and two cartilaginous endplates connecting the NP and AF to the vertebral body on both sides. Discogenic pain possibly caused by degenerative intervertebral disc disease (DDD) and disc herniations has been identified as a major problem in our modern society. To study possible mechanisms of IVD degeneration, in vitro organ culture systems with live disc cells are highly appealing. The in vitro culture of intact bovine coccygeal IVDs has advanced to a relevant model system, which allows the study of mechano-biological aspects in a well-controlled physiological and mechanical environment. Bovine tail IVDs can be obtained relatively easy in higher numbers and are very similar to the human lumbar IVDs with respect to cell density, cell population and dimensions. However, previous bovine caudal IVD harvesting techniques retaining cartilaginous endplates and bony endplates failed after 1-2 days of culture since the nutrition pathways were obviously blocked by clotted blood. IVDs are the biggest avascular organs, thus, the nutrients to the cells in the NP are solely dependent on diffusion via the capillary buds from the adjacent vertebral body. Presence of bone debris and clotted blood on the endplate surfaces can hinder nutrient diffusion into the center of the disc and compromise cell viability. Our group established a relatively quick protocol to "crack"-out the IVDs from the tail with a low risk for contamination. We are able to permeabilize the freshly-cut bony endplate surfaces by using a surgical jet lavage system, which removes the blood clots and cutting debris and very efficiently reopens the nutrition diffusion pathway to the center of the IVD. The presence of growth plates on both sides of the vertebral bone has to be avoided and to be removed prior to culture. In this video, we outline the crucial steps during preparation and demonstrate the key to a successful organ culture maintaining high cell viability for 14 days under free swelling culture. The culture time could be extended when appropriate mechanical environment can be maintained by using mechanical loading bioreactor. The technique demonstrated here can be extended to other animal species such as porcine, ovine and leporine caudal and lumbar IVD isolation.     

Load-displacement properties of lower cervical spine motion segments

The load-displacement behavior of 35 fresh adult cervical spine motion segments was measured in compression, shear, flexion, extension, lateral bending and axial torsion tests. Motion segments were tested both intact and with posterior elements removed. Applied forces ranged to 73.6 N in compression and to 39 N in shear, while applied moments ranged to 2.16 Nm. For each mode of loading, principal and coupled motions were measured and stiffnesses were calculated. The effect of disc degeneration on motion segment stiffnesses and the moments required for motion segment failure were also measured. In compression, the stiffnesses of the cervical motion segments were similar to those of thoracic and lumbar motion segments. In other modes of loading, cervical stiffnesses were considerably smaller than thoracic or lumbar stiffnesses. Removal of the posterior elements decreased cervical motion segment stiffnesses by as much as 50%. Degenerated cervical discs were less stiff in compression and stiffer in shear than less degenerated discs, but in bending or axial torsion, no statistically significant differences were evident. Bending moments causing failure were an order of magnitude lower than those for lumbar segments.     

Neural arch load-bearing in old and degenerated spines

We validate a technique for measuring neural arch load-bearing in cadaveric spines, and use it to test the hypothesis that such load-bearing rises to high levels in old and degenerated spines. Fifty-nine cadaveric lumbar motion segments, aged 19-92 yr, were subjected to compressive creep loading to reduce intervertebral disc water content and height to in vivo levels. The distribution of compressive "stress" within the disc was then measured by pulling a miniature pressure transducer, side-mounted in a 1.3mm-diameter needle, along its mid-sagittal diameter. During these measurements, the motion segment was subjected to a compressive load of 2 kN, and positioned in 2 degrees of extension to simulate erect standing. Measurements of compressive "stress" were integrated over disc area, and this force subtracted from the applied 2 kN to give the force resisted by the neural arch. An empirical calibration factor was applied to normalise results from each disc to values obtained under conditions when all of the compressive force could be assumed to pass through the disc. Disc degeneration was graded macroscopically on a scale of 1-4. Validation tests showed that calculated values of disc loading were proportional to actual applied load (r(2)>0.96) and predicted it with errors of 2-8%. Neural arch load-bearing in non-degenerated specimens was generally less than 20%, but averaged 49% for specimens aged over 70 yr. Multiple regression showed that neural arch load bearing (%)=14.4 x disc degeneration score+0.46 x age-35. These results indicate a substantial shift in vertebral load-bearing with increasing age and degeneration.     

Intervertebral disc creep behavior assessment through an open source finite element solver

Degenerative Disc Disease (DDD) is one of the largest health problems faced worldwide, based on lost working time and associated costs. By means of this motivation, this work aims to evaluate a biomimetic Finite Element (FE) model of the Intervertebral Disc (IVD). Recent studies have emphasized the importance of an accurate biomechanical modeling of the IVD, as it is a highly complex multiphasic medium. Poroelastic models of the disc are mostly implemented in commercial finite element packages with limited access to the algorithms. Therefore, a novel poroelastic formulation implemented on a home-developed open source FE solver is briefly addressed throughout this paper. The combination of this formulation with biphasic osmotic swelling behavior is also taken into account. Numerical simulations were devoted to the analysis of the non-degenerated human lumbar IVD time-dependent behavior. The results of the tests performed for creep assessment were inside the scope of the experimental data, with a remarkable improvement of the numerical accuracy when compared with previously published results obtained with ABAQUS^®. In brief, this in-development open-source FE solver was validated with literature experimental data and aims to be a valuable tool to study the IVD biomechanics and DDD mechanisms.     

Moderately degenerated lumbar motion segments: Are they truly unstable?

The two main load bearing tissues of the intervertebral disc are the nucleus pulposus and the annulus fibrosus. Both tissues are composed of the same basic components, but differ in their organization and relative amounts. With degeneration, the clear distinction between the two tissues disappears. The changes in biochemical content lead to changes in mechanical behaviour of the intervertebral disc. The aim of the current study was to investigate if well-documented moderate degeneration at the biochemical and fibre structure level leads to instability of the lumbar spine. By taking into account biochemical and ultrastructural changes to the extracellular matrix of degenerating discs, a set of constitutive material parameters were determined that described the individual tissue behaviour. These tissue biomechanical models were then used to simulate dynamic behaviour of the degenerated spinal motion segment, which showed instability in axial rotation, while a stabilizing effect in the other two principle bending directions. When a shear load was applied to the degenerated spinal motion segment, no sign of instability was found. This study found that reported changes to the nucleus pulposus and annulus fibrosus matrix during moderate degeneration lead to a more stable spinal motion segment and that such biomechanical considerations should be incorporated into the general pathophysiological understanding of disc degeneration and how its progress could affect low back pain and its treatments thereof.     

A biomechanical assessment to evaluate breed differences in normal porcine medial collateral ligaments

Little information is available on the role of genetic factors and heredity in normal ligament behaviour and their ability to heal. Assessing these factors is challenging because of the lack of suitable animal models. Therefore, the purpose of this study was to develop a porcine model in order to evaluate and compare the biomechanical differences of normal medial collateral ligaments (MCLs) between Yorkshire (YK) and red Duroc (RD) breeds. It was hypothesized that biomechanical differences would not exist between normal YK and RD MCLs. Comparisons between porcine and human MCL were also made. A biomechanical testing apparatus and protocol specific to pig MCL were developed. Ligaments were subjected to cyclic and static creep tests and then elongated to failure. Pig MCL morphology, geometry, and low- and high-load mechanical behaviour were assessed. The custom-designed apparatus and protocol were sufficiently sensitive to detect mechanical property differences between breeds as well as inter-leg differences. The results reveal that porcine MCL is comparable in both shape and size to human MCL and exhibits similar structural and material failure properties, thus making it a feasible model. Comparisons between RD and YK breeds revealed that age-matched RD pigs weigh more, have larger MCL cross-sectional area, and have lower MCL failure stress than YK pigs. The effect of weight may have influenced MCL geometrical and biomechanical properties, and consequently, the differences observed may be due to breed type and/or animal weight. In conclusion, the pig serves as a suitable large animal model for genetic-related connective tissue studies.     

Simulating the response of a standing operator to vibration stress by means of a biomechanical model

Under vibration stress the compressive forces transmitted in the joints of a standing operator are composed of nearly static and oscillating force parts. Because these forces can hardly be measured they were assessed by means of a biomechanical model. In the model, 27 rigid bodies with 103 degrees of freedom represent the segments of the human body. 106 force elements imitate the muscles of the trunk and the legs. At first, the model parameter were varied so that for the simulated sitting posture the model fits the seat-to-head transmissibility given in the literature and in ISO/CD 5982. For the standing posture, the transfer functions between the ground acceleration and the oscillating forces in the ankle, the knee, the hip, and the motion segment L3-L4 were computed. According to the moduli of these functions the forces in the ankles are higher than those in the knees or the hips and they nearly come up to the forces in the lumbar spine. Further the results of the simulation indicate that under equal vibration stress in the standing and the sitting posture the differences between the compressive forces in the lumbar spine are small.     

Effect of the intra-abdominal pressure and the center of segmental body mass on the lumbar spine mechanics - a computational parametric study

Determination of physiological loads in human lumbar spine is critical for understanding the mechanisms of lumbar diseases and for designing surgical treatments. Computational models have been used widely to estimate the physiological loads of the spine during simulated functional activities. However, various assumptions on physiological factors such as the intra-abdominal pressure (IAP), centers of mass (COMs) of the upper body and lumbar segments, and vertebral centers of rotation (CORs) have been made in modeling techniques. Systematic knowledge of how these assumptions will affect the predicted spinal biomechanics is important for improving the simulation accuracy. In this paper, we developed a 3D subject-specific numerical model of the lumbosacral spine including T12 and 90 muscles. The effects of the IAP magnitude and COMs locations on the COR of each motion segment and on the joint/muscle forces were investigated using a global convergence optimization procedure when the subject was in a weight bearing standing position. The data indicated that the line connecting the CORs showed a smaller curvature than the lordosis of the lumbar spine in standing posture when the IAP was 0?kPa and the COMs were 10?mm anterior to the geometric center of the T12 vertebra. Increasing the IAP from 0 kPa to 10 kPa shifted the location of CORs toward the posterior direction (from 1.4?±?8.9 mm anterior to intervertebral disc (IVD) centers to 40.5?±?3.1 mm posterior to the IVD centers) and reduced the average joint force (from 0.78?±?0.11 Body weight (BW) to 0.31?±?0.07 BW) and overall muscle force (from 349.3?±?57.7 N to 221.5?±?84.2 N). Anterior movement of the COMs from -30 mm to 70 mm relative to the geometric center of T12 vertebra caused an anterior shift of the CORs (from 25.1?±?8.3 mm posterior to IVD centers to 7.8?±?6.2 mm anterior to IVD centers) and increases of average joint forces (from 0.78?±?0.1 BW to 0.93?±?0.1 BW) and muscle force (from 348.9?±?47.7 N to 452.9?±?58.6 N). Therefore, it is important to consider the IAP and correct COMs in order to accurately simulate human spine biomechanics. The method and results of this study could be useful for designing prevention strategies of spinal injuries and recurrences, and for enhancing rehabilitation efficiency.