由於計算機水平和軟體能力的迅速提高,數值模擬技術取得很大的進展。人們已經開發出頭部、大腦和血管等模型來研究頭部損傷機理和損傷防護。
有限元法是研究人體組織損傷機理的重要方法,由於人體組織結構的複雜性,如何建立高生物模擬度的有限元模型是有限元計算中必須解決的首要問題;其次,如何解決數值模擬中的高度非線性、流固耦合問題,是有限元模擬的關鍵。
Mimics與ANSYS聯合模擬技術
Mimics軟體是一套整合逆向工程、材料工程、生物力學工程、有限元分析等 多工業、多學科領域的統一解決方案的專業軟體,它透過對三維體圖或二維平面序列影象(CT/MRI)進行提取、分割、合併等操作,生成表面重建後的三維影象。
ANSYS是國際上最著名的通用有限元分析軟體之一,在計算理論和求解問題的廣泛性方面處於全球領先的地位,尤其針對結構非線性、流/固耦合等複雜問題的求解具有強大優勢,被業內人士認為是非線性有限元發展方向的先導,被廣泛應用於各個工業領域的工程模擬計算,包括土木建築、交通運輸、石油化工、機械製造、航空航天、汽車、國防軍工、船舶、生物力學以及科學研究等各個領域。
運用Mimics與ANSYS聯合模擬技術,可以將人體組織的CT/MRI二維斷層掃描圖片資料建立三維模型。ANSYS強大的非線性力學計算功能和流固耦合求解能力,高效並精確的模擬組織應力、應變、傳熱、流體動力學等問題。
分析流程
腦脊髓的流體動力學計算
基於MRI影象,運用逆向工程軟體進行三維模型重建,並建立了位於中樞神經系統中的腦脊髓三維流體動力學有限元模型。本例中需要計算腦脊髓的流動和迴圈系統中腦組織的變形,因此計算中需要呼叫求解器(流固耦合)
從MRI影象到三維有限元模型
模擬的腦脊髓流動動力學模擬:壓力、速度分佈
腦積水病人腦部應力分佈和速度場分佈
腦動脈瘤的流固耦合計算
腦動脈瘤是動脈病態的膨脹,一般發現於韋利斯氏環的前部和後部。腦動脈瘤的破裂會引起蛛網膜下出血以及嚴重的併發症。因此,在腦動脈瘤破裂的分析中,血液動力學發揮著很重要的作用。下面的例子模擬了血液的流動以及流固耦合問題。模型的建立運用了三維影象處理技術,並選用有限元分析軟體中的非牛頓不可壓縮流體、流固耦合求解器。
影象處理並輸出有限元網格
圖示分別為血液壓力分佈、組織剪下應力、有效應力、位移分佈
心臟的流固耦合模擬
大多數的生物組織運動過程都與流固耦合有關。本例就使用了流固耦合計算方法,模擬了某位病人左右心房和膜片的模型,目的是用來最佳化心臟肺動脈瓣的外科手術。
心臟的流固耦合模型
心臟的主應力和主應變雲圖
主動脈瓣的流固耦合模擬
本例使用流體分析軟體建立了主動脈根和主動脈瓣的二維平面應變模型,弱可壓縮流體,並在入口給定14mmHg的壓力。
不同時刻的速度場分佈
多孔彈性有限元模型預測腰椎間盤的漸進破壞
腰椎間盤突出症是較為常見的疾患之一,主要是因為腰椎間盤各部分(髓核、纖維環及軟骨板),尤其是髓核,有不同程度的退行性改變後,在外力因素的作用下,椎間盤的纖維環破裂,髓核組織從破裂之處突出(或脫出)於後方或椎管內,導致相鄰脊神經根遭受刺激或壓迫,從而產生腰部疼痛,一側下肢或雙下肢麻木、疼痛等一系列臨床症狀。腰椎間盤複雜的模型建立,是有限元分析的重點。
本例中使用有限元分析軟體進行迴圈載入,來模擬腰椎間盤漸進破壞。模擬的結果可以用來研究如何避免由於重複的負重而引起的疼痛。
使用多孔介質材料來模擬椎間盤,並且考慮固相和液相。下圖為腰椎間盤的有限元模型。
腰椎間盤網格圖
一個迴圈載入後,各部位所受的力與時間的關係如下圖所示。
受力與時間的關係
同時,我們也能獲得在不同載入力和迴圈次數條件下的椎間盤有效應力雲圖。
有效應力雲圖