奚如如 王興松 韓亞麗
(1 東南大學(xué)機(jī)械工程學(xué)院,南京211189)
(2 南京工程學(xué)院機(jī)械工程學(xué)院,南京211167)
被動雙足行走指的是雙足行走機(jī)構(gòu)只依靠本身結(jié)構(gòu)、重力、足與地面的摩擦力及碰撞力而實(shí)現(xiàn)雙足行走[1-2].被動行走機(jī)構(gòu)不包含可以將系統(tǒng)非機(jī)械能轉(zhuǎn)換為機(jī)械能的部件(如電機(jī)等),但可以包含如彈簧阻尼之類的儲能元件.基于被動行走理論的雙足機(jī)器人是欠驅(qū)動的,實(shí)現(xiàn)自然狀態(tài)下的行走運(yùn)動,無附加驅(qū)動力而產(chǎn)生的能量損耗.被動行走研究的主要思路是首先研究無驅(qū)動的純被動機(jī)器人[3].現(xiàn)有簡單類雙足機(jī)構(gòu)二維無邊輪輻模型,可在斜坡上滾動,通過間歇性碰撞實(shí)現(xiàn)相對持續(xù)的運(yùn)動.簡單無驅(qū)動被動雙足行走模型[1-2]的支撐腿與擺動腿通過自由鉸鏈連接,支撐腿呈現(xiàn)倒立擺行走特性,擺動腿以自由鉸鏈為節(jié)點(diǎn)自由擺動.輪輻支撐輻條的運(yùn)動及其與斜坡間的碰撞分別對應(yīng)雙足行走模型支撐腿的運(yùn)動及其擺動腿與斜面的碰撞,隨著碰撞的發(fā)生,輪輻兩相鄰支撐輻條完成交替,且對應(yīng)于雙足模型支撐腿與擺動腿的交替.輪輻的運(yùn)動具有純被動特性,呈現(xiàn)了人類行走過程中足與地面的接觸及支撐腿與地面的類似倒立擺的相對運(yùn)動,對其動力學(xué)及相鄰輻條間碰撞交替的研究對被動行走本質(zhì)的探究更具直觀性.
McGeer[1-2]最先基于無邊輪輻的運(yùn)動特質(zhì)提出被動雙足行走理論,并對無邊輪輻運(yùn)動及其穩(wěn)定性進(jìn)行了初步理論分析.文獻(xiàn)[4-9]基于McGeer提出的具有被動行走特性的無邊輪輻模型,對被動雙足行走機(jī)器人進(jìn)行了深入的研究.文獻(xiàn)[5-13]指出無邊輪輻系統(tǒng)是包含連續(xù)動力學(xué)與離散動力學(xué)的復(fù)合系統(tǒng).文獻(xiàn)[14-19]雖對無邊輪輻運(yùn)動及其被動特性進(jìn)行了基本描述和運(yùn)動仿真,但對動力學(xué)特性的數(shù)值分析、輪輻運(yùn)動極限環(huán)、初始條件及機(jī)構(gòu)參數(shù)之間的相互依存關(guān)系沒有進(jìn)行直接深入的研究,從而限制了對其動力學(xué)特性的直接應(yīng)用.
本文對無邊輪輻模型進(jìn)行動力學(xué)特性分析,揭示輪輻實(shí)現(xiàn)穩(wěn)定被動行走的可能性、影響因素及其所依據(jù)的數(shù)值關(guān)系.基于倒立擺模型及角動量矩守恒原理建立表征無邊輪輻模型的運(yùn)動學(xué)方程及碰撞瞬間切換方程,利用龐加萊截面及映射理論對斜坡上的無邊輪輻進(jìn)行運(yùn)動穩(wěn)定性分析.結(jié)果表明,具有復(fù)合動力學(xué)特性的輪輻系統(tǒng)在特定條件下可實(shí)現(xiàn)穩(wěn)定持續(xù)的下坡滾動,不動點(diǎn)與極限環(huán)的存在是輪輻運(yùn)動穩(wěn)定性的主要表征參數(shù).不動點(diǎn)與極限環(huán)的存在與輪輻的輻條數(shù)目、慣性矩參數(shù)、運(yùn)動初始條件和斜坡角度有關(guān),驗(yàn)證了一定的斜坡角度和初始角速度下的輪輻可實(shí)現(xiàn)穩(wěn)定持續(xù)的下坡滾動.本文對輪輻運(yùn)動不動點(diǎn)及極限環(huán)的數(shù)值研究方法具有普遍適用性,可直接指導(dǎo)無驅(qū)動雙足被動行走模型的研究.若將輪輻下坡時(shí)重力利用等同于雙足機(jī)器人的驅(qū)動力,則對以固有頻率行走、減少驅(qū)動能耗及簡化控制理論為目標(biāo)的主動雙足行走及下肢外骨骼的步態(tài)規(guī)劃及機(jī)構(gòu)研制具有借鑒意義.
如圖1所示,無邊輪輻模型成軸對稱分布,質(zhì)量為m,質(zhì)心點(diǎn)C,轉(zhuǎn)動慣量為IC,輻條分布均勻,輻條總數(shù)為n,輻條長度為l,相鄰輻條之間的角度為β=2π/n,斜坡角度α.滾動過程中當(dāng)前接觸地面的輻條定義為支撐輻條,輻條角度θ 定義為輻條相對于地面垂直線的角度,逆時(shí)針方向?yàn)檎?
圖1 無邊輪輻系統(tǒng)模型
為了更好地研究無邊輪輻運(yùn)動特性,作如下假設(shè):①當(dāng)前支撐輻條端點(diǎn)與坡面的接觸點(diǎn)相當(dāng)于鉸鏈鏈接點(diǎn);②輻條與坡面間的碰撞為非彈性且瞬時(shí)完成,無滑移和彈起現(xiàn)象;③輻條與坡面碰撞瞬間無沖擊,兩相鄰輻條分別平穩(wěn)地接觸和離開斜坡;④輻條端點(diǎn)與坡面接觸點(diǎn)之間的摩擦力損失忽略不計(jì)[5-7].
鑒于上坡和下坡現(xiàn)象在無邊輪輻運(yùn)動過程中可能交替出現(xiàn),需對輪輻運(yùn)動參數(shù)標(biāo)示明確.
1)為避免模運(yùn)算,輻條表征參數(shù)k 采用逐一計(jì)數(shù)模式,k=0 定義為無邊輪輻與地面的初始狀態(tài).無邊輪輻從k=0 輻條開始運(yùn)動,逆時(shí)針運(yùn)動時(shí)首個(gè)支撐輻條記為k=1,順時(shí)針運(yùn)動時(shí)首個(gè)支撐輻條記為k= -1,依次加減1 類推.
2)i 表征碰撞次數(shù),i=0 為初始狀態(tài),i 采用逐一計(jì)數(shù)方法,每發(fā)生一次碰撞,i 加1,0≤i≤+∞.j(i)表征第i 次碰撞后的支撐輻條,支撐輻條隨著每次碰撞的發(fā)生而改變.輪輻歷經(jīng)i+1 次碰撞后,若其為逆時(shí)針下坡運(yùn)動,支撐輻條j(i +1)=j(i)+1;反之,j(i+1)=j(i)-1.
3)上標(biāo)(-)和(+)表征碰撞前后瞬間.記ti為第i 次碰撞發(fā)生時(shí)刻,t+i為第i 次碰撞后瞬間,為第i 次碰撞前瞬間,θk(t)為輻條k 在t 時(shí)刻的角度,θk(ti)為輻條k 在ti時(shí)刻的角度.
4)鑒于無邊輪輻與雙足行走特性的研究相關(guān)性,僅考慮當(dāng)前支撐輻條而不是所有輻條的運(yùn)動狀態(tài),輪輻當(dāng)前支撐輻條運(yùn)動特性對應(yīng)于雙足行走機(jī)器人的支撐腿.記(t)為輪輻第i 次碰撞后支撐輻條j(i)在任意時(shí)刻t 的角度,為第i 次碰撞后瞬間支撐輻條j(i)的角度,其與意義相同.類似地,
如圖2所示,無邊輪輻在角度為α 的斜坡上作下坡運(yùn)動,支撐輻條以逆時(shí)針方向滾動.如圖2(a)、(b)、(e)、(f)所示,碰撞發(fā)生后瞬間,無邊輪輻分別以(iθ+,i+)與(i+1θ+,i+1+)作為初始條件進(jìn)入無碰撞階段,運(yùn)動至下一次碰撞發(fā)生.如圖2(d)所示,輪輻支撐輻條呈豎直狀態(tài),即支撐輻條角度為零,整個(gè)輪輻機(jī)構(gòu)關(guān)于支撐輻條對稱.如圖2(c)、(d)、(e)所示,無邊輪輻處于無碰撞階段,整個(gè)輪輻的運(yùn)動可看作是以角度iθ+=α -π/n、角速度i+作為初始條件的非完整周期的倒立擺運(yùn)動,終止于邊界條件i+1θ-=α +π/n.由圖2(a)、(b)、(c)可知,無邊輪輻整體參數(shù)配置在碰撞前后無顯著變化,但支撐輻條在碰撞前后完成j(i)與j(i +1)兩相鄰輻條的交替,支撐輻條角度發(fā)生幅度為±2π/n 的躍變,且因碰撞引起輪輻整體動能損耗,導(dǎo)致碰撞前后輪輻角速度變化很大.顯然,碰撞前后的支撐輻條角度與角速度均不連續(xù).
當(dāng)無邊輪輻作上坡運(yùn)動時(shí),支撐輻條沿順時(shí)針方向滾動,可從左向右依次逆解讀圖2中輪輻運(yùn)動狀態(tài)來說明.圖2(b)、(f)為碰撞瞬間,碰撞點(diǎn)分別為點(diǎn)A 和點(diǎn)D,圖2(a)、(e)為碰撞前瞬間,圖2(c)為碰撞后瞬間,圖2(d)為支撐輻條呈豎直狀態(tài).圖2(c)、(d)、(e)中的輪輻處于無碰撞階段,其參數(shù)配置方式的遵循原則與下坡運(yùn)動時(shí)保持一致.
無邊輪輻無論處于何種初始運(yùn)動趨勢,在運(yùn)動過程中其運(yùn)動趨勢均可能出現(xiàn)翻轉(zhuǎn).輪輻到達(dá)如圖2(d)支撐角度為零狀態(tài)時(shí),角速度將決定輪輻是否可以通過支撐輻條來達(dá)到豎直狀態(tài),從而影響輪輻是否繼續(xù)初始運(yùn)動或出現(xiàn)搖擺運(yùn)動.
因此,本文定義無邊輪輻一個(gè)完整的運(yùn)動周期為一個(gè)完整的無碰撞運(yùn)動過程和接下來的一次碰撞運(yùn)動,并在不改變模型參數(shù)意義的基礎(chǔ)上簡化角度及角速度的符號標(biāo)示,以當(dāng)前支撐輻條角度θ(t)、角速度(t)作為研究對象,研究無邊輪輻的運(yùn)動特性.一個(gè)完整的輪輻運(yùn)動周期可表現(xiàn)為上坡運(yùn)動、下坡運(yùn)動或組合運(yùn)動,即在一個(gè)運(yùn)動周期內(nèi),輪輻運(yùn)動趨勢可能發(fā)生改變,或出現(xiàn)搖擺現(xiàn)象.
圖2 無邊輪輻模型運(yùn)動狀態(tài)示意簡圖
如圖2(c)、(d)、(e)所示,無論輪輻初始運(yùn)動趨勢為上坡或下坡運(yùn)動,輪輻都處于無碰撞運(yùn)動階段,整個(gè)輪輻關(guān)于支撐點(diǎn)A 的角動量矩平衡,即
式中,MA為輪輻關(guān)于點(diǎn)A 的所受力矩總和;˙HA為輪輻關(guān)于點(diǎn)A 的角動量變化率.
由式(1)可得無碰撞階段的無量綱運(yùn)動方程為
輪輻碰撞前后瞬間參數(shù)配置如圖2(a)~(c)所示,支撐輻條角度在碰撞前后有±2π/n 的躍變.無邊輪輻在碰撞發(fā)生前為逆時(shí)針下坡運(yùn)動,輪輻碰撞后瞬間相對于碰撞前瞬間轉(zhuǎn)過了-2π/n 角度;反之,若輪輻在碰撞發(fā)生前呈順時(shí)針上坡運(yùn)動趨勢,輪輻則轉(zhuǎn)過2π/n 角度.第i 次碰撞前后瞬間角度躍變公式為
無邊輪輻在碰撞發(fā)生前為逆時(shí)針下坡運(yùn)動,iθ-=α+π/n,iθ+=α -π/n;反之,若無邊輪輻在碰撞發(fā)生前為順時(shí)針上坡運(yùn)動,則iθ-=α -π/n,iθ+=α+ π/n.無論何種運(yùn)動趨勢下發(fā)生碰撞運(yùn)動,碰撞前后含符號函數(shù)的角度躍變公式(3)均可改寫為如下簡單形式:
無邊輪輻的角速度隨著碰撞的發(fā)生而改變,忽略重力及碰撞瞬間對輪輻引起的沖擊,則整個(gè)輪輻關(guān)于點(diǎn)B 角動量守恒,即
式中,H-B 為第i 次碰撞前瞬間輪輻關(guān)于點(diǎn)B 的角動量;H+B為第i 次碰撞后瞬間輪輻關(guān)于點(diǎn)B 的角動量.
由式(5)可得第i 次碰撞前后瞬間角速度躍變方程為
式中,碰撞參數(shù)λ=1 + μ2(cos(2π/n)-1),0 ≤λ <1.
第i 次碰撞前后瞬間角度、角速度躍變公式(4)、(6)表征了碰撞發(fā)生前后輪輻支撐輻條運(yùn)動參數(shù)的改變,式(4)、(6)可統(tǒng)稱為碰撞瞬間切換模型.
引進(jìn)臨界角速度參數(shù)dn和up,獲取依據(jù)為無邊輪輻介于2 次碰撞之間的運(yùn)動遵循能量守恒原則.在一個(gè)運(yùn)動周期的無碰撞階段,輪輻以一定的初始角速度運(yùn)動,輪輻的初始動能可根據(jù)圖2(c)、(d)、(e)中的動能和勢能表征條轉(zhuǎn)換.在θ=0 處完全轉(zhuǎn)化為勢能,即在θ=0 處,輪輻角速度亦為零,此時(shí)輪輻處于基于1 個(gè)支撐輻條的不穩(wěn)定豎直平衡狀態(tài).
當(dāng)無邊輪輻呈下坡運(yùn)動趨勢時(shí),輪輻以臨界角速度dn作初始運(yùn)動,可根據(jù)圖2(d)與(e)間的能量表征條進(jìn)行能量轉(zhuǎn)換;反之,當(dāng)輪輻呈上坡運(yùn)動趨勢時(shí),輪輻以臨界角速度up作為初始速度運(yùn)動,根據(jù)圖2(c)與(d)間能量表征條進(jìn)行能量轉(zhuǎn)換,計(jì)算可得
無邊輪輻呈下坡運(yùn)動趨勢且碰撞后瞬間角速度為dn或上坡運(yùn)動趨勢且碰撞后瞬間角速度為時(shí),輪輻將耗時(shí)無限長到達(dá)并停留在位置θ=0處,輪輻基于1 個(gè)支撐輻條豎直靜止,是一種臨界不穩(wěn)定平衡狀態(tài).
無碰撞階段的運(yùn)動方程式(2)可改寫成如下一階微分方程組:
積分后可得介于2 次碰撞之間的相空間軌跡流,即
鑒于支撐輻條與斜坡之間的空間限制,支撐輻條并不能進(jìn)行完整的倒立擺運(yùn)動,因此無邊輪輻的軌跡變化需在相空間U 中進(jìn)行研究,即
無碰撞階段運(yùn)動方程式(8)及碰撞瞬間切換模型(4)、(6),均可表達(dá)為如下形式:
2)圖3中相空間U 內(nèi)的連續(xù)運(yùn)動軌跡流,根據(jù)獨(dú)立向量場v 發(fā)生連續(xù)變化,表明輪輻在無碰撞階段的動力學(xué)特性.相空間邊界?U 上的碰撞切換軌跡流,依據(jù)獨(dú)立向量場h 發(fā)生離散性質(zhì)的跳躍變化,表明輪輻碰撞發(fā)生前后支撐輻條角度與角速度的躍變.一個(gè)完整的運(yùn)動周期軌跡流為連續(xù)運(yùn)動軌跡流與相應(yīng)的碰撞切換軌跡流的組合,具有連續(xù)動力學(xué)與離散動力學(xué)并存的復(fù)合動力學(xué)特性.
3)以所屬相空間邊界?U 內(nèi)含臨界值的邊界點(diǎn)x1(α-π/n,dn)與x2(α +π/n,up)為初始條件的臨界運(yùn)動軌跡流終止于不穩(wěn)定平衡位置θ=0處,輪輻不再出現(xiàn)后續(xù)運(yùn)動軌跡流.當(dāng)輪輻連續(xù)滾動時(shí),輪輻上一次碰撞后瞬間(θ,)即為輪輻下一個(gè)運(yùn)動周期的初始條件.若輪輻歷經(jīng)數(shù)次碰撞后,其運(yùn)動參數(shù)達(dá)到臨界邊界點(diǎn),則輪輻將到達(dá)并停留在θ=0 位置,輪輻將不會發(fā)生下一次碰撞.
4)如果斜坡角度α >n/π,輪輻角度θ 始終大于零,則相空間內(nèi)不包含不穩(wěn)定豎直平衡狀態(tài)θ=0.
圖3 兩次碰撞間的相空間軌跡流圖
以滿足式(13)的邊界點(diǎn)為初始運(yùn)動條件,圖3中連續(xù)運(yùn)動軌跡流表明,輪輻均未通過θ=0 位置,而是逆向終止于相空間邊界?U 上,輪輻出現(xiàn)運(yùn)動趨勢發(fā)生逆轉(zhuǎn)的搖擺運(yùn)動.反之,若(θ,)滿足
以滿足式(14)的邊界點(diǎn)為初始運(yùn)動條件,圖3中相空間U 內(nèi)連續(xù)運(yùn)動軌跡流均都正向終止于相空間邊界?U 上,輪輻運(yùn)動趨勢保持不變.
6)無論輪輻以何種初始條件運(yùn)動,歷經(jīng)數(shù)次碰撞、上坡、搖擺及下坡運(yùn)動后,輪輻可能處于如下3 種狀態(tài):①基于單個(gè)支撐輻條上的豎直狀態(tài),即運(yùn)動軌跡流終止于不穩(wěn)定平衡位置;②基于兩相鄰支撐輻條上的穩(wěn)定靜止?fàn)顟B(tài),即運(yùn)動軌跡流終止于相空間邊界?U 內(nèi)的=0 邊界點(diǎn)上;③輪輻以一定速度作持續(xù)下坡滾動.如圖4所示,設(shè)n=6,α=π/16,μ=0.75,輪輻以θ0=α -n/π,0=1 為初始條件開始運(yùn)動,歷經(jīng)足夠多次碰撞后,連續(xù)運(yùn)動軌跡流與碰撞切換軌跡流逐漸成為閉合曲線,說明極限環(huán)在特定條件下是存在的,輪輻可以實(shí)現(xiàn)穩(wěn)定持續(xù)的下坡滾動.
圖4 無邊輪輻相空間軌跡流及極限環(huán)
利用龐加萊回歸映射方法研究輪輻系統(tǒng)不動點(diǎn)、極限環(huán)的影響因素及存在條件,從而判定輪輻是否可實(shí)現(xiàn)穩(wěn)定持續(xù)的下坡滾動[5-7].龐加萊截面Γ 定義為
建立如下龐加萊映射關(guān)系:
式中,iq+=(iθ+,i+)∈Γ,映射函數(shù)f:?!#琭=(h ?v)是基于離散向量場h 和連續(xù)向量場v 的復(fù)合映射,向量場h 與v 有如下關(guān)系:
映射函數(shù)P 表示輪輻第i 次碰撞后與第i +1次碰撞后的角速度之間的關(guān)系.如式(19)所示,當(dāng)0≤α≤π/n 時(shí),函數(shù)P 是受臨界值dn,up和參數(shù)集合ψ={μ2,n,α}影響的分段連續(xù)函數(shù),從而輪輻在運(yùn)動過程中會出現(xiàn)交替上坡運(yùn)動、下坡運(yùn)動及包含2 種運(yùn)動的搖擺運(yùn)動;如式(20)所示,當(dāng)π/n<α≤π/2 時(shí),函數(shù)P 是受參數(shù)集合ψ={μ2,n,α}影響的連續(xù)線性函數(shù),若將無邊輪輻兩相鄰輻條支撐于斜坡上,即使給輪輻零初始角速度或上坡運(yùn)動趨勢,輪輻最終仍會自主持續(xù)下坡滾動.
映射函數(shù)P 表征了輪輻相鄰2 次碰撞后瞬間角速度之間的關(guān)系,函數(shù)P 的不動點(diǎn)的存在意味著輪輻歷經(jīng)數(shù)次碰撞后,其后每次碰撞后的瞬間角速度值均相同,這表明輪輻運(yùn)動存在極限環(huán),由式(19)、(20)可得函數(shù)P 的不動點(diǎn)ω*.
式中,ω0*為靜止不動點(diǎn),表明輪輻最終可實(shí)現(xiàn)基于兩相鄰輻條支撐的穩(wěn)定靜止?fàn)顟B(tài);為極限環(huán)不動點(diǎn),表明極限環(huán)的存在及輪輻最終可實(shí)現(xiàn)持續(xù)穩(wěn)定的下坡滾動,即輪輻歷經(jīng)數(shù)次碰撞后角速度值趨向且最終等于,其后每次碰撞后的瞬間角速度值均等于ω*1.
由不動點(diǎn)表達(dá)式(21)、(22)可知,不動點(diǎn)的存在與輪輻自身模型參數(shù)μ、輪輻數(shù)目n 及斜坡角度α 有關(guān).當(dāng)輪輻模型參數(shù)一定時(shí),函數(shù)P 的2 個(gè)不動點(diǎn)ω0*與ω*1與斜坡角度的關(guān)系為
式中,αe為臨界斜坡角度,是決定極限環(huán)不動點(diǎn)ω1*存在與否的臨界值,由式(21)可推出αe滿足
由圖5可知,不同的斜坡角度對應(yīng)不同的一維回轉(zhuǎn)圖,不動點(diǎn)性質(zhì)及存在個(gè)數(shù)也不同.
如圖5(a)所示,當(dāng)0≤α≤αe時(shí),輪輻存在靜止不動點(diǎn)ω*0及表征不穩(wěn)定平衡位置的臨界值up和dn,不存在極限環(huán)不動點(diǎn)ω*1 ,輪輻最終處于穩(wěn)定靜止?fàn)顟B(tài)或不穩(wěn)定平衡狀態(tài).
如圖5(b)所示,當(dāng)αe<α≤π/n 時(shí),輪輻存在靜止不動點(diǎn)ω0*、極限環(huán)不動點(diǎn)ω1*及臨界值和等.輪輻以不同的初始條件運(yùn)動,歷經(jīng)數(shù)次碰條件取決于角速度臨界值撞,輪輻角速度逐漸趨向于ω0*,ω1*,up或dn等4種角速度值,此時(shí)輪輻最終運(yùn)動狀態(tài)可能有3 種:①輪輻基于兩相鄰輻條搖擺數(shù)次,角速度在有限時(shí)間內(nèi)衰減至靜止不動點(diǎn)ω*0,輪輻最終靜止于兩相鄰輻條上;②輪輻碰撞后瞬間角速度值以增速和減速2 種方式逐漸趨向且最終等于極限環(huán)不動點(diǎn)ω1*,輪輻實(shí)現(xiàn)以ω1*為運(yùn)動角速度的穩(wěn)定持續(xù)的下坡滾動;③碰撞后瞬間輪輻角速度達(dá)到臨界值up和dn,然后輪輻將耗時(shí)無限長到達(dá)并停留在豎直位置,即不穩(wěn)定平衡狀態(tài)θ=0 位置.
圖5 不同斜坡角度對應(yīng)的不動點(diǎn)與相應(yīng)的一維回轉(zhuǎn)圖(n=6,μ=0.75)
不動點(diǎn)的穩(wěn)定性由映射函數(shù)在不動點(diǎn)處的一階導(dǎo)數(shù)決定.函數(shù)P 在不動點(diǎn)ω*0與ω*1處的一階導(dǎo)數(shù)為
由于0≤λ <1,映射函數(shù)P 在不動點(diǎn)ω*0與ω1*處的一階導(dǎo)數(shù)DP(ω*)滿足因此,映射函數(shù)P 的不動點(diǎn)ω*0與ω*1是漸進(jìn)穩(wěn)定的.
當(dāng)α=π/15,n=6,μ=0.75,回轉(zhuǎn)圖如圖6所示.圖6表示了輪輻靜止不動點(diǎn)ω*0與極限環(huán)不動點(diǎn)ω*1及臨界值dn和up等均存在時(shí)的多種運(yùn)動狀態(tài).靜止不動點(diǎn)ω*0與極限環(huán)不動點(diǎn)ω*1交替吸引相應(yīng)區(qū)域內(nèi)的初始條件,從而決定輪輻運(yùn)動狀態(tài).吸引盆在回轉(zhuǎn)圖中表現(xiàn)為不同區(qū)域內(nèi)的線段,吸引盆的邊界由角速度臨界值來界定.圖6中奇數(shù)值區(qū)域內(nèi)線段為極限環(huán)不動點(diǎn)的吸引盆,偶數(shù)值區(qū)域內(nèi)線段為靜止不動點(diǎn)的吸引盆.由此可知,滿足極限環(huán)存在的初始所界定的區(qū)域,且不動點(diǎn)ω*0與ω*1相對應(yīng)的吸引盆交替出現(xiàn).
圖6 臨界角速度值與相應(yīng)不動點(diǎn)的回轉(zhuǎn)圖
當(dāng)輪輻以式(26)中臨界值系列為初始條件開始運(yùn)動時(shí),其終止于θ=0 位置而無法到達(dá)任一不動點(diǎn),由映射函數(shù)P 可推導(dǎo)出m≥1 時(shí)的系列臨界值,其表達(dá)式為
綜上所述,輪輻最終所處運(yùn)動狀態(tài)依賴于斜坡角度及由臨界角速度值界定的不動點(diǎn)ω0*與ω1*吸引盆內(nèi)的初始角速度.斜坡角度、初始角速度與不動點(diǎn)類型的具體依存關(guān)系如表1所示,輪輻最終運(yùn)動狀態(tài)與初始運(yùn)動趨勢如圖7所示.①若輪輻初始角速度或歷經(jīng)數(shù)次碰撞后的角速度為臨界角速度值,無論斜坡角度如何,輪輻終將處于不穩(wěn)定豎直狀態(tài);②當(dāng)斜坡角度為0≤α≤αe時(shí),初始角速度為非臨界角速度值,所有角速度值都位于靜止不動點(diǎn)ω*0的吸引盆內(nèi),即無論輪輻為何種初始運(yùn)動趨勢,輪輻終將歷經(jīng)M 次搖擺后基于兩相鄰支撐輻條處于穩(wěn)定靜止?fàn)顟B(tài);③當(dāng)斜坡角度滿足αe<α≤π/n,且角速度滿足或則所有角速度值都位于靜止不動點(diǎn)ω*0的吸引盆內(nèi),即輪輻歷經(jīng)M 次搖擺后將基于兩相鄰支撐輻條處于穩(wěn)定靜止?fàn)顟B(tài);④當(dāng)斜坡角度滿足αe<α≤π/n 且初始角速度滿足則角速度值均位于極限環(huán)不動點(diǎn)ω1*的吸引盆內(nèi),輪輻由最初的上坡滾動通過一次搖擺逆轉(zhuǎn)為下坡滾動,歷經(jīng)數(shù)次碰撞,輪輻角速度最終增速趨向并等于極限環(huán)不動點(diǎn)ω*1,輪輻呈穩(wěn)定持續(xù)下坡滾動狀態(tài);⑤當(dāng)斜坡角度與初始角速度滿足關(guān)系αe<α≤π/n,>dn或π/n <α≤π/2,≥0 時(shí),輪輻會以增速或減速2 種方式趨向且最終等于極限環(huán)不動點(diǎn)ω1*,輪輻實(shí)現(xiàn)純被動性質(zhì)的穩(wěn)定持續(xù)下坡滾動.
表1 不動點(diǎn)與相應(yīng)吸引盆及斜坡角度的依存關(guān)系
圖7 輪輻初始條件與運(yùn)動過程及最終運(yùn)動狀態(tài)關(guān)系圖
1)輪輻模型參數(shù)一定時(shí),不同的斜坡角度對應(yīng)不同的龐加萊映射回轉(zhuǎn)圖及不動點(diǎn)的存在狀態(tài),是輪輻運(yùn)動極限環(huán)存在的重要影響因素.決定極限環(huán)不動點(diǎn)存在的臨界斜坡角度值依賴于輪輻輻條數(shù)目及慣性矩參數(shù).
2)輪輻自身模型參數(shù)及所在斜坡角度一定時(shí),靜止不動點(diǎn)與極限環(huán)不動點(diǎn)的吸引盆區(qū)域確定,輪輻初始運(yùn)動條件所屬吸引盆類別直接決定輪輻最終運(yùn)動狀態(tài).
3)臨界角速度系列為不動點(diǎn)的吸引盆邊界值,由輪輻自身模型參數(shù)及斜坡角度決定.歷經(jīng)數(shù)次碰撞后,以此系列臨界角速度為初始條件開始運(yùn)動的輪輻將耗時(shí)無限長,趨向并停留于不穩(wěn)定豎直平衡位置,而無法到達(dá)靜止不動點(diǎn)或極限環(huán)不動點(diǎn).
4)輪輻只在重力作用下發(fā)生滾動,無附加驅(qū)動力,具有純被動性質(zhì),輪輻自身模型參數(shù)、斜坡角度、初始運(yùn)動條件將決定輪輻運(yùn)動過程及最終運(yùn)動狀態(tài).
5)輪輻自身模型參數(shù)即輻條數(shù)目、慣性矩參數(shù)與臨界斜坡角度、臨界角速度值、不動點(diǎn)、不動點(diǎn)的吸引盆相互依賴并相互制約.
無邊輪輻輻條數(shù)目與慣性矩參數(shù)類比對應(yīng)于簡單被動雙足行走模型運(yùn)動過程中雙腿間最大角度、腿部與髖關(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].機(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ū)動步行機(jī)器人運(yùn)動學(xué)機(jī)理與控制策略研究[D].長春:吉林大學(xué)通信工程學(xué)院,2007.
[15] 劉振澤,周長久,田彥濤.Compass-like 無動力雙足行走機(jī)器人的運(yùn)動狀態(tài)[J].吉林大學(xué)學(xué)報(bào):工學(xué)版,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)行走雙足機(jī)器人的穩(wěn)定性分析與控制研究[D].長春:吉林大學(xué)通信工程學(xué)院,2011.
[17] 胡運(yùn)富.無膝雙足被動機(jī)器人的運(yùn)動特性和穩(wěn)定性研究[D].哈爾濱:哈爾濱工業(yè)大學(xué)機(jī)電工程學(xué)院,2009.
[18] 郭欣.雙足被動機(jī)器人行走模式的設(shè)計(jì)構(gòu)想[D].北京:北京郵電大學(xué)自動化學(xué)院,2009.
[19] 胡運(yùn)富,朱延河,吳曉光,等.被動無邊輪輻運(yùn)動特性的Adams 仿真分析[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào),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)