Biomechanical analysis of the anterior cervical fusion

This paper presents a biomechanical analysis of the cervical C5-C6 functional spine unit before and after the anterior cervical discectomy and fusion. The aim of this work is to study the influence of the medical procedure and its instrumentation on range of motion and stress distribution. First, a three-dimensional finite element model of the lower cervical spine is obtained from computed tomography images using a pipeline of image processing, geometric modelling and mesh generation software. Then, a finite element study of parameters' influence on motion and a stress analysis at physiological and different post-operative scenarios were made for the basic movements of the cervical spine. It was confirmed that the results were very sensitive to intervertebral disc properties. The insertion of an anterior cervical plate influenced the stress distribution at the vertebral level as well as in the bone graft. Additionally, stress values in the graft decreased when it is used together with a cage.     

Computational analysis of bone remodeling during an anterior cervical fusion

The anterior cervical fusion is an established surgical procedure for spine stabilization after the removal of an intervertebral disc. However, it is not yet clear which bone graft represents the best choice and whether surgical devices can be efficient and beneficial for fusion. The aim of this work is to study the influence of the spine instrumentation on bone remodeling after a cervical spine surgery and, consequently, on the fusion process. A finite element model of the cervical spine was developed, having computed tomography images as input. Bone was modeled as a porous material characterized by the relative density at each point and the bone remodeling law was derived assuming that bone self-adapts in order to achieve the stiffest structure for the supported loads, with the total bone mass regulated by the metabolic cost of maintaining bone tissue. Apart from the analysis of healthy cervical spine, different surgical scenarios were tested: bone graft with or without a cage and the use of a stabilization plate system. Results showed that the anterior and posterior regions of the disc space are more important to stress transmission and that spinal devices reduce bone growth within bone grafts, being plate systems the most interfering elements. The material of the interbody cages plays a major role in fusion and, therefore, it should be carefully chosen.     

前路椎间盘减压融合与前路椎体次全切除减压融合治疗多节段颈椎病的疗效比较

目的:探索前路椎间盘减压融合与前路椎体次全切除减压治疗多节段颈椎病的疗效,为临床手术方式的选择提供依据。方法:收集我院骨科2008年6月到2014年6月收治的多节段颈椎病患者26例,按照患者手术方式分为研究组(13例)和对照组(13例),研究组给予前路椎间盘减压融合治疗,对照组给予前路椎体次全切除减压治疗,对比两组手术时间、术中出血量、术后住院时间,记录并分析两组术前和术后3月、6月、12个月JOA评分、颈椎总活动度、颈椎曲度、颈椎节段高度。结果:研究组手术时间、术中出血量低于对照组(P0.05);两组JOA评分术前、术后3月、6月、12个月逐渐升高(P0.05),术后12月组间差异有统计学意义(P0.05)。两组颈椎总活动度术前、术后3月、6月、12个月逐渐降低(P0.05),但是术后同时期组间差异无统计学意义(P0.05)。两组颈椎曲度与颈椎节段高度术后3月、6月、12个月差异有统计学意义(P0.05)。结论:前路椎间盘减压融合治疗多节段颈椎病较前路椎体次全切除减压治疗效果好,手术时间短、术中出血量少,并且颈椎曲度和节段高度恢复好。     

Bryan颈椎间盘假体置换术对脊髓型颈椎病患者疗效及颈椎生物力学的影响

目的:探讨Bryan颈椎间盘假体置换术对脊髓型颈椎病患者疗效及颈椎生物力学的影响。方法:选取2015年1月到2016年12月期间在我院接受治疗的脊髓型颈椎病患者48例,根据手术方式的不同将其分为植骨融合组(25例)和假体置换组(23例),其中植骨融合组采用颈椎前路减压植骨融合术进行治疗,假体置换组采用Bryan颈椎间盘假体置换术进行治疗。比较两组患者的日本骨科协会(JOA)颈椎评分、颈椎功能障碍指数(NDI)评分、视觉模拟疼痛量表(VAS)评分、颈椎生理曲度、颈椎活动度、手术节段活动度、上邻近节段活动度、下邻近节段活动度,并比较两组患者的并发症情况。结果:术后12个月假体置换组的NDI评分明显低于植骨融合组(P0.05);术后6个月、术后12个月植骨融合组的颈椎活动度低于假体置换组(P0.05);术后1个月、术后3个月、术后6个月、术后12个月假体置换组的手术节段活动度高于植骨融合组(P0.05);术后12个月植骨融合组的上邻近节段活动度、下邻近节段活动度高于假体置换组(P0.05);两组患者随访期间颈部轴性症状发生率比较差异有统计学意义(P0.05)。结论:与颈椎前路减压植骨融合术比较,Bryan颈椎间盘假体置换术对脊髓型颈椎病患者的远期疗效更佳,可更好的改善患者的颈椎生物力学,降低颈部轴性症状发生率,值得临床推广应用。     

Influence of cervical disc degeneration after posterior surgical techniques in combined flexion-extension--a nonlinear analytical study

Laminectomy and facetectomy are surgical techniques used for decompression of the cervical spinal stenosis. Recent in vitro and finite element studies have shown significant cervical spinal instability after performing these surgical techniques. However, the influence of degenerated cervical disk on the biomechanical responses of the cervical spine after these surgical techniques remains unknown. Therefore, a three-dimensional nonlinear finite element model of the human cervical spine (C2-C7) was created. Two types of disk degeneration grades were simulated. For each grade of disk degeneration, the intact as well as the two surgically altered models simulating C5 laminectomy with or without C5-C6 total facetectomies were exercised under flexion and extension. Intersegmental rotational motions, internal disk annulus, cancellous and cortical bone stresses were obtained and compared to the normal intact model. Results showed that the cervical rotational motion decreases with progressive disk degeneration. Decreases in the rotational motion due to disk degeneration were accompanied by higher cancellous and cortical bone stress. The surgically altered model showed significant increases in the rotational motions after laminectomies and facetectomies when compared to the intact model. However, the percentage increases in the rotational motions after various surgical techniques were reduced with progressive disk degeneration.     

Biomechanical investigation of post-operative C5 palsy due to ossification of the posterior longitudinal ligament in different types of cervical spinal alignment

Post-operative C5 palsies are among the most common complications seen after cervical surgery for ossification of the posterior longitudinal ligament (OPLL). Although C5 palsy is a well-known complication of cervical spine surgery, its pathogenesis is poorly understood and depends on many other factors. In this study, a finite element model of the cervical spine and spinal cord-nerve roots complex structures was developed. The changes in stress in the cord and nerve roots, posterior shift of the spinal cord, and displacement and elongation of the nerve roots after laminectomy for cervical OPLL were analyzed for three different cervical sagittal alignments (lordosis, straight, and kyphosis). The results suggest that high stress concentrated on the nerve roots after laminectomy could be the main cause of C5 palsy because ossification of ligaments increases spinal cord shifting and root displacement. The type of sagittal alignment had no influence on changes in cord stress after laminectomy, although cases of kyphosis with a high degree of occupying ratio resulted in greater increases in nerve root stress after laminectomy. Therefore, kyphosis with a high OPLL occupying ratio could be a risk factor for poor surgical outcomes or post-operative complications and should be carefully considered for surgical treatment.     

Predictive modelling of cervical disc implant wear

This study presents a chain of simulations aimed at estimating the wear in a cervical disc implant and providing insight into the in vivo biomechanical performance of the implant. The simulation chain can start with determining a representative maximum range of motion (ROM) of a person's head. The ROM is used as motion input to a kinematic simulation of the cervical spine containing a disc implant. The cervical spine geometry is obtained from computed tomography (CT) scans and converted to STL format using reverse engineering software. The time histories of the loads imposed by the adjacent vertebrae on the implant, as well as the vertebral relative rotations can be extracted from the kinematic simulation. Alternatively, force and motion profiles prescribed by wear test protocols (e.g. ISO 18192-1 and ASTM F2423-05) can be used. The force and motion profiles are applied as boundary conditions to a non-linear finite element model (FEM) of the implant to determine the time-varying contact stress and slip velocity distributions at the interface between the two halves of the implant. The stresses and slip velocities are used in a linear wear model to estimate the wear rate distribution at the FEM's nodal points where contact occurs. Reverse engineering software is used to triangulate the contact surface so that the total wear volume can be calculated. The simulation chain's predicted wear rate shows good agreement with in vitro results in the literature. The simulation chain is thereby demonstrated to be suitable for comparative pre-experimental studies of spinal implant designs.     

Parametric and subject-specific finite element modelling of the lower cervical spine. Influence of geometrical parameters on the motion patterns

Morphometrical and postural features of the cervical spine are supposed to significantly influence its biomechanical behaviour. However, the effects of these geometrical parameters are quite difficult to evaluate. An original numerical method is proposed in order to automatically generate parametric and subject-specific meshes of the lower cervical spine. Sixteen finite element models have been built from cadaver specimens using low dose biplanar X-rays. All the generated meshes fulfilled the quality criteria. A preliminary evaluation was performed on the C5–C6 functional units using a database of previous experimental tests. The principal and coupled motions were simulated. The responses of the numerical models were within the experimental standard deviation corridors in most cases. Rotation–moment relationships were then compared to assess the influence of geometry on the mechanical response. Geometry was found to play a significant role in the motion patterns.     

Artificial Cervical Vertebra and Intervertebral Complex Replacement through the Anterior Approach in Animal Model: A Biomechanical and In Vivo Evaluation of a Successful Goat Model

This was an in vitro and in vivo study to develop a novel artificial cervical vertebra and intervertebral complex (ACVC) joint in a goat model to provide a new method for treating degenerative disc disease in the cervical spine. The objectives of this study were to test the safety, validity, and effectiveness of ACVC by goat model and to provide preclinical data for a clinical trial in humans in future. We designed the ACVC based on the radiological and anatomical data on goat and human cervical spines, established an animal model by implanting the ACVC into goat cervical spines in vitro prior to in vivo implantation through the anterior approach, and evaluated clinical, radiological, biomechanical parameters after implantation. The X-ray radiological data revealed similarities between goat and human intervertebral angles at the levels of C2-3, C3-4, and C4-5, and between goat and human lordosis angles at the levels of C3-4 and C4-5. In the in vivo implantation, the goats successfully endured the entire experimental procedure and recovered well after the surgery. The radiological results showed that there was no dislocation of the ACVC and that the ACVC successfully restored the intervertebral disc height after the surgery. The biomechanical data showed that there was no significant difference in range of motion (ROM) or neural zone (NZ) between the control group and the ACVC group in flexion-extension and lateral bending before or after the fatigue test. The ROM and NZ of the ACVC group were greater than those of the control group for rotation. In conclusion, the goat provides an excellent animal model for the biomechanical study of the cervical spine. The ACVC is able to provide instant stability after surgery and to preserve normal motion in the cervical spine.     

Including surrounding tissue improves ultrasound-based 3D mechanical characterization of abdominal aortic aneurysms

ObjectivesIn this study the influence of surrounding tissues including the presence of the spine on wall stress analysis and mechanical characterization of abdominal aortic aneurysms using ultrasound imaging has been investigated.MethodsGeometries of 7 AAA patients and 11 healthy volunteers were acquired using 3-D ultrasound and converted to finite element based models. Model complexity of externally unsupported (aorta-only) models was complemented with inclusion of both soft tissue around the aorta and a spine support dorsal to the aorta. Computed 3-D motion of the aortic wall was verified by means of ultrasound speckle tracking. Resulting stress, strain, and estimated shear moduli were analyzed to quantify the effect of adding surrounding material supports.ResultsAn improved agreement was shown between the ultrasound measurements and the finite element tissue and spine models compared to the aorta-only models. Peak and 99-percentile Von Mises stress showed an overall decrease of 23–30%, while estimated shear modulus decreased with 12–20% after addition of the soft tissue. Shear strains in the aortic wall were higher in areas close to the spine compared to the anterior region.ConclusionsImproving model complexity with surrounding tissue and spine showed a homogenization of wall stresses, reduction in homogeneity of shear strain at the posterior side of the AAA, and a decrease in estimated aortic wall shear modulus. Future research will focus on the importance of a patient-specific spine geometry and location.     

四个节段以上伴椎间不稳颈椎病手术治疗的分析

目的:分析四个节段以上同时伴有椎间不稳的多节段颈椎病的手术入路及手术方法。方法:回顾性分析2008年10月-2012年12月收治的符合入选标准的颈椎病患者64例,其中A组33例,采用传统颈前路分节段开窗减压植骨内固定术;B组31例,采用改良颈后路单开门椎管减压轴侧植骨Arch钛板内固定术。采用日本矫形外科学会(JOA)评分标准和疼痛视觉模拟(VAS)评分标准对患者术后疗效进行评价,并对两组患者的术中出血量、手术时间、住院日数、颈椎活动度、颈椎曲度进行比较。结果:所有患者手术均顺利,A组患者的手术时间为(150.7±30.3)min,B组为(90.8±22.2)min,较A组明显缩短,差异具有统计学意义(P0.05)。A组患者的术中失血量为(320±50)m L,B组为(180±45)m L,较A组明显减少,差异具有统计学意义(P0.05)。此外,B组患者的住院时间显著短于短于A组(P0.05)。两组患者术后切口均I期愈合,出院时JOA评分及VAS评分均较术前明显改善(P0.05)。出院后6个月时,A组患者的Ishihara指数较术前显著改善(2.2±1.6),而B组无明显改善,A、B两组比较差异具有统计学意义(P0.05);A组患者的活动度丢失(4.2±3.3)°,B组活动度丢失(4.0±2.9)°,两组比较差异无统计学意义(t=0.26,P0.05)。结论:颈前路及颈后路手术方式治疗四个节段以上伴有椎间不稳颈椎病的患者均可获得理想的临床疗效,但颈后路手术方式的手术时间短,术中出血量少,住院周期短,安全性高,适应症广,是治疗四个节段以上颈椎病伴椎间不稳首选的手术方式。     

青年颈椎生理曲度变直全节段有限元建模与分析

目的建立青年颈椎生理曲度变直和正常的全节段有限元模型,进行对比与分析。方法采集青年颈椎曲度变直志愿者CT数据,利用生物力学有限元软件构建高质量颈椎全节段模型,然后采用基于离散微分属性的体网格变形技术,将变直模型映射生成曲度正常模型,在相同边界条件进行下进行对比分析。结果在前屈、后伸、侧弯、和旋转工况下,生理曲度变直模型的活动度范围比正常值要小,并且出现再分配。在应力分布方面,C3-C4,C4-C5之间的小关节、钩突关节,和C3-C4之间的椎间盘出现应力集中。结论通过使用新型建模软件和体网格变形技术,能方便构建颈椎生理曲度变直和正常模型,分析结果对青年颈椎病的临床诊治有指导意义。     

Biomechanical contributions of upper cervical ligamentous structures in Type II odontoid fractures

Fractures of the odontoid present frequently in spinal trauma, and Type II odontoid fractures, occurring at the junction of the odontoid process and C2 vertebrae, represent the bulk of all traumatic odontoid fractures. It is currently unclear what soft-tissue stabilizers contribute to upper cervical motion in the setting of a Type II odontoid fracture, and evaluation of how concomitant injury contributes to cervical stability may inform surgical decision-making as well as allow for the creation of future, accurate, biomechanical models of the upper cervical spine. The objective of the current study was to determine the contribution of soft-tissue stabilizers in the upper cervical spine following a Type II odontoid fracture. Eight cadaveric C0-C2 specimens were evaluated using a robotic testing system with motion tracking. The unilateral facet capsule (UFC) and anterior longitudinal ligament (ALL) were serially resected to determine their biomechanical role following odontoid fracture. Range of motion (ROM) and moment at the end of intact specimen replay were the primary outcomes. We determined that fracture of the odontoid significantly increases motion and decreases resistance to intact motion for flexion–extension (FE), axial rotation (AR), and lateral bending (LB). Injury to the UFC increased AR by 3.2° and FE by 3.2°. ALL resection did not significantly increase ROM or decrease end-point moment. The UFC was determined to contribute to 19% of intact flexion resistance and 24% of intact AR resistance. Overall, we determined that Type II fracture of the odontoid is a significant biomechanical destabilizer and that concurrent injury to the UFC further increases upper cervical ROM and decreases resistance to motion in a cadaveric model of traumatic Type II odontoid fractures.     

Effect of bone graft characteristics on the mechanical behavior of the lumbar spine

There is little information about the influence of bone graft size, position and elasticity on the mechanical behavior of the lumbar spine. Intersegmental motion, intradiscal pressure and stresses in the lumbar spine were calculated using a three-dimensional, nonlinear finite element model which included an internal spinal fixation device and a bone graft. Cross-sectional area, position, and elastic modulus of the graft were varied in this study. Bone grafts, especially very stiff ones, increase stresses on adjacent endplates. Though larger grafts lead to less contact pressure, it is difficult to judge the quality of different bone graft positions. In general, ventral flexion results in lower maximum contact pressure than lateral bending. There is always little intersegmental rotation in the bridged region compared with that of an intact spine.A larger graft with low stiffness should be favored from a mechanical point of view. Patients should avoid lateral bending of the upper body shortly after surgery.     

Stress amplifications in dental non-carious cervical lesions

This study aims to investigate the influence of the presence, shape and depth of NCCLs on the mechanical response of a maxillary second premolar subjected to functional and non-functional occlusal loadings using 3-D finite element (FE) analysis. A three-dimensional model of a maxillary second premolar and its supporting bone was constructed based on the contours of their cross-sections. From the sound model, cervical defects having either V- or U-shapes, as found clinically, were subtracted in three different depths. The models were loaded with 105 N to simulate normal chewing forces according to a functional occlusal loading (F1) vertically applied and two non-functional loadings (F2 and F3) obliquely oriented. Two alveolar bone crest heights were tested. Ansys™ FE software was used to compute stress distributions and maximum principal stress for each of the models. The presence of a lesion had no effect on the overall stress distribution of the system, but affected local stress concentrations. Non-functional loadings exhibited tensile stresses concentrating at the cervical areas and root surfaces, while the functional loading resulted in homogeneous stress distributions within the tooth. V-shaped lesions showed higher stress levels concentrated at the zenith of the lesion, whereas in U-shaped defect stresses concentrated over a wider area. As the lesions advanced in depth, the stress was amplified at their deepest part. A trend of stress amplification was observed with decreasing bone height. These results suggest a non-linear lesion progression with time, with the progression rate increasing with patient's age (deeper lesions and lower bone support).     

Implementation and validation of probabilistic models of the anterior longitudinal ligament and posterior longitudinal ligament of the cervical spine

The objective of this investigation was to develop probabilistic finite element (FE) models of the anterior longitudinal ligament (ALL) and posterior longitudinal ligament (PLL) of the cervical spine that incorporate the natural variability of biological specimens. In addition to the model development, a rigorous validation methodology was developed to quantify model performance. Experimental data for the geometry and dynamic properties of the ALL and PLL were used to create probabilistic FE models capable of predicting not only the mean dynamic relaxation response but also the observed experimental variation of that response. The probabilistic FE model uses a quasilinear viscoelastic material constitutive model to capture the time-dependent behaviour of the ligaments. The probabilistic analysis approach yields a statistical distribution for the model-predicted response at each time point rather than a single deterministic quantity (e.g. ligament force) and that response can be statistically compared to experimental data for validation. A quantitative metric that compares the cumulative distribution functions of the experimental data and model response is computed for both the ALL and PLL throughout the time histories and is used to quantify model performance.     

Influence of different mucosal resiliency and denture reline on stress distribution in peri-implant bone tissue during osseointegration. A three-dimensional finite element analysis

doi: 10.1111/j.1741‐2358.2011.00569.x Influence of different mucosal resiliency and denture reline on stress distribution in peri‐implant bone tissue during osseointegration. A three‐dimensional finite element analysis Objective: The aim of this study was to evaluate the influence of mucosal properties and relining material on the stress distribution in peri‐implant bone tissue during masticatory function with a conventional complete denture during the healing period through finite element analysis. Materials and Methods: Three‐dimensional models of a severely resorbed mandible with two recently placed implants in the anterior region were created and divided into the following situations: (i) conventional complete denture and (ii) relined denture with soft lining material. The mucosal tissue properties were divided into soft, resilient and hard. The models were exported to mechanical simulation software; two simulations were carried out with a load at the lower right canine (35 N) and the lower right first molar (50 N). Data were qualitatively evaluated using Maximum Principal Stress, in MPa, given by the software. Results: All models showed stress concentrations in the cortical bone corresponding to the cervical part of the implant. The mucosal properties influenced the stress in peri‐implant bone tissue showing a different performance according to the denture base material. The simulations with relined dentures showed lower values of stress concentration than conventional ones. Conclusions: It seems that the mucosal properties and denture reline have a high influence on the stress distribution in the peri‐implant bone during the healing period.     

Comparison of the anatomical morphology of cervical vertebrae between humans and macaques: related to a spinal cord injury model

Non-human primates are most suitable for generating cervical experimental models, and it is necessary to study the anatomy of the cervical spine in non-human primates when generating the models. The purpose of this study was to provide the anatomical parameters of the cervical spine and spinal cord in long-tailed macaques (Macaca fascicularis) as a basis for cervical spine-related experimental studies. Cervical spine specimens from 8 male adult subjects were scanned by micro-computed tomography, and an additional 10 live male subjects were scanned by magnetic resonance imaging. The measurements and parameters from them were compared to those of 12 male adult human subjects. Additionally, 10 live male subjects were scanned by magnetic resonance imaging, and the width and depth of the spinal cord and spinal canal and the thickness of the anterior and posterior cerebrospinal fluid were measured and compared to the relevant parameters of 10 male adult human subjects. The tendency of cervical parameters to change with segmental changes was similar between species. The vertebral body, spinal canal, and spinal cord were significantly flatter in the human subjects than in the long-tailed macaques. The cerebrospinal fluid space in the long-tailed macaques was smaller than that in the human subjects. The anatomical features of the cervical vertebrae of long-tailed macaques provide a reference for establishing a preclinical model of cervical spinal cord injury.     

Finite element analysis of moment-rotation relationships for human cervical spine

A comprehensive, geometrically accurate, nonlinear C0-C7 FE model of head and cervical spine based on the actual geometry of a human cadaver specimen was developed. The motions of each cervical vertebral level under pure moment loading of 1.0 Nm applied incrementally on the skull to simulate the movements of the head and cervical spine under flexion, tension, axial rotation and lateral bending with the inferior surface of the C7 vertebral body fully constrained were analysed. The predicted range of motion (ROM) for each motion segment were computed and compared with published experimental data. The model predicted the nonlinear moment-rotation relationship of human cervical spine. Under the same loading magnitude, the model predicted the largest rotation in extension, followed by flexion and axial rotation, and least ROM in lateral bending. The upper cervical spines are more flexible than the lower cervical levels. The motions of the two uppermost motion segments account for half (or even higher) of the whole cervical spine motion under rotational loadings. The differences in the ROMs among the lower cervical spines (C3-C7) were relatively small. The FE predicted segmental motions effectively reflect the behavior of human cervical spine and were in agreement with the experimental data. The C0-C7 FE model offers potentials for biomedical and injury studies.     

A three-dimensional finite element model of the human anterior cruciate ligament: a computational analysis with experimental validation

In this study, the force and stress distribution within the anteromedial (AM) and posterolateral (PL) bundles of the anterior cruciate ligament (ACL) in response to an anterior tibial load with the knee at full extension was calculated using a validated three-dimensional finite element model (FEM) of a human ACL. The interaction between the AM and PL bundles, as well as the contact and friction caused by the ACL wrapping around the bone during knee motion, were included in the model. The AM and PL bundles of the ACL were simulated as incompressible homogeneous and isotropic hyperelastic materials. The multiple-degrees-of-freedom (DOF) knee kinematics of a cadaveric knee were first obtained using a robotic/universal force-moment sensor testing system. These data were used as the boundary conditions for the FEM of the ACL to calculate the forces in the ACL. The calculated forces were compared to the in situ force in the ACL, determined experimentally, to validate the model. The validated FEM was then used to calculate the force and stress distribution within the ACL under an anterior tibial load at full extension. The AM and PL bundles shared the force, and the stress distribution was non-uniform within both bundles with the highest stress localized near the femoral insertion site. The contact and friction caused by the ACL wrapping around the bone during knee motion played the role of transferring the force from the ACL to the bone, and had a direct effect on the force and stress distribution of the ACL. This validated model will enable the analysis of force and stress distribution in the ACL in response to more complex loading conditions and has the potential to help design improved surgical procedures following ACL injuries.