楔狀缺損深度與角度對(duì)下頜第一前磨牙牙體應(yīng)力分布影響的有限元分析

杜珊珊 童錦發(fā) 范中凱

牙體楔狀缺損（wedge-shaped defect，WSD）是一種口腔臨床工作中常見的一種非齲性牙體硬組織疾?。?］，其通常發(fā)生在牙齒的頰側(cè)牙頸部，舌側(cè)較少發(fā)生，中老年人楔型缺損的發(fā)生率較高，＞50歲人群的發(fā)病率約為85%，且男性患者多于女性［2］。楔狀缺損在口腔中并非均勻分布，通常以第一前磨牙最為多見，其次為第二前磨牙［1］。牙體組織的應(yīng)力分布會(huì)因?yàn)樾钊睋p的出現(xiàn)而改變，使應(yīng)力無法沿著牙長(zhǎng)軸均勻連續(xù)傳導(dǎo)，牙折的風(fēng)險(xiǎn)因此增加［3］。臨床上常見的楔狀缺損的患牙可有不同的角度與深度。而不同的缺損深度與角度會(huì)引起牙體應(yīng)力的不同變化。本資料通過三維有限元法建立人工楔狀缺損的下頜前磨牙三維有限元模型，模擬咬合運(yùn)動(dòng)中牙的受力情況，分析多種深度與角度楔狀缺損對(duì)牙體組織應(yīng)力分布的影響，通過研究各種變量之間的相互影響，進(jìn)一步探討楔狀缺損的發(fā)生和變化機(jī)理。

1 材料與方法

1.1 主要設(shè)備及條件 （1）設(shè)備：電腦設(shè)備：Thinkpad CPU，酷睿15-480M，內(nèi)存：4G，硬盤：1T，顯卡：1G。SOMATOM Sensation 16層螺旋CT機(jī)（德國Siemens公司）。（2）相關(guān)軟件：醫(yī)學(xué)建模軟件mimics 20（比利時(shí)Materialise公司），逆向工程軟件Geomagic Studio 2014（美國3D Systems公司），三維機(jī)械制圖專用軟件NX10（德國Siemens公司），有限元分析專用軟件Ansys workbench 18.2（美國ANSYS公司）。

1.2 模型的建立 （1）下頜第一前磨牙三維實(shí)體模型的建立：以一名健康成年女性為觀察對(duì)象，要求其牙列完整，咬合關(guān)系基本正常，且無牙體牙髓牙周疾患，其下頜第一前磨牙符合王惠蕓統(tǒng)計(jì)的我國人牙體解剖標(biāo)準(zhǔn)［4］。采用口腔錐形束CT掃描志愿者的上下頜頜骨，層厚設(shè)置為0.2mm，共計(jì)400張斷層影像，后將文件導(dǎo)出DICOM格式備用。之后，將文件導(dǎo)入至MIMICS20軟件中，根據(jù)不同材質(zhì)灰度的區(qū)別，調(diào)整閾值數(shù)值，提取出右側(cè)下頜第一前磨牙牙本質(zhì)，牙釉質(zhì)及周邊牙槽骨區(qū)域，并且使用編輯工具，分別對(duì)各個(gè)模型進(jìn)行優(yōu)化。利用偏移工具，將牙本質(zhì)向外偏移0.2mm厚，形成牙周膜，將牙周膜向外偏移0.5mm得到固有牙槽骨的模型，將牙槽骨向內(nèi)偏移1.5mm，得到骨松質(zhì)模型，使用曲面工具建立牙本質(zhì)、牙釉質(zhì)、牙周膜、固有牙槽骨、骨皮質(zhì)和骨松質(zhì)的曲面模型，分別保存為STP格式，并導(dǎo)入至NX三維建模軟件中，使用NX軟件的草圖和實(shí)體建模工具，在牙頸部建立楔狀缺損模型。（2）建立下頜第一前磨牙不同深度與角度楔形缺損的三維模型：在NX軟件中，在已建立的釉牙骨質(zhì)界附近做一橫截面命名為H，使平面H與牙體長(zhǎng)軸垂直，并以該平面作為楔狀缺損尖端所在的平面，從H的中心至該平面上牙齒頰側(cè)的外形高點(diǎn)做一條直線L，并在L上距牙體表面1mm處定點(diǎn)，作為楔狀缺損的頂點(diǎn)，再以L作為楔狀缺損的角平分線，建立了角度為30°、45°和60°的牙齒楔狀缺損實(shí)體模型，并采用相同的方法建立深度為2mm、夾角為30°、45°、60°的牙體缺損實(shí)體模型，并將楔狀缺損的頂點(diǎn)改為極小的弧形，避免此處應(yīng)力過度集中（見圖1）。（3）建立下頜第一前磨牙不同深度與不同角度楔形缺損的有限元模型：將建立好的不同深度和角度的楔狀缺損模型分別導(dǎo)入AnsysWorkbench有限元分析軟件中，使用軟件的前處理工具，對(duì)各個(gè)結(jié)構(gòu)進(jìn)行網(wǎng)格劃分。分別對(duì)各個(gè)結(jié)構(gòu)賦材質(zhì)，所用材料屬性［5-8］見表1，設(shè)定各個(gè)結(jié)構(gòu)之間為bonded接觸關(guān)系，相對(duì)無滑動(dòng)，約束牙槽骨近遠(yuǎn)中面［9］。將所有材料設(shè)置為均質(zhì)、各向同性和線彈性材料。本實(shí)驗(yàn)為了模擬下頜第一前磨牙的正常咬合情況，采用靜態(tài)載荷，加載方向與牙體長(zhǎng)軸呈30°偏頰側(cè)，加載部位取位于頰尖頂偏頰面區(qū)域的一點(diǎn)，加載部位的面積為0.88mm2，載荷力為100N［10］（見圖1），完成下頜第一前磨牙不同深度與角度楔形缺損有限元模型的建立。

圖1 缺損深度角度及缺損尖端設(shè)計(jì)

1.3 材料性能設(shè)定［6，11］見表 1。

1.4 模型分組 根據(jù)缺損深度的不同將模型分為兩組，缺損深度1mm為A組，缺損深度2mm為B組，再根據(jù)缺損角度的不同分為A30°、A45°、A60°、B30°、B45°、B60°共計(jì)6種工況。兩組之間進(jìn)行組間對(duì)照，為臨床上楔狀缺損的充填治療提供依據(jù)。

表1 材料屬性

1.5 主要應(yīng)力分析指標(biāo) 本實(shí)驗(yàn)研究的重點(diǎn)為剩余牙體組織的應(yīng)力。選用Von mises應(yīng)力與最大主應(yīng)力（S1）作為主要應(yīng)力分析指標(biāo)。Von mises應(yīng)力是基于剪切應(yīng)變能的一種等效應(yīng)力值，常用于表示某種材料承受的總體應(yīng)力情況。最大主應(yīng)力（S1）可以反映材料內(nèi)部的某一點(diǎn)在不同方向上的最大拉應(yīng)力［12］。由于牙體組織是脆性材料，抗壓不抗拉，其抗壓強(qiáng)度是抗拉強(qiáng)度的5～6倍，故選用Von mises應(yīng)力與最大主應(yīng)力（S1）為主要應(yīng)力分析指標(biāo)。

2 結(jié)果

2.1 楔狀缺損對(duì)牙體應(yīng)力分布的影響 在100N的側(cè)向載荷條件下，在A、B組中，缺損角度分別為30°、45°、60°時(shí)，剩余牙體的最大主應(yīng)力S1受力分布云圖基本一致：應(yīng)力峰值集中區(qū)域均出現(xiàn)在楔狀缺損的尖端部位，而牙體的冠部與根部受力均較均勻且所受應(yīng)力值也較低（見圖2），6種工況下剩余牙體組織的Von mises應(yīng)力分布云圖見圖3。

圖2 6種工況下剩余牙體組織的最大主應(yīng)力(S1)分布云圖

圖3 6種工況下剩余牙體組織的Von mises應(yīng)力分布云圖

2.2 楔狀缺損深度對(duì)余留牙體受力的影響 楔狀缺損的深度對(duì)余留牙體組織的受力具有一定的影響，表現(xiàn)為：在力的載荷恒定條件下，缺損角度一定時(shí)，缺損深度2mm（B組）組相較缺損深度為1mm（A組），牙體的Von mises應(yīng)力峰值與最大主應(yīng)力S1峰值均有所升高。B30°比A30°的Von mises應(yīng)力峰值升高67.8%，最大主應(yīng)力S1峰值升高70.42%；B45°比A45°的Von mises應(yīng)力峰值升高67.12%，最大主應(yīng)力S1峰值升高59.68%；B60°比A60°的Von mises應(yīng)力峰值升高49.08%，最大主應(yīng)力S1峰值升高48.83%（見表2）。

2.3 楔狀缺損角度對(duì)余留牙體受力的影響 楔狀缺損的角度也同樣影響牙體組織的Von mises應(yīng)力與最大主應(yīng)力S1峰值，具體如下：在缺損深度與力的載荷條件一致時(shí)：隨著楔狀缺損角度的增加，余留牙體組織的Von mises應(yīng)力峰值與最大主應(yīng)力S1峰值呈下降趨勢(shì)。在A組中，A45°與A30°相比，牙體組織的 Von mises應(yīng)力峰值下降了16.18%，最大主應(yīng)力S1峰值下降了15.66%；A60°與A45°相比，牙體組織的 Von mises應(yīng)力峰值下降了2.28%，最大主應(yīng)力S1峰值下降了4.34%。在B組中，B45°與B30°相比，牙體組織的 Von mises應(yīng)力峰值下降了16.52%，最大主應(yīng)力S1峰值下降了20.97%；B60°與B45相比，牙體組織的 Von mises應(yīng)力峰值下降了12.82%，最大主應(yīng)力S1峰值下降了10.84%（見表2）。

表2 6種工況下剩余牙體組織的Von mises與最大主應(yīng)力峰值應(yīng)力峰值

3 討論

楔狀缺損是一種臨床上常見的牙體硬組織非齲性疾病，輕型癥狀可不明顯，當(dāng)疾病進(jìn)一步發(fā)展可出現(xiàn)冷熱及機(jī)械刺激敏感癥狀，較深的楔狀缺損有時(shí)會(huì)累及牙髓，從而導(dǎo)致牙髓炎與根尖周病變，嚴(yán)重時(shí)甚至?xí)?dǎo)致根折［13］，影響美觀和咀嚼功能。第一前磨牙楔狀缺損的發(fā)病率最高、缺損也最為嚴(yán)重，因此，對(duì)楔狀缺損的研究尤為重要。本資料中采用三維有限元分析法（FEM）建立下頜第一前磨牙多類型楔狀缺損的模型，構(gòu)建的模型單元?jiǎng)澐旨?xì)致、圖像清晰完整，建立的模型形態(tài)良好、圖像清晰完整，與天然牙幾何相似度高，模型力的加載和計(jì)算可重復(fù)，能夠較好地模擬下頜第一前磨牙的受力情況。

在本實(shí)驗(yàn)的6組工況中：楔狀缺損的下頜第一磨牙牙體組織的最大主應(yīng)力和Von Mises應(yīng)力峰值均集中于楔狀缺損的尖端。這可能是由于楔狀缺損中斷了應(yīng)力沿著牙體長(zhǎng)軸從牙冠向牙根部的傳遞所導(dǎo)致的［3］。由于牙體硬組織為脆性材料，應(yīng)力的過度集中必然會(huì)對(duì)牙齒結(jié)構(gòu)造成損害。因此，臨床上對(duì)于楔狀缺損的牙齒應(yīng)及時(shí)治療。

本資料結(jié)果顯示，剩余牙體組織的應(yīng)力峰值隨著楔狀缺損深度的增加呈上升趨勢(shì)。楔狀缺損的角度也影響牙體硬組織的應(yīng)力，本實(shí)驗(yàn)數(shù)據(jù)證實(shí)，楔狀缺損的患牙余留牙體組織的應(yīng)力峰值均隨著缺損角度的增加呈下降趨勢(shì)，但本資料中采用的缺損角度僅3種，當(dāng)缺損角度更大時(shí)，余留牙體組織的應(yīng)力分布規(guī)律有何變化還需進(jìn)一步探討。

三維有限元分析結(jié)果表明，楔狀缺損的患牙在楔缺的尖端會(huì)出現(xiàn)應(yīng)力集中現(xiàn)象，且缺損深度越大，這種應(yīng)力集中就越明顯，并且缺損角度越小，牙體組織承受的拉應(yīng)力愈大，提示臨床遇到楔狀缺損的患牙，應(yīng)及早進(jìn)行有效治療，可有效降低患牙頸部的折裂風(fēng)險(xiǎn)。本實(shí)驗(yàn)采用靜力加載，只是在體外一定程度上模擬牙的受力情況，但是口腔內(nèi)的咀嚼運(yùn)動(dòng)是一個(gè)復(fù)雜的過程，而且材料的疲勞問題以及口腔內(nèi)的溫度、口腔內(nèi)腺體分泌等環(huán)境變化也會(huì)對(duì)實(shí)驗(yàn)結(jié)果造成影響。