柴川頁 ,顏國正 *,韓 玎 ,趙 凱 ,華芳芳
(1. 上海交通大學電子信息與電氣工程學院,上海200240;2. 上海交通大學 醫療機器人研究院,上海 200240)
大便失禁指肛管和括約肌對糞便失去蓄控力,嚴重影響患者的生活。人工肛門括約肌(Ar?tificial Anal Sphincter,AAS)是可治療大便失禁的新興手段之一,通過原位植入腹腔代替功能缺失的括約肌實現控便與排便動作。夾持機構是AAS 的關鍵組成部分,通常環繞直腸肛管并與腸道外壁緊密接觸,其運動學和力學特點直接影響到患者體驗。
Yan 團隊在AAS 動物實驗中發現夾持機構被大量組織增生包裹,而且植入2 個月后機構出現嚴重的機械故障[1],這些現象說明夾持機構的生物相容性不足。AAS 的材料選擇直接影響生物相容性,長期植入式醫療設備廣泛使用的材料有醫用硅橡膠、鈦合金、尼龍等。Artificial Bowel Sphincter(ABS)、Soft Anal Band(SAB)和 Ger?many Artificial Sphincter(GAS)等就是硅橡膠環形袖套式AAS[2-6]。硅橡膠輕軟、生物相容性好,但因強度不足僅用于形狀簡單的結構。Magnet?ic Anal Sphincter(MAS)是醫用鈦合金磁珠式AAS[7],其強度高于硅橡膠,但磁珠間被動收緊力較小導致最大控便量不足。上述AAS 采用環形肛門括約肌擠壓腸道的仿生學控便機理,受夾持腸壁將形成大量褶皺,致使腸道局部高壓產生血供安全隱患,80%的患者須再次手術治療[8]。Marziale 等人通過數值模擬和實驗評估了周向與鉗夾式人工尿道括約肌的力學特點[9],Zan 等人也通過數值模擬的方法分析了一種夾鉗式AAS的生物力學相容性[10],但僅停留在理論研究階段而未經過工程驗證。Zhou 等研發了擺臂式仿恥骨 AAS[11],Wang 等研發了恒力彈壓式 AAS[12],但由于腸道的粘彈性及機構并非完全封閉,夾持腸道時常出現移位、膨出的現象,影響控便性能。
除了生物相容性,控排便響應時間也是AAS夾持機構的關鍵指標之一。傳統AAS 響應時間較長,執行效率低下,液壓驅動的ABS 耗時5~8 min[2],基于形狀記憶合金的 Anal Sphincter-Shape Memory Alloy 開閉 時間約 5 min[13]。采用 電機驅動后,Doll 等設計的 GAS 響應時長達 8~9 min[3],Ke 等設計的機械式假體響應時間為 80 s[14]。
上述分析可見,現有AAS 夾持機構生物相容性較低、響應時間較長,本文在人體排便機制及其力學特點研究的基礎上,提出一種新型夾持機構,并設計出由三組交錯疊放的夾持臂構成的封閉式括約肌假體(Novel Closed-link Artificial Anal Sphincter,NCAAS)。基于新型假體結構,本文進行了運動學與力學理論分析,確定設計參數,并通過有限元仿真說明了機構控便有效性和安全性。樣機的制作依照理論計算及仿真結果,并實施了控排便時長實驗與夾持性能實驗,實驗結果亦驗證了夾持機構的有效性及安全性。NCAAS 兼有肛門括約肌封閉式收縮與恥骨直腸肌形成肛腸角的仿生學特點,可滿足人體日常所需控便量,具備良好的生物相容性,而且與傳統AAS 相比較大縮短了響應時間。
人體在靜息控便時,隨著糞便積累使直腸膨脹、直腸內壁受體對糞便量采樣,判斷是否需要排便。排便行為的執行取決于個體腹壓、結腸動力、直腸擴張性、肛腸角以及盆底肌之間的協調[15]。恥骨直腸肌包繞直腸肛管,另一端連接于恥骨末端。如需排便,恥骨直腸肌、括約肌放松,直腸伸直,肛腸角為100°~110°,肛門打開排出糞便。待排便結束后,內括約肌、恥骨直腸肌收縮,肛門關閉,肛腸角減小至 80°~90°恢復控便[16]。文獻[17]表明盡管各種哺乳動物的直腸長度在4~40 cm 內變化,但排便時間維持在12±7 s。
腸道作為生物組織因膠原特性具有非線性機械行為、明顯的粘彈性。為簡化腸道模型,研究中通常忽略腸道的黏膜皺襞,將其看作柱狀壓力管,故因糞便積累產生的壓力不僅使腸道管道在圓周方向擴展,同時引起管道長軸方向伸展、徑向收縮,并致使管壁變薄[18]。因此,研究中采用三維Ogden 不可壓縮超彈模型對腸道組織材料建模,密度1 020 kg·m-3,應變勢能為式(1):

其中:λ1,λ2,λ3分別為主應變伸長率,μi,αi為材料常數,取值如下[19]:μ1=8.3 kPa,μ2=0.3 kPa,μ3=6.2 kPa,α1=7.625,α2=13.875,α3=7.625。
NCAAS 夾持機構主要包括夾持臂和傳動機構,傳動機構用于驅動夾持臂運動。設計夾持機構時主要考慮以下要求:
(1)AAS 長期植入性要求夾持機構在滿足功能的條件下,各組件集成小型、緊湊;夾持臂材料輕便、生物相容性良好;傳動機構具有自鎖性,即保持某一夾持狀態的能力。
(2)為順利包繞直腸,夾持臂完全張開時的內切圓直徑應不小于松弛直腸的平均直徑D=20 mm;要減少夾緊腸道造成的大量褶皺,夾持臂完全閉合時的長度應不小于空癟腸道的寬度
(3)將夾持機構設置勻速夾持,排便時間設計為哺乳動物的排便時長;否則,會損傷腹腔內生物組織或導致排便效率低下。
(4)要求能至少控制4.6 kPa的糞便量[14],夾持臂封閉夾緊腸道時肛腸角不大于 90°[16];同時,為保證腸道血供安全閾值,夾持壓強不應大于10 kPa[10]。
如圖1 所示,力學建模中將夾持臂簡化為曲柄滑塊機構進行剛體運動學分析[20]。

圖1 運動學與力學分析Fig.1 Kinematics and mechanical analysis
模型以驅動器輸出軸O為原點創建直角坐標系Oxy,定義夾持臂各鉸接點坐標A(x1,y1),C(x2,y2),D(x3,y3)。建立鉸接點坐標矩陣的數學模型,如式(2):

由 于D點 沿y軸 運 動 ,可 得β和λ的 關 系 ,如式(3):

夾持臂根部A為固定點,設|x1|=c,y1=e;c,e為常數。各鉸接點坐標矩陣數學模型,如式(4):

夾持臂運動時角速度ωA恒定,且與機構夾持臂張角λ滿足式(5):

夾持機構長度d以及D點的速度vD和加速度aD,如式(6)~式(8):


3.2.1 夾持力方程
因夾持臂運動速度平穩緩慢,故可將夾持臂夾持直腸肛管作為靜力系統進行力學分析。取圖1(b)中心線左側機構為研究對象,根據虛功原理,驅動夾持臂運動的扭矩TA經微小時間增量Δt轉過微小角度增量 Δθ=ωAΔt,該過程所作的虛功為WA。考慮到腸道生物組織具有柔性與粘彈性,夾持臂與腸道外壁作用時受力可視為均布載荷,并作用在位于s 區域的內。扭矩TA的虛功WA等于 桿AC 所 受 夾持 反 力q1與 桿CD 所受 夾持反力q2在Δθ內所作虛功之和W1+W2。根據虛功原理,可得式(9):

桿 AC 所作虛功W1,如式(10):

以下說明如何計算桿CD 所作的虛功W2:
將 C,D 的微位移 ΔlC,ΔlD沿平行于桿 CD 和垂直于桿CD 的方向分解,如式(11)~式(12):

由于桿件不會被拉長,故滿足式(13):

桿 CD 上任意點P的虛功 ΔWP,如式(14):


將式(16)代入虛功方程(10),有:

將式(17)代入(1)、(3)、(4)、(7),故夾持力載荷方程如式(18):

3.2.2 夾持臂有限元分析
實際上夾持臂具有一定曲度,故前文將夾持臂視為連桿機構分析夾持腸道的過程存在一定誤差,因此,本節基于Ansys 有限元分析對夾持機構夾持腸道的動作做瞬時動力學仿真。醫用剛性材料中鈦合金密度4 500 kg·m-3是醫用尼龍的4 倍,建模過程中考慮到夾持機構的輕量性,因此將夾持臂材料設置為醫用尼龍PA6,密度1 130 kg·m-3,楊氏模量 2.32 GPa,泊松比 0.34,三層夾持臂間距設置為4 mm;腸道簡化為外徑20 mm 厚度 2 mm 的圓柱管,材料根據公式(1)設置。如圖2 所示,網格劃分方法為Tetrahedrons,夾持臂與腸道接觸面網格單元尺寸設置為0.002 m,節點數為25 220,單元數為12 908。考慮真實情況下機構夾持腸道的方式與腸道變形情況,本文設置了6 處操作臂-腸道接觸面,1 處腸道-腸道接觸面,摩擦系數均設為0.1;操作臂之間通過joint 連接,在夾持臂根部設置3 處Body-Ground 轉動副,余下5 處鉸鏈設置為Body-Body轉動副。固定腸道上、下兩個端面,在Body-Ground 處施加 3(°)/s 的恒定角速度載荷,運動12 s 后停止,求解時間設置為15 s。夾持過程腸道和夾持臂位移云圖如圖3 所示,夾持機構完好,直腸變形由夾持臂控制,在t=9.09 s 時腸道基本閉合。

圖2 網格劃分與施加荷載示意圖Fig.2 Schematic diagram of meshing and loading

圖3 夾持過程位移云圖Fig.3 Displacement during clamping
為了驗證張角范圍內是否能滿足控便性能,本文對夾持過程中腸道形成肛腸角的過程也進行了分析。如圖4 所示,t=9.09 s 時,肛腸角為95.2°達到近似控便狀態;t=11.667 s 時,肛腸角為62.2°,滿足生理控便的要求。

圖4 仿真肛腸角Fig. 4 Simulated anorectal angle
如圖5 所示,腸道處于夾緊狀態時,高壓區僅集中出現于腸道夾持區域兩側,腸道表面平滑。此外,如圖6 所示,夾持臂與腸道的接觸力隨夾持過程逐漸增大,在θ=32°時存在較大抖動。原因在于此時腸道幾乎為夾緊狀態,而腸道組織盡管存在粘滯性,其體積在受壓過程中近乎不變。仿真結果還顯示,接觸力在極限位置達到最大,約為1.6 N;腸道血供壓強閾值10 kPa 與接觸面積200 mm2之 積 為 血 供 壓 力 閾 值 2 N[10],1.6 N<2 N 說明夾持機構能夠保障腸道血供安全。

圖5 腸道應變云圖Fig.5 Colon strain

圖6 操作臂-腸道接觸力與擺桿轉角的關系Fig.6 Relationship between arm-colon contact force and swing rod angle
已知醫用尼龍PA6 的屈服強度σs=30 Mpa,設安全系數[s]=1.1,則安全應力σ[s]=σs/[s]=27.27 MPa。夾持臂應力仿真得到最大應力σmax=6.298 9 MPa<σ[s],因此,醫用尼龍 PA6材料符合夾持臂安全應力強度要求。
3.2.3 傳動機構設計
已知夾持臂旋轉點位置固定,跟據張角范圍所設計的雙軸驅動傳動機構包含多級行星齒輪減速器以及基于溝槽凸輪擺桿結構。為增大原始電機的輸出扭矩,減速器由總傳動比為3 802的5 級行星齒輪構成,每級行星齒輪有三個行星輪和一個太陽輪。此外,多級齒輪還增大了摩擦力,使反行程時驅動力不足以克服摩擦力矩,從而實現“自鎖”。基于溝槽凸輪擺桿的擺動轉盤是傳動機構的關鍵:減速器輸出軸與擺動轉盤中心重合,圓端銷在外部約束為溝槽的凸輪輪廓里滑動,通過雙軸驅動兩側夾持臂相向運動。因此需要兩個溝槽凸輪,且二者組合在平面擺動轉盤上。
記擺桿初始擺角θ0,凸輪初始角為φ0,擺桿長度l,凸輪擺桿中心距a,凸輪升程角為φ,擺桿轉角θ,可根據幾何關系通過解析法計算凸輪的理論輪廓線。因機構勻速夾持,傳動機構的擺桿轉角θ與凸輪轉角φ滿足式(19),其中C為常數:

3.2.4 設計結果
表1 給出夾持機構的設計參數。由式(7)、(8)得頂點D的速度和加速度,如圖 7 所示。D的速度與張角近乎呈線性從2.65 mm·s-1減少至0.15 mm·s-1,R2=0.9958;加速度aD不超過 6.3 mm·s-2,不至于損傷體內生物組織。NCAAS 夾持機構設計結構如圖8 所示。

圖8 夾持機構Fig.8 Clamping mechanism

表1 設計參數表Tab.1 Design Parameters

圖7 頂點速度,加速度與夾持臂張角的關系Fig.7 Relationship between vertex velocity,acceleration and opening angle of the clamping arm
NCAAS 樣機如圖9 所示,機構主體由醫用尼龍PA6 材料經3D 打印制成,高度42.7 mm,長度最大為68.2 mm,夾持臂張口內徑為10.1~23.8 mm。執行機構質量為55.19 g,表面覆蓋一層2 mm 厚生物相容性良好的硅膠薄膜,以減少植入體內導致的組織增生,并能提升夾持腸道的柔性,減少醫用尼龍剛性材料直接接觸腸道造成的損傷。

圖9 NCAAS 樣機Fig.9 NCAAS prototype
夾持機構由從上至下依次設置三組交錯疊放的夾持臂和基于溝槽凸輪擺桿的傳動機構構成,每組夾持臂包含通過鈦合金螺釘鉸接的兩個操作臂。夾持臂近端分別固定于兩根軸上,上、下層夾持臂連于同一根軸,中層夾持臂連接另一根軸并置于另側,夾持臂間縱向間隔4 mm。為實現封閉功能,三組夾持臂遠端通過醫用鈦合金螺釘一并鉸接。經實驗驗證,NCAAS 機構只能由凸輪驅動擺桿,無法反向驅動;此外,與僅依靠齒輪傳動的方法[21]相比,其受力面、受力強度更大。
使用秒表記錄3 次執行機構的控便、排便時長,二者分別對應夾持機構的夾持、張開時間,結果如表2 所示,控排便平均時間為7.25 s。

表2 排便和控便時長記錄表Tab.2 Defecation and defecation control time record
為驗證NCAAS 的夾持性能,本文實施了機構夾持豬大腸的離體實驗。考慮到人體直腸肛管的特征,離體實驗裝置如圖10(a)所示,將一段直徑約20 mm、厚度3 mm、長度200 mm 的新鮮豬大腸的上端套入波紋管并用扎口帶密封,波紋管豎直固定在鐵架臺上,鐵架臺下方有一塊中間有圓孔的支撐臺。NCAAS 樣機放置于圓孔上方,夾持豬大腸下端保持夾緊控便的狀態。實驗時,向波紋管頂注入自來水使腸腔膨脹、增壓以模擬糞便積累,糞便量以100 g 為間隔通過注水量模擬,等待腸道及機構狀態穩定時記錄注水量與液體滲漏情況。

圖10 夾持性能實驗圖Fig. 10 Clamping performance experimental
文獻[14]表明糞便積累到 110 g、腸腔壓力約2.5 kPa 時人體將產生便意;而當糞便增至220g、腸壓增至4.6 kPa 時內括約肌將失去協助控便能力,人體必須排便。表3 結果表明NCAAS 樣機在糞便量達700 g 時發生少量滲漏,因此最大控便量約為700 g,滿足人體日常生理需求。如圖10(b)所示,NCAAS 樣機閉合時腸道表面保持平滑、無明顯褶皺,該現象說明夾持機構尺寸選擇合理、安全;如圖 10(c)所示,NCAAS 控便時形成的肛腸角不大于90°。

表3 夾持性能實驗記錄表Tab.3 Clamping performance test record
本文提出了一種AAS 夾持機構優化設計的方法:首先建立夾持機構的幾何模型和數學模型,對相關力學性能進行理論推導;然后通過有限元分析對機構優化設計,驗證理論分析可行性;最后通過實驗測試機構的有效性。該方法對機械結構的設計與優化具有普適性。
基于該方法,本文設計了一種能封閉包繞直腸和肛管的NCAAS 夾持機構,機構由基于溝槽凸輪擺桿傳動的三層醫用尼龍夾持臂組成,表面覆一層醫用硅膠膜,具有較好的生物相容性和較短的響應時間。夾持腸道過程中,本文基于虛功原理推導了夾持力方程,基于有限元進行瞬時動力學仿真,仿真結果顯示近似控便狀態時肛腸角為95.2°,夾緊時肛腸角為62.2°,說明了該夾持機構控便的可行性;而且夾持臂在腸道上的作用力1.6 N 小于腸道破裂壓力,具備一定生物力學相容性。質量55.19 g 的NCAAS 樣機高度42.7 mm,長度不超過68.2 mm,控排便平均時長7.25 s,較大縮短了AAS 的響應時間。離體實驗結果表明NCAAS 控便量達700 g,肛腸角小于90°,腸道夾持區域平滑,因此能夠滿足人體控便要求,達到了設計優化的目的。
但本文所述NCAAS 樣機所實施的離體實驗采用自來水模擬糞便,而實際糞便多為半固態。實驗過程中,夾持力數據收集不足,因此實驗結果具有一定局限性。NCAAS 最終將原位植入于潮濕、復雜的體內環境,今后工作將聚焦于動物實驗從而驗證機構長期植入的有效性。