膝關(guān)節(jié)有限元模型屬性設(shè)置

膝關(guān)節(jié)有限元模型：憑仗有限元解析法，將膝關(guān)節(jié)解剖模型離散化為由有限個(gè)單元網(wǎng)格構(gòu)成的模型，膝關(guān)節(jié)有限元模型是研發(fā)膝關(guān)節(jié)生物 力學(xué)的一類(lèi)有效途徑。膝關(guān)節(jié)模型屬性設(shè)置：膝關(guān)節(jié)模型中最首要的4個(gè)方面的屬性分別為模型如果、資料賦值、網(wǎng)格區(qū)分并且載荷加以。

摘要

布景：緊隨膝關(guān)節(jié)有關(guān)研發(fā)的深入，對(duì)于膝關(guān)節(jié)有限元模型創(chuàng)建時(shí)屬性選用的確切性及應(yīng)用性仍存在絕對(duì)的爭(zhēng)論，亟待研發(fā)。

目標(biāo)：綜合膝關(guān)節(jié)有限元模型的文獻(xiàn)，為模型建立和解析供應(yīng)更合理的屬性設(shè)置，并提高其應(yīng)用性。

方式：對(duì)于9組海內(nèi)近5年膝關(guān)節(jié)有限元模型文獻(xiàn)中骨性與非骨性組織的模型如果、資料賦值、網(wǎng)格區(qū)分并且載荷加以4個(gè)方面進(jìn)行對(duì)照， 結(jié)合外文文獻(xiàn)及海外研發(fā)，驗(yàn)證其合感性。

結(jié)果與論斷：①模型如果：骨骼設(shè)剛體，韌帶設(shè)超彈性，半月板設(shè)橫向各向異性，軟骨設(shè)單相各向異性線(xiàn)彈性；②資料賦值：骨骼無(wú)需賦 值，韌帶賦剪切模量，半月板辨別軸向、徑向和圓周方向賦值，軟骨賦彈性模量；③網(wǎng)格區(qū)分：骨骼用殼單元，非骨性組織用8節(jié)點(diǎn)六面 體或二階四面體單元；④載荷加以：在股骨考慮點(diǎn)加以1倍體品質(zhì)軸向載荷；⑤通過(guò)對(duì)照獲得的屬性設(shè)置能讓膝關(guān)節(jié)有限元模型的應(yīng)用性 更強(qiáng)，為以后進(jìn)一步研發(fā)供應(yīng)愈加確切的模仿結(jié)果。

膝關(guān)節(jié)；有限元模型；資料賦值；模型如果；網(wǎng)格區(qū)分；載荷加以

引言 Introduction 膝關(guān)節(jié)是身體最高難的關(guān)節(jié)之一，位于髖關(guān)節(jié)與踝關(guān)節(jié) 之間，是下肢活動(dòng)的樞紐。緊隨社會(huì)人口增加，老齡化形勢(shì) 逐漸突出，膝關(guān)節(jié)有關(guān)重病如關(guān)節(jié)炎、半月板傷害等的發(fā)病 率也逐年提升。據(jù)預(yù)計(jì)，美國(guó)有 5 000 萬(wàn)膝骨關(guān)節(jié)炎患者， 在超越60歲人群中每年有5%的人因患此病而喪失舉措本領(lǐng)， 而據(jù)海內(nèi)統(tǒng)計(jì)顯現(xiàn)，海內(nèi)約有 3% 的人患有骨關(guān)節(jié)炎，膝骨 關(guān)節(jié)炎占大一些比率 [1]。

膝關(guān)節(jié)有限元解析成為一類(lèi)可以直觀(guān)反映膝關(guān)節(jié)內(nèi)部應(yīng) 力散布的方式，從 1972 年 BREKELMANS 等 [2] 和 RYBICKI 等 [3] 將限元法初次應(yīng)用于股骨的應(yīng)力解析此后，距今已有近半個(gè) 時(shí)代的時(shí)間。緊隨有限元解析的全面推行，其在研究膝關(guān)節(jié) 病患生物力學(xué)體制中的應(yīng)用也呈現(xiàn)出了廣大的遠(yuǎn)景。近些年 來(lái)，海內(nèi)對(duì)于膝關(guān)節(jié)有限元模型的研發(fā)在骨骼、韌帶、半月 板并且關(guān)節(jié)軟骨的資料賦值、模型如果、網(wǎng)格區(qū)分并且載荷 加以這 4 個(gè)方面存在著絕對(duì)的差別。此文對(duì)于 2015 至 2020 年 5 年間文獻(xiàn)報(bào)導(dǎo)的 9 組膝關(guān)節(jié)有限元模型，選取站立位狀 態(tài)即屈膝 0° 時(shí)的屬性進(jìn)行以上 4 個(gè)方面的橫向?qū)φ?，結(jié)合 有關(guān)資料屬性的外文文獻(xiàn)出處，通過(guò)對(duì)照得出能提高建模效 率、減小計(jì)算量并且提升膝關(guān)節(jié)模型應(yīng)用性的建模方式，為 膝關(guān)節(jié)有限元模型的建立并且后期研發(fā)供應(yīng)絕對(duì)的根據(jù)。

1 材料和方式 Data and methods

1.1 設(shè)計(jì) 文獻(xiàn)屬性對(duì)照解析，膝關(guān)節(jié)有限元模型創(chuàng)建及參 數(shù)設(shè)置。

1.2 文獻(xiàn)選取 查詢(xún) 2015 至 2020 年間，下載量靠前的中文 文獻(xiàn)，通過(guò)篩選獲得膝關(guān)節(jié)模型創(chuàng)建相對(duì)完好的 9 篇文獻(xiàn)。通過(guò)中文文獻(xiàn)中引用的外文文獻(xiàn)屬性進(jìn)行查找，對(duì)照有關(guān)外 文文獻(xiàn)中引用源文獻(xiàn)的差別；對(duì)照在膝關(guān)節(jié)有限元模型中模 型如果、資料賦值、網(wǎng)格區(qū)分并且載荷加以 4 個(gè)方面的屬性。1.3 膝關(guān)節(jié)有限元模型案例對(duì)照

1.3.1 骨骼、韌帶、半月板和關(guān)節(jié)軟骨的模型如果 骨骼的 模型如果首要?jiǎng)潪閯傮w和線(xiàn)彈性?xún)深?lèi) [4-12]。在膝關(guān)節(jié)有限 元解析中，將骨骼簡(jiǎn)化為各向異性線(xiàn)彈性資料能在保留一 定扭曲量的同時(shí)簡(jiǎn)化后期的計(jì)算量?？墒?，骨骼硬度宏大 于非骨性組織 [13]，況且其在全部膝關(guān)節(jié)模型中占比最大。DONAHUE 等 [14] 在研發(fā)中證明將骨骼如果為剛體時(shí)碰觸變量 的改變小過(guò) 2%，單元數(shù)目能減小約 5 000 個(gè)，計(jì)算時(shí)間減小 50%。表 1 中陳文棟等 [9] 為了不讓解析結(jié)果構(gòu)成受力彌散現(xiàn) 象，將骨骼如果為彈性后，還單獨(dú)分隔股骨殘端與脛骨殘端 一小片層將其設(shè)為剛體。

韌帶的模型如果首要?jiǎng)潪槌瑥椥院途€(xiàn)彈性?xún)深?lèi) [4-12]。韌帶是軟組織，其擁有非線(xiàn)性、各向同性并且黏彈性的特征 [15]， 將韌帶如果成各向異性線(xiàn)彈性資料能夠減小運(yùn)算量減低建模 難度，但因?yàn)檫M(jìn)行了較大簡(jiǎn)化而使結(jié)果確切性不高。相比之 下，超彈性資料具有很好的伸縮性和恢復(fù)性以及資料特征和 幾何特點(diǎn)都呈非線(xiàn)性改變 [16]，將韌帶如果為超彈性更適合 其力學(xué)功能，結(jié)果也更牢靠。與此同時(shí)，超彈性資料應(yīng)變能 密度函數(shù)本構(gòu)模型劃為 Neo-Hookean 模型、多項(xiàng)式模型、 Mooney-Rivlin 模型等 [17]，特別是 Neo-Hookean 模型能夠在 確保模仿確切性的同時(shí)顯著簡(jiǎn)化須要設(shè)置的資料屬性。半月板的模型如果首要?jiǎng)潪闄M向各向異性和線(xiàn)彈性?xún)?種 [4-12]。

半月板和關(guān)節(jié)軟骨都屬于含水軟骨組織，首要由固相 和液相構(gòu)成 [18]，其成為膝關(guān)節(jié)模仿中首要應(yīng)力散布領(lǐng)域須要 采取小尺寸網(wǎng)格區(qū)分以確保模仿精度，采取各向異性線(xiàn)彈性假 設(shè)能夠在網(wǎng)格數(shù)目加大的條件下減小運(yùn)算量。相比之下，使 用橫向各向異性如果能表現(xiàn)半月板固相基質(zhì)的各向同性 [19]， 以及 HAUT 等 [20] 證明了壓縮載荷下要在脛骨平臺(tái)上實(shí)行常態(tài) 碰觸壓力散布，采取橫向各向異性如果是必須的。

關(guān)節(jié)軟骨模型如果首要?jiǎng)潪閱蜗喔飨虍愋跃€(xiàn)彈性和各 向異性線(xiàn)彈性?xún)深?lèi) [4-12]。因?yàn)殛P(guān)節(jié)軟骨自身也擁有雙相性， MOW 等 [21] 認(rèn)定雙相性理論是唯獨(dú)可以描繪大多數(shù)能夠觀(guān) 察到的關(guān)節(jié)軟骨扭曲言行的流變模型；而 GARCIA 等 [22] 認(rèn) 為關(guān)節(jié)軟骨簡(jiǎn)化為單相性在短加載時(shí)間內(nèi)所構(gòu)成的邊緣牽引 力和雙相性如果是相近；DONZELLI 等 [23] 也發(fā)掘在載荷加以 后的很短時(shí)間內(nèi)，關(guān)節(jié)軟骨的碰觸反應(yīng)沒(méi)有顯著的改變。DONAHUE 等 [14] 認(rèn)定關(guān)節(jié)軟骨組織內(nèi)的流體在較短的加載時(shí) 間內(nèi)沒(méi)有產(chǎn)生流動(dòng)，因此關(guān)節(jié)軟骨線(xiàn)彈性如果要優(yōu)于雙相性。LI 等 [24] 并且 PE?A 等 [25] 在關(guān)節(jié)軟骨的模型如果中應(yīng)用的是 單相各向異性線(xiàn)彈性如果。相較于線(xiàn)彈性和雙相性，在近期 呼應(yīng)狀況下，采取單相各向異性如果能夠進(jìn)行有效的碰觸應(yīng) 力模仿。

1.3.2 骨骼、韌帶、半月板和關(guān)節(jié)軟骨的資料賦值 見(jiàn)表 2， 3。

骨骼的資料賦值有 3 種：第一類(lèi)是剛體 [4-5，7，10-11]，其 不須要任意賦值；第二種是線(xiàn)彈性 [6，8，12]，ASHMAN 等 [26] 通 過(guò)對(duì) 60 具身體股骨樣件進(jìn)行實(shí)驗(yàn)測(cè)驗(yàn)和公式計(jì)算，獲得了 股骨的彈性模量平均值分別為：E1=12 GPa，E2=13.4 GPa，E3= 20 GPa，泊松比：v12=0.376，v13=0.222，v23=0.235。DONAHUE 等[14]在模型中，對(duì)小梁骨賦值彈性模量0.4 GPa和泊松比0.3?；诤Ｍ鈱W(xué)者的測(cè)驗(yàn)和模仿中屬性的選取能夠看出針對(duì)骨骼 賦值的彈性模量值比非骨性組織要高于不少。

韌帶的賦值有 3 種：第 1 種是依據(jù)位子和功效辨別其 賦值的數(shù)值 [5-6，12]；第 2 種是將一切韌帶賦以相近的彈性模量 [7-9]，但如此獲得的模仿結(jié)果會(huì)缺失 確切性；第 3 種是應(yīng)用 Neo-Hookean 模 型的超彈性如果韌帶 [4，10-11]，其須要設(shè) 置初始剪切模量“C1”和資料不能壓縮 屬性“D1”2個(gè)屬性，若資料采取不能 壓縮如果，則只要初始剪切模量“C1”。

針對(duì)剪切模量的賦值在表 3 中首要 能夠差別 2 種：①一類(lèi)是源自 PE?A 等 [25] 基于此前研發(fā)者測(cè)驗(yàn)的信息將韌帶劃為 外側(cè)副韌帶、內(nèi)側(cè)副韌帶、前交叉及后 交叉韌帶并分隔賦值 6.06，6.43，5.83 和 6.06 MPa[27]；②另一類(lèi)是源自下列國(guó) 外學(xué)者的研發(fā)結(jié)果，BUTLER 等 [28] 通過(guò) 對(duì)一具 38 歲男人的髕韌帶、前交叉及 后交叉韌帶樣件進(jìn)行實(shí)驗(yàn)解析獲得了 它們的單軸應(yīng)力應(yīng)變曲線(xiàn)；GARDINER 等 [29] 通過(guò)對(duì) 8 具男人膝關(guān)節(jié)進(jìn)舉措態(tài) 測(cè)驗(yàn)并且有限元模型解析，獲得內(nèi)側(cè) 副韌帶的 C1 平均值為 1.44 MPa；PE?A 等 [13] 擬合了單軸應(yīng)力 - 應(yīng)變曲線(xiàn)和相 關(guān)常數(shù)獲得了內(nèi)外側(cè)副韌帶、前交叉韌 帶及后交叉韌帶的 C1 值為 1.44，1.95 及 3.25 MPa[28，30]，并且相應(yīng)的資料不能壓 縮屬性為0.001 26，0.006 83并且0.004 1。

將韌帶簡(jiǎn)化為線(xiàn)性彈簧進(jìn)行賦值 也能絕對(duì)水平提高模仿的效益，部分國(guó) 外學(xué)者在進(jìn)行膝關(guān)節(jié)力學(xué)仿生時(shí)應(yīng)用了 線(xiàn)性彈簧來(lái)取代韌帶 [5，31-32]，如此能在 確保韌帶力學(xué)功能的同時(shí)盡也許讓模型 的屬性設(shè)置獲得簡(jiǎn)化。

半月板的賦值首要有 2 種：①第 一類(lèi)是彈性模量 59 MPa[4-9，11-12]、泊松 比 0.49[4-8，11-12]，賦值源自 LEROUX 等 [19] 通過(guò)對(duì)10具狗的半月板進(jìn)行拉伸實(shí)驗(yàn)， 獲得了半月板在纖維方向及垂直于纖維方向的楊氏模量為 67.8 和 11.1 MPa， 泊松比在各向異性的限定下不宜該超越 0.5。PE?A 等 [25] 基 于 LEROUX 等 [19] 研 究的根基上，在模型中將半月板如果為 了各向異性線(xiàn)彈性，并賦值平均參數(shù)彈 性模量 59 MPa 并且泊松比 0.49。②第 二種是將半月板劃為軸向、徑向和圓周 方向分隔賦值 [10]。因?yàn)榘朐掳逯心z原 纖維 ( Ⅰ型 ) 首要沿圓周方向擺列，從 而使半月板在該方向上比在徑向上更 硬 [33]。TISSAKHT 等 [34] 對(duì) 31 具身體半 月板樣件進(jìn)行拉伸測(cè)驗(yàn)，獲得徑向和周 向樣件的前部、中部、后部的彈性模 量分別為 7.82，11.49，13.04 MPa 并且 99.75，90.22，102.12 MPa。DONAHUE 等 [14] 和 YANG 等 [35] 在進(jìn)行模仿時(shí)將半月板的軸向和徑向賦值 20 MPa，圓周 方向辨別賦值 140 和 120 MPa。針對(duì)泊 松比的賦值，KEMPTON[36] 通過(guò)測(cè)驗(yàn)和 對(duì)照信息發(fā)掘半月板的泊松比在 0.4- 1.0 之間，以及泊松比在樣品表層到達(dá) 最大值。ASPDEN[37] 在模仿壓縮載荷對(duì) 半月板的傷害狀況時(shí)，用到的泊松比也 是 0.2 和 0.4。

關(guān)節(jié)軟骨的賦值首要有 2 種：第 一類(lèi)是彈性模量 5 MPa[4，7-8，10-12]，泊松 比 0.46[4，7-8，10-12]。LI 等 [24] 在膝關(guān)節(jié)模 擬中對(duì)關(guān)節(jié)軟骨賦值彈性模量 5 MPa 和 泊松比 0.45 以代表常態(tài)行走流程中的 軟骨言行。第二種是彈性模量 15 MPa 及泊松比 0.3[5]。ARMSTRONG 等 [38] 提 出關(guān)節(jié)軟骨成為雙相資料其緊隨加載時(shí)

https://mp.weixin.qq.com/s/9mvsqOQ-WHHM67Sk-RCX1g