膝關(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ū)φ眨Y(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?;诤M鈱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