螺旋機器人的結構參數對人體腸道壁損傷影響的優化研究

梁 亮 唐 勇 龐佑霞 胡冠昱

長沙學院,長沙,410003

0 引言

微創外科醫療技術是當今國際上的一個研究熱點。目前,主要從兩方面對微創外科醫療手術進行研究:一是改進傳統的醫用內窺鏡系統,將系統進一步微型化;二是改進醫用內窺鏡系統進入人體內腔的驅動方式。例如,研制與開發能自動把醫用內窺鏡系統送到人體內腔各個部位的微型機器人[1-5]。但這些機構作為醫用機器人在進入人體內腔后都直接與內腔壁接觸,或多或少會給患者帶來不適和痛苦。另外,以色列開發的胃腸道微創內窺鏡膠囊不具備驅動結構[6],也只能隨消化道蠕動做被動運行。

現有的內窺鏡機器人均或多或少地存在一些問題和缺陷?；谏鲜龇治?文獻[7-8]介紹了螺旋式無損傷醫用微型機器人,并對機器人運行速度和懸浮的高度等進行了研究;張永順等[9]設計了螺旋式腸道膠囊微型機器人,并對驅動原理和螺旋參數等進行了研究。利用螺旋外殼,以達到無損傷驅動目的時,螺旋外殼高速旋轉運行所產生的動壓力,人體是否能夠承受,是否會損傷腸道,什么樣的機器人結構參數能夠減小動壓力,都需要進一步分析研究和優化。

筆者設計了一種結構簡單的單節螺旋式無線無創微型機器人。建立了腸道黏液的動力學控制方程,并借助計算流體力學(CFD)軟件計算了機器人高速旋轉運行時對人體腸道壁產生的壓力,比較了不同機器人結構參數對腸道壁所受壓力的影響,優化了結構參數。

1 螺旋式微型機器人的驅動機構

圖1為螺旋式微型機器人內部驅動機構示意圖。該機器人是一個帶螺旋槽外殼的圓柱體,內部裝有微電池、微電機、4根形狀記憶合金SMA構成的轉向機構、圖像傳輸模塊、無線收發與控制模塊。微電機外表面與外殼內表面之間為緊配合,微電機外表面與外殼內表面之間無摩擦;當微電機轉軸旋轉時,將帶動帶螺旋槽的圓柱體外殼旋轉。當正向接通微電機電源時,帶螺旋槽外殼產生的軸向摩擦牽引力帶動微型機器人前進;當反向接通微電機電源時,微型機器人后退。連接形狀記憶合金(SMA)的兩端可相對移動,通過控制4根SMA的伸縮狀態和幅度實現任意轉向。人體腸道中存在的黏液可看成是一層動壓潤滑黏液膜,它使微型機器人處于懸浮狀態,該潤滑膜能避免機器人與腸道壁發生直接接觸。這樣就達到了無創或微創驅動的目的。同時機器人內部帶有電池和無線收發及控制裝置,又可實現對機器人的無線控制。

2 腸道黏液的動力學控制方程

要研究機器人旋轉運行對人體腸道壁產生的損傷,也就是要研究機器人和腸道之間黏液對腸道壁的壓力,因此首先要建立腸道黏液的動力學控制方程。

雖然人體腸道內的黏液在低切變速率下多為黏彈性流體,但在中高切變速率下則可把它看作為牛頓流體[10],并且假設黏液的密度ρ為常數,以黏液中的某一微元體為研究對象,可得質量守恒方程:

根據動量守恒,可得方程:

式中,u、v、w分別為黏液微元體速度矢量在x、y、z方向的分量;μ為黏液動力黏度;p為黏液微元體上的壓力;Fx、Fy、Fz為黏液微元體上的體力。

因為體力只有重力,加上重力方向與y軸相反,則Fx=Fz=0,Fy=—ρ g 。式(1)和式(2)即為腸道黏液的動力學控制方程,它是腸道黏液流場數值計算的數學模型。

3 影響腸道壁所受壓力的機器人結構參數的優化研究

圖2和圖3是螺旋機器人、腸道和黏液三維簡圖。影響腸道壁所受壓力的機器人結構參數主要有以下幾個(圖4):傾角α;槽面寬度a與螺距t之比β1;槽底寬度b與螺距t之比β2;螺旋槽槽深h與半徑間隙c之比γ;螺紋升角 Φ;螺紋線數n。在數值計算中,假定機器人與腸道間隙中充滿黏液,黏液密度為1.3g/cm3,動力黏度為1Pa˙s,腸道半徑R為6mm,螺旋機器人半徑r為4.5mm,正向和反向旋轉速度均為1000r/min,機器人軸向長度為 15mm,初設參數 α=60°,β1=β2=0.2,γ=1,Φ=45°,n=6。

利用Fluent6.3計算流體動力學軟件對機器人、腸道及黏液組成的系統進行三維數值模擬。首先,運用Pro/E wildfire5.0對機器人系統進行三維繪圖,然后輸出igs格式文件。利用CFD前置處理軟件Gambit2.2導入剛剛輸出的igs格式文件,對機器人系統進行網格劃分和邊界條件設置:機器人系統采用了非結構化四面體網格,機器人附近區域采用加密的計算網格;并設置腸道兩端為壓力進口和壓力出口,設置機器人鄰近區域流體類型為第一流體類型,設置腸道機器人附近區域之外的流體類型為第二流體類型,最后生成網格文件。利用Fluent6.3導入生成的網格文件,將單位轉換為mm,對網格進行光滑和交換,數值分析將雷諾平均的NS方程作為控制方程,湍流模型為標準k—ε模型,近壁處流動采用標準壁面函數處理,壓力和速度耦合方程采用標準SIMPLE算法求解。考慮到腸道壁的變形只引起腸道黏液流場區域的微小變化,并且機器人結構參數對腸壁損傷的影響也很小,因此忽略腸壁的柔性。為了模擬機器人鄰近區域流體的運動,采用滑移網格方法進行處理,給定轉速為1000 r/min(正轉)或—1000r/min(反轉)。解算收斂精度如下:腸道兩端進出口的流量誤差小于0.0001,x、y、z方向速度、k和ε的殘差收斂精度均為0.0001。

圖5所示為當螺旋機器人正向和反向旋轉運行時,傾角α對腸道壁所受最大正壓和最大負壓的影響曲線?？梢钥闯?隨著α的變化,機器人正轉和反轉時腸道壁所受最大正壓和最大負壓的絕對值相差不大,但變化規律相反。圖5中,最大正壓力和最大負壓力分別為 244.4Pa和—226.7Pa,這都在人體能夠承受的范圍內(腸道能夠承受的最大負壓力約為—3.7MPa[11])。

圖6所示為當螺旋機器人正向和反向旋轉運行時,參數β1對腸道壁所受最大正壓力和最大負壓力的影響曲線?？梢钥闯?隨著β1的增大,機器人正轉和反轉時腸道壁所受最大正壓力和最大負壓力的絕對值相差不大,不論正轉和反轉,腸道壁所受最大正壓力和最大負壓力的絕對值基本都隨之減小。圖6中,最大正壓力和最大負壓力分別為244.0Pa和—228.2Pa,這對人體腸道壁的影響是很小的。

圖7所示為當螺旋機器人正向和反向旋轉運行時,參數β2對腸道壁所受最大正壓力和最大負壓力的影響曲線?？梢钥闯?隨著β2的增大,機器人正轉和反轉時腸道壁所受最大正壓力和最大負壓力的絕對值相差不大,不論正轉和反轉,腸道壁所受最大正壓力和最大負壓力變化都很小。圖7中,最大正壓力和最大負壓力分別為245.5Pa和—229.2Pa,這也在人體能夠承受的范圍內。

圖8所示為當螺旋機器人正向和反向旋轉運行時,參數γ對腸道壁所受最大正壓力和最大負壓力的影響曲線。可以看出,隨著γ的增大,機器人正轉和反轉時腸道壁所受最大正壓力和最大負壓力的絕對值相差不大,不論正轉和反轉,腸道壁所受最大正壓力和最大負壓力的絕對值基本都隨之增大。圖8中,最大正壓力和最大負壓力分別為253.7Pa和—235.5Pa,這對人體腸道的影響是不大的。

圖9所示為當螺旋機器人正向和反向旋轉運行時,螺旋角Φ對腸道壁所受最大正壓力和最大負壓力的影響曲線?？梢钥闯?隨著Φ的增大,機器人正轉和反轉時腸道壁所受最大正壓力和最大負壓力的絕對值相差不大,不論正轉和反轉,腸道壁所受最大正壓力和最大負壓力的絕對值基本都隨之先增大后減小。螺旋角 Φ為45°時,腸道壁所受最大正壓力和最大負壓力都最大,分別為243.3Pa和—227.1Pa,這都在人體能夠承受的范圍內。

圖10所示為當螺旋機器人正向和反向旋轉運行時,螺紋線數n對腸道壁所受最大正壓和最大負壓的影響曲線?？梢钥闯?隨著n的增大,機器人正轉和反轉時腸道壁所受最大正壓力和最大負壓力的絕對值相差不大,不論正轉和反轉,腸道壁所受最大正壓力和最大負壓力的絕對值基本都隨之減小。螺紋線數n為6時,腸道壁所受最大正壓力和最大負壓力都最大,分別為243.3Pa和—227.1Pa。

4 結論

(1)所提出的螺旋微型機器人在人體腸道內懸浮運行時,腸道壁所受壓力都是人體能夠承受的,所以該機器人可對人體腸道實施無創或微創手術。

(2)不論正轉和反轉,腸道壁所受壓力的大小相差不大,并且機器人結構參數 α、β1、β2和 γ對腸道壁所受最大正壓力和最大負壓力的影響都不大,螺旋角 Φ和螺紋線數n對腸道壁所受最大正壓力和最大負壓力的影響較大。機器人正轉時,對于減小腸道壁所受最大正壓力來說,參數α、β1取較大的值,β2、γ、Φ取較小的值,n取8;對于減小腸道壁所受最大負壓力來說,參數α、γ、Φ取較小的值,β1取較大的值 ,β2取 0.3,n取 8。機器人反轉時,對于減小腸道壁所受最大正壓力來說,參數α、γ、Φ取較小的值,β1 取較大的值,β2 取 0.4,n取8;對于減小腸道壁所受最大負壓力來說,參數 α、β1取較大的值,β2、γ、Φ取較小的值,n取 8。

[1]Kassim I,Phee L,Wan S,et al.Locomotion Techniques for Robotic Colonoscopy[J].IEEE Engineering in Medicine and Biology Magazine,2006,25(3):49-56.

[2]Byungkyu K,Moon G L,Young P L,et al.An Earthworm—like Micro Robot Using Shape Memory Alloy Actuator[J].Sensors and Actuators,2006,125(2):429-437.

[3]Cheung E,Mustafa E K,Sukho P,et al.A New Endoscopic Microcapsule Robot Using Beetle Inspired Micro Fibrillar Adhesives[C]//IEEE/ASME International Conference on Advanced Intelligent Mechatronics.Monterey,2005:551-557.

[4]Byungkyu K,Sunghak L,Jong H P,et al.Design and Fabrication of a Locomotive Mechanism for Capsule—type Endoscopes Using Shape Memory Alloys(SMAs)[J].IEEE/ASM E Transactions on Mechatronics,2005,10(1):77-86.

[5]李佳峻,顏國正,王坤東,等.胃腸道內窺鏡機器人驅動機構研究[J].北京生物醫學工程,2009,28(2):179-183.

[6]Meron G D.The Development of the Swallowable Video—capsule(M2A)[J].Gastrointest Endosc,2000,52(6):812-819.

[7]周銀生,賀惠農,全永昕.無損傷腸道機器人運行速度的研究[J].摩擦學學報,1999,19(4):299-303.

[8]陳柏,蔣景榮,陳筍,等.螺旋式單鏈介入微機器人的研究[J].中國機械工程,2009,20(17):2100-2103.

[9]張永順,于宏海,阮曉燕,等.新型腸道膠囊式微型機器人的運動特性[J].機械工程學報,2009,45(8):18-23.

[10]Powell R L,Aharonson E F,Schwarz W H,et al.Rheological Behavior of Normal Tracheobronchial Mucus of Canines[J].Journal of Applied Physiology,1974,37(3):447-451.

[11]劉金華,王樹卿,劉丹陽.醫用微型機器人對腸道損傷的探討[J].世界華人消化雜志,2006,14(6):618-620.