王亞翔 張俊俊
(西南科技大學(xué)特殊環(huán)境機(jī)器人技術(shù)四川省重點(diǎn)實(shí)驗(yàn)室 四川綿陽(yáng) 621010)
履帶是履帶機(jī)器人中的重要部件之一。其功能是保證機(jī)器人在無(wú)路地面上的通過(guò)性,降低車(chē)輛的行駛阻力。而履帶機(jī)器人的通過(guò)穩(wěn)定性可以通過(guò)履帶張緊力來(lái)評(píng)估,履帶的張緊力對(duì)于履帶行駛系統(tǒng)的性能有深刻影響。張緊力過(guò)大,履帶剛度太大,起不到緩沖作用,會(huì)增加履帶和底盤(pán)部件摩擦,能耗嚴(yán)重,并有斷帶危險(xiǎn);張緊力過(guò)小,使得履帶松弛,容易脫帶、耙齒,造成履帶失效,車(chē)輛癱瘓[1-4]。為了使機(jī)器人更好地適應(yīng)環(huán)境、提高機(jī)器人在遙操作過(guò)程中的穩(wěn)定性和高效性,有必要對(duì)其在典型環(huán)境中工作時(shí)張緊力的變化進(jìn)行研究,以提高自適應(yīng)履帶機(jī)器人的綜合性能,方便實(shí)際應(yīng)用中設(shè)計(jì)合理的機(jī)器人結(jié)構(gòu)并明確對(duì)機(jī)器人的操作和控制[5],因此張緊力的計(jì)算和優(yōu)化成為研究的重點(diǎn)。
國(guó)內(nèi)外許多學(xué)者對(duì)履帶張緊力進(jìn)行了分析,Park等[6]將低速履帶環(huán)看成柔性帶,通過(guò)建立履帶的數(shù)學(xué)模型分析履帶車(chē)輛的行駛性能,Wong等[7]以NTVPM模型預(yù)測(cè)小型、輕型履帶車(chē)輛在可變形地形上的越野性能。劉漢光等[8]研究了履帶液壓挖掘機(jī)行走裝置的合理預(yù)張緊力。馬興國(guó)等[9]通過(guò)建立數(shù)學(xué)模型得到了履帶系統(tǒng)各部位的張緊力以及地面對(duì)負(fù)重輪法向力的計(jì)算公式。本文根據(jù)自行設(shè)計(jì)的自適應(yīng)履帶式機(jī)器人地盤(pán)為模型,通過(guò)理論分析研究履帶機(jī)器人張緊力分布,得出了主動(dòng)輪附近張緊力計(jì)算公式,并通過(guò)Recurdyn軟件的Track(LM)工具建立移動(dòng)平臺(tái)的虛擬樣機(jī)模型驗(yàn)證了計(jì)算公式的可行性。同時(shí)通過(guò)正交實(shí)驗(yàn)得出履帶機(jī)器人設(shè)計(jì)中的主要參數(shù)對(duì)張緊力的影響顯著性,為履帶式機(jī)器人設(shè)計(jì)和優(yōu)化提供依據(jù)。
本文設(shè)計(jì)的履帶機(jī)器人機(jī)構(gòu)如圖1所示,其中為單邊履帶系統(tǒng),包括了主動(dòng)輪1、誘導(dǎo)輪3、負(fù)重輪10、電機(jī)與減速器8、電源7和主控6,履帶單邊可以單獨(dú)提供數(shù)據(jù)傳輸與控制,獨(dú)立的電源與電機(jī)可以是履帶單邊獨(dú)立完成行走任務(wù)。其中為了提高越障性能,前進(jìn)方向處的負(fù)重輪處設(shè)計(jì)了被動(dòng)自適應(yīng)變形結(jié)構(gòu),包括了變形板4和負(fù)重輪10以及變形減震彈簧11。變形板在越障受力后被動(dòng)變形,減小前角大小從而增加越障性能。
1.主動(dòng)輪, 2.托帶輪, 3.誘導(dǎo)輪, 4.變形板, 5.驅(qū)動(dòng)器,
6.主控通信模塊, 7.電池, 8.動(dòng)力模塊(電機(jī)與減速器),
9.履帶, 10.負(fù)重輪, 11.彈簧
圖1履帶機(jī)器人結(jié)構(gòu)
Fig.1 Structure of tracked robot
履帶機(jī)器人的履帶系統(tǒng)由主動(dòng)輪、托帶輪、誘導(dǎo)輪、負(fù)重輪組成,如圖2所示。由于設(shè)計(jì)的履帶機(jī)器人其托帶輪與主動(dòng)輪和誘導(dǎo)輪相切履帶為直線,因此忽略主動(dòng)輪和誘導(dǎo)輪之間的履帶板的重力和振動(dòng)以及托帶輪的影響,分析時(shí)將履帶系統(tǒng)分為主動(dòng)輪、負(fù)重輪、誘導(dǎo)輪三部分,則可以認(rèn)為T(mén)d1≈Tc2。
圖2 履帶張緊力系統(tǒng)模型Fig.2 Mechanical model for tension system
取主動(dòng)輪為隔離體,得到主動(dòng)輪受力如圖3(a)所示,可以得出主動(dòng)輪處張緊力為:
(1)
式中:M為主動(dòng)輪轉(zhuǎn)矩;Td1,Td2為主動(dòng)輪處松邊張緊力和主動(dòng)輪處緊邊張緊力。
圖3 主從輪系受力分析Fig.3 Mechanical model for driving wheel and induction wheel
取誘導(dǎo)輪為隔離體,得到誘導(dǎo)輪輪受力如圖3(b)所示,可以得出誘導(dǎo)輪處張緊力動(dòng)力學(xué)方程為:
(2)
式中:mc為從動(dòng)輪質(zhì)量;Fct為履帶預(yù)張緊力,F(xiàn)c為誘導(dǎo)輪周?chē)膸D(zhuǎn)的離心合力;γ為離心合力與水平夾角;Tc1,Tc2為誘導(dǎo)輪處履帶受張緊力松邊、誘導(dǎo)輪處履帶受張緊力緊邊;α為履帶機(jī)器人前角,Gc為誘導(dǎo)輪受到的重力;Jc為履帶誘導(dǎo)輪轉(zhuǎn)動(dòng)慣量;Rc為誘導(dǎo)輪半徑;ωc為履帶機(jī)器人勻速前進(jìn)誘導(dǎo)輪角速度;Fcn為誘導(dǎo)輪y方向支撐力;ρ為履帶單位質(zhì)量。
取履帶負(fù)重輪系為隔離,不變形負(fù)重輪系受力如圖4(a)所示,變形負(fù)重輪系受力如圖4(b)所示。
圖4 負(fù)重輪系分析Fig.4 Mechanical model for load wheel
得出不變形負(fù)重輪處履帶張緊力動(dòng)力學(xué)方程為:
(3)
(4)
被動(dòng)自適應(yīng)變形處的負(fù)重輪履帶張緊力動(dòng)力學(xué)方程為:
(5)
(6)
式中:ml為負(fù)重輪質(zhì)量;TLi1,TLi2為各負(fù)重輪處張緊力;β為履帶后角;α為履帶前角;RL為負(fù)重輪半徑;Jl為負(fù)重輪轉(zhuǎn)動(dòng)慣量;fLi為各負(fù)重輪與履帶間的摩擦力;FLin,F(xiàn)Lit為負(fù)重輪y,x方向支撐力;ωLi為各負(fù)重輪角速度;Gl為負(fù)重輪受到的重力;NLi為各負(fù)重輪與履帶間的法線接觸力,u·NLi=fLi;u為履帶與負(fù)重輪間摩擦系數(shù),u1為履帶與負(fù)重輪間滾動(dòng)摩擦系數(shù)。
綜合式(1)-式(6)可以得出主動(dòng)輪處兩端的張緊力分別為:
(7)
對(duì)于變形負(fù)重輪處的變形板有:
(8)
其中:
(9)
式中:W為履帶機(jī)器人對(duì)在變形板上的均布?jí)毫?;mB為變形板質(zhì)量;ml為負(fù)重輪質(zhì)量;φ為變形板斜邊與法線夾角;LB為變形板邊長(zhǎng);TB為變形板受到的轉(zhuǎn)矩;k為彈簧剛度;θB為變形板微變形后變化角度。
(10)
從式中可以看出,自適應(yīng)履帶主動(dòng)輪附近張緊力與主動(dòng)輪、從動(dòng)輪半徑、質(zhì)量、自適應(yīng)負(fù)重輪處彈簧剛度、轉(zhuǎn)動(dòng)慣量、前角和后角角度相關(guān)。
為了驗(yàn)證張緊力推導(dǎo)公式的可靠性,在多體動(dòng)力學(xué)軟件Recurdyn中使用Track模塊建立了被動(dòng)自適應(yīng)履帶機(jī)器人的動(dòng)力學(xué)模型,當(dāng)選用表1參數(shù)進(jìn)行理論計(jì)算和仿真實(shí)驗(yàn)時(shí)(履帶預(yù)張緊力以總質(zhì)量百分比為單位),得到的仿真曲線對(duì)比圖如圖5所示。通過(guò)文中推導(dǎo)公式計(jì)算張緊力Td1平均值為1 134.3 N,仿真計(jì)算張緊力Td1平均值為1 092.1 N,通過(guò)圖5可以看出兩個(gè)曲線變化規(guī)律相符合,因此公式準(zhǔn)確性較好。
表1 仿真參數(shù)Table 1 Parameters of simulation
圖5 張緊力Td1公式計(jì)算與仿真對(duì)比Fig.5 Comparison between calculation of tension force Td1 and that of simulation
從式(10)可以得出當(dāng)機(jī)器人質(zhì)量、尺寸參數(shù)確定時(shí),改變其中任何一個(gè)設(shè)計(jì)參數(shù)都會(huì)引起對(duì)主動(dòng)輪附近張緊力的影響。為了得出多個(gè)設(shè)計(jì)參數(shù)對(duì)運(yùn)行過(guò)程中履帶張緊力變化影響的顯著性,本文將前角α、后角β、預(yù)緊力Fct、變形彈簧剛度k4個(gè)重要履帶設(shè)計(jì)參數(shù)作為正交實(shí)驗(yàn)因素,每個(gè)因素選取4個(gè)水平,選取L(45)正交表進(jìn)行實(shí)驗(yàn),正交實(shí)驗(yàn)表如表2所示,16次正交實(shí)驗(yàn)得出由理論計(jì)算與仿真實(shí)驗(yàn)的張緊力均值變化圖如圖6所示。從表3和圖6可以看出正交實(shí)驗(yàn)中理論公式計(jì)算的張緊力Td1的變化趨勢(shì)與仿真相一致且計(jì)算得出的張緊力均值相近,因此可以進(jìn)一步驗(yàn)證理論公式的有效性,同時(shí)在極差與方差分析中將會(huì)采用正交表中的仿真均值作為分析數(shù)據(jù)。根據(jù)表2數(shù)據(jù)得出的4因素方差分析結(jié)果如表3所示??梢愿鶕?jù)表2的極差和表3的方差分析看出前角α、后角β、預(yù)緊力Fct對(duì)實(shí)驗(yàn)結(jié)果有顯著性影響,其數(shù)值的改變會(huì)引起張緊力數(shù)值的顯著性變化,其中預(yù)緊力Fct對(duì)實(shí)驗(yàn)結(jié)果影響極其顯著。
表2 正交實(shí)驗(yàn)設(shè)計(jì)Table 2 Orthogonal test design
圖6 正交實(shí)驗(yàn)中張緊力變化Fig.6 The change of tension in orthogonal test
因素偏差平方和自由度方差FProb>Fα1.11e433 732.1490.621.56e-4β321.83107.314.120.028Fct3.62e531.21e515 850.198.5e-7k9.633.20.420.75誤差e22.837.6總和3.73e515
通過(guò)利用正交實(shí)驗(yàn)中的數(shù)據(jù)在Recurdyn中建立了動(dòng)力學(xué)模型,仿真得出了16次試驗(yàn)中履帶機(jī)器人在法線方向的加速的變化,如圖7所示。從表3和圖7可以看出,履帶整體尺寸確定時(shí),前角α與后角β過(guò)小時(shí),履帶法線加速度及其均方差較大,穩(wěn)定性較差;當(dāng)前角α與后角β增大后,預(yù)緊力Fct的變化對(duì)其法線方向加速度的均方差影響較大;其中第11次實(shí)驗(yàn)中履帶機(jī)器人的法線加速度和均方差最小,初始預(yù)張緊力為整體重力的50%,這時(shí)履帶主動(dòng)輪處張緊力均值為710.2 N,相比正交實(shí)驗(yàn)的其他實(shí)驗(yàn)中履帶張緊力處于中間水平,可以同時(shí)綜合脫帶與減震的情況。
圖7 履帶機(jī)器人法線方向加速度變化Fig.7 Acceleration change of crawler robot in normal direction
通過(guò)建立數(shù)學(xué)模型、理論分析、仿真以及正交實(shí)驗(yàn),得出以下結(jié)論:(1)建立了被動(dòng)自適應(yīng)履帶機(jī)器人負(fù)重輪、主動(dòng)輪、誘導(dǎo)輪附近的履帶張緊力數(shù)學(xué)模型,得到了被動(dòng)自適應(yīng)履帶機(jī)器人張緊力的計(jì)算公式,并通過(guò)用公式計(jì)算得到的曲線與仿真結(jié)果相比較,計(jì)算公式與仿真實(shí)驗(yàn)得出的張緊力曲線變化趨勢(shì)相符合,驗(yàn)證了論文建立的被動(dòng)自適應(yīng)履帶張緊力公式的可行性。(2)通過(guò)建立正交實(shí)驗(yàn),得到了被動(dòng)自適應(yīng)履帶機(jī)器人多個(gè)結(jié)構(gòu)參數(shù)對(duì)履帶張緊力的影響顯著性。得出了在多個(gè)結(jié)構(gòu)參數(shù)同時(shí)影響的情況下,預(yù)緊力對(duì)自適應(yīng)履帶運(yùn)行中張緊力影響顯著性最大。并得到了履帶機(jī)器人以前角45°、后角45°,預(yù)緊力為總重力的50%~60%時(shí)設(shè)計(jì)時(shí),自適應(yīng)履帶機(jī)器人直線行走穩(wěn)定性最優(yōu)。