阮朝陽 何永清 曾強(qiáng)華 朱群威 付強(qiáng)

腰椎間孔成形幅度對腰椎生物力學(xué)影響的三維有限元分析

阮朝陽 何永清 曾強(qiáng)華 朱群威 付強(qiáng)

目的 應(yīng)用腰椎L3～5節(jié)段三維有限元模型，分析椎間盤摘除+不同幅度椎間孔成形前后對腰椎生物力學(xué)的影響。方法 選取1位20歲既往無腰椎疾病史中國男性志愿者的CT數(shù)據(jù)，在Mimics 15.0軟件中建立L3～5節(jié)段三維有限元模型，在Geomagic 12.0中修補(bǔ)、降噪及曲面化，在Pro/E 5.0中行椎間盤的曲面建模，修改髓核和纖維環(huán)的材料參數(shù)模擬L4～5椎間盤中度退變（模型M1），在Hypermesh 12.0中進(jìn)行有限元網(wǎng)格處理，并導(dǎo)入ABAQUS軟件進(jìn)行分析。模擬經(jīng)皮椎間孔脊柱內(nèi)鏡技術(shù)側(cè)后入路，以上關(guān)節(jié)突尖部到下位椎體后上緣中點為穿刺基線建立通道，去除椎間盤左后側(cè)約1/4的纖維環(huán)中部及1/4的髓核，以模擬腰椎間盤摘除手術(shù)，構(gòu)建模型M2，以圓柱體代替環(huán)鋸模擬切除上關(guān)節(jié)突部分骨質(zhì)行椎間孔擴(kuò)大成形術(shù)，構(gòu)建一級椎間孔成形模型M3（環(huán)鋸直徑5mm）、二級椎間孔成形模型M4（環(huán)鋸直徑6.5mm）和三級椎間孔成形模型M5（環(huán)鋸直徑7.5mm），給予特定加載條件，比較其在前屈、后伸、左右側(cè)彎、左右旋轉(zhuǎn)6種工況下的生物力學(xué)特征。 結(jié)果 建立了有效的L3～5三維有限元模型，除了在后伸工況下L4椎體的位移無明顯變化外，其他5種工況下椎間孔成形后L4椎體的位移均有所增加，但成形的幅度對L4椎體的位移無明顯影響。M1模型在各種工況下L4～5關(guān)節(jié)突關(guān)節(jié)的應(yīng)力左側(cè)較右側(cè)大。M2模型在腰椎前屈及右旋工況下L4～5左側(cè)關(guān)節(jié)突關(guān)節(jié)的應(yīng)力下降，右側(cè)的應(yīng)力上升；腰椎后伸下L4～5左側(cè)關(guān)節(jié)突關(guān)節(jié)的應(yīng)力上升，右側(cè)的應(yīng)力下降；腰椎左側(cè)屈工況下L4～5關(guān)節(jié)突關(guān)節(jié)的左右兩側(cè)應(yīng)力均上升；腰椎右側(cè)屈工況下L4～5關(guān)節(jié)突關(guān)節(jié)的左右兩側(cè)應(yīng)力均下降。除腰椎左旋工況下M5模型的L4～5右側(cè)關(guān)節(jié)突關(guān)節(jié)應(yīng)力增加明顯外，其他各種工況下M3、M4、M5模型L4～5關(guān)節(jié)突關(guān)節(jié)應(yīng)力變化不明顯。 結(jié)論 在精確穿刺的指引下，1～3級椎間孔成形對術(shù)后即刻的腰椎穩(wěn)定及關(guān)節(jié)突關(guān)節(jié)的應(yīng)力分布無明顯影響。

腰椎 生物力學(xué) 椎間盤移位 應(yīng)力 物理 有限元分析

近年來，腰椎間盤突出癥（lumbar disc herniation，LDH）的治療越來越趨于精細(xì)和微創(chuàng)。經(jīng)皮椎間孔脊柱內(nèi)鏡（transforaminal endoscopic spine system，TESSYS）技術(shù)與傳統(tǒng)手術(shù)療效相當(dāng)，但因其具有創(chuàng)傷小、出血少、快速康復(fù)等優(yōu)點，受到越來越多的關(guān)注。工作套管從側(cè)后方進(jìn)入椎管內(nèi)硬膜囊前方間隙的過程中，關(guān)節(jié)突是主要障礙，因此，對椎間孔擴(kuò)大成形，有時是必需的舉措。已有大量系統(tǒng)性的生物力學(xué)研究結(jié)果證實傳統(tǒng)后路逐級切除關(guān)節(jié)突關(guān)節(jié)對腰椎穩(wěn)定性有一定影響[1-3]，但TESSYS技術(shù)行腰椎間孔擴(kuò)大成形對腰椎穩(wěn)定性有何影響，如何避免腰椎間孔擴(kuò)大成形后對椎骨造成的不適當(dāng)?shù)膽?yīng)力、應(yīng)變？尸體標(biāo)本研究[4]并未給出滿意的答案。本研究利用三維有限元模型進(jìn)行模擬實驗，分析腰椎椎間孔成形幅度對腰椎生物力學(xué)的影響。

1 材料和方法

1.1 三維有限元模型的建立 選取1位20歲既往無腰椎疾病史的中國籍男性志愿者，采用64排螺旋CT機(jī)（德國西門子公司）對其L1～S1節(jié)段連續(xù)掃描，獲得層厚為0.625mm的連續(xù)斷面圖像，以DICOM格式保存。將DICOM數(shù)據(jù)導(dǎo)入Mimics 15.0軟件（比利時Materialise公司），建立L3～5節(jié)段三維模型，導(dǎo)出STL文件，在Geomagic 12.0（美國Geomagic公司）中修補(bǔ)、降噪及曲面化，導(dǎo)出STP格式，依據(jù)椎間盤解剖學(xué)形態(tài)，在Pro/E 5.0（美國參數(shù)技術(shù)公司）中行椎間盤的曲面建模，最后整體導(dǎo)出IGES格式，在Hypermesh 12.0（美國Altair公司）中進(jìn)行有限元網(wǎng)格處理，并導(dǎo)入ABAQUS（美國ABAQUS軟件公司）軟件進(jìn)行分析。有限元模型各結(jié)構(gòu)的材料參數(shù)參考相關(guān)文獻(xiàn)[3-9]。骨皮質(zhì)、骨松質(zhì)及終板單獨建模，修改髓核和纖維環(huán)的材料參數(shù)模擬L4～5椎間盤中度退變[10]。韌帶（包括棘上韌帶、棘間韌帶、橫突間韌帶、黃韌帶、后縱韌帶和前縱韌帶等）采用Truss單元。具體模型見圖1。

1.2 椎間孔成形模型的構(gòu)建 在L4～5椎間盤中度退變模型（model1，M1）的基礎(chǔ)上，模擬TESSYS技術(shù)側(cè)后入路，以上關(guān)節(jié)突尖部到下位椎體后上緣中點為穿刺基線建立通道，去除椎間盤左后側(cè)約1/4的纖維環(huán)中部及1/4的髓核，以模擬腰椎間盤摘除手術(shù)，構(gòu)建模型M2，以圓柱體代替環(huán)鋸模擬切除上關(guān)節(jié)突部分骨質(zhì)行椎間孔擴(kuò)大成形術(shù)，構(gòu)建一級椎間孔成形模型模型M3（環(huán)鋸直徑5mm）、二級椎間孔成形模型M4（環(huán)鋸直徑6.5mm）和三級椎間孔成形模型M5（環(huán)鋸直徑7.5mm）。

圖1 正常腰椎椎體L3～5有限元模型圖

1.3 設(shè)定模型邊界和載荷條件 約束模型L5的下表面，在L3椎體上表面向終板軸施加負(fù)荷為500N的正壓力，在矢狀面、冠狀面和水平面上分別施加10N·m的純扭矩。模擬L4～5節(jié)段在左右旋轉(zhuǎn)、左右側(cè)屈、前屈和后伸等6種工況下的運動狀態(tài)。具體模型見圖2（插頁）。

圖2 L3～5運動節(jié)段模型在6種工況下的綜合應(yīng)力云圖（a：左旋；b：右旋；c：左側(cè)屈；d：右側(cè)屈；e：前屈；f：后伸）

1.4 觀察項目 給予上述加載條件，取L4椎體上表面前、后、左右側(cè)緣的中點及上表面的中心點5個節(jié)點在不同幅度椎間孔成形后的位移改變，取關(guān)節(jié)面凸面中心的節(jié)點分析L4～5關(guān)節(jié)突關(guān)節(jié)面的應(yīng)力分布變化。

2 結(jié)果

2.1 有限元模型驗證 將正常腰椎三維有限元模型經(jīng)過上述載荷和邊界條件后，進(jìn)行力學(xué)計算，將各節(jié)段的活動度與Shim等[11]的標(biāo)本生物力學(xué)測試結(jié)果作比較，結(jié)果表明本研究的正常腰椎L3～5三維有限元模型在載荷工況下，其活動度均處在合理結(jié)果范圍內(nèi)，證明了本模型的有效性，見表1。

表1 正常腰椎模型不同工況下的活動度與文獻(xiàn)報告數(shù)據(jù)比較（°）

2.2 5組模型在6種工況下的位移變化 L4椎體上表面前、后、左右側(cè)緣及L4椎體上表面中心5個固定點在6種工況下的位移見圖3。

由圖3可見，除了在后伸工況下L4椎體的位移無明顯變化外，其他5種工況下椎間孔成形后L4椎體的位移均有所增加，但成形的幅度對L4椎體的位移無明顯影響。

2.3 5組模型L4～5關(guān)節(jié)突關(guān)節(jié)的應(yīng)力分布比較 M1模型在各種工況下L4～5關(guān)節(jié)突關(guān)節(jié)的應(yīng)力左側(cè)均較右側(cè)大。M2模型在腰椎前屈及右旋工況下L4～5左側(cè)關(guān)節(jié)突關(guān)節(jié)的應(yīng)力下降，右側(cè)的應(yīng)力上升；腰椎后伸工況下L4～5左側(cè)關(guān)節(jié)突關(guān)節(jié)的應(yīng)力上升，右側(cè)的應(yīng)力下降；腰椎左側(cè)屈工況下L4～5關(guān)節(jié)突關(guān)節(jié)的左右兩側(cè)應(yīng)力均上升；腰椎右側(cè)屈工況下L4～5關(guān)節(jié)突關(guān)節(jié)的左右兩側(cè)應(yīng)力均下降。除腰椎左旋工況下M5模型的L4～5右側(cè)關(guān)節(jié)突關(guān)節(jié)應(yīng)力增加明顯外，其他各種工況下椎間孔1、2、3的成形后（M3、M4、M5）L4～5關(guān)節(jié)突關(guān)節(jié)應(yīng)力折線幾乎平直，說明環(huán)鋸直徑5、6.5、7mm的椎間孔成形對L4～5關(guān)節(jié)突關(guān)節(jié)應(yīng)力變化幾乎無明顯影響，見圖4。

圖3 L4椎體上表面前、后、左、右側(cè)緣及L4椎體上表面中心5個固定點在6種工況下的位移（a：L4椎體上表面后緣中點；b：L4椎體上表面前緣中點；c：L4椎體上表面左側(cè)緣中點；d：L4椎體上表面右側(cè)緣中點；e：L4椎體上表面中心點）

圖4 5組模型L4～5關(guān)節(jié)突關(guān)節(jié)的應(yīng)力分布比較（a：L4～5左側(cè)關(guān)節(jié)突關(guān)節(jié)的應(yīng)力分布；b：L4～5右側(cè)關(guān)節(jié)突關(guān)節(jié)的應(yīng)力分布）

3 討論

本研究建立了有效的L3～5三維有限元模型，通過改變椎間盤材料性質(zhì)建立了L4～5椎間盤中度退變模型，并在此基礎(chǔ)上模擬手術(shù)建立了椎間盤摘除、椎間孔1～3級成形+椎間盤摘除的模型，觀察了6種工況下腰椎節(jié)段位移及關(guān)節(jié)突關(guān)節(jié)的應(yīng)力分布。

3.1 髓核摘除術(shù)及髓核摘除術(shù)+不同幅度椎間孔成形后對腰椎穩(wěn)定性的影響 本研究發(fā)現(xiàn)，除了在后伸工況下L4椎體的位移無明顯變化外，其他5種工況下椎間盤摘除+1～3級椎間孔成形后L4椎體的位移均有所增加，但1～3級椎間孔成形幅度對L4椎體的位移影響甚小。換而言之，不同幅度椎間孔成形對腰椎位移無明顯影響，而椎間盤的摘除對腰椎位移的改變相對要大。說明髓核的正常形態(tài)、結(jié)構(gòu)在維持脊柱功能單元的穩(wěn)定性有重要意義，宋升等[12]、李海波等[13]的研究也證實了這一點。

3.2 髓核摘除術(shù)及髓核摘除術(shù)+不同幅度椎間孔成形后對腰椎關(guān)節(jié)突關(guān)節(jié)應(yīng)力分布的影響 在各種工況、各模型下L4～5關(guān)節(jié)突關(guān)節(jié)的應(yīng)力左側(cè)較右側(cè)大，這可能與模型的腰椎間盤突出及椎間盤摘除的模擬手術(shù)在左側(cè)有關(guān)，而模型在不同工況下L4～5左、右側(cè)關(guān)節(jié)突關(guān)節(jié)應(yīng)力變化規(guī)律可能與腰椎的關(guān)節(jié)突關(guān)節(jié)呈冠狀面排列方式有關(guān)，M5模型在腰椎左旋工況下L4～5右側(cè)關(guān)節(jié)突關(guān)節(jié)應(yīng)力增加明顯，可能是椎間孔成形幅度增大破壞左側(cè)的關(guān)節(jié)突關(guān)節(jié)的關(guān)節(jié)囊，對左側(cè)關(guān)節(jié)突關(guān)節(jié)的約束力下降導(dǎo)致右側(cè)關(guān)節(jié)突關(guān)節(jié)的應(yīng)力反常增大，而其他工況下M2、M3、M4、M5各模型間變化不大，這提示穿刺定位要求精確無誤即以上關(guān)節(jié)突尖部到下位椎體后上緣中點為穿刺基線建立通道，過多地去除上關(guān)節(jié)突可能影響腰椎的穩(wěn)定性。李振宙等[4]的結(jié)果表明大幅度椎間孔成形時上關(guān)節(jié)突尖部及部分關(guān)節(jié)面、腹側(cè)關(guān)節(jié)囊被切除，腰椎側(cè)屈穩(wěn)定性顯著下降。

髓核摘除后瘢痕組織修復(fù)后對腰椎節(jié)段剛度和后部結(jié)構(gòu)應(yīng)力的影響較髓核摘除后即刻明顯增大[13]，它對椎間盤的功能和應(yīng)力分布有不利影響[14]。本研究為一次性載荷分析，沒有充分考慮疲勞載荷的影響，髓核摘除術(shù)+椎間孔成形術(shù)后，相應(yīng)椎體的位移及關(guān)節(jié)突關(guān)節(jié)的應(yīng)力較術(shù)前有所改變，雖然統(tǒng)計學(xué)上無明顯差異，但腰椎反復(fù)活動的累積效應(yīng)可能引發(fā)相應(yīng)節(jié)段的退變。

[1] 趙凡,劉正,王炳強(qiáng),等.有限元模擬單節(jié)段腰椎小關(guān)節(jié)分級切除對腰椎穩(wěn)定性的影響[J].中華醫(yī)學(xué)雜志,2015,95(13):973-977.

[2] 靳安民,袁野,曹虹,等.單雙側(cè)小關(guān)節(jié)分級切除的有限元分析[J].中國臨床解剖學(xué)雜志,2002,20(6):473-474.

[3] 趙勇,李玉茂,李平生,等.單側(cè)小關(guān)節(jié)分級切除對腰椎穩(wěn)定性影響的三維有限元分析[J].實用骨科雜志,2009,15(10):764-767.

[4] 李振宙,侯樹勛,吳聞文,等.經(jīng)皮側(cè)后路腰椎間孔成形術(shù)對腰椎解剖及生物力學(xué)影響的實驗研究[J].中國骨與關(guān)節(jié)雜志,2010,9(6): 503-508.

[5] 蘇晉,趙文志,陳秉智,等.建立全腰椎有限元接觸模型[J].醫(yī)用生物力學(xué),2010,25(3):200-205.

[6] Pitzen T,Geisler F,Matthis D,et al.A finite element model for predicting the biomechanical behaviour of the human lumbar spine[J].Control Engineering Practice,2002,10(1):83-90.doi: org/10.1016/s0967-0661(01)00129-0.

[7] Bresnahan L,Ogden A T,Natarajan R N,et al.A biomechanical evaluation of graded posterior element removal for treatment of lumbar stenosis:comparison of a minimally invasive approach with two standard laminectomy techniques[J].Spine,2009,34 (34):17-23.doi:10.1097/BRS.0b013e318191438b.

[8] Serhan H,Mhatre D,Defossez H,et al.Motion-preserving technologies for degenerative lumbar spine:the past,present,and future horizons[J].Sas Journal,2011,5(3):75-89.doi:10.1016/j. esas.2011.05.001.eCollection 2011.

[9] 黃菊英,李海云,吳浩.腰椎間盤突出癥力學(xué)特征的仿真計算方法[J].醫(yī)用生物力學(xué),2012,27(1):96-101.

[10] Ruberté L M,Natarajan R N,Andersson G B.Influence of single-level lumbar degenerative disc disease on the behavior of the adjacent segments-a finite element model study[J].Journal of Biomechanics,2009,42(3):341-348.doi:10.1016/j.jbiomech. 2008.11.024.

[11] Shim C S,Park S W,Lee S H,et al.Biomechanical evaluation of an interspinous stabilizing device,locker[J].Spine(Phila Pa 1976),2008,33(22):E820-827.doi:10.1097/BRS.0b013e31818-94fb1.

[12] 宋升,芮永軍,孫振,等.模擬腰椎間盤髓核摘除后疲勞載荷對腰椎穩(wěn)定性的影響[J].中國組織工程研究與臨床康復(fù),2011,15(48): 9040-9043.

[13] 李海波,方杰,陳其昕,等.髓核摘除術(shù)對腰椎生物力學(xué)特性影響的有限元研究[J].中國臨床解剖學(xué)雜志,2011,29(1):89-93.

[14] Hayashi S,Taira A,Inoue G,et al.TNF-alpha in nucleus pulposus induces sensory nerve growth:a study of the mechanism of discogenic low back pain using TNF-alpha-deficient mice [J].Spine,2008,33(14):1542-1546.doi:10.1097/BRS.0b013e3-18178e5ea.

Lumbar biomechanical characteristics after foraminoplasty in different amplitudes:a finite element analysis

RUAN Zhaoyang,HE Yongqing,ZENG Qianghua,et al.Department of Orthopedics,Haining People's Hospital,Haining 314400,China

Objective To analyze lumbar biomechanical characteristics after foraminoplasty in different amplitudes based on a three-dimensional finite element model of L3～5segments of the lumbar spine. Methods A 20-year-old Chinese male volunteer who had no history of lumbar disease was selected.Using the obtained CT data,the three-dimensional model of L3～5segment was established by Mimics15.0 software,and curved,repaired and denoised by Geomagic 12.0 software,the degeneration of L4～5(model M1)was simulated by modifying the material parameters of the nucleus pulposus and annulus fibrosus of intervertebral disc and grided by Hypermesh 12,and then was analyzed by ABAQUS software.The lumbar discectomy was simulated via percutaneous transforaminal endoscopic posterior approach to establish a channel between the midpoint of the tip of the protruding points and the top of the vertebral body,which removed the superior facet bone for foraminal plasty,the whole layer of middle of the fiber ring at the 1/4 right of the inter-vertebral and about a quarter of the nucleus pulposus to establish model M2.The model of the first-stage foraminoplasty(model M3,diameter 5mm),the second-stage foraminoplasty(model M4, diameter 6.5mm)and three-stage foraminoplasty(model M5,diameter 7.5mm)were established by using the cylinder instead of the trepan remove the superior facet bone for foraminoplasty.Under specific loading conditions,biomechanical characteristics of 6 kinds of loads (the forward bends,stretch,left and right side bend,rotate)were compared. Results An effective L3～5three-dimensional finite element model was established.In addition to the displacement of the L4vertebral body during extension,the forming of the lumbar intervertebral hole had no significant difference on the displacement of intervertebral lumbar.Different forming amplitude of intervertebral lumbar had no significant difference on the stress of the protruding joints.In the M1 model,the stress of the L4～5left facet joints was larger than the right side under various conditions.In the M2 model,the stress L4～5left facet joints decreased in lumbar flexion and dextral,the stress increased on the right side.The stress on the right side of L4～5facet joints increased in left lumbar flexion,but decreased on the right lumbar flexion.In addition to the stress of the L4～5right facet joint increased significantly,and the stress changes of L4～5facet joints of M3,M4,M5 models were not significant under other working conditions. Conclusion Under the guidance of the precise puncture,the hole forming of grade 1 to grade 3 lumbar intervertebral hole has no significant effect on the stabilization of the lumbar and the stress distribution of the joints.

Lumbar vertebrae Biomechanics Intervertebral disk displacement Stress MechanicalFinite element analysis analysis

2017-01-22）

（本文編輯：陳麗）

10.12056/j.issn.1006-2785.2017.39.10.2017-164

浙江省醫(yī)藥衛(wèi)生科技計劃項目（2015KYB395）

314400 海寧市人民醫(yī)院骨二科（阮朝陽、何永清、曾強(qiáng)華、朱群威）；上海長海醫(yī)院脊柱外科（付強(qiáng)）

阮朝陽，E-mail：zhyangr3061@126.com