Internal load-sharing in the human passive lumbar spine: Review of in vitro and finite element model studies

Nucleus pulposus cell network modelling in the intervertebral disc

Intervertebral disc degeneration (IDD) results from an imbalance between anabolic and catabolic processes in the extracellular matrix (ECM). Due to complex biochemical interactions, a comprehensive understanding is needed. This study presents a regulatory network model (RNM) for nucleus pulposus cells (NPC), representing normal intervertebral disc (IVD) conditions. The RNM includes 33 proteins, and 153 interactions based on literature, incorporating key NPC regulatory mechanisms. A semi-quantitative approach calculates the basal steady state, accurately reflecting normal NPC activity. Model validation through published studies replicated pro-catabolic and pro-anabolic shifts, emphasizing the roles of transforming growth factor beta (TGF-β) and interleukin-1 receptor antagonist (IL-1Ra) in ECM regulation. This IVD RNM is a valuable tool for predicting IDD progression, offering insights into ECM degradation mechanisms and guiding experimental research on IVD health and degeneration.

Study on the biodynamic characteristics and internal vibration behaviors of a seated human body under biomechanical characteristics

To provide reference and theoretical guidance for establishing human body dynamics models and studying biomechanical vibration behavior, this study aimed to develop and verify a computational model of a three-dimensional seated human body with detailed anatomical structure under complex biomechanical characteristics to investigate dynamic characteristics and internal vibration behaviors of the human body. Fifty modes of a seated human body were extracted by modal method. The intervertebral disc and head motions under uniaxial white noise excitation (between 0 and 20 Hz at 1.0, 0.5 and 0.5 m/s2 r.m.s. for vertical, fore-aft and lateral direction, respectively) were computed by random response analysis method. It was found that there were many modes of the seated human body in the low-frequency range, and the modes that had a great impact on seated human vibration were mainly distributed below 13 Hz. The responses of different positions of the spine varied greatly under the fore-aft and lateral excitation, but the maximum stress was distributed in the lumbar under different excitations, which could explain why drivers were prone to lower back pain after prolonged driving. Moreover, there was a large vibration coupling between the vertical and fore-aft direction of an upright seated human body, while the vibration couplings between the lateral and other directions were very small. Overall, the study could provide new insights into not only the overall dynamic characteristics of the human body, but also the internal local motion and biomechanical characteristics under different excitations.

Biomechanics of annulus fibrosus: Elastic fiber simplification and degenerative impact on damage initiation and propagation

This study addresses three primary objectives related to lumbar intervertebral disc (IVD) biomechanics under ramping quasi-static loading conditions. First, we explore the conditions justifying the simplification of axisymmetric elastic fiber families into single fiber bundles through discretized strain energy functions. Simulations reveal that a concentration factor exceeding 10 allows for a consistent deviation below 10% between simplified and non-simplified responses. Second, we investigate the impact of elastic fibers on the physiological stiffness in IVDs, revealing minimal influence on biological motions but significant effects on degeneration. Lastly, we examine the initiation and progression of annulus fibrosus (AF) damage. Our findings confirm the validity of simplifying elastic fiber families and underscore the necessity of considering elastic fiber damage in biomechanical studies of AF tissues. Elastic fibers contribute to increased biaxial stretch stiffness, and their damage significantly affects the loading capacity of the inner AF. Additionally, degeneration significantly alters the susceptibility to damage in the AF, with specific regions exhibiting higher vulnerability. Damage tends to extend circumferentially and radially, emphasizing the regional variations in collagen and elastic fiber properties. This study offers useful insights for refining biomechanical models, paving the way for a more comprehensive understanding of IVD responses and potential clinical implications.

Influence of foot excitation and shin posture on the vibration behavior of the entire spine inside a seated human body

Due to ethical issues and simplification of traditional biomechanical models, experimental methods and traditional computer methods were difficult to quantify the effects of foot excitation and shin posture on vibration behavior of the entire spine inside a seated human body under vertical whole-body vibration. This study developed and verified different three-dimensional (3D) finite element (FE) models of seated human body with detailed anatomical structure under the biomechanical characteristics to predict vibration behavior of the entire spine inside a seated human body with different foot excitation (with and without vibration) and shin posture (vertical and tilt posture). Random response analysis was performed to study the transmissibility of the entire spine to seat under vertical white noise excitation between 0 and 20 Hz at 0.5 m/s2 r.m.s. The results showed that although the foot excitation could reduce the fore-aft transmissibility in the cervical spine (23% reduction), it could significantly increase that in the lumbar spine (52% increase), which resulted in complex alternating stresses at lumbar spine and made the lumbar spine more vulnerable to injury in long-term vibration environment. Moreover, the shin tilt posture made the maximum fore-aft transmissibility in the lumbar spine move to the upper lumbar spine. The study provided new insights into the influence of foot excitation and shin posture on the vibration behavior of the entire spine inside a seated human body. Foot excitation exposed the lumbar spine to complex alternating stresses and made it more vulnerable to injury in long-term whole body vibration.

Marked differences between continuous long-term and clinical snapshot examinations: is the current standard of back pain diagnostics outdated?

Current clinical examination of low back pain (LBP) patients primarily relies on static clinical examinations, which rarely represent the dynamic postures patients adopt during daily activities. To gain an overview on the dynamic kinematics-kinetics changes over a day, the lumbar back kinematics of asymptomatic individuals and LBP patients were measured over 24 h, and the passively resisted bending and torsional moments were estimated. 208 asymptomatic subjects (115 females) and 116 LBP patients (71 females) were analysed. Compared to static upright standing, the mean lumbar lordosis of asymptomatic subjects drops significantly by 21° during everyday life (p < 0.01). Maximum bending moments of 44.0-50.6 Nm were estimated at the L2-L3. LBP patients showed significantly lower (p < 0.01) lumbar flattening during daily life of about 16°. Maximum bending moments of 27-52 Nm were found at the L3-L4. The initial static upright lumbar lordosis was significantly lower in LBP population (by 6°) resulting in almost similar average lumbar shapes during daily activities in both groups. The torsional movements were with 2.2° greatest in L1-L2 independent of sex (p = 0.19) and LBP (p = 0.54) with moments of 6-16 Nm. The lumbar profile and associated internal moments during daily life differ substantially from those recorded during clinical examinations. LBP patients demonstrates significantly lower lordosis at the snapshot assessment and significantly lower movement variations and internal moments during daily life. Only the dynamic long-term assessments unravelled a less flexed posture in LBP population. Apparently, such a reduced dynamic flexed posture indicates a compensatory habit for pain relief.

Impact of osteoporosis and Cement-Augmented fusion on adjacent spinal levels Post-Fusion Surgery: Patient-Specific finite element analysis

Cement-augmentation is a technique commonly used during posterior lumbar instrumented fusion (PLIF) to reinforce compromised osteoporotic vertebral bone, minimize the risk of loosening screws, enhance stability, and improve overall surgical outcomes. In this study, we introduce a novel segmented vertebral body regional modeling approach to investigate the effects of osteoporosis and cement-augmented lumbar fusion on disc biomechanics at spinal levels adjacent to the fused vertebrae. Using our previously validated personalized-poroelastic-osteoligamentous FE model of the spine, fusion was simulated at L4-L5, and the biomechanics of adjacent levels were studied for 30 patients (non-osteoporotic patients (N = 15), osteoporotic patients (N = 15)). PLIF models, with and without cement-augmentation, were developed and compared after an 8 h-rest period (200 N), following a 16 h-cyclic compressive loading of 500-1000 N (40 and 20 min, respectively). Movement in different directions (flexion/ extension/ lateral bending/ axial rotation) was simulated using 10Nm moment before and after cyclic loading. The material mapping algorithm was validated by comparing the results of voxel-based and parametric models. The FE cement-augmented models, subject to daily activity loading, demonstrated significant differences in disc height loss and fluid loss as compared to non-cemented models. The calculated axial stress and fiber strain values were also significantly higher for these models. This work demonstrates that although osteoporosis does not significantly alter the time-dependent characteristics of adjacent IVDs post-surgery, cement-augmentation increases the risk of adjacent segment disease (ASD) incidence. A holistic understanding of the trade-offs and long-term complex interplay between structural reinforcement modalities, including cement augmentation, and altered biomechanics warrants further investigation.

Comprehensive modeling of annulus fibrosus: From biphasic refined characterization to damage accumulation under viscous loading

The annulus fibrosus (AF), a permeable, hydrated, and fiber-reinforced soft tissue, exhibits complex responses influenced by fluid pressure, osmotic pressure, and structural mechanics. Existing models struggle to comprehensively represent these intricate interactions and the heterogeneous solid responses within the AF. Additionally, the mechanisms driving differential damage accumulation between non-degenerative and degenerative intervertebral discs remain poorly understood. In this study, we introduce a biphasic-swelling damage model for the AF. We conceptually develop and rigorously validate this model through tissue-level tests employing various loading modes, consistently aligning model predictions with experimental data. Leveraging parametric geometric algorithms and custom Python scripts, we construct models simulating both non-degenerative and degenerative discs. Following calibration, we subject these models to viscous loading protocols. Our findings reveal the posterior AF's susceptibility to damage, contingent upon loading rate and water content. We elucidate the underlying mechanisms by examining the temporal evolution of fluid pressure, osmotic pressure, and the regionally dependent fiber network. This research presents a highly accurate model of the AF, providing valuable insights into disc damage. Future research endeavors should expand this model to incorporate ionic transport and diffusion, enabling a more profound exploration of intervertebral disc mechanobiology. This comprehensive model contributes to a better understanding of AF behavior and may inform therapeutic strategies for disc-related pathologies. STATEMENT OF SIGNIFICANCE: This research presents a comprehensive model of the annulus fibrosus (AF), a crucial component of the intervertebral disc that provides structural support and resists deformation. The study introduces a biphasic-swelling damage model for the AF and validates it through tissue-level tests. The model accounts for fluid pressure, osmotic pressure, and matrix mechanics, providing a more accurate representation of the AF's behavior. The study also investigates the differential damage accumulation between non-degenerative and degenerative discs, shedding light on the mechanisms driving disc degeneration. The findings have significant implications for medical treatments and interventions, as they highlight the posterior AF's susceptibility to damage. This research is of great interest to readers interested in biomechanics, tissue engineering, and medical treatments for disc degeneration.

On the identification of the ultra-structural organization of elastic fibers and their effects on the integrity of annulus fibrosus

Due to the complicated structure of the elastic fiber network in annulus fibrosus, existing in-silico studies either simplified or just overlooked its distribution pattern. Nonetheless, experimental and simulation results have proven that elastic fibers are of great importance to maintaining the structural integrity of annulus fibrosus and therefore to ensuring the load-bearing ability of intervertebral discs. Such needs call for a fine model. This work aims at developing a biphasic annulus fibrosus model by incorporating the accurate distribution pattern of collagen and elastic fibers. Both the structural parameters and intrinsic mechanical parameters were successfully identified using single lamella and inter-lamella microscopy anatomy and micromechanical testing data. The proposed model was then used to implement finite element simulations on various anterior and posterolateral multi-lamellae annulus fibrosus specimens. In general, simulation results agree well with available experimental and simulation data. On this basis, the effects of elastic fibers on the integrity of annulus fibrosus were further investigated. It was found that elastic fibers significantly influence the free swelling, radial stretching and circumferential shear performances of annulus fibrosus. Nonetheless, no significant effects were found for the circumferential stretching capability. The proposed biphasic model considers for the first time the distribution characteristics of elastic fibers at two scales, including both the principal orientations of all fiber families and the detailed distribution pattern within each family. Better understandings on the functions of collagen and elastic fibers can therefore be realized. To further enhance its prediction capability, the current model can be extended in the future by taking the fiber-matrix interaction as well as progressive damages into consideration.

A New Interlaminar/Interspinous and Facet-Joint Stabilization System in Lumbar Degenerative Disk Disease: 2 Years of Results

INTRODUCTION

In lumbar degenerative disk diseases (DDDs), we include a wide range of lumbar pathologies. Lumbar spinal stenosis with or without spondylolisthesis is a common cause of lower-limb pain in elderly patients. The surgical treatment of lumbar DDD consists of the decompression of the neural structures or the decompression and fusion of the involved motion segment. Unfortunately, rigid spinal implants followed by fusion cause increased stresses on the neighboring spinal segments, often leading to adjacent segment degeneration. The aim of this paper is to present a new system for interlaminar/interspinous and facet-joint stabilization and fusion.

MATERIAL AND METHODS

From March 2018 to June 2021, 175 patients with severe lumbar back and/or leg pain were operated on with this device after the failure of conservative treatment for a minimum of 6 months. For this study, we considered 75 available patients with a minimum follow-up time of 24 months. Patients rated their back pain and leg pain on a visual analog scale (VAS) after surgery and at the last follow-up; also, the postoperative consumption of analgesic drugs was investigated. Finally, patients were asked whether they would undergo this kind of surgery again or recommend it. The device is designed to promote a better and more efficient intervertebral fusion: Bone Ingrowth technology allows the system to achieve the best response from the bone tissue.

RESULTS

All procedures were performed without any complications. A reoperation was required in five cases (6.7%) because of implant failure. In one case, a spinous process and unilateral articular process fracture occurred. In four cases, a subcutaneous ISCHIA (one case) or FILICUDI (three cases) pullout was observed. Significant improvements following lumbar surgery were observed when evaluating the postoperative VAS values, analgesic drug consumption levels, and patients' satisfaction ratings. Evident fusion was seen in 58 of 75 patients (77.3%).

CONCLUSIONS

Our interspinous/interlaminar and facet-joint implant solution, associated with bone grafting, provided vertebral fusion in most stenotic patients with Grade I DS undergoing bilateral microdecompression (BMD) or bilateral decompression via the unilateral approach (BDUA). A higher number of patients and a longer follow-up will certainly be required to completely validate these new devices, but this minimally invasive surgery (MIS) is currently very encouraging and satisfactory.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

Discussion on the possibility of multi-layer intelligent technologies to achieve the best recover of musculoskeletal injuries: Smart materials, variable structures, and intelligent therapeutic planning

Although intelligent technologies has facilitated the development of precise orthopaedic, simple internal fixation, ligament reconstruction or arthroplasty can only relieve pain of patients in short-term. To achieve the best recover of musculoskeletal injuries, three bottlenecks must be broken through, which includes scientific path planning, bioactive implants and personalized surgical channels building. As scientific surgical path can be planned and built by through AI technology, 4D printing technology can make more bioactive implants be manufactured, and variable structures can establish personalized channels precisely, it is possible to achieve satisfied and effective musculoskeletal injury recovery with the progress of multi-layer intelligent technologies (MLIT).

Immuno-Modulatory Effects of Intervertebral Disc Cells

Low back pain is a highly prevalent, chronic, and costly medical condition predominantly triggered by intervertebral disc degeneration (IDD). IDD is often caused by structural and biochemical changes in intervertebral discs (IVD) that prompt a pathologic shift from an anabolic to catabolic state, affecting extracellular matrix (ECM) production, enzyme generation, cytokine and chemokine production, neurotrophic and angiogenic factor production. The IVD is an immune-privileged organ. However, during degeneration immune cells and inflammatory factors can infiltrate through defects in the cartilage endplate and annulus fibrosus fissures, further accelerating the catabolic environment. Remarkably, though, catabolic ECM disruption also occurs in the absence of immune cell infiltration, largely due to native disc cell production of catabolic enzymes and cytokines. An unbalanced metabolism could be induced by many different factors, including a harsh microenvironment, biomechanical cues, genetics, and infection. The complex, multifactorial nature of IDD brings the challenge of identifying key factors which initiate the degenerative cascade, eventually leading to back pain. These factors are often investigated through methods including animal models, 3D cell culture, bioreactors, and computational models. However, the crosstalk between the IVD, immune system, and shifted metabolism is frequently misconstrued, often with the assumption that the presence of cytokines and chemokines is synonymous to inflammation or an immune response, which is not true for the intact disc. Therefore, this review will tackle immunomodulatory and IVD cell roles in IDD, clarifying the differences between cellular involvements and implications for therapeutic development and assessing models used to explore inflammatory or catabolic IVD environments.

Biomechanical Properties of Paraspinal Muscles Influence Spinal Loading—A Musculoskeletal Simulation Study

Paraspinal muscles are vital to the functioning of the spine. Changes in muscle physiological cross-sectional area significantly affect spinal loading, but the importance of other muscle biomechanical properties remains unclear. This study explored the changes in spinal loading due to variation in five muscle biomechanical properties: passive stiffness, slack sarcomere length (SSL), in situ sarcomere length, specific tension, and pennation angle. An enhanced version of a musculoskeletal simulation model of the thoracolumbar spine with 210 muscle fascicles was used for this study and its predictions were validated for several tasks and multiple postures. Ranges of physiologically realistic values were selected for all five muscle parameters and their influence on L4-L5 intradiscal pressure (IDP) was investigated in standing and 36° flexion. We observed large changes in IDP due to changes in passive stiffness, SSL, in situ sarcomere length, and specific tension, often with interesting interplays between the parameters. For example, for upright standing, a change in stiffness value from one tenth to 10 times the baseline value increased the IDP only by 91% for the baseline model but by 945% when SSL was 0.4 μm shorter. Shorter SSL values and higher stiffnesses led to the largest increases in IDP. More changes were evident in flexion, as sarcomere lengths were longer in that posture and thus the passive curve is more influential. Our results highlight the importance of the muscle force-length curve and the parameters associated with it and motivate further experimental studies on in vivo measurement of those properties.

Computational modelling of the scoliotic spine: A literature review

Scoliosis is a deformity of the spine that in severe cases requires surgical treatment. There is still disagreement among clinicians as to what the aim of such treatment is as well as the optimal surgical technique. Numerical models can aid clinical decision-making by estimating the outcome of a given surgical intervention. This paper provided some background information on the modelling of the healthy spine and a review of the literature on scoliotic spine models, their validation, and their application. An overview of the methods and techniques used to construct scoliotic finite element and multibody models was given as well as the boundary conditions used in the simulations. The current limitations of the models were discussed as well as how such limitations are addressed in non-scoliotic spine models. Finally, future directions for the numerical modelling of scoliosis were addressed.

Numerical analysis of the mechanical behaviour of intact and implanted lumbar functional spinal units: Effects of loading and boundary conditions

The objective of this study was to develop an improved finite element (FE) model of a lumbar functional spinal unit (FSU) and to subsequently analyse the deviations in load transfer owing to implantation. The effects of loading and boundary conditions on load transfer in intact and implanted FSUs and its relationship with the potential risk of vertebral fracture were investigated. The FE models of L1-L5 and L3-L4 FSUs, intact and implanted, were developed using patient-specific CT-scan dataset and segmentation of cortical and cancellous bone regions. The effect of submodelling technique, as compared to artificial boundary conditions, on the elastic behaviour of lumbar spine was examined. Applied forces and moments, corresponding to physiologic movements, were used as loading conditions. Results indicated that the loading and boundary conditions considerably affect stress-strain distributions within a FSU. This study, based on an improved FE model of a vertebra, highlights the importance of using the submodelling technique to adequately evaluate the mechanical behaviour of a FSU. In the intact FSU, strains of 200-400 µε were observed in the cancellous bone of vertebral body and pedicles. High equivalent stresses of 10-25 MPa and 1-5 MPa were generated around the pars interarticularis for cortical and cancellous regions, respectively. Implantation caused reductions of 85%-92% in the range of motion for all movements. Insertion of the intervertebral cage resulted in major deviations in load transfer across a FSU for all movements. The cancellous bone around cage experienced pronounced increase in stresses of 10-15 MPa, which indicated potential risk of failure initiation in the vertebra.

A novel coupled musculoskeletal finite element model of the spine - Critical evaluation of trunk models in some tasks

Spine musculoskeletal (MS) models make simplifying assumptions on the intervertebral joint degrees-of-freedom (rotational and/or translational), representation (spherical or beam-like joints), and properties (linear or nonlinear). They also generally neglect the realistic structure of the joints with disc nuclei/annuli, facets, and ligaments. We aim to develop a novel MS model where trunk muscles are incorporated into a detailed finite element (FE) model of the ligamentous T12-S1 spine thus constructing a gold standard coupled MS-FE model. Model predictions are compared under some tasks with those of our earlier spherical joints, beam joints, and hybrid (uncoupled) MS-FE models. The coupled model predicted L4-L5 intradiscal pressures (R² ≅ 0.97, RMSE ≅ 0.27 MPa) and L1-S1 centers of rotation (CoRs) in agreement to in vivo data. Differences in model predictions grew at larger trunk flexion angles; at the peak (80°) flexion the coupled model predicted, compared to the hybrid model, much smaller global/local muscle forces (~38%), segmental (~44%) and disc (~22%) compression forces but larger segmental (~9%) and disc (~17%) shear loads, ligament forces at the lower lumbar levels (by up to 57%) and facet forces at all levels. The spherical/beam joints models predicted much greater muscle forces and segmental loads under larger flexion angles. Unlike the spherical joints model with fixed CoRs, the beam joints model predicted CoRs closer (RMSE = 2.3 mm in flexion tasks) to those of the coupled model. The coupled model offers a great potential for future studies towards improvement of surgical techniques, management of musculoskeletal injuries and subject-specific simulations.

Computational Challenges in Tissue Engineering for the Spine

This paper deals with a brief review of the recent developments in computational modelling applied to innovative treatments of spine diseases. Additionally, it provides a perspective on the research directions expected for the forthcoming years. The spine is composed of distinct and complex tissues that require specific modelling approaches. With the advent of additive manufacturing and increasing computational power, patient-specific treatments have moved from being a research trend to a reality in clinical practice, but there are many issues to be addressed before such approaches become universal. Here, it is identified that the major setback resides in validation of these computational techniques prior to approval by regulatory agencies. Nevertheless, there are very promising indicators in terms of optimised scaffold modelling for both disc arthroplasty and vertebroplasty, powered by a decisive contribution from imaging methods.

Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis

There is an increased interest in studying the biomechanics of the facet joints. For in silico studies, it is therefore important to understand the level of reliability of models for outputs of interest related to the facet joints. In this work, a systematic review of finite element models of multi-level spinal section with facet joints output of interest was performed. The review focused on the methodology used to model the facet joints and its associated validation. From the 110 papers analysed, 18 presented some validation of the facet joints outputs. Validation was done by comparing outputs to literature data, either computational or experimental values; with the major drawback that, when comparing to computational values, the baseline data was rarely validated. Analysis of the modelling methodology showed that there seems to be a compromise made between accuracy of the geometry and nonlinearity of the cartilage behaviour in compression. Most models either used a soft contact representation of the cartilage layer at the joint or included a cartilage layer which was linear elastic. Most concerning, soft contact models usually did not contain much information on the pressure-overclosure law. This review shows that to increase the reliability of in silico model of the spine for facet joints outputs, more needs to be done regarding the description of the methods used to model the facet joints, and the validation for specific outputs of interest needs to be more thorough, with recommendation to systematically share input and output data of validation studies.

Development of a multiscale model of the human lumbar spine for investigation of tissue loads in people with and without a transtibial amputation during sit-to-stand

Quantification of lumbar spine load transfer is important for understanding low back pain, especially among persons with a lower limb amputation. Computational modeling provides a helpful solution for obtaining estimates of in vivo loads. A multiscale model was constructed by combining musculoskeletal and finite element (FE) models of the lumbar spine to determine tissue loading during daily activities. Three-dimensional kinematic and ground reaction force data were collected from participants with ([Formula: see text]) and without ([Formula: see text]) a unilateral transtibial amputation (TTA) during 5 sit-to-stand trials. We estimated tissue-level load transfer from the multiscale model by controlling the FE model with intervertebral kinematics and muscle forces predicted by the musculoskeletal model. Annulus fibrosis stress, intradiscal pressure (IDP), and facet contact forces were calculated using the FE model. Differences in whole-body kinematics, muscle forces, and tissue-level loads were found between participant groups. Notably, participants with TTA had greater axial rotation toward their intact limb ([Formula: see text]), greater abdominal muscle activity ([Formula: see text]), and greater overall tissue loading throughout sit-to-stand ([Formula: see text]) compared to able-bodied participants. Both normalized (to upright standing) and absolute estimates of L4-L5 IDP were close to in vivo values reported in the literature. The multiscale model can be used to estimate the distribution of loads within different lumbar spine tissue structures and can be adapted for use with different activities, populations, and spinal geometries.

Estimation of Trunk Muscle Forces Using a Bio-Inspired Control Strategy Implemented in a Neuro-Osteo-Ligamentous Finite Element Model of the Lumbar Spine

Low back pain (LBP), the leading cause of disability worldwide, remains one of the most common and challenging problems in occupational musculoskeletal disorders. The effective assessment of LBP injury risk, and the design of appropriate treatment modalities and rehabilitation protocols, require accurate estimation of the mechanical spinal loads during different activities. This study aimed to: (1) develop a novel 2D beam-column finite element control-based model of the lumbar spine and compare its predictions for muscle forces and spinal loads to those resulting from a geometrically matched equilibrium-based model; (2) test, using the foregoing control-based finite element model, the validity of the follower load (FL) concept suggested in the geometrically matched model; and (3) investigate the effect of change in the magnitude of the external load on trunk muscle activation patterns. A simple 2D continuous beam-column model of the human lumbar spine, incorporating five pairs of Hill's muscle models, was developed in the frontal plane. Bio-inspired fuzzy neuro-controllers were used to maintain a laterally bent posture under five different external loading conditions. Muscle forces were assigned based on minimizing the kinematic error between target and actual postures, while imposing a penalty on muscular activation levels. As compared to the geometrically matched model, our control-based model predicted similar patterns for muscle forces, but at considerably lower values. Moreover, irrespective of the external loading conditions, a near (<3°) optimal FL on the spine was generated by the control-based predicted muscle forces. The variation of the muscle forces with the magnitude of the external load within the simulated range at the L1 level was found linear. This work presents a novel methodology, based on a bio-inspired control strategy, that can be used to estimate trunk muscle forces for various clinical and occupational applications toward shedding light on the ever-elusive LBP etiology.