Computational Challenges in Tissue Engineering: Comparison
Please note this is a comparison between Version 2 by Vivi Li and Version 1 by André Castro.

This entry 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.

  • spine
  • intervertebral disc
  • vertebrae
  • tissue engineering
  • finite element modelling
  • biomechanics

1. Introduction

Lumbar spine conditions related to degeneration are a serious source of disability worldwide, as they affect a large portion of the population. This results in high care costs for therapy and treatment, especially in Western societies. Currently, this problem is one of the major reasons for work absenteeism and productivity decrease, denoting itself to be an alert for society [1,2][1][2]. Recent statistics from the current decade point out that 70% to 85% of the Western population has, or is expected to have at some point, a record of low back pain (LBP), the global prevalence of chronic LBP being estimated between 7% and 9% [3,4][3][4].

The high frequency of these conditions can be explained by the ageing of the population and by the prominence of sedentary lifestyles, but its actual roots are far from being fully understood [5]. Intervertebral disc (IVD) degeneration or ligament hypertrophy, which in the most severe cases can greatly impact people’s quality of life, are among the known risk factors [2,6][2][6].

The IVD works with vertebral bodies (VB), ligaments and muscles to provide flexibility, shock absorbance and motion to the spine [7]. It has three components: (i) nucleus pulposus (NP), a highly hydrated biphasic material with dispersed collagen fibres; (ii) annulus fibrosus (AF), composed of successive lamellae of highly anisotropic collagen fibres, responsible for containing the NP; and (iii) cartilaginous endplate (CEP), a thin layer of hyaline cartilage, responsible for most of the nutrient exchange [8,9,10][8][9][10]. The complexity of each component and their connection is a challenge for both healthiness and degeneration studies. Each VB–IVD–VB complex is denoted as a functional spinal unit (FSU), the basic unit of the spine [11,12][11][12].

Furthermore, ligaments and muscles play an important role in stability as well as in motion proprioception. There are seven major ligaments constraining the spine components and limiting the range of motion (ROM) in all directions. Due to their elasticity, they passively help return the spine back to the neutral position [13,14][13][14]. These ligaments of the spine are the ligamentum flavum (LF), attached between the lamina of each VB; the anterior longitudinal ligament (ALL); the posterior longitudinal ligament (PLL), which works with the ALL along the consecutive VB; the intertransverse ligament (ITL), which connects the transverse processes; the interspinous ligament (ISL), which connects the opposing edges of the spinous processes; the supraspinous ligament (SSL), which connects the peaks of the spinous processes; and the capsular ligament (CL), which connects the circumferences of the joining articular facet joints [15,16][15][16]. The first three ligaments allow for flexion and extension of the spine, while keeping the bony parts in their correct alignment. The remaining four ligaments connect all the posterior elements of the VB [13,17][13][17].

Degenerative spine diseases (DSD) associated with LBP include degenerative disc disease (DDD), spondylolisthesis (wherein a vertebra suffers a displacement relative to the adjacent one) and stenosis (narrowing of the spinal canal) [18,19][18][19]. A lack of widely accepted explanatory models limits the understanding of such diseases, retarding the development of effective therapies [20,21][20][21]. The most common treatments of DSD are either conservative (no surgery) or surgical (invasive). However, these methods may act only to relieve symptoms; further complications frequently arise from massively invasive procedures, such as interbody fusion, disc arthroplasty or vertebroplasty [22,23,24,25,26][22][23][24][25][26]. Therefore, minimally invasive surgical approaches, based on tissue engineering (TE) methods, should be applied, whenever possible, to recover healthy functionality of the spine [27,28,29][27][28][29]. Additionally, to better diagnose and treat spinal instability, as well as to create models to support clinical decisions and surgical procedures, it is of utmost relevance to understand the kinematics and constitutive mechanics of the spine. Since in vitro and in vivo studies have some inherent limitations, computational models are useful tools in the investigation of these issues [30,31][30][31].

2. Computational Biomechanics

Since the first application of the finite element (FE) method in biomechanics by Brekelmans et al. in 1972 [32] and the first reported spine model developed by Belytschko et al. in 1974 [33], the number of studies on the human spine without resorting to animal and cadaveric experiments has increased significantly [34]. This application has countless benefits to health and biomedical research. As the human spine is a complex system, the intrinsic easiness with FE modelling in changing the geometry, constitutive laws, material properties and boundary conditions can lead to many important advances regarding spinal health, diseases, degeneration, trauma, surgery procedures and spinal instrumentation. Other computational methods, such as multibody systems modelling (usually dedicated to movement studies), can also be linked with FE modelling to establish multiscale approaches to spine biomechanics, providing macromechanical movement data (e.g., load, pressure or velocity) as simulation inputs [35,36][35][36].

In fact, FE modelling has already led to relevant conclusions and supported many pre-clinical studies, namely in what concerns to IVD behaviour or implant development [8,37,38,39,40][8][37][38][39][40]. Additionally, it has also motivated the development of new computational tools (imaging reconstruction methods, FE solvers or even kinematics simulators) that have later led to wider achievements in biomechanics and tissue engineering [41[41][42][43][44][45][46][47].

2.1. Macromechanical Level

From the works of Breau et al. in 1991 [48] or Chen et al. in 2001 [49], in which FE models of the lumbar spine were built from computed tomography (CT) scans, a more accurate perspective of the spine components was potentiated. The first group developed an L1–S1 model, including the seven relevant spine ligaments (LF, ALL, PLL, ITL, ISL, SSL and CL) plus the iliolumbar and the fascia,42 simulated by uniaxial elements. The AF was modelled as a ground substance with collagenous fibres embedded,43 while the NP was represented as an inviscid fluid. The fact joints were designed as a general moving contact problem [48]. The second group included a similar FSU modelling approach,44 with the seven ligaments represented by cable elements,45 along with the fibres of the AF. The facet articulations were modelled as 3D contact elements. This model was employed in simulating a lumbar fusion intervention,46 in which the IVD of the affected level was replaced by an interbody bone graft. It was possible to study the ROM and the stress distribution on the adjacent discs after surgery,47] which would not be directly possible without these numerical methods [49].

2.1. Macromechanical Level

From the works of Breau et al. in 1991 [48] or Chen et al. in 2001 [49], in which FE models of the lumbar spine were built from computed tomography (CT) scans, a more accurate perspective of the spine components was potentiated. The first group developed an L1–S1 model, including the seven relevant spine ligaments (LF, ALL, PLL, ITL, ISL, SSL and CL) plus the iliolumbar and the fascia, simulated by uniaxial elements. The AF was modelled as a ground substance with collagenous fibres embedded, while the NP was represented as an inviscid fluid. The fact joints were designed as a general moving contact problem [48]. The second group included a similar FSU modelling approach, with the seven ligaments represented by cable elements, along with the fibres of the AF. The facet articulations were modelled as 3D contact elements. This model was employed in simulating a lumbar fusion intervention, in which the IVD of the affected level was replaced by an interbody bone graft. It was possible to study the ROM and the stress distribution on the adjacent discs after surgery, which would not be directly possible without these numerical methods [49]. Figure 1

 shows a common patient-specific pathway for 3D model reconstruction from imaging examinations (CT scans), considering the design of the IVD on a given computer aided design (CAD) software and then the inclusion of the ligaments within the FE modelling software.

Figure 1.

 Example of a common pathway for patient-specific finite element (FE) modelling of the spine, going from medical imaging to the final FE model.

More recent studies were conducted to investigate new possible designs for spine implants. Choi et al. in 2016 [50] studied the biomechanical effects in surgical and adjacent FSUs of a newly interspinous process compressor and compared it with the traditional pedicle screw fixation. The different FE models of the lumbar spine were evaluated in terms of the ROM, facet contact forces and stress in the adjacent IVD. Guo and Yin in 2019 [51] proposed a new device interspinous stabilisation process using topology optimisation, which is a very important numerical technique to achieve better designs and more adequate material properties for biomedical devices [52,53]. A new device with reduced volume was fabricated from their FE modelling strategy, and the results showed that the innovative device provided stability in all motion directions at the affected segment, while decreasing the stress of the implant structure [51]. Lu and Lu in 2019 [54] built an L3–L5 FE model comprising the NP modelled as a linear elastic fluid, the AF represented by a ground substance with layers of diagonal fibres of varying strength and the seven major ligaments. Different interventions were simulated at the L4–L5 level, leading to the conclusion that lumbar interbody fusion techniques permitted greater stability, having shown a smaller ROM and stress peaks on the posterior instrumentation than posterolateral fusion. Demir et al. in 2020 [40] compared three different hybrid stabilisation systems (combination of dynamic and static fixation), starting from an L1–L5 FE model from CT scans. Hybrid stabilisation has been a focus of discussion as regards lumbar interbody fusion, and this study was able to demonstrate that the movement supported by the dynamic fixation in the adjacent level could reduce fracture risk and the overall intradiscal pressure (IDP). Such insight into tissue behaviour and mechanical responses would also not be possible with other in vitro or ex vivo techniques [8,54]. Zhou et al. in 2020 [55] combined in vivo imaging with FE simulation to evaluate the dynamics of adjacent-level degeneration after interbody fusion, creating a framework for pre- and post-operative assessment of the patients. This study is important to show the way for simulation integration in clinical practice.

More recent studies were conducted to investigate new possible designs for spine implants. Choi et al. in 2016 [50] studied the biomechanical effects in surgical and adjacent FSUs of a newly interspinous process compressor and compared it with the traditional pedicle screw fixation. The different FE models of the lumbar spine were evaluated in terms of the ROM, facet contact forces and stress in the adjacent IVD. Guo and Yin in 2019 [51] proposed a new device interspinous stabilisation process using topology optimisation, which is a very important numerical technique to achieve better designs and more adequate material properties for biomedical devices [52][53]. A new device with reduced volume was fabricated from their FE modelling strategy, and the results showed that the innovative device provided stability in all motion directions at the affected segment, while decreasing the stress of the implant structure [51]. Lu and Lu in 2019 [54] built an L3–L5 FE model comprising the NP modelled as a linear elastic fluid, the AF represented by a ground substance with layers of diagonal fibres of varying strength and the seven major ligaments. Different interventions were simulated at the L4–L5 level, leading to the conclusion that lumbar interbody fusion techniques permitted greater stability, having shown a smaller ROM and stress peaks on the posterior instrumentation than posterolateral fusion. Demir et al. in 2020 [40] compared three different hybrid stabilisation systems (combination of dynamic and static fixation), starting from an L1–L5 FE model from CT scans. Hybrid stabilisation has been a focus of discussion as regards lumbar interbody fusion, and this study was able to demonstrate that the movement supported by the dynamic fixation in the adjacent level could reduce fracture risk and the overall intradiscal pressure (IDP). Such insight into tissue behaviour and mechanical responses would also not be possible with other in vitro or ex vivo techniques [8][54]. Zhou et al. in 2020 [55] combined in vivo imaging with FE simulation to evaluate the dynamics of adjacent-level degeneration after interbody fusion, creating a framework for pre- and post-operative assessment of the patients. This study is important to show the way for simulation integration in clinical practice.

2.2. Tissue Level

The development of implants or new surgical approaches is definitely within the top priorities for human spine studies. However, there are still questions to be answered as regards the biomechanics and mechanobiology of native spine components. Experimental characterisation and numerical modelling of the IVD have been a challenge throughout the years, particularly as regards its multiphasic charge-driven behaviour and the anisotropic behaviour of its fibres [9,56,57,58,59,60]. CT scans are one of the sources of imaging data, given their high resolution, although magnetic resonance imaging (MRI) is more accurate for IVD degeneration detection. X-ray scans are more limited due to their 2D nature, but they are also employed as they permit less radiation exposure than CT equipment when considering in vivo evaluations [44,61,62,63].

The development of implants or new surgical approaches is definitely within the top priorities for human spine studies. However, there are still questions to be answered as regards the biomechanics and mechanobiology of native spine components. Experimental characterisation and numerical modelling of the IVD have been a challenge throughout the years, particularly as regards its multiphasic charge-driven behaviour and the anisotropic behaviour of its fibres [9][56][57][58][59][60]. CT scans are one of the sources of imaging data, given their high resolution, although magnetic resonance imaging (MRI) is more accurate for IVD degeneration detection. X-ray scans are more limited due to their 2D nature, but they are also employed as they permit less radiation exposure than CT equipment when considering in vivo evaluations [44][61][62][63].

Additionally, known DDD triggers such as osmotic pressure variation, CEP calcification, NP internal cracks or AF bulging are directly associated with intricate constitutive properties [27,68]. Consequently, the characterisation of FSU tissues (not only the IVD but also the VB, muscles and ligaments) before and after initiation of the degenerative process is the most accurate strategy to tackle the alterations associated with DDD and find potential treatments [44,69,70,71], i.e., it is not realistic to study a DDD symptom such as IVD height decrease or dehydration without analysing each adjacent component.

Additionally, known DDD triggers such as osmotic pressure variation, CEP calcification, NP internal cracks or AF bulging are directly associated with intricate constitutive properties [27][64]. Consequently, the characterisation of FSU tissues (not only the IVD but also the VB, muscles and ligaments) before and after initiation of the degenerative process is the most accurate strategy to tackle the alterations associated with DDD and find potential treatments [44][65][66][67], i.e., it is not realistic to study a DDD symptom such as IVD height decrease or dehydration without analysing each adjacent component.

The work of Schroeder et al. in 2010 [64] was one of the first and most complete approaches to the real mechanical behaviour of the IVD by implementing a 3D osmo-poro-visco-hyper-elastic FE model. In this study, the behaviour of the IVD was simulated as a function of NP and AF tissue biochemical composition, organisation and constituent properties, bridging chemistry, biology and physics. Fluid flow or the AF bulging response was consistent with experimental studies, demonstrating that a comprehensive approach is necessary to understand the overall IVD response. The idea of studying the influence of each specific parameter in order to distinguish between native and degenerated conditions was then followed in many other works [8,72,73,74].

The work of Schroeder et al. in 2010 [68] was one of the first and most complete approaches to the real mechanical behaviour of the IVD by implementing a 3D osmo-poro-visco-hyper-elastic FE model. In this study, the behaviour of the IVD was simulated as a function of NP and AF tissue biochemical composition, organisation and constituent properties, bridging chemistry, biology and physics. Fluid flow or the AF bulging response was consistent with experimental studies, demonstrating that a comprehensive approach is necessary to understand the overall IVD response. The idea of studying the influence of each specific parameter in order to distinguish between native and degenerated conditions was then followed in many other works [8][69][70][71].

The most recent generation of IVD/FSU studies is focused on patient-specific modelling, i.e., to associate what is known about IVD or VB biomechanics with the specific characteristics of each patient [75]. However, this strategy will only be effective with large databases of cases, provided that such databases contain complete medical data: top-quality medical imaging examinations are not enough, as the clinical background of the patient is also necessary to achieve a meaningful patient-specific FE model of the spine [76]. Fattor et al. in 2016 [77] employed CT scans to tailor total disc replacement devices and hence ensure correct geometrical fitting. The study considered not only pre-operative examinations but also follow-up post-operative data to establish a complete evaluation of the expected device performance and intraoperative alignment.

The most recent generation of IVD/FSU studies is focused on patient-specific modelling, i.e., to associate what is known about IVD or VB biomechanics with the specific characteristics of each patient [72]. However, this strategy will only be effective with large databases of cases, provided that such databases contain complete medical data: top-quality medical imaging examinations are not enough, as the clinical background o[73] employed CT scans to tailor total disc replacement devices and hence ensure correct geometrical fitting. The study considered not only pre-operative examinations but also follow-up post-operative data to establish a complete evaluation of the expected device performance and intraoperative alignment.

The works of Lavecchia et al. in 2018 [47] or Bojairami et al. in 2020 [78] show dedicated tools for the generation of geometrically accurate models for patient-specific FE studies, based on either age-height correlation for a given population or imaging examinations. This eliminates the need for a full 3D reconstruction from imaging, while ensuring good-quality FE meshes. Still, it must be highlighted that geometry is only one of the vectors to be considered in patient-specific spine modelling. The work of van Rijsbergen et al. in 2018 [79] is one of the most complete approaches to understanding how to go from imaging to actual clinical decision support. The developed framework to predict treatment outcomes is supported in high-level numerical modelling of FSU tissues, from mechanoregulation of the IVD environment to adjacent bone remodelling.

The works of Lavecchia et al. in 2018 [47] or Bojairami et al. in 2020 [74] show dedicated tools for the generation of geometrically accurate models for patient-specific FE studies, based on either age-height correlation for a given population or imaging examinations. This eliminates the need for a full 3D reconstruction from imaging, while ensuring good-quality FE meshes. Still, it must be highlighted that geometry is only one of the vectors to be considered in patient-specific spine modelling. The work of van Rijsbergen et al. in 2018 [75] is one of the most complete approaches to understanding how to go from imaging to actual clinical decision support. The developed framework to predict treatment outcomes is supported in high-level numerical modelling of FSU tissues, from mechanoregulation of the IVD environment to adjacent bone remodelling.

A different trend of FE-based studies is to emphasise on the modelling of specific FSU tissues, namely the cortical/trabecular layers of the VB [62,66,80,81], the layered organisation of the AF [58,65,82,83] or the fluid exchange and osmotic swelling properties of the NP [84,85,86,87]. These (and other) studies are contributing to identifying the complexity of FSU tissues and to understanding their interaction, ranging from vertebroplasty support to the importance of the AF criss-cross angles in the herniation process, despite the frequent simplifications on the modelling of surrounding tissues (i.e., the FSU components that are not the focus of such study).

A different trend of FE-based studies is to emphasise on the modelling of specific FSU tissues, namely the cortical/trabecular layers of the VB [62][76][77][78], the layered organisation of the AF [58][79][80][81] or the fluid exchange and osmotic swelling properties of the NP [82][83][84][85]. These (and other) studies are contributing to identifying the complexity of FSU tissues and to understanding their interaction, ranging from vertebroplasty support to the importance of the AF criss-cross angles in the herniation process, despite the frequent simplifications on the modelling of surrounding tissues (i.e., the FSU components that are not the focus of such study).

The common ground is that although numerical methods are increasingly established for tissue characterisation, the experimental support (through cadaveric samples or in vivo imaging) are more important than ever to validate these achievements and to obtain the approval by regulatory agencies for their use in clinical practice [47,56,88,89,90].

The common ground is that although numerical methods are increasingly established for tissue characterisation, the experimental support (through cadaveric samples or in vivo imaging) are more important than ever to validate these achievements and to obtain the approval by regulatory agencies for their use in clinical practice [47][56][86][87][88].

3. Challenges

Similar to tissue characterisation, FE modelling is now an essential tool for the development of TE solutions. This review intends to look at the application of FE modelling in TE studies, instead of looking at the multiple TE works dedicated to human spine tissue treatment or substitution.

TE substitutes to the IVD must mimic its complex environment and the interactions between the major components under native and degenerated conditions, which leads to a considerable number of design requisites. The IVD is an osmo-poro-hyper-viscoelastic organ, with anisotropic fibre reinforcement, so its mechanical functionality in terms of load support and flexibility needs to be ensured in arthroplasty or partial regeneration (e.g., NP replacement with hydrogel) [29,91,92].

TE substitutes to the IVD must mimic its complex environment and the interactions between the major components under native and degenerated conditions, which leads to a considerable number of design requisites. The IVD is an osmo-poro-hyper-viscoelastic organ, with anisotropic fibre reinforcement, so its mechanical functionality in terms of load support and flexibility needs to be ensured in arthroplasty or partial regeneration (e.g., NP replacement with hydrogel) [29][89][90].

The work of Silva et al. in 2005 [93] was one of the first to include FE simulation in the evaluation of a TE hydrogel for IVD replacement, the hydrogel being supported by an annular ring, i.e., mimicking the native NP–AF structure. Poro-elastic simulation allowed for the verification of stress patterns generated by this substitute under physiological fluid exchange. Strange et al. in 2010 [94] developed an FE model to assess the behaviour of an elastomeric NP replacement, assuming that such implant would have different mechanical properties than the native IVD, namely a different compressive stiffness. The focus was again on the stress distribution patterns: the NP natively distributes the pressure to the AF, which is then responsible for supporting the circadian variation. Different nucleotomy procedures/implant sizes were studied towards achieving the most appropriate solution, before advancing to actual implant production and approval. Schmidt et al. in 2014 [95] also investigated the effect of geometry/volumetric properties of an NP substitute, but here mimicking the poro-elastic medium of the IVD and performing an overall assessment of the implant behaviour, i.e., evaluating its effects on the adjacent FSU components. Kang et al. in 2013 [96] and Zhang et al. in 2018 [97] evaluated the design or porosity of interbody fusion cages towards achieving an optimal solution for each specific spinal level. Lim et al. in 2019 [98] addressed the issue of cage subsidence by combining additive manufacturing techniques and FE simulation to compare the behaviour of solid vs. porous devices built with PEEK–titanium composites, highlighting the biomechanical advantages of using a porous structure that could be built only with titanium. Their study is important to look at the multiple possibilities of computational simulations in TE, both for design and/or for testing of solutions.

The work of Silva et al. in 2005 [91] was one of the first to include FE simulation in the evaluation of a TE hydrogel for IVD replacement, the hydrogel being supported by an annular ring, i.e., mimicking the native NP–AF structure. Poro-elastic simulation allowed for the verification of stress patterns generated by this substitute under physiological fluid exchange. Strange et al. in 2010 [92] developed an FE model to assess the behaviour of an elastomeric NP replacement, assuming that such implant would have different mechanical properties than the native IVD, namely a different compressive stiffness. The focus was again on the stress distribution patterns: the NP natively distributes the pressure to the AF, which is then responsible for supporting the circadian variation. Different nucleotomy procedures/implant sizes were studied towards achieving the most appropriate solution, before advancing to actual implant production and approval. Schmidt et al. in 2014 [93] also investigated the effect of geometry/volumetric properties of an NP substitute, but here mimicking the poro-elastic medium of the IVD and performing an overall assessment of the implant behaviour, i.e., evaluating its effects on the adjacent FSU components. Kang et al. in 2013 [94] and Zhang et al. in 2018 [95] evaluated the design or porosity of interbody fusion cages towards achieving an optimal solution for each specific spinal level. Lim et al. in 2019 [96] addressed the issue of cage subsidence by combining additive manufacturing techniques and FE simulation to compare the behaviour of solid vs. porous devices built with PEEK–titanium composites, highlighting the biomechanical advantages of using a porous structure that could be built only with titanium. Their study is important to look at the multiple possibilities of computational simulations in TE, both for design and/or for testing of solutions.

Still, most of these studies have focussed on achieving a coherent geometry to the intervertebral space available after discectomy or nucleotomy, which is of upmost importance to ensure soft integration of the implant and an appropriate response of the surrounding tissues. The step forward from generic to dedicated solutions is another important aspect of the most recent computationally based TE approaches as regards structural optimisation methods [52,96]. The advances in additive manufacturing have potentiated the applicability of these tailored solutions. As seen above, the other level of research is the material properties of the substitute and their configuration, which will also be pivotal for the success of the procedure. Additive manufacturing has a main role here, as the use of biocompatible materials in these technologies is increasingly frequent [38,99]. However, as happens with the modelling (and subsequent substitution) of every biological tissue, the dependence of the constitutive material parameters on the patient’s age, physical condition and clinical history has to be addressed in such TE strategies [100,101,102]. This means that it may not be a good strategy to replace an aged IVD with a TE substitute based on the characteristics of a young and barely degenerated IVD, which will most likely induce a different biomechanical response. If the option is for interbody fusion instead of arthroplasty, the replacement cage should be adapted to the environment of the affected level and not to its native geometry/material properties as regards rehabilitating the correct sagittal balance for that specific patient/situation [103,104].

Still, most of these studies have focussed on achieving a coherent geometry to the intervertebral space available after discectomy or nucleotomy, which is of upmost importance to ensure soft integration of the implant and an appropriate response of the surrounding tissues. The step forward from generic to dedicated solutions is another important aspect of the most recent computationally based TE approaches as regards structural optimisation methods [52][94]. The advances in additive manufacturing have potentiated the applicability of these tailored solutions. As seen above, the other level of research is the material properties of the substitute and their configuration, which will also be pivotal for the success of the procedure. Additive manufacturing has a main role here, as the use of biocompatible materials in these technologies is increasingly frequent [38][97]. However, as happens with the modelling (and subsequent substitution) of every biological tissue, the dependence of the constitutive material parameters on the patient’s age, physical condition and clinical history has to be addressed in such TE strategies [98][99][100]. This means that it may not be a good strategy to replace an aged IVD with a TE substitute based on the characteristics of a young and barely degenerated IVD, which will most likely induce a different biomechanical response. If the option is for interbody fusion instead of arthroplasty, the replacement cage should be adapted to the environment of the affected level and not to its native geometry/material properties as regards rehabilitating the correct sagittal balance for that specific patient/situation [101][102].

In other words, a given device may have been modelled according to the best practices and produced with expensive techniques and materials, but it may still fail if careful pre-analysis of the individual situation is not performed. FE modelling then appears as the key technique to pre-screen the mechanobiological environment prior to the clinical decision, despite the challenging conditions to achieve significant simulations.

References

  1. Whatley, B.R.; Wen, X. Intervertebral disc (IVD): Structure, degeneration, repair and regeneration. Mater. Sci. Eng. C 2012, 32, 61–77.
  2. Foster, N.E.; Anema, J.R.; Cherkin, D.; Chou, R.; Cohen, S.P.; Gross, D.P.; Ferreira, P.H.; Fritz, J.M.; Koes, B.W.; Peul, W.; et al. Prevention and treatment of low back pain: Evidence, challenges, and promising directions. Lancet 2018, 391, 2368–2383.
  3. Hartvigsen, J.; Hancock, M.J.; Kongsted, A.; Louw, Q.; Ferreira, M.L.; Genevay, S.; Hoy, D.; Karppinen, J.; Pransky, G.; Sieper, J.; et al. What low back pain is and why we need to pay attention. Lancet 2018, 391, 2356–2367.
  4. Geurts, J.W.; Willems, P.C.; Kallewaard, J.W.; Van Kleef, M.; Dirksen, C. The Impact of Chronic Discogenic Low Back Pain: Costs and Patients’ Burden. Pain Res. Manag. 2018, 2018.
  5. Urban, J.P.G.; Fairbank, J.C.T. Current perspectives on the role of biomechanical loading and genetics in development of disc degeneration and low back pain; a narrative review. J. Biomech. 2020, 102, 109573.
  6. Khan, A.N.; Jacobsen, H.E.; Khan, J.; Filippi, C.G.; Levine, M.; Lehman, R.A.; Riew, K.D.; Lenke, L.G.; Chahine, N.O. Inflammatory biomarkers of low back pain and disc degeneration: A review. Ann. N. Y. Acad. Sci. 2017, 1410, 68–84.
  7. Bayoglu, R.; Geeraedts, L.; Groenen, K.H.J.; Verdonschot, N.; Koopman, B.; Homminga, J. Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the thoracic and cervical regions of the human spine. J. Biomech. 2017, 58, 52–63.
  8. Castro, A.P.G.; Paul, C.P.L.; Detiger, S.E.L.; Smit, T.H.; van Royen, B.J.; Pimenta Claro, J.C.; Mullender, M.G.; Alves, J.L. Long-Term Creep Behavior of the Intervertebral Disk: Comparison between Bioreactor Data and Numerical Results. Front. Bioeng. Biotechnol. 2014, 2, 56.
  9. Ghezelbash, F.; Shirazi-Adl, A.; Baghani, M.; Eskandari, A.H. On the modeling of human intervertebral disc annulus fibrosus: Elastic, permanent deformation and failure responses. J. Biomech. 2020, 102, 109463.
  10. Walter, B.A.; Torre, O.M.; Laudier, D.; Naidich, T.P.; Hecht, A.C.; Iatridis, J.C. Form and function of the intervertebral disc in health and disease: A morphological and stain comparison study. J. Anat. 2015, 227, 707–716.
  11. Zirbel, S.A.; Stolworthy, D.K.; Howell, L.L.; Bowden, A.E. Intervertebral disc degeneration alters lumbar spine segmental stiffness in all modes of loading under a compressive follower load. Spine J. 2013, 13, 1134–1147.
  12. Oxland, T.R. Fundamental biomechanics of the spine--What we have learned in the past 25 years and future directions. J. Biomech. 2016, 49, 817–832.
  13. Widmer, J.; Cornaz, F.; Scheibler, G.; Spirig, J.M.; Snedeker, J.G.; Farshad, M. Biomechanical contribution of spinal structures to stability of the lumbar spine—novel biomechanical insights. Spine J. 2020.
  14. Palanca, M.; Ruspi, M.L.; Cristofolini, L.; Liebsch, C.; Villa, T.; Brayda-Bruno, M.; Galbusera, F.; Wilke, H.J.; Barbera, L. La The strain distribution in the lumbar anterior longitudinal ligament is affected by the loading condition and bony features: An in vitro full-field analysis. PLoS ONE 2020, 15, 1–21.
  15. Naserkhaki, S.; Arjmand, N.; Shirazi-Adl, A.; Farahmand, F.; El-Rich, M. Effects of eight different ligament property datasets on biomechanics of a lumbar L4-L5 finite element model. J. Biomech. 2018, 70, 33–42.
  16. Ebraheim, N.A.; Hassan, A.; Lee, M.; Xu, R. Functional anatomy of the lumbar spine. Semin. Pain Med. 2004, 2, 131–137.
  17. Loenen, A.C.Y.; Noriega, D.C.; Wills, C.R.; Noailly, J.; Nunley, P.D.; Kirchner, R.; Ito, K.; van Rietbergen, B. Misaligned spinal rods can induce high internal forces consistent with those observed to cause screw pullout and disc degeneration. Spine J. 2020, 21, 528–537.
  18. Ito, K.; Creemers, L. Mechanisms of intervertebral disk degeneration/injury and pain: A review. Glob. spine J. 2013, 3, 145–152.
  19. Mok, F.P.S.; Samartzis, D.; Karppinen, J.; Fong, D.Y.T.; Luk, K.D.K.; Cheung, K.M.C. Modic changes of the lumbar spine: Prevalence, risk factors, and association with disc degeneration and low back pain in a large-scale population-based cohort. Spine J. 2016, 16, 32–41.
  20. Vergroesen, P.-P.A.; Kingma, I.; Emanuel, K.S.; Hoogendoorn, R.J.W.; Welting, T.J.; van Royen, B.J.; van Dieën, J.H.; Smit, T.H. Mechanics and biology in intervertebral disc degeneration: A vicious circle. Osteoarthr. Cartil. 2015, 23, 1057–1070.
  21. Maher, C.; Underwood, M.; Buchbinder, R. Non-specific low back pain. Lancet 2017, 389, 736–747.
  22. Andreula, C.; Muto, M.; Leonardi, M. Interventional spinal procedures. Eur. J. Radiol. 2004, 50, 112–119.
  23. Sasso, R.C.; Best, N.M.; Mummaneni, P.V.; Reilly, T.M.; Hussain, S.M. Analysis of Operative Complications in a Series of 471 Anterior Lumbar Interbody Fusion Procedures. Spine (Phila. Pa. 1976) 2005, 30, 670–674.
  24. Teng, I.; Han, J.; Phan, K.; Mobbs, R. A meta-analysis comparing ALIF, PLIF, TLIF and LLIF. J. Clin. Neurosci. 2017, 44, 11–17.
  25. Mobbs, R.; Phan, K.; Malham, G. Lumbar interbody fusion: Techniques, indications and comparison of interbody fusion options including PLIF, TLIF, MI-TLIF, OLIF/ATP, LLIF and ALIF. J. Spine Surg. 2015, 1, 2–18.
  26. Lazennec, J.Y. Lumbar and cervical viscoelastic disc replacement: Concepts and current experience. World J. Orthop. 2020, 11, 345–356.
  27. Weber, K.T.; Jacobsen, T.D.; Maidhof, R.; Virojanapa, J.; Overby, C.; Bloom, O.; Quraishi, S.; Levine, M.; Chahine, N.O. Developments in intervertebral disc disease research: Pathophysiology, mechanobiology, and therapeutics. Curr. Rev. Musculoskelet. Med. 2015, 8, 18–31.
  28. Singha, K.; Singha, M. Biomechanism Profile of Intervertebral Discs (IVD): Strategies to Successful Tissue Engineering for Spinal Healing by Reinforced Composite Structure. J. Tissue Sci. Eng. 2012, 3, 1–14.
  29. Bron, J.L.; Helder, M.N.; Meisel, H.J.; Van Royen, B.J.; Smit, T.H. Repair, regenerative and supportive therapies of the annulus fibrosus: Achievements and challenges. Eur. Spine J. 2009, 18, 301–313.
  30. Weisse, B.; Aiyangar, A.K.; Affolter, C.; Gander, R.; Terrasi, G.P.; Ploeg, H. Determination of the translational and rotational stiffnesses of an L4-L5 functional spinal unit using a specimen-specific finite element model. J. Mech. Behav. Biomed. Mater. 2012, 13, 45–61.
  31. Ghezelbash, F.; Schmidt, H.; Shirazi-Adl, A.; El-Rich, M. Internal load-sharing in the human passive lumbar spine: Review of in vitro and finite element model studies. J. Biomech. 2020, 102, 109441.
  32. Brekelmans, W.A.M.; Poort, H.W.; Slooff, T.J.J.H. A New Method to Analyse the Mechanical Behaviour of Skeletal Parts. Acta Orthop. Scand. 1972, 43, 301–317.
  33. Belytschko, T.; Kulak, R.F.; Schultz, A.B.; Galante, J.O. Finite element stress analysis of an intervertebral disc. J. Biomech. 1974, 7, 277–285.
  34. Schmidt, H.; Galbusera, F.; Rohlmann, A.; Shirazi-Adl, A. What have we learned from finite element model studies of lumbar intervertebral discs in the past four decades? J. Biomech. 2013, 46, 2342–2355.
  35. Fradet, L.; Wang, X.; Lenke, L.G.; Aubin, C.E. Biomechanical analysis of proximal junctional failure following adult spinal instrumentation using a comprehensive hybrid modeling approach. Clin. Biomech. 2016, 39, 122–128.
  36. Monteiro, N.M.B.; Silva, M.P.T.; Folgado, J.O.M.G.; Melancia, J.P.L. Structural analysis of the intervertebral discs adjacent to an interbody fusion using multibody dynamics and finite element cosimulation. Multibody Syst. Dyn. 2011, 25, 245–270.
  37. Ambati, D.V.; Wright, E.K.; Lehman, R.A.; Kang, D.G.; Wagner, S.C.; Dmitriev, A.E. Bilateral pedicle screw fixation provides superior biomechanical stability in transforaminal lumbar interbody fusion: A finite element study. Spine J. 2015, 15, 1812–1822.
  38. Mobbs, R.J.; Coughlan, M.; Thompson, R.; Sutterlin, C.E.; Phan, K. The utility of 3D printing for surgical planning and patient-specific implant design for complex spinal pathologies: Case report. J. Neurosurg. Spine 2017, 26, 513–518.
  39. Dreischarf, M.; Schmidt, H.; Putzier, M.; Zander, T. Biomechanics of the L5-S1 motion segment after total disc replacement—Influence of iatrogenic distraction, implant positioning and preoperative disc height on the range of motion and loading of facet joints. J. Biomech. 2015, 48, 3283–3291.
  40. Demir, E.; Eltes, P.; Castro, A.P.G.; Lacroix, D.; Toktaş, İ. Finite element modelling of hybrid stabilization systems for the human lumbar spine. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2020.
  41. Castro, A.P.G.P.G.; Alves, J.L.L. Numerical implementation of an osmo-poro-visco-hyperelastic finite element solver: Application to the intervertebral disc. Comput. Methods Biomech. Biomed. Engin. 2020, 0, 1–13.
  42. Huynh, K.T.; Gao, Z.; Gibson, I.; Lu, W.F. Haptically integrated simulation of a finite element model of thoracolumbar spine combining offline biomechanical response analysis of intervertebral discs. CAD Comput. Aided Des. 2010, 42, 1151–1166.
  43. Chen, Y.; Dall׳Ara, E.; Sales, E.; Manda, K.; Wallace, R.; Pankaj, P.; Viceconti, M. Micro-CT based finite element models of cancellous bone predict accurately displacement once the boundary condition is well replicated: A validation study. J. Mech. Behav. Biomed. Mater. 2017, 65, 644–651.
  44. Disney, C.M.; Eckersley, A.; McConnell, J.C.; Geng, H.; Bodey, A.J.; Hoyland, J.A.; Lee, P.D.; Sherratt, M.J.; Bay, B.K. Synchrotron tomography of intervertebral disc deformation quantified by digital volume correlation reveals microstructural influence on strain patterns. Acta Biomater. 2019, 92, 290–304.
  45. Wang, S.; Park, W.M.; Gadikota, H.R.; Miao, J.; Kim, Y.H.; Wood, K.B.; Li, G. A combined numerical and experimental technique for estimation of the forces and moments in the lumbar intervertebral disc. Comput. Methods Biomech. Biomed. Eng. 2013, 16, 1278–1286.
  46. Finley, S.M.; Brodke, D.S.; Spina, N.T.; DeDen, C.A.; Ellis, B.J. FEBio finite element models of the human lumbar spine. Comput. Methods Biomech. Biomed. Eng. 2018, 21, 444–452.
  47. Lavecchia, C.E.; Espino, D.M.; Moerman, K.M.; Tse, K.M.; Robinson, D.; Lee, P.V.S.; Shepherd, D.E.T. Lumbar model generator: A tool for the automated generation of a parametric scalable model of the lumbar spine. J. R. Soc. Interface 2018, 15.
  48. Breau, C.; Shirazi-Adl, A.; de Guise, J. Reconstruction of a human ligamentous lumbar spine using CT images—A three-dimensional finite element mesh generation. Ann. Biomed. Eng. 1991, 19, 291–302.
  49. Chen, C.S.; Cheng, C.K.; Liu, C.L.; Lo, W.H. Stress analysis of the disc adjacent to interbody fusion in lumbar spine. Med. Eng. Phys. 2001, 23, 485–493.
  50. Choi, J.; Kim, S.; Shin, D.A. Biomechanical comparison of spinal fusion methods using interspinous process compressor and pedicle screw fixation system based on finite element method. J. Korean Neurosurg. Soc. 2016, 59, 91–97.
  51. Guo, L.; Yin, J. Finite element analysis and design of an interspinous device using topology optimization. Med. Biol. Eng. Comput. 2019, 57, 89–98.
  52. Coelho, P.G.; Fernandes, P.C.; Folgado, J.; Fernandes, P.R. Development of a spinal fusion cage by multiscale modelling: Application to the human cervical spine. Procedia Eng. 2015, 110, 183–190.
  53. Zhong, Z.-C.; Wei, S.-H.; Wang, J.-P.; Feng, C.-K.; Chen, C.-S.; Yu, C. Finite element analysis of the lumbar spine with a new cage using a topology optimization method. Med. Eng. Phys. 2006, 28, 90–98.
  54. Lu, T.; Lu, Y. Comparison of Biomechanical Performance Among Posterolateral Fusion and Transforaminal, Extreme, and Oblique Lumbar Interbody Fusion: A Finite Element Analysis. World Neurosurg. 2019, 129, e890–e899.
  55. Zhou, C.; Cha, T.; Wang, W.; Guo, R.; Li, G. Investigation of Alterations in the Lumbar Disc Biomechanics at the Adjacent Segments After Spinal Fusion Using a Combined In Vivo and In Silico Approach. Ann. Biomed. Eng. 2020, 49, 601–616.
  56. Stokes, I.A.F.; Gardner-Morse, M. A database of lumbar spinal mechanical behavior for validation of spinal analytical models. J. Biomech. 2016, 49, 780–785.
  57. Hollingsworth, N.T.; Wagner, D.R. Modeling shear behavior of the annulus fibrosus. J. Mech. Behav. Biomed. Mater. 2011, 4, 1103–1114.
  58. Newell, N.; Carpanen, D.; Grigoriadis, G.; Little, J.P.; Masouros, S.D. Material properties of human lumbar intervertebral discs across strain rates. Spine J. 2019, 19, 2013–2024.
  59. Roriz, P.; Abe, I.; Schiller, M.; Gabriel, J.; Simões, J. Ex Vivo Intervertebral Disc Bulging Measurement Using a Fibre Bragg Grating Sensor. Exp. Mech. 2011, 51, 1573–1577.
  60. Jacobs, N.T.; Cortes, D.H.; Peloquin, J.M.; Vresilovic, E.J.; Elliott, D.M. Validation and application of an intervertebral disc finite element model utilizing independently constructed tissue-level constitutive formulations that are nonlinear, anisotropic, and time-dependent. J. Biomech. 2014, 47, 2540–2546.
  61. Yu, L.P.; Qian, W.W.; Yin, G.Y.; Ren, Y.X.; Hu, Z.Y. MRI Assessment of Lumbar Intervertebral Disc Degeneration with Lumbar Degenerative Disease Using the Pfirrmann Grading Systems. PLoS ONE 2012, 7, 1–7.
  62. Costa, M.C.; Eltes, P.; Lazary, A.; Varga, P.P.; Viceconti, M.; Dall’Ara, E. Biomechanical assessment of vertebrae with lytic metastases with subject-specific finite element models. J. Mech. Behav. Biomed. Mater. 2019, 98, 268–290.
  63. O’Reilly, M.A.; Whyne, C.M. Comparison of computed tomography based parametric and patient-specific finite element models of the healthy and metastatic spine using a mesh-morphing algorithm. Spine (Phila. Pa. 1976) 2008, 33, 1876–1881.
  64. Wognum, S.; Huyghe, J.M.; Baaijens, F.P.T. Influence of osmotic pressure changes on the opening of existing cracks in 2 intervertebral disc models. Spine (Phila. Pa. 1976) 2006, 31, 1783–1788.
  65. Vergroesen, P.-P.A.; Emanuel, K.S.; Peeters, M.; Kingma, I.; Smit, T.H. Are axial intervertebral disc biomechanics determined by osmosis? J. Biomech. 2018, 70, 4–9.
  66. Paul, C.P.L.; Emanuel, K.S.; Kingma, I.; van der Veen, A.J.; Holewijn, R.M.; Vergroesen, P.-P.A.; van de Ven, P.M.; Mullender, M.G.; Helder, M.N.; Smit, T.H. Changes in Intervertebral Disk Mechanical Behavior During Early Degeneration. J. Biomech. Eng. 2018, 140, 091008.
  67. Qasim, M.; Natarajan, R.N.; An, H.S.; Andersson, G.B.J. Initiation and progression of mechanical damage in the intervertebral disc under cyclic loading using continuum damage mechanics methodology: A finite element study. J. Biomech. 2012, 45, 1934–1940.
  68. Schroeder, Y.; Huyghe, J.M.; Van Donkelaar, C.C.; Ito, K. A biochemical/biophysical 3D FE intervertebral disc model. Biomech. Model. Mechanobiol. 2010, 9, 641–650.
  69. Barthelemy, V.M.P.; van Rijsbergen, M.M.; Wilson, W.; Huyghe, J.M.; van Rietbergen, B.; Ito, K. A computational spinal motion segment model incorporating a matrix composition-based model of the intervertebral disc. J. Mech. Behav. Biomed. Mater. 2016, 54, 194–204.
  70. Velísková, P.; Bashkuev, M.; Shirazi-Adl, A.; Schmidt, H. Computational study of the role of fluid content and flow on the lumbar disc response in cyclic compression: Replication of in vitro and in vivo conditions. J. Biomech. 2017.
  71. Massey, C.J.; Van Donkelaar, C.C.; Vresilovic, E.; Zavaliangos, A.; Marcolongo, M. Effects of aging and degeneration on the human intervertebral disc during the diurnal cycle: A finite element study. J. Orthop. Res. 2012, 30, 122–128.
  72. Malandrino, A.; Pozo, J.M.; Castro-Mateos, I.; Frangi, A.F.; van Rijsbergen, M.M.; Ito, K.; Wilke, H.-J.; Dao, T.T.; Ho Ba Tho, M.-C.; Noailly, J. On the Relative Relevance of Subject-Specific Geometries and Degeneration-Specific Mechanical Properties for the Study of Cell Death in Human Intervertebral Disk Models. Front. Bioeng. Biotechnol. 2015, 3, 1–15.
  73. Fattor, J.A.; Hollenbeck, J.F.M.; Laz, P.J.; Rullkoetter, P.J.; Burger, E.L.; Patel, V.V.; Cain, C.M.J. Patient-specific templating of lumbar total disk replacement to restore normal anatomy and function. Orthopedics 2016, 39, 97–102.
  74. El Bojairami, I.; El-Monajjed, K.; Driscoll, M. Development and validation of a timely and representative finite element human spine model for biomechanical simulations. Sci. Rep. 2020, 10, 1–15.
  75. Van Rijsbergen, M.; Van Rietbergen, B.; Barthelemy, V.; Eltes, P.; Lazáry, Á.; Lacroix, D.; Noailly, J.; Tho, M.C.H.B.; Wilson, W.; Ito, K. Comparison of patient-specific computational models vs. clinical follow-up, for adjacent segment disc degeneration and bone remodelling after spinal fusion. PLoS ONE 2018, 13, 1–24.
  76. Dall’Ara, E.; Karl, C.; Mazza, G.; Franzoso, G.; Vena, P.; Pretterklieber, M.; Pahr, D.; Zysset, P. Tissue properties of the human vertebral body sub-structures evaluated by means of microindentation. J. Mech. Behav. Biomed. Mater. 2013, 25, 23–32.
  77. Calvo-Echenique, A.; Cegoñino, J.; Chueca, R.; Pérez-del Palomar, A. Stand-alone lumbar cage subsidence: A biomechanical sensitivity study of cage design and placement. Comput. Methods Programs Biomed. 2018, 162, 211–219.
  78. Jackman, T.M.; DelMonaco, A.M.; Morgan, E.F. Accuracy of finite element analyses of CT scans in predictions of vertebral failure patterns under axial compression and anterior flexion. J. Biomech. 2016, 49, 267–275.
  79. Eberlein, R.; Holzapfel, G.A.; Schulze-Bauer, C.A.J. An Anisotropic Model for Annulus Tissue and Enhanced Finite Element Analyses of Intact Lumbar Disc Bodies. Comput. Methods Biomech. Biomed. Engin. 2001, 4, 209–229.
  80. Vergari, C.; Chan, D.; Clarke, A.; Mansfield, J.C.; Meakin, J.R.; Winlove, P.C. Bovine and degenerated human annulus fibrosus: A microstructural and micromechanical comparison. Biomech. Model. Mechanobiol. 2017, 1–10.
  81. Noailly, J.; Planell, J.A.; Lacroix, D. On the collagen criss-cross angles in the annuli fibrosi of lumbar spine finite element models. Biomech. Model. Mechanobiol. 2011, 10, 203–219.
  82. Ruiz, C.; Noailly, J.; Lacroix, D. Material property discontinuities in intervertebral disc porohyperelastic finite element models generate numerical instabilities due to volumetric strain variations. J. Mech. Behav. Biomed. Mater. 2013, 26, 1–10.
  83. Galbusera, F.; Schmidt, H.; Noailly, J.; Malandrino, A.; Lacroix, D.; Wilke, H.J.; Shirazi-Adl, A. Comparison of four methods to simulate swelling in poroelastic finite element models of intervertebral discs. J. Mech. Behav. Biomed. Mater. 2011, 4, 1234–1241.
  84. Marini, G.; Studer, H.; Huber, G.; Püschel, K.; Ferguson, S.J. Geometrical aspects of patient-specific modelling of the intervertebral disc: Collagen fibre orientation and residual stress distribution. Biomech. Model. Mechanobiol. 2016, 15, 543–560.
  85. Bezci, S.E.; O’Connell, G.D. Osmotic Pressure Alters Time-dependent Recovery Behavior of the Intervertebral Disc. Spine (Phila. Pa. 1976) 2018, 43, E334–E340.
  86. Berton, A.; Salvatore, G.; Giambini, H.; Ciuffreda, M.; Longo, U.G.; Denaro, V.; Thoreson, A.; An, K.-N. A 3D finite element model of prophylactic vertebroplasty in the metastatic spine: Vertebral stability and stress distribution on adjacent vertebrae. J. Spinal Cord Med. 2020, 43, 39–45.
  87. Viceconti, M. Biomechanics-based in silico medicine: The manifesto of a new science. J. Biomech. 2015, 48, 193–194.
  88. Agarwal, A.; Kodigudla, M.; Kelkar, A.; Jayaswal, D.; Goel, V. North American Spine Society Journal (NASSJ) Towards a validated patient-specific computational modeling framework to identify failure regions in traditional growing rods in patients with early onset scoliosis. North Am. Spine Soc. J. 2021, 5, 100043.
  89. Huang, J.; Yan, H.; Jian, F.; Wang, X.; Li, H. Numerical analysis of the influence of nucleus pulposus removal on the biomechanical behavior of a lumbar motion segment. Comput. Methods Biomech. Biomed. Engin. 2014, 5842, 1–9.
  90. Reitmaier, S.; Shirazi-Adl, A.; Bashkuev, M.; Wilke, H.-J.; Gloria, A.; Schmidt, H. In vitro and in silico investigations of disc nucleus replacement. J. R. Soc. Interface 2012, 9, 1869–1879.
  91. Silva, P.; Crozier, S.; Veidt, M.; Pearcy, M.J. An experimental and finite element poroelastic creep response analysis of an intervertebral hydrogel disc model in axial compression. J. Mater. Sci. Mater. Med. 2005, 16, 663–669.
  92. Strange, D.G.T.; Fisher, S.T.; Boughton, P.C.; Kishen, T.J.; Diwan, A.D. Restoration of compressive loading properties of lumbar discs with a nucleus implant-a finite element analysis study. Spine J. 2010, 10, 602–609.
  93. Schmidt, H.; Bashkuev, M.; Galbusera, F.; Wilke, H.-J.; Shirazi-Adl, A. Finite element study of human lumbar disc nucleus replacements. Comput. Methods Biomech. Biomed. Engin. 2014, 17, 1762–1776.
  94. Kang, H.; Hollister, S.J.; La Marca, F.; Park, P.; Lin, C.-Y. Porous Biodegradable Lumbar Interbody Fusion Cage Design and Fabrication Using Integrated Global-Local Topology Optimization With Laser Sintering. J. Biomech. Eng. 2013, 135, 101013.
  95. Zhang, Z.; Fogel, G.R.; Liao, Z.; Sun, Y.; Liu, W. Biomechanical analysis of lumbar interbody fusion cages with various lordotic angles: A finite element study. Comput. Methods Biomech. Biomed. Engin. 2018, 21, 247–254.
  96. Lim, K.M.; Park, T.H.; Lee, S.J.; Park, S.J. Design and biomechanical verification of additive manufactured composite spinal cage composed of porous titanium cover and PEEK body. Appl. Sci. 2019, 9, 4258.
  97. Liang, H.; Li, X.; Shimer, A.L.; Balian, G.; Shen, F.H. A novel strategy of spine defect repair with a degradable bioactive scaffold preloaded with adipose-derived stromal cells. Spine J. 2014, 14, 445–454.
  98. Reutlinger, C.; Bürki, A.; Brandejsky, V.; Ebert, L.; Büchler, P. Specimen specific parameter identification of ovine lumbar intervertebral discs: On the influence of fibre-matrix and fibre-fibre shear interactions. J. Mech. Behav. Biomed. Mater. 2014, 30, 279–289.
  99. O’Connell, G.D.; Sen, S.; Elliott, D.M. Human annulus fibrosus material properties from biaxial testing and constitutive modeling are altered with degeneration. Biomech. Model. Mechanobiol. 2012, 11, 493–503.
  100. Bezci, S.E.; Lim, S.; O’Connell, G.D. Nonlinear stress-dependent recovery behavior of the intervertebral disc. J. Mech. Behav. Biomed. Mater. 2020, 110, 103881.
  101. Caprara, S.; Moschini, G.; Snedeker, J.G.; Farshad, M.; Senteler, M. Spinal sagittal alignment goals based on statistical modelling and musculoskeletal simulations. J. Biomech. 2020, 102, 109621.
  102. Galbusera, F.; Wilke, H.J.; Brayda-Bruno, M.; Costa, F.; Fornari, M. Influence of sagittal balance on spinal lumbar loads: A numerical approach. Clin. Biomech. 2013, 28, 370–377.
More
ScholarVision Creations