仵健磊 彭 偉# 董輝躍 姜獻(xiàn)峰 劉云峰*

1(浙江工業(yè)大學(xué)特種裝備制造與先進(jìn)加工技術(shù)教育部/浙江省重點(diǎn)實(shí)驗(yàn)室，杭州 310014)2(浙江大學(xué)浙江省先進(jìn)制造技術(shù)重點(diǎn)實(shí)驗(yàn)室，杭州 310027)

引言

牙周膜是連接在牙齒與牙槽骨之間的一層結(jié)締纖維組織，具有支持牙齒，傳遞、吸收和緩沖咬合力的作用[1-2]。臨床正畸治療中，對(duì)于錯(cuò)位牙齒的矯正，牙周膜起著至關(guān)重要的作用。一般來(lái)說(shuō)，由矯治器產(chǎn)生的矯治力首先傳遞給牙齒，再經(jīng)牙周膜緩沖之后作用在牙槽骨上，從而引起牙槽骨改建等一系列生物反應(yīng)，并最終在矯治力的作用下，牙齒產(chǎn)生相應(yīng)的移動(dòng)和轉(zhuǎn)動(dòng)，達(dá)到正畸的目的。在此矯治過(guò)程中，牙周膜起著重要作用，越來(lái)越多的學(xué)者認(rèn)為牙槽骨的改建與牙周膜的應(yīng)力和應(yīng)變息息相關(guān)[3-4]。因此，需要針對(duì)牙周膜開(kāi)展更加深入的研究。

牙周膜作為人體生物組織，實(shí)驗(yàn)材料來(lái)源有限，只能從捐獻(xiàn)尸體的口腔中獲取；另一方面，牙周膜為活性組織且極易遭到破壞，正常人體牙周膜厚度為0.2～0.3 mm[5]；牙周膜與其他人體口腔組織分離之后將失去生物活性，所以牙周膜的保存條件也極為苛刻。在實(shí)驗(yàn)條件有限的情況下，更多的研究人員采用有限元的方法對(duì)牙周膜進(jìn)行相關(guān)研究。目前，有限元法已經(jīng)在口腔正畸乃至整個(gè)醫(yī)療領(lǐng)域得到了廣泛應(yīng)用，具有準(zhǔn)確度高、重復(fù)性好、可操作性強(qiáng)等特點(diǎn)，越來(lái)越受到正畸領(lǐng)域人員的青睞[6-9]。

為了更好地研究矯治力-牙周膜應(yīng)力應(yīng)變-牙槽骨改建-牙齒移動(dòng)之間的復(fù)雜作用關(guān)系，需要建立合適的牙周膜本構(gòu)模型。在之前的許多有限元工作中，牙周膜多采用線彈性材料模型[10-11]，這顯然與實(shí)際不符。研究表明，牙周膜的本構(gòu)關(guān)系曲線表現(xiàn)出很強(qiáng)的指數(shù)形式[12]。為此，黃輝祥等構(gòu)建了牙周膜的超彈性模型[1-2]，董晶等將牙周膜設(shè)定為雙線性模型進(jìn)行有限元分析[13]，魏志剛等分別將牙周膜設(shè)定為線彈性、超彈性及次彈性模型進(jìn)行研究[14]。但是，相比于彈性材料，牙周膜更貼近黏彈性材料，其黏性作用不能忽視[15]。魏志剛等同時(shí)又構(gòu)建了牙周膜的黏彈性模型，但并沒(méi)有對(duì)黏彈性模型的有效性進(jìn)行合理的驗(yàn)證[16]；國(guó)外學(xué)者NataliA等和Iman Z等分別構(gòu)建了牙周膜的黏彈性模型，但是其模型數(shù)據(jù)來(lái)源為豬和牛，并沒(méi)有對(duì)人體牙周膜進(jìn)行直接研究，其本構(gòu)模型參數(shù)并不適用于人體牙周膜[4,17]。

為了構(gòu)建一種準(zhǔn)確度高、實(shí)用性強(qiáng)的牙周膜本構(gòu)模型，以人體牙周膜平面剪切和應(yīng)力松弛實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ)，利用有限元ABAQUS中自帶的非線性材料模型，通過(guò)實(shí)驗(yàn)數(shù)據(jù)擬合確定符合牙周膜特性的超-黏彈性本構(gòu)模型及參數(shù)，并通過(guò)對(duì)牙周膜平面剪切實(shí)驗(yàn)過(guò)程的模擬來(lái)驗(yàn)證牙周膜本構(gòu)模型的正確性，最后通過(guò)對(duì)比牙周膜線彈性與超-黏彈性模型，分析不同牙周膜本構(gòu)模型力學(xué)響應(yīng)的差異，為口腔正畸生物力學(xué)研究提供有力的理論支持。

1 材料和方法

1.1 人體牙周膜實(shí)驗(yàn)數(shù)據(jù)

為了構(gòu)建更加準(zhǔn)確的牙周膜本構(gòu)模型，首先需要可靠性好的牙周膜實(shí)驗(yàn)數(shù)據(jù)。牙周膜作為人體生物活性組織，其實(shí)驗(yàn)材料的獲取受到很大程度的限制，經(jīng)過(guò)查閱大量文獻(xiàn)發(fā)現(xiàn)，Toms等的牙周膜實(shí)驗(yàn)數(shù)據(jù)具有較高的可信度，這里將以此實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ)進(jìn)行牙周膜模型的構(gòu)建[18]。

Toms等的實(shí)驗(yàn)牙周膜對(duì)象取自一位男性(年齡78歲)尸體的下頜前磨牙，針對(duì)前磨牙牙根進(jìn)行0.85 mm厚度的組織切片處理，最終獲得了A1～A5共5組實(shí)驗(yàn)試樣。牙根組織切片試樣包括牙根、牙周膜以及牙槽骨，試樣的橫斷面垂直于牙長(zhǎng)軸方向，如圖1所示。

針對(duì)牙根組織切片試樣，分別進(jìn)行了剪切實(shí)驗(yàn)和應(yīng)力松弛實(shí)驗(yàn)，實(shí)驗(yàn)測(cè)試簡(jiǎn)圖如圖 2所示，實(shí)驗(yàn)過(guò)程中牙槽骨固定，對(duì)切片的中間牙根施加位移載荷，壓頭下壓速度為0.2 mm/min。實(shí)驗(yàn)中以牙周膜的應(yīng)變角γ作為應(yīng)變，以壓頭受到的反作用力F與牙周膜圓周面積As的比值作為應(yīng)力，試樣厚度為d，牙周膜的寬度為WL，壓頭的位移量為Δl。

經(jīng)過(guò)對(duì)5組牙根組織切片的實(shí)驗(yàn)測(cè)試，獲得了牙周膜剪切實(shí)驗(yàn)應(yīng)力應(yīng)變曲線和應(yīng)力松弛曲線，Toms等通過(guò)對(duì)實(shí)驗(yàn)曲線的處理進(jìn)一步獲得了牙周膜剪切應(yīng)力應(yīng)變曲線方程和歸一化應(yīng)力松弛曲線方程，如表1所示。接下來(lái)，將以該實(shí)驗(yàn)數(shù)據(jù)作為牙周膜本構(gòu)模型及參數(shù)構(gòu)建的數(shù)據(jù)基礎(chǔ)。

圖2 實(shí)驗(yàn)測(cè)試簡(jiǎn)圖

表1 5組牙周膜試樣的剪切應(yīng)力-應(yīng)變曲線方程和歸一化應(yīng)力松弛曲線方程

1.2 牙周膜本構(gòu)模型

牙周膜既不屬于符合胡克定律的彈性體，也不屬于非牛頓流體，而是介于兩者之間的一種高黏彈性體，其力學(xué)模型及參數(shù)與彈性體和流體均不相同。黏彈性材料的流動(dòng)變形包括兩個(gè)階段：在初始階段首先表現(xiàn)為一個(gè)瞬間的彈性變形，隨后是持續(xù)一段時(shí)間的黏性變形，在此期間保持應(yīng)變不變的情況下，應(yīng)力會(huì)隨著時(shí)間而逐步衰減。牙周膜作為一種黏彈性材料，在受載荷的一瞬間表現(xiàn)為超彈性，隨后又表現(xiàn)為黏彈性，因此需要分別定義牙周膜的超彈性力學(xué)行為和黏彈性力學(xué)行為，即超-黏彈性行為[1]，其中超彈性模型的模量時(shí)間尺度為瞬態(tài)。

1.2.1超彈性本構(gòu)模型

基于連續(xù)介質(zhì)力學(xué)理論，多項(xiàng)式形式模型是最早提出的一種經(jīng)典本構(gòu)模型。對(duì)于各向同性材料，應(yīng)變能密度可以分解為應(yīng)變偏量能和體積應(yīng)變能兩部分，形式如下：

(1)

式中，f和g分別表示應(yīng)變偏量能和體積應(yīng)變能。

將式(1)進(jìn)行泰勒展開(kāi)，可得

(2)

(3)

一般來(lái)說(shuō)，縮減多項(xiàng)式模型的階數(shù)N越高，其適用性越廣。在有限元軟件ABAQUS中縮減多項(xiàng)式模型的最高階數(shù)為六階。這里，將采用縮減多項(xiàng)式模型描述牙周膜的超彈性行為。

1.2.2黏彈性本構(gòu)模型

(4)

如果將多個(gè)Maxwell模型并聯(lián)或多個(gè)Kelvin模型串聯(lián)起來(lái)，就得到了廣義Maxwell模型和廣義Kelvin模型，后者又稱(chēng)為Kelvin鏈。對(duì)于Kelvin鏈而言，其本構(gòu)模型為

(5)

即

(6)

式(5)、(6)即為線性黏彈性本構(gòu)方程的一般表達(dá)。對(duì)于不可壓縮各向同性材料，體積應(yīng)力只改變體積，應(yīng)力偏量導(dǎo)致等體積的形狀畸變，而且兩種情形下的黏彈特性與效應(yīng)可分別考慮。因而其本構(gòu)關(guān)系可以表示為

P′·Sij=Q′·eij，P″·σkk=Q″·εkk

(7)

對(duì)于大變形黏彈性，又可將其黏彈性本構(gòu)模型改寫(xiě)成應(yīng)變能形式，其Prony形式為

(8)

式中，W0為瞬時(shí)應(yīng)變能函數(shù)，W∞為穩(wěn)態(tài)應(yīng)變能函數(shù)，δn為能量函數(shù)乘子，λn為松弛時(shí)間常數(shù)，N為階數(shù)，t為時(shí)間。

由式(8)，可得名義應(yīng)力為

(9)

式中，xI為變形狀態(tài)時(shí)的坐標(biāo)，XJ為初始狀態(tài)時(shí)的坐標(biāo)。

假設(shè)W∞=0，則

(10)

由此可得牙周膜黏彈性的歸一化應(yīng)力松弛函數(shù)為

(11)

式(11)為目前有限元軟件ABAQUS中用來(lái)描述材料黏彈性特性的表達(dá)式。

1.3 牙周膜模型數(shù)據(jù)擬合

為了確定牙周膜本構(gòu)模型的具體形式以及模型參數(shù)，需要基于ABAQUS中現(xiàn)有縮減多項(xiàng)式模型和歸一化應(yīng)力松弛函數(shù)模型進(jìn)行人體牙周膜實(shí)驗(yàn)數(shù)據(jù)擬合，其中牙周膜平面剪切實(shí)驗(yàn)的擬合范圍是0～0.6 rad，應(yīng)力松弛實(shí)驗(yàn)的擬合范圍是0～60 s，分別針對(duì)試樣A1～A5進(jìn)行配準(zhǔn)擬合，選擇擬合度最高的材料模型作為牙周膜的本構(gòu)模型。

1.4 牙周膜剪切實(shí)驗(yàn)的有限元模擬

為驗(yàn)證數(shù)據(jù)擬合獲得的牙周膜本構(gòu)模型及參數(shù)，將針對(duì)牙周膜平面剪切實(shí)驗(yàn)過(guò)程進(jìn)行有限元模擬，通過(guò)對(duì)比模擬曲線與實(shí)驗(yàn)曲線驗(yàn)證模型的正確性。為簡(jiǎn)化分析，假設(shè)牙根組織切片試樣近似為軸對(duì)稱(chēng)圓盤(pán)結(jié)構(gòu)，取其圓周向截面作為研究對(duì)象，并以此建立二維平面軸對(duì)稱(chēng)有限元模型，如圖3所示(牙槽骨模型未畫(huà)出)。牙周膜的外側(cè)設(shè)置為固定約束，牙周膜的最大應(yīng)變角為0.6 rad，屬性定義采用本文第1.3節(jié)擬合獲得的牙周膜模型及參數(shù)；牙根設(shè)定為彈性體，彈性模量為20 000 MPa，泊松比0.3[10]，其上端面以0.2 mm/min的速度施加位移載荷，牙根與牙周膜之間設(shè)定為黏結(jié)約束。因模型為軸對(duì)稱(chēng)模型，同時(shí)為了保證計(jì)算精度，牙齒和牙周膜單元類(lèi)型均采用四邊形雜交軸對(duì)稱(chēng)單元CAX8H。試樣A1～A5的建模數(shù)據(jù)如表2所示，網(wǎng)格單元總數(shù)分別為900，1 100，1 100，1 100和1 100。有限元仿真中，以牙周膜的應(yīng)變角γ作為剪切應(yīng)變，剪切應(yīng)力則采用牙根上端面反作用力RF之和與牙根切片厚度d的比值。

圖3 實(shí)驗(yàn)?zāi)M的有限元模型

1.5 不同牙周膜模型力學(xué)響應(yīng)的對(duì)比

為探究牙周膜超-黏彈性模型在正畸治應(yīng)用的力學(xué)響應(yīng)，以牙周膜線彈性模型為對(duì)比，對(duì)不同牙周膜模型對(duì)正畸力的傳遞、緩沖效果進(jìn)行有限元分析。為方便有限元分析計(jì)算，簡(jiǎn)化牙齒、牙周膜為二維平面模型(牙槽骨省略)，如圖4所示。牙周膜

表2 實(shí)驗(yàn)?zāi)M中試樣A1～A5的有限元參數(shù)

與牙齒為黏結(jié)約束，牙周膜的一側(cè)固定約束，牙根的一側(cè)施加位移載荷，總位移量為0.1 mm，加載時(shí)長(zhǎng)1 s，為獲得較高的計(jì)算精度，牙齒與牙周膜網(wǎng)格單元均采用八結(jié)點(diǎn)雙向二次平面應(yīng)力四邊形單元CPS8R，其中牙根網(wǎng)格單元數(shù)為400，牙周膜網(wǎng)格單元數(shù)為120。牙周膜厚度取0.25 mm，牙周膜線彈性模型的彈性模量取0.68 MPa，泊松比取0.49[18]；牙周膜超-黏彈性模型取自上述研究結(jié)果，材料參數(shù)參照A1組試樣。

圖4 用于牙周膜力學(xué)響應(yīng)對(duì)比的有限元模型

2 結(jié)果

2.1 牙周膜實(shí)驗(yàn)數(shù)據(jù)擬合及模型參數(shù)

以人體牙周膜平面剪切數(shù)據(jù)為基礎(chǔ)，選擇縮減多項(xiàng)式模型對(duì)牙周膜超彈性特性進(jìn)行數(shù)據(jù)擬合，其中試樣A1不同階數(shù)縮減多項(xiàng)式模型的擬合結(jié)果如圖5所示；以牙周膜平面剪切數(shù)據(jù)和應(yīng)力松弛實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ)，以歸一化的應(yīng)力松弛函數(shù)對(duì)牙周膜的黏彈性特性進(jìn)行擬合，試樣A1的黏彈性模型擬合結(jié)果如圖6所示(其他試樣A2～A5擬合曲線與A1基本相同，不再列出)。

從牙周膜超彈性模型的擬合結(jié)果來(lái)看，縮減多項(xiàng)式的階數(shù)越高，擬合曲線越是貼近于實(shí)驗(yàn)曲線，其中六階縮減多項(xiàng)式數(shù)據(jù)擬合獲得的模型參數(shù)列于表3；牙周膜黏彈性模型的擬合曲線與實(shí)驗(yàn)曲線完全吻合，擬合獲得的參數(shù)如表3所示。

圖5 超彈性模型擬合曲線

圖6 黏彈性模型擬合曲線

2.2 牙周膜剪切實(shí)驗(yàn)?zāi)M結(jié)果

通過(guò)不同牙周膜試樣(包括試樣A1～A5)平面剪切實(shí)驗(yàn)過(guò)程的有限元模擬，經(jīng)后處理獲得了牙周膜剪切實(shí)驗(yàn)的應(yīng)力應(yīng)變曲線，如圖7所示。其中，實(shí)驗(yàn)曲線是指Toms等在牙周膜平面剪切實(shí)驗(yàn)中獲得的剪切應(yīng)力應(yīng)變曲線，模擬曲線是指利用本文第2.1節(jié)中構(gòu)建的牙周膜本構(gòu)模型及參數(shù)，通過(guò)有限元對(duì)Toms等的實(shí)驗(yàn)過(guò)程進(jìn)行模擬，獲得的牙周膜剪切應(yīng)力應(yīng)變曲線。

2.3 不同牙周膜模型的力學(xué)響應(yīng)

在相同位移載荷作用下，牙周膜線彈性模型與超-黏彈性模型所產(chǎn)生的力學(xué)響應(yīng)云圖(mises應(yīng)力)如圖8所示；利用后處理手段提取牙周膜底側(cè)

表3 牙周膜本構(gòu)模型曲線擬合參數(shù)

圖7 5組試樣實(shí)驗(yàn)曲線與模擬曲線對(duì)比結(jié)果。(a)試樣A1；(b)試樣A2；(c)試樣A3；(d)試樣A4；(e)試樣A5

圖8 不同牙周膜模型的有限元模擬結(jié)果。(a)線彈性模型；(b)超-黏彈性模型

的反作用力，繪制牙周膜反作用力-牙根位移曲線如圖9所示。

圖9 不同牙周膜模型的力學(xué)響應(yīng)

4 討論

從牙周膜試樣A1平面剪切和應(yīng)力松弛實(shí)驗(yàn)數(shù)據(jù)的擬合結(jié)果圖5可以看出，對(duì)于縮減多項(xiàng)式模型，模型的階數(shù)越高，擬合曲線越是貼近于實(shí)驗(yàn)曲線。在擬合的初始階段，牙周膜應(yīng)變角較小，基本所有階數(shù)的縮減多項(xiàng)式模型均能很好地吻合實(shí)驗(yàn)曲線，但是當(dāng)應(yīng)變角持續(xù)增大到0.3 rad時(shí)，一階和二階縮減多項(xiàng)式模型的擬合效果不夠理想，隨著應(yīng)變角的持續(xù)增加，一階和二階模型不再適用于牙周膜剪切實(shí)驗(yàn)的數(shù)據(jù)擬合。可以看出，低階縮減多項(xiàng)式模型適用于應(yīng)變角較小時(shí)的數(shù)據(jù)擬合，而當(dāng)應(yīng)變角較大時(shí)擬合精度較差。相比低階縮減多項(xiàng)式模型，高階模型與實(shí)驗(yàn)曲線具有更高的吻合度，當(dāng)應(yīng)變角增大到0.6 rad時(shí)，四階、五階和六階縮減多項(xiàng)式模型仍能夠保持與實(shí)驗(yàn)曲線的吻合。但是，在正畸有限元研究工作中，牙周膜的最大應(yīng)變角不可能只有0.6 rad，選擇更高階數(shù)的縮減多項(xiàng)式模型更加有利于準(zhǔn)確描述牙周膜的力學(xué)行為，為此選擇六階縮減多項(xiàng)式模型描述牙周膜的超彈性特性。

對(duì)于牙周膜應(yīng)力松弛實(shí)驗(yàn)的數(shù)據(jù)擬合，從結(jié)果圖6來(lái)看，歸一化的應(yīng)力松弛函數(shù)能夠與松弛實(shí)驗(yàn)曲線吻合，從而準(zhǔn)確描述牙周膜的黏彈性力學(xué)行為。從應(yīng)力松弛的實(shí)驗(yàn)曲線以及模擬曲線可以看出，在整個(gè)應(yīng)力松弛過(guò)程中，當(dāng)應(yīng)變保持不變時(shí)，牙周膜的應(yīng)力是隨著時(shí)間逐漸減小的。在初始階段，應(yīng)力下降的速度較快，隨著時(shí)間推移，應(yīng)力的下降趨勢(shì)逐漸平緩，并逐步趨向于一個(gè)穩(wěn)定值，這完全符合牙周膜的黏彈性特性[18]。

為了驗(yàn)證牙周膜本構(gòu)模型的準(zhǔn)確性，針對(duì)牙周膜平面剪切的實(shí)驗(yàn)過(guò)程進(jìn)行了有限元模擬，從結(jié)果圖7來(lái)看，5組試樣的模擬結(jié)果與實(shí)驗(yàn)結(jié)果基本一致，從而驗(yàn)證了本研究所構(gòu)建的牙周膜是否黏彈性本構(gòu)模型的正確性。另外，從模擬結(jié)果可以看出， A1、A2、A3、A5試樣具有更高的吻合度，A4組試樣的擬合度效果不夠理想，但模擬結(jié)果和實(shí)驗(yàn)結(jié)果均在同一個(gè)數(shù)量級(jí)，且整體趨勢(shì)一致。實(shí)驗(yàn)曲線與擬合曲線之間存在誤差的原因可能是

1)真實(shí)的牙根組織切片有一定的錐度，而在有限元建模時(shí)將其簡(jiǎn)化為軸對(duì)稱(chēng)的圓盤(pán)結(jié)構(gòu)，這種簡(jiǎn)化可能會(huì)對(duì)分析結(jié)果產(chǎn)生不同程度的影響。

2)實(shí)驗(yàn)數(shù)據(jù)取自一位78歲的男性捐獻(xiàn)者，其牙周組織狀況相比于年輕人較差，牙周膜的韌性和抗剪切能力較弱。對(duì)于A4組試樣，其牙周膜厚度為(0.14±0.03)mm，而正常人體牙周膜的組織厚度在0.2～0.3 mm附近[1, 5]，由于其牙周膜薄弱，易受環(huán)境因素的干擾，導(dǎo)致實(shí)驗(yàn)結(jié)果與模擬結(jié)果偏差相對(duì)較大。

3)為簡(jiǎn)化有限元分析，假設(shè)牙根組織切片近似為圓盤(pán)結(jié)構(gòu)，為此只構(gòu)建了實(shí)驗(yàn)試樣的二維模型，實(shí)際上牙根切片更近似于橢圓柱形，構(gòu)建牙根切片的三維橢圓模型可能更加貼近于實(shí)際情況。

以上情況都可能對(duì)模擬結(jié)果產(chǎn)生一定的影響，但總體來(lái)說(shuō)，構(gòu)建的牙周膜超一黏彈性本構(gòu)模型仍具有較高的準(zhǔn)確性。

如圖9所示，針對(duì)不同牙周膜本構(gòu)模型對(duì)載荷力學(xué)響應(yīng)的模擬結(jié)果可知，線彈性模型與超-黏彈性模型存在明顯的差異性，線彈性模型的力學(xué)響應(yīng)近似為線性曲線，超-黏彈性模型的響應(yīng)曲線近似為指數(shù)形式。在位移量為0～0.06 mm時(shí)，兩種模型的力學(xué)響應(yīng)差異較小，超-黏彈性模型曲線略低于線彈性模型曲線，但在位移量超過(guò)0.06 mm時(shí)，超-黏彈性曲線上升趨勢(shì)明顯，線彈性曲線基本保持原來(lái)的曲線斜率上升，兩種模型差異越來(lái)越大；當(dāng)位移量為0.1 mm時(shí)，線彈性模型的牙周膜作用力大小為2.98 N，超-黏彈性模型的為29.92 N，數(shù)量級(jí)差異達(dá)到10倍。由此可知，在牙周膜變形量較小、計(jì)算結(jié)果要求不精確時(shí)，采用線彈性模型可近似替代牙周膜超-黏彈性模型；但在牙周膜變形量稍大時(shí)，線彈性模型力學(xué)響應(yīng)與超-黏彈性模型差異較大，需要構(gòu)建更符合牙周膜特性的超-黏彈性模型[16]。符合牙周膜材料屬性的超-黏彈性本構(gòu)模型，對(duì)于臨床正畸的有限元研究則更為準(zhǔn)確可靠。

此外，牙周膜作為人體生物組織，其材料來(lái)源較少，另一方面牙周膜組織極其薄弱，易受外界因素干擾，這導(dǎo)致人體牙周膜的實(shí)驗(yàn)數(shù)據(jù)較為匱乏。由于只依靠牙周膜的平面剪切數(shù)據(jù)和應(yīng)力松弛數(shù)據(jù)進(jìn)行牙周膜本構(gòu)模型的重建，而缺少牙周膜單/雙向拉壓和體積實(shí)驗(yàn)數(shù)據(jù)，更多類(lèi)型實(shí)驗(yàn)數(shù)據(jù)的結(jié)合更加有利于構(gòu)建準(zhǔn)確度高的牙周膜本構(gòu)模型[2]，因此對(duì)人體牙周膜的實(shí)驗(yàn)以及仿真研究仍有許多工作要做。

5 結(jié)論

1)由六階縮減多項(xiàng)式模型和歸一化應(yīng)力松弛函數(shù)模型組合成的牙周膜超-黏彈性模型，能夠準(zhǔn)確擬合牙周膜的平面剪切和應(yīng)力松弛實(shí)驗(yàn)數(shù)據(jù)。

2)通過(guò)對(duì)牙周膜平面剪切實(shí)驗(yàn)過(guò)程的有限元模擬，可以看出所構(gòu)建的牙周膜本構(gòu)模型具有很高的準(zhǔn)確性。

3)從不同牙周膜本構(gòu)模型的力學(xué)響應(yīng)，可以看出牙周膜應(yīng)變較小時(shí)可假設(shè)為線彈性材料，但應(yīng)變較大時(shí)必須采用更加符合其特性的超-黏彈性模型。

[1] 黃輝祥，湯文成，吳斌，等. 基于超彈性模型的牙周膜生物力學(xué)響應(yīng) [J]. 東南大學(xué)學(xué)報(bào)(自然科學(xué)版)2013, 43(2): 340-344.

[2] 黃輝祥，湯文成，吳斌，等. 基于超彈性模型的牙周膜力學(xué)行為數(shù)值模擬 [J]. 上海交通大學(xué)學(xué)報(bào), 2014, 48(9): 1263-1273.

[3] Kawarizadeh A, Bourauel C, Jager A. Experimental and numerical determination of initial tooth mobility and material properties of the periodontal ligament in rat molar specimens [J]. European Journal of Orthodontics, 2003, 25(6): 569-578.

[4] Natali AN, Pavan PG, Carniel EL, et al. Viscoelastic response of the periodontal ligament: an experimental-numerical analysis [J]. Connective Tissue Research, 2004, 45(4-5): 222-230.

[5] 王洪寧，秦行林，劉東旭，等. 基于micro-CT成像的仿真牙周膜有限元建模探討 [J]. 口腔疾病防治, 2016, 24(5): 283-287.

[6] 仵健磊，劉云峰，彭偉，等. 基于有限元仿真的形狀記憶聚合物弓絲初始正畸力分析 [J]. 中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào), 2016, 35(2), 202-210.

[7] Liao Zhipeng, Chen Junning, Li Wei, et al. Biomechanical investigation into the role of the periodontal ligament in optimising orthodontic force: a finite element case study [J]. Archives of Oral Biology, 2016, 66: 98-107.

[8] Masakazu H, Taiji A, Teruko TY. Computer simulation of orthodontic tooth movement using CT image-based voxel finite element models with the level set method [J]. Computer Methods in Biomechanics and Biomedical Engineering, 2016, 19(5): 474-483.

[9] Anshul C, Maninder SS, Girish C, et al. Evaluation of stress changes in the mandible with a fixed functional appliance: a finite element study [J]. American Journal of Orthodontics and Dentofacial Orthopedics, 2015, 147(2): 226-234.

[10] Qian Yinglian, Fan Yubo, Liu Zhan, et al. Numerical simulation of tooth movement in a therapy period [J]. Clinical Biomechanics, 2008, 23(1): S48-S52.

[11] Hohmann A, Kober C, Young P, et al. Influence of different modeling strategies for the periodontal ligament on finite element simulation result [J]. American Journal of Orthodontics and Dentofacial Orthopedics, 2011, 139(6): 775-783.

[12] Fill TS, Carey JP, Toogood RW, et al. Experimentally determined mechanical properties of, and models for, the periodontal ligament: critical review of current literature [J]. Journal of Dental Biomechanics, 2011, 10:1-9.

[13] 董晶，張哲湛，周?chē)?guó)良. 上頜第一磨牙遠(yuǎn)中移動(dòng)時(shí)牙周膜應(yīng)力分布的非線性三維有限元分析 [J]. 上海口腔醫(yī)學(xué), 2015, 24(3): 315-320.

[14] 魏志剛，湯文成，嚴(yán)斌，等. 基于有限元法的牙周膜本構(gòu)模型研究 [J]. 工程力學(xué), 2009, 26(10): 211-216.

[15] Papadopoulou K, Hasan I, Keilig L, et al. Biomechanical time dependency of the periodontal ligament: a combined experimental and numerical approach [J]. European Journal of Orthodontics, 2013, 35(6): 811-819.

[16] 魏志剛，湯文成，嚴(yán)斌，等. 基于黏彈性模型的牙周膜生物力學(xué)研究 [J]. 東南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 39(3): 484-489.

[17] Iman Z, Oskui, Ata H. Dynamic tensile properties of bovine periodontal ligament: a nonlinear viscoelastic model [J]. Journal of Biomechanics, 2016, 49(5), 756-764.

[18] Toms SR, Dakin GJ, Lemons JE, et al. Quasi-linear viscoelastic behavior of the human periodontal ligament [J]. Journal of Biomechanics, 2002, 35(10): 1411-1415.