考慮材料參數(shù)不確定性的皮膚傷口縫合力預(yù)測模型*

溫廣全, 紀(jì)小剛,2, 段玉順, 鄧 霖

(1. 江南大學(xué) 機(jī)械工程學(xué)院, 江蘇 無錫 214122;2. 江蘇省食品先進(jìn)制造裝備技術(shù)重點(diǎn)實(shí)驗(yàn)室, 江蘇 無錫 214122)

0 引 言

皮膚傷口縫合作為臨床醫(yī)生的必備技能,是外科手術(shù)治療的重要環(huán)節(jié)．隨著科學(xué)技術(shù)的進(jìn)步,機(jī)器人手術(shù)縫合技術(shù)得到快速發(fā)展[1],拓展了外科醫(yī)生的診斷治療能力．對(duì)于機(jī)器人手術(shù)縫合過程,力覺信息對(duì)其縫合效果的影響至關(guān)重要[2],了解閉合傷口所需縫合力對(duì)機(jī)器人縫合算法的開發(fā)具有重要的指導(dǎo)作用[3]．

近年來,數(shù)值仿真技術(shù)已成為研究皮膚等生物軟組織力學(xué)響應(yīng)的強(qiáng)大工具[4],可有效地用于預(yù)測閉合傷口所需的力[5]．然而,皮膚等生物組織材料具有一定的特殊性及復(fù)雜性[6],材料模型往往難以與皮膚真實(shí)的力學(xué)響應(yīng)完全匹配[7]．盡管現(xiàn)在各國學(xué)者提出了一些更為復(fù)雜,針對(duì)性更強(qiáng)的材料模型[8-9],大大提高了模型的準(zhǔn)確性和預(yù)測能力,但材料參數(shù)的確定依賴實(shí)驗(yàn)數(shù)據(jù)的準(zhǔn)確性與多樣性,且材料特性受年齡、性別、身體部位等因素影響具有較大的差異性[10],這就造成數(shù)值計(jì)算過程存在輸入材料參數(shù)的不確定性,如不加以考慮,可能會(huì)導(dǎo)致計(jì)算結(jié)果產(chǎn)生波動(dòng)甚至偏離預(yù)期[11]．為了嚴(yán)格評(píng)估材料參數(shù)的不確定性對(duì)傷口縫合力預(yù)測結(jié)果的影響,亟需一套不確定性傳播分析方法．

決策和優(yōu)化過程中,不確定性傳播分析在減少不確定性的影響方面起著關(guān)鍵作用,已被應(yīng)用于解決科學(xué)和工程領(lǐng)域的各種現(xiàn)實(shí)問題[12]．然而,很少有針對(duì)皮膚組織數(shù)值計(jì)算不確定性問題的相關(guān)研究．為了得到可靠的傷口縫合力預(yù)測結(jié)果,如何在充分了解皮膚力學(xué)特性的基礎(chǔ)上,量化影響皮膚傷口縫合數(shù)值計(jì)算結(jié)果的不確定性材料參數(shù),進(jìn)一步獲取材料參數(shù)不確定性影響下皮膚傷口縫合力的概率統(tǒng)計(jì)信息,是亟需解決的問題．Monte-Carlo(MC)方法[13]結(jié)合代理模型技術(shù)[14]為此問題的解決提供了有效的思路．MC方法因構(gòu)造簡單,在不確定性傳播分析中應(yīng)用廣泛[15]．而代理模型通過構(gòu)建一個(gè)數(shù)學(xué)表達(dá)式來描述輸入變量和輸出響應(yīng)之間的關(guān)系,只需要運(yùn)行有限數(shù)量的仿真計(jì)算,便可得到寬泛參數(shù)空間下的預(yù)測結(jié)果,用于解決MC方法對(duì)計(jì)算資源需求大的問題[16]．同時(shí)通過使用光滑連續(xù)的插值函數(shù),代理模型可以有效減少數(shù)值計(jì)算過程產(chǎn)生的噪聲．常用的代理模型方法有響應(yīng)面模型、Gauss過程回歸和人工神經(jīng)網(wǎng)絡(luò)[17]等．其中,人工神經(jīng)網(wǎng)絡(luò)方法具有近似復(fù)雜非線性函數(shù)的能力、極好的泛化能力,在預(yù)測材料的力學(xué)響應(yīng)方面具有一定的優(yōu)勢(shì)[18]．

本文基于皮膚組織結(jié)構(gòu)及力學(xué)特性,通過文獻(xiàn)歸納及相關(guān)皮膚力學(xué)試驗(yàn),獲取了大量構(gòu)建不確定性皮膚傷口縫合力預(yù)測模型所必需的力學(xué)參數(shù)信息．然后將橢球基(elliptical basis functions, EBF)神經(jīng)網(wǎng)絡(luò)代理模型方法與有限元方法、MC方法相結(jié)合,構(gòu)建出皮膚傷口縫合力預(yù)測模型．最后分析了皮膚材料參數(shù)的不確定性對(duì)計(jì)算結(jié)果的影響,并通過實(shí)驗(yàn)對(duì)預(yù)測結(jié)果進(jìn)行驗(yàn)證,以期為機(jī)器人手術(shù)縫合提供可靠的力學(xué)參考數(shù)據(jù)．

1 皮膚傷口縫合力數(shù)值計(jì)算

1.1 皮膚材料本構(gòu)模型

皮膚結(jié)構(gòu)從上到下分為表皮層、真皮層和皮下組織,其中真皮層是主要承載層,由嵌在基質(zhì)中的纖維網(wǎng)絡(luò)組成．在力學(xué)上,皮膚表現(xiàn)出非線性的應(yīng)力-應(yīng)變響應(yīng)特性,且在很小的外力下即可產(chǎn)生較大的變形,Holzapfel超彈性本構(gòu)模型可以很好地描述其力學(xué)行為[8-9]．其應(yīng)變勢(shì)能函數(shù)為

(1)

(2)

應(yīng)變能函數(shù)中C10,k1,k2,D,κ是有限元分析的輸入?yún)?shù)．其中,C10與材料的剪切模量相關(guān),反應(yīng)皮膚非膠原纖維部分的組織剛度;k1代表皮膚纖維的剛度模量;k2為無量綱參數(shù),反應(yīng)組織的非線性力學(xué)性能;D與材料的體積模量相關(guān),參數(shù)κ(0<κ<1/3)描述了皮膚纖維方向的分散程度,κ= 0纖維完全對(duì)齊,κ= 1/3纖維隨機(jī)分布,材料呈各向同性．

1.2 皮膚傷口有限元模型

使用ABAQUS有限元分析軟件對(duì)皮膚傷口進(jìn)行建模計(jì)算．傷口輪廓如圖1所示．皮膚組織被簡化為100 mm×100 mm的平面,傷口形狀被定義為紡錘形,尺寸由長短軸確定．長軸a固定為40 mm,短軸b分別取5 mm,10 mm,15 mm,20 mm,模型厚度t參考人體皮膚厚度范圍[4]分別取1.5 mm,2 mm,2.5 mm,3 mm,3.5 mm,4 mm．縫合間距取5 mm,傷口共計(jì)需縫合7次,縫合點(diǎn)位置如圖1所示．

圖1 皮膚傷口輪廓 圖2 皮膚傷口有限元簡化模型 Fig. 1 The skin wound profileFig. 2 The simplified finite element model for the skin wound

考慮到模型的對(duì)稱性, 取傷口的二分之一進(jìn)行分析．短軸b與厚度t排列組合, 24種傷口尺寸模型被創(chuàng)建．皮膚組織采用Holzapfel超彈性本構(gòu)模型,C10,k1,k2,D,κ是需要輸入的材料參數(shù)．鑒于皮膚軟組織通常被認(rèn)為是近似不可壓縮的[19],同時(shí)為了保證所有輸入?yún)?shù)的收斂,D被固定為0.2．在分析中,假設(shè)κ= 1/3,皮膚表現(xiàn)為各向同性．雖然皮膚的確表現(xiàn)出各向異性,但纖維色散相對(duì)較大[20],隨機(jī)纖維色散是一種合理的近似．以往關(guān)于皮膚力學(xué)特性的研究[10, 20],為人體皮膚本構(gòu)模型參數(shù)C10,k1,k2提供了合理的取值范圍,歸納如表1所示．此外,由于活體皮膚正常情況下處于張力狀態(tài)[4],為得到皮膚組織中力學(xué)響應(yīng)更真實(shí)的預(yù)測,分別在平面模型X,Y方向橫截面上施加預(yù)應(yīng)力σX,σY,取值參考文獻(xiàn)[5],歸納如表1所示．為了模擬手術(shù)縫合過程,減少建模過程對(duì)分析結(jié)果的影響,模型統(tǒng)一采用八節(jié)點(diǎn)線性六面體單元,根據(jù)縫合間距統(tǒng)一單元尺寸為5 mm．通過設(shè)置7個(gè)靜力通用分析步,在傷口邊緣節(jié)點(diǎn)處逐次施加線性位移約束,以模擬傷口的閉合．有限元簡化模型如圖2所示．有限元計(jì)算結(jié)果為7個(gè)分析步中,施加位移約束節(jié)點(diǎn)的反力,記為F1～F7,以此作為傷口縫合過程中縫合線對(duì)傷口邊緣的縫合力．

表1 人體皮膚材料參數(shù)取值范圍

2 皮膚傷口縫合力預(yù)測模型構(gòu)建

2.1 樣本數(shù)據(jù)的獲取

以尺寸及材料參數(shù)的不同組合得到的皮膚傷口縫合力數(shù)值計(jì)算結(jié)果,作為構(gòu)建皮膚傷口縫合力預(yù)測模型的樣本數(shù)據(jù)．考慮表1獲取的材料參數(shù)可能存在誤差,本文將表1的參數(shù)范圍進(jìn)一步擴(kuò)展20%,以便在較大參數(shù)范圍內(nèi)訓(xùn)練代理模型后,用來傳播更小范圍的、真實(shí)的參數(shù)聯(lián)合概率分布信息．

為充分探索整個(gè)材料參數(shù)范圍對(duì)輸出結(jié)果產(chǎn)生的響應(yīng),采用優(yōu)化拉丁超立方采樣方法,創(chuàng)建了基于參數(shù)b,t,σX,σY,C10,k1和k2的2 000組不同組合,總共進(jìn)行了2 000次仿真計(jì)算．本文首先隨機(jī)抽取1 400組仿真數(shù)據(jù)用來構(gòu)建代理模型,剩余600組數(shù)據(jù)用作驗(yàn)證數(shù)據(jù)集,以校驗(yàn)代理模型的擬合精度．

2.2 EBF神經(jīng)網(wǎng)絡(luò)模型的構(gòu)建

EBF神經(jīng)網(wǎng)絡(luò)模型采用橢圓單元隱層和線性單元輸出層,來描述輸入變量和輸出響應(yīng)之間的關(guān)系．本文所構(gòu)建的EBF神經(jīng)網(wǎng)絡(luò)模型基本結(jié)構(gòu)如圖3所示．其中輸入層的x1～x7分別代表傷口尺寸b,t和材料參數(shù)σX,σY,C10,k1,k2,輸出層的y1～y7分別代表有限元輸出結(jié)果F1～F7．

圖3 EBF神經(jīng)網(wǎng)絡(luò)基本結(jié)構(gòu)Fig. 3 The basic structure of the EBF neural network

設(shè)輸入層中變量個(gè)數(shù)為N,隱含層樣本數(shù)為n,輸出層響應(yīng)變量個(gè)數(shù)為M．以待測點(diǎn)與樣本點(diǎn)之間的Mahalanobis距離為自變量,通過線性疊加構(gòu)造出EBF神經(jīng)網(wǎng)絡(luò)模型,其響應(yīng)函數(shù)表達(dá)式如下:

(3)

式中,gi(x)為基函數(shù),由式(4)給出

gi(x)≡g(‖x-xi‖m),i=1,2,…,N,

(4)

‖x-xi‖m=(x-xi)TS-1(x-xi),

(5)

(6)

式中,‖x-xi‖m是待測點(diǎn)與樣本點(diǎn)之間的Mahalanobis距離,S為多維隨機(jī)變量的協(xié)方差矩陣,μ為樣本中心點(diǎn)．

式(3)中αi為權(quán)重系數(shù),通過求解式(7)線性方程組求得

(7)

2.3 不確定性傳播分析

構(gòu)建好代理模型后,除了無需在特定輸入上重新進(jìn)行有限元計(jì)算便可得到預(yù)測結(jié)果外,還可以進(jìn)行不確定性傳播分析．本文將采用MC方法進(jìn)行皮膚傷口縫合力預(yù)測中的材料參數(shù)不確定性問題分析．

MC方法也稱統(tǒng)計(jì)試驗(yàn)方法, 是以概率統(tǒng)計(jì)理論為基礎(chǔ)的一種數(shù)值計(jì)算方法．其基本思想是, 當(dāng)想要獲取某個(gè)隨機(jī)變量的期望值時(shí), 可以采用抽樣試驗(yàn)或隨機(jī)模擬的方法, 得到這個(gè)隨機(jī)變量的平均值, 以此作為問題的解．

本文實(shí)現(xiàn)步驟為:

1) 假設(shè)皮膚材料參數(shù)圍繞實(shí)驗(yàn)測量或經(jīng)驗(yàn)評(píng)估結(jié)果呈正態(tài)分布,根據(jù)獲取材料參數(shù)的可靠程度,給出參數(shù)的不確定度．本文統(tǒng)一設(shè)為±10%．

2) 采用描述抽樣方法,抽取10 000個(gè)樣本點(diǎn),通過代理模型快速得到縫合力輸出響應(yīng)結(jié)果．需要強(qiáng)調(diào)的是,MC方法統(tǒng)計(jì)結(jié)果的準(zhǔn)確性非常依賴于樣本數(shù)據(jù)的規(guī)模,規(guī)模越大樣本統(tǒng)計(jì)特性越優(yōu)越,因此,在計(jì)算資源足夠的情況下,應(yīng)抽取盡可能多的樣本數(shù)據(jù)．

(8)

(9)

其中,Q為抽樣容量,yi(i=1,2,…,Q)為樣本輸出變量的響應(yīng)值．

4) 計(jì)算輸出結(jié)果總體均值的置信區(qū)間:

(10)

式中,zα/2為標(biāo)準(zhǔn)分?jǐn)?shù),置信水平為100(1-α),α通常取值0.1或0.05．

氣象信息化建設(shè)是一項(xiàng)持續(xù)性的工作,對(duì)于基層氣象部門來說,信息化基礎(chǔ)設(shè)施構(gòu)建要按照上級(jí)氣象部門的統(tǒng)一部署,符合中國氣象局《氣象信息化發(fā)展規(guī)劃(2018—2022年)》構(gòu)建統(tǒng)籌集約、協(xié)同高效、開放共享、安全可靠的氣象信息化體系的總體目標(biāo),在此基礎(chǔ)上廣泛應(yīng)用成熟、先進(jìn)的信息化技術(shù)、同時(shí)參考本地區(qū)氣象部門經(jīng)過實(shí)踐檢驗(yàn),廣泛應(yīng)用的技術(shù)規(guī)范,以便和本地區(qū)氣象部門順利對(duì)接,逐步探索符合當(dāng)?shù)貧庀髽I(yè)務(wù)特點(diǎn)的氣象信息化建設(shè)思路和建設(shè)模式,以期為提高本地區(qū)氣象災(zāi)害監(jiān)測防御能力,預(yù)防和減輕氣象災(zāi)害損失,保障人民生命財(cái)產(chǎn)安全,促進(jìn)當(dāng)?shù)亟?jīng)濟(jì)社會(huì)健康發(fā)展提供強(qiáng)有力的保障。

3 案 例 分 析

文獻(xiàn)表明,豬皮膚生理結(jié)構(gòu)和力學(xué)特性與人體相似[21]．本文以普通家豬腹部皮膚組織為研究對(duì)象,通過實(shí)驗(yàn)測量縫合傷口所需的力,并與預(yù)測模型進(jìn)行比較,以驗(yàn)證模型的可靠性．實(shí)驗(yàn)均在同一環(huán)境下進(jìn)行,未考慮溫度、濕度等外在因素的影響．

3.1 皮膚材料參數(shù)獲取

由于缺乏家豬皮膚的材料參數(shù)信息,本文設(shè)計(jì)了兩個(gè)實(shí)驗(yàn)來進(jìn)行評(píng)估．需要獲取的材料參數(shù)有:皮膚本構(gòu)模型參數(shù)和皮膚預(yù)應(yīng)力參數(shù)．

圖4給出了通過試驗(yàn)最終得到的真實(shí)應(yīng)力-應(yīng)變曲線．為獲取皮膚本構(gòu)模型參數(shù),我們制作了圖4所示的厚度為3 mm的板狀皮膚拉伸試樣,試樣取自五個(gè)月大、重約120 kg的普通家豬腹部．使用CTM2500萬能材料試驗(yàn)機(jī)(最大荷重為10 kN,荷重精度為±0.01%,精度等級(jí)為0.5級(jí),位移分辨率為0.03 μm)對(duì)試樣進(jìn)行單軸拉伸試驗(yàn)．試驗(yàn)采用楔形夾具,夾具與皮膚組織接觸處有許多整齊的小齒,可以有效增大夾爪與皮膚組織之間的摩擦因數(shù),避免試樣在拉伸的過程中發(fā)生滑移．應(yīng)變測量采用標(biāo)距為25 mm,量程為25 mm的引伸計(jì)．根據(jù)大變形后體積V不變的假設(shè),即A0L0=AL,真實(shí)應(yīng)力σ=F/A=FL/(A0L0);真實(shí)應(yīng)變?chǔ)?ln(L/L0)．其中,F為拉伸載荷,A0,A分別為試樣初始及任一瞬時(shí)的橫截面積,L0,L(L=L0+ΔL)分別為初始標(biāo)距和任一瞬時(shí)的標(biāo)距．由真實(shí)應(yīng)力-應(yīng)變曲線擬合得到皮膚的本構(gòu)模型參數(shù)為C10=24.22 kPa,k1=24 421 kPa,k2=369.7．

圖4 豬皮單軸拉伸試樣尺寸、試驗(yàn)設(shè)備及應(yīng)力-應(yīng)變曲線Fig. 4 Sizes, the test equipment and the stress-strain curve of the pig skin uniaxial tensile test specimen

為獲取皮膚預(yù)應(yīng)力信息,參考金屬構(gòu)件殘余應(yīng)力測量方法[22],從普通家豬腹部切取100 mm×100 mm皮膚組織試樣進(jìn)行實(shí)驗(yàn)．由于預(yù)應(yīng)力釋放,切取的皮膚組織輪廓產(chǎn)生尺寸收縮．假設(shè)尺寸收縮是由皮膚預(yù)應(yīng)力彈性釋放造成的,并且皮膚橫截面在預(yù)應(yīng)力釋放過程中產(chǎn)生等變形收縮．如果施加外力將變形后的皮膚組織恢復(fù)到切割前的狀態(tài),所得到的截面應(yīng)力就等效于切割前該截面的預(yù)應(yīng)力．基于此,采用彈性有限元方法,計(jì)算皮膚輪廓恢復(fù)到切割前的應(yīng)力狀態(tài)．具體步驟如下:

1) 在ABAQUS中以收縮后的試樣尺寸建模,鑒于模型的對(duì)稱性,取皮膚組織四分之一進(jìn)行建模,并采用八節(jié)點(diǎn)線性六面體單元?jiǎng)澐志W(wǎng)格;

2) 本構(gòu)模型參數(shù)采用上述單軸拉伸試驗(yàn)結(jié)果;

3) 模型邊界施加對(duì)稱約束并限制底面Z向位移,切割橫截面上施加位移約束使試樣尺寸變?yōu)槌跏汲叽?00 mm×100 mm,如圖5(a)所示;

3.2 傷口縫合力測量

在普通家豬腹部分別設(shè)計(jì)了40 mm×10 mm和40 mm×14 mm、厚度為3 mm的各3組皮膚傷口;使用精度為0.01 N,型號(hào)為AIPLI-SF-10的電子拉力計(jì)測量閉合傷口所需的力;采用單純間斷縫合法閉合傷口,縫合線兩端水平對(duì)齊,同時(shí)拉動(dòng)兩端縫合線,使傷口邊緣沿中軸線緊密對(duì)齊;重復(fù)測量,記錄均值與標(biāo)準(zhǔn)差,傷口尺寸及縫合力測量方法如圖6所示,兩種傷口縫合力測量結(jié)果如圖7所示．其中橫坐標(biāo)1～7代表傷口的7個(gè)縫合點(diǎn),其位置如圖1所示,縱坐標(biāo)為該點(diǎn)縫合力測量結(jié)果．

(a) 邊界條件 (b) 橫向預(yù)應(yīng)力 (c) 縱向預(yù)應(yīng)力 (a) Boundary conditions(b) Transverse prestresses (c) Longitudinal prestresses圖5 皮膚組織預(yù)應(yīng)力釋放的反向模擬邊界條件及計(jì)算結(jié)果Fig. 5 Boundary conditions and calculation results of the reverse simulation of prestress release in the skin tissue

圖6 傷口尺寸及縫合力測量方法Fig. 6 Measurement of wound sizes and suture forces

圖7 傷口縫合力測量結(jié)果Fig. 7 Wound suture force measurement results

3.3 模型預(yù)測結(jié)果統(tǒng)計(jì)分析

材料參數(shù)取實(shí)驗(yàn)估計(jì)值時(shí),40 mm×10 mm傷口縫合力有限元計(jì)算結(jié)果如圖8所示．其中圖8(a)為傷口施加預(yù)應(yīng)力后的初始狀態(tài),圖8(b)—(h)為縫合過程．可以看出施加位移約束節(jié)點(diǎn)處,即縫合點(diǎn)處力最大,并且7個(gè)縫合點(diǎn)處縫合力大小按縫合針次依次呈先增高后逐漸降低趨勢(shì),縫合力大小在0～1.5 N之間,與圖7實(shí)驗(yàn)測量結(jié)果接近．

圖8 有限元模擬縫合過程Fig. 8 The finite element simulation of the suture process

構(gòu)建的EBF神經(jīng)網(wǎng)絡(luò)模型需要進(jìn)行校驗(yàn),本文采用式(11)的決定系數(shù)R2來評(píng)估所構(gòu)建代理模型的擬合精度:

(11)

圖9 3種代理模型擬合精度 圖10 輸出響應(yīng)概率分布圖 Fig. 9 The fitting precision of 3 proxy modelsFig. 10 The output the response probability distribution graph

輸入測得的豬皮材料參數(shù)及傷口尺寸參數(shù),通過構(gòu)建的皮膚傷口縫合力預(yù)測模型,計(jì)算得到各輸出變量的概率密度分布、均值以及標(biāo)準(zhǔn)差．以40 mm×10 mm傷口中心節(jié)點(diǎn)處縫合力F4的統(tǒng)計(jì)結(jié)果為例,得到輸出響應(yīng)的概率分布,如圖10所示．可以看出,結(jié)果呈現(xiàn)出良好的正態(tài)分布,表明參數(shù)樣本具有良好的統(tǒng)計(jì)特性,同時(shí)受材料參數(shù)不確定性影響,傷口縫合力預(yù)測結(jié)果表現(xiàn)出較高的離散性．

模型預(yù)測結(jié)果與實(shí)驗(yàn)測量結(jié)果的對(duì)比如圖11所示．結(jié)果顯示,采用間斷縫合橢圓形皮膚傷口,所需縫合力呈先增后減趨勢(shì),峰值力發(fā)生在傷口中線前．40 mm×10 mm傷口縫合力峰值約為1.7 N;40 mm×14 mm傷口縫合力峰值約為2.5 N．受材料參數(shù)不確定性影響,縫合力預(yù)測結(jié)果最大有±0.6 N的波動(dòng),隨著所取置信水平的增大,模型預(yù)測結(jié)果的置信區(qū)間變大,絕大部分實(shí)驗(yàn)測量結(jié)果落在95%置信區(qū)間范圍,少數(shù)結(jié)果超出這一范圍,所有測量結(jié)果均落在90%置信區(qū)間范圍,相比于單純地進(jìn)行一次仿真計(jì)算,考慮了參數(shù)的不確定性后得到的結(jié)果更加穩(wěn)健可靠．從40 mm×14 mm傷口對(duì)比結(jié)果可以看出,對(duì)于建模尺寸外的傷口,預(yù)測模型亦有較好的預(yù)測效果．決策者可以根據(jù)不同的精度要求,選擇相信不同置信水平對(duì)應(yīng)的縫合力置信區(qū)間,作為機(jī)器人縫合算法的參考輸入．

(a) 40 mm×10 mm傷口 (b) 40 mm×14 mm傷口 (a) 40 mm×10 mm wound (b) 40 mm×14 mm wound圖11 傷口縫合力預(yù)測結(jié)果與實(shí)驗(yàn)測量結(jié)果對(duì)比Fig. 11 The prediction results of wound suture forces compared with the experimental results

4 結(jié) 論

本文應(yīng)用有限元理論對(duì)皮膚傷口縫合過程進(jìn)行力學(xué)仿真計(jì)算,并將代理模型方法及MC方法結(jié)合,得到了考慮皮膚材料參數(shù)不確定性的皮膚傷口縫合力預(yù)測模型．通過與普通家豬皮膚傷口縫合實(shí)驗(yàn)結(jié)果相比較,驗(yàn)證了模型的可靠性．

此外,皮膚材料實(shí)際上還表現(xiàn)出黏彈性、可壓縮性、各向異性等材料力學(xué)特性,本文方法為分析皮膚等生物軟組織數(shù)值仿真不確定傳播問題提供了一般思路,額外考慮皮膚的黏彈性、可壓縮性、各向異性,并與皮膚力學(xué)特性在體測試技術(shù)相結(jié)合,將得到更加準(zhǔn)確的預(yù)測結(jié)果．同時(shí)該方法還可以擴(kuò)展到更廣泛的傷口尺寸及形狀,期望能為機(jī)器人手術(shù)縫合算法的開發(fā)及術(shù)前指導(dǎo)提供參考依據(jù)．