劉 巖 張恒巖 王 君 張 嘉, 郇 勇,??
?(中國(guó)科學(xué)院力學(xué)研究所,非線性力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100190)
?(中國(guó)科學(xué)院大學(xué)工程科學(xué)學(xué)院,北京100049)
??(中國(guó)醫(yī)學(xué)科學(xué)院北京協(xié)和醫(yī)學(xué)院,北京協(xié)和醫(yī)院骨科,北京100730)
??(中國(guó)科學(xué)院力學(xué)研究所,工程化構(gòu)建與力學(xué)生物學(xué)北京市重點(diǎn)實(shí)驗(yàn)室,北京100190)
全膝關(guān)節(jié)置換術(shù)(total knee arthroplasty,TKA)是目前治療晚期膝關(guān)節(jié)疾病的常用有效手段[1]。近年來(lái),TKA數(shù)量明顯增多[2],假體周?chē)钦垡苍桨l(fā)常見(jiàn),其發(fā)生的概率為0.2%~2.5%,翻修后的發(fā)生率為2%~4%[3-4]。受到人口老齡化、患者術(shù)后活動(dòng)量增加、假體使用壽命變長(zhǎng)等因素的影響,發(fā)生假體周?chē)钦鄣幕颊邤?shù)量明顯上升,且呈年輕化趨勢(shì)[5-6]。臨床發(fā)現(xiàn),術(shù)后的老年人患者在正常行走狀態(tài)下也會(huì)有骨折的發(fā)生。相關(guān)研究報(bào)道稱(chēng):術(shù)后骨折多為類(lèi)似的低能量損傷[7],如跌倒、崴腳、摔倒、站立位屈膝等情況所致,少數(shù)高能量損傷病例發(fā)生于活動(dòng)量大的年輕人。TKA術(shù)后假體周?chē)钦壑校晒羌袤w周?chē)钦圩顬槌R?jiàn),其中股骨遠(yuǎn)端骨折患者30天、6個(gè)月和1年的死亡率分別為8%,24%和27%[8],脛骨假體周?chē)钦叟c髕骨假體周?chē)钦垡矔r(shí)有發(fā)生[9]。
臨床發(fā)現(xiàn)一種特殊類(lèi)型的股骨假體周?chē)钦郏鐖D1所示,骨折線長(zhǎng)度雖不完全相同,但骨折線走向均為股骨遠(yuǎn)端外側(cè)下方至內(nèi)側(cè)上方。據(jù)臨床醫(yī)生粗略統(tǒng)計(jì),該種骨折的發(fā)生概率約為T(mén)KA術(shù)后股骨假體周?chē)钦鄣?0%,是術(shù)后常見(jiàn)的一種股骨遠(yuǎn)端骨折,甚至有些患者雙腿同時(shí)發(fā)生該種骨折。由于骨折發(fā)生在股骨髁附近,常伴有骨質(zhì)疏松等癥狀,且周?chē)延屑袤w存在,處理骨折的傳統(tǒng)方法已不適用,給后期翻修治療造成很大困難。因而探究該種骨折發(fā)生的原因?qū)︻A(yù)防骨折有重要的意義。

圖1 臨床骨折CT圖
為此,一些學(xué)者利用有限元方法對(duì)假體及周?chē)墙M織的應(yīng)力狀態(tài)進(jìn)行了研究。文獻(xiàn)[10-13]對(duì)行走、上樓梯等工況下的假體接觸應(yīng)力進(jìn)行了分析,為假體性能優(yōu)化提供了指導(dǎo);文獻(xiàn)[14-15]模擬對(duì)比了使用不同功能類(lèi)型的假體后膝關(guān)節(jié)應(yīng)力分布情況,為新型假體的設(shè)計(jì)提供了參考依據(jù);文獻(xiàn)[16-18]對(duì)比了不同材料、不同構(gòu)型的假體周?chē)墙M織應(yīng)力分布情況,對(duì)理解、減輕術(shù)后應(yīng)力遮擋問(wèn)題有所幫助;周建華[19]的研究表明,股骨遠(yuǎn)端手術(shù)切跡的深度和厚度與術(shù)后髁上骨折的發(fā)生有一定關(guān)聯(lián),切跡深度應(yīng)避免超過(guò)皮質(zhì)骨與松質(zhì)骨的交界處以減少骨折發(fā)生的可能性;Conlisk等[20]發(fā)現(xiàn)使用有柄假體有助于減少假體周?chē)墙M織應(yīng)力,而骨柄長(zhǎng)度是影響老年骨質(zhì)疏松患者術(shù)后骨折風(fēng)險(xiǎn)的重要因素;Kim等[21]對(duì)計(jì)算機(jī)導(dǎo)航膝關(guān)節(jié)置換手術(shù)的安全性進(jìn)行了研究,發(fā)現(xiàn)定位孔的直徑與孔周?chē)膽?yīng)力呈正相關(guān),是引發(fā)術(shù)后骨折的一個(gè)重要影響因素。目前的研究,從假體材料和構(gòu)型、手術(shù)操作參數(shù)等方面對(duì)假體周?chē)鷳?yīng)力分布情況進(jìn)行了研究,分析了應(yīng)力遮擋情況及骨折發(fā)生風(fēng)險(xiǎn)的影響因素,對(duì)提高假體使用壽命、降低術(shù)后骨折發(fā)生風(fēng)險(xiǎn)具有普遍的指導(dǎo)意義。然而,對(duì)于本文所提到的特殊類(lèi)型的假體周?chē)钦鄣陌l(fā)生原因,尚未有針對(duì)性研究。
為此,本文利用有限元方法針對(duì)站立、屈膝、行走三種典型工況對(duì)TKA術(shù)后股骨遠(yuǎn)端的應(yīng)力狀態(tài)進(jìn)行了模擬分析,尋找造成這種特殊類(lèi)型假體骨折的原因,用于幫助醫(yī)生給TKA術(shù)后患者提供更科學(xué)的預(yù)防建議,也可為假體優(yōu)化設(shè)計(jì)提供參考。
以志愿者的膝關(guān)節(jié)CT圖像數(shù)據(jù)(女性,59歲,左側(cè))為研究對(duì)象。CT圖像的實(shí)質(zhì)是通過(guò)不同角度的射線掃描獲得被掃描體各點(diǎn)的衰減系數(shù)值,膝關(guān)節(jié)CT如圖2(a)所示。骨、肌肉、軟組織等組織結(jié)構(gòu)不同,各組織的衰減系數(shù)值范圍不同。通過(guò)閾值分割、區(qū)域增長(zhǎng)、填充等工作,獨(dú)立出完整的骨組織后進(jìn)行三維模型重建。骨模型建立過(guò)程如圖2(b)和圖2(c)所示。建立膝關(guān)節(jié)假體模型,通過(guò)布爾運(yùn)算進(jìn)行“模擬手術(shù)”,得到裝配后的股骨假體計(jì)算模型如圖2(d)所示。骨材料的泊松比設(shè)為0.3。采用彈塑性模型來(lái)描述骨材料的力學(xué)行為,設(shè)置材料在發(fā)生0.7%的彈性變形時(shí)出現(xiàn)塑性應(yīng)變[22]。通過(guò)選定CT閾值對(duì)骨進(jìn)行分區(qū)域材料賦值,骨組織按照骨髓腔、松質(zhì)骨、皮質(zhì)骨區(qū)分開(kāi),再細(xì)分為十個(gè)區(qū)域。依據(jù)股骨的CT值、骨密度及彈性模量之間的關(guān)系式(1)和式(2),得到每個(gè)骨質(zhì)區(qū)域的彈性模量等相關(guān)參數(shù),賦值后骨模型如圖2(e)所示。

圖2 模型建立過(guò)程
骨密度與CT值的關(guān)系式為

彈性模量與密度的關(guān)系式為

式中,密度ρ單位為kg/m3,CT代表CT值,單位為Hu,彈性模量E單位為MPa。
市面上的膝關(guān)節(jié)假體多為鈷、鉬、鉻合金等材質(zhì),模擬中假體材料設(shè)置為鈷鉻合金,線彈性材料,不考慮塑性變形的發(fā)生,彈性模量為209 GPa,泊松比為0.31,密度為8.9×103kg/m3。股骨與假體之間摩擦系數(shù)設(shè)置為0.3。
人體在不同的運(yùn)動(dòng)狀態(tài)下,膝關(guān)節(jié)的受力情況都會(huì)有相應(yīng)的改變。雙腿站立時(shí),關(guān)節(jié)面受力合力的方向沿著下肢力線的方向,作用力在股骨假體面上均勻施加。進(jìn)行屈膝等彎曲動(dòng)作時(shí),膝關(guān)節(jié)處的受力會(huì)明顯增加,假體上的力主要作用在偏后位置,均勻分布。行走狀態(tài)下,單腿承擔(dān)整體身體重量,且因單腿站立,股骨關(guān)節(jié)面上整體受力的方向向內(nèi)傾斜,外側(cè)髁較內(nèi)側(cè)髁受到更多力的作用。對(duì)以上三種工況進(jìn)行模擬,限制股骨上端約束股骨頭上節(jié)點(diǎn)的6個(gè)自由度為0,模型受力圖如圖3所示。

圖3 模型受力示意圖
骨在結(jié)構(gòu)上主要分為皮質(zhì)骨與松質(zhì)骨。股骨髁多為松質(zhì)骨,松質(zhì)骨的破壞與骨小梁的微損傷有關(guān),多用應(yīng)變來(lái)表征[23]。皮質(zhì)骨分布于骨干,為承力的主要組織,應(yīng)力超過(guò)組織強(qiáng)度時(shí),皮質(zhì)骨可發(fā)生骨折[22];因此本文工作將分別考察松質(zhì)骨的應(yīng)變、皮質(zhì)骨的應(yīng)力狀態(tài)來(lái)分析評(píng)估假體周?chē)钦郯l(fā)生的原因。
圖4為雙腿站立時(shí)的股骨遠(yuǎn)端力學(xué)模擬結(jié)果圖。對(duì)于以松質(zhì)骨為主的股骨髁來(lái)說(shuō),如圖4(a)和圖4(b)所示,塑性應(yīng)變主要出現(xiàn)在髁間窩和假體內(nèi)側(cè)邊緣線附近。參考文獻(xiàn)[22]結(jié)果(松質(zhì)骨的極限應(yīng)變?yōu)?%~2.5%),本文定義應(yīng)變超過(guò)1%的區(qū)域?yàn)楣钦畚kU(xiǎn)區(qū)域,統(tǒng)計(jì)分析發(fā)現(xiàn)外側(cè)髁與內(nèi)側(cè)髁的骨折危險(xiǎn)區(qū)域體積比為1.09:1,無(wú)明顯差異。而對(duì)股骨干部位來(lái)說(shuō),如圖4(c)所示,股骨干上呈現(xiàn)出近乎對(duì)稱(chēng)的兩條應(yīng)力集中線。為量化分析股骨干的應(yīng)力分布情況,以10 MPa為例,統(tǒng)計(jì)發(fā)現(xiàn)股骨干外側(cè)與內(nèi)側(cè)應(yīng)力超過(guò)10 MPa的單元體積比為1.32:1,外側(cè)應(yīng)力較大的單元略多于內(nèi)側(cè)。

圖4 雙腿站立時(shí)股骨遠(yuǎn)端力學(xué)模擬圖
圖5所示為屈膝狀態(tài)下的股骨遠(yuǎn)端力學(xué)模擬結(jié)果。對(duì)于松質(zhì)骨為主的股骨髁,如圖5(a)和圖5(b)所示,塑性應(yīng)變主要出現(xiàn)在髁間窩位置。統(tǒng)計(jì)發(fā)現(xiàn)外側(cè)與內(nèi)側(cè)的骨折危險(xiǎn)區(qū)域體積比為1.46:1,比站立狀態(tài)差異明顯。股骨干的應(yīng)力如圖5(c)所示,分布相對(duì)均勻,統(tǒng)計(jì)股骨干遠(yuǎn)端股骨干外側(cè)、內(nèi)側(cè)應(yīng)力超過(guò)10 MPa的單元體積,比例為1.31:1,外側(cè)應(yīng)力較大的單元略多于內(nèi)側(cè)。

圖5 屈膝狀態(tài)股骨遠(yuǎn)端力學(xué)模擬結(jié)果圖
圖6所示為行走狀態(tài)下的股骨遠(yuǎn)端力學(xué)模擬結(jié)果。對(duì)于以松質(zhì)骨為主的股骨髁,如圖6(a)和圖6(b)所示,塑性應(yīng)變主要集中在髁間窩外側(cè)和假體外側(cè)邊緣線附近,假體外側(cè)邊緣線附近的塑性應(yīng)變區(qū)呈現(xiàn)向內(nèi)側(cè)上方擴(kuò)展的走勢(shì),這和臨床上發(fā)現(xiàn)的骨折線呈現(xiàn)由外側(cè)下方向內(nèi)側(cè)上方走勢(shì)吻合。統(tǒng)計(jì)發(fā)現(xiàn),外側(cè)與內(nèi)側(cè)的骨折危險(xiǎn)區(qū)域體積比為4.33:1,遠(yuǎn)高于站立狀態(tài)和屈膝狀態(tài),說(shuō)明行走時(shí)股骨髁外側(cè)發(fā)生骨折的風(fēng)險(xiǎn)明顯增大。而對(duì)于股骨干部位來(lái)說(shuō),外側(cè)應(yīng)力明顯大于內(nèi)側(cè),統(tǒng)計(jì)遠(yuǎn)端股骨干外側(cè)、內(nèi)側(cè)應(yīng)力超過(guò)10 MPa的單元體積比為3.40:1,有顯著差異,且大應(yīng)力區(qū)集中在外側(cè)靠近股骨髁的部位,這會(huì)增大假體周?chē)钦塾赏鈧?cè)股骨髁部位發(fā)生的風(fēng)險(xiǎn)。

圖6 行走態(tài)股骨遠(yuǎn)端力學(xué)模擬結(jié)果圖
本文對(duì)站立、屈膝、行走工況下股骨遠(yuǎn)端假體周?chē)鷳?yīng)力和應(yīng)變情況進(jìn)行了分析,結(jié)果表明:行走工況下,股骨髁部位的塑性應(yīng)變主要集中在髁間窩外側(cè)和假體外側(cè)邊緣線附近,假體外側(cè)邊緣線附近的塑性應(yīng)變區(qū)呈現(xiàn)向內(nèi)側(cè)上方擴(kuò)展的走勢(shì),這和臨床上發(fā)現(xiàn)的骨折線呈現(xiàn)由外側(cè)下方向內(nèi)側(cè)上方走勢(shì)吻合,是誘發(fā)該類(lèi)型骨折的重要原因。統(tǒng)計(jì)發(fā)現(xiàn),外側(cè)與內(nèi)側(cè)的骨折危險(xiǎn)區(qū)域體積比遠(yuǎn)高于站立狀態(tài)和屈膝狀態(tài),說(shuō)明行走時(shí)股骨髁外側(cè)發(fā)生骨折的風(fēng)險(xiǎn)明顯增大。
而對(duì)于股骨干部位來(lái)說(shuō),行走時(shí)外側(cè)應(yīng)力明顯大于內(nèi)側(cè),且大應(yīng)力區(qū)集中的外側(cè)靠近股骨髁部位,這是導(dǎo)致假體周?chē)钦塾赏鈧?cè)股骨髁部位發(fā)生的另一個(gè)重要原因。