Next Article in Journal
An Open-Source COVID-19 CT Dataset with Automatic Lung Tissue Classification for Radiomics
Previous Article in Journal
Biocompatible Electrospun Polycaprolactone-Polyaniline Scaffold Treated with Atmospheric Plasma to Improve Hydrophilicity
Previous Article in Special Issue
Mineralization of Biomaterials for Bone Tissue Engineering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Computational Challenges in Tissue Engineering for the Spine

IDMEC, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal
Submission received: 31 December 2020 / Revised: 4 February 2021 / Accepted: 13 February 2021 / Published: 15 February 2021
(This article belongs to the Special Issue Biomaterials for Bone Tissue Engineering)

Abstract

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

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]. 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].
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].
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]. 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].
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]. 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]. 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].
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]. A lack of widely accepted explanatory models limits the understanding of such diseases, retarding the development of effective therapies [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]. 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]. 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].
Therefore, this brief review intends to demonstrate the growing importance of computational methods in the study of the healthy, degenerated and implanted human spine, as well as to identify the next generation of simulation-based TE studies towards the development of increasingly adequate treatment solutions for DSD and LBP.

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].
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]. 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,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, 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.
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].
Table 1 lists the most common constitutive models and data sources for different FSU tissues. It must be highlighted that despite the complex morphology of bone, cortical bone is usually modelled as linear elastic and isotropic; the local trabeculae of trabecular bone are markedly anisotropic, but such behaviour is not frequently taken into consideration in FSU modelling.
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.
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 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 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.
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).
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].

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].
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.
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].
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.

4. Conclusions

To conclude, it has been demonstrated here that FE modelling is nowadays an essential tool to pre-screen the specificities of a patient suffering from LBP, particularly the targeted tissues. It is also important to evaluate the potential TE solutions to be employed. Still, with growing computational power and developing imaging methods, the contribution of medical imaging examinations is decisive, along with the also extremely important contribution of clinical essays. Therefore, the major challenge in computational modelling for TE solutions to tackle DSD is to achieve comprehensive strategies that include the tissue properties and the functionality target. TE developments still have a way to go before full clinical practice implementation and fast regulatory approval, but the association of FE modelling is definitely the most promising approach to achieve successful patient-specific solutions.

Funding

This work was supported by Fundação para a Ciência e Tecnologia, through IDMEC, under the LAETA project UIDB/50022/2020.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Whatley, B.R.; Wen, X. Intervertebral disc (IVD): Structure, degeneration, repair and regeneration. Mater. Sci. Eng. C 2012, 32, 61–77. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  16. Ebraheim, N.A.; Hassan, A.; Lee, M.; Xu, R. Functional anatomy of the lumbar spine. Semin. Pain Med. 2004, 2, 131–137. [Google Scholar] [CrossRef]
  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. [Google Scholar]
  18. Ito, K.; Creemers, L. Mechanisms of intervertebral disk degeneration/injury and pain: A review. Glob. spine J. 2013, 3, 145–152. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  21. Maher, C.; Underwood, M.; Buchbinder, R. Non-specific low back pain. Lancet 2017, 389, 736–747. [Google Scholar] [CrossRef] [Green Version]
  22. Andreula, C.; Muto, M.; Leonardi, M. Interventional spinal procedures. Eur. J. Radiol. 2004, 50, 112–119. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [PubMed]
  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. [Google Scholar]
  26. Lazennec, J.Y. Lumbar and cervical viscoelastic disc replacement: Concepts and current experience. World J. Orthop. 2020, 11, 345–356. [Google Scholar] [PubMed]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hollingsworth, N.T.; Wagner, D.R. Modeling shear behavior of the annulus fibrosus. J. Mech. Behav. Biomed. Mater. 2011, 4, 1103–1114. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  64. 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. [Google Scholar] [CrossRef] [Green Version]
  65. 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. [Google Scholar] [CrossRef]
  66. 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. [Google Scholar] [CrossRef]
  67. Toumanidou, T.; Noailly, J. Musculoskeletal Modeling of the Lumbar Spine to Explore Functional Interactions between Back Muscle Loads and Intervertebral Disk Multiphysics. Front. Bioeng. Biotechnol. 2015, 3, 1–13. [Google Scholar]
  68. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. 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. [Google Scholar] [CrossRef] [PubMed]
  70. 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. [Google Scholar] [CrossRef] [Green Version]
  71. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. 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. [Google Scholar] [CrossRef] [PubMed]
  73. 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. [Google Scholar] [CrossRef]
  74. 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. [Google Scholar] [CrossRef]
  75. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Zanjani-Pour, S.; Winlove, C.P.; Smith, C.W.; Meakin, J.R. Image driven subject-specific finite element models of spinal biomechanics. J. Biomech. 2016, 49, 919–925. [Google Scholar] [CrossRef]
  77. 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. [Google Scholar] [CrossRef] [Green Version]
  78. 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. [Google Scholar] [CrossRef]
  79. 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. [Google Scholar] [CrossRef]
  80. 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. [Google Scholar] [CrossRef] [Green Version]
  81. 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. [Google Scholar] [CrossRef] [Green Version]
  82. 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. [Google Scholar] [CrossRef]
  83. 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. [Google Scholar] [CrossRef]
  84. 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. [Google Scholar] [CrossRef] [PubMed]
  85. 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. [Google Scholar] [CrossRef] [PubMed]
  86. 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. [Google Scholar] [CrossRef]
  87. 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. [Google Scholar] [CrossRef]
  88. 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. [Google Scholar] [CrossRef] [PubMed]
  89. Viceconti, M. Biomechanics-based in silico medicine: The manifesto of a new science. J. Biomech. 2015, 48, 193–194. [Google Scholar] [CrossRef]
  90. 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. [Google Scholar]
  91. 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. [Google Scholar] [CrossRef] [PubMed]
  92. 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. [Google Scholar] [CrossRef]
  93. 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. [Google Scholar] [CrossRef] [PubMed]
  94. 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. [Google Scholar] [CrossRef]
  95. 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. [Google Scholar] [CrossRef] [PubMed]
  96. 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. [Google Scholar] [CrossRef] [Green Version]
  97. 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. [Google Scholar] [CrossRef]
  98. 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. [Google Scholar] [CrossRef] [Green Version]
  99. 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. [Google Scholar]
  100. 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. [Google Scholar] [CrossRef]
  101. 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. [Google Scholar] [CrossRef] [Green Version]
  102. 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. [Google Scholar] [CrossRef] [PubMed]
  103. 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. [Google Scholar] [CrossRef] [PubMed]
  104. 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. [Google Scholar] [CrossRef] [PubMed]
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.
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.
Bioengineering 08 00025 g001
Table 1. Common constitutive modelling options and preferential imaging methods to characterise functional spinal unit (FSU) tissues.
Table 1. Common constitutive modelling options and preferential imaging methods to characterise functional spinal unit (FSU) tissues.
ComponentConstitutive ModelImaging Method
NPHyper-elastic, poro-elastic with osmotic swelling [41,64]MRI
AFHyper-elastic, anisotropic (fibre reinforcement), poro-elastic with osmotic swelling [41,65]MRI
CEPHyper-elastic, poro-elastic [41]MRI
Cortical boneLinear elastic, orthotropic [66]CT
Trabecular boneHyper-elastic, anisotropic [66]CT
LigamentsNon-linear, tension-only [17] MRI
MusclesNon-linear [67]MRI
MRI: Magnetic resonance imaging; CT: Computed tomography.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Castro, A.P.G. Computational Challenges in Tissue Engineering for the Spine. Bioengineering 2021, 8, 25. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering8020025

AMA Style

Castro APG. Computational Challenges in Tissue Engineering for the Spine. Bioengineering. 2021; 8(2):25. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering8020025

Chicago/Turabian Style

Castro, André P. G. 2021. "Computational Challenges in Tissue Engineering for the Spine" Bioengineering 8, no. 2: 25. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering8020025

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop