王沫楠+高鍵+徐雯雯+張凱
摘要:為了避免在股骨頸骨折手術(shù)中人為選擇內(nèi)固定方式產(chǎn)生的不良后果,滿足用PC機(jī)優(yōu)選內(nèi)固定方式的需要,從生物力學(xué)的角度出發(fā),基于有限元法開發(fā)了集合前處理、求解、后處理于一體的股骨頸骨折內(nèi)固定方式優(yōu)選程序,目的是將最優(yōu)固定效果下的緊固螺釘數(shù)量、內(nèi)固定角度及組合方式提供給醫(yī)生施術(shù),該系統(tǒng)在分析了股骨與載荷之間的力學(xué)關(guān)系的基礎(chǔ)上,運(yùn)用參數(shù)化的思想,建立個(gè)體化的股骨有限元模型和骨折內(nèi)固定模型,井應(yīng)用有限元法進(jìn)行模型解算,得到相應(yīng)內(nèi)固定方式下的股骨頭下沉位移、股骨頭橫向位移、股骨頭扭轉(zhuǎn)位移、股骨頭總位移及骨折面之間最大斷裂位移的解算結(jié)果,并根據(jù)評(píng)價(jià)算法優(yōu)選出最合理的內(nèi)固定方案.為臨床選擇內(nèi)固定方式提供有效的定量參考。
關(guān)鍵詞:股骨頸骨折;有限元法;面向?qū)ο?;參?shù)化
DOI:10.15938/j.hust.2016.01.001
中圖分類號(hào):TP301 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1007-2683(2016)01-0001-05
0引言
股骨頸骨折的治療分為保守治療和手術(shù)治療兩種,鑒于保守治療的骨不合率約為48%,發(fā)生股骨頭壞死的概率約為34%,因此現(xiàn)代骨科多提倡積極的手術(shù)治療.手術(shù)治療通常分為骨折復(fù)位、內(nèi)固定方式選擇、穿釘施術(shù)這三個(gè)階段.整個(gè)手術(shù)過(guò)程中除了操作手法由醫(yī)生把握外,內(nèi)固定方式的選擇也是完全憑借醫(yī)生的經(jīng)驗(yàn),常常忽略力學(xué)結(jié)構(gòu)的合理性且不可重復(fù).股骨頸骨折手術(shù)治療過(guò)程中,內(nèi)固定方式的選擇是其中唯一的可控變量,也是在醫(yī)生施術(shù)手法嫻熟的前提下,唯一制約術(shù)后效果的關(guān)鍵因素.所以,選擇最合理的內(nèi)固定方式是穿釘施術(shù)成功的根本保障,本文以股骨頸骨折內(nèi)固定方式優(yōu)選為目標(biāo),將醫(yī)生的操作與理論分析相結(jié)合,從工程學(xué)的角度提出一個(gè)可供醫(yī)生操作的優(yōu)選系統(tǒng).該系統(tǒng)基于有限元法解算出個(gè)體化骨折情況下,內(nèi)固定方式的術(shù)后骨折部位的應(yīng)力、應(yīng)變數(shù)據(jù)結(jié)果,輔助醫(yī)生優(yōu)選出最合理的內(nèi)固定手術(shù)方案.
1股骨頸骨折固定方式優(yōu)選系統(tǒng)方案設(shè)計(jì)
1.1設(shè)置骨折類型與內(nèi)固定方式
緊固螺釘數(shù)量、穿釘角度、螺釘?shù)慕M合方式這三者排列組合后形成了多種內(nèi)固定方式.本系統(tǒng)以股骨頸骨折情況中最常見的經(jīng)頸型骨折為研究對(duì)象,按照臨床治療方法,將骨折角度(0)設(shè)置從30°~70°公差為10°的5種情況;將穿釘角度(B)設(shè)置從30°~70°公差為10°的5種情況;緊固螺釘數(shù)量的選擇為1~3;將組合方式設(shè)置為單釘、2釘斜放、2釘豎放、2釘平放、3釘正放、3釘?shù)狗诺?組合方式.
1.2參數(shù)化的有限元三維建模方式
將骨視為彈性體,并遵守線彈性理論.將CT數(shù)據(jù)讀入醫(yī)學(xué)軟件,建立精確的股骨三維模型,應(yīng)用圣維南定理對(duì)研究對(duì)象進(jìn)行有效的簡(jiǎn)化,確保應(yīng)用彈性力學(xué)的方法來(lái)解決實(shí)際問(wèn)題時(shí)的有效性.對(duì)于股骨來(lái)說(shuō),由于實(shí)際受力環(huán)境復(fù)雜,根據(jù)Scige等提出的較精確的股骨力學(xué)模型,可知在典型的單腿站立位時(shí),股骨受到的力合成后劃分為3種:關(guān)節(jié)反力R、臀肌肌群力M、骼脛束肌力T.
為了實(shí)現(xiàn)參數(shù)化的建模方式,在有限元語(yǔ)言中將頸干角b、股骨頭的直徑d、股骨頸的長(zhǎng)度1,股骨干的直徑d1等患者的解剖學(xué)參數(shù)應(yīng)用*SET函數(shù)設(shè)置為可變參數(shù),具體的值根據(jù)患者的CT圖像測(cè)量得到,由界面輸入,傳遞至有限元代碼中,形成參數(shù)化的有限元代碼文件.
1.3有限元解算
應(yīng)用有限元解算時(shí)本文提出下列假設(shè):
1)解算起始點(diǎn)設(shè)定為骨折部位剛剛完成復(fù)位和內(nèi)固定,骨的愈合過(guò)程還未開始;
2)應(yīng)用圓柱形螺釘固定骨折部位,并假設(shè)螺釘與骨質(zhì)之間是緊密結(jié)合的;
3)假設(shè)圓柱螺釘不會(huì)發(fā)生斷裂,其強(qiáng)度足夠;
4)根據(jù)骨本身的粗糙情況,設(shè)定螺釘與骨之間的摩擦系數(shù)為0.8.
1.4設(shè)置判定準(zhǔn)則
因?yàn)榛颊呤┬g(shù)后只能進(jìn)行輕微的活動(dòng)如緩慢行走等,這是一個(gè)近似靜力分析的過(guò)程,本文考慮其中的單腿站立位.根據(jù)生物力學(xué)原理和AO準(zhǔn)則,從下列參數(shù)人手:
1)股骨頭的位移變化:在單腿站立的情況下,股骨頭上將受到載荷的作用,產(chǎn)生位移,有X方向上反映骨折面抗張力性能和骨折面結(jié)合牢固性的橫向位移UX、Y方向上的反映結(jié)構(gòu)穩(wěn)定性的下沉位移UY、Z方向上的反映抗扭性的扭轉(zhuǎn)位移UZ.這三個(gè)方向的位移值越小表明固定方式越好,結(jié)構(gòu)的抗變形能力越強(qiáng),固定的越牢固。為了綜合評(píng)價(jià)這3個(gè)方向上的位移情況,設(shè)置一個(gè)判斷標(biāo)準(zhǔn)US表示3個(gè)方向上的合成位移,U3越小穩(wěn)定性越好.
2)骨折面的最大斷裂位移:股骨近端是一種類似懸臂梁的結(jié)構(gòu),在外載荷的作用下上側(cè)受拉應(yīng)力,下側(cè)受壓應(yīng)力.假設(shè)骨折面完全斷裂,因此在載荷的作用下會(huì)產(chǎn)生一定量的斷裂位移,這個(gè)位移DS通過(guò)測(cè)量骨折處的兩個(gè)節(jié)點(diǎn)之間的距離得到.DS是反映固定結(jié)構(gòu)的抗彎性,其值越小抗彎性越好。
首先按照U3排序選擇出最優(yōu)U5結(jié)果對(duì)應(yīng)的內(nèi)固定方式,其余內(nèi)固定方式的US值與其的比例關(guān)系是1.047以內(nèi),先選出來(lái),比較每種情況下的DS值與最優(yōu)DS值之間的比例是否在1.17以內(nèi),若在,那么以U3最優(yōu)結(jié)果內(nèi)固定方式為最優(yōu)選項(xiàng),若不在則排除該情況,然后將該種比較方法形成C++算法,供系統(tǒng)優(yōu)選,最終形成的評(píng)價(jià)準(zhǔn)則見表1.
2股骨頸骨折內(nèi)固定方式優(yōu)選系統(tǒng)建立
2.1股骨頸骨折內(nèi)固定方式優(yōu)選系統(tǒng)程序設(shè)計(jì)
步驟1 編寫Visual C++與有限元參數(shù)化語(yǔ)言之間的通信代碼,將用戶界面上搜集的信息傳遞至有限元代碼中實(shí)現(xiàn)建模、解算、結(jié)果實(shí)時(shí)顯示、后處理結(jié)果顯示的功能.
步驟2 根據(jù)獲取的患者信息,如前傾角、頸干角、股骨頭直徑等,將這些信息通過(guò)步驟1的代碼傳遞至有限元代碼中對(duì)應(yīng)的可變參數(shù),建立參數(shù)化的股骨三維模型.
步驟3 依據(jù)患者的病史信息,考慮到患者可能對(duì)內(nèi)固定方式產(chǎn)生影響的疾病,如骨質(zhì)疏松癥.若患者有這種骨質(zhì)疏松的情況,根據(jù)文獻(xiàn)提醒用戶在誤差10′的范圍內(nèi)優(yōu)先選擇單釘?shù)墓潭ǚ绞?