
==== Front
Global Spine J
Global Spine J
spgsj
GSJ
Global Spine Journal
2192-5682
2192-5690
SAGE Publications Sage CA: Los Angeles, CA

37132375
10.1177_21925682231166612
10.1177/21925682231166612
Original Articles
Biomechanical Analysis of Lumbar Interbody Fusion Cages With Various Elastic Moduli in Osteoporotic and Non-osteoporotic Lumbar Spine: A Finite Element Analysis
Zou Da MD 123*
Yue Lihao BS 4*
Fan Zheyu BS 4*
https://orcid.org/0000-0002-3460-5425
Zhao Yi BS 4
Leng Huijie PhD 123
Sun Zhuoran MD 123
https://orcid.org/0000-0001-9512-5436
Li Weishi MD 123
1 Orthopaedics Department, 66482 Peking University Third Hospital , China
2 Ministry of Education, Engineering Research Center of Bone and Joint Precision Medicine , China
3 Orthopaedics Department, Beijing Key Laboratory of Spinal Disease Research , China
4 Orthopaedics Department, 33133 Peking University Health Science Center , China
Weishi Li, Orthopaedic Department, Peking University Third Hospital, No. 49 North Garden Road, Haidian District, Beijing 100191, China. Email: puh3liweishi@163.com
* Da Zou, Lihao Yue, and Zheyu Fan contributed equally to this article.

3 5 2023
9 2024
14 7 20532061
© The Author(s) 2023
2023
AO Spine, unless otherwise noted. Manuscript content on this site is licensed under Creative Commons Licenses
https://creativecommons.org/licenses/by-nc-nd/4.0/ This article is distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 License (https://creativecommons.org/licenses/by-nc-nd/4.0/) which permits non-commercial use, reproduction and distribution of the work as published without adaptation or alteration, without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/open-access-at-sage).

Study Design

Finite element analysis (FEA).

Objective

This study aimed to explore the effects of cage elastic modulus (Cage-E) on the endplate stress in different bone conditions: osteoporosis (OP) and non-osteoporosis (non-OP). We also explored the correlation between endplate thickness and endplate stress.

Methods

The FEA models of L4-L5 with lumbar interbody fusion were designed to access the effects of Cage-E on the endplate stress in different bone conditions. Two groups of the Young’s moduli of bony structure were assigned to simulate the conditions of OP and non-OP, and the bony endplates were analyzed in 2 kinds of thicknesses: .5 mm and 1.0 mm, with the insertion of cages with different Young’s moduli including .5, 1.5, 3, 5, 10, and 20 GPa. After model validation, an axial compressive load of 400 N and a flexion/extension moment of 7.5Nm was performed on the superior surface of L4 vertebral body in order to analyze the distribution of stress.

Results

The maximum Von Mises stress in the endplates increased by up to 100% in the OP model compared with non-OP model under the same condition of cage-E and endplate thickness. In both OP and non-OP models, the maximum endplate stress decreased as the cage-E decreased, but the maximum stress in the lumbar posterior fixation increased as the cage-E decreased. Thinner endplate thickness was associated with increased endplate stress.

Conclusion

The endplate stress is higher in osteoporotic bone than non-osteoporotic bone, which explains part of the mechanism of OP-related cage subsidence. It is reasonable to reduce the endplate stress by reducing the cage-E, but we should balance the risk of fixation failure. Endplate thickness is also important when evaluating the cage subsidence risk.

osteoporosis
cage
finite element analysis
elastic modulus
endplate
typesetterts10
==== Body
pmcIntroduction

Lumbar degenerative diseases (LDD) are considered to be a series of age-related diseases with low back pain and sciatica as the main clinical manifestations caused by the degenerative changes of the lumbar intervertebral, facet joints and the surrounding ligaments. They are the leading cause of worldwide productivity loss and the top cause of years lived with disability in 126 countries. 1 With the aging of world population, an increasing number of elderly patients with LDD suffer from poor quality of life and dysfunction, therefore need to receive lumber fusion surgery. Osteoporosis (OP) is also a common disease characterized by low bone mineral density (BMD) in the aging population, which affects nearly one-third of people aged 50 years and older. 2 It can also reduce the strength of the lumbar spine, and frequently cause the vertebral compression fractures.3,4 OP was very common in patients aged≥50 years undergoing lumbar fusion surgery, with the prevalence of up to 53.4%. 5 Accordingly, it’s essential to give attention to postoperative complications caused by OP.

Lumbar interbody fusion routinely requires the implantation of interbody cages to provide axial support for the intervertebral space. Cage subsidence is a common complication after lumbar interbody fusion, it may affect the clinical outcomes of lumbar fusion surgery, leading to reduced fusion rate and reoperation.6,7 OP is a well-known risk factor for cage subsidence.8,9 It has been reported that low BMD leads to compromised bone strength and decreased failure loads of the bony endplates.4,10 The elastic modulus of fusion cages is also believed to affect the subsidence risk. For example, the titanium cage tends to have higher rate of subsidence than polyetheretherketone (PEEK) cages because of the high Young’s modulus of titanium cages. 11 The Young’s modulus of human bone usually decreases as the increasing of age and bone loss. 12 In this case, another possible reason for the high risk of cage subsidence in osteoporotic spine is the mismatch in biomechanical characteristics between fusion cages and osteoporotic vertebral bone. The fusion cages used in clinical practice are always designed according to the elastic modulus of normal bone. We hypothesized that the fusion cages with low elastic modulus can better match the elastic modulus of osteoporotic lumbar spine, and decrease the stress on cage-endplate interface.

Finite element analysis (FEA) allows for the alteration of various parameters to systematically understand the role of each structure and condition on the spine’s biomechanical response. 13 In our study, we established 2 FE models according to different bone conditions: OP and non-osteoporosis (non-OP). We aim to access the effects of cage elastic modulus (Cage-E) on the endplate stress in different bone conditions, to guide the design of low elastic modulus cages for osteoporotic patients in our future work.

Methods

Images Data Acquisition and Reconstruction

This study was approved by Peking University Third Hospital Medical Science Research Ethics Committee (M2022083). The informed consent was waived because this was a retrospective study using patients’ computed tomography (CT) images. The lumbar CT scan of a 39-year-old man was used to develop the model geometry. The data were imported into the Mimic 23.0 software (Materialise Inc, Leuven, Belgium) in the DICOM format. According to the different gray-scale value (HU) range of the vertebral bone, the 3D geometric models of the bony part of the lumbar spine were reconstructed by using threshold selection, mask, 3D region growing and calculation.

The 3D models of L4-L5 were then imported into the Geomagic Studio 2015 (3D Systems, Ltd, USA) in the STL format. According to the original CT images, error correction and surface smoothing were performed with the STL data. The parametric surface of the entity model was further optimized to build the L4-L5 geometric model. The 3D geometric models of cage, rod and screws were created by the UG 12.0 software (Siemens PLM Software, USA). Finally, the 3D geometric models of L4-L5, cage, rods and screws were exported in the STP format.

FEA Model Establishment

The geometric model STP file was meshed in Hypermesh 14.0 (Altair, USA). The bony part, cage, screws and rods were divided into tetrahedral mesh elements with the size of 1 mm. The ligaments were divided into beam elements. The meshed model was saved as INP form and imported into Abaqus 6.13 software (SIMULIA, French) for further static analysis. Each vertebral body was divided into 3 parts: cortical bone, cancellous bone and bony endplate. Especially, the cortical bone was 1.0 mm thick while the endplate was analyzed in 2 kinds of thickness: .5 mm and 1.0 mm.14,15 A total of 7 kinds of ligaments were included: anterior longitudinal ligament (ALL), posterior longitudinal ligament (PLL), ligament flava (LF), interspinous ligament (ISL), supraspinous ligament (SSL), intertransverse ligament (ITL), and capsular ligament (CL). The bullet-shaped cage (29 mm in length, 12 mm in width, 11 mm in height) was transversely placed into the middle of intervertebral space. A total of 4 simulated pedicle screws (6 mm in diameter and 45 mm in length) were inserted through the pedicles of L4 and L5, and connected by 2 rods (5.5 mm in diameter and 50 mm in length). The cage and endplates were assigned to be bonded, ensuring there is no relative displacement between the cage-endplate interfaces, which could simulate the condition of complete fusion. Similarly, the connection between the screws and the bone was set up as binding contact. The FEA model was showed in Figure 1. The material properties for each component in the FEA model are listed in Table 1. In particular, the Young’s moduli of cortical bone and bony endplates in OP and non-OP models were determined by reference to the work of Lu et al. 12 The Young’s modulus of cancellous bone in the 2 models were determined by reference to the BMD cut-off values of quantitative computed tomography (QCT) for identifying OP and non-OP patients, and the correlation between QCT BMD and Young’s modulus. 4 Figure 1. The FEA model of L4-5 lumbar interbody fusion with pedicle screw fixation.

Table 1. Material Properties of the Components in FEA Model.

Components	Young’s modulus (MPa)	Poisson ratio	Cross-sectional area mm2	Reference	
Normal spine	Cortical bone	12 000	0.3	/	Lu et al 12	
Cancellous bone	400	0.3	/	Tawara et al 4	
Bony endplate	10 000	0.3	/	Lu et al 12	
Osteoporotic	Cortical bone	8000	0.3	/	Lu et al 12	
Cancellous bone	100	0.3	/	Tawara et al 4	
Bony endplate	6700	0.3	/	Lu et al 12	
	ALL	20	0.3	38	Polikeit et al 16	
	PLL	70	0.3	20	Polikeit et al 16	
	LF	50	0.3	60	Polikeit et al 16	
	ISL	28	0.3	35.5	Polikeit et al 16	
	SSL	28	0.3	35.5	Polikeit et al 16	
	TL	50	0.3	10	Polikeit et al 16	
	CL	20	0.3	40	Polikeit et al 16	
	Rod	110 000	.28	/	Wang et al 17	
	Screws	110 000	.28	/	Wang et al 17	
	Cage	500, 1500, 3000, 5000, 10 000, 20 000	0.3	/	/	

Validation and Load Condition

This FE model was validated by comparing the range of motion (ROM) results with those of the cadaveric study conducted by Shim et al under similar loads.18,19 In the model for validation, the intervertebral disc was complete and the material properties for annulus fibers, nucleus pulposus, and annulus fibrous substance were similar to previous research.20,21 The range-of-motion (ROM) results from this study were compared with those from the cadaveric study conducted by Shim et al after applying similar loads to current model. A good agreement between the results and the reported data was shown in Figure 2. Thus, our model can offer effective information for further analysis.Figure 2. Comparison of the ROM results between this study and the cadaveric study conducted by Shim et al. L, left; R, right.

In this study, we only analyzed the lumbar flexion and extension conditions. A reference point was created and constrained to the surface nodes of the L4. 14 To simulate the upper torso of a normal adult and flexion-extension movement, an axial compressive load of 400 N and a flexion or extension moment of 7.5Nm were assigned to the superior surface of L4 vertebral body.14,22 The inferior surface of L5 vertebral body in the model was fixed rigidly. The von Mises stresses in L4/L5 vertebral body, endplates, cage and screws were calculated.

Results

Effect of OP on the Endplate Stress

When cage-E was set to be 3 GPa (PEEK cage-E≈3 GPa), the maximum von Mises stress in the endplates was found to be higher in OP model than that in non-OP model (Table 2). This phenomenon was more evident under the flexion condition than the extension condition, and more evident in L4 endplate than L5 endplate. The maximum stress in the endplates increased by up to 100% in the OP model compared with non-OP model (Tables 2–4).Table 2. The Maximum Von Mises Stress in the Endplates When Cage-E = 3 GPa.

	Maximum von mises stress (MPa)	
Lumbar flexion	Lumbar extension	
Non-OP	OP	Non-OP	OP	
L4 endplate (0.5 mm thick)	21.84	40.75	8.763	12.74	
L5 endplate (0.5 mm thick)	21.64	37.05	7.151	8.358	
L4 endplate (1 mm thick)	16.83	31.48	6.862	10.05	
L5 endplate (1 mm thick)	18.07	26.93	5.675	5.814	

Table 3. The Maximum von Mises Stress in the Endplates in OP Model.

	Maximum von mises stress (MPa)	
	Lumbar flexion	Lumbar extension	
Cage-E (GPa)	L4 endplate (.5 mm)	L5 endplate (.5 mm)	L4 endplate (1 mm)	L5 endplate (1 mm)	L4 endplate (.5 mm)	L5 endplate (.5 mm)	L4 endplate (1 mm)	L5 endplate (1 mm)	
0.5	29.05	27.12	24.32	20.69	9.71	7.313	7.168	5.657	
1.5	38.73	33.52	28.25	24.71	11.83	7.918	8.936	5.741	
3	40.75	37.05	31.48	26.93	12.74	8.358	10.05	5.814	
5	42.22	39.83	34.72	28.39	15.01	8.713	10.91	5.891	
10	45.86	42.46	34.76	29.97	17.14	9.2	11.75	6.047	
20	47.83	44.00	35.62	31.83	18.92	9.662	12.3	6.224	

Table 4. The Maximum Von Mises Stress in the Endplates in Non-OP Model.

	Maximum von mises stress (MPa)	
	Lumbar flexion	Lumbar extension	
Cage-E (GPa)	L4 endplate (.5 mm)	L5 endplate (.5 mm)	L4 endplate (1 mm)	L5 endplate (1 mm)	L4 endplate (.5 mm)	L5 endplate (.5 mm)	L4 endplate (1 mm)	L5 endplate (1 mm)	
0.5	16.13	14.17	12.86	12.37	6.076	6.117	4.422	5.167	
1.5	20.41	19.04	15.39	15.29	6.976	6.627	5.838	5.49	
3	21.84	21.64	16.83	18.07	8.763	7.151	6.862	5.675	
5	23.72	23.36	18.60	19.84	9.263	7.593	7.658	5.827	
10	24.96	26.48	22.18	19.3	10.24	8.259	8.878	6.108	
20	27.27	28.62	23.25	20.74	13.55	8.951	9.807	6.401	

Effect of Cage-E on Endplate Stress

The results in Table 3 and Table 4 showed the maximum von Mises stress in the endplates decreased as the cage-E decreased in both OP and non-OP model. As illustrated in Figure 3, OP and thinner endplates are related to higher endplate stress under the conditions of different levels of cage-E.Figure 3. The maximum von Mises stress in the endplates of L4 and L5 with different Young’s modulus of cages. “OP/Non-OP; L4/L5; .5/1.0 mm” means “L4/L5 endplate with .5/1 mm-thick in the OP/Non-OP model”.

Effect of Cage-E on the Stress in the Lumbar Posterior Fixation

Table 5 shows the maximum von Mises stress in the lumbar posterior fixation was higher in OP model than that of non-OP model. The maximum stress in the lumbar posterior fixation decreased as the cage-E increased from .5 GPa to 20 GPa both in the non-OP and OP model. The endplate thickness slightly affected the stress (Figure 4).Table 5. The maximum von Mises stress in the lumbar posterior fixation.

	Maximum von mises stress (MPa)	
	Lumbar flexion	Lumbar extension	
Cage-E (GPa)	Non-OP (.5 mm)	Non-OP (1 mm)	OP (.5 mm)	OP (1 mm)	Non-OP (.5 mm)	Non-OP (1 mm)	OP (.5 mm)	OP (1 mm)	
0.5	51.71	51.62	66.46	64.77	40.24	39.43	50.95	48.89	
1.5	47.05	46.26	60.04	57.80	35.98	34.93	47.02	44.59	
3	45.13	44.15	58.02	55.51	34.5	34.16	45.66	43.13	
5	44.07	42.94	57.20	54.41	34.49	34.1	44.97	42.4	
10	43.05	41.67	56.45	53.32	34.44	33.94	44.32	41.67	
20	42.39	40.79	55.99	52.57	34.37	33.78	43.9	41.15	

Figure 4. The maximum von Mises stress in the lumbar posterior fixation with different cage-E, different thicknesses of endplate and different bone conditions in the FE model. “OP/Non-OP; .5/1.0 mm” means “OP/Non-OP model with .5/1 mm-thick endplate”.

Discussion

FEA has been widely used to simulate the mechanical environment of lumbar interbody fusion with the insertion of interbody cages. However, previous studies have mostly established the FEA models of lumbar spine with normal BMD, few studies have focused on the condition of interbody fusion in osteoporotic lumbar spine.12,23 Cage subsidence is 1 of the serious complications worthy of attention in lumbar spine surgery. Various studies have reported that OP is an important risk factor of cage subsidence, and associates with poor clinical outcomes.24,25 Considering the worldwide rapid aging process, it is urgent to put more emphasis on the elderly patients with high risk of OP-related cage subsidence. This study explored the influence of OP, cage-E and endplate thickness on biomechanical behaviors of the bony endplates and posterior fixation system. The results can help spine surgeons understand the mechanism of cage subsidence in osteoporotic patients, and provide evidence for designing cages with low elastic modulus which matches the mechanical properties of osteoporotic bone.

For osteoporotic patients, the Young’s modulus generally decreased by 66% for cancellous bone and by 33% for cortical bone and bony endplates.12,23 In the 1990s, Mizrahi et al first proved that the decrease of bone modulus could significantly increase peak stress in endplates. 26 Lu et al also compared the FEA model of osteoporotic vertebral bone with normal bone, and found that the cortical shell stress was sensitive to the decrease of bone modulus during aging process. 12 However, their FEA models had no surgical implants like fusion cages. The results of current study revealed that endplate maximal stress increased by up to 100% in OP model compared to non-OP model, especially for L4 endplates under the condition of lumbar flexion. The cage-E was set to be 3 GPa when we performed the comparison between OP and non-OP model, because Young’s modulus of the most common cages used in clinical practice was about 3 GPa. 27 Therefore, the increase of endplate stress in osteoporotic bone explains part of the mechanism of OP-related cage subsidence.

The high modulus of certain kinds of fusion cages, such as titanium cages, is also the risk factor of cage subsidence. 11 PEEK cages are preferred by most spine surgeons because the modulus of PEEK cages is similar to that of normal bone, but they may not match the biomechanical properties of osteoporotic bone. The mismatch between the modulus of cages and that of the bony tissue may cause stress shielding effects and bone remodeling, thus increasing subsidence rate.28,29 A few previous studies found that endplate maximal stress significantly decreased as cage-E decreased,29,30 which was in agreement with our results. However, there studies only analyzed the non-OP model. Our work indicates that the decrease of cage-E can also reduce endplate stress in osteoporotic bone. Especially, the endplate maximal stress decreased sharply when the cage-E was of ≤1.5 GPa. Of note, we found that the side effect of reducing cage-E was increasing the stress in posterior fixation system, which was also observed by Lu et al. 29 The high stress in fixation system is the risk factor of postoperative complications like pedicle screw loosening. Thus, it is important to keep the balance between reducing cage subsidence risk and maintaining the fixation strength.

The bony endplates have direct contact with fusion cages and therefore play a significant role in supporting the cages. Apart from the elastic modulus of bony endplates, the endplate thickness also decreases in aging and bone loss process. 31 Furthermore, the endplate thickness varies among different positions of the same endplate and different endplates of different vertebral levels. 32 In literature of FEA, the endplate thickness ranges from .5 mm to 1 mm.12,22,29 The effect of endplate thickness on endplate stress in interbody fusion system is still unknown. In this regard, this study found that thinner endplate was associated with higher endplate stress. Therefore, the endplate thickness in the planned position of cage also needs to be assessed when we evaluate the cage subsidence risk.

There are several limitations of our work. First, the shape and size of cages and fixation implants may also influence the stress in cage-endplate interface, we did not include such variables in the analysis. Second, there is no standard FEA model for osteoporotic lumbar spine, the FEA models in our study are unable to represent the patients with various degrees of osteoporosis. Third, considering that lumbar flexion-extension is the most common daily actions, and the endplate stress reaches a maximum in flexion position, 19 we only simulated the lumbar flexion-extension postures. The loading conditions of lateral bending and axial rotation were not analyzed. Forth, the FEA model is only a simplified system used to simulate the human body environment, more in vitro biomechanical tests and clinical research are needed to verify the effect of using cages with low elastic modulus. Lastly, many other material properties of cages are related to their performance in interbody fusion, such as the compressive strength. We should guarantee that the cages have enough compressive strength to prevent cage broken when reducing the cage-E.

In conclusion, this study reveals that higher endplate stress in osteoporotic condition is the possible mechanism of OP-related cage subsidence. Lower cage-E is associated with lower endplate stress. In this case, it is reasonable to design BMD-matched or modulus-matched cages for osteoporotic patients. One kind of BMD-matched cage has been reported in recent literature. 33 Such cages are designed to match the BMD T-score measured by dual energy x-ray absorptiometry (DXA), but it has been reported that DXA results are less accurate than BMD measured by CT Hounsfield units (HU) in older patients requiring lumbar fusion.5,34 Instead, we propose to design customized cages matching the biomechanical properties of patients with various degrees of osteoporosis measured in CT HU (mild OP:≤110HU; severe OP: 50-80 HU; extremely severe OP:≤50HU).35,36 In future, the biomechanical-matched cages may benefit more older patients. Besides, thinner endplate can also lead to higher risk of cage subsidence, we should pay attention to endplate thickness assessment when we design customized cages and choose where to place the cages.

Supplemental Material

Supplemental Material - Biomechanical Analysis of Lumbar Interbody Fusion Cages With Various Elastic Moduli in Osteoporotic and Non-osteoporotic Lumbar Spine: A Finite Element Analysis

Supplemental Material for Biomechanical Analysis of Lumbar Interbody Fusion Cages With Various Elastic Moduli in Osteoporotic and Non-osteoporotic Lumbar Spine: A Finite Element Analysis by Da Zou, MD, Lihao Yue, BS, Zheyu Fan, BS, Yi Zhao, BS, Huijie Leng, PhD, Zhuoran Sun, MD, and Weishi Li, MD in Global Spine Journal

Supplemental Material - Biomechanical Analysis of Lumbar Interbody Fusion Cages With Various Elastic Moduli in Osteoporotic and Non-osteoporotic Lumbar Spine: A Finite Element Analysis

Supplemental Material for Biomechanical Analysis of Lumbar Interbody Fusion Cages With Various Elastic Moduli in Osteoporotic and Non-osteoporotic Lumbar Spine: A Finite Element Analysis by Da Zou, MD, Lihao Yue, BS, Zheyu Fan, BS, Yi Zhao, BS, Huijie Leng, PhD, Zhuoran Sun, MD, and Weishi Li, MD in Global Spine Journal

Supplemental Material - Biomechanical Analysis of Lumbar Interbody Fusion Cages With Various Elastic Moduli in Osteoporotic and Non-osteoporotic Lumbar Spine: A Finite Element Analysis

Supplemental Material for Biomechanical Analysis of Lumbar Interbody Fusion Cages With Various Elastic Moduli in Osteoporotic and Non-osteoporotic Lumbar Spine: A Finite Element Analysis by Da Zou, MD, Lihao Yue, BS, Zheyu Fan, BS, Yi Zhao, BS, Huijie Leng, PhD, Zhuoran Sun, MD, and Weishi Li, MD in Global Spine Journal

Acknowledgment

We thank Xinyu Li and Shengqiang Liu at Beijing AKEC Medical Company Ltd For their valuable input on the FEA model for this study.

ORCID iDs

Yi Zhao https://orcid.org/0000-0002-3460-5425

Weishi Li https://orcid.org/0000-0001-9512-5436

The authors declared no potential conflicts of interest with respect to the research, authorship, and publication of this article.

Funding: The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was funded by the National Natural Science Foundation of China (Grant No.82202779), and the Research and Application of Clinical Diagnosis and Treatment Technology Foundation of Beijing (Grant No. Z201100005520073).

Supplemental Material: Supplemental material for this article is available online.
==== Refs
References

1 GBD 2017 Disease and Injury Incidence and Prevalence Collaborators . Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392 (10159 ):1789-1858. doi:10.1016/s0140-6736(18)32279-7 30496104
2 Zeng Q Li N Wang Q , et al. The prevalence of osteoporosis in China, a nationwide, multicenter DXA survey. J Bone Miner Res. 2019;34 (10 ):1789-1797. doi:10.1002/jbmr.3757 31067339
3 Kim G Cole JH Boskey AL Baker SP van der Meulen MCH . Reduced tissue-level stiffness and mineralization in osteoporotic cancellous bone. Calcif Tissue Int. 2014;95 (2 ):125-131. doi:10.1007/s00223-014-9873-4 24888692
4 Tawara D Sakamoto J Murakami H Kawahara N Oda J Tomita K . Mechanical evaluation by patient-specific finite element analyses demonstrates therapeutic effects for osteoporotic vertebrae. J Mech Behav Biomed Mater. 2010;3 (1 ):31-40. DOI: 10.1016/j.jmbbm.2009.03.001 19878900
5 Zou D Jiang S Zhou S , et al. Prevalence of osteoporosis in patients undergoing lumbar fusion for lumbar degenerative diseases: A combination of DXA and hounsfield units. Spine(Phila Pa 1976). 2020;45 (7 ):E406-e410. doi:10.1097/brs.0000000000003284 31725127
6 Le TV Baaj AA Dakwar E , et al. Subsidence of polyetheretherketone intervertebral cages in minimally invasive lateral retroperitoneal transpsoas lumbar interbody fusion. Spine(Phila Pa 1976). 2012;37 (14 ):1268-1273.22695245
7 Dua K Kepler CK Huang RC Marchenko A . Vertebral body fracture after anterolateral instrumentation and interbody fusion in two osteoporotic patients. Spine J. 2010;10 (9 ):e11-e15. doi:10.1016/j.spinee.2010.07.007
8 Cho JH Hwang CJ Kim H Joo Y-S Lee D-H Lee CS . Effect of osteoporosis on the clinical and radiological outcomes following one-level posterior lumbar interbody fusion. J Orthop Sci. 2018;23 (6 ):870-877. doi:10.1016/j.jos.2018.06.009 30431006
9 Formby PM Kang DG Helgeson MD Wagner SC . Clinical and radiographic outcomes of transforaminal lumbar interbody fusion in patients with osteoporosis. Global Spine J. 2016;6 (7 ):660-664. doi:10.1055/s-0036-1578804 27781185
10 Hou Y Yuan W . Influences of disc degeneration and bone mineral density on the structural properties of lumbar end plates. Spine J. 2012;12 (3 ):249-256. DOI: 10.1016/j.spinee.2012.01.021 22366078
11 Seaman S Kerezoudis P Bydon M Torner JC Hitchon PW . Titanium vs. polyetheretherketone (PEEK) interbody fusion: Meta-analysis and review of the literature. J Clin Neurosci. 2017;44 :23-29. doi:10.1016/j.jocn.2017.06.062 28736113
12 Lu Y Rosenau E Paetzold H , et al. Strain changes on the cortical shell of vertebral bodies due to spine ageing: A parametric study using a finite element model evaluated by strain measurements. Proc Inst Mech Eng H. 2013;227 (12 ):1265-1274. doi:10.1177/0954411913501293 23990044
13 Ellingson AM Shaw MN Giambini H An KN . Comparative role of disc degeneration and ligament failure on functional mechanics of the lumbar spine. Comput Methods Biomech Biomed Eng. 2016;19 (9 ):1009-1018. doi:10.1080/10255842.2015.1088524
14 Sin DA Heo DH . Comparative finite element analysis of lumbar cortical screws and pedicle screws in transforaminal and posterior lumbar interbody fusion. Neurospine. 2019;16 (2 ):298-304. doi:10.14245/ns.1836030.015 31154694
15 Wang Y Battié MC Boyd SK Videman T . The osseous endplates in lumbar vertebrae: Thickness, bone mineral density and their associations with age and disk degeneration. Bone. 2011;48 (4 ):804-809. DOI: 10.1016/j.bone.2010.12.005 21168539
16 Polikeit A Ferguson SJ Nolte LP Orr TE . Factors influencing stresses in the lumbar spine after the insertion of intervertebral cages: finite element analysis. Eur Spine J . 2003;12 (4 ):413-420. 10.1007/s00586-002-0505-8 12955610
17 Wang W Pei B Pei Y , et al. Biomechanical effects of posterior pedicle fixation techniques on the adjacent segment for the treatment of thoracolumbar burst fractures: A biomechanical analysis. Comput Methods Biomech Biomed Eng. 2019;22 (13 ):1083-1092. doi:10.1080/10255842.2019.1631286
18 Shim CS Park SW Lee SH Lim TJ Chun K Kim DH . Biomechanical evaluation of an interspinous stabilizing device, Locker. Spine(Phila Pa 1976). 2008;33 (22 ):E820-E827. doi:10.1097/BRS.0b013e3181894fb1.18923305
19 Xu H Tang H Guan X , et al. Biomechanical comparison of posterior lumbar interbody fusion and transforaminal lumbar interbody fusion by finite element analysis. Neurosurgery. 2013;72 (1 suppl Operative ):21-26. doi:10.1227/NEU.0b013e3182742a69 23037820
20 Chen SH Tai CL Lin CY Hsieh PH Chen WP . Biomechanical comparison of a new stand-alone anterior lumbar interbody fusion cage with established fixation techniques - a three-dimensional finite element analysis. BMC Muscoskel Disord. 2008;9 :88.
21 Wu TK Meng Y Wang BY , et al. Biomechanics following skip-level cervical disc arthroplasty versus skip-level cervical discectomy and fusion: a finite element-based study. BMC Muscoskel Disord. 2019;20 (1 ):49. doi:10.1186/s12891-019-2425-3
22 Zhu H Zhong W Zhang P , et al. Biomechanical evaluation of autologous bone-cage in posterior lumbar interbody fusion: A finite element analysis. BMC Muscoskel Disord. 2020;21 (1 ):379. doi:10.1186/s12891-020-03411-1
23 Liu ZX Gao ZW Chen C , et al. Effects of osteoporosis on the biomechanics of various supplemental fixations co-applied with oblique lumbar interbody fusion (OLIF): a finite element analysis. BMC Muscoskel Disord. 2022;23 (1 ):794. doi:10.1186/s12891-022-05645-7
24 Amorim-Barbosa T Pereira C Catelas D , et al. Risk factors for cage subsidence and clinical outcomes after transforaminal and posterior lumbar interbody fusion. Eur J Orthop Surg Traumatol. 2022;32 (7 ):1291-1299. doi:10.1007/s00590-021-03103-z 34462820
25 Yao Y-C Chou P-H Lin H-H Wang S-T Liu C-L Chang M-C . Risk factors of cage subsidence in patients received minimally invasive transforaminal lumbar interbody fusion. Spine(Phila Pa 1976) 2020;45 (19 ):E1279-E1285. doi:10.1097/brs.0000000000003557 32472823
26 Mizrahi J Silva MJ Keaveny TM Edwards WT Hayes WC . Finite-element stress analysis of the normal and osteoporotic lumbar vertebral body. Spine(Phila Pa 1976). 1993;18 (14 ):2088-2096. doi:10.1097/00007632-199310001-00028 8272965
27 Liao C Li Y Tjong SC . Polyetheretherketone and its composites for bone replacement and regeneration. Polymers (Basel). 2020;12 (12 ):2858. doi:10.3390/polym12122858 33260490
28 Wang Ln Yuan B Chen F , et al. Ability of a novel biomimetic titanium alloy cage in avoiding subsidence and promoting fusion: A goat spine model study. Mater Des. 2022;213 :110361. DOI: 10.1016/j.matdes.2021.110361
29 Lu T Ren J Sun Z , et al. Relationship between the elastic modulus of the cage material and the biomechanical properties of transforaminal lumbar interbody fusion: A logarithmic regression analysis based on parametric finite element simulations. Comput Methods Progr Biomed. 2022;214 :106570. doi:10.1016/j.cmpb.2021.106570
30 Zhang Z Li H Fogel GR Xiang D Liao Z Liu W . Finite element model predicts the biomechanical performance of transforaminal lumbar interbody fusion with various porous additive manufactured cages. Comput Biol Med. 2018;95 :167-174. doi:10.1016/j.compbiomed.2018.02.016 29501735
31 Ritzel H Amling M Pösl M Hahn M Delling G . The thickness of human vertebral cortical bone and its changes in aging and osteoporosis: A histomorphometric analysis of the complete spinal column from thirty-seven autopsy specimens. J Bone Miner Res. 1997;12 (1 ):89-95. doi:10.1359/jbmr.1997.12.1.89 9240730
32 Palepu V Rayaprolu SD Nagaraja S . Differences in trabecular bone, cortical shell, and endplate microstructure across the lumbar spine. Internet J Spine Surg. 2019;13 (4 ):361-370.
33 Falowski SM Koga SF Northcutt T Garamszegi L Leasure J Block JE . Improving the management of patients with osteoporosis undergoing spinal fusion: Tthe need for a bone mineral density-matched interbody cage. Orthop Res Rev. 2021;13 :281-288. doi:10.2147/ORR.S339222 34934366
34 Zou D Li W Deng C Du G Xu N . The use of CT Hounsfield unit values to identify the undiagnosed spinal osteoporosis in patients with lumbar degenerative diseases. Eur Spine J. 2019;28 (8 ):1758-1766. doi:10.1007/s00586-018-5776-9 30306332
35 Zou D Ye K Tian Y , et al. Characteristics of vertebral CT Hounsfield units in elderly patients with acute vertebral fragility fractures. Eur Spine J. 2020;29 (5 ):1092-1097. doi:10.1007/s00586-020-06363-1 32157387
36 Ye K Zou D Zhou F Li W Tian Y . Low vertebral CT Hounsfield units: A risk factor for new osteoporotic vertebral fractures after the treatment of percutaneous kyphoplasty. Arch Osteoporosis. 2022;17 (1 ):137. doi:10.1007/s11657-022-01177-8
