奚如如 王興松 韓亞麗

(1 東南大學機械工程學院，南京211189)

(2 南京工程學院機械工程學院，南京211167)

被動雙足行走指的是雙足行走機構(gòu)只依靠本身結(jié)構(gòu)、重力、足與地面的摩擦力及碰撞力而實現(xiàn)雙足行走［1-2］.被動行走機構(gòu)不包含可以將系統(tǒng)非機械能轉(zhuǎn)換為機械能的部件(如電機等)，但可以包含如彈簧阻尼之類的儲能元件.基于被動行走理論的雙足機器人是欠驅(qū)動的，實現(xiàn)自然狀態(tài)下的行走運動，無附加驅(qū)動力而產(chǎn)生的能量損耗.被動行走研究的主要思路是首先研究無驅(qū)動的純被動機器人［3］.現(xiàn)有簡單類雙足機構(gòu)二維無邊輪輻模型，可在斜坡上滾動，通過間歇性碰撞實現(xiàn)相對持續(xù)的運動.簡單無驅(qū)動被動雙足行走模型［1-2］的支撐腿與擺動腿通過自由鉸鏈連接，支撐腿呈現(xiàn)倒立擺行走特性，擺動腿以自由鉸鏈為節(jié)點自由擺動.輪輻支撐輻條的運動及其與斜坡間的碰撞分別對應(yīng)雙足行走模型支撐腿的運動及其擺動腿與斜面的碰撞，隨著碰撞的發(fā)生，輪輻兩相鄰支撐輻條完成交替，且對應(yīng)于雙足模型支撐腿與擺動腿的交替.輪輻的運動具有純被動特性，呈現(xiàn)了人類行走過程中足與地面的接觸及支撐腿與地面的類似倒立擺的相對運動，對其動力學及相鄰輻條間碰撞交替的研究對被動行走本質(zhì)的探究更具直觀性.

McGeer［1-2］最先基于無邊輪輻的運動特質(zhì)提出被動雙足行走理論，并對無邊輪輻運動及其穩(wěn)定性進行了初步理論分析.文獻［4-9］基于McGeer提出的具有被動行走特性的無邊輪輻模型，對被動雙足行走機器人進行了深入的研究.文獻［5-13］指出無邊輪輻系統(tǒng)是包含連續(xù)動力學與離散動力學的復(fù)合系統(tǒng).文獻［14-19］雖對無邊輪輻運動及其被動特性進行了基本描述和運動仿真，但對動力學特性的數(shù)值分析、輪輻運動極限環(huán)、初始條件及機構(gòu)參數(shù)之間的相互依存關(guān)系沒有進行直接深入的研究，從而限制了對其動力學特性的直接應(yīng)用.

本文對無邊輪輻模型進行動力學特性分析，揭示輪輻實現(xiàn)穩(wěn)定被動行走的可能性、影響因素及其所依據(jù)的數(shù)值關(guān)系.基于倒立擺模型及角動量矩守恒原理建立表征無邊輪輻模型的運動學方程及碰撞瞬間切換方程，利用龐加萊截面及映射理論對斜坡上的無邊輪輻進行運動穩(wěn)定性分析.結(jié)果表明，具有復(fù)合動力學特性的輪輻系統(tǒng)在特定條件下可實現(xiàn)穩(wěn)定持續(xù)的下坡滾動，不動點與極限環(huán)的存在是輪輻運動穩(wěn)定性的主要表征參數(shù).不動點與極限環(huán)的存在與輪輻的輻條數(shù)目、慣性矩參數(shù)、運動初始條件和斜坡角度有關(guān)，驗證了一定的斜坡角度和初始角速度下的輪輻可實現(xiàn)穩(wěn)定持續(xù)的下坡滾動.本文對輪輻運動不動點及極限環(huán)的數(shù)值研究方法具有普遍適用性，可直接指導無驅(qū)動雙足被動行走模型的研究.若將輪輻下坡時重力利用等同于雙足機器人的驅(qū)動力，則對以固有頻率行走、減少驅(qū)動能耗及簡化控制理論為目標的主動雙足行走及下肢外骨骼的步態(tài)規(guī)劃及機構(gòu)研制具有借鑒意義.

1 無邊輪輻系統(tǒng)模型

1.1 系統(tǒng)模型及運動參數(shù)配置

如圖1所示，無邊輪輻模型成軸對稱分布，質(zhì)量為m，質(zhì)心點C，轉(zhuǎn)動慣量為IC，輻條分布均勻，輻條總數(shù)為n，輻條長度為l，相鄰輻條之間的角度為β=2π/n，斜坡角度α.滾動過程中當前接觸地面的輻條定義為支撐輻條，輻條角度θ 定義為輻條相對于地面垂直線的角度，逆時針方向為正.

圖1 無邊輪輻系統(tǒng)模型

為了更好地研究無邊輪輻運動特性，作如下假設(shè):①當前支撐輻條端點與坡面的接觸點相當于鉸鏈鏈接點;②輻條與坡面間的碰撞為非彈性且瞬時完成，無滑移和彈起現(xiàn)象;③輻條與坡面碰撞瞬間無沖擊，兩相鄰輻條分別平穩(wěn)地接觸和離開斜坡;④輻條端點與坡面接觸點之間的摩擦力損失忽略不計［5-7］.

鑒于上坡和下坡現(xiàn)象在無邊輪輻運動過程中可能交替出現(xiàn)，需對輪輻運動參數(shù)標示明確.

1)為避免模運算，輻條表征參數(shù)k 采用逐一計數(shù)模式，k=0 定義為無邊輪輻與地面的初始狀態(tài).無邊輪輻從k=0 輻條開始運動，逆時針運動時首個支撐輻條記為k=1，順時針運動時首個支撐輻條記為k= －1，依次加減1 類推.

2)i 表征碰撞次數(shù)，i=0 為初始狀態(tài)，i 采用逐一計數(shù)方法，每發(fā)生一次碰撞，i 加1，0≤i≤+∞.j(i)表征第i 次碰撞后的支撐輻條，支撐輻條隨著每次碰撞的發(fā)生而改變.輪輻歷經(jīng)i+1 次碰撞后，若其為逆時針下坡運動，支撐輻條j(i +1)=j(i)+1;反之，j(i+1)=j(i)－1.

3)上標(－)和(+)表征碰撞前后瞬間.記ti為第i 次碰撞發(fā)生時刻，t+i為第i 次碰撞后瞬間，為第i 次碰撞前瞬間，θk(t)為輻條k 在t 時刻的角度，θk(ti)為輻條k 在ti時刻的角度.

4)鑒于無邊輪輻與雙足行走特性的研究相關(guān)性，僅考慮當前支撐輻條而不是所有輻條的運動狀態(tài)，輪輻當前支撐輻條運動特性對應(yīng)于雙足行走機器人的支撐腿.記(t)為輪輻第i 次碰撞后支撐輻條j(i)在任意時刻t 的角度，為第i 次碰撞后瞬間支撐輻條j(i)的角度，其與意義相同.類似地，

1.2 模型運動周期及分析

如圖2所示，無邊輪輻在角度為α 的斜坡上作下坡運動，支撐輻條以逆時針方向滾動.如圖2(a)、(b)、(e)、(f)所示，碰撞發(fā)生后瞬間，無邊輪輻分別以(iθ+，i+)與(i+1θ+，i+1+)作為初始條件進入無碰撞階段，運動至下一次碰撞發(fā)生.如圖2(d)所示，輪輻支撐輻條呈豎直狀態(tài)，即支撐輻條角度為零，整個輪輻機構(gòu)關(guān)于支撐輻條對稱.如圖2(c)、(d)、(e)所示，無邊輪輻處于無碰撞階段，整個輪輻的運動可看作是以角度iθ+=α －π/n、角速度i+作為初始條件的非完整周期的倒立擺運動，終止于邊界條件i+1θ－=α +π/n.由圖2(a)、(b)、(c)可知，無邊輪輻整體參數(shù)配置在碰撞前后無顯著變化，但支撐輻條在碰撞前后完成j(i)與j(i +1)兩相鄰輻條的交替，支撐輻條角度發(fā)生幅度為±2π/n 的躍變，且因碰撞引起輪輻整體動能損耗，導致碰撞前后輪輻角速度變化很大.顯然，碰撞前后的支撐輻條角度與角速度均不連續(xù).

當無邊輪輻作上坡運動時，支撐輻條沿順時針方向滾動，可從左向右依次逆解讀圖2中輪輻運動狀態(tài)來說明.圖2(b)、(f)為碰撞瞬間，碰撞點分別為點A 和點D，圖2(a)、(e)為碰撞前瞬間，圖2(c)為碰撞后瞬間，圖2(d)為支撐輻條呈豎直狀態(tài).圖2(c)、(d)、(e)中的輪輻處于無碰撞階段，其參數(shù)配置方式的遵循原則與下坡運動時保持一致.

無邊輪輻無論處于何種初始運動趨勢，在運動過程中其運動趨勢均可能出現(xiàn)翻轉(zhuǎn).輪輻到達如圖2(d)支撐角度為零狀態(tài)時，角速度將決定輪輻是否可以通過支撐輻條來達到豎直狀態(tài)，從而影響輪輻是否繼續(xù)初始運動或出現(xiàn)搖擺運動.

因此，本文定義無邊輪輻一個完整的運動周期為一個完整的無碰撞運動過程和接下來的一次碰撞運動，并在不改變模型參數(shù)意義的基礎(chǔ)上簡化角度及角速度的符號標示，以當前支撐輻條角度θ(t)、角速度(t)作為研究對象，研究無邊輪輻的運動特性.一個完整的輪輻運動周期可表現(xiàn)為上坡運動、下坡運動或組合運動，即在一個運動周期內(nèi)，輪輻運動趨勢可能發(fā)生改變，或出現(xiàn)搖擺現(xiàn)象.

圖2 無邊輪輻模型運動狀態(tài)示意簡圖

2 系統(tǒng)動力學分析

2.1 運動方程

如圖2(c)、(d)、(e)所示，無論輪輻初始運動趨勢為上坡或下坡運動，輪輻都處于無碰撞運動階段，整個輪輻關(guān)于支撐點A 的角動量矩平衡，即

式中，MA為輪輻關(guān)于點A 的所受力矩總和;˙HA為輪輻關(guān)于點A 的角動量變化率.

由式(1)可得無碰撞階段的無量綱運動方程為

2.2 碰撞瞬間切換方程

輪輻碰撞前后瞬間參數(shù)配置如圖2(a)～(c)所示，支撐輻條角度在碰撞前后有±2π/n 的躍變.無邊輪輻在碰撞發(fā)生前為逆時針下坡運動，輪輻碰撞后瞬間相對于碰撞前瞬間轉(zhuǎn)過了－2π/n 角度;反之，若輪輻在碰撞發(fā)生前呈順時針上坡運動趨勢，輪輻則轉(zhuǎn)過2π/n 角度.第i 次碰撞前后瞬間角度躍變公式為

無邊輪輻在碰撞發(fā)生前為逆時針下坡運動，iθ－=α+π/n，iθ+=α －π/n;反之，若無邊輪輻在碰撞發(fā)生前為順時針上坡運動，則iθ－=α －π/n，iθ+=α+ π/n.無論何種運動趨勢下發(fā)生碰撞運動，碰撞前后含符號函數(shù)的角度躍變公式(3)均可改寫為如下簡單形式:

無邊輪輻的角速度隨著碰撞的發(fā)生而改變，忽略重力及碰撞瞬間對輪輻引起的沖擊，則整個輪輻關(guān)于點B 角動量守恒，即

式中，H－B 為第i 次碰撞前瞬間輪輻關(guān)于點B 的角動量;H+B為第i 次碰撞后瞬間輪輻關(guān)于點B 的角動量.

由式(5)可得第i 次碰撞前后瞬間角速度躍變方程為

式中，碰撞參數(shù)λ=1 + μ2(cos(2π/n)－1)，0 ≤λ ＜1.

第i 次碰撞前后瞬間角度、角速度躍變公式(4)、(6)表征了碰撞發(fā)生前后輪輻支撐輻條運動參數(shù)的改變，式(4)、(6)可統(tǒng)稱為碰撞瞬間切換模型.

2.3 關(guān)鍵角速度臨界值

引進臨界角速度參數(shù)dn和up，獲取依據(jù)為無邊輪輻介于2 次碰撞之間的運動遵循能量守恒原則.在一個運動周期的無碰撞階段，輪輻以一定的初始角速度運動，輪輻的初始動能可根據(jù)圖2(c)、(d)、(e)中的動能和勢能表征條轉(zhuǎn)換.在θ=0 處完全轉(zhuǎn)化為勢能，即在θ=0 處，輪輻角速度亦為零，此時輪輻處于基于1 個支撐輻條的不穩(wěn)定豎直平衡狀態(tài).

當無邊輪輻呈下坡運動趨勢時，輪輻以臨界角速度dn作初始運動，可根據(jù)圖2(d)與(e)間的能量表征條進行能量轉(zhuǎn)換;反之，當輪輻呈上坡運動趨勢時，輪輻以臨界角速度up作為初始速度運動，根據(jù)圖2(c)與(d)間能量表征條進行能量轉(zhuǎn)換，計算可得

無邊輪輻呈下坡運動趨勢且碰撞后瞬間角速度為dn或上坡運動趨勢且碰撞后瞬間角速度為時，輪輻將耗時無限長到達并停留在位置θ=0處，輪輻基于1 個支撐輻條豎直靜止，是一種臨界不穩(wěn)定平衡狀態(tài).

2.4 相圖及運動軌跡流

無碰撞階段的運動方程式(2)可改寫成如下一階微分方程組:

積分后可得介于2 次碰撞之間的相空間軌跡流，即

鑒于支撐輻條與斜坡之間的空間限制，支撐輻條并不能進行完整的倒立擺運動，因此無邊輪輻的軌跡變化需在相空間U 中進行研究，即

無碰撞階段運動方程式(8)及碰撞瞬間切換模型(4)、(6)，均可表達為如下形式:

2)圖3中相空間U 內(nèi)的連續(xù)運動軌跡流，根據(jù)獨立向量場v 發(fā)生連續(xù)變化，表明輪輻在無碰撞階段的動力學特性.相空間邊界?U 上的碰撞切換軌跡流，依據(jù)獨立向量場h 發(fā)生離散性質(zhì)的跳躍變化，表明輪輻碰撞發(fā)生前后支撐輻條角度與角速度的躍變.一個完整的運動周期軌跡流為連續(xù)運動軌跡流與相應(yīng)的碰撞切換軌跡流的組合，具有連續(xù)動力學與離散動力學并存的復(fù)合動力學特性.

3)以所屬相空間邊界?U 內(nèi)含臨界值的邊界點x1(α－π/n，dn)與x2(α +π/n，up)為初始條件的臨界運動軌跡流終止于不穩(wěn)定平衡位置θ=0處，輪輻不再出現(xiàn)后續(xù)運動軌跡流.當輪輻連續(xù)滾動時，輪輻上一次碰撞后瞬間(θ，)即為輪輻下一個運動周期的初始條件.若輪輻歷經(jīng)數(shù)次碰撞后，其運動參數(shù)達到臨界邊界點，則輪輻將到達并停留在θ=0 位置，輪輻將不會發(fā)生下一次碰撞.

4)如果斜坡角度α ＞n/π，輪輻角度θ 始終大于零，則相空間內(nèi)不包含不穩(wěn)定豎直平衡狀態(tài)θ=0.

圖3 兩次碰撞間的相空間軌跡流圖

以滿足式(13)的邊界點為初始運動條件，圖3中連續(xù)運動軌跡流表明，輪輻均未通過θ=0 位置，而是逆向終止于相空間邊界?U 上，輪輻出現(xiàn)運動趨勢發(fā)生逆轉(zhuǎn)的搖擺運動.反之，若(θ，)滿足

以滿足式(14)的邊界點為初始運動條件，圖3中相空間U 內(nèi)連續(xù)運動軌跡流均都正向終止于相空間邊界?U 上，輪輻運動趨勢保持不變.

6)無論輪輻以何種初始條件運動，歷經(jīng)數(shù)次碰撞、上坡、搖擺及下坡運動后，輪輻可能處于如下3 種狀態(tài):①基于單個支撐輻條上的豎直狀態(tài)，即運動軌跡流終止于不穩(wěn)定平衡位置;②基于兩相鄰支撐輻條上的穩(wěn)定靜止狀態(tài)，即運動軌跡流終止于相空間邊界?U 內(nèi)的=0 邊界點上;③輪輻以一定速度作持續(xù)下坡滾動.如圖4所示，設(shè)n=6，α=π/16，μ=0.75，輪輻以θ0=α －n/π，0=1 為初始條件開始運動，歷經(jīng)足夠多次碰撞后，連續(xù)運動軌跡流與碰撞切換軌跡流逐漸成為閉合曲線，說明極限環(huán)在特定條件下是存在的，輪輻可以實現(xiàn)穩(wěn)定持續(xù)的下坡滾動.

圖4 無邊輪輻相空間軌跡流及極限環(huán)

3 系統(tǒng)不動點及穩(wěn)定性分析

3.1 龐加萊截面

利用龐加萊回歸映射方法研究輪輻系統(tǒng)不動點、極限環(huán)的影響因素及存在條件，從而判定輪輻是否可實現(xiàn)穩(wěn)定持續(xù)的下坡滾動［5-7］.龐加萊截面Γ 定義為

3.2 龐加萊回歸映射

建立如下龐加萊映射關(guān)系:

式中，iq+=(iθ+，i+)∈Γ，映射函數(shù)f:Γ→Γ，f=(h ?v)是基于離散向量場h 和連續(xù)向量場v 的復(fù)合映射，向量場h 與v 有如下關(guān)系:

映射函數(shù)P 表示輪輻第i 次碰撞后與第i +1次碰撞后的角速度之間的關(guān)系.如式(19)所示，當0≤α≤π/n 時，函數(shù)P 是受臨界值dn，up和參數(shù)集合ψ={μ2，n，α}影響的分段連續(xù)函數(shù)，從而輪輻在運動過程中會出現(xiàn)交替上坡運動、下坡運動及包含2 種運動的搖擺運動;如式(20)所示，當π/n＜α≤π/2 時，函數(shù)P 是受參數(shù)集合ψ={μ2，n，α}影響的連續(xù)線性函數(shù)，若將無邊輪輻兩相鄰輻條支撐于斜坡上，即使給輪輻零初始角速度或上坡運動趨勢，輪輻最終仍會自主持續(xù)下坡滾動.

3.3 映射函數(shù)的不動點

映射函數(shù)P 表征了輪輻相鄰2 次碰撞后瞬間角速度之間的關(guān)系，函數(shù)P 的不動點的存在意味著輪輻歷經(jīng)數(shù)次碰撞后，其后每次碰撞后的瞬間角速度值均相同，這表明輪輻運動存在極限環(huán)，由式(19)、(20)可得函數(shù)P 的不動點ω*.

式中，ω0*為靜止不動點，表明輪輻最終可實現(xiàn)基于兩相鄰輻條支撐的穩(wěn)定靜止狀態(tài);為極限環(huán)不動點，表明極限環(huán)的存在及輪輻最終可實現(xiàn)持續(xù)穩(wěn)定的下坡滾動，即輪輻歷經(jīng)數(shù)次碰撞后角速度值趨向且最終等于，其后每次碰撞后的瞬間角速度值均等于ω*1.

3.4 不動點存在的影響因素

由不動點表達式(21)、(22)可知，不動點的存在與輪輻自身模型參數(shù)μ、輪輻數(shù)目n 及斜坡角度α 有關(guān).當輪輻模型參數(shù)一定時，函數(shù)P 的2 個不動點ω0*與ω*1與斜坡角度的關(guān)系為

式中，αe為臨界斜坡角度，是決定極限環(huán)不動點ω1*存在與否的臨界值，由式(21)可推出αe滿足

由圖5可知，不同的斜坡角度對應(yīng)不同的一維回轉(zhuǎn)圖，不動點性質(zhì)及存在個數(shù)也不同.

如圖5(a)所示，當0≤α≤αe時，輪輻存在靜止不動點ω*0及表征不穩(wěn)定平衡位置的臨界值up和dn，不存在極限環(huán)不動點ω*1 ，輪輻最終處于穩(wěn)定靜止狀態(tài)或不穩(wěn)定平衡狀態(tài).

如圖5(b)所示，當αe＜α≤π/n 時，輪輻存在靜止不動點ω0*、極限環(huán)不動點ω1*及臨界值和等.輪輻以不同的初始條件運動，歷經(jīng)數(shù)次碰條件取決于角速度臨界值撞，輪輻角速度逐漸趨向于ω0*，ω1*，up或dn等4種角速度值，此時輪輻最終運動狀態(tài)可能有3 種:①輪輻基于兩相鄰輻條搖擺數(shù)次，角速度在有限時間內(nèi)衰減至靜止不動點ω*0，輪輻最終靜止于兩相鄰輻條上;②輪輻碰撞后瞬間角速度值以增速和減速2 種方式逐漸趨向且最終等于極限環(huán)不動點ω1*，輪輻實現(xiàn)以ω1*為運動角速度的穩(wěn)定持續(xù)的下坡滾動;③碰撞后瞬間輪輻角速度達到臨界值up和dn，然后輪輻將耗時無限長到達并停留在豎直位置，即不穩(wěn)定平衡狀態(tài)θ=0 位置.

圖5 不同斜坡角度對應(yīng)的不動點與相應(yīng)的一維回轉(zhuǎn)圖(n=6，μ=0.75)

3.5 不動點穩(wěn)定性判據(jù)

不動點的穩(wěn)定性由映射函數(shù)在不動點處的一階導數(shù)決定.函數(shù)P 在不動點ω*0與ω*1處的一階導數(shù)為

由于0≤λ ＜1，映射函數(shù)P 在不動點ω*0與ω1*處的一階導數(shù)DP(ω*)滿足因此，映射函數(shù)P 的不動點ω*0與ω*1是漸進穩(wěn)定的.

3.6 不動點的吸引盆

當α=π/15，n=6，μ=0.75，回轉(zhuǎn)圖如圖6所示.圖6表示了輪輻靜止不動點ω*0與極限環(huán)不動點ω*1及臨界值dn和up等均存在時的多種運動狀態(tài).靜止不動點ω*0與極限環(huán)不動點ω*1交替吸引相應(yīng)區(qū)域內(nèi)的初始條件，從而決定輪輻運動狀態(tài).吸引盆在回轉(zhuǎn)圖中表現(xiàn)為不同區(qū)域內(nèi)的線段，吸引盆的邊界由角速度臨界值來界定.圖6中奇數(shù)值區(qū)域內(nèi)線段為極限環(huán)不動點的吸引盆，偶數(shù)值區(qū)域內(nèi)線段為靜止不動點的吸引盆.由此可知，滿足極限環(huán)存在的初始所界定的區(qū)域，且不動點ω*0與ω*1相對應(yīng)的吸引盆交替出現(xiàn).

圖6 臨界角速度值與相應(yīng)不動點的回轉(zhuǎn)圖

當輪輻以式(26)中臨界值系列為初始條件開始運動時，其終止于θ=0 位置而無法到達任一不動點，由映射函數(shù)P 可推導出m≥1 時的系列臨界值，其表達式為

綜上所述，輪輻最終所處運動狀態(tài)依賴于斜坡角度及由臨界角速度值界定的不動點ω0*與ω1*吸引盆內(nèi)的初始角速度.斜坡角度、初始角速度與不動點類型的具體依存關(guān)系如表1所示，輪輻最終運動狀態(tài)與初始運動趨勢如圖7所示.①若輪輻初始角速度或歷經(jīng)數(shù)次碰撞后的角速度為臨界角速度值，無論斜坡角度如何，輪輻終將處于不穩(wěn)定豎直狀態(tài);②當斜坡角度為0≤α≤αe時，初始角速度為非臨界角速度值，所有角速度值都位于靜止不動點ω*0的吸引盆內(nèi)，即無論輪輻為何種初始運動趨勢，輪輻終將歷經(jīng)M 次搖擺后基于兩相鄰支撐輻條處于穩(wěn)定靜止狀態(tài);③當斜坡角度滿足αe＜α≤π/n，且角速度滿足或則所有角速度值都位于靜止不動點ω*0的吸引盆內(nèi)，即輪輻歷經(jīng)M 次搖擺后將基于兩相鄰支撐輻條處于穩(wěn)定靜止狀態(tài);④當斜坡角度滿足αe＜α≤π/n 且初始角速度滿足則角速度值均位于極限環(huán)不動點ω1*的吸引盆內(nèi)，輪輻由最初的上坡滾動通過一次搖擺逆轉(zhuǎn)為下坡滾動，歷經(jīng)數(shù)次碰撞，輪輻角速度最終增速趨向并等于極限環(huán)不動點ω*1，輪輻呈穩(wěn)定持續(xù)下坡滾動狀態(tài);⑤當斜坡角度與初始角速度滿足關(guān)系αe＜α≤π/n，＞dn或π/n ＜α≤π/2，≥0 時，輪輻會以增速或減速2 種方式趨向且最終等于極限環(huán)不動點ω1*，輪輻實現(xiàn)純被動性質(zhì)的穩(wěn)定持續(xù)下坡滾動.

表1 不動點與相應(yīng)吸引盆及斜坡角度的依存關(guān)系

圖7 輪輻初始條件與運動過程及最終運動狀態(tài)關(guān)系圖

4 結(jié)論

1)輪輻模型參數(shù)一定時，不同的斜坡角度對應(yīng)不同的龐加萊映射回轉(zhuǎn)圖及不動點的存在狀態(tài)，是輪輻運動極限環(huán)存在的重要影響因素.決定極限環(huán)不動點存在的臨界斜坡角度值依賴于輪輻輻條數(shù)目及慣性矩參數(shù).

2)輪輻自身模型參數(shù)及所在斜坡角度一定時，靜止不動點與極限環(huán)不動點的吸引盆區(qū)域確定，輪輻初始運動條件所屬吸引盆類別直接決定輪輻最終運動狀態(tài).

3)臨界角速度系列為不動點的吸引盆邊界值，由輪輻自身模型參數(shù)及斜坡角度決定.歷經(jīng)數(shù)次碰撞后，以此系列臨界角速度為初始條件開始運動的輪輻將耗時無限長，趨向并停留于不穩(wěn)定豎直平衡位置，而無法到達靜止不動點或極限環(huán)不動點.

4)輪輻只在重力作用下發(fā)生滾動，無附加驅(qū)動力，具有純被動性質(zhì)，輪輻自身模型參數(shù)、斜坡角度、初始運動條件將決定輪輻運動過程及最終運動狀態(tài).

5)輪輻自身模型參數(shù)即輻條數(shù)目、慣性矩參數(shù)與臨界斜坡角度、臨界角速度值、不動點、不動點的吸引盆相互依賴并相互制約.

無邊輪輻輻條數(shù)目與慣性矩參數(shù)類比對應(yīng)于簡單被動雙足行走模型運動過程中雙腿間最大角度、腿部與髖關(guān)節(jié)處質(zhì)量分布，斜坡角度決定無邊輪輻或被動雙足行走模型行走過程中重力驅(qū)動力.因此，兩相鄰輻條間角度的表征參數(shù)輻條數(shù)目、慣性矩參數(shù)與斜坡角度是無邊輪輻、被動雙足行走模型的穩(wěn)定周期狀態(tài)的決定性因素.

References)

［1］ McGeer T.Passive dynamic walking［J］.The International Journal of Robotics Research，1990，9(2):62-82.

［2］ McGeer T.Passive dynamic biped catalogue ［C］//Proceedings of the Experimental Robotics Ⅱ:the 2nd International Symposium.Toulouse，F(xiàn)rance，1991，190:463-490.

［3］ 毛勇，王家廞，賈培發(fā)，等.雙足被動步行研究綜述［J］.機器人，2007，29(3):274-280.Mao Yong，Wang Jiaxin，Jia Peifa，et al.Passive dynamic biped walking:a survey ［J］.Robot，2007，29(3):274-280.(in Chinese)

［4］ Collins S，Ruina A，Tedrake R，et al.Efficient bipedal robots based on passive-dynamic walkers［J］.Science，2005，307(5712):1082-1085.

［5］ Coleman M J，Chatterjee A，Ruina A.Motions of a rimless spoked wheel:a simple three-dimensional system with impacts［J］.Dynamics and Stability of Systems，1997，12(3):139-159.

［6］ Coleman M J.A stability study of a three-dimensional passive-dynamic model of human gait ［D］.Ithaca，NY，USA:Cornell University，1998.

［7］ Coleman M J.Dynamics and stability of a rimless spoked wheel:a simple 2D system with impacts［J］.Dynamical Systems an International Journal ，2010，25(2):215-238.

［8］ Hodgins J K.Bipedal gait transitions ［C］//Proceedings of the IEEE International Conference on Robotics and Automation.Sacramento，CA，USA，1991:2092-2097.

［9］ Anderson S O，Wisse M，Atkeson C G，et al.Powered bipeds based on passive dynamic principles［C］//Proceedings of the 5th IEEE-RAS International Conference on Humanoid Robots.Tsukuba，Japan，2005:110-116.

［10］ Asano F，Kawamoto J.Modeling and analysis of passive viscoelastic-legged rimless wheel that generates measurable period of double-limb support［J］.Multibody System Dynamics，2014，31(2):111-126.

［11］ Beigzadeh B，Meghdari A，Sohrabpour S.Passive dynamic object manipulation:a framework for passive walking systems［J］.Journal of Multi-Body Dynamics，2013，227(2):185-198.

［12］ Sabaaapour M R，Hairi-Yazdi M R，Beigzadeh B.Towards passive turning in biped walkers［J］.Procedia Technology，2014，12(1):98-104.

［13］ Gamus B，Or Y.Dynamic bipedal walking under stick-slip transitions ［J］.Society for Industrial and Applied Mathematics:Journal of Applied Dynamical Systems，2015，14(2):609-642.

［14］ 劉振澤.欠驅(qū)動步行機器人運動學機理與控制策略研究［D］.長春:吉林大學通信工程學院，2007.

［15］ 劉振澤，周長久，田彥濤.Compass-like 無動力雙足行走機器人的運動狀態(tài)［J］.吉林大學學報:工學版，2007，37(5):1175-1180.Liu Zhenze，Zhou Changjiu，Tian Yantao.Analysis of compass-like passive biped robot locomotion ［J］.Journal of Jilin University:Engineering and Technology Edition，2007，37(5):1175-1180.(in Chinese)

［16］ 劉麗梅.被動動態(tài)行走雙足機器人的穩(wěn)定性分析與控制研究［D］.長春:吉林大學通信工程學院，2011.

［17］ 胡運富.無膝雙足被動機器人的運動特性和穩(wěn)定性研究［D］.哈爾濱:哈爾濱工業(yè)大學機電工程學院，2009.

［18］ 郭欣.雙足被動機器人行走模式的設(shè)計構(gòu)想［D］.北京:北京郵電大學自動化學院，2009.

［19］ 胡運富，朱延河，吳曉光，等.被動無邊輪輻運動特性的Adams 仿真分析［J］.哈爾濱工業(yè)大學學報，2010，42 (7):1076-1079.Hu Yunfu，Zhu Yanhe，Wu Xiaoguang，et al.Adams simulation analysis of motion characteristic of passive rimless-wheel ［J］.Journal of Harbin Institute of Technology，2010，42(7):1076-1079.(in Chinese)