Nonlinear numerical analysis of the structural response of the intervertebral disc to impact loading

3D-Printed PEEK/Silicon Nitride Scaffolds with a Triply Periodic Minimal Surface Structure for Spinal Fusion Implants

The issue of spine-related disorders is a global healthcare concern that requires effective solutions to restore normal spine functioning. Spinal fusion implants have become a standard approach for this purpose, making it crucial to develop biomaterials and structures that possess high osteogenic capacities and exhibit mechanical properties and dynamic responses similar to those of the host bone. This study focused on the fabrication of 3D-printed polyether ether ketone/silicon nitride (PEEK/SiN) scaffolds with a triply periodic minimal surface (TPMS) structure, which offers several advantages, such as a large surface area and uniform stress distribution under load. The mechanical properties and dynamic response of PEEK/SiN scaffolds with varying porosities were evaluated through mechanical testing and finite element analysis. The scaffold with 30% porosity exhibited a compressive strength (34.56 ± 1.91 MPa) and elastic modulus (734 ± 64 MPa) similar to those of trabecular bone. In addition, the scaffold demonstrated favorable damping properties. The biological data revealed that incorporating silicon nitride into the PEEK scaffold stimulated osteogenic differentiation. In light of these findings, it can be inferred that PEEK/SiN TPMS scaffolds exhibit significant potential for use in bone tissue engineering and represent a promising option as candidates for spinal fusion implants.

A comparative study on dynamic stiffness in typical finite element model and multi-body model of C6-C7 cervical spine segment

In contrast to numerous researches on static or quasi-static stiffness of cervical spine segments, very few investigations on their dynamic stiffness were published. Currently, scale factors and estimated coefficients were usually used in multi-body models for including viscoelastic properties and damping effects, meanwhile viscoelastic properties of some tissues were unavailable for establishing finite element models. Because dynamic stiffness of cervical spine segments in these models were difficult to validate because of lacking in experimental data, we tried to gain some insights on current modeling methods through studying dynamic stiffness differences between these models. A finite element model and a multi-body model of C6-C7 segment were developed through using available material data and typical modeling technologies. These two models were validated with quasi-static response data of the C6-C7 cervical spine segment. Dynamic stiffness differences were investigated through controlling motions of C6 vertebrae at different rates and then comparing their reaction forces or moments. Validation results showed that both the finite element model and the multi-body model could generate reasonable responses under quasi-static loads, but the finite element segment model exhibited more nonlinear characters. Dynamic response investigations indicated that dynamic stiffness of this finite element model might be underestimated because of the absence of dynamic stiffen effect and damping effects of annulus fibrous, while representation of these effects also need to be improved in current multi-body model. Copyright © 2015 John Wiley & Sons, Ltd.

NUMERICAL STUDY OF BIDIRECTIONAL GLENN WITH UNILATERAL PULMONARY ARTERY STENOSIS

Purpose: Hypoplastic left heart syndrome (HLHS) is a congenital heart disease and is usually associated with pulmonary artery stenosis. The superior vena cava-to-pulmonary artery (bidirectional Glenn) shunt is used primarily as a staging procedure to the total cava-to-pulmonary connection for single-ventricle complex. When HLHS coexists with pulmonary artery stenosis, the surgeons then face a multiple problem. This leads to high demand of optimized structure of Glenn surgery. The objective of this article is to investigate the influence of various anastomotic structures and the direction of superior vena cava (SVC) in Glenn on hemodynamics under pulse inflow conditions and try to find an optimal structure of SVC in Glenn surgery with unilateral pulmonary artery stenosis.

Method: First, 3D patient-specific models were constructed from medical images of a HLHS patient before any surgery by using the commercial software Mimics, and another software Free-form was used to deform the reconstructed models in the computer. Four 3D patient-specific Glenn models were constructed: model-1 (normal Glenn), model-2 (lean the SVC back to the stenotic pulmonary artery), model-3 (lean the SVC towards the stenotic pulmonary artery), model-4 (add patch at junction of the SVC toward stenosis at pulmonary artery). Second, a lumped parameter model (LPM) was established to predict boundary conditions for computational fluid dynamics (CFD). In addition, numerical simulations were conducted using CFD through the finite volume method. Finally, hemodynamic parameters were obtained and evaluated.

Results: It was showed that model-4 have relatively balanced vena cava blood perfusion into the left pulmonary artery (LPA) and right pulmonary artery (RPA), this may be due to less helical flow and the patch at junction of the SVC. Near stenosis of pulmonary artery, model-4 performed with the higher wall shear stress (WSS), which would benefit endothelial cell function and gene expression. In addition, results showed that model-4 performed with the lower oscillatory shear index (OSI) and wall shear stress gradient (WSSG), which would decrease the opportunity of vascular intimal hyperplasia.

Conclusion: It is benefited that surgeons adds patch at junction of the SVC towards stenosis at pulmonary artery. These results can impact the surgical design and planning of the Glenn surgery with unilateral pulmonary artery stenosis.

What have we learned from finite element model studies of lumbar intervertebral discs in the past four decades?

Finite element analysis is a powerful tool routinely used to study complex biological systems. For the last four decades, the lumbar intervertebral disc has been the focus of many such investigations. To understand the disc functional biomechanics, a precise knowledge of the disc mechanical, structural and biochemical environments at the microscopic and macroscopic levels is essential. In response to this need, finite element model studies have proven themselves as reliable and robust tools when combined with in vitro and in vivo measurements. This paper aims to review and discuss some salient findings of reported finite element simulations of lumbar intervertebral discs with special focus on their relevance and implications in disc functional biomechanics. Towards this goal, the earlier investigations are presented, discussed and summarized separately in three distinct groups of elastic, multi-phasic transient and transport model studies. The disc overall response as well as the relative role of its constituents are markedly influenced by loading rate, magnitude, combinations/preloads and posture. The nucleus fluid content and pressurizing capacity affect the disc compliance, annulus strains and failure sites/modes. Biodynamics of the disc is affected by not only the excitation characteristics but also preloads, existing mass and nucleus condition. The role of fluid pressurization and collagen fiber stiffening diminish with time during diurnal loading. The endplates permeability influences the time-dependent response of the disc in both loaded and unloaded recovery phases. The transport of solutes is substantially influenced by the disc size, tissue diffusivity and endplates permeability.