顏旖萌，歐陽志強，江麗青，劉薇

(1.南昌大學附屬口腔醫(yī)院正畸科，江西 南昌 330000；2.九江學院附屬口腔醫(yī)院正畸科，江西 九江 332000)

有限元分析(finite element analysis，F(xiàn)EA)是正畸生物力學分析中被廣泛應用的一種研究方法。它可以通過建立結構模型，并對模型施加受力條件，從而分析其力學性能，模擬臨床效果。20 世紀 70年代 Thresher等[1]首次將有限元法引入口腔醫(yī)學領域，因其具有非侵入性，實驗數(shù)據(jù)可重復使用，應力分析無可比擬等優(yōu)點在口腔領域各個方向得到廣泛應用。隨著循證醫(yī)學的發(fā)展，口腔正畸學也逐漸從基于經(jīng)驗的實踐主義轉變?yōu)榛谧C據(jù)的實踐主義[2]。任何治療方式都要有科學的理論基礎，但由于正畸加力的獨特方式，細胞實驗或動物實驗并不能完全模擬臨床上實際的正畸牙齒移動和頜骨的生長改良，因此，有限元分析在這方面展現(xiàn)了獨特的優(yōu)勢，它可以根據(jù)數(shù)據(jù)條件分析應力和位移。其中包括3個基本步驟：前處理、處理和后處理。前處理包括幾何模型的構建、模型材料數(shù)據(jù)的轉換、設定接觸關系、構建坐標系、定義邊界條件和設置載荷。處理即為有限元模型中的數(shù)據(jù)分析過程，包括求解線性代數(shù)方程組。后處理包括對結果的解釋及臨床應用[3]。

1 前處理——有限元模型的構建和轉換

切片、磨片法是最早建立牙頜組織數(shù)字模型的方法，但這種建模方法要求使用特殊的實驗設備，人力物力消耗大，效率低，并且屬于破壞性的建模方法，可能會改變原有的生物學結構，對于復雜細微的結構顯示常不清晰，目前已經(jīng)較少采用[4]。隨著數(shù)字成像系統(tǒng)(CT和MRI)的發(fā)展，可以將個體特定的骨骼幾何數(shù)據(jù)輸入到有限元模型中，并且不損壞原物體結構，可根據(jù)需要調節(jié)掃描間距，掃描獲得的數(shù)據(jù)信息可重復應用，每個斷面的解剖結構清晰可辨，能較真實地代表原物體的結構，很好地改善了切片、磨片法的不足之處，而其中近年來出現(xiàn)的錐體束計算機斷層掃描技術(cone beam computer tomography，CBCT)，無論在速度或是精確度上都有明顯的提高，分辨率甚至可以達到微米(μm)級別，以往的測量是在兩條直線之間進行的，而CBCT成像提供了在兩個平面之間進行測量的可能性，從而有助于重建和理解三維模型中的骨骼、牙齒和軟組織皺襞之間的關系，因此逐漸被廣泛應用[5]。在獲得CBCT掃描得到的數(shù)據(jù)后，使用美國放射學會和國家電子制造商協(xié)會聯(lián)合頒布的醫(yī)學數(shù)字成像和通訊的標準——DICOM格式對數(shù)據(jù)進行存儲，有研究表明，以DICOM格式作為存儲格式的圖像處理軟件，除了可以直接讀入CT機輸出的數(shù)據(jù)，同時可以改善CT圖像的質量[6]。常見進行數(shù)據(jù)讀取與處理的軟件有Mimics Medical、Geomagic Studio、Unigraphics NX、Solidworks等。Mimics可以讀入DICOM格式的數(shù)據(jù)，并轉化為初始的三維結構，Geomagic Studio通過對模型進行加工處理、質量提升，生成清晰完整的三維有限元模型，并通過Unigraphics NX以及Solidworks等進行模型的組合與裝配。

2 處理——有限元模型中的數(shù)據(jù)分析

大致分為如下3個步驟[7]：(1)網(wǎng)格劃分，將預求解的模型區(qū)域劃分為有限個單元網(wǎng)格，在三維模型中，單元網(wǎng)格往往為四面體或者多面體，單元間的連接點為節(jié)點，也即為單元的頂點，單元、節(jié)點和節(jié)點連線構成整體的模型網(wǎng)格。(2)單元分析，也就是建立每個單元中節(jié)點位移和節(jié)點力之間的關系，其中節(jié)點位移為基本變量，先通過給單元內部節(jié)點位移確定一個近似表達式，然后建立節(jié)點位移和節(jié)點力之間的關系式。(3)整體分析，對模型區(qū)域整體進行分析，將之前得到每個單元中位移與力之間的關系式，進行數(shù)據(jù)代入，從而揭示整體負載與整個模型之間的關系。較常用于有限元模型數(shù)據(jù)分析的軟件有Ansys、Abaqus等。

3 后處理——有限元分析在口腔正畸生物力學中的應用

3.1 用有限元方法構建正畸矯治器模型

3.1.1 固定正畸矯治器

嚴擁慶等[8]通過有限元方法獲得下頜牙列、托槽和方絲模型，初步建立了用三維有限元法對唇側固定矯治技術進行全牙列分析的基礎。直絲弓矯治器用托槽定位牙齒，簡化了臨床操作，使牙齒定位更加精準迅速，因而得到發(fā)展。相亞寧等[9]建立了顳下頜關節(jié)、上下頜骨及帶有直絲弓矯治器牙列的三維有限元模型，從而有利于研究直絲弓矯治器的力學特性。隨著患者對美學的要求越來越高，越來越多的患者選擇定制的舌側正畸矯治器，因此作為臨床醫(yī)生，有必要了解舌側正畸中有關牙齒運動的生物力學[10]。Hernández-Vázquez等[11]利用三維有限元法對舌側正畸矯治技術中所用的2種不同類型的正畸弓絲做了力學分析，從而發(fā)現(xiàn)蘑菇形正畸弓絲所帶來的應力更大。

3.1.2 功能矯治器

除了傳統(tǒng)的固定正畸矯治器，臨床上還常使用功能矯治器，早期的有限元分析法主要是在二維層面開展，Singh等[12]利用有限元尺度分析法證明了Twin-block 矯治器對軟組織側貌改善的有效性。這之后，林銘等[13]建立了導下頜向前的 Twin-block 矯治器三維有限元模型，從而幫助推進對 Twin-block 矯治機制的研究。Chaudhry等[14]通過三維有限元應力分析法評估固定功能矯治器——Forsus抗疲勞裝置對下頜骨的影響，結果顯示安裝有Forsus抗疲勞裝置的下頜骨模型有更多的應力集中區(qū)，其中下頜皮質骨和髁突區(qū)域的最大主應力及von Mises應力均增加了2倍以上。

3.1.3 隱形矯治器

近年來，隱形矯治技術也逐漸成熟，由于它所具有的美學效果，吸引了越來越多尋求正畸治療的成年患者。起初，隱形矯治器常用于治療輕度或中度前牙擁擠的簡單病例[15-16]，但隨著計算機技術的發(fā)展和人們對矯治器材料生物力學特性的認識加深，隱形矯治技術已顯示出治療更復雜病例的能力，例如拔牙病例。通過三維有限元法可以模擬隱形矯治中牙齒位移模式，從而有利于分析隱形矯治技術的有效性。Jiang等[17]設計了拔除第一前磨牙后采用隱形矯治的三維有限元模型，模擬前后牙的位移模式，結果顯示，前后牙均表現(xiàn)出不同程度的傾斜角度，從而意識到在隱形矯治中牙齒移動的復雜性。曾紅等[18]利用三維有限元法研究了隱形矯治器中附件的不同放置位點對矯治效果的影響。

3.2 用有限元方法研究牙體組織的負載和位移

在有限元方法中，牙體組織常因非均質各向異性的生物材料特性而使得分析過程十分復雜，經(jīng)過研究表明，在進行牙體組織的靜力分析和形變量小時，可以對材料特性進行假設，使其簡化為均勻、連續(xù)、各向同性的線彈性材料[19-21]。利用此方法，李志華等[22]通過有限元方法建立了上頜第一磨牙共4 983個節(jié)點、4 281個單元的三維有限元模型，為上頜第一磨牙在正畸作用下的生物力學分析提供了一個良好的數(shù)字模型。

3.2.1 上頜牙列

阻力中心被視作預測牙齒移動的基本參考點，Luu等[23]對患者CBCT圖像進行處理，以提取與測定相關的三維生物結構，并將其轉化為一個由四面體組成的虛擬網(wǎng)格，然后施加載荷條件，模擬系統(tǒng)上的應力和應變，幫助定位上頜牙齒的阻抗中心，從而為準確預測牙齒移動打下基礎。Kawamura等[24]研究上頜牙列在不同受力角度作用下的牙齒移動情況時發(fā)現(xiàn)，牙齒移動方式與施力的角度有關，當對上頜第二磨牙施加3N的近中力，受力角度與平面成28°時，作用力剛好通過上頜全牙列的阻力中心。

3.2.2 下頜牙列

對于下頜牙列，Chae等[25]通過三維有限元模擬，闡明不同施力角度對下頜牙齒移動模式影響的力學機制，得出結論：力的角度可以影響牙列的運動模式，當力的作用線在阻力中心下方或上方時，整個牙列是順時針或逆時針旋轉的，選擇性地使用不同的施力角度可以實現(xiàn)下頜全牙列的遠移。Lombardo等[26]通過三維有限元建立了下頜中切牙模型，并對該模型施加不同載荷，從而分析不同載荷條件下下頜中切牙受到的生物力學效應，研究顯示，當對下頜中切牙施加舌向壓低力時可以使下頜中切牙接近整體移動，而當施加唇向壓低力時，會使下頜中切牙發(fā)生唇側傾斜，這種改變，對于設計治療計劃至關重要。除此之外，下頜第三磨牙阻生常常影響正畸治療療效的穩(wěn)定性，Oenning等[27]建立了有限元模型來評估阻生的第三磨牙對第二磨牙牙根的生物力學影響。

3.3 用有限元方法模擬頜面部的生物力學行為

3.3.1 顳下頜關節(jié)

顳下頜關節(jié)是一對高度復雜并且可移動的滑膜關節(jié)，同時具有轉動和滑動的雙重運動特點，在日?；顒又?，顳下頜關節(jié)總是通過分散載荷來減少應力峰值，促進下頜運動。同時，顳下頜關節(jié)又是一個承載關節(jié)，在咀嚼過程中需要傳遞力，因此，研究整個口腔功能必須考慮到顳下頜關節(jié)的生物力學[28]。顳下頜關節(jié)紊亂是臨床中一種常見的疾病，雖然該疾病病因尚不明確，但微創(chuàng)傷和關節(jié)內應力被視為是顳下頜關節(jié)紊亂的一些主要原因，以往由于關節(jié)結構復雜，傳統(tǒng)的生物力學研究往往會對關節(jié)造成損傷且實驗不可重復，從而限制了顳下頜關節(jié)生物力學的研究[29]。有限元分析可以用于研究不同的口腔運動、不同載荷條件下顳下頜關節(jié)的生物學行為。Shu等[30]根據(jù)CBCT圖像數(shù)據(jù)重建下頜骨、關節(jié)盤和上頜骨的三維打印模型，然后利用具有高度相似性的三維打印模型對有限元模型進行驗證，發(fā)現(xiàn)用有限元模型來模擬顳下頜關節(jié)和上下牙列之間的相互作用是準確的。Lai等[31]通過有限元模型發(fā)現(xiàn)在顳下頜關節(jié)紊亂患者中，髁突、關節(jié)盤和下頜骨的應力值均高于正常側關節(jié)。

3.3.2 上下頜骨

4 發(fā)展前景

正畸治療主要是通過各種不同類型的矯正裝置來協(xié)調口腔頜面部骨骼、牙齒、神經(jīng)肌肉之間的關系以達到口頜系統(tǒng)的平衡、穩(wěn)定和美觀。因此，在研究牙頜組織受力后的應力分布和移動趨勢方面，三維有限元分析有其獨特的優(yōu)勢，在正畸生物力學研究中有良好的應用前景。然而，三維有限元分析也有一定的局限性。

4.1 需要進行臨床驗證

由計算機分析得出的力學數(shù)值目前只能為臨床操作及結果預測提供一定的參考，并不能完全模擬口腔內的實際臨床狀況，其實際應用效果需要更多的臨床驗證。

4.2 優(yōu)化建模的材料屬性

由于牙頜組織非均質各向異性、粘彈性等的生物材料特性復雜，使得目前三維有限元實現(xiàn)精確建模難度較大，分析時往往使其簡化為均勻、連續(xù)、各向同性的線彈性材料，這也就意味著一定程度上的信息丟失，從而有可能影響到整個模型的幾何及力學特性。

4.3 實現(xiàn)動態(tài)建模

現(xiàn)有的三維有限元模型大多還是靜態(tài)模型，而人類的生理活動本身是動態(tài)的，如矯治過程中發(fā)生的牙槽骨改建等。并且隨著矯治技術的提高，靜磁場和生物學電傳導等作用會隨著時間的推移而發(fā)生力值的改變。在未來，三維有限元建模應該要致力實現(xiàn)從靜態(tài)模型到動態(tài)模型的轉化，以實現(xiàn)真正的生物仿真。

可以預見，通過更多的醫(yī)學研究探索，在不久的將來，三維有限元分析一定可以逐步克服自身的瓶頸，實現(xiàn)生物力學的精確模擬，從而推動整個正畸事業(yè)的發(fā)展。