Effects of surgical joint destabilization on load sharing between ligamentous structures in the thoracic spine: a finite element investigation

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.

Finite element analysis of the cervical spine: dynamic characteristics and material property sensitivity study

The study aimed to investigate the dynamic characteristics of the cervical spine and determine the effect of the material properties of the cervical spinal components on it. A finite element model of the head-cervical spine was developed based on CT scan data, and the first six orders of modes (e.g. flexion-extension, lateral bending, and vertical, etc.) were verified by experimental and simulation studies. The material sensitivity study was conducted by varying elasticity modulus of cervical hard tissues (cortical bone, cancellous bone, endplates, and posterior elements) and soft tissues (intervertebral disc and ligaments). The results showed that increasing the elastic modulus of ligaments by 4 times increased the natural frequency by 77%, while increasing that of cancellous bone by 4 times only increased the natural frequency by 6%. In the axial mode, the cervical spine had not only axial deformation but also anterior-posterior deformation, with the largest deformation located at the intervertebral disc C6-C7. Decreasing the elastic modulus of a component in soft tissues by 80% increased modal displacement by up to 62%. The material properties of the intervertebral discs and ligaments had opposite effects on the modal displacement and deformation of the cervical spine. Low cervical discs were more susceptible to injury in a vertical vibration environment. Cervical spine dynamics were more sensitive to soft tissue material properties than to hard tissue material properties. Disc degeneration could reduce the range of vibratory motion of the cervical spine, thereby reducing the ability of the cervical spine to cushion head impacts.

A finite element method study of the effect of vibration on the dynamic biomechanical response of the lumbar spine

BACKGROUND

Studies focusing on lumbar spine biomechanics are very limited, and the mechanism of the effect of vibration on lumbar spine biodynamics is unclear. To provide guidance and reference for lumbar spine biodynamics research and vibration safety assessment, this study aims to investigate the effects of different vibrations on lumbar spine biodynamics.

METHODS

A validated finite element model of the lumbosacral spine was utilized. The model incorporated a 40 kg mass on the upper side and a 400 N follower preload. As a comparison, another model without a coupled mass was also employed. A sinusoidal acceleration with an amplitude of 1 m/s2 and a frequency of 5 Hz was applied to the upper and lower sides of the model respectively.

FINDINGS

When the coupled mass point is not introduced: in the case of upper-side excitation, the lumbar spine shows a significantly larger response in the x-direction than in the z-direction, while in the case of lower-side excitation, the lumbar spine experiences rigid body displacement in the z-direction without any movement, deformation, rotation, or stress changes in the x-direction. When the coupled mass point is introduced: both upper and lower-side excitations result in significant differences in z-directional displacement, with relatively small differences in vertebral rotation angle, disc deformation, and stress. Under upper excitation, low-frequency oscillations occur in the x-direction. In both types of excitations, the anterior-posterior deformation of the L2-L3 and L4-L5 intervertebral discs is greater than the vertical deformation. The peak (maximum) disc stress exceeds the average stress and stress amplitude across the entire disc. Regardless of the excitation type, the stress distribution within the disc at the moment of peak displacement remains nearly identical, with the maximum stress consistently localized on the anterior side of the L4-L5 disc.

INTERPRETATION

Accurately simulating lumbar spine biodynamics requires the inclusion of the upper body mass in the lumbosacral spine model. The physiological curvature of the lumbar spine could escalate the risk of lumbar spine vibration injuries. It is more instructive to apply local high stress in the disc as a lumbar spine vibration safety evaluation parameter.

The Impact of Spine Pathology on Posterior Ligamentous Complex Structure and Function

PURPOSE OF REVIEW

Spinal ligament is an important component of the spinal column in mitigating biomechanical stress. Particularly the posterior ligamentous complex, which is composed of the ligamentum flavum, interspinous, and supraspinous ligaments. However, research characterizing the biomechanics and role of ligament health in spinal pathology and clinical context are scarce. This article provides a comprehensive review of the implications of spinal pathology on the structure, function, and biomechanical properties of the posterior ligamentous complex.

RECENT FINDINGS

Current research characterizing biomechanical properties of the posterior ligamentous complex is primarily composed of cadaveric studies and finite element modeling, and more recently incorporating patient-specific anatomy into finite element models. The ultimate goal of current research is to understand the relative contributions of these ligamentous structures in healthy and pathological spine, and whether preserving ligaments may play an important role in spinal surgical techniques. At baseline, posterior ligamentous complex structures account for 30-40% of spinal stability, which is highly dependent on the intrinsic biomechanical properties of each ligament. Biomechanics vary widely with pathology and following rigid surgical fixation techniques and are generally maladaptive. Often secondary to morphological changes in the setting of spinal pathology, but morphological changes in ligament may also serve as a primary pathology. Biomechanical maladaptations of the spinal ligament adversely influence overall spinal column integrity and ultimately predispose to increased risk for surgical failure and poor clinical outcomes. Future research is needed, particularly in living subjects, to better characterize adaptations in ligaments that can provide targets for improved treatment of spinal pathology.

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.

A comparison of intervertebral ligament properties utilized in a thoracic spine functional unit through kinematic evaluation

Ligament properties in the literature are variable, yet scarce, but needed to calibrate computational models for spine clinical research applications. A comparison of ligament stiffness properties and their effect on the kinematic behavior of a thoracic functional spinal unit (FSU) is examined in this paper. Six unique ligament property sets were utilized within a volumetric T7-T8 finite element (FE) model developed using computer-aided design (CAD) spinal geometry. A 7.5 Nm moment was applied along three anatomical planes both with and without costovertebral (CV) joints present. Range of Motion (RoM) was assessed for each property set and compared to published experimental data. Intact and serial ligament removal procedures were implemented in accordance with experimental protocol. The variance in both kinematic behavior and comparability with experimental data among property sets emphasizes the role nonlinear characterization plays in determining proper kinematic behavior in spinal FE models. Additionally, a decrease in RoM variation among property sets was exhibited when the model setup incorporated the CV joint. With proper assessment of the source and size of each ligament, the material properties considered here could be expanded and justified for implementation into thoracic spine clinical studies.

Influence of configuration and anchor in ligamentous augmentation to prevent proximal junctional kyphosis: A finite element study

Ligament augmentation has been applied during spinal surgery to prevent proximal junctional kyphosis (PJK), but the configuration and distal anchor strategies are diverse and inconsistent. The biomechanics of different ligament augmentation strategies are, therefore, unclear. We aimed to create a finite element model of the spine for segments T6–S1. Model Intact was the native form, and Model IF was instrumented with a pedicle screw from segments T10 to S1. The remaining models were based on Model IF, with ligament augmentation configurations as common (CM), chained (CH), common and chained (CHM); and distal anchors to the spinous process (SP), crosslink (CL), and pedicle screw (PS), creating SP-CH, PS-CHM, PS-CH, PS-CM, CL-CHM, CL-CH, and CL-CM models. The range of motion (ROM) and maximum stress on the intervertebral disc (IVD), PS, and interspinous and supraspinous ligaments (ISL/SSL) was measured. In the PS-CH model, the ROM for segments T9–T10 was 73% (of Model Intact). In the CL-CHM, CL-CH, CL-CM, PS-CM, and PS-CHM models, the ROM was 8%, 17%, 7%, 13%, and 30%, respectively. The PS-CH method had the highest maximum stress on IVD and ISL/SSL, at 80% and 72%, respectively. The crosslink was more preferable as the distal anchor. In the uppermost instrumented vertebrae (UIV) + 1/UIV segment, the CM was the most effective configuration. The PS-CH model had the highest flexion load on the UIV + 1/UIV segment and the CL-CM model provided the greatest reduction. The CL-CM model should be verified in a clinical trial. The influence of configuration and anchor in ligament augmentation is important for the choice of surgical strategy and improvement of technique.

How Does the Rib Cage Affect the Biomechanical Properties of the Thoracic Spine? A Systematic Literature Review

The vast majority of previous experimental studies on the thoracic spine were performed without the entire rib cage, while significant contributive aspects regarding stability and motion behavior were shown in several other studies. The aim of this literature review was to pool and increase evidence on the effect of the rib cage on human thoracic spinal biomechanical characteristics by collating and interrelating previous experimental findings in order to support interpretations of in vitro and in silico studies disregarding the rib cage to create comparability and reproducibility for all studies including the rib cage and provide combined comparative data for future biomechanical studies on the thoracic spine. After a systematic literature search corresponding to PRISMA guidelines, eleven studies were included and quantitatively evaluated in this review. The combined data exhibited that the rib cage increases the thoracic spinal stability in all motion planes, primarily in axial rotation and predominantly in the upper thorax half, reducing thoracic spinal range of motion, neutral zone, and intradiscal pressure, while increasing thoracic spinal neutral and elastic zone stiffness, compression resistance, and, in a neutral position, the intradiscal pressure. In particular, the costosternal connection was found to be the primary stabilizer and an essential determinant for the kinematics of the overall thoracic spine, while the costotransverse and costovertebral joints predominantly reinforce the stability of the single thoracic spinal segments but do not alter thoracic spinal kinematics. Neutral zone and neutral zone stiffness were more affected by rib cage removal than the range of motion and elastic zone stiffness, thus also representing the essential parameters for destabilization of the thoracic spine. As a result, the rib cage and thoracic spine form a biomechanical entity that should not be separated. Therefore, usage of entire human non-degenerated thoracic spine and rib cage specimens together with pure moment application and sagittal curvature determination is recommended for future in vitro testing in order to ensure comparability, reproducibility, and quasi-physiological validity.

Preliminary evaluation of sagittal and transverse plane cross-sectional variations of the trunk during quiet and deep breathing by optical reflective motion analysis in patients with scoliosis

Trunk movements during quiet and deep breathing in untreated and in operated patients with scoliosis are not well defined. To evaluate sagittal and transverse plane cross-sectional variations of the trunk during quiet and deep breathing by optical reflective motion analysis (ORMA) in children with scoliosis. Twenty-one patients were divided into three groups: normal subjects (A; n = 6), subjects with untreated scoliosis >50° (B; n = 7) and operated patients (C; n = 8). Standing and sitting height, T1-T12 and L1-L5 length, arm span, chest perimeter, weight and BMI were recorded. Trunk movements of all patients, during quiet and deep breathing, were measured with a 10-camera 3D ORMA system (82 markers) with the subjects in a standard standing position. Groups were comparable with respect to age, sex, height, arm span and weight (P > 0.05). Significant differences were found in Cobb angle, chest perimeter and BMI (P < 0.05). Trunk sagittal and transverse plane movements during quiet and deep breathing decreased significantly in group B and group C when compared to group A (P < 0.05). Surgery does not completely eliminate sagittal and transverse plane kinematics, although the spine is rendered more rigid. This preliminary study involving a relatively limited number of patients outlines ORMA is a useful tool for analyzing sagittal and transverse plane motion abnormalities of the trunk. Trunk kinematics was altered during breathing in unoperated patients with scoliosis >50°. Operated subjects had sagittal plane trunk kinematics closer to normal subjects, although changes were less pronounced at the level of the convex side.

Morphometric analysis of the costal facet of the thoracic vertebrae

Various studies have examined morphometric features of the vertebrae to understand the functional aspects of the spine. Geometric analysis of vertebral zygapophyseal facets has also been related to functional and clinical aspects of the spine, but no quantitative investigation of the costotransverse joint facet is found in the literature. The costal facet geometry may partly determine the mechanical interaction between the rib cage and spine for trunk stabilization during functional tasks and during breathing. Therefore, the present study proposes a method for estimating the 3D geometric features of the costal facets of the first 10 thoracic vertebrae (Th1-Th10). Series of landmarks (95 ± 43) were placed on 258 costal facets from a sample of 14 asymptomatic individuals to determine their 3D location and orientation. The relative location of the costal facet was used to investigate symmetry and asymmetry components of the overall vertebrae shape variation among thoracic levels using 3D geometric morphometric methods. Results showed significant variation in sagittal orientation (inclination angle) between levels with a gradual cephalic orientation in the lower levels. No significant difference was observed on transverse orientation (declination angle). The shape of the costal facet was flatter at Th1 and from Th5 to Th10 and more concave from Th2 to Th4. An average difference of 7° between right and left facet orientation in both sagittal and transverse plane was demonstrated. Asymmetry of costal facet relative location was also detected and significantly influenced by the thoracic level. Nevertheless, location and orientation of the costal facets seem to be independent features of vertebrae morphology.

Ligaments of the Costovertebral Joints including Biomechanics, Innervations, and Clinical Applications: A Comprehensive Review with Application to Approaches to the Thoracic Spine

Few studies have examined the costovertebral joint and its ligaments in detail. Therefore, the following review was performed to better elucidate their anatomy, function and involvement in pathology. Standard search engines were used to find studies concerning the costovertebral joints and ligaments. These often-overlooked ligaments of the body serve important functions in maintaining appropriate alignment between the ribs and spine. With an increasing interest in minimally invasive approaches to the thoracic spine and an improved understanding of the function and innervation of these ligaments, surgeons and clinicians should have a good working knowledge of these structures.

Mechanism of right thoracic adolescent idiopathic scoliosis at risk for progression; a unifying pathway of development by normal growth and imbalance

Adolescent idiopathic scoliosis is regarded as a multifactorial disease and none of the many suggested causal etiologies have yet prevailed. I will suggest that adolescent idiopathic scoliosis has one common denominator, namely that initial curve development is mediated through one common normal physiological pathway of thoracic rotational instability. This is a consequence of gender specific natural growth of the passive structural components of thoracic spinal tissues for the adolescent female. This causes an unbalanced mechanical situation, which progresses if the paravertebral muscles cannot maintain spinal alignment. The alteration in the coronal plane with the lateral curve deformity is an uncoupling effect due to a culmination of a secondary, temporary sagittal plane thoracic flattening and of a primary, temporary transverse plane rotational instability for the adolescent female. Treatment of adolescent idiopathic scoliosis should address this physiological pathway and the overall treatment strategy is early intervention with strengthening of thoracic rotational stability for small curve adolescent idiopathic scoliosis.

Canine thoracic costovertebral and costotransverse joints: three case reports of dysfunction and manual therapy guidelines for assessment and treatment of these structures

The costovertebral and costotransverse joints receive little attention in research. However, pain associated with rib articulation dysfunction is reported to occur in human patients. The anatomic structures of the canine rib joints and thoracic spine are similar to those of humans. As such, it is proposed that extrapolation from human physical therapy practice could be used for the assessment and treatment of the canine patient with presumed rib joint pain. This article presents 3 case studies that demonstrate signs of rib dysfunction and successful treatment using primarily physical therapy manual techniques. General assessment and select treatment techniques are described.

Engineering-aided treatment of chest deformities to improve the process of breathing

In this paper, the application of new visualization technologies in correction of funnel chest is presented. Often, such defects cause problems with breathing, and therefore, to improve the cardiorespiratory efficiency, a correction is required. The virtual model of chest was formulated to determine the strains and stresses after correction of deformation with using the stabilizing plate. From biomechanical point of view, the knowledge of the load-affecting stabilizer is necessary to select optimal parameters of the plate. The force-affecting plate that was used in the correction of chest deformation by Nuss method was determined on the basis of finite element model of funnel chest. The main aim of this work was to select the optimal thickness of plate. Calculations were conducted for three alternative constraints.

An FE investigation simulating intra-operative corrective forces applied to correct scoliosis deformity

Background

Adolescent idiopathic scoliosis (AIS) is a deformity of the spine, which may require surgical correction by attaching a rod to the patient’s spine using screws implanted in the vertebral bodies. Surgeons achieve an intra-operative reduction in the deformity by applying compressive forces across the intervertebral disc spaces while they secure the rod to the vertebra. We were interested to understand how the deformity correction is influenced by increasing magnitudes of surgical corrective forces and what tissue level stresses are predicted at the vertebral endplates due to the surgical correction.

Methods

Patient-specific finite element models of the osseoligamentous spine and ribcage of eight AIS patients who underwent single rod anterior scoliosis surgery were created using pre-operative computed tomography (CT) scans. The surgically altered spine, including titanium rod and vertebral screws, was simulated. The models were analysed using data for intra-operatively measured compressive forces – three load profiles representing the mean and upper and lower standard deviation of this data were analysed. Data for the clinically observed deformity correction (Cobb angle) were compared with the model-predicted correction and the model results investigated to better understand the influence of increased compressive forces on the biomechanics of the instrumented joints.

Results

The predicted corrected Cobb angle for seven of the eight FE models were within the 5° clinical Cobb measurement variability for at least one of the force profiles. The largest portion of overall correction was predicted at or near the apical intervertebral disc for all load profiles. Model predictions for four of the eight patients showed endplate-to-endplate contact was occurring on adjacent endplates of one or more intervertebral disc spaces in the instrumented curve following the surgical loading steps.

Conclusion

This study demonstrated there is a direct relationship between intra-operative joint compressive forces and the degree of deformity correction achieved. The majority of the deformity correction will occur at or in adjacent spinal levels to the apex of the deformity. This study highlighted the importance of the intervertebral disc space anatomy in governing the coronal plane deformity correction and the limit of this correction will be when bone-to-bone contact of the opposing vertebral endplates occurs.

The rib head as a landmark in the anterolateral approach to the thoracic spine: a computed tomography-based morphometric study

OBJECT

The rib head is an important landmark in the anterolateral approach to the thoracic spine. Resection of the rib head is typically the first step in gaining access to the underlying pedicle and ultimately the spinal canal. The goal of this work is to quantify the relationship of the rib head to the spinal canal and adjacent aorta at each thoracic level using CT-based morphometric measurements.

METHODS

One hundred thoracic spine CT scans (obtained in 50 male and 50 female subjects) were evaluated in this study. The width and depth of each vertebra body were measured from T-1 to T-12. In addition, the distance of each rib head to the spinal canal was determined by drawing a line connecting the rib heads bilaterally and measuring the distance to this line from the most ventral aspect of the spinal canal. Finally, the distance of the left rib head to the thoracic aorta was measured at each thoracic level below the aortic arch.

RESULTS

The vertebral body depth progressively increased in a rostral to caudal direction. The vertebral body width was at its minimum at T-4 and progressively increased to T-12. The rib head extended beyond the spinal canal maximally at T-1. This distance incrementally decreased toward the caudal levels, with the tip of the rib head lying approximately even with the ventral canal at T-11 and T-12. The distance between the aorta and the left rib head increased in a rostral to caudal direction as well.

CONCLUSIONS

The rib head is an important landmark in the anterolateral approach to the thoracic spine. At more cephalad levels, a larger portion of rib head requires resection to gain access to the spinal canal. At more caudad levels, there is a safer working distance between the rib head and aorta.