<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Bioeng. Biotechnol.</journal-id>
<journal-title>Frontiers in Bioengineering and Biotechnology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Bioeng. Biotechnol.</abbrev-journal-title>
<issn pub-type="epub">2296-4185</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1294658</article-id>
<article-id pub-id-type="doi">10.3389/fbioe.2024.1294658</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Bioengineering and Biotechnology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Analyzing isolated degeneration of lumbar facet joints: implications for degenerative instability and lumbar biomechanics using finite element analysis</article-title>
<alt-title alt-title-type="left-running-head">Park et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fbioe.2024.1294658">10.3389/fbioe.2024.1294658</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Park</surname>
<given-names>Jun Sung</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>&#x2020;</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Goh</surname>
<given-names>Tae Sik</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1461215/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Lee</surname>
<given-names>Jung Sub</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Lee</surname>
<given-names>Chiseung</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1867162/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Biomedical Engineering</institution>, <institution>Graduate School</institution>, <institution>Pusan National University</institution>, <addr-line>Busan</addr-line>, <country>Republic of Korea</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Department of Orthopaedic Surgery</institution>, <institution>School of Medicine</institution>, <institution>Pusan National University</institution>, <addr-line>Busan</addr-line>, <country>Republic of Korea</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Biomedical Research Institute</institution>, <institution>Pusan National University Hospital</institution>, <addr-line>Busan</addr-line>, <country>Republic of Korea</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Department of Biomedical Engineering</institution>, <institution>School of Medicine</institution>, <institution>Pusan National University</institution>, <addr-line>Busan</addr-line>, <country>Republic of Korea</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/418863/overview">Alexandros E. Tsouknidas</ext-link>, University of Western Macedonia, Greece</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/234327/overview">Francesco Travascio</ext-link>, University of Miami, United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1504517/overview">Zhonghai Li</ext-link>, First Affiliated Hospital of Dalian Medical University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Jung Sub Lee, <email>jungsublee@pusan.ac.kr</email>; Chiseung Lee, <email>victorich@pusan.ac.kr</email>
</corresp>
<fn fn-type="equal" id="fn001">
<label>
<sup>&#x2020;</sup>
</label>
<p>These authors have contributed equally to this work</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>27</day>
<month>03</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1294658</elocation-id>
<history>
<date date-type="received">
<day>15</day>
<month>09</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>26</day>
<month>02</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Park, Goh, Lee and Lee.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Park, Goh, Lee and Lee</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The facet joint contributes to lumbar spine stability as it supports the weight of body along with the intervertebral discs. However, most studies on the causes of degenerative lumbar diseases focus on the intervertebral discs and often overlook the facet joints. This study aimed to investigate the impact of facet joint degeneration on the degenerative changes and diseases of the lumbar spine. A finite element model of the lumbar spine (L1&#x2013;S1) was fabricated and validated to study the biomechanical characteristics of the facet joints. To simulate degeneration of the facet joint, the model was divided into four grades based on the number of degenerative segments (L4&#x2013;L5 or L4&#x2013;S1) and the contact condition between the facet joint surfaces. Finite element analysis was performed on four spine motions: flexion, extension, lateral bending, and axial torsion, by applying a pure moment to the upper surface of L1. Important parameters that could be used to confirm the effect of facet joint degeneration on the lumbar spine were calculated, including the range of motion (ROM) of the lumbar segments, maximum von Mises stress on the intervertebral discs, and reaction force at the facet joint. Facet joint degeneration affected the biomechanical characteristics of the lumbar spine depending on the movements of the spine. When analyzed by dividing it into degenerative onset and onset-adjacent segments, lumbar ROM and the maximum von Mises stress of the intervertebral discs decreased as the degree of degeneration increased in the degenerative onset segments. The reaction force at the facet joint decreased with flexion and increased with lateral bending and axial torsion. In contrast, lumbar ROM of the onset-adjacent segments remained almost unchanged despite severe degeneration of the facet joint, and the maximum von Mises stress of the intervertebral discs increased with flexion and extension but decreased with lateral bending and axial torsion. Additionally, the facet joint reaction force increased with extension, lateral bending, and axial rotation. This analysis, which combined the ROM of the lumbar segment, maximum von Mises stress on the intervertebral disc, and facet joint reaction force, confirmed the biomechanical changes in the lumbar spine due to the degeneration of isolated facet joints under the load of spinal motion. In the degenerative onset segment, spinal instability decreased, whereas in the onset-adjacent segment, a greater load was applied than in the intact state. When conducting biomechanical studies on the lumbar spine, considering facet joint degeneration is important since it can lead to degenerative spinal diseases, including adjacent segment diseases.</p>
</abstract>
<kwd-group>
<kwd>lumbar facet joint</kwd>
<kwd>isolated degeneration</kwd>
<kwd>degenerative instability</kwd>
<kwd>lumbar biomechanics</kwd>
<kwd>finite element analysis</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Biomechanics</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1 Introduction</title>
<p>The facet joint is a functional joint that consists of a pair of zygapophyseal joints located at the posterior part of the vertebrae (<xref ref-type="bibr" rid="B18">Gellhorn et al., 2013</xref>; <xref ref-type="bibr" rid="B49">Yin et al., 2020</xref>). It maintains the stability of spinal motion in the functional spinal unit (FSU) of each vertebral segment and supports approximately 6%&#x2013;30% of the axial compressive load; therefore, degenerative changes can occur in the facet joint when subjected to excessive loads (<xref ref-type="bibr" rid="B43">Takigawa et al., 2010</xref>; <xref ref-type="bibr" rid="B20">Iorio et al., 2015</xref>; <xref ref-type="bibr" rid="B49">Yin et al., 2020</xref>). Degeneration of the facet joint mainly occurs at the L4&#x2013;L5 and L5&#x2013;S1 spinal segments; moreover, pathological changes, such as narrowing of the facet joint space, hypertrophy of the articular process, subchondral cysts, osteophyte formation, and subarticular bone erosions, can be observed on medical images (<xref ref-type="bibr" rid="B48">Weishaupt et al., 1999</xref>; <xref ref-type="bibr" rid="B24">Kalichman et al., 2008</xref>). The incidence and severity of degeneration tend to increase with age; furthermore, the severity of degeneration is evaluated using the grading system for facet joint degeneration (FJD) developed by <xref ref-type="bibr" rid="B48">Weishaupt et al. (1999)</xref>, <xref ref-type="bibr" rid="B13">Eubanks et al. (2007)</xref>.</p>
<p>Degenerative changes in the facet joint can also be accelerated by other spinal disorders, such as degenerative disc disease (DDD) (<xref ref-type="bibr" rid="B16">Fujiwara et al., 1999</xref>; <xref ref-type="bibr" rid="B22">Jaumard et al., 2011</xref>; <xref ref-type="bibr" rid="B44">Varlotta et al., 2011</xref>; <xref ref-type="bibr" rid="B18">Gellhorn et al., 2013</xref>). DDD causes a reduction in the height of the intervertebral discs, which increases the compressive load and pressure transmitted to the facet joint, and this eventually damages facet joints (<xref ref-type="bibr" rid="B9">Dunlop et al., 1984</xref>; <xref ref-type="bibr" rid="B34">Panjabi et al., 1984</xref>; <xref ref-type="bibr" rid="B28">Li et al., 2011</xref>). Despite the high incidence of degenerative facet joint disease due to various reasons, relatively few patients with FJD are diagnosed without DDD. Previous studies on FJD have also included patients with degenerative changes in the intervertebral discs. Therefore, conducting studies that focus solely on FJD is difficult.</p>
<p>Recently, a significant increase in the use of computational biomechanics techniques to model and simulate the spine (<xref ref-type="bibr" rid="B15">Fujiwara et al., 2000</xref>; <xref ref-type="bibr" rid="B8">Du et al., 2016</xref>; <xref ref-type="bibr" rid="B47">Wang et al., 2018</xref>; <xref ref-type="bibr" rid="B46">Wang et al., 2019</xref>; <xref ref-type="bibr" rid="B42">Song et al., 2021</xref>; <xref ref-type="bibr" rid="B25">Kim et al., 2022</xref>; <xref ref-type="bibr" rid="B41">Sim et al., 2022</xref>) have been reported. Studies focusing on the lumbar facet joints have also been conducted with increasing frequency. <xref ref-type="bibr" rid="B8">Du et al. (2016)</xref> investigated the changes in the biomechanical characteristics of the facet joint when various levels of follower preload were applied to the lumbar L1&#x2013;L5 finite element (FE) model. <xref ref-type="bibr" rid="B46">Wang et al. (2019)</xref> developed an FE model of the C5&#x2013;C6 cervical segment and implemented FJD by modeling the stiffness of the capsular ligament and the angle of the zygapophyseal joint. They also evaluated the biomechanical characteristics of the cervical spine by performing simulations using these parameters and comparing and analyzing the range of motion (ROM) and intradiscal pressure (IDP) of the intervertebral disc. However, confirming the biomechanical characteristics of FJD at each stage was difficult since they were not classified into grades. Regarding the grades of degeneration, <xref ref-type="bibr" rid="B47">Wang et al. (2018)</xref> developed a lumbar L1&#x2013;S1 FE model and performed a series of simulations by adding anterior osteophytes, which are degenerative changes, at each segment. They evaluated the effects on the lumbar spine by grading the anterior osteophytes. <xref ref-type="bibr" rid="B15">Fujiwara et al. (2000)</xref> evaluated the biomechanical characteristics of the spine by analyzing the spine motion according to the degeneration grade of the facet joint and intervertebral disc on medical images. Many previous studies have evaluated the effects of degenerative changes in the spine. Nonetheless, there is a paucity of research that effectively differentiates between the grades of degeneration specifically within the lumbar facet joint, excluding considerations of intervertebral disc degeneration. Moreover, there is a dearth of studies assessing the biomechanical repercussions of lumbar facet joint degeneration on the lumbar spine.</p>
<p>Therefore, this study focused only on lumbar FJD. Using finite element analysis, we aimed to confirm the effects of various levels of lumbar FJD on lumbar activity. In order to simulate diverse degrees of degeneration, we utilized the grading system for FJD as a reference and assumed that the gap and friction coefficient of the articular surface were key variables (<xref ref-type="bibr" rid="B48">Weishaupt et al., 1999</xref>). We classified the lesions into five levels based on this grading system: intact state (grade 0; G0), mild (Grade 1; G1), mild-moderate (Grade 2; G2), moderate (Grade 3; G3), and severe generation (Grade 4; G4). In addition, we considered two cases with different numbers of degenerated segments, considering the location of the segments where degeneration mainly occurs in the facet joint: a single segment at the L4&#x2013;L5 and double segments at the L4&#x2013;S1 segments. Both of these correspond to the onset segment. And the onset adjacent segments are L3-L4 (upper segment), L5-S1 (lower segment) for single segment facet joint degeneration and L3-L4 (upper segment) for double segments facet joint degeneration. In this study, we aimed to strategically evaluate the impact of the degeneration grade and the number of onset segments on onset and onset-adjacent segments.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Materials and methods</title>
<p>The normalized shape of the lumbosacral spine (L1&#x2013;S1) was extracted from computed tomography (CT) images of the spine of a 28-year old male (80&#xa0;kg, 178&#xa0;cm) without metabolic bone disease using Mimics 21.0 (Materialise, Leuven, Belgium). The shape of the vertebral body (L1&#x2013;S1) was extracted using Geomagic Design X software (3D Systems, Rock Hill, SC, United States). Afterwards, noise and defects that occurred during shape extraction were removed from the medical images. These images were then exported as an STP file for additional detailed modeling using the Inventor 2018 (Autodesk, Mill Valley, CA, United States). The tasks carried out with Inventor included simplification of the vertebral body and additional modeling of intervertebral discs (annulus fibrosus and nucleus pulposus) and endplates to create a 3D geometric model, as shown in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Finite element model of the lumbar spine including the L1&#x2013;S1 spinal segment with vertebral bodies, intervertebral discs, and ligaments.</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g001.tif"/>
</fig>
<p>The generated 3D model was subjected to FEA preprocessing on an ANSYS Workbench 2022 R2 (ANSYS Inc. Canonsburg. PA. United States), a commercial finite element analysis program. The material properties of each component of the 3D model are listed in <xref ref-type="table" rid="T1">Table 1</xref> (<xref ref-type="bibr" rid="B19">Goel et al., 1995</xref>; <xref ref-type="bibr" rid="B31">Lu et al., 1996</xref>; <xref ref-type="bibr" rid="B38">Rohlmann et al., 2006</xref>; <xref ref-type="bibr" rid="B39">Rohlmann et al., 2009</xref>; <xref ref-type="bibr" rid="B29">Liu et al., 2011</xref>; <xref ref-type="bibr" rid="B6">Dreischarf et al., 2014</xref>; <xref ref-type="bibr" rid="B7">Du et al., 2014</xref>; <xref ref-type="bibr" rid="B33">Nikkhoo et al., 2020</xref>; <xref ref-type="bibr" rid="B32">Mengoni, 2021</xref>). Among the components of the FE model, vertebral bodies (cortical bone, cancellous bone, and posterior elements) and endplates were adopted as isotropic linear material properties. The cortical bone, cancellous bone, and endplate were hexahedral elements, and the posterior elements were tetrahedral elements. The annulus fibrosus and incompressible nucleus pulposus of the intervertebral disc were adopted as hyperelastic material models using Mooney Rivlin with tetrahedral elements. Additionally, seven representative ligaments were defined as nonlinear materials using tension-only spring elements: the anterior longitudinal ligament, posterior longitudinal ligament, interspinal ligament, supraspinal ligament, intertransverse ligament, ligamentum flavum, and capsular ligament (<xref ref-type="table" rid="T2">Table 2</xref>). Based on previous studies, the modeling of the intact facet joint assumes certain parameters, including a frictionless interface, an initial gap of 0.5&#xa0;mm, an exponential formulation, and a pressure of 120&#xa0;MPa at a zero gap (<xref ref-type="bibr" rid="B30">Liu et al., 2020</xref>; <xref ref-type="bibr" rid="B32">Mengoni, 2021</xref>). A finite element model (Model0; M0) of the intact state of lumbar L1&#x2013;S1 was constructed, as shown in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Material properties used in the developed FE model.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Components</th>
<th align="center">Young&#x2019;s modulus (MPa)</th>
<th align="center">Poisson&#x2019;s ratio</th>
<th align="center">References</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Cortical bone</td>
<td align="center">10,000</td>
<td align="center">0.3</td>
<td align="center">
<xref ref-type="bibr" rid="B39">Rohlmann et al. (2009)</xref>
</td>
</tr>
<tr>
<td align="center">Cancellous bone</td>
<td align="center">100</td>
<td align="center">0.2</td>
<td align="center">
<xref ref-type="bibr" rid="B19">Goel et al. (1995)</xref>
</td>
</tr>
<tr>
<td align="center">Posterior bone</td>
<td align="center">3500</td>
<td align="center">0.25</td>
<td align="center">
<xref ref-type="bibr" rid="B19">Goel et al. (1995)</xref>
</td>
</tr>
<tr>
<td align="center">Sacrum</td>
<td align="center">5000</td>
<td align="center">0.2</td>
<td align="center">
<xref ref-type="bibr" rid="B7">Du et al. (2014)</xref>
</td>
</tr>
<tr>
<td align="center">Endplate</td>
<td align="center">23.8</td>
<td align="center">0.4</td>
<td align="center">
<xref ref-type="bibr" rid="B31">Lu et al. (1996)</xref>
</td>
</tr>
<tr>
<td align="center">Ground substance</td>
<td colspan="2" align="center">Mooney-Rivlin C10 &#x3d; 0.42, C01 &#x3d; 0.105</td>
<td align="center">
<xref ref-type="bibr" rid="B29">Liu et al. (2011)</xref>
</td>
</tr>
<tr>
<td align="center">Nucleus pulposus</td>
<td colspan="2" align="center">Mooney-Rivlin C10 &#x3d; 0.12, C01 &#x3d; 0.03</td>
<td align="center">
<xref ref-type="bibr" rid="B7">Du et al. (2014)</xref>
</td>
</tr>
<tr>
<td align="center">Ligaments</td>
<td colspan="2" align="center">Nonlinear stress-strain curve</td>
<td align="center">
<xref ref-type="bibr" rid="B38">Rohlmann et al. (2006)</xref>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Properties of the ligaments (<xref ref-type="bibr" rid="B38">Rohlmann et al., 2006</xref>).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Ligament</th>
<th align="center">Stiffness (K1) (N/mm)</th>
<th align="center">Strains (&#x3b5;1) (&#x2212;)</th>
<th align="center">Stiffness (K2) (N/mm)</th>
<th align="center">Strains (&#x3b5;2) (&#x2212;)</th>
<th align="center">Stiffness (K3) (N/mm)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">ALL</td>
<td align="center">347</td>
<td align="center">0.122</td>
<td align="center">787</td>
<td align="center">0.203</td>
<td align="center">1864</td>
</tr>
<tr>
<td align="center">PLL</td>
<td align="center">29.5</td>
<td align="center">0.111</td>
<td align="center">61.7</td>
<td align="center">0.230</td>
<td align="center">236</td>
</tr>
<tr>
<td align="center">ISL</td>
<td align="center">1.4</td>
<td align="center">0.139</td>
<td align="center">1.5</td>
<td align="center">0.200</td>
<td align="center">14.7</td>
</tr>
<tr>
<td align="center">SSL</td>
<td align="center">2.5</td>
<td align="center">0.200</td>
<td align="center">5.3</td>
<td align="center">0.250</td>
<td align="center">34</td>
</tr>
<tr>
<td align="center">ITL</td>
<td align="center">0.3</td>
<td align="center">0.182</td>
<td align="center">1.8</td>
<td align="center">0.233</td>
<td align="center">10.7</td>
</tr>
<tr>
<td align="center">LF</td>
<td align="center">7.7</td>
<td align="center">0.059</td>
<td align="center">9.6</td>
<td align="center">0.490</td>
<td align="center">58.2</td>
</tr>
<tr>
<td align="center">CL</td>
<td align="center">36</td>
<td align="center">0.250</td>
<td align="center">159</td>
<td align="center">0.300</td>
<td align="center">384</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>To validate the FE model constructed in this study, we compared and analyzed the ROM of L1&#x2013;S1 under the following conditions: flexion of 8 Nm, extension of 6 Nm, lateral bending of &#xb1;6 Nm, and axial torsion of &#xb1;4 Nm, which were the same as the spine motion in previous studies. We also compared and analyzed the IDP under a compressive follower load and the height of the intervertebral disc between the L4&#x2013;L5 vertebrae at a compressive follower load of 1,200&#xa0;N (<xref ref-type="bibr" rid="B37">Renner et al., 2007</xref>; <xref ref-type="bibr" rid="B6">Dreischarf et al., 2014</xref>).</p>
<p>The location and number of degenerative segments of facet joints were the main variables in this study and were represented by L4/L5 FJD (Model1; M1) and L4/S1 FJD (Model2; M2), respectively. The degeneration grade of the facet joint was classified into five stages, including the intact state, based on the severity of degeneration. The gap and friction coefficient of the contacting surface of the facet joint were assumed to be graded based on the stage of degeneration, where Grade 0 corresponds to 100% and Grade 4 corresponds to 10%, with intermediate grades of 75, 50, and 25% corresponding to Grade 1, 2, and 3, respectively (<xref ref-type="table" rid="T3">Table 3</xref>). Finally, the gap in the contacting surface of the facet joint in L4&#x2013;L5 FSU and L5&#x2013;S1 FSU according to the degenerative location and degree of the facet joint is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Primary variables associated with lumbar facet joint degeneration and the values of each primary variable based on the degree of degeneration.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Degree of FJD</th>
<th colspan="2" align="center">Lumbar facet joint contact parameters</th>
</tr>
<tr>
<th align="center">Facet joint gap size</th>
<th align="center">Friction coefficient</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Grade0 (G0)</td>
<td align="right">0.5&#xa0;mm (100%)</td>
<td align="right">0 (100%)</td>
</tr>
<tr>
<td align="center">Grade1 (G1)</td>
<td align="right">0.375&#xa0;mm (75%)</td>
<td align="right">0.25 (75%)</td>
</tr>
<tr>
<td align="center">Grade2 (G2)</td>
<td align="right">0.25&#xa0;mm (50%)</td>
<td align="right">0.5 (50%)</td>
</tr>
<tr>
<td align="center">Grade3 (G3)</td>
<td align="right">0.125&#xa0;mm (25%)</td>
<td align="right">0.75 (25%)</td>
</tr>
<tr>
<td align="center">Grade4 (G4)</td>
<td align="right">0.05&#xa0;mm (10%)</td>
<td align="right">0.9 (10%)</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Localization of facet joint degeneration and modified facet joint space narrowing method at L4&#x2013;L5 and L5&#x2013;S1 functional spinal units (FSUs).</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g002.tif"/>
</fig>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<sec id="s3-1">
<title>3.1 Validation of the finite element (FE) model</title>
<p>The FE model was validated before conducting the study. We compared the measured ROM during spine movements (flexion, extension, lateral bending, and axial torsion) in <xref ref-type="fig" rid="F3">Figures 3</xref>, <xref ref-type="fig" rid="F4">4</xref>, and checked whether they were within the range of results from previous studies. We also compared the IDP at the L4&#x2013;L5 FSU and the height of the intervertebral disc at the L4&#x2013;L5 FSU under follower load with the results of previous studies. The ROM and intervertebral disk height measured in the FE model developed in this study were compared with the results of the study by <xref ref-type="bibr" rid="B37">Renner et al. (2007)</xref> on <italic>ex vivo</italic> spinal structures and spinal FE model behavior simulations. We confirmed that the ROM of the entire segment (L1&#x2013;S1) and each segment measured in this study were within the range of results reported by Renner et al. Additionally, when comparing IDP at the L4&#x2013;L5 FSU, our model was also within the range of IDP from seven previous studies on spine motion simulations and <italic>ex vivo</italic> spinal structure experiments (<xref ref-type="bibr" rid="B6">Dreischarf et al., 2014</xref>). Based on this, we concluded that the lumbar L1&#x2013;S1 FE model developed in this study can be utilized.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Comparison of the range of motion (ROM) based on spinal motion between the model used in this study and that of a study by <xref ref-type="bibr" rid="B37">Renner et al. (2007)</xref>. <bold>(A)</bold> Flexion &#x2b; Extension, <bold>(B)</bold> Right&#x2013;Left lateral bending, <bold>(C)</bold> Right&#x2013;Left axial torsion, <bold>(D)</bold> Compression.</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g003.tif"/>
</fig>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Comparison of intradiscal pressure (IDP) findings at the L4&#x2013;L5 functional spinal unit (FSU) in this study and those in previous studies (<xref ref-type="bibr" rid="B6">Dreischarf et al., 2014</xref>).</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g004.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>3.2 Range of motion (ROM)</title>
<p>
<xref ref-type="fig" rid="F5">Figures 5</xref>, <xref ref-type="fig" rid="F6">6</xref> show the results of the lumbar segmental ROM simulations performed under the loading conditions described in <xref ref-type="table" rid="T4">Table 4</xref>. Specifically, <xref ref-type="fig" rid="F5">Figure 5</xref> shows the total ROM of lumbar L1&#x2013;S1, and <xref ref-type="fig" rid="F6">Figure 6</xref> shows the ROM of the degenerated and adjacent segments of the facet joint.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Comparison of the total range of motion (ROM) based on the degree of lumbar facet joint degeneration. <bold>(A)</bold> Flexion, <bold>(B)</bold> Extension, <bold>(C)</bold> Lateral bending, <bold>(D)</bold> Axial torsion.</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Comparison of the range of motion (ROM) in onset and onset-adjacent segments (L3&#x2013;L4, L4&#x2013;L5, and L5&#x2013;S1) according to the grade of lumbar facet joint degeneration. <bold>(A)</bold> Flexion, <bold>(B)</bold> Extension, <bold>(C)</bold> Lateral bending, <bold>(D)</bold> Axial torsion.</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g006.tif"/>
</fig>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Loading conditions applied to L1 vertebra in this study.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Spinal motion</th>
<th align="center">Flexion</th>
<th align="center">Extension</th>
<th align="center">Lateral bending</th>
<th align="center">Axial torsion</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Moment (Nm)</td>
<td align="center">7.5</td>
<td align="center">7.5</td>
<td align="center">7.5</td>
<td align="center">7.5</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>As shown in <xref ref-type="fig" rid="F5">Figure 5</xref>, as the degeneration of the facet joint worsened in all spine motions, the ROM of the L1&#x2013;S1 vertebrae gradually decreased. Furthermore, when degeneration occurred in two segments (L4/S1 FJD; M2) rather than one segment (L4/L5 FJD; M1), the ROM of the L1&#x2013;S1 vertebrae decreased. During flexion movements, Model 1 did not differ significantly from Model 0 in terms of ROM, even as the degeneration grade increased. Model 2 decreased by approximately 1.2%&#x2013;4.5% compared to Model 0 (<xref ref-type="fig" rid="F5">Figure 5A</xref>). During extension movements, both Models 1 and 2 decreased the ROM compared to Model 0, while Model 2 showed a maximum decrease of 11.9% (Grade 4) (<xref ref-type="fig" rid="F5">Figure 5B</xref>). During lateral bending movements, both Models 1 and 2 decreased the ROM compared to Model 0, whereas Model 2 decreased by 6.4% (Grade 4) compared to Model 0 (<xref ref-type="fig" rid="F5">Figure 5C</xref>). Moreover, during axial torsion movements, Model 1 increased by 1.1%&#x2013;8.9% in Grade 1, and Model 3 increased by Grade 1&#x2013;2 compared to Model 0. However, as the degree of FJD increased, the ROM of Model 0 decreased compared with that of Model 1 in Grade 1 and Model 3 in Grade 1&#x2013;2 (<xref ref-type="fig" rid="F5">Figure 5D</xref>).</p>
<p>The ROM of each segment (L4&#x2013;L5 FSU, L5&#x2013;S1 FSU) decreased gradually as the degeneration of the facet joint worsened in all spinal motions, as shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. However, the ROM of each segment (L3&#x2013;L4 FSU, L5&#x2013;S1 FSU) in adjacent segments was not significantly different from that of Model 0 in an intact state. In the initial analysis of the ROM for each segment (L4&#x2013;L5 FSU, L5&#x2013;S1 FSU) in the onset segment, the ROMs of Models 1 and 2 remained nearly identical during the flexion movement within the L4&#x2013;L5 FSU. Furthermore, the ROM did not change significantly even as the grade of degeneration increased (<xref ref-type="fig" rid="F6">Figure 6A</xref>). The ROM difference between the two models was almost negligible during other movements, such as extension, lateral bending, and axial torsion; however, the ROM gradually decreased depending on the degeneration grade, unlike during flexion movement. Compared to the ROM of Model 0, the ROM reduction rates of the two models were 14.62%&#x2013;39.28% (M1G1&#x2013;M1G4) and 14.59%&#x2013;39.51% (M2G1&#x2013;M2G4) during extension, and 5%&#x2013;18.81% (M1G1&#x2013;M1G4) and 5%&#x2013;19.71% (M2G1&#x2013;M2G4) during lateral bending (<xref ref-type="fig" rid="F6">Figures 6B,C</xref>). Furthermore, during axial torsion, the ROM slightly increased in Grade 1 but showed a reduction rate of 26.43% (M1G4) and 28.33% (M2G4) in Grade 4 (<xref ref-type="fig" rid="F6">Figure 6D</xref>). In the L5&#x2013;S1 FSU as shown in <xref ref-type="fig" rid="F6">Figure 6</xref>, Model 2 showed a gradual decrease in ROM depending on the grade of degeneration in all spine motions (flexion, extension, lateral bending, and axial torsion), and the maximum reduction rates in each movement observed in Grade 4 were 22.24 (flexion), 15.43 (extension), 21.53 (lateral bending), and 12.87% (axial torsion). Finally, regarding the ROM of each segment (L3&#x2013;L4 FSU, L5&#x2013;S1 FSU) in adjacent segments, the ROMs of Models 1 and 2 showed the same values, regardless of the degeneration grade, and were almost identical to the ROM of Model 0.</p>
</sec>
<sec id="s3-3">
<title>3.3 Maximum von mises stress of the intervertebral disc</title>
<p>
<xref ref-type="fig" rid="F7">Figures 7</xref>, <xref ref-type="fig" rid="F8">8</xref> show the maximum von Mises stress results of the simulations performed under the load conditions listed in <xref ref-type="table" rid="T4">Table 4</xref> for the intervertebral discs. Specifically, <xref ref-type="fig" rid="F7">Figure 7</xref> shows the maximum von Mises stress of the intervertebral disc at the onset segment, whereas <xref ref-type="fig" rid="F8">Figure 8</xref> shows the maximum von Mises stress of the intervertebral disc in the adjacent upper and lower segments at the onset; both are presented as a percentage change from the intact state (M0).</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Comparison of the maximum von Mises stress in intervertebral discs at the onset segment based on the degree of lumbar facet joint degeneration and spinal motion (flexion, extension, lateral bending, and axial torsion). <bold>(A)</bold> Flexion, <bold>(B)</bold> Extension, <bold>(C)</bold> Lateral bending, <bold>(D)</bold> Axial torsion.</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g007.tif"/>
</fig>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Comparative analysis of the percentage changes in maximum von Mises stress at the intervertebral discs of adjacent levels compared to that of the intact state, based on the degree of lumbar facet joint degeneration and spine motion (flexion, extension, lateral bending, and axial torsion). <bold>(A)</bold> Flexion, <bold>(B)</bold> Extension, <bold>(C)</bold> Lateral bending, <bold>(D)</bold> Axial torsion. Percentage change &#x3d; (data of Model1 or Model2&#x2013;data of Model0)/data of Model0&#xd7;100%.</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g008.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F7">Figure 7</xref>, the maximum von Mises stress on the intervertebral disc decreased as the degeneration of the facet joint increased in all spinal motions. Moreover, in Model 2, which comprises two segments, there was a notably higher maximum von Mises stress observed on the L5&#x2013;S1 FSU during axial rotation when compared to that of Model 0. However, the trend of gradually decreasing maximum von Mises stress on the intervertebral disc according to the degeneration grade was consistent. Firstly, regarding the maximum von Mises stress of the intervertebral disc at onset, it was observed that Model 1 and Model 2 increased slightly or remained the same as Model 0 during flexion movement in L4&#x2013;L5 FSU, and that Model 2 increased by 4.04 kPa&#x2013;520.9&#xa0;kPa (M2G4) in Grade 4. Furthermore, regarding the L5&#x2013;S1 FSU, Model 2 gradually decreased as degeneration increased, showing a decrease rate of 0.74%&#x2013;8.23% (M2G1&#x2013;M2G4) compared to Model 0 (<xref ref-type="fig" rid="F7">Figure 7A</xref>). Furthermore, considering extension, lateral bending, and axial torsion movements (<xref ref-type="fig" rid="F7">Figures 7B&#x2013;D</xref>), when compared to Model 0, the decreased rate in Model 2 in the L4&#x2013;L5 FSU was 21.29%&#x2013;54.60% (M1G1&#x2013;M1G4) during extension, 21.18%&#x2013;54.72% (M2G1&#x2013;M2G4) during lateral bending, and 9.34%&#x2013;32.49% (M1G1&#x2013;M1G4) and 11.10%&#x2013;31.95% (M2G1&#x2013;M2G4) during axial torsion. In the L5&#x2013;S1 FSU, the decreased rate in Model 2 was 12.69%&#x2013;34.75% (M2G1&#x2013;M2G4) during extension and 0.00%&#x2013;16.15% (M2G1&#x2013;M2G4) during lateral bending. However, it increased compared to Model 0 during axial torsion, showing an increase rate of 27%&#x2013;76% (M2G4&#x2013;M2G1).</p>
<p>As shown in <xref ref-type="fig" rid="F8">Figure 8</xref>, the percentage change in the maximum von Mises stress of the intervertebral disc in a single segment adjacent to the affected area gradually increased during flexion and extension movements as the degeneration of the facet joint worsened compared to the intact state (M0) (<xref ref-type="fig" rid="F8">Figures 8A, B</xref>). Additionally, the percentage change in the maximum von Mises stress of the intervertebral disc gradually decreased during lateral bending and axial torsion movements (<xref ref-type="fig" rid="F8">Figures 8C, D</xref>). First, when observing the adjacent upper segment, there was no change in Model 1 as the grade of degeneration of the facet joint increased during flexion movements; furthermore, in Model 2, it only slightly increased, with a maximum increase rate of 0.18% (M2G4), which was not statistically significant (<xref ref-type="fig" rid="F8">Figure 8A</xref>). During extension, the maximum increase rates of Model 1 and Model 2 were slightly increased to 1.16% (M1G4) and 1.18% (M2G4), respectively (<xref ref-type="fig" rid="F8">Figure 8B</xref>). In the lateral bending and axial torsion movements, the maximum decrease rates of Model 1 and Model 2 were 0.46% (M1G4) and 1.00% (M2G4) for lateral bending and 2.95% (M1G4) and 2.25% (M2G4) for axial torsion, respectively (<xref ref-type="fig" rid="F8">Figures 8C,D</xref>). Furthermore, when observing the adjacent lower segment, which was the L5&#x2013;S1 FSU, only Model 1 was observed. During the extension movements, it increased to approximately 4.00% up to Grade 3, depending on the degeneration grade of the facet joint, and then decreased to approximately 3.34% at Grade 4 (<xref ref-type="fig" rid="F8">Figure 8B</xref>). The maximum increase and decrease rates during the lateral bending and axial torsion movements were 0.90% (M1G4) and 0.65% (M1G4), respectively, depending on the degeneration grade of the facet joint (<xref ref-type="fig" rid="F8">Figures 8C,D</xref>).</p>
</sec>
<sec id="s3-4">
<title>3.4 Reaction force on the facet joint</title>
<p>
<xref ref-type="fig" rid="F9">Figures 9</xref>, <xref ref-type="fig" rid="F10">10</xref> present the reaction force results of the simulations performed under the load conditions of the facet joints listed in <xref ref-type="table" rid="T4">Table 4</xref>. Specifically, <xref ref-type="fig" rid="F9">Figure 9</xref> shows the reaction force of the facet joint at the onset segment, whereas <xref ref-type="fig" rid="F10">Figure 10</xref> shows the reaction force of the facet joint in the adjacent upper and lower segments at the onset. Both are presented as a percentage change from the intact state (M0). However, flexion movements that did not cause a facet joint reaction force because the articular surfaces of the zygapophyseal joints that did not come into contact during spinal motion were excluded (Cai et al.). In addition, only the right direction was implemented for lateral bending and axial torsion movements, and the right zygapophyseal joint was excluded when the gap between the zygapophyseal joints widened during right lateral bending and right axial torsion movements.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Comparison of the reaction force in the facet joint at the onset segment based on the degree of lumbar facet joint degeneration and spine motion. <bold>(A)</bold> Extension, <bold>(B)</bold> Lateral bending, <bold>(C)</bold> Axial torsion.</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g009.tif"/>
</fig>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Comparative analysis of the percentage changes in the reaction force at the facet joints of adjacent levels compared to the intact state based on the degree of lumbar facet joint degeneration and spine motion. <bold>(A)</bold> Right of facet joint during extension, <bold>(B)</bold> Left of facet joint during extension, <bold>(C)</bold> Right of facet joint during lateral bending, <bold>(D)</bold> Left of facet joint during axial torsion. Percentage change &#x3d; (data of Model1 or Model2&#x2013;data of Model0)/data of Model0&#xd7;100%.</p>
</caption>
<graphic xlink:href="fbioe-12-1294658-g010.tif"/>
</fig>
<p>First, regarding the facet joint reaction force in the onset segment, the forces of both the left and right facet joints during spinal motion showed a similar trend depending on the degeneration grade of the facet joint. During the simulated extension movements (<xref ref-type="fig" rid="F9">Figure 9A</xref>), the facet joint reaction force in both Models 1 and 2 decreased and then increased depending on the degeneration of the onset segments. In Model 1, at L4&#x2013;L5 FSU, the facet joint reaction force decreased from 69.496 N to 66.68&#xa0;N (M1G1) to 63.398 N and 61.585&#xa0;N (M1G3), and increased to 64.389 N and 62.263&#xa0;N (M1G4). In Model 2, at the L4&#x2013;L5 FSU, the facet joint reaction force decreased from 69.566 N to 67.343&#xa0;N (M2G1) to 63.887 N and 62.378&#xa0;N (M2G3), and increased to 64.623 N and 63.028&#xa0;N (M2G4), while at the L5&#x2013;S1 FSU, it decreased from 51.152 N to 46.467&#xa0;N (M2G1) to 33.191 N and 28.509&#xa0;N (M2G2) and increased to 39.927 N and 34.801&#xa0;N (M2G4). During the simulated lateral bending movements (<xref ref-type="fig" rid="F9">Figure 9B</xref>), the facet joint reaction force showed an increasing trend depending on the degeneration of the onset segments in the L4&#x2013;L5 FSU and a decreasing trend in the L5&#x2013;S1 FSU. In Model 1, at the L4&#x2013;L5 FSU, the facet joint reaction force gradually increased from 72.516&#xa0;N (M1G1) to 85.46&#xa0;N (M1G4). In Model 2, at L4&#x2013;L5 FSU, the facet joint reaction force gradually increased from 72.256&#xa0;N (M2G1) to 87.673&#xa0;N (M2G4), and at L5&#x2013;S1 FSU, it increased from 119.93&#xa0;N (M2G1) to 142.26&#xa0;N (M2G2) and then decreased to 125.1&#xa0;N (M2G4). Owing to the degeneration of the facet joint during the axial torsion movement (<xref ref-type="fig" rid="F9">Figure 9C</xref>), the facet joint reaction force gradually increased in the L4&#x2013;L5 and L5&#x2013;S1 FSUs. Model 1 gradually increases from 73.83&#xa0;N (M1G1) to 112.81&#xa0;N (M1G4) in the L4&#x2013;L5 FSU. Model 2 showed an increase from 75.189&#xa0;N (M2G1) to 111.67&#xa0;N (M2G4) in the L4&#x2013;L5 FSU and an increase from 71.47&#xa0;N (M2G1) to 98.564&#xa0;N (M2G4) in the L5&#x2013;S1 FSU.</p>
<p>We confirmed that the percentage change in the facet joint reaction force gradually increased during all spinal motions compared with that of the intact facet joint (M0) in the adjacent upper and lower segments where the onset occurred (<xref ref-type="fig" rid="F10">Figure 10</xref>). Moreover, in the case of onset in both segments (M2), the increase was greater than that in the case of onset in one segment (M1). During extension movements (<xref ref-type="fig" rid="F10">Figures 10A, B</xref>), the percentage change in the reaction force compared to the intact (M0) left and right facet joints in the adjacent upper segment gradually increased as the degeneration grade of the facet joint increased, reaching a maximum of 1.82% and 1.56% (M1G4), 0.48% and 0.53% (M2G4), and 1.81% and 1.78% (M3G4). In addition, during adjacent lower segment movements the percentage change in the reaction force of the left and right facet joints gradually increased, similar to that of the adjacent upper segment, ranging from 2.9% to 4.52% and 0.8%&#x2013;3.39% (M1G1&#x2013;M1G4). During lateral bending movements (<xref ref-type="fig" rid="F10">Figure 10C</xref>), the percentage change in reaction force compared with that of the intact facet joint (M0) in the adjacent upper segment gradually increased as the degeneration grade of the facet joint increased, reaching a maximum of 3.7% (M1G4) and 6.5% (M2G4). In addition, even in the facet joint of the adjacent lower segment, a maximum of 1.45% (M1G4) was observed. During axial torsion movements (<xref ref-type="fig" rid="F10">Figure 10D</xref>), the percentage change in the reaction force compared with that of the intact (M0) facet joint in the adjacent upper segment gradually increased as the degeneration grade of the facet joint increased, reaching a maximum of 2.43% (M1G4) and 3.99% (M2G4). The change in percentage compared to the intact right facet joint of the adjacent lower segment was minimal, with a maximum of 0.35% (M1G4).</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<p>In the present study, the biomechanical effects of FJD on the lumbar spine were investigated through finite element analysis. Finite element models with varying degrees of FJD and numbers of onset segments were fabricated to simulate four common physiological loads experienced during daily activities. The degree of FJD affected the ROM of the lumbar spine, von Mises stress of the intervertebral discs, and facet joint reaction force, demonstrating how FJD affects onset and onset-adjacent segments. As the grade of FJD increased during spinal motion, the ROM gradually decreased, particularly in the onset segments. Model 2, where FJD occurred in two segments, showed a much greater decrease in ROM than did Model 1. Similarly, in each segment with FJD, the decrease in ROM was similar to or slightly greater in Model 2 than in Model 1, demonstrating the influence of the number of onset segments on the decrease in ROM. However, the number of onset segments did not have a significant effect on the maximum von Mises stress of the intervertebral discs or facet joint reaction force.</p>
<p>In a recent study (<xref ref-type="bibr" rid="B49">Yin et al., 2020</xref>), medical images of patients with degenerative lumbar facet joint disease who experienced spine motion after maximum bending of the body during flexion-extension, lateral bending, and axial torsion were captured using a dual fluoroscopic imaging system. Based on this, a 3D model of the L3&#x2013;S1 spinal segment was developed, and the ROM of each was measured. As a result, the ROM was found to be highest in the moderate stage of spinal motion but decreased in the severe stage. These results differ from the tendency of ROM to gradually decrease with FJD, as reported in previous studies. This is attributed to difficulty in completely controlling for variables, such as the possibility that patients with degenerative facet joint disease may also experience degeneration of the IVDs in previous studies. However, both studies showed that ROM decreases during spinal motion in patients with severe FJD, indicating that degenerative facet joint disease limits spinal motion.</p>
<p>In a study conducted by <xref ref-type="bibr" rid="B15">Fujiwara et al. (2000)</xref>, intervertebral discs and facet joints of patients were graded according to the degeneration, and their ROM was compared based on the motion of their spines. Specifically, degenerative changes in the facet joints were categorized into cartilage degeneration, osteophyte formation, and osteoarthritis, each of which was further divided into four grades to compare the ROM. Osteoarthritis, one of the three symptoms of facet joint degeneration, exhibits a gradual decrease in ROM as degeneration worsens, which is similar to the findings of the current study. Since the main parameters of this study were assumed to be the gap and friction coefficient of the posterior joint, we believe that the results are similar to those of osteoarthritis, where the articular surface of the facet joint becomes rougher as the joint space narrows (<xref ref-type="bibr" rid="B21">Jarraya et al., 2018</xref>).</p>
<p>In addition, it has been reported that degenerative changes in the facet joint due to osteoarthritis can cause friction between the bones of the facet joint (<xref ref-type="bibr" rid="B21">Jarraya et al., 2018</xref>; <xref ref-type="bibr" rid="B50">Anastasia et al., 2022</xref>). This means that as the grade of degeneration of the facet joint increases, a higher coefficient of friction can occur between the bones. However, it has been difficult to formalize this. In this study, we assumed four grades of degeneration and a constant coefficient of friction for each grade. This way, we aimed to emulate the interosseous friction that occurs in clinical degeneration of the facet joint. </p>
<p>Spinal degeneration occurs naturally with aging and is further accelerated by mechanical weakening of tissues, such as ligaments, due to repetitive injuries or trauma (<xref ref-type="bibr" rid="B22">Jaumard et al., 2011</xref>; <xref ref-type="bibr" rid="B46">Wang et al., 2019</xref>). This degenerative process is characterized by temporary functional impairment, instability, and restabilization (<xref ref-type="bibr" rid="B27">Kirkaldy-Willis and Farfan, 1982</xref>). Additionally, it has been reported that the role of facet joints changes as lumbar degeneration progresses (<xref ref-type="bibr" rid="B26">Kim et al., 2019</xref>), leading to degenerative changes in these joints due to excessive load, which increases lumbar instability. The role of facet joints is considered to be as critical as that of intervertebral discs in lumbar degeneration. Biomechanically, facet joints also play a role in limiting axial torsion and are involved in rotational kinematic mechanisms. Consequently, the stages of FJD are thought to be observable in the axial torsion behavior (<xref ref-type="bibr" rid="B1">Adams and Hutton, 1983</xref>; <xref ref-type="bibr" rid="B23">Kalichman and Hunter, 2007</xref>).</p>
<p>In this study, we observed instability and restabilization of the lumbar spine during axial torsion. Specifically, in the case of the L5&#x2013;S1 FSU, which is the onset segment, Model 2 exhibited a higher ROM than Model 0 in Grade 1, indicating lumbar instability. As Grade 4 was reached, the ROM decreased compared to that of Model 0, indicating lumbar restabilization. We confirmed that as the ROM of the onset segment decreased during the restabilization stage, the load on the intervertebral disc caused by the onset segment was significantly reduced. Notably, the average decrease rate in the L4&#x2013;L5 FSU of both models showed a reduction of 54% (flexion), 35% (lateral bending), and 32% (axial torsion) in Grade 4 compared to those of the intact state. However, in the facet joints of the onset segment, the load gradually increased as the degeneration grade increased. For Grade 4, the average increase rates of the two models were 12.5% (lateral bending) and 26% (axial torsion) compared with the intact state. Through this research, it is possible to quantitatively check the information that is difficult to check in medical images, such as the stress on the intervertebral disc and the reaction force of the facet joint, in addition to the information that can be checked in medical images during the development of posterior joint degeneration. Therefore, we believe that the results of this research, such as the stress on the intervertebral disc of the adjacent segment and the reaction force of the facet joint, can be used as a predictor of the development of potential lumbar degenerative diseases.</p>
<p>The ROM of a healthy adjacent segment increases to compensate for the decreased ROM of a degenerated segment (<xref ref-type="bibr" rid="B35">Park et al., 2015</xref>). <xref ref-type="bibr" rid="B47">Wang et al. (2018)</xref>, introduced anterior vertebral osteophytes and decreased disc height, demonstrating that as the degeneration grade increased, the ROM of the adjacent segment also increased. Similarly, <xref ref-type="bibr" rid="B35">Park et al. (2015)</xref> conducted a study on cervical segments, considering both facet joint and disc degeneration, and found that as the degeneration grade increased, the ROM of the adjacent segment gradually increased. However, in our study, we observed almost no change in the ROM of the adjacent segment, which differs from the findings of previous studies. This discrepancy is attributed to the exclusion of disc degeneration and height reduction because our study focused solely on FJD. Therefore, we believe that disc degeneration may have a more significant impact on ROM changes in the adjacent segment compared to FJD.</p>
<p>In this study, we found no significant differences in spinal ROM, the maximum von Mises stress on the intervertebral disc, and the reaction force of the facet joint values between intact, L4/5 FJD, and L4/S1 FJD for the upper segments corresponding to L1-L3 FSU. As a result, we did not include the upper segments in the results of the present study.</p>
<p>As mentioned earlier, although little change was observed in the ROM of the adjacent segment in this study, we observed an increase in the load on both the disc and facet joint of the adjacent segment. Notably, the facet joint had a greater impact in both segments. This finding aligns with the effects of disc or FJD in the adjacent segments due to spinal degeneration. Consequently, in clinical practice, if disc and FJD occur simultaneously, the impact on the adjacent segment is expected to be greater, potentially accelerating additional degeneration. Previous studies have predominantly focused on lumbar intervertebral discs in the context of spinal degeneration, resulting in numerous research efforts related to disc degeneration (<xref ref-type="bibr" rid="B36">Pfirrmann et al., 2001</xref>; <xref ref-type="bibr" rid="B40">Schmidt et al., 2009</xref>; <xref ref-type="bibr" rid="B14">Fan et al., 2015</xref>; <xref ref-type="bibr" rid="B10">Ellingson et al., 2016</xref>). However, considering that spinal load is distributed between the lumbar intervertebral discs and facet joints is important. Therefore, research on facet joints is necessary to comprehensively understand spinal degeneration. This study adopted a novel approach as it exclusively examined lumbar FJD.</p>
<p>Several simplifications and assumptions were made in this study. Specifically, the investigation focused solely on degenerative changes in facet joints during spinal degeneration. Although various shape changes related to FJD, such as osteophytes and asymmetry, can occur, this study employed the number of onset segments and the gap and friction coefficients of the facet joint as parameters to explore changes in lumbar characteristics across the four grades of degeneration. The location of the onset segment with FJD was assumed to be limited to the most commonly affected L4&#x2013;L5 FSU and L5&#x2013;S1 FSU. Moreover, real patients with FJD often experience disc height reduction; however, for the purposes of this study, it was assumed that the disc height remained constant regardless of the facet joint gap. This decision was made to exclusively investigate the impact of FJD. Therefore, this study offers insights into the biomechanical characteristics of the lumbar region by considering FJD. However, further observational studies are necessary to validate the clinical changes caused by isolated lumbar FJD, as examined in this research.</p>
<p>There are also limitations in the modeling of the lumbar spine and the material properties applied. First of all, the author&#x2019;s vertebrae are divided into cortical bone, cancellous bone, and posterior body, and each is assumed to have simple isotropic linear elastic material properties. However, when we look at the cervical vertebrae, we see that each part has different bone density and consequently different mechanical properties (<xref ref-type="bibr" rid="B17">Garay et al., 2022</xref>). And we see a non-uniform distribution of bone density throughout the vertebrae (<xref ref-type="bibr" rid="B2">Al-Barghouthi et al., 2020</xref>).</p>
<p>In addition, the cartilage and intervertebral discs of the facet joints are limited by the approach of implementing the facet joints under contact conditions among the various facet joint implementations in this study. However, the actual cartilage of the facet joints is a porous elastic model, and its mechanical properties change as the water content decreases and become stiffer during facet joint degeneration (<xref ref-type="bibr" rid="B11">Elmasry et al., 2017</xref>; <xref ref-type="bibr" rid="B12">Elmasry et al., 2018</xref>). The intervertebral disc model developed in this study is a hyperelastic material model, which has the limitation of simplicity. In recent years, various modeling methods have been used in intervertebral disc studies. One of them is an analysis using a porous elastic model (<xref ref-type="bibr" rid="B45">Volz et al., 2022</xref>). The porous elastic model models the biomechanical properties of intervertebral discs by considering the porosity and fluid flow in the disc. The model takes into account the interaction of liquids and solids inside the intervertebral disc and helps to understand the mechanisms of intervertebral disc degeneration. Cappetti et al. demonstrate the strong influence of geometric parameters in intervertebral disc modeling (<xref ref-type="bibr" rid="B5">Cappetti et al., 2016</xref>). Geometric parameters include the radius, thickness, status, curvature, position, and orientation of the intervertebral disc, which are said to affect the biomechanical properties of the intervertebral disc. To evaluate how each parameter affects the output, they used a sensitivity analysis method using Taguchi Orthogonal Array, which allows them to consider many parameters while minimizing the number of experiments, and quickly identify the parameters that have a significant impact on the output. This allows us to improve the accuracy of our modeling with less time and cost. In future research, we believe that by applying the advanced techniques in lumbar spine modeling and material properties mentioned above, we can simulate spinal behavior closer to the human body.</p>
<p>Another factor that was not monitored in this study but contributes significantly to pain is the disc bulging (<xref ref-type="bibr" rid="B3">Amirouche et al., 2015</xref>). Amirouche et al. describe segmental stiffness due to disc degeneration and the resulting degree of the disc bulging in cervical spine subjects. Further studies should consider not only facet joint degeneration but also intervertebral disc degeneration, so that the degree of the disc bulging can be utilized as a measure of spinal pain.</p>
<p>Among the various methods used to implement the FJs in this study, contact conditions were utilized (<xref ref-type="bibr" rid="B32">Mengoni, 2021</xref>). However, contact conditions cannot be applied when the contact surfaces are not in contact, such as during flexion (<xref ref-type="bibr" rid="B4">Cai et al., 2020</xref>). Additionally, this study did not consider changes in the properties of capsular ligament during FJD. Consequently, it does not have a significant impact on the ROM and intervertebral disc. Therefore, the FJD implementation method in this study was found to have no significant impact on the ROM and biomechanical characteristics of the intervertebral disc in flexion due to its simplification. In this study, four levels of degeneration were arbitrarily set by adjusting the gap and friction coefficient of the facet joint to 75, 50, 25, and 10% of the intact state (100%) for implementing FJD. However, considering that spinal load is distributed between the lumbar intervertebral discs and facet joints is important. Therefore, research on facet joints is necessary to comprehensively understand spinal degeneration. Recently, studies have been reported on the effects of facet joint parameters such as facet orientation (FO) and facet tropism (FT) on the lumbar spine (<xref ref-type="bibr" rid="B54">Teo et al., 2003</xref>; <xref ref-type="bibr" rid="B52">Kim et al., 2013</xref>; <xref ref-type="bibr" rid="B53">Li et al., 2020</xref>; <xref ref-type="bibr" rid="B51">Ke et al., 2021</xref>). The left and right facet joint angles of the vertebral body in the sagittal plane are referred to as FO, and the difference between the left and right angles is FT. The biomechanical effects of FT and FO on lumbar segmental stresses have been investigated, and it is reported that FT has the greatest effect on increasing intradiscal pressure and facet joint pressure (<xref ref-type="bibr" rid="B52">Kim et al., 2013</xref>; <xref ref-type="bibr" rid="B51">Ke et al., 2021</xref>). Studies investigating the association between FO and FT and recurrent lumbar disc herniation (rLDH) have found that FO and FT are associated with the development of rLDH, with decreased FO reported to be associated with an increased risk of rLDH and increased FT reported to be associated with an increased incidence of rLDH (<xref ref-type="bibr" rid="B53">Li et al., 2020</xref>). This suggests that facet joint parameters that were not considered in this study may also influence the lumbar spine. Therefore, this study is original in that it considered the parameter (gap size, friction coefficient) of lumbar FJD.</p>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>Overall, the results of this study showed that the decrease in the ROM of the onset segment due to FJD was caused by a decrease in the ROM of the entire segment and a decrease in the maximum von Mises stress of the intervertebral disc of the onset segment. However, the maximum von Mises stress on the IVD of the onset-adjacent segment during flexion and extension gradually increases as the level of degeneration becomes more severe. The facet joint reaction force decreased during the extension motion in the onset segment and then increased, whereas it increased during lateral bending and axial torsion. In contrast, the facet joint reaction force gradually increased in the adjacent segments.</p>
<p>The onset segment enters the stage of restabilization at Grade 4 degeneration and its instability is reduced further. However, the increased load on the intervertebral disc and facet joint in the adjacent segment may contribute to its instability. Therefore, inadequate neural control and compensatory lumbar muscles are thought to contribute to the degeneration of the adjacent segments in clinical scenarios.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>JP: Formal Analysis, Investigation, Methodology, Software, Validation, Writing&#x2013;original draft, Writing&#x2013;review and editing. TG: Conceptualization, Data curation, Funding acquisition, Resources, Validation, Writing&#x2013;original draft. JL: Funding acquisition, Project administration, Supervision, Writing&#x2013;review and editing. CL: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Validation, Visualization, Writing&#x2013;original draft, Writing&#x2013;review and editing.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the National Research Foundation of Korea (NRF) and funded by the Korean Government (Ministry of Science and ICT) (NRF-2021R1F1A1064056, 2021R1A2B5B01087440, and 2022R1F1A1066509). This research was supported by a grant (RS-2023-00215638) from Ministry of Food and Drug Safety in 2023.</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Adams</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Hutton</surname>
<given-names>W. C.</given-names>
</name>
</person-group> (<year>1983</year>). <article-title>The mechanical function of the lumbar apophyseal joints</article-title>. <source>Spine</source> <volume>8</volume> (<issue>3</issue>), <fpage>327</fpage>&#x2013;<lpage>330</lpage>. <pub-id pub-id-type="doi">10.1097/00007632-198304000-00017</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Al-Barghouthi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Solitro</surname>
<given-names>G. F.</given-names>
</name>
<name>
<surname>Latta</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Travascio</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Relationships among bone morphological parameters and mechanical properties of cadaveric human vertebral cancellous bone</article-title>. <source>JBMR plus</source> <volume>4</volume> (<issue>5</issue>), <fpage>e10351</fpage>. <pub-id pub-id-type="doi">10.1002/jbm4.10351</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Amirouche</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Solitro</surname>
<given-names>G. F.</given-names>
</name>
<name>
<surname>Siemionow</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Drucker</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Upadhyay</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Patel</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Role of posterior elements in the disc bulging of a degenerated cervical spine</article-title>. <source>Int. J. Spine Surg.</source> <volume>9</volume>, <fpage>13</fpage>. <pub-id pub-id-type="doi">10.14444/2013</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Anastasia</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Sukmaningtyas</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Priambodo</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Setiawati</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Analysis of risk factors affecting lumbal facet joint osteoarthritis in mri scan</article-title>. <source>Diponegoro Medical J. (Jurnal Kedokteran Diponegoro)</source> <volume>11</volume> (<issue>1</issue>), <fpage>19</fpage>&#x2013;<lpage>23</lpage>. <pub-id pub-id-type="doi">10.14710/dmj.v11i1.32588</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cai</surname>
<given-names>X.-Y.</given-names>
</name>
<name>
<surname>Sang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Yuchi</surname>
<given-names>C.-X.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>C.-F.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Using finite element analysis to determine effects of the motion loading method on facet joint forces after cervical disc degeneration</article-title>. <source>Comput. Biol. Med.</source> <volume>116</volume>, <fpage>103519</fpage>. <pub-id pub-id-type="doi">10.1016/j.compbiomed.2019.103519</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cappetti</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Naddeo</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Naddeo</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Solitro</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Finite elements/Taguchi method based procedure for the identification of the geometrical parameters significantly affecting the biomechanical behavior of a lumbar disc</article-title>. <source>Comput. methods biomechanics Biomed. Eng.</source> <volume>19</volume> (<issue>12</issue>), <fpage>1278</fpage>&#x2013;<lpage>1285</lpage>. <pub-id pub-id-type="doi">10.1080/10255842.2015.1128529</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dreischarf</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zander</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Shirazi-Adl</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Puttlitz</surname>
<given-names>C. M.</given-names>
</name>
<name>
<surname>Adam</surname>
<given-names>C. J.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>C. S.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Comparison of eight published static finite element models of the intact lumbar spine: predictive power of models improves when combined together</article-title>. <source>J. Biomechanics</source> <volume>47</volume> (<issue>8</issue>), <fpage>1757</fpage>&#x2013;<lpage>1766</lpage>. <pub-id pub-id-type="doi">10.1016/j.jbiomech.2014.04.002</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Du</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Mo</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Biomechanical investigation of thoracolumbar spine in different postures during ejection using a combined finite element and multi&#x2010;body approach</article-title>. <source>Int. J. Numer. Methods Biomed. Eng.</source> <volume>30</volume> (<issue>11</issue>), <fpage>1121</fpage>&#x2013;<lpage>1131</lpage>. <pub-id pub-id-type="doi">10.1002/cnm.2647</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Du</surname>
<given-names>C.-F.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>J.-C.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>Y.-P.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Biomechanical response of lumbar facet joints under follower preload: a finite element study</article-title>. <source>BMC Musculoskelet. Disord.</source> <volume>17</volume> (<issue>1</issue>), <fpage>126</fpage>. <pub-id pub-id-type="doi">10.1186/s12891-016-0980-4</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dunlop</surname>
<given-names>R. B.</given-names>
</name>
<name>
<surname>Adams</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Hutton</surname>
<given-names>W. C.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>Disc space narrowing and the lumbar facet joints</article-title>. <source>J. bone Jt. Surg. Br. volume</source> <volume>66</volume> (<issue>5</issue>), <fpage>706</fpage>&#x2013;<lpage>710</lpage>. <pub-id pub-id-type="doi">10.1302/0301-620x.66b5.6501365</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ellingson</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Shaw</surname>
<given-names>M. N.</given-names>
</name>
<name>
<surname>Giambini</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>K.-N.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Comparative role of disc degeneration and ligament failure on functional mechanics of the lumbar spine</article-title>. <source>Comput. Methods Biomechanics Biomed. Eng.</source> <volume>19</volume> (<issue>9</issue>), <fpage>1009</fpage>&#x2013;<lpage>1018</lpage>. <pub-id pub-id-type="doi">10.1080/10255842.2015.1088524</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elmasry</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Asfour</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Travascio</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Effectiveness of pedicle screw inclusion at the fracture level in short-segment fixation constructs for the treatment of thoracolumbar burst fractures: a computational biomechanics analysis</article-title>. <source>Comput. methods biomechanics Biomed. Eng.</source> <volume>20</volume> (<issue>13</issue>), <fpage>1412</fpage>&#x2013;<lpage>1420</lpage>. <pub-id pub-id-type="doi">10.1080/10255842.2017.1366995</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elmasry</surname>
<given-names>S. S.</given-names>
</name>
<name>
<surname>Asfour</surname>
<given-names>S. S.</given-names>
</name>
<name>
<surname>Travascio</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Finite element study to evaluate the biomechanical performance of the spine after augmenting percutaneous pedicle screw fixation with kyphoplasty in the treatment of burst fractures</article-title>. <source>J. biomechanical Eng.</source> <volume>140</volume> (<issue>6</issue>), <fpage>061005</fpage>. <pub-id pub-id-type="doi">10.1115/1.4039174</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eubanks</surname>
<given-names>J. D.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>M. J.</given-names>
</name>
<name>
<surname>Cassinelli</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Ahn</surname>
<given-names>N. U.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Prevalence of lumbar facet arthrosis and its relationship to age, sex, and race: an anatomic study of cadaveric specimens</article-title>. <source>Spine</source> <volume>32</volume> (<issue>19</issue>), <fpage>2058</fpage>&#x2013;<lpage>2062</lpage>. <pub-id pub-id-type="doi">10.1097/brs.0b013e318145a3a9</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fan</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Gong</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Qiu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Effects of resting modes on human lumbar spines with different levels of degenerated intervertebral discs: a finite element investigation</article-title>. <source>BMC Musculoskelet. Disord.</source> <volume>16</volume> (<issue>1</issue>), <fpage>221</fpage>. <pub-id pub-id-type="doi">10.1186/s12891-015-0686-z</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fujiwara</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Lim</surname>
<given-names>T.-H.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>H. S.</given-names>
</name>
<name>
<surname>Tanaka</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Jeon</surname>
<given-names>C.-H.</given-names>
</name>
<name>
<surname>Andersson</surname>
<given-names>G. B. J.</given-names>
</name>
<etal/>
</person-group> (<year>2000</year>). <article-title>The effect of disc degeneration and facet joint osteoarthritis on the segmental flexibility of the lumbar spine</article-title>. <source>Spine</source> <volume>25</volume> (<issue>23</issue>), <fpage>3036</fpage>&#x2013;<lpage>3044</lpage>. <pub-id pub-id-type="doi">10.1097/00007632-200012010-00011</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fujiwara</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Tamai</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yamato</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>H. S.</given-names>
</name>
<name>
<surname>Yoshida</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Saotome</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>1999</year>). <article-title>The relationship between facet joint osteoarthritis and disc degeneration of the lumbar spine: an MRI study</article-title>. <source>Eur. Spine J.</source> <volume>8</volume> (<issue>5</issue>), <fpage>396</fpage>&#x2013;<lpage>401</lpage>. <pub-id pub-id-type="doi">10.1007/s005860050193</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garay</surname>
<given-names>R. S.</given-names>
</name>
<name>
<surname>Solitro</surname>
<given-names>G. F.</given-names>
</name>
<name>
<surname>Lam</surname>
<given-names>K. C.</given-names>
</name>
<name>
<surname>Morris</surname>
<given-names>R. P.</given-names>
</name>
<name>
<surname>Albarghouthi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Lindsey</surname>
<given-names>R. W.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Characterization of regional variation of bone mineral density in the geriatric human cervical spine by quantitative computed tomography</article-title>. <source>PloS one</source> <volume>17</volume> (<issue>7</issue>), <fpage>e0271187</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0271187</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gellhorn</surname>
<given-names>A. C.</given-names>
</name>
<name>
<surname>Katz</surname>
<given-names>J. N.</given-names>
</name>
<name>
<surname>Suri</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Osteoarthritis of the spine: the facet joints</article-title>. <source>Nat. Rev. Rheumatol.</source> <volume>9</volume> (<issue>4</issue>), <fpage>216</fpage>&#x2013;<lpage>224</lpage>. <pub-id pub-id-type="doi">10.1038/nrrheum.2012.199</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Goel</surname>
<given-names>V. K.</given-names>
</name>
<name>
<surname>Monroe</surname>
<given-names>B. T.</given-names>
</name>
<name>
<surname>Gilbertson</surname>
<given-names>L. G.</given-names>
</name>
<name>
<surname>Brinckmann</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>Interlaminar shear stresses and laminae separation in a disc. Finite element analysis of the L3-L4 motion segment subjected to axial compressive loads</article-title>. <source>Spine</source> <volume>20</volume> (<issue>6</issue>), <fpage>689</fpage>&#x2013;<lpage>698</lpage>. <pub-id pub-id-type="doi">10.1097/00007632-199503150-00010</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Iorio</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Jakoi</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Singla</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Biomechanics of degenerative spinal disorders</article-title>. <source>Asian Spine J.</source> <volume>10</volume> (<issue>2</issue>), <fpage>377</fpage>&#x2013;<lpage>384</lpage>. <pub-id pub-id-type="doi">10.4184/asj.2016.10.2.377</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jarraya</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Guermazi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Lorbergs</surname>
<given-names>A. L.</given-names>
</name>
<name>
<surname>Brochin</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Kiel</surname>
<given-names>D. P.</given-names>
</name>
<name>
<surname>Bouxsein</surname>
<given-names>M. L.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>A longitudinal study of disc height narrowing and facet joint osteoarthritis at the thoracic and lumbar spine, evaluated by computed tomography: the Framingham Study</article-title>. <source>Spine J.</source> <volume>18</volume> (<issue>11</issue>), <fpage>2065</fpage>&#x2013;<lpage>2073</lpage>. <pub-id pub-id-type="doi">10.1016/j.spinee.2018.04.010</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jaumard</surname>
<given-names>N. V.</given-names>
</name>
<name>
<surname>Welch</surname>
<given-names>W. C.</given-names>
</name>
<name>
<surname>Winkelstein</surname>
<given-names>B. A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Spinal facet joint biomechanics and mechanotransduction in normal, injury and degenerative conditions</article-title>. <source>J. Biomechanical Eng.</source> <volume>133</volume> (<issue>7</issue>), <fpage>071010</fpage>. <pub-id pub-id-type="doi">10.1115/1.4004493</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kalichman</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Hunter</surname>
<given-names>D. J.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Lumbar facet joint osteoarthritis: a review</article-title>. <source>Seminars Arthritis Rheumatism</source> <volume>37</volume> (<issue>2</issue>), <fpage>69</fpage>&#x2013;<lpage>80</lpage>. <pub-id pub-id-type="doi">10.1016/j.semarthrit.2007.01.007</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kalichman</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>D. H.</given-names>
</name>
<name>
<surname>Guermazi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Berkin</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>O&#x27;Donnell</surname>
<given-names>C. J.</given-names>
</name>
<etal/>
</person-group> (<year>2008</year>). <article-title>Facet joint osteoarthritis and low back pain in the community-based population</article-title>. <source>Spine</source> <volume>33</volume> (<issue>23</issue>), <fpage>2560</fpage>&#x2013;<lpage>2565</lpage>. <pub-id pub-id-type="doi">10.1097/brs.0b013e318184ef95</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ke</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>The biomechanical influence of facet joint parameters on corresponding segment in the lumbar spine: a new visualization method</article-title>. <source>Spine</source> <volume>21</volume> (<issue>12</issue>), <fpage>2112</fpage>&#x2013;<lpage>2121</lpage>. <pub-id pub-id-type="doi">10.1016/j.spinee.2021.05.024</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>C.-J.</given-names>
</name>
<name>
<surname>Son</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Ryu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Spinal stability analysis of lumbar interbody fusion according to pelvic type and cage angle based on simplified spinal model with various pelvic indices</article-title>. <source>Front. Bioeng. Biotechnol.</source> <volume>10</volume>, <fpage>1002276</fpage>. <pub-id pub-id-type="doi">10.3389/fbioe.2022.1002276</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Sharan</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Cho</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Emam</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Hagen</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>S. Y.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>The prevalence of asymptomatic cervical and lumbar facet arthropathy: a computed tomography study</article-title>. <source>Asian Spine J.</source> <volume>13</volume> (<issue>3</issue>), <fpage>417</fpage>&#x2013;<lpage>422</lpage>. <pub-id pub-id-type="doi">10.31616/asj.2018.0235</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>H. -J.</given-names>
</name>
<name>
<surname>Chun</surname>
<given-names>H. -J.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>H. -M.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>K. -T.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>C. -K.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>B. -S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>The biomechanical influence of the facet joint orientation and the facet tropism in the lumbar spine</article-title>. <source>Spine J.</source> <volume>13</volume> (<issue>10</issue>), <fpage>1301</fpage>&#x2013;<lpage>1308</lpage>. <pub-id pub-id-type="doi">10.1016/j.spinee.2013.06.025</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kirkaldy-Willis</surname>
<given-names>W. H.</given-names>
</name>
<name>
<surname>Farfan</surname>
<given-names>H. F.</given-names>
</name>
</person-group> (<year>1982</year>). <article-title>Instability of the lumbar spine</article-title>. <source>Clin. Orthop. Relat. Res.</source> <volume>165</volume>, <fpage>110</fpage>&#x2013;<lpage>123</lpage>. <pub-id pub-id-type="doi">10.1097/00003086-198205000-00015</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Passias</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Kozanek</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wood</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>Lumbar facet joint motion in patients with degenerative disc disease at affected and adjacent levels</article-title>. <source>Spine</source> <volume>36</volume> (<issue>10</issue>), <fpage>E629</fpage>&#x2013;<lpage>E637</lpage>. <pub-id pub-id-type="doi">10.1097/brs.0b013e3181faaef7</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Gui</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Are facet joint parameters risk factors for recurrent lumbar disc herniation? a pilot study in a Chinese population</article-title>. <source>J. Clin. Neurosci.</source> <volume>77</volume>, <fpage>36</fpage>&#x2013;<lpage>40</lpage>. <pub-id pub-id-type="doi">10.1016/j.jocn.2020.05.048</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>C.-L.</given-names>
</name>
<name>
<surname>Zhong</surname>
<given-names>Z.-C.</given-names>
</name>
<name>
<surname>Hsu</surname>
<given-names>H.-W.</given-names>
</name>
<name>
<surname>Shih</surname>
<given-names>S.-L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.-T.</given-names>
</name>
<name>
<surname>Hung</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>Effect of the cord pretension of the Dynesys dynamic stabilisation system on the biomechanics of the lumbar spine: a finite element analysis</article-title>. <source>Eur. Spine J.</source> <volume>20</volume> (<issue>11</issue>), <fpage>1850</fpage>&#x2013;<lpage>1858</lpage>. <pub-id pub-id-type="doi">10.1007/s00586-011-1817-3</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>C.-W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.-L.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.-K.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>C.-M.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.-C.</given-names>
</name>
<name>
<surname>Tsai</surname>
<given-names>W.-T.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Traditional and cortical trajectory screws of static and dynamic lumbar fixation-a finite element study</article-title>. <source>BMC Musculoskelet. Disord.</source> <volume>21</volume> (<issue>1</issue>), <fpage>463</fpage>. <pub-id pub-id-type="doi">10.1186/s12891-020-03437-5</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>Y. M.</given-names>
</name>
<name>
<surname>Hutton</surname>
<given-names>W. C.</given-names>
</name>
<name>
<surname>Gharpuray</surname>
<given-names>V. M.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>Do bending, twisting, and diurnal fluid changes in the disc affect the propensity to prolapse? A viscoelastic finite element model</article-title>. <source>Spine</source> <volume>21</volume> (<issue>22</issue>), <fpage>2570</fpage>&#x2013;<lpage>2579</lpage>. <pub-id pub-id-type="doi">10.1097/00007632-199611150-00006</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mengoni</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis</article-title>. <source>Biomechanics Model. Mechanobiol.</source> <volume>20</volume> (<issue>2</issue>), <fpage>389</fpage>&#x2013;<lpage>401</lpage>. <pub-id pub-id-type="doi">10.1007/s10237-020-01403-7</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nikkhoo</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Khoz</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>C.-H.</given-names>
</name>
<name>
<surname>Niu</surname>
<given-names>C.-C.</given-names>
</name>
<name>
<surname>El-Rich</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Khalaf</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Development of a novel geometrically-parametric patient-specific finite element model to investigate the effects of the lumbar lordosis angle on fusion surgery</article-title>. <source>J. Biomechanics</source> <volume>102</volume>, <fpage>109722</fpage>. <pub-id pub-id-type="doi">10.1016/j.jbiomech.2020.109722</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Panjabi</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Krag</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Chung</surname>
<given-names>T. Q.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>Effects of disc injury on mechanical behavior of the human spine</article-title>. <source>Spine</source> <volume>9</volume> (<issue>7</issue>), <fpage>707</fpage>&#x2013;<lpage>713</lpage>. <pub-id pub-id-type="doi">10.1097/00007632-198410000-00010</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Park</surname>
<given-names>W. M.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>Y. H.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Changes in range of motion, intradiscal pressure, and facet joint force after intervertebral disc and facet joint degeneration in the cervical spine</article-title>. <source>J. Mech. Sci. Technol.</source> <volume>29</volume> (<issue>7</issue>), <fpage>3031</fpage>&#x2013;<lpage>3038</lpage>. <pub-id pub-id-type="doi">10.1007/s12206-015-0633-9</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pfirrmann</surname>
<given-names>C. W. A.</given-names>
</name>
<name>
<surname>Metzdorf</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Zanetti</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Hodler</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Boos</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Magnetic resonance classification of lumbar intervertebral disc degeneration</article-title>. <source>Spine</source> <volume>26</volume> (<issue>17</issue>), <fpage>1873</fpage>&#x2013;<lpage>1878</lpage>. <pub-id pub-id-type="doi">10.1097/00007632-200109010-00011</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Renner</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Natarajan</surname>
<given-names>R. N.</given-names>
</name>
<name>
<surname>Patwardhan</surname>
<given-names>A. G.</given-names>
</name>
<name>
<surname>Havey</surname>
<given-names>R. M.</given-names>
</name>
<name>
<surname>Voronov</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>B. Y.</given-names>
</name>
<etal/>
</person-group> (<year>2007</year>). <article-title>Novel model to analyze the effect of a large compressive follower pre-load on range of motions in a lumbar spine</article-title>. <source>J. Biomechanics</source> <volume>40</volume> (<issue>6</issue>), <fpage>1326</fpage>&#x2013;<lpage>1332</lpage>. <pub-id pub-id-type="doi">10.1016/j.jbiomech.2006.05.019</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rohlmann</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Bauer</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zander</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Bergmann</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Wilke</surname>
<given-names>H.-J.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Determination of trunk muscle forces for flexion and extension by using a validated finite element model of the lumbar spine and measured <italic>in vivo</italic> data</article-title>. <source>J. Biomechanics</source> <volume>39</volume> (<issue>6</issue>), <fpage>981</fpage>&#x2013;<lpage>989</lpage>. <pub-id pub-id-type="doi">10.1016/j.jbiomech.2005.02.019</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rohlmann</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Zander</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Rao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bergmann</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Realistic loading conditions for upper body bending</article-title>. <source>J. Biomechanics</source> <volume>42</volume> (<issue>7</issue>), <fpage>884</fpage>&#x2013;<lpage>890</lpage>. <pub-id pub-id-type="doi">10.1016/j.jbiomech.2009.01.017</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schmidt</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Heuer</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wilke</surname>
<given-names>H.-J.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Dependency of disc degeneration on shear and tensile strains between annular fiber layers for complex loads</article-title>. <source>Med. Eng. Phys.</source> <volume>31</volume> (<issue>6</issue>), <fpage>642</fpage>&#x2013;<lpage>649</lpage>. <pub-id pub-id-type="doi">10.1016/j.medengphy.2008.12.004</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sim</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Ryu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Stress distribution on spinal cord according to type of laminectomy for large focal cervical ossification of posterior longitudinal ligament based on finite element method</article-title>. <source>Bioengineering</source> <volume>9</volume> (<issue>10</issue>), <fpage>519</fpage>. <pub-id pub-id-type="doi">10.3390/bioengineering9100519</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname>
<given-names>C.-H.</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>J.-S.</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>B.-W.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>J. S.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>C.-S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Computational investigation for biomechanical characteristics of lumbar spine with various porous Ti&#x2013;6Al&#x2013;4V implant systems</article-title>. <source>Appl. Sci.</source> <volume>11</volume> (<issue>17</issue>), <fpage>8023</fpage>. <pub-id pub-id-type="doi">10.3390/app11178023</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Takigawa</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Or&#xed;as</surname>
<given-names>A. A. E.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>H. S.</given-names>
</name>
<name>
<surname>Gohgi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Udayakumar</surname>
<given-names>R. K.</given-names>
</name>
<name>
<surname>Sugisaki</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2010</year>). <article-title>Spinal kinematics and facet load transmission after total disc replacement</article-title>. <source>Spine</source> <volume>35</volume> (<issue>22</issue>), <fpage>E1160</fpage>&#x2013;<lpage>E1166</lpage>. <pub-id pub-id-type="doi">10.1097/brs.0b013e3181e5352d</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Teo</surname>
<given-names>E. C.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>K. K.</given-names>
</name>
<name>
<surname>Ng</surname>
<given-names>H. W.</given-names>
</name>
<name>
<surname>Qiu</surname>
<given-names>T. X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Determination of load transmission and contact force at facet joints of L2&#x2013;L3 motion segment using Fe method</article-title>. <source>J. Musculoskelet. Res.</source> <volume>7</volume> (<issue>02</issue>), <fpage>97</fpage>&#x2013;<lpage>109</lpage>. <pub-id pub-id-type="doi">10.1142/s0218957703001034</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Varlotta</surname>
<given-names>G. P.</given-names>
</name>
<name>
<surname>Lefkowitz</surname>
<given-names>T. R.</given-names>
</name>
<name>
<surname>Schweitzer</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Errico</surname>
<given-names>T. J.</given-names>
</name>
<name>
<surname>Spivak</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Bendo</surname>
<given-names>J. A.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>The lumbar facet joint: a review of current knowledge: part 1: anatomy, biomechanics, and grading</article-title>. <source>Skelet. Radiol.</source> <volume>40</volume> (<issue>1</issue>), <fpage>13</fpage>&#x2013;<lpage>23</lpage>. <pub-id pub-id-type="doi">10.1007/s00256-010-0983-4</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Volz</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Elmasry</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Jackson</surname>
<given-names>A. R.</given-names>
</name>
<name>
<surname>Travascio</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Computational modeling intervertebral disc pathophysiology: a review</article-title>. <source>Front. physiology</source> <volume>12</volume>, <fpage>750668</fpage>. <pub-id pub-id-type="doi">10.3389/fphys.2021.750668</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>Y.-X.</given-names>
</name>
<name>
<surname>Zhan</surname>
<given-names>H.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Effects of facet joint degeneration on stress alterations in cervical spine C5&#x2013;C6: a finite element analysis</article-title>. <source>Math. Biosci. Eng.</source> <volume>16</volume> (<issue>6</issue>), <fpage>7447</fpage>&#x2013;<lpage>7457</lpage>. <pub-id pub-id-type="doi">10.3934/mbe.2019373</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Niu</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The biomechanical influence of anterior vertebral body osteophytes on the lumbar spine: a finite element study</article-title>. <source>Spine J.</source> <volume>18</volume> (<issue>12</issue>), <fpage>2288</fpage>&#x2013;<lpage>2296</lpage>. <pub-id pub-id-type="doi">10.1016/j.spinee.2018.07.001</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Weishaupt</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zanetti</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Boos</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Hodler</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>MR imaging and CT in osteoarthritis of the lumbar facet joints</article-title>. <source>Skelet. Radiol.</source> <volume>28</volume> (<issue>4</issue>), <fpage>215</fpage>&#x2013;<lpage>219</lpage>. <pub-id pub-id-type="doi">10.1007/s002560050503</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lai</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Effect of facet-joint degeneration on the <italic>in vivo</italic> motion of the lower lumbar spine</article-title>. <source>J. Orthop. Surg. Res.</source> <volume>15</volume>, <fpage>340</fpage>&#x2013;<lpage>349</lpage>. <pub-id pub-id-type="doi">10.1186/s13018-020-01826-z</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>