朱佳鵬， 廖茂林， 朱志強， 曾子衿， 李 智

(北京科技大學(xué) 機械工程學(xué)院， 北京 100083)

隨著腸胃病患者基數(shù)不斷增加，消化道的醫(yī)療檢測需求量大幅增加，然而傳統(tǒng)的內(nèi)窺鏡檢查的管線探頭結(jié)構(gòu)，在檢測過程中易對患者造成強烈的不適感，并存在對病變區(qū)域造成二次傷害的安全隱患。為有效解決此類問題，無線式膠囊內(nèi)窺鏡逐漸引起關(guān)注，目前可作臨床使用的無線式膠囊內(nèi)窺鏡主要分為磁控式與被動式。其中，磁控式無線內(nèi)窺鏡膠囊依靠大型的磁控設(shè)備與膠囊表面的磁控衣作用，控制膠囊在消化道當(dāng)中的運動，并利用膠囊自身的圖像采集功能完成檢測[1]，但其費用昂貴，操作復(fù)雜。而被動式膠囊內(nèi)窺鏡雖然成本較低，操作簡單，但由于膠囊在消化道內(nèi)的移動僅依靠腸胃的蠕動，無法主動靠近或遠離病變區(qū)域[2]，給檢測帶來諸多不便的同時還存在無法排出體外的隱患。

為更好地解決無線式膠囊內(nèi)窺鏡的驅(qū)動問題，激振沖擊方式被引入對膠囊進行結(jié)構(gòu)改造。激振驅(qū)動膠囊系統(tǒng)由Liu等[3-4]在研究偏遠地區(qū)管道檢測問題時被提出，隨后被考慮用于消化道的無線檢測。其核心思想是在膠囊內(nèi)部設(shè)置一個激振沖擊結(jié)構(gòu)，通過內(nèi)部激振體與膠囊外殼的周期性沖擊碰撞驅(qū)動膠囊以一定速度進行自推進移動。在此基礎(chǔ)上，對具有單邊和雙邊約束結(jié)構(gòu)的自推進膠囊進行了對比分析[5-6]，發(fā)現(xiàn)單邊約束結(jié)構(gòu)更容易實現(xiàn)膠囊的高速移動，而雙邊約束結(jié)構(gòu)則在實現(xiàn)膠囊的雙向移動方面更具優(yōu)勢。

然而，膠囊內(nèi)窺鏡的主要運行環(huán)境為消化道，其小腸部分在變形時具有指數(shù)增強的強非線性應(yīng)力-應(yīng)變響應(yīng)[7]。因此，膠囊內(nèi)窺鏡在通過狹隘的小腸時，除了庫倫摩擦阻力，還會受到其他多種腸道阻力的影響。針對此類問題，Ciarletta等[8]利用超彈性理論對腸壁進行了分層分析，并提出了五應(yīng)力超彈性模型，該模型對腸道的力學(xué)模型搭建具有重要意義。隨后，Kim 等[9-10]基于Ciarletta 模型建立了膠囊在小腸內(nèi)運動的摩擦阻力模型，并利用該模型開展試驗，測量了不同尺寸的膠囊通過豬小腸內(nèi)部與表面時的摩擦阻力，最終得到其范圍在[10 mN，50 mN]。譚人嘉等[11]利用Ciarletta 模型，分別建立了膠囊機器人頭部壓力、中段摩擦力和頭部摩擦力的阻力模型，并通過試驗發(fā)現(xiàn)膠囊半徑與腸道半徑之比和腸道摩擦因數(shù)的增大都會導(dǎo)致膠囊前進阻力增大。Tian等[12]則從膠囊的安全性出發(fā)，用不同形狀與速度的膠囊依次通過平坦、塌陷、收縮和彎曲4種不同的小腸結(jié)構(gòu)，發(fā)現(xiàn)當(dāng)速度處于[4 mm/s，10 mm/s]，膠囊柱段的半徑等于膠囊末端的長度時，對腸道的損傷最小。此外，也有學(xué)者[13]從膠囊的移動速度出發(fā)，得到膠囊受到的摩擦阻力范圍在 [10 mN，250 mN]；而針對小腸收縮率[14]進行試驗得到了膠囊受到的收縮壓力峰值約為0.24 kPa。

因此，目前的問題主要集中在如何將建立好的腸道阻力模型與自推進膠囊內(nèi)窺鏡的動力學(xué)模型耦合，從而更好的對膠囊在小腸中的運動進行分析。為此，Yan等[15]建立了膠囊在小腸中運動的動力學(xué)模型，并以低能耗與高速度作為參考指標制定了膠囊的最優(yōu)參數(shù)，這也為本文阻力模型的搭建提供了參考思路。從腸道環(huán)境的不穩(wěn)定性出發(fā)，Liao等[16]利用組合優(yōu)化算法，得到了自推進膠囊在摩擦因數(shù)變化的腸腔內(nèi)運動時，5種不同運動的控制參數(shù)，包括前進(快慢)、后退(快慢)和懸停，并在此基礎(chǔ)上，Zhu等[17]對自推進膠囊在腸道中運動的能耗、速度與外部驅(qū)動力三者進行了多目標優(yōu)化。此外，Guo等[18-19]利用試驗對3種不同小腸環(huán)境下自推進膠囊的運動進行優(yōu)化分析，得到膠囊在腸道中最大前進速度為8.49 mm/s，最大的反向運動速度為-4.9 mm/s。

雖然國內(nèi)外學(xué)者們已經(jīng)開展了大量相關(guān)研究，但是大多并未考慮腸道的蠕動以及腸道內(nèi)流體環(huán)境對膠囊運動的影響。鑒于此，本文將搭建小腸-膠囊耦合模型，在小腸阻力分析的基礎(chǔ)上充分考慮小腸的能動性，分析腸道蠕動下自推進膠囊在通過狹窄小腸段時的動力學(xué)特性，討論膠囊的激振參數(shù)和結(jié)構(gòu)參數(shù)對自推進膠囊的非線性動力學(xué)行為演化規(guī)律的影響，以及分析流體環(huán)境對膠囊運動的影響。從而為實現(xiàn)膠囊在小腸中的最優(yōu)移動，并為自推進膠囊內(nèi)窺鏡的研發(fā)方案提供支撐。

1 小腸-膠囊動力學(xué)模型建立

1.1 自推進膠囊的動力學(xué)模型

為實現(xiàn)在不需要增加外部大型配套設(shè)備的情況下，膠囊在消化道中具有可控的移動能力，考慮在膠囊內(nèi)部設(shè)置磁致伸縮式的激振沖擊結(jié)構(gòu)來賦予膠囊自推進能力。自推激振膠囊的物理模型如圖1所示，質(zhì)量分別為mc，mm的膠囊外部殼體和內(nèi)部激振體通過一個Kelvin-Voigt(k1-c)模型相連，且膠囊內(nèi)部右側(cè)設(shè)置代表膠囊殼體剛度的彈簧k2。激振體在簡諧激勵Fa作用下發(fā)生振動時，會產(chǎn)生與外部殼體的相對位移Xm-Xc，當(dāng)相對位移大于兩者間的間隙e1時，激振體會與彈簧k2發(fā)生碰撞，這使得膠囊具有了接觸非線性特征。鑒于此，根據(jù)膠囊內(nèi)部激振體與彈簧k2是否發(fā)生碰撞，可以將膠囊分成以下兩個移動階段。

圖1 自推進膠囊基礎(chǔ)物理模型Fig.1 Physical model of self-propelled capsule

(1) 當(dāng)激振體與膠囊發(fā)生碰撞：Xm-Xc≥e1

此時，激振體mm的振動方程如式(1)所示

(1)

式中：Pd為激振振幅；f為激振頻率；Fc為膠囊外殼受到的內(nèi)部作用力，當(dāng)Fc大于外界阻力Ff時，膠囊整體開始運動，具體表達式為

(2)

(2) 當(dāng)激振體與膠囊未發(fā)生碰撞時：Xm-Xc

激振體的振動方程可用式(3)表示為

(3)

當(dāng)膠囊外殼受到的內(nèi)部作用力Fc大于外界阻力Ff時，膠囊整體開始運動，具體表達式為

(4)

在此基礎(chǔ)上，進一步引入Heaviside判定函數(shù)，對膠囊系統(tǒng)的各個移動階段的動力學(xué)方程進行整合

(5)

綜上，自推進膠囊的動力學(xué)模型可以表示為

(6)

其中

1.2 小腸模型

膠囊在吞咽后，會經(jīng)過人體的食道、胃、小腸和結(jié)腸，最后經(jīng)肛門排出；其中小腸作為主要的消化道器官，整體長度達4～6 m，其結(jié)構(gòu)最為復(fù)雜，對膠囊移動限制也最大。因此考慮對狹窄的小腸進行建模，以對自推進膠囊在其內(nèi)部的移動規(guī)律進行分析。

小腸的橫截面如圖2(a)所示，其結(jié)構(gòu)包括了漿膜、縱向肌肉層、環(huán)形肌肉層、粘膜層(含粘膜下層與粘膜)與上皮細胞層。從其功能方面來看，小腸依靠收縮環(huán)的不斷推進形成蠕動波完成對食物的運輸，過程中肌肉層與粘膜層對小腸的運動進行協(xié)調(diào)，如圖2(b)所示。由此可知，當(dāng)膠囊在腸腔內(nèi)運動時，會受到小腸結(jié)構(gòu)的應(yīng)變應(yīng)力與蠕動的作用，為此需要建立小腸的應(yīng)力應(yīng)變模型并且確定其蠕動速度，從而實現(xiàn)對小腸的建模分析。

(a) 小腸橫截面

(b) 小腸蠕動示意圖圖2 小腸結(jié)構(gòu)圖Fig.2 Structure of small intestine

1.2.1 應(yīng)力應(yīng)變模型

小腸結(jié)構(gòu)同時具備黏彈性與彈性兩種力學(xué)特征，其中粘膜層和上皮層是黏彈性的，而彈性主要來自肌肉層。但相比于彈性，黏彈性在描述小腸組織隨時間變化的力學(xué)特性方面更具優(yōu)勢，因而通常被用來代表小腸的特性。

因此，選擇小腸的黏彈性作為研究對象，通過引進三元麥克斯韋模型對小腸的黏彈性進行建模，如圖3所示，并得到小腸的應(yīng)力應(yīng)變關(guān)系表達式為

圖3 三元麥克斯韋模型Fig.3 The ternary Maxwell model

(7)

式中：ε(t)為小腸的應(yīng)變；E1，E2為模型中涉及的兩項楊氏模量；η為阻尼黏性系數(shù)；t為時間。

1.2.2 小腸蠕動速度

圖4 收縮環(huán)速度變化曲線Fig.4 Speed variation of shrink ring

1.3 小腸阻力模型

基于自推進膠囊的動力學(xué)模型與小腸模型，進一步分析兩者耦合作用下，膠囊在通過狹窄的小腸段時受到的外界阻力Ff，分別從包裹阻力，流體阻力與摩擦阻力3個方面構(gòu)建小腸-膠囊阻力模型。

1.3.1 包裹阻力模型

考慮到膠囊在通過小腸時會持續(xù)受到收縮環(huán)的作用，且膠囊直徑與小腸腸徑接近，因此膠囊在小腸內(nèi)運動時，全程基本處于整體包裹型的內(nèi)壁環(huán)境中，如圖5所示。針對小腸段的包裹阻力，建立了分段力學(xué)模型，將膠囊受到的腸道阻力分解為頭、中、尾3個區(qū)域，其中頭部受到來自小腸水平方向擠壓力形成的阻力Fhp，與垂直方向擠壓力形成的摩擦阻力Fhf。同理，尾部受到水平擠壓阻力Ftp，摩擦阻力Ftf，而膠囊的中部僅受到小腸的垂直方向擠壓力產(chǎn)生的摩擦阻力Fmf，如圖5所示。

圖5 包裹式內(nèi)壁環(huán)境下膠囊整體受力Fig.5 Force analysis of capsule in wrapped inner wall

(a)

(b)圖6 膠囊頭部受力分析Fig.6 Force analysis of capsule head

(8)

并通過引入式(7)可以得到相應(yīng)的小腸應(yīng)力為

(9)

基于小腸應(yīng)力的分析，進一步推算膠囊受到的擠壓力，如圖6所示，首先分析小腸與膠囊擠壓后的變化，厚度由初始值hi變?yōu)閔ie，由式(10)表示為

(10)

結(jié)合式(9)與式(10)，可以得到膠囊頭部區(qū)域受到的壓強表達式

(11)

引入dS(x)表示為某一位置的環(huán)向受力面積，從而計算出膠囊頭部任意區(qū)域受到的擠壓力dp(x)，表達式為

dp(x)=q(x)dS(x)=q(x)2πR(x)·

(12)

式中，R′(x)為R(x)的導(dǎo)數(shù)。

得到作用于膠囊壁面的擠壓力后，將其分解為水平與垂直方向的力，開始討論膠囊受到的外界阻力，其中水平方向的擠壓阻力Fhp為

(13)

式中，cosφ(x)=(xc-x)/Rc，垂直方向擠壓力引起的水平摩擦阻力Fhf由式(14)表示

(14)

式中，μ為小腸環(huán)境的摩擦因數(shù)。

同理可得，膠囊尾部的壓力導(dǎo)致的水平擠壓阻力為

(15)

此時的cosφ(x)=(xc+L-x)/Rc<0，L為膠囊中部柱段長度，x∈[xc+L,2xc+L]，為膠囊尾部受到的水平擠壓阻力與頭部相反，膠囊尾部垂直方向的壓力引起的Ftf水平摩擦阻力為

(16)

膠囊中部柱段受到的壓力與x方向垂直，因此只產(chǎn)生水平方向的摩擦阻力，其表達式為

(17)

綜合式(13)～式(17)，其中水平方向的頭部擠壓阻力Fhp與尾部擠壓阻力Ftp相互作用下結(jié)果為零，最終可以得到包裹阻力Ffi

Ffi=Ftf+Fmf

(18)

1.3.2 流體阻力模型

除去小腸包裹作用引起的阻力，膠囊自身的移動速度與小腸的蠕動速度會形成一定的相對速度，加上腸腔表層流體環(huán)境的影響，如圖7所示，膠囊在小腸中移動時全程會受到流體阻力Ffq作用，具體表達式[21]為

圖7 流體阻力與庫倫摩擦阻力Fig.7 Fluid resistance and Coulomb friction resistance

(19)

式中，d為流體厚度。

1.3.3 庫倫摩擦阻力模型

此外膠囊在小腸環(huán)境中還受到自身重力帶來的庫倫摩擦阻力Ffc，如圖7所示，表達式為

Ffc=Mmμg

(20)

式中:g為重力加速度；Mm為膠囊系統(tǒng)的總質(zhì)量。

最終通過以上3個模型分析，可得到膠囊在小腸中運動時受到的外界阻力Ff可以表示為

Ff=Ffi+Ffq+Ffc

(21)

1.4 自推進膠囊的動力學(xué)模型

圖8 小腸-膠囊物理模型Fig.8 Small intestine-capsule physical model

可見，相比于初始的膠囊模型，小腸-膠囊模型不僅要考慮碰撞導(dǎo)致的接觸非線性特征，還需要對小腸蠕動速度與膠囊自推進速度之間的相對速度變化引起的阻力Ff方向進行判定。為此，需要在激振膠囊模型的基礎(chǔ)上引入兩個新的Heaviside函數(shù)，以考慮小腸與膠囊相對速度的影響

(22)

結(jié)合腸壁包裹環(huán)境下對應(yīng)的Ff表達式(21)，最終建立小腸-膠囊的耦合動力學(xué)模型

(23)

此外，為更好地了解自推進膠囊在小腸環(huán)境下移動速度變化規(guī)律，計算膠囊在單位時間內(nèi)的平均速度

(24)

式中：N為周期數(shù)；T=1/f為單個激勵周期。

2 非線性動力學(xué)分析

利用建立的小腸-膠囊耦合動力學(xué)模型，進一步分析自推進膠囊在小腸環(huán)境中的非線性動力學(xué)行為特性。根據(jù)前人[22]對小腸結(jié)構(gòu)的試驗研究，首先確定了小腸參數(shù)，如表1所示。并從膠囊的應(yīng)用要求出發(fā)，初步確定了膠囊的結(jié)構(gòu)參數(shù)與激振參數(shù)，激振參數(shù)主要為式(1)中激振力Fa中的頻率、振幅，如表2所示。選擇膠囊的激振參數(shù)(頻率、振幅)與部分結(jié)構(gòu)參數(shù)(激振體的質(zhì)量與剛度)作為研究對象，通過大范圍調(diào)節(jié)各項參數(shù)，得到膠囊在小腸中運動的非線性動力學(xué)行為演化規(guī)律。具體而言，對膠囊進行分岔分析、相平面分析，以及速度變化對比，并通過確定自推進膠囊最優(yōu)速度出現(xiàn)時的系統(tǒng)動態(tài)響應(yīng)，為制定合理的激振控制方案與膠囊結(jié)構(gòu)設(shè)計提供支持。

表1 小腸參數(shù)表Tab.1 Parameters of small intestine

表2 膠囊參數(shù)表Tab.2 Parameters of capsule

2.1 激振參數(shù)的影響

采用四階龍格庫塔算法對建立的小腸-膠囊耦合動力學(xué)模型進行數(shù)值求解，計算200個激勵周期，并取最后10個周期作為參考，以確保自推進膠囊處于穩(wěn)態(tài)響應(yīng)狀態(tài)。

為進一步分析流體阻力對膠囊運動帶來的影響，分別在圖9(a)、圖9(b)中分別比較了有無流體阻力膠囊的平均速度曲線與分岔圖，如圖10所示。無流體阻力時膠囊的速度得到小幅提升，在振幅為60 mN時，兩者速度差最大，且膠囊系統(tǒng)的分岔區(qū)間由之前[21 mN,29 mN]縮減至[23 mN,26 mN]。

圖9 關(guān)于振幅Pd的動力學(xué)分析(圖9(c)～圖9(h)中虛線表示膠囊外殼，實線表示激振體；圖9(i)～圖9(n)中直線為撞擊面位置，圓點為龐加萊截面截取的點)Fig.9 Dynamic analysis of amplitude

接著，以頻率f作為變量，參數(shù)區(qū)間為[0,50 Hz]，參數(shù)間隔取0.05 Hz，得到膠囊的平均速度變化曲線圖圖10(a)與分岔圖圖10(b)，與Pd不同，在f的變化區(qū)間內(nèi)，膠囊系統(tǒng)始終保持單周期運動，其中，當(dāng)頻率達到10 Hz時，激振體開始與殼體碰撞，見圖10(i)，此時膠囊的速度為29.9 mm/s，與小腸蠕動速度基本保持一致，見圖10(a)。當(dāng)頻率增大到15 Hz時，見10(j)，單個激勵周期下激振體將與殼體發(fā)生兩次碰撞，值得注意的是，此時膠囊平均速度vavg下降到21.3 mm/s，見圖10(d)，相對于小腸蠕動速度30 mm/s，膠囊在往反方向移動，這表明選擇合理的激振頻率可以控制膠囊的雙向移動。繼續(xù)增大頻率至17 Hz時，見圖10(k)，激勵周期內(nèi)的碰撞次數(shù)重新回到一次，且由圖10(c)可知，此時膠囊的平均速度達到最大為39.3 mm/s。當(dāng)頻率處于25～35 Hz時，以30 Hz，35 Hz為例，vavg基本穩(wěn)定，速度保持在36 mm/s附近，且碰撞次數(shù)保持為一次，見圖10(f)、圖10(g)、圖10(l)、圖10(m)。一直到f=36 Hz時，見圖10(n)，激振體與殼體之間碰撞消失，膠囊靜止僅隨腸道蠕動前進，速度為30 mm/s，見圖10(h)；由此可知，當(dāng)激勵頻率過高時，反而無法推動膠囊高效移動。

2.2 結(jié)構(gòu)參數(shù)的影響

為保證自推進膠囊在小腸內(nèi)的高效穩(wěn)定移動，不僅要對激振參數(shù)進行調(diào)節(jié)，還需合理設(shè)計膠囊的結(jié)構(gòu)參數(shù)；因此，本文以膠囊激振體的質(zhì)量mm與剛度k1作為研究對象，分析結(jié)構(gòu)參數(shù)對自推進膠囊的動力學(xué)行為演化規(guī)律的影響。

首先，分析膠囊激振體的質(zhì)量mm，其參數(shù)區(qū)間設(shè)為[1 g,10 g]，參數(shù)間隔為0.05 g，結(jié)果如圖11所示。觀察發(fā)現(xiàn)，m1在區(qū)間內(nèi)變化時，激振膠囊始終保持單周期運動，見圖11(b)。初始階段，當(dāng)mm=1.1 g時，膠囊速度為22.82 mm/s，相對于小腸往后運動，系統(tǒng)單個激勵周期內(nèi)發(fā)生了兩次碰撞，見圖11(c)、圖11(i)。當(dāng)mm=1.4 g時，由圖11(d)、圖11(j)可知，系統(tǒng)變?yōu)閱沃芷趩闻鲎矤顟B(tài)，膠囊平均速度速度達到38.29 mm/s，移動方向轉(zhuǎn)為相對小腸向前運動。持續(xù)增大mm至7 g，vavg呈現(xiàn)先減小再增大至穩(wěn)定的變化趨勢，過程中整體變化幅度較小，膠囊平均速度最終維持在41 mm/s附近，系統(tǒng)維持單周期單碰撞響應(yīng)狀態(tài)，見圖11(e)、圖11(g)與圖11(k)、圖11(m)。當(dāng)mm超過7 g時，見圖11(h)、圖11(n)，vavg變?yōu)榧s30 mm/s，此時，激振體與殼體不再發(fā)生碰撞，膠囊只隨腸道蠕動；主要原因是當(dāng)激振體質(zhì)量持續(xù)增大時，將導(dǎo)致庫倫摩擦力過大，激振沖擊力小于外界阻力，膠囊無法相對小腸移動。另外，對激振質(zhì)量塊質(zhì)量變化時，膠囊的多穩(wěn)態(tài)進行了分析，見圖11(a)、圖11(b)，令mm由高到低進行變化，由圖11(f)、圖11(l)可知系統(tǒng)的階躍點由之前的7 g變?yōu)?.9 g，且存在明顯的遲滯回線。

接著，分析激振體的剛度k1對膠囊系統(tǒng)動態(tài)響應(yīng)的影響，k1的變化區(qū)間設(shè)置為[10 N/m,120 N/m]，參數(shù)間隔為0.05 N/m，結(jié)果如圖12所示。特別地，剛度k1過小時系統(tǒng)會出現(xiàn)混沌現(xiàn)象，見圖12(i)，此時膠囊的平均速度為32.5 mm/s，見圖12(c)。當(dāng)k1=12 N/m時，系統(tǒng)響應(yīng)為的單周期單碰撞狀態(tài)，速度達到最大為44.9 mm/s，見圖12(d)、圖12(j)。隨著k1繼續(xù)增大至68 N/m，vavg逐步下降至30 mm/s，見圖12(e)～圖12(g)，響應(yīng)保持在單周期單碰撞，但碰撞逐漸減弱，見圖12(k)～圖12(m)。當(dāng)k1>68 N/m時，碰撞消失，膠囊的移動速度穩(wěn)定在30 mm/s，即僅依靠腸道蠕動實現(xiàn)前進，見圖12(h)～圖12(n)。

圖12 關(guān)于激振塊剛度k1的動力學(xué)分析(圖12(c)～圖12(b)中虛線表示膠囊外殼，實線表示激振體；圖12(i)～圖12(n)中豎線為撞擊面位置，圓點為龐加萊截面截取的點)Fig.12 Dynamic analysis of oscillator stiffness

4 結(jié) 論

(1) 本文通過考慮小腸的應(yīng)力應(yīng)變特征與其蠕動規(guī)律，建立了小腸-膠囊耦合動力學(xué)模型，重點討論了膠囊的振幅、頻率、激振體剛度與質(zhì)量4項參數(shù)對膠囊在小腸內(nèi)移動的影響。由于所建立的動力學(xué)模型耦合了自推進膠囊自身的接觸非線性以及小腸阻力的非線性特征，因而在系統(tǒng)的非線性動力學(xué)行為規(guī)律研究中發(fā)現(xiàn)了倍周期分岔、逆倍周期分岔、混沌等現(xiàn)象，它們對膠囊的自推進移動影響顯著。

(2) 當(dāng)調(diào)節(jié)激振參數(shù)時，膠囊的振動形式主要表現(xiàn)出單周期單碰撞、單周期雙碰撞、雙周期單碰撞、雙周期雙碰撞4種狀態(tài)；從膠囊的平均移動速度來看，最優(yōu)的前進速度都出現(xiàn)在單周期單碰撞的振動狀態(tài)下，原因在于單周期單碰撞能更高效地實現(xiàn)能量集中并通過激振體對殼體的唯一一次碰撞，最大程度地將激勵能量轉(zhuǎn)化為驅(qū)動沖擊力。此外，在給定的參數(shù)范圍內(nèi)，當(dāng)振幅增大時，膠囊始終保持相對小腸正向運動，膠囊的絕對速度持續(xù)增加并最終穩(wěn)定在44 mm/s左右；而頻率變化過程中會出現(xiàn)膠囊相對小腸反方向移動的情況(f處于15 Hz附近時)，此時膠囊處于單周期雙碰撞狀態(tài)。此外，頻率過高時，會導(dǎo)致激勵驅(qū)動失效，膠囊表現(xiàn)為相對靜止，只隨小腸蠕動前進。綜上，在調(diào)節(jié)膠囊激振參數(shù)時，應(yīng)保證系統(tǒng)在穩(wěn)定的單周期單碰撞狀態(tài)；而要使膠囊實現(xiàn)反向移動，則需要對頻率進行調(diào)整，使膠囊處于單周期雙碰撞狀態(tài)；此外，還應(yīng)注意避免由于激勵頻率過高而導(dǎo)致的驅(qū)動失效。

(3) 相比于激振參數(shù)，改變結(jié)構(gòu)參數(shù)時膠囊基本保持在單周期單碰撞、單周期雙碰撞與混沌3種狀態(tài)下，其中混沌與單周期雙碰撞分別只出現(xiàn)在k1與mm極小時，膠囊的最優(yōu)前進速度分別出現(xiàn)在mm=7 g與k1=12 N/m附近時，系統(tǒng)狀態(tài)都為單周期單碰撞狀態(tài)。此外，膠囊的反向移動僅出現(xiàn)在mm極小時，此時系統(tǒng)為單周期雙碰撞狀態(tài)，與調(diào)節(jié)激振頻率時發(fā)現(xiàn)的規(guī)律一致。此外，在對膠囊結(jié)構(gòu)進行設(shè)計時，應(yīng)避免激振體質(zhì)量過大或剛度過大而導(dǎo)致的驅(qū)動失效。

(4) 本文現(xiàn)有的分析與結(jié)論，都建立在對膠囊的單參數(shù)分析之上，接下來將考慮對膠囊多個參數(shù)同時進行分析并進行多目標優(yōu)化，在平均速度的基礎(chǔ)上增加能耗、以及腸道損傷度等優(yōu)化目標；并進行可靠性分析，以此來制定更為合理的激振控制方案與膠囊結(jié)構(gòu)設(shè)計，從而更好地指導(dǎo)后續(xù)的試驗研究。