一種流體驅(qū)動式模塊化軟體仿生象鼻關(guān)節(jié)變剛度結(jié)構(gòu)設(shè)計及仿真分析

劉思遠(yuǎn)， 張京濤， 王 杰， 鄧亦杰， 辛闖庭

(燕山大學(xué) 河北省重型機械流體動力傳輸與控制重點實驗室， 河北 秦皇島 066004)

引言

軟體抓持裝置是兼具位姿變換功能的運動軟體夾持器，實現(xiàn)了自身抓持結(jié)構(gòu)特性與運動功能的巧妙融合[1]。由于結(jié)構(gòu)靈活、與環(huán)境交互性好，在國防軍事、航空航天、消防救援、人機交互、醫(yī)療外科、康復(fù)助力器等多個關(guān)鍵應(yīng)用領(lǐng)域發(fā)揮著重要作用，是目前軟體機器人研究的熱門方向[2-3]。

國內(nèi)外眾多科研學(xué)者[4-5]針對軟體抓持裝置多自由度變形的順應(yīng)性和承載特性開展了大量的研究工作。在順應(yīng)性方面，CONNOLLY F等[6]提出通過改變纖維角度來調(diào)整制動器以實現(xiàn)延伸、徑向膨脹和扭轉(zhuǎn)等運動。YAN J等[7]基于氣動驅(qū)動螺旋腔的原理，設(shè)計了一種新型的扭轉(zhuǎn)執(zhí)行器模塊可在復(fù)雜和狹窄的空間中快速實現(xiàn)位置和姿態(tài)的調(diào)整。李曙光[8]提出一種重量輕、真空驅(qū)動的軟機器人手爪折紙“魔術(shù)球”和柔性薄膜，可以舉起各種各樣的物體。在承載能力方面，LIU等[9]提出了一種基于形狀記憶合金(SMA)的具有可變剛度的軟體機器人抓手，關(guān)節(jié)中的可變剛度材料可以通過加熱或冷卻來改變機器人手指的剛度。張海榮等[10]提出了一個嵌入低熔點合金(LMPA)的變剛度軟體仿人手指，通過改變溫度，使合金在固液相之間轉(zhuǎn)變而進行整體剛度的轉(zhuǎn)變。

綜上以上文獻發(fā)現(xiàn)，同時兼具良好順應(yīng)性和較強承載能力的軟體抓持裝置的研究成果尚不多見，其剛度變化連續(xù)性差和響應(yīng)速度慢的問題限制了軟體抓持裝置的應(yīng)用范圍。

為此本研究提出了一種可變剛度的流體驅(qū)動式模塊化軟體仿生象鼻機構(gòu)。通過PNS—嵌入式氣動網(wǎng)絡(luò)各通道內(nèi)氣體的正負(fù)壓協(xié)同控制，實現(xiàn)關(guān)節(jié)多自由度變形；利用電磁場控制磁流變液流變特性實現(xiàn)剛度連續(xù)調(diào)節(jié)，增強仿生象鼻的承載能力；建立仿生象鼻關(guān)節(jié)變形量和剛度理論數(shù)學(xué)模型,基于Yeoh模型利用有限元方法對仿生象鼻關(guān)節(jié)模型進行變形仿真分析，優(yōu)選出變形量最大的模型結(jié)構(gòu)和最佳輸入氣體壓力值；對關(guān)節(jié)變形后兩端電磁鐵的疊加磁場分布進行了有限元仿真，通過分析磁場分布均勻性初步確定了兩端電磁鐵形成磁場的最佳距離即關(guān)節(jié)長度；最后對不同長度關(guān)節(jié)進行了變形仿真分析，綜合考慮磁場分布均勻性和變形量，通過兩者的最優(yōu)組合確定最終關(guān)節(jié)長度，驗證本研究所設(shè)計的仿生象鼻的多自由度變形和可變剛度特性。預(yù)期效果可拓寬軟體機器人在軍事偵察、災(zāi)難救援以及科學(xué)探測等重要領(lǐng)域的應(yīng)用范圍，具有重要的應(yīng)用價值。

1 仿生象鼻結(jié)構(gòu)設(shè)計

所設(shè)計的仿生象鼻機構(gòu)由柔性控制系統(tǒng)、剛度控制系統(tǒng)和輔助支承系統(tǒng)組成，如圖1所示。

圖1 仿生象鼻整體圖

柔性控制系統(tǒng)包括氣動系統(tǒng)和四個軟體關(guān)節(jié)的PNS通道。關(guān)節(jié)的4個通道可通入不同壓力的正負(fù)壓氣體變化組合，實現(xiàn)多個自由度的變形控制，其部分變形狀態(tài)如表1所示，狀態(tài)1為象鼻向Y軸正方向彎曲；狀態(tài)2為象鼻向Y軸正方向增強彎曲；狀態(tài)3為象鼻向X軸正方向彎曲；狀態(tài)4為象鼻向第一象限方向彎曲，根據(jù)1通道和2通道通入氣壓的大小決定偏向于Y軸正方向還是X軸正方向；狀態(tài)5為象鼻向第二象限方向彎曲；狀態(tài)6為象鼻向第三象限方向彎曲；狀態(tài)7為象鼻向第四象限方向彎曲。

表1 變形情況分析

剛度控制系統(tǒng)由關(guān)節(jié)兩端電磁鐵和填充磁流變液的聚氨酯泡沫墊組成，集中安放在每個關(guān)節(jié)的軸心位置。電磁鐵采用直流電源供電，通過控制回路電流大小來調(diào)節(jié)外部磁場。由于磁流變液是由磁性顆粒和硅油組成的混合物，其物理狀態(tài)可以隨著外界磁場發(fā)生變化。當(dāng)磁場迅速增強時，在數(shù)毫秒內(nèi)磁流變液將從低黏度流體轉(zhuǎn)變成高黏度Bingham流體，呈現(xiàn)半固體狀態(tài)，并具有一定的抗剪切能力；當(dāng)磁場迅速減弱時，其狀態(tài)在數(shù)毫秒內(nèi)會立刻恢復(fù)為一般流體，且這樣的可逆變化是連續(xù)的，且狀態(tài)穩(wěn)定，此效應(yīng)稱之為磁流變液的磁流變效應(yīng)[11]。因此可以通過控制磁場的強弱來改變磁流變液黏度，從而實現(xiàn)關(guān)節(jié)剛度的動態(tài)連續(xù)可調(diào)。

輔助支承系統(tǒng)可實現(xiàn)對象鼻的支承和彎曲變形前后的靈活控制。

注：正壓為“+”，負(fù)壓為“-”，不通為“0”。

2 關(guān)節(jié)性能參數(shù)數(shù)學(xué)模型

2.1 變形量

假設(shè)關(guān)節(jié)變形時中性軸長度不變，且彎曲變形在長度方向上均勻變化，由此建立關(guān)節(jié)結(jié)構(gòu)的理想物理模型，如圖2所示。

圖2 關(guān)節(jié)結(jié)構(gòu)理想物理模型

向關(guān)節(jié)內(nèi)充入一定壓力的氣體，其變形量和彎曲角度的關(guān)系式為：

(1)

由(1)式可得：

(2)

式中，R—— 關(guān)節(jié)的邊緣距離中性軸的距離

L0—— 關(guān)節(jié)的初始長度，m

L1—— 關(guān)節(jié)變形后的長度，m

θ—— 關(guān)節(jié)的彎曲角度，°

R1—— 關(guān)節(jié)的彎曲半徑，m。

2.2 剛度

關(guān)節(jié)變形后需要產(chǎn)生一定的剛度來承載，剛度越大承載能力越強。引入力與運動的變換，建立仿生象鼻關(guān)節(jié)的可變剛度模型。在外力f作用下，關(guān)節(jié)發(fā)生轉(zhuǎn)動量為θ，當(dāng)外力和關(guān)節(jié)轉(zhuǎn)動量足夠小時，二者可以看做近似的線性關(guān)系[12]：

(3)

式中，k—— 關(guān)節(jié)的剛度

機器人的剛度包含多種，如靜剛度、動剛度、伺服剛度和機械剛度。考慮到剛度的調(diào)節(jié)是在象鼻完全抓牢物體后發(fā)生，因此，本研究只分析靜剛度的變化。

對于施加在關(guān)節(jié)上的扭矩τ，當(dāng)扭矩和關(guān)節(jié)位移較小時象鼻關(guān)節(jié)靜剛度將滿足：

(4)

磁流變液在無磁場作用時表現(xiàn)為牛頓流體的一般流動特性，此時產(chǎn)生的阻尼力矩很小；而在外加磁場作用下，磁流變液中的磁性微顆粒相互吸引形成鏈狀或者柱狀，使阻尼力矩增大，從而作用在關(guān)節(jié)的扭矩也增大。

本研究采用的磁流變液磁化性能曲線如圖3所示。可通過電磁場磁感應(yīng)強度分析來反映磁流變液的變剛度特性。

圖3 磁化性能曲線

采用直流電源供電時關(guān)節(jié)變剛度一般工作在磁感應(yīng)穩(wěn)定區(qū)，磁致屈服應(yīng)力σ與外加磁場強度H存在近似的線性關(guān)系：

σ(H)=αH

(5)

根據(jù)B-H曲線進一步地：

(6)

上式表示的是磁流變液的磁致屈服應(yīng)力、外加磁場和磁感應(yīng)強度之間的關(guān)系，磁流變液的磁致屈服應(yīng)力與外加磁場成正相關(guān)。因為磁流變液與關(guān)節(jié)是直接接觸的，所以磁致屈服應(yīng)力近似等于磁流變液作用在關(guān)節(jié)的扭矩，因此有：

(7)

由式(7)可知：關(guān)節(jié)的剛度大小由外加磁場強度大小決定。

3 關(guān)節(jié)的仿真分析

3.1 變形量仿真分析

仿生象鼻關(guān)節(jié)選取硅膠作為軟體材料，采用Yeoh模型建立硅膠材料的本構(gòu)模型，典型的二項參數(shù)形式密度函數(shù)W為：

W=C10(I1-3)+C20(I1-3)2

(8)

式中，I1=λ1+λ2+λ3，λ1、λ2、λ3分別為軸向、徑向和周向主拉伸比；C10、C20為硅膠材料的 Yeoh 模型材料常系數(shù)。

令材料常系數(shù)C10=0.11、C20=0.02，在75 kPa氣體壓力下分析關(guān)鍵結(jié)構(gòu)參數(shù)對關(guān)節(jié)變形量的影響。圖4為關(guān)節(jié)的關(guān)鍵結(jié)構(gòu)參數(shù)示意圖，其中a為內(nèi)壁壁厚，b為外壁與內(nèi)壁通道寬度，c為外壁通道高度，d為外壁厚度，e為外壁通道寬度。圖5為改變內(nèi)壁壁厚a的部分仿真圖，由仿真結(jié)果可知，隨著內(nèi)壁壁厚a的增大，其彎曲角度逐漸減小，由此可確定不同內(nèi)壁壁厚a對應(yīng)的彎曲角度的大小。進而通過式(2)求出相應(yīng)變形量L0-L1。

圖4 關(guān)鍵結(jié)構(gòu)參數(shù)示意圖

圖5 改變內(nèi)壁壁厚a的部分變形仿真

同理可求出改變關(guān)鍵結(jié)構(gòu)參數(shù)b、c、d、e對應(yīng)的變形量L0-L1，其變形量隨以上g變化規(guī)律，如圖6所示。

由圖6可知，變形量隨著內(nèi)壁壁厚a的增大而減小，隨著外壁與內(nèi)壁通道寬度b的增大而增大，隨著外壁通道高度c的增大而增大，隨著外壁壁厚d的增大而減小，隨著外壁通道寬度e的增大而減小，且外壁與內(nèi)壁通道寬度b、外壁通道高度c和外壁厚度d對變形量影響比較大，內(nèi)壁壁厚a和外壁通道寬度e對變形量影響比較小。為確定最優(yōu)關(guān)節(jié)參數(shù)尺寸，還需要考慮3個影響比較大的結(jié)構(gòu)參數(shù)，對變形量的影響，因此，在這3個影響比較大的結(jié)構(gòu)參數(shù)中，每個結(jié)構(gòu)參數(shù)選出3個最優(yōu)數(shù)值；在剩下的2個影響比較小的結(jié)構(gòu)參數(shù)中，選出一個最優(yōu)數(shù)值，使用正交試驗法設(shè)計仿真實驗，用以上同樣的方法進行仿真分析，其結(jié)果如圖7所示。

圖6 改變關(guān)鍵結(jié)構(gòu)參數(shù)與關(guān)節(jié)變形量關(guān)系

圖7 結(jié)構(gòu)參數(shù)組合與變形量關(guān)系

通過27組仿真分析，確定第25組結(jié)構(gòu)參數(shù)組合即a=4,b=6.5,c=14,d=4,e=5為最優(yōu)組合。

用以上確定最優(yōu)關(guān)節(jié)參數(shù)的方法，確定仿生象鼻變形的最優(yōu)氣體壓力。本研究在不同氣體壓力下進行了變形仿真分析，其變形量隨氣體壓力的變化規(guī)律如圖8所示。隨著氣體壓力的增加，變形量逐漸增大，當(dāng)氣體壓力為100 kPa時，單個關(guān)節(jié)變形量達(dá)到151.42 mm，仿生象鼻4個關(guān)節(jié)配合變形可以將直徑為90 mm、容積為500 ml的燒杯卷起，軟體結(jié)構(gòu)模型達(dá)到變形要求，綜合考慮壓力損失等影響因素，取氣動系統(tǒng)額定壓力為150 kPa。

圖8 氣體壓力與變形量關(guān)系

3.2 變剛度仿真分析

本研究選用的的圓盤型電磁鐵結(jié)構(gòu)參數(shù)如表2所示。

表2 電磁鐵結(jié)構(gòu)參數(shù)

通過改變回路電流，來改變電磁鐵產(chǎn)生磁場的大小，在任何工況下都可以等效為變化的靜磁場來進行分析，于是選用靜磁場分析模式對電磁鐵進行模擬。

用磁場大小分布的均勻性來衡量關(guān)節(jié)的變剛度特性[13]，由此來評價仿生象鼻承載能力的提升效果。因此，本研究以變形前關(guān)節(jié)長度L0，即聚氨酯填充的磁流變液兩端電磁鐵距離為變量進行磁場分布的仿真分析。由于該仿生象鼻是發(fā)生完變形之后變剛度控制，因此對變形后的關(guān)節(jié)進行磁場大小分布均勻性分析。

以最優(yōu)模型在100 kPa的變形仿真圖來放置電磁鐵，如圖9所示。

圖9 電磁鐵相對位置

相對位置的具體數(shù)值可由變形前兩端電磁鐵的距離求出，以變形前關(guān)節(jié)長度L0為變量進行磁場分布的仿真分析，部分仿真結(jié)果，如圖10所示。

圖10 磁場分布仿真圖

上圖設(shè)置的變形前關(guān)節(jié)長度分別為L0=20 cm、L0=30 cm、L0=40 cm、L0=50 cm、L0=60 cm。取均勻分布的5個點，查出磁場強度大小，求其方差，得出變形前關(guān)節(jié)長度與磁場強度方差的關(guān)系圖，如圖11所示。用方差ζ的大小來描述磁場分布的均勻性，方差越小均勻性越好。

由圖11可知，L0=40 cm時，其磁場強度方差最小，磁場分布均勻性最好，由2.2節(jié)式(7)可知，關(guān)節(jié)的磁場均勻性越好，關(guān)節(jié)的剛度均勻性越好，因此確定變形前關(guān)節(jié)長度為40厘米。

圖11 變形前關(guān)節(jié)長度與磁場強度方差的關(guān)系

3.3 變形量和變剛度綜合仿真分析

由3.2節(jié)可知，關(guān)節(jié)的最優(yōu)長度為40 cm，但是關(guān)節(jié)的長度也會影響其變形效果，因此，需要綜合考慮剛度和變形效果來確定關(guān)節(jié)的最優(yōu)長度。所以，還需要分析不同關(guān)節(jié)長度L0的變形效果，因此選用與3.1節(jié)同樣的方法求出變形量L0-L1，做出關(guān)節(jié)長度L0與變形量L0-L1關(guān)系圖，如圖12所示。

圖12 變形前關(guān)節(jié)長度與變形量關(guān)系

由上圖12可知，關(guān)節(jié)長度越長，其變形量越大，變形效果越好。

因此，綜合分析磁場分布均勻性和變形效果，確定關(guān)節(jié)的最優(yōu)長度為40 cm，2個電磁鐵間產(chǎn)生的磁場分布較為均勻、變形效果顯著，較好地解決了順應(yīng)性和承載能力不能兼顧和剛度變化不連續(xù)的問題。

4 結(jié)論

本研究對一種新型軟體抓持裝置開展了變形量與變剛度特性的理論和仿真研究，具體得到以下結(jié)論：

(1) 針對軟體抓持裝置順應(yīng)性和承載能力難于兼顧的問題，提出了一種可變剛度的流體驅(qū)動式模塊化軟體仿生象鼻機構(gòu)；

(2) 建立了仿生象鼻關(guān)節(jié)變形量和剛度參數(shù)數(shù)學(xué)模型；

(3) 利用有限元分析方法，基于Yeoh模型對不同結(jié)構(gòu)參數(shù)的仿生象鼻關(guān)節(jié)進行變形量的仿真分析，得到了75 kPa氣體壓力下的關(guān)鍵結(jié)構(gòu)參數(shù)；

(4) 對電磁鐵的磁場分布均勻性進行了仿真分析，初步確定了兩端電磁鐵形成磁場的最佳距離即關(guān)節(jié)長度；

(5) 對不同長度關(guān)節(jié)進行了變形仿真分析，綜合考慮磁場分布均勻性和變形量，通過兩者的最優(yōu)組合確定最終關(guān)節(jié)長度，驗證了該仿生象鼻的多自由度變形和可變剛度特性。