• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    一種諧振式壓電爬行機(jī)器人的設(shè)計與實驗1)

    2023-10-29 10:15:50何世界吳醫(yī)博周生喜
    力學(xué)學(xué)報 2023年9期
    關(guān)鍵詞:壓電驅(qū)動頻率

    何世界 吳醫(yī)博 周生喜,2)

    * (西北工業(yè)大學(xué)航空學(xué)院,西安 710072)

    ? (洛陽船舶材料研究所,河南洛陽 471023)

    引言

    近些年來,機(jī)器人技術(shù)在國內(nèi)越來越受到重視,同時也在各領(lǐng)域的生產(chǎn)實踐中得到發(fā)展與應(yīng)用.爬行機(jī)器人[1]作為機(jī)器人中的一種,將代替人類執(zhí)行偵察、探測、工程作業(yè)及顯微操作等多項任務(wù),它們的高精確性、高效率、高耐久度證明機(jī)器人可以更好地完成人類的工作任務(wù).微小型機(jī)器人在如今高科技信息時代扮演著至關(guān)重要的角色,是近些年來智能機(jī)器人技術(shù)發(fā)展研究的重點(diǎn)方向,基于其體積較小、靈敏度高、抗干擾強(qiáng)、運(yùn)動靈活等優(yōu)點(diǎn),可以應(yīng)用于災(zāi)后救援搜索、極端環(huán)境探測和醫(yī)療手術(shù)等諸多領(lǐng)域.其中,在管道、峭壁、沙漠以及醫(yī)療過程中一些人體部位等人類不易工作的地點(diǎn),微小型爬行機(jī)器人能夠代替人類完成高精度任務(wù),進(jìn)而可有效降低人工作業(yè)的危險性并提高工作效率.

    由于微小型爬行機(jī)器人的基體較小,承載能力較低,當(dāng)我們需要將機(jī)器人體積縮小到cm 或mm級別時,傳統(tǒng)的機(jī)器人結(jié)構(gòu)模型與驅(qū)動方式將無法適用,這就要求爬行機(jī)器人利用較小的驅(qū)動腿結(jié)構(gòu)來產(chǎn)生較大的驅(qū)動力.與此同時,電池與傳感器的重量也必須輕量化,這些問題對機(jī)器人的微小型化來說都是巨大的挑戰(zhàn).

    哈佛大學(xué)的微型機(jī)器人實驗室(Microrobotics Laboratory)研制了一種微型多節(jié)爬行機(jī)器人[2],該機(jī)器人重量為2.4 g,采用PZT-5H 壓電陶瓷作為驅(qū)動器,在200 V 電壓的驅(qū)動下能達(dá)到0.5 Hz 的腿部擺動頻率.斯坦福大學(xué)的工程師們設(shè)計出了一款叫做MicroTugs 的微型爬行機(jī)器人[3-4],該機(jī)器人僅重12 g,在平面上的投影與一枚硬幣相當(dāng).爬行機(jī)器人的腿上安裝有受到外界壓力時會自動彎曲的橡膠釘,增加了其與接觸面之間的相對面積從而加強(qiáng)了機(jī)器人腿部與接觸面之間的黏附性.2010 年,美國哈佛大學(xué)的Baisch 等[5]研制了一種仿蟑螂運(yùn)動的微小型六足機(jī)器人HAMR,每條支腿可實現(xiàn)兩個自由度的運(yùn)動輸出,并通過撓性接頭與基體相連.2011年,美國哈佛大學(xué)的Baisch 等[6]研制了新一代仿蟑螂運(yùn)動的微小型六足機(jī)器人 HAMR-3,配備了專用的小型電源,實現(xiàn)了無源運(yùn)動.2014 年,Baisch 等[7]在前期工作的基礎(chǔ)上設(shè)計制造出了一種可組裝的微型四足機(jī)器人HAMR-VP,與HAMR 相比,前者降低了制造復(fù)雜性且同時實現(xiàn)了準(zhǔn)靜態(tài)和動態(tài)操作.

    在微小型機(jī)器人研究中,利用智能材料進(jìn)行機(jī)器人的驅(qū)動可以使運(yùn)動結(jié)構(gòu)更加簡單,制造更加便捷.其中,利用壓電陶瓷逆壓電效應(yīng)[8]的新型驅(qū)動方式具有占用體積小、傳動效率高、響應(yīng)速度快且不易受干擾等諸多優(yōu)點(diǎn),對提高機(jī)器人的可操控性、運(yùn)動分辨力和響應(yīng)速度等十分有利,在極端環(huán)境的探測和物品運(yùn)輸?shù)阮I(lǐng)域有著潛在的應(yīng)用前景.由于以上各種優(yōu)勢,壓電陶瓷材料驅(qū)動為微小型爬行機(jī)器人的未來發(fā)展提供了新方向,因此,對于復(fù)雜環(huán)境下的作業(yè),利用逆壓電效應(yīng)和摩擦驅(qū)動以及黏滑運(yùn)動原理設(shè)計出的各種新型結(jié)構(gòu)的壓電爬行機(jī)器人具有非常廣的研究前景和實用價值.

    2014 年,Takato 等[9]提出了一種壓電沖擊式旋轉(zhuǎn)驅(qū)動器,致動器由多層壓電元件組成,該毫米大小的機(jī)器人有6 條腿,可以像昆蟲一樣在不平整的表面行走,每側(cè)3 支腿通過連接機(jī)構(gòu)連接.2016 年,Rios 等[10-11]設(shè)計了一款壓電驅(qū)動多足式機(jī)器人,利用機(jī)器人足部在不同驅(qū)動電壓下的不同模態(tài)實現(xiàn)機(jī)器人的運(yùn)動.2016 年,Graule 等[12]利用壓電驅(qū)動器驅(qū)動仿蜻蜓機(jī)器人的翅膀,設(shè)計并制造了一個撲翼型機(jī)器人.2017 年,Hariri 等[13]設(shè)計了一款駐波驅(qū)動的腿式壓電微型機(jī)器人并建立了其運(yùn)動學(xué)模型,設(shè)計的機(jī)器人利用壓電梁的不同模態(tài)進(jìn)行運(yùn)動.2017 年,Rios 等[14]提出并設(shè)計了一種每條驅(qū)動腿的運(yùn)動由兩個壓電雙晶梁來控制的微小型六足機(jī)器人MinRAR V1,當(dāng)受到350 Hz 頻率的信號激勵時可實現(xiàn)約為520 mm/s 的最大運(yùn)動速度.2018 年,Rios 等[10]研制出了新一代的微小型六足機(jī)器人MinRAR V2,與之前相比,該機(jī)器人裝配了驅(qū)動電源與微型控制系統(tǒng),可以實現(xiàn)兩種驅(qū)動模式的運(yùn)動,即利用外部電源激勵運(yùn)動和用自身的電源激勵運(yùn)動.

    在國內(nèi),蔣振宇等[15-16]研制了一種基于雙壓電膜驅(qū)動的微小型足式機(jī)器人,該機(jī)器人每條腿都是先在平面彈性板上加工出平面結(jié)構(gòu)再彎折而成,結(jié)構(gòu)簡單.2015 年,李魁等[17]設(shè)計了一種有3 條不同固有頻率的圓弧曲梁足的微小型三足機(jī)器人,通過施加每個曲梁足所對應(yīng)固有頻率的激勵信號,可以使機(jī)器人沿3 個方向運(yùn)動.2017 年,鄭龍龍等[18-19]設(shè)計了一種由壓電雙晶片驅(qū)動的微小型八足機(jī)器人,其4 條內(nèi)腿和4 條外腿足端與地面均有大小相同的一個傾斜角,且每一條腿只有一個自由度.2018 年,陳暢等[20-21]設(shè)計并研制了一種六足爬行機(jī)器人,利用多層材料疊合的壓電驅(qū)動器驅(qū)動,機(jī)身結(jié)構(gòu)由3 個驅(qū)動單元節(jié)組成并利用碳纖維連桿傳動機(jī)構(gòu)放大位移并傳動.2019 年,李一帆等[22]和李京等[23]研究了一種輪腿式的足腿一體化六足壓電機(jī)器人,每條腿利用4 個壓電陶瓷片諧振驅(qū)動,并且在實驗中該機(jī)器人可以在不同粗糙度的表面行走,展現(xiàn)出了較強(qiáng)的環(huán)境適應(yīng)性.2019 年,Peng 等[24]設(shè)計了一種微小型雙足壓電機(jī)器人Milli-Walker,主要由一片壓電陶瓷和一個雙足基體組成,體長僅為6 mm,當(dāng)施加的激勵信號達(dá)到機(jī)器人共振頻率時,可以實現(xiàn)較快的運(yùn)動速度.2020 年,盧鵬輝[25]設(shè)計了一種小型防傾翻結(jié)構(gòu)的壓電機(jī)器人,該結(jié)構(gòu)在傾翻后仍可以各個姿態(tài)運(yùn)動,并利用LabVIEW 軟件構(gòu)建了控制系統(tǒng),能使機(jī)器人按既定路徑運(yùn)動.2021 年,高煜斐等[8]設(shè)計并制作了一種微型結(jié)構(gòu)的三足壓電爬行機(jī)器人,通過驅(qū)動腿的諧振原理來使其行走,利用左右兩側(cè)不同的驅(qū)動電壓使機(jī)器人實現(xiàn)了較快的左轉(zhuǎn)和右轉(zhuǎn)運(yùn)動,并實現(xiàn)了設(shè)計的3 個方向的運(yùn)動,該結(jié)構(gòu)可為未來微型壓電機(jī)器人的設(shè)計與研究提供參考.為了提高壓電機(jī)器人的移動靈活性,人們對提高機(jī)器人的移動速度和運(yùn)動分辨力產(chǎn)生了很大興趣.2023年,Li 等[26]在前期研究的基礎(chǔ)上研制了一種小型四足壓電機(jī)器人,該結(jié)構(gòu)模擬動物運(yùn)動時的姿態(tài)產(chǎn)生形變,可以分別實現(xiàn)納米分辨率和高速移動兩種運(yùn)動狀態(tài),并可實現(xiàn)上述兩種運(yùn)動機(jī)制的有效配合.2023 年,Deng 等[27]提出了一種小型諧振蠕動壓電機(jī)器人,該結(jié)構(gòu)具有沿圓周均勻設(shè)置的6 個驅(qū)動腳,通過兩種振動模式的疊加使驅(qū)動腳產(chǎn)生橢圓運(yùn)動,該機(jī)器人最大速度可達(dá)200 mm/s,位移分辨力為0.71 μm.同年,Chen 等[28]研究了一種可快速移動和翻轉(zhuǎn)的軟壓電機(jī)器人,長10 mm,體重僅0.058 g.受獵豹四肢在奔馳過程中收縮和擴(kuò)張的運(yùn)動步態(tài)的啟發(fā),利用雙螺旋的共振用于放大運(yùn)動性能,該軟體機(jī)器人可在復(fù)雜地形上移動、爬坡,并承受重負(fù)載以及溫度傳感器.

    目前關(guān)于小型壓電機(jī)器人的實驗研究在現(xiàn)有的依靠逆壓電效應(yīng)的執(zhí)行機(jī)構(gòu)驅(qū)動位移較小以及不能控制其進(jìn)行復(fù)雜環(huán)境下移動與作業(yè)等方面尚存在不足,因此本文在上述研究的基礎(chǔ)上,設(shè)計了一種結(jié)構(gòu)簡單、響應(yīng)快速的四足爬行機(jī)器人.基于歐拉-伯努利梁理論計算了其驅(qū)動腿部的固有頻率,并推導(dǎo)了機(jī)器人基體結(jié)構(gòu)的力學(xué)模型,設(shè)計并制作了壓電機(jī)器人實物,通過實驗測試分析了不同驅(qū)動頻率、不同負(fù)載、不同電壓、不同驅(qū)動足以及不同接觸面對運(yùn)動方向及運(yùn)動速度的影響.與此同時,本文還設(shè)計了一種機(jī)器人的控制程序,可以使其執(zhí)行期望的運(yùn)動,以期為后續(xù)進(jìn)一步優(yōu)化機(jī)器人結(jié)構(gòu)以用于人類無法進(jìn)入的環(huán)境或者探險救災(zāi)等提供參考數(shù)據(jù).

    1 機(jī)器人結(jié)構(gòu)設(shè)計與驅(qū)動原理

    1.1 機(jī)器人的模型設(shè)計

    在對稱結(jié)構(gòu)中,要想實現(xiàn)運(yùn)動必須實現(xiàn)不對稱的驅(qū)動力,借鑒于不同變化角度的傾斜面與地面接觸時的前后摩擦系數(shù)不同,設(shè)想出了如圖1 所示的3 種驅(qū)動足,分別為圓弧體、圓弧面(鏟式)和傾斜面驅(qū)動足,驅(qū)動力實現(xiàn)方式主要為利用不對稱的形狀和表面傾斜結(jié)構(gòu)形成不對稱摩擦力[29].選擇在材料表面加工傾斜的微結(jié)構(gòu)實現(xiàn)各向異性的摩擦力,該驅(qū)動足采用PLA 材料使用30%填充度打印而成,且該結(jié)構(gòu)中間為鏤空結(jié)構(gòu),實現(xiàn)輕量化.壓電片采用PZT-5H 壓電陶瓷,經(jīng)極化處理后再與黃銅片粘在一起.

    圖1 驅(qū)動足的結(jié)構(gòu)Fig.1 The structure of the driving foot

    本文設(shè)計了一種基于逆壓電效應(yīng)的足腿一體化摩擦驅(qū)動結(jié)構(gòu),實驗中選用的壓電片參數(shù)如表1 所示,機(jī)器人基體為軸對稱結(jié)構(gòu),中間基板為半徑18 mm 的圓形鈹銅板,中間右側(cè)第2 段鈹銅基板長為5 mm,寬為12 mm,連接在圓形基底右側(cè),與腿部梁成135°夾角,第3 段鈹銅基板長為45 mm,寬為12 mm,且與足部連接板成135°夾角,第4 段鈹銅基板長為20 mm,寬為20 mm,用來連接驅(qū)動足,厚度均為0.2 mm.

    表1 材料參數(shù)Table 1 Material parameters

    其中該機(jī)器人總質(zhì)量為10.8 g,摩擦力驅(qū)動足通過雙面膠粘貼在足部連接板下端,借鑒已有的菱形、杠桿等位移放大機(jī)械裝置,將機(jī)器人結(jié)構(gòu)設(shè)計成近似梯形形狀使其頂部振動幅值得到放大.將該結(jié)構(gòu)看做一個兩足的運(yùn)動單元節(jié)以方便后續(xù)利用其來進(jìn)行多足機(jī)器人的制作.設(shè)計的兩足機(jī)器人實物模型參數(shù)如表2 所示,機(jī)器人基體展開俯視圖如圖2 所示,中間圓形基底為粘貼壓電片的部分.

    表2 機(jī)器人一個單元節(jié)的參數(shù)Table 2 Parameters of a robot unit

    圖2 機(jī)器人基體展開圖及尺寸 (單位: mm)Fig.2 Unfolded drawing and dimensions of the robot substrate (unit: mm)

    1.2 壓電機(jī)器人驅(qū)動原理

    在Comsol 中施加與兩足機(jī)器人結(jié)構(gòu)諧振頻率一致的電壓信號到壓電片后,它對基體產(chǎn)生壓電力,同時帶動基體產(chǎn)生相應(yīng)的應(yīng)力應(yīng)變,機(jī)器人的形變?nèi)鐖D3 所示.

    圖3 機(jī)器人在諧振頻率下的變形Fig.3 Deformation of the robot at a resonant frequency

    首先,分析機(jī)器人單側(cè)足腿的受力變形,在Comsol中設(shè)置壓電陶瓷片受到第一階振動模態(tài)的激勵頻率,對模型施加正弦波電壓激勵載荷,并對其一個周期內(nèi)進(jìn)行瞬態(tài)求解可得出機(jī)器人在一個運(yùn)動周期內(nèi)的形態(tài)變化.如圖4 所示,狀態(tài)(1)為模型上表面由平衡位置向上彎曲的狀態(tài),到正弦波峰值時機(jī)器人上表面達(dá)到最大彎曲狀態(tài),狀態(tài)(2)為彎曲達(dá)到最高點(diǎn)并即將向下彎曲的變形狀態(tài),狀態(tài)(3)為機(jī)器人由平衡位置開始向下彎曲達(dá)到最低點(diǎn),狀態(tài)(4)為模型運(yùn)動到最低點(diǎn)并開始向上運(yùn)動,最終再次回到平衡位置;在這個周期的瞬態(tài)解中可以發(fā)現(xiàn)機(jī)器人兩端的足部存在兩種運(yùn)動狀態(tài): 當(dāng)模型上基底向上彎曲時,機(jī)器人兩足端相互靠近;當(dāng)模型上基底向下彎曲時,機(jī)器人兩足端相互遠(yuǎn)離.

    圖4 結(jié)構(gòu)在正弦波激勵下的變形Fig.4 Structural deformation under sinusoidal wave excitation

    本文設(shè)計了一種基于逆壓電效應(yīng)的足腿一體化摩擦驅(qū)動結(jié)構(gòu),當(dāng)確定壓電片的極化方向后,在d31工作模式下(施加電壓方向與壓電片形變方向相垂直)的壓電陶瓷片在周期性電壓的驅(qū)動下會產(chǎn)生周期性彎曲,與之連接的彈性腿也將產(chǎn)生對應(yīng)的彎曲振動.當(dāng)驅(qū)動壓電片的激勵頻率達(dá)到機(jī)器人基體的共振頻率時,它的振動幅值會顯著增大,彈性腿連接的摩擦力驅(qū)動足與工作表面發(fā)生摩擦,將產(chǎn)生左右兩側(cè)不對稱的摩擦力來驅(qū)動結(jié)構(gòu)行走[30].當(dāng)施加連續(xù)正弦激勵和脈沖正弦激勵信號時,壓電機(jī)器人可分別實現(xiàn)較快速度和較高分辨力的運(yùn)動.基于以上原理設(shè)計出的機(jī)器人具有體積較小、傳動效率高、響應(yīng)速度快的優(yōu)點(diǎn).

    在機(jī)器人兩端足部設(shè)置不對稱的結(jié)構(gòu),可以引入左右方向不等的摩擦力,它可以在不對稱摩擦力的驅(qū)動下向摩擦力較小的一側(cè)運(yùn)動.將壓電片粘貼在兩足機(jī)器人單元體的圓形基底上,經(jīng)分析可得到具體的運(yùn)動過程如圖5 所示: 狀態(tài)(a)機(jī)器人的平衡位置,狀態(tài)(b) 為圓形基底處于最高點(diǎn)位置,狀態(tài)(c)為壓電片彎曲到最低點(diǎn)位置.當(dāng)壓電片由最高點(diǎn)向最低點(diǎn)運(yùn)動時,機(jī)器人兩端腿部相互遠(yuǎn)離,此時機(jī)器人右端足部受到正方向的大值摩擦力作用,左端受到負(fù)方向的小值摩擦力作用,受到的合力為正向力,驅(qū)使機(jī)器人向正方向運(yùn)動;當(dāng)壓電片由最低點(diǎn)向最高點(diǎn)運(yùn)動時,機(jī)器人兩端腿部相互靠近,此時機(jī)器人右端足部受到負(fù)方向的小值摩擦力作用,左端受到正方向的大值摩擦力作用,受到的合力驅(qū)動機(jī)器人繼續(xù)向正向運(yùn)動.

    圖5 移動機(jī)理Fig.5 The movement mechanism

    2 受力分析及動力學(xué)建模

    2.1 驅(qū)動腿的共振頻率分析

    如圖6 所示,利用材料力學(xué)知識[31]對驅(qū)動腿進(jìn)行理論分析,設(shè)橫向振動位移為w,軸向位移為u,則壓電作動器的應(yīng)力應(yīng)變關(guān)系為

    圖6 驅(qū)動腿的簡化模型Fig.6 Simplified model of the driving leg

    式中,z為縱坐標(biāo),S為梁上的應(yīng)變.

    將圓形壓電片和鈹銅基底簡化為一個質(zhì)量點(diǎn),下文將具體分析壓電片產(chǎn)生的激勵.將與其連接的機(jī)器人壓電腿部和足部連接板簡化為歐拉-伯努利梁[32],并對每段結(jié)構(gòu)進(jìn)行分析,w1,w2,w3分別代表3 段結(jié)構(gòu)的橫向位移,僅考慮彎曲振動和軸向振動,當(dāng)x∈(0,l1) 時

    當(dāng)x∈(0,l2) 時

    當(dāng)x∈(0,l3) 時

    式中,mk為各分段單位長度的質(zhì)量,g/mm;Ek為各分段等效楊氏模量,GPa;Ik為各分段的慣性矩,m4;ck為黏性阻尼系數(shù).

    式中,Wn為每段的寬度,m;tc為鈹銅梁厚度,m;ρc為梁體積密度,g/mm3.

    忽略黏性阻尼系數(shù)c時,設(shè)梁的橫向振動方程[33]的解為

    將式(14)帶入式(2)中,可得

    進(jìn)一步整理可得

    該方程左端為x的函數(shù),右端為t的函數(shù),且x與t相互獨(dú)立,因此上式應(yīng)該等于一個常數(shù),且該常數(shù)非負(fù),記為ω2≥ 0,因此上式可以分離為兩個獨(dú)立的微分方程

    由于各個分段的寬度和厚度也不同,因此彎曲梁的模態(tài)方程由3 個不同的部分組成.各段的橫向坐標(biāo)分別表示為x1,x2和x3.

    當(dāng)x1=0 時,邊界條件為

    當(dāng)xn=ln,xn+1=0,n=1,2 時,連續(xù)性邊界條件為

    當(dāng)x3=l3時,邊界條件為

    彎曲梁的橫向振型函數(shù)為

    設(shè)縱向振動方程的解為

    其中,An,Bn,Cn,Dn,En,Fn,Gn均為系數(shù)

    由式(14)~式(34)可以得出壓電驅(qū)動腿各階的固有頻率fn,驅(qū)動腿每段結(jié)構(gòu)簡化模型具有一致的共振頻率,可求得f1=242.52 Hz,f2=542.04 Hz.

    在Comsol 中建立機(jī)器人對稱單元節(jié)一半的結(jié)構(gòu)模型,中間半圓板橫截面處施加滾動鉸支座約束,有限元仿真分析的結(jié)果中可得到低階固有頻率為249.67 Hz,可以看到數(shù)值解和仿真結(jié)果相差并不大,說明分析結(jié)果較為準(zhǔn)確.

    2.2 圓形壓電片的等效激勵

    本文所選用的壓電片結(jié)構(gòu)如圖7 所示,下端為黃銅金屬層,厚度為t1;上端為PZT-5H 壓電陶瓷材料,厚度為t2;上層的壓電陶瓷焊接有導(dǎo)線接正極,下層的金屬接負(fù)極;以下端金屬層的中心位置為坐標(biāo)原點(diǎn)建立柱坐標(biāo)系,壓電陶瓷與金屬層為同圓心的層合圓板結(jié)構(gòu),壓電板的半徑為a,令z軸垂直平面向上,徑向坐標(biāo)r與環(huán)向坐標(biāo)θ服從右手定則.將壓電層的激勵等效為關(guān)節(jié)處的彎矩M(t),當(dāng)在壓電層上下兩端之間施加電壓時,由于d31效應(yīng)壓電片一端伸長另一端縮短,將產(chǎn)生彎曲變形,設(shè)金屬層在ro-z平面內(nèi)的轉(zhuǎn)角為φ(r,t),當(dāng)r=a時將其表示為φ(t).設(shè)粘貼有壓電片的黃銅板的彎曲剛度為K2,以產(chǎn)生等量轉(zhuǎn)角為原則有M(t)=K2φ(t).

    圖7 壓電片的柱坐標(biāo)系Fig.7 The cylindrical coordinate system of piezoelectric patch

    由以上基本假設(shè)可得徑向應(yīng)變、環(huán)向應(yīng)變與徑向坐標(biāo)、轉(zhuǎn)角之間的關(guān)系如下[34]

    式中,εr為徑向應(yīng)變,εθ為環(huán)向應(yīng)變,r為徑向坐標(biāo),φ為繞中心面的轉(zhuǎn)角.

    采用第3 類g 型壓電方程

    其中,σr為徑向應(yīng)力,σθ為環(huán)向應(yīng)力,Ez為z方向電場強(qiáng)度,D為電位移向量,g為壓電電壓常數(shù)矩陣,β為應(yīng)力不變條件下介質(zhì)隔離率矩陣,sD為恒電流密度條件下的柔性常數(shù)矩陣.

    本文選用的壓電材料為PZT-5H,其壓電介電常數(shù)矩陣、柔度系數(shù)矩陣、電位移向量以及壓電電壓常數(shù)矩陣參照標(biāo)定標(biāo)準(zhǔn)(z方向極化).其g 型壓電方程中的材料參數(shù)sD,g,β分別為

    將式(37)代入式(38),解得上壓電層的本構(gòu)方程[35]如下

    壓電層中r-o-z平面內(nèi)的力矩為

    將式(39)代入式(40)得

    由圓形薄板的軸對稱彎曲理論[35],可知黃銅金屬板中r-o-z平面內(nèi)的力矩為

    其中,M為壓電層中的力矩;Mm為金屬層中的彎矩;Em為黃銅的彈性模量;υ 為黃銅的泊松比.

    將Ez在z方向積分有

    其中,V為壓電層的上表面電壓值.

    將式(39)代入式(43)可得

    由力矩平衡關(guān)系得

    將上式代入可得

    2.3 對機(jī)器人右側(cè)腿的受力分析

    將計算得到的壓電層合板的激勵等效為關(guān)節(jié)處的彎矩M(t),如圖8 所示.ki為柔性足的等效關(guān)節(jié)剛度,zi為柔性足每段梁的等效質(zhì)量,φi為關(guān)節(jié)轉(zhuǎn)角,li為梁的長度,i=1,2,3.自由狀態(tài)下l1呈水平狀態(tài),l2與y軸方向的夾角為θ2.取x軸與地面平行,y軸正方向向上,建立直角坐標(biāo)系o-xy.在o-xy中,每個節(jié)點(diǎn)處的位置為(ai,bi),i=1,2,3,4,初始時刻a1與x軸上的o點(diǎn)重合,b4與y軸上的0 點(diǎn)重合.根據(jù)幾何關(guān)系,可得ai,bi與φi之間的關(guān)系如下

    圖8 一側(cè)腿部的受力分析Fig.8 Force analysis of one leg

    式中,θ2為初始狀態(tài)第2 段關(guān)節(jié)與y軸方向的夾角;θ3為彎曲后第3 段關(guān)節(jié)與x軸的夾角.

    根據(jù)上圖的結(jié)構(gòu)模型,可得系統(tǒng)的動力學(xué)方程式(59)~式(64)

    計算與地面未接觸時驅(qū)動足受激勵時y方向上的位移變化

    當(dāng)V=80sin(2πft) 時,y方向上的振幅為

    壓電片與圓板彎曲時,橫向位移為 Δa,隨時間變化的驅(qū)動足連接板末端x軸上的位置為

    未彎曲變形時,驅(qū)動足連接板末端x軸上的位置a4=74.42 mm,可計算x方向上的位移變化為

    將信號發(fā)生器頻率調(diào)為機(jī)器人共振頻率240 Hz,電壓為2 V,經(jīng)功率放大器放大40 倍后輸出到壓電片上,此時將機(jī)器人左側(cè)一半固定,右邊一半懸空放置,如圖9 所示.機(jī)器人腿部受到激勵帶動足端振動起來,由于達(dá)到了結(jié)構(gòu)的共振頻率,此時足端振幅最大,可采用合肥富煌君達(dá)高科信息技術(shù)有限公司的ISP502 千里狼高速攝像機(jī)(Revealer)記錄下機(jī)器人足部x方向和y方向隨時間變化的位移,設(shè)備如下圖10 所示.

    圖9 壓電力作用時驅(qū)動足的位移Fig.9 Displacement of the driving foot under the piezoelectric force

    圖10 高速攝像機(jī)拍攝足端振動Fig.10 High speed camera shooting foot end vibration

    由圖11 可以看到高速攝像機(jī)測得的y方向上的最大位移為0.299 62 mm,x方向上的最大位移為0.163 62 mm,而理論上的y方向和x方向上的振幅分別為 Δb=0.273 mm 和 Δa=0.128 mm.由于理論計算時忽略了連接到腿部連接板上驅(qū)動足的質(zhì)量,在實驗過程中,驅(qū)動足振動時的慣性會帶動腿部的搖擺,因此會導(dǎo)致機(jī)器人足端位移的增加.可以看出實驗測得的數(shù)據(jù)與理論計算所得的相差并不大,證明實驗結(jié)果較為可靠.

    圖11 正弦波激勵下機(jī)器人足部的振幅Fig.11 Amplitude of robot foot under sine wave excitation

    3 實驗驗證

    首先對機(jī)器人的單個運(yùn)動單元節(jié)進(jìn)行研究分析,實驗需要儀器如圖12 所示,通過信號發(fā)生器產(chǎn)生一定頻率和電壓的信號源,然后將此信號輸入功率放大器的輸入端,經(jīng)過放大后輸出到壓電陶瓷銅片使其彎曲振動;當(dāng)達(dá)到機(jī)器人的共振頻率時,此時機(jī)器人驅(qū)動足的位移變化最大,由于與地面接觸的足部前后兩個方向的摩擦力大小不同,機(jī)器人會按照摩擦力小的方向開始運(yùn)動.同時通過仿真和實驗測試機(jī)器人不同的模態(tài)頻率,當(dāng)達(dá)到機(jī)器人的其他模態(tài)頻率時,金屬基體連接驅(qū)動足會產(chǎn)生相應(yīng)的變形,從而帶動機(jī)器人與地面發(fā)生相對運(yùn)動,實現(xiàn)轉(zhuǎn)彎功能.

    圖12 機(jī)器人實驗流程圖Fig.12 Robot experiment flowchart

    將壓電片下端銅片焊接線連接負(fù)極且設(shè)置為零電壓,壓電陶瓷上方焊接線接正極并輸入周期性正弦電壓信號.為了驗證上述設(shè)計加工的壓電機(jī)器人工作性能和相關(guān)參數(shù),本文設(shè)計了相關(guān)實驗來研究機(jī)器人的頻率-速度關(guān)系、電壓-速度關(guān)系以及不同承載位置下的負(fù)載能力等.實驗中通過導(dǎo)軌來約束機(jī)器人單個運(yùn)動單元節(jié)的直行,驅(qū)動足及直線導(dǎo)軌均為3D 打印PLA 材料.

    通過實驗從大范圍來測試機(jī)器人在3 種不同足的驅(qū)動下速度較快的頻率點(diǎn),如圖13 所示,可以看到在200 Hz 和500 Hz 處機(jī)器人速度值較大,然后將信號源的電壓的頻率設(shè)置在共振頻率附近,驅(qū)動電壓幅值恒定為80 V,驅(qū)動頻率分別設(shè)置為190~270 Hz 和480~570 Hz.下面對不同類型的驅(qū)動足對在不同驅(qū)動頻率下的機(jī)器人運(yùn)動速度造成差異的原因進(jìn)行分析.由圖1 可得,斜面足與地面成一個45°的傾斜角度,相比于其他兩種驅(qū)動足,它與地面之間的摩擦力顯然更大,而且可以測得其所帶來的前后方向摩擦力的差異相較于圓弧面足更小,因此在斜面足驅(qū)動下機(jī)器人運(yùn)動速度最慢但最平穩(wěn),圓弧面與四圓弧體足與地面接觸部分在建模時均由一段圓弧線掃描拉伸而成,由于圓弧面足在回退時其內(nèi)表面更平滑,因此所帶來的與運(yùn)動方向相反的摩擦力較小,因此相同驅(qū)動頻率下運(yùn)動速度要比斜面足更快,四圓弧體足相比于其他兩結(jié)構(gòu)而言,其與地面具有更小的接觸面,在驅(qū)動腿帶動足部振動時,會受到摩擦、碰撞等多種工況作用,由于其與地面的接觸面積較小,因此其能量的損耗也較小,綜上可得,在圓弧體足驅(qū)動下機(jī)器人運(yùn)動速度最快,但運(yùn)動時穩(wěn)定性相較于其他兩足較差.

    圖13 3 種驅(qū)動足下機(jī)器人的頻率-速度圖Fig.13 Frequency-velocity diagram of three types of robots driving feet

    之前對比了3 種不同驅(qū)動足下的頻率-速度參數(shù),得到結(jié)果為四圓弧體足的驅(qū)動效果最好,黑色尼龍足的結(jié)構(gòu)形態(tài)也為圓弧體足,且尼龍的材料性能較 PLA 材料具有結(jié)構(gòu)更輕,韌性和抗壓強(qiáng)度更好,耐磨,表面更光滑且能夠吸收沖擊和振動等優(yōu)點(diǎn),為了形成與上述3 種結(jié)構(gòu)之間的對比,因此將機(jī)器人的驅(qū)動足更換為材料性能更好的黑色尼龍驅(qū)動足進(jìn)行實驗,并且在共振頻率附近處得到了在尼龍足驅(qū)動下的的機(jī)器人頻率-速度曲線,可以看出機(jī)器人的速度特性確有提升.圖14 為實驗得到的頻率-速度圖,可以看到速度曲線在240 Hz 附近存在著一個峰值,這個峰值即對應(yīng)著機(jī)器人驅(qū)動腿對應(yīng)的低階固有頻率;與此同時頻率-速度曲線在540 Hz 附近存在著另一個小的峰值,這個峰值即對應(yīng)著機(jī)器人的高階固有頻率.實驗曲線中的速度最高峰值出現(xiàn)在230~240 Hz 附近,這也與Comsol 當(dāng)中的機(jī)器人結(jié)構(gòu)固有模態(tài)的仿真結(jié)果能夠?qū)?yīng)上,在驅(qū)動頻率為240 Hz 附近的速度為92.9 mm/s,而驅(qū)動頻率為540 Hz 附近的速度接近60.2 mm/s.下面分析造成兩諧振點(diǎn)處的速度值差距較大的原因,通過高速攝像機(jī)拍攝機(jī)器人驅(qū)動足部在不同諧振點(diǎn)振動的位移,可以看到在240 Hz 低階固有頻率下驅(qū)動足水平方向振動幅值更大,與地面之間的作用力也會更大,所以移動速度會更快,而在高階諧振點(diǎn)時,足端在左右方向振動位移較小,因此速度較慢.

    在確定機(jī)器人固有振動頻率后,分別在設(shè)置為低階固有頻率和高階固有頻率的正弦波信號源的驅(qū)動下,來研究機(jī)器人移動速度與驅(qū)動電壓之間的關(guān)系,驅(qū)動頻率分別設(shè)置為240 Hz 和540 Hz.機(jī)器人的驅(qū)動足有3 種,分別為兩圓弧面、四圓弧體和黑色尼龍材料驅(qū)動足.由于材料加工和人工裝配誤差等原因會導(dǎo)致機(jī)器人結(jié)構(gòu)并不完全對稱,造成機(jī)器人無法在不加約束的條件下直行,因此加工了長度為40 cm 的導(dǎo)軌(3D 打印PLA 材料),記錄其在軌道上行走前20 cm 的時間,便可以計算得到機(jī)器人不同電壓下的移動速度.

    實驗結(jié)果如圖15 所示,可以看到,在240 Hz 驅(qū)動頻率下,機(jī)器人的速度隨電壓的增加而增加,在尼龍足的驅(qū)動下,機(jī)器人的速度遠(yuǎn)遠(yuǎn)高于其他兩足的速度,在110 V 時最快速度能達(dá)到14.2 cm/s,圓弧體足的速度高于弧面足,且隨著電壓的增加,機(jī)器人速度增加的趨勢變緩,這是由于電壓過大時,機(jī)器人足端振幅會過大,引起它在行進(jìn)過程中跳躍,與地面接觸不充分;在540 Hz 驅(qū)動頻率下,機(jī)器人的速度也是隨電壓的增加而增加,在尼龍足的驅(qū)動下機(jī)器人速度最快,其他兩足驅(qū)動下的速度差距不大,在實驗中可以觀察到在540 Hz 驅(qū)動頻率下,機(jī)器人雖然在速度上略小于240 Hz 時的速度,但就運(yùn)動平穩(wěn)性上來說,在高階固有頻率下機(jī)器人更占優(yōu),且速度隨電壓變化較緩,在低電壓驅(qū)動時,速度也不會過低,相比低階固有頻率具有一定優(yōu)勢;通過線性擬合可得到在尼龍足驅(qū)動下當(dāng)激勵信號頻率為240 Hz 時,機(jī)器人電壓-速度的一次方程為v1=0.208 8V-7.580 6,當(dāng)在540 Hz 激勵頻率下尼龍足驅(qū)動時,機(jī)器人電壓-速度的一次方程為v2=0.132 3V-3.659 1,其中驅(qū)動電壓V的單位為V,速度單位為cm/s.

    圖15 不同驅(qū)動足下機(jī)器人移動速度隨電壓的變化Fig.15 The variation of movement speed versus voltage for robots with different driving foots

    為了測試機(jī)器人在負(fù)載情況下的實驗速度,在機(jī)器人的頂部圓形基底和兩足端連接板處分別以圓形塑料片和鐵片來模擬機(jī)器人在實際中的負(fù)載,實驗照片如圖16 所示,實驗速度曲線如圖17 所示.驅(qū)動頻率分別設(shè)置為240 Hz 和540 Hz,施加在壓電片上的電壓為80 V,從圖17 中可以看出,隨著負(fù)載重量的增加,在大體的趨勢上機(jī)器人的速度是不斷下降的,當(dāng)負(fù)載在兩足端連接板上時,隨著塑料片數(shù)量的增加,在負(fù)重1.6 g 時,機(jī)器人速度會有一個微小的提升,這是因為機(jī)器人承受一定負(fù)重時,豎直方向上的跳躍會減少,減少能量損失的同時會增加機(jī)器人與接觸面的摩擦力,使得機(jī)器人的速度會增加一些;在更換負(fù)重為鐵片時,在x=3.8 g 機(jī)器人速度有一個階躍的提升是因為塑料片堆疊過高會導(dǎo)致機(jī)器人運(yùn)動不平穩(wěn);當(dāng)頂部圓形基底上的負(fù)重更換為鐵片時,機(jī)器人速度會增加也是由于類似的原因,當(dāng)堆疊過高會導(dǎo)致機(jī)器人運(yùn)動時重心不穩(wěn),對運(yùn)動時的干擾較大,負(fù)重更換為鐵片時就解決了這一問題;在頂部負(fù)載0.8 g 時,應(yīng)該是機(jī)器人在240 Hz 的驅(qū)動頻率下的最佳負(fù)載,此時不僅豎直方向上的位移會減少,與接觸面的摩擦力增大,而且在對比空載時的速度沒有較大下降的前提下會使機(jī)器人運(yùn)動的更加平穩(wěn).由此可以看出負(fù)載質(zhì)量對機(jī)器人運(yùn)動速度有較大影響,通過實驗驗證了機(jī)器人在240 Hz 的驅(qū)動頻率下單個單元節(jié)無負(fù)載條件下速度可達(dá)92.9 mm/s,當(dāng)在頂部加鐵片負(fù)載2.5 g 時,移動速度可達(dá)58.7 mm/s.雖然機(jī)器人運(yùn)動單元節(jié)最大負(fù)載能力有限,但在兩足端連接板上負(fù)載超過自身重量15.2 g 時速度也可以達(dá)到21.2 mm/s.因此,在實際應(yīng)用中機(jī)器人可以承擔(dān)一些運(yùn)輸小型物品的功能;此外,將負(fù)重更換為微型傳感器或攝像頭時,也可以為人類探測一些較危險的地帶,還能為機(jī)器人的實時運(yùn)動提供反饋,使機(jī)器人的行進(jìn)路線更準(zhǔn)確.

    圖16 在不同位置負(fù)載的機(jī)器人Fig.16 Robots loaded at different positions

    圖17 不同負(fù)載位置下機(jī)器人的速度Fig.17 The speed of robots under different load positions

    4 爬行機(jī)器人的控制

    4.1 不同構(gòu)型設(shè)計的爬行機(jī)器人

    在驅(qū)動原理確定的基礎(chǔ)上,可以將此兩足壓電結(jié)構(gòu)當(dāng)作一個模塊單元,通過不同的連接方式和構(gòu)型設(shè)計利用該單元節(jié)構(gòu)造出多種不同結(jié)構(gòu)構(gòu)型的多足爬行機(jī)器人,以上實驗研究主要是針對機(jī)器人的一個基礎(chǔ)運(yùn)動單元節(jié)來進(jìn)行的,實驗數(shù)據(jù)可以為之后更復(fù)雜構(gòu)型的機(jī)器人結(jié)構(gòu)提供參考,根據(jù)此結(jié)構(gòu)的布局變化可以令爬行機(jī)器人有多種形態(tài),如下圖18所示.

    圖18 壓電機(jī)器人的不同構(gòu)型Fig.18 Experimental equipment for piezoelectric robots

    之后本文設(shè)計制作了一種小型四足爬行機(jī)器人實物如圖19 所示,整體結(jié)構(gòu)由壓電陶瓷片、機(jī)器人基體本身和粘貼的摩擦力驅(qū)動足構(gòu)成.機(jī)器人基體結(jié)構(gòu)為一體化設(shè)計,由0.2 mm 厚的鈹銅板切割加工彎折而成,摩擦力驅(qū)動足被粘貼在足部連接板上,兩片圓形蜂鳴壓電片分別粘貼在機(jī)器人連接兩彈性腿的圓形基底的上表面.將驅(qū)動部件與傳動和執(zhí)行元件整合為一體,結(jié)構(gòu)設(shè)計簡單,傳動效率提高.

    圖19 壓電驅(qū)動的四足爬行機(jī)器人實物圖Fig.19 Manufactured prototype of the piezoelectric driven quadruped crawling robot

    4.2 四足爬行機(jī)器人直走和轉(zhuǎn)彎頻率的實驗驗證

    參考上面對兩足機(jī)器人的實驗研究,下文將對四足爬行機(jī)器人的性能參數(shù)進(jìn)行研究.通過對四足機(jī)器人整體進(jìn)行大量實驗時可以得出: 當(dāng)信號發(fā)生器產(chǎn)生頻率為560 Hz 電壓為2 V 的信號源,經(jīng)過功率放大器放大后輸出到壓電陶瓷銅片上,機(jī)器人可以直行且此時移動速度較快,而且得到在PLA 板上運(yùn)動時機(jī)器人的電壓-速度關(guān)系如下圖20 所示,通過線性擬合可得到四足機(jī)器人直線行走時的電壓-速度關(guān)系約為v=0.058V-0.483 9,其中電壓V的單位為V,速度單位為cm/s,電壓范圍為60~120 V;當(dāng)調(diào)節(jié)機(jī)器人頻率為520 Hz 時,可以帶動機(jī)器人向左轉(zhuǎn)彎,且轉(zhuǎn)彎順暢;當(dāng)信號發(fā)生器輸出頻率為940 Hz 左右電壓不變的信號源時,機(jī)器人可以實現(xiàn)右轉(zhuǎn)的功能.此外,當(dāng)信號發(fā)生器頻率為560 Hz時,也可以調(diào)節(jié)兩個單元節(jié)上壓電陶瓷片上的電壓,當(dāng)左側(cè)電壓較大時,機(jī)器人左側(cè)受到的驅(qū)動力更大,會使其向右轉(zhuǎn)動,同理可得,當(dāng)右側(cè)壓電片上的電壓調(diào)大時,機(jī)器人會向左轉(zhuǎn)彎.與此同時,可以在Simulink中搭建一定的程序,搭配不同頻率和電壓的信號源,以及一些控制輸入的模塊,將此程序連接到Quanser半物理仿真平臺;然后通過電壓放大到機(jī)器人壓電片上,即可通過更換不同的指令讓機(jī)器人進(jìn)行直行或轉(zhuǎn)彎等一系列動作,實驗裝置圖如圖21 所示.

    圖20 在560 Hz 時機(jī)器人的電壓-速度圖Fig.20 Voltage-velocity diagram of the robot at 560 Hz

    圖21 壓電機(jī)器人的實驗設(shè)備Fig.21 Experimental equipment for piezoelectric robots

    對機(jī)器人在560 Hz 驅(qū)動頻率下在不同接觸面上的電壓-速度進(jìn)行研究,圖22 為機(jī)器人在電壓60~110 V 區(qū)間內(nèi)在不同面上運(yùn)動時的速度曲線,可以看出隨著接觸面粗糙度的增加,同一驅(qū)動電壓下機(jī)器人的速度是逐漸降低的,從木板到辦公桌面隨著粗糙度的減小,運(yùn)動速度增加,但太過于光滑的表面也會導(dǎo)致機(jī)器人運(yùn)動時足部與接觸面打滑(比如玻璃板),使得運(yùn)動速度下降,因此當(dāng)機(jī)器人在適中的粗糙度接觸面(塑膠地)上的運(yùn)動速度最快.實驗研究表明,設(shè)計的四足機(jī)器人可以很好地運(yùn)用于各種環(huán)境以及不同粗糙度的接觸面.

    圖22 不同接觸面上的速度Fig.22 Velocity of the robot on different contact surfaces

    4.3 展望

    Simulink 是Matlab 中的可視化仿真工具,它支持系統(tǒng)仿真、程序設(shè)計、系統(tǒng)控制和信號處理等方面.Quanser 半物理仿真平臺[36]、與Simulink 接口板卡和Windows 在線控制程序Quarc 等,能夠直接將計算機(jī)與受控對象相連,構(gòu)成閉環(huán)控制結(jié)構(gòu).

    通過在Simulink 中編譯程序來對一些結(jié)構(gòu)進(jìn)行智能化控制,可以將該工具運(yùn)用到壓電機(jī)器人的控制當(dāng)中讓機(jī)器人實現(xiàn)預(yù)期化的運(yùn)動,從工程應(yīng)用的角度來說,不僅需要提供一個穩(wěn)定的電壓驅(qū)動壓電機(jī)器人運(yùn)動,還需要讓壓電片驅(qū)動機(jī)器人結(jié)構(gòu)在給定速度下進(jìn)行直線運(yùn)動,才能適應(yīng)各種實際需求.可以基于PID (proportional,integral,differential)控制器,在Simulink 中設(shè)計一種控制壓電驅(qū)動爬行機(jī)器人速度的模型,來使機(jī)器人實現(xiàn)給定速度下的穩(wěn)定運(yùn)動,然后連接到Quanser 半物理仿真平臺對壓電驅(qū)動四足爬行機(jī)器人的直行和轉(zhuǎn)彎運(yùn)動進(jìn)行控制實驗,在線控制程序Quarc 能夠啟動Simulink 模型生成的代碼,向板卡發(fā)送命令或從板卡輸出端收集反饋回的數(shù)據(jù),達(dá)到實時控制的目的,而且可以通過反饋的實時速度來改變程序中的參數(shù)使機(jī)器人實現(xiàn)預(yù)期下的運(yùn)動.分析對比實驗與理論的誤差,并進(jìn)一步優(yōu)化控制模型,為以后更精確控制壓電爬行機(jī)器人的速度模型提供參考.

    首先,對機(jī)器人單個模塊單元即兩足壓電結(jié)構(gòu)進(jìn)行驅(qū)動和控制,通過之前實驗測得的電壓-速度數(shù)據(jù)以及擬合方程,將其輸入到程序中建立兩足機(jī)器人速度與驅(qū)動電壓的關(guān)系模型,然后在程序中設(shè)置期望速度與行走時間,最后啟動仿真控制程序連接到Quanser 板卡上并通過激光位移傳感器傳輸回采集的數(shù)據(jù)從而實現(xiàn)反饋控制,圖23 為機(jī)器人在軌道上的行走.通過實驗可以發(fā)現(xiàn)兩足機(jī)器人在程序驅(qū)動控制下,行走并不穩(wěn)定且機(jī)身搖晃較劇烈,因此開始對四足機(jī)器人進(jìn)行實時控制實驗,將特定程序代碼輸出到板卡連接到功率放大器然后輸出到兩個壓電片上,研究結(jié)果發(fā)現(xiàn)驅(qū)動效果較好,比兩足結(jié)構(gòu)運(yùn)動更穩(wěn)定.圖24 為通過Simulink 搭建模型然后連接到板卡,設(shè)置機(jī)器人在程序啟動1 s 后開始運(yùn)動到8 s 時停止運(yùn)動,在此期間實現(xiàn)期望速度為5 cm/s 的直線運(yùn)動,通過激光位移傳感器實時采集得到位移數(shù)據(jù)微分后得到速度數(shù)據(jù),由于激光位移傳感器采集傳輸回的位移數(shù)據(jù)比較雜亂,干擾較多,因此在模型中加入濾波器模塊,將濾波器過濾掉雜波后的位移數(shù)據(jù)經(jīng)過限幅和微分后反饋到PID 模型中,讓程序自動調(diào)節(jié)輸出到功率放大器上的電壓參數(shù),使得機(jī)器人實現(xiàn)期望下的運(yùn)動,圖25 即為設(shè)置的目標(biāo)速度以及加上反饋后經(jīng)PID 調(diào)節(jié)后的擬合曲線,此時比例參數(shù)P值為2,積分參數(shù)I值為0.3,D值為1.

    圖23 兩種機(jī)器人在軌道上的運(yùn)動控制Fig.23 Motion control of two types of robots on orbits

    圖24 通過Simulink 搭建模型連接到板卡Fig.24 A Simulink model connected to the card

    圖25 期望速度與經(jīng)PID 調(diào)整后的速度曲線Fig.25 Expected speed and PID adjusted speed curve

    未來可以用Simscape 模塊來對機(jī)器人的運(yùn)動進(jìn)行路徑規(guī)劃,通過一系列數(shù)據(jù)點(diǎn)的設(shè)置,設(shè)定程序來計算出兩端驅(qū)動腿所需要達(dá)到的速度,進(jìn)而得出壓電片上所需的電壓,再通過功率放大器放大到壓電片上,從而使機(jī)器人完成規(guī)劃下的運(yùn)動.

    將本文的壓電機(jī)器人與近年來其他文獻(xiàn)中壓電機(jī)器人在重量、運(yùn)動速度、驅(qū)動器末端位移、最大負(fù)載、移動靈活性和可控性等方面的參數(shù)作對比,可以發(fā)現(xiàn)現(xiàn)有壓電機(jī)器人存在執(zhí)行機(jī)構(gòu)的驅(qū)動位移較小以及在復(fù)雜環(huán)境下移動與作業(yè)的能力不佳等問題,本文設(shè)計的機(jī)器人相比于現(xiàn)有的一些壓電機(jī)器人具有較快的移動速度和良好的移動靈活性,且具有較好的負(fù)載能力,可承載微型攝像頭和傳感器以勘測環(huán)境,且可通過仿真程序?qū)ζ溥M(jìn)行實時控制來完成未知環(huán)境下的作業(yè).

    5 結(jié)論

    本文設(shè)計了一種可以在各種平面上穩(wěn)定運(yùn)動的諧振式壓電爬行機(jī)器人,并利用振動力學(xué)和材料力學(xué)方法對該機(jī)器人建立了整體受力分析方程,得到了以下結(jié)論.

    (1)利用壓電片逆效應(yīng)設(shè)計出了一種新型的機(jī)器人結(jié)構(gòu),期望能按照預(yù)定方向去行走,其具有速度快、易控制等特點(diǎn).

    (2)通過Comsol 對機(jī)器人單元節(jié)的分析,得出了兩足機(jī)器人受到交變壓電力時的應(yīng)力-應(yīng)變和運(yùn)動原理,通過大量實驗分析總結(jié)出驅(qū)動電壓和機(jī)器人行進(jìn)速度的關(guān)系,并分析不同種類驅(qū)動足對機(jī)器人移動速度的影響,然后通過實驗驗證了機(jī)器人的低階和高階固有頻率與仿真結(jié)果較為吻合.此外,還對兩足機(jī)器人在不同位置負(fù)重時進(jìn)行了實驗,得出當(dāng)頂部負(fù)載2.5 g 時,機(jī)器人的峰值速度接近58.7 mm/s,在兩足端連接板上負(fù)載15.2 g 時速度也可以達(dá)到21.2 mm/s.證明該機(jī)器人在工程探測、搶險救援等方面具有潛在應(yīng)用價值.

    (3)當(dāng)通過Simulink 程序連接半物理仿真平臺Quancer 板卡利用不同頻率和幅值的驅(qū)動電壓控制四足爬行機(jī)器人時,實現(xiàn)了左轉(zhuǎn)、右轉(zhuǎn)和繞自身圓心自轉(zhuǎn)以及期望速度下的近似直線運(yùn)動.未來可以對機(jī)器人的運(yùn)動進(jìn)行路徑規(guī)劃,使機(jī)器人行進(jìn)過程更精確,更加完美地完成期望的路線.

    猜你喜歡
    壓電驅(qū)動頻率
    基于模糊PI控制的驅(qū)動防滑仿真系統(tǒng)分析
    屈宏斌:未來五年,雙輪驅(qū)動,砥礪前行
    振動與頻率
    軌旁ATC系統(tǒng)門控柜接收/驅(qū)動板改造
    《壓電與聲光》征稿啟事
    壓電與聲光(2019年1期)2019-02-22 09:46:06
    新型壓電疊堆泵設(shè)計及仿真
    極限頻率
    基于S3C6410的Wi-Fi驅(qū)動移植實現(xiàn)
    導(dǎo)航頻率源的同步與控制
    基于壓電激振的彈性模量測量方法
    成熟少妇高潮喷水视频| 又粗又硬又长又爽又黄的视频 | 亚洲国产精品久久男人天堂| 久久久久久久久久成人| 日韩精品青青久久久久久| 又粗又硬又长又爽又黄的视频 | 麻豆成人av视频| 一进一出抽搐gif免费好疼| 婷婷六月久久综合丁香| 亚洲婷婷狠狠爱综合网| 天堂影院成人在线观看| 免费无遮挡裸体视频| 日韩人妻高清精品专区| 色综合站精品国产| 久久精品国产亚洲网站| 国产精品人妻久久久久久| 两个人的视频大全免费| 亚洲内射少妇av| 毛片女人毛片| 国产探花在线观看一区二区| 天天躁夜夜躁狠狠久久av| 亚洲国产欧洲综合997久久,| 国产成人aa在线观看| 成人美女网站在线观看视频| ponron亚洲| 国产国拍精品亚洲av在线观看| 毛片一级片免费看久久久久| 禁无遮挡网站| 美女高潮的动态| 久久人妻av系列| 特级一级黄色大片| 99久久久亚洲精品蜜臀av| 久久这里只有精品中国| 亚洲18禁久久av| 国产国拍精品亚洲av在线观看| 91久久精品国产一区二区三区| 欧美日韩乱码在线| 亚洲va在线va天堂va国产| 中国美女看黄片| 亚洲国产色片| 中文字幕精品亚洲无线码一区| 国产日本99.免费观看| 晚上一个人看的免费电影| 男女那种视频在线观看| 少妇的逼水好多| 色综合色国产| 成人毛片a级毛片在线播放| 欧美日韩在线观看h| 爱豆传媒免费全集在线观看| 欧美性猛交╳xxx乱大交人| 男插女下体视频免费在线播放| a级毛色黄片| 级片在线观看| 人人妻人人澡人人爽人人夜夜 | 欧美日韩国产亚洲二区| 日韩中字成人| 最近手机中文字幕大全| 成人午夜精彩视频在线观看| 亚洲av中文av极速乱| 日韩强制内射视频| 亚洲av不卡在线观看| 国产精品一区二区三区四区免费观看| 久久久久久久久大av| 性色avwww在线观看| 老女人水多毛片| 在线天堂最新版资源| 日日撸夜夜添| 欧美色视频一区免费| 午夜精品一区二区三区免费看| 日本黄色视频三级网站网址| a级毛片免费高清观看在线播放| 欧美色视频一区免费| 一级毛片aaaaaa免费看小| 欧美一区二区精品小视频在线| 哪里可以看免费的av片| 国产精品一二三区在线看| 国产伦在线观看视频一区| 日产精品乱码卡一卡2卡三| 美女xxoo啪啪120秒动态图| 一级毛片久久久久久久久女| 直男gayav资源| 成人亚洲欧美一区二区av| www.色视频.com| 一边摸一边抽搐一进一小说| 天堂av国产一区二区熟女人妻| 久久热精品热| 青春草亚洲视频在线观看| 亚洲一区高清亚洲精品| 日本熟妇午夜| 美女黄网站色视频| 中文字幕精品亚洲无线码一区| 日韩成人伦理影院| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久久久久av不卡| 精品国内亚洲2022精品成人| 深夜a级毛片| 亚洲国产精品成人综合色| 国内久久婷婷六月综合欲色啪| 国产老妇女一区| 深夜精品福利| 桃色一区二区三区在线观看| 久久久久久久久久久免费av| 国产国拍精品亚洲av在线观看| 最近2019中文字幕mv第一页| 国内精品一区二区在线观看| 亚洲国产欧洲综合997久久,| 啦啦啦啦在线视频资源| 人人妻人人澡人人爽人人夜夜 | 国产极品精品免费视频能看的| 波多野结衣高清作品| 狂野欧美白嫩少妇大欣赏| 日韩成人伦理影院| 久久久久久久久大av| 久久午夜亚洲精品久久| 天堂av国产一区二区熟女人妻| av在线观看视频网站免费| 69人妻影院| 国内久久婷婷六月综合欲色啪| 久久精品91蜜桃| 国产一区二区激情短视频| 99热这里只有是精品在线观看| www.av在线官网国产| 色综合亚洲欧美另类图片| 亚洲欧美清纯卡通| 久久亚洲国产成人精品v| 免费观看a级毛片全部| 亚洲婷婷狠狠爱综合网| 国产午夜精品一二区理论片| 午夜a级毛片| 18禁裸乳无遮挡免费网站照片| 国产在线精品亚洲第一网站| 一区二区三区高清视频在线| 精品免费久久久久久久清纯| 99久久成人亚洲精品观看| 午夜福利在线观看免费完整高清在 | 啦啦啦啦在线视频资源| 女同久久另类99精品国产91| 少妇被粗大猛烈的视频| 精品久久久久久久久av| a级毛片a级免费在线| 超碰av人人做人人爽久久| 国产美女午夜福利| 婷婷亚洲欧美| 国产一区二区亚洲精品在线观看| 蜜桃久久精品国产亚洲av| 久久精品91蜜桃| a级毛片a级免费在线| 国产精品久久久久久久久免| 欧美三级亚洲精品| 日本爱情动作片www.在线观看| 尤物成人国产欧美一区二区三区| 国产亚洲5aaaaa淫片| 欧美zozozo另类| 国产大屁股一区二区在线视频| 成人毛片60女人毛片免费| 99九九线精品视频在线观看视频| 国产又黄又爽又无遮挡在线| 国产精品无大码| videossex国产| 可以在线观看毛片的网站| 国产一区二区三区av在线 | 22中文网久久字幕| 久久韩国三级中文字幕| 成人午夜精彩视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲一区高清亚洲精品| 国内精品美女久久久久久| 听说在线观看完整版免费高清| 中文字幕制服av| 成人特级黄色片久久久久久久| 丰满乱子伦码专区| 中文在线观看免费www的网站| 在线观看av片永久免费下载| 毛片一级片免费看久久久久| 国产激情偷乱视频一区二区| 婷婷六月久久综合丁香| 国产视频内射| 身体一侧抽搐| 国产黄色小视频在线观看| 一本精品99久久精品77| 国产 一区精品| 一级毛片我不卡| 看片在线看免费视频| 成人三级黄色视频| 亚洲国产精品成人久久小说 | 男人的好看免费观看在线视频| 成人一区二区视频在线观看| 久久精品国产亚洲av香蕉五月| 熟女人妻精品中文字幕| 网址你懂的国产日韩在线| 亚洲内射少妇av| 国产一级毛片在线| 不卡一级毛片| 99久久无色码亚洲精品果冻| 国产伦理片在线播放av一区 | 精品日产1卡2卡| 免费观看在线日韩| 亚洲国产精品成人综合色| 联通29元200g的流量卡| 男女下面进入的视频免费午夜| 国产片特级美女逼逼视频| 精品少妇黑人巨大在线播放 | 男人和女人高潮做爰伦理| 99热精品在线国产| 高清日韩中文字幕在线| 欧美成人免费av一区二区三区| 日本黄大片高清| 嫩草影院新地址| 国产在线精品亚洲第一网站| 亚洲欧美成人综合另类久久久 | kizo精华| 国产成人精品婷婷| 国产在视频线在精品| 2022亚洲国产成人精品| 精品一区二区三区视频在线| 在线a可以看的网站| 欧美性猛交黑人性爽| 亚洲久久久久久中文字幕| 中文字幕av成人在线电影| 亚洲欧美日韩卡通动漫| 中文字幕免费在线视频6| 久久久久免费精品人妻一区二区| 国产精品福利在线免费观看| 天堂影院成人在线观看| 91av网一区二区| 1024手机看黄色片| 欧美不卡视频在线免费观看| 我要搜黄色片| 在线国产一区二区在线| 黑人高潮一二区| 久久草成人影院| 99久久久亚洲精品蜜臀av| 免费人成视频x8x8入口观看| 亚洲,欧美,日韩| 国产成人精品婷婷| 老师上课跳d突然被开到最大视频| 日本与韩国留学比较| 国内精品宾馆在线| 国产精品人妻久久久久久| 精品一区二区三区人妻视频| 校园春色视频在线观看| 亚洲人成网站在线观看播放| 午夜福利在线观看免费完整高清在 | 国产精品久久久久久亚洲av鲁大| 日韩欧美三级三区| 亚洲av免费高清在线观看| 国内精品久久久久精免费| 12—13女人毛片做爰片一| 精品久久国产蜜桃| 欧美xxxx性猛交bbbb| 97超碰精品成人国产| 村上凉子中文字幕在线| 国产黄片视频在线免费观看| 国产久久久一区二区三区| 亚洲电影在线观看av| 五月玫瑰六月丁香| 国产精品永久免费网站| 五月伊人婷婷丁香| 美女 人体艺术 gogo| 全区人妻精品视频| 国产精品久久视频播放| 大香蕉久久网| 99九九线精品视频在线观看视频| 看黄色毛片网站| 99久久中文字幕三级久久日本| 国产精品伦人一区二区| 亚洲乱码一区二区免费版| 成人性生交大片免费视频hd| 人体艺术视频欧美日本| 精品久久久久久久末码| 国产精品.久久久| 听说在线观看完整版免费高清| 国内精品久久久久精免费| 男人舔女人下体高潮全视频| 极品教师在线视频| 91aial.com中文字幕在线观看| 精品一区二区三区视频在线| 亚洲一区二区三区色噜噜| 性插视频无遮挡在线免费观看| 国产精品不卡视频一区二区| 国产精品爽爽va在线观看网站| av在线老鸭窝| 美女黄网站色视频| 国产亚洲av片在线观看秒播厂 | 天天一区二区日本电影三级| 深夜精品福利| 午夜激情福利司机影院| 欧洲精品卡2卡3卡4卡5卡区| 老司机影院成人| 99久久成人亚洲精品观看| 全区人妻精品视频| 97在线视频观看| 嘟嘟电影网在线观看| 夫妻性生交免费视频一级片| 久久人人爽人人爽人人片va| 久久久国产成人精品二区| 网址你懂的国产日韩在线| 美女黄网站色视频| 亚洲人与动物交配视频| 免费看光身美女| 亚洲国产欧美在线一区| 国产精品一二三区在线看| 一个人看的www免费观看视频| 爱豆传媒免费全集在线观看| 高清毛片免费看| 亚洲精品日韩av片在线观看| 亚洲精品456在线播放app| 日日啪夜夜撸| 久久久久性生活片| 此物有八面人人有两片| 久久精品夜色国产| 波多野结衣高清无吗| 一本一本综合久久| 麻豆久久精品国产亚洲av| 性色avwww在线观看| 国语自产精品视频在线第100页| 一夜夜www| 中国美白少妇内射xxxbb| 麻豆国产97在线/欧美| 国产在线精品亚洲第一网站| 国产一区亚洲一区在线观看| 亚洲在线观看片| 99热网站在线观看| 在线播放国产精品三级| 丰满的人妻完整版| 在线播放国产精品三级| 日韩制服骚丝袜av| av黄色大香蕉| 久久精品国产清高在天天线| 丝袜喷水一区| 国产成年人精品一区二区| 亚洲天堂国产精品一区在线| 日韩欧美一区二区三区在线观看| 日本在线视频免费播放| 久久婷婷人人爽人人干人人爱| 免费看日本二区| 亚洲乱码一区二区免费版| 韩国av在线不卡| 22中文网久久字幕| 国产中年淑女户外野战色| 精品人妻偷拍中文字幕| 美女 人体艺术 gogo| 国产乱人视频| 午夜激情欧美在线| av在线观看视频网站免费| 秋霞在线观看毛片| 国产av在哪里看| 久久人人爽人人片av| 亚洲激情五月婷婷啪啪| 乱系列少妇在线播放| 欧美一区二区精品小视频在线| 久久久国产成人免费| 国产成人精品久久久久久| 美女cb高潮喷水在线观看| 中文亚洲av片在线观看爽| 久久午夜福利片| 亚洲经典国产精华液单| 日韩高清综合在线| 亚洲一区二区三区色噜噜| 日本欧美国产在线视频| 成人毛片60女人毛片免费| 好男人在线观看高清免费视频| 寂寞人妻少妇视频99o| 99久久精品国产国产毛片| 在线天堂最新版资源| 午夜福利在线观看吧| av天堂中文字幕网| 国产亚洲av嫩草精品影院| 熟女电影av网| 午夜精品在线福利| 精品人妻偷拍中文字幕| 青春草国产在线视频 | 99国产精品一区二区蜜桃av| 免费av不卡在线播放| 午夜福利在线观看吧| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久末码| 成年女人看的毛片在线观看| 成人av在线播放网站| 12—13女人毛片做爰片一| 成人av在线播放网站| 日韩人妻高清精品专区| 免费在线观看成人毛片| 村上凉子中文字幕在线| 欧美另类亚洲清纯唯美| 国产精品一及| 联通29元200g的流量卡| 中文字幕熟女人妻在线| 亚洲最大成人中文| 精品久久久久久久末码| 国产极品精品免费视频能看的| 男人和女人高潮做爰伦理| 我的老师免费观看完整版| 久99久视频精品免费| 又粗又爽又猛毛片免费看| 久久精品国产亚洲av涩爱 | 在线天堂最新版资源| 免费观看的影片在线观看| 国产精品乱码一区二三区的特点| 成年女人永久免费观看视频| 97超碰精品成人国产| 97在线视频观看| 国产又黄又爽又无遮挡在线| 久久久成人免费电影| 精品人妻一区二区三区麻豆| 长腿黑丝高跟| 看黄色毛片网站| 美女xxoo啪啪120秒动态图| 不卡视频在线观看欧美| 国模一区二区三区四区视频| 人妻夜夜爽99麻豆av| 国产美女午夜福利| 91aial.com中文字幕在线观看| 91精品国产九色| 午夜a级毛片| 午夜久久久久精精品| 精品午夜福利在线看| av专区在线播放| 欧美性猛交黑人性爽| 国产午夜精品论理片| 一区二区三区四区激情视频 | 国产成人aa在线观看| 亚洲精品久久久久久婷婷小说 | 国产在线精品亚洲第一网站| 亚洲国产欧美人成| 成人永久免费在线观看视频| 99久久无色码亚洲精品果冻| 永久网站在线| 韩国av在线不卡| 精华霜和精华液先用哪个| 国产久久久一区二区三区| 97超视频在线观看视频| 赤兔流量卡办理| 亚洲av一区综合| 人妻制服诱惑在线中文字幕| 成人无遮挡网站| 国产一级毛片七仙女欲春2| 国产精品美女特级片免费视频播放器| 日本与韩国留学比较| 午夜精品国产一区二区电影 | 蜜桃久久精品国产亚洲av| 91麻豆精品激情在线观看国产| 国产色婷婷99| 欧美性猛交╳xxx乱大交人| 黄色配什么色好看| 禁无遮挡网站| av天堂在线播放| 亚洲av电影不卡..在线观看| 成人亚洲精品av一区二区| 亚洲欧洲国产日韩| 欧美一级a爱片免费观看看| 丰满乱子伦码专区| 特级一级黄色大片| 一级毛片我不卡| 国产高潮美女av| 青春草国产在线视频 | 国产亚洲欧美98| 成人国产麻豆网| 亚洲人成网站高清观看| 乱人视频在线观看| 波野结衣二区三区在线| 人妻少妇偷人精品九色| av免费观看日本| 国产亚洲av片在线观看秒播厂 | 高清毛片免费观看视频网站| 人妻久久中文字幕网| 一区二区三区高清视频在线| 高清毛片免费观看视频网站| 人妻制服诱惑在线中文字幕| 久久久精品欧美日韩精品| av黄色大香蕉| 热99在线观看视频| 国产精品久久电影中文字幕| 婷婷六月久久综合丁香| 尤物成人国产欧美一区二区三区| 久久精品久久久久久噜噜老黄 | 成人鲁丝片一二三区免费| 久久欧美精品欧美久久欧美| 亚洲人与动物交配视频| 1000部很黄的大片| 男人和女人高潮做爰伦理| 国产成人精品久久久久久| 成人性生交大片免费视频hd| 色哟哟哟哟哟哟| 插阴视频在线观看视频| 亚洲一区二区三区色噜噜| 噜噜噜噜噜久久久久久91| 亚洲成av人片在线播放无| 少妇裸体淫交视频免费看高清| 天天一区二区日本电影三级| 91久久精品国产一区二区成人| 在线观看66精品国产| 免费观看的影片在线观看| 我的老师免费观看完整版| 国产精品一区www在线观看| 天堂中文最新版在线下载 | 女人被狂操c到高潮| 日韩欧美精品v在线| 色视频www国产| 亚洲乱码一区二区免费版| 成熟少妇高潮喷水视频| 少妇丰满av| 欧美成人精品欧美一级黄| 一区福利在线观看| 中文字幕制服av| 日日撸夜夜添| 最近最新中文字幕大全电影3| 精品一区二区免费观看| 亚洲人成网站高清观看| 变态另类丝袜制服| 日韩欧美在线乱码| 偷拍熟女少妇极品色| 国产伦在线观看视频一区| 国产高清有码在线观看视频| 精品不卡国产一区二区三区| 18禁在线无遮挡免费观看视频| 国产成人一区二区在线| 91久久精品国产一区二区成人| 女同久久另类99精品国产91| 久久久国产成人免费| 国产三级中文精品| 欧美色欧美亚洲另类二区| 亚洲国产精品sss在线观看| 中文字幕免费在线视频6| 91麻豆精品激情在线观看国产| 亚洲在线自拍视频| 久久精品综合一区二区三区| av在线亚洲专区| 两个人视频免费观看高清| 人妻制服诱惑在线中文字幕| 欧美一区二区国产精品久久精品| 99久久精品一区二区三区| 99在线人妻在线中文字幕| 免费观看精品视频网站| 国产黄a三级三级三级人| 久久久久久久久中文| 国产精品三级大全| 亚洲在线观看片| 国产成年人精品一区二区| 国产又黄又爽又无遮挡在线| 丰满的人妻完整版| 久久精品国产鲁丝片午夜精品| 精品国产三级普通话版| 精品久久久久久成人av| 色噜噜av男人的天堂激情| 女的被弄到高潮叫床怎么办| 97人妻精品一区二区三区麻豆| 日本黄色片子视频| 免费看a级黄色片| 麻豆精品久久久久久蜜桃| 久久精品人妻少妇| 热99re8久久精品国产| 91午夜精品亚洲一区二区三区| 国产真实乱freesex| 99久久久亚洲精品蜜臀av| 麻豆成人午夜福利视频| 一区二区三区高清视频在线| 国产人妻一区二区三区在| 99久久久亚洲精品蜜臀av| 国产片特级美女逼逼视频| 91久久精品电影网| 国产黄色视频一区二区在线观看 | 国产精品一区二区性色av| 亚洲欧美日韩高清专用| 成人毛片60女人毛片免费| 久久九九热精品免费| 91狼人影院| 国产视频首页在线观看| 亚洲国产日韩欧美精品在线观看| 国产一区二区亚洲精品在线观看| av在线亚洲专区| 97热精品久久久久久| 村上凉子中文字幕在线| 青春草视频在线免费观看| 一个人观看的视频www高清免费观看| 午夜福利高清视频| 亚洲不卡免费看| 国产一区二区三区在线臀色熟女| 国产老妇女一区| 一级av片app| 九色成人免费人妻av| 老熟妇乱子伦视频在线观看| 一级黄片播放器| 国产成年人精品一区二区| 大香蕉久久网| 久久精品国产亚洲网站| av又黄又爽大尺度在线免费看 | 成人毛片60女人毛片免费| 亚洲美女搞黄在线观看| 欧美性感艳星| 12—13女人毛片做爰片一| 国产av一区在线观看免费| 超碰av人人做人人爽久久| 日韩欧美 国产精品| 国产av一区在线观看免费| 免费观看人在逋| 日本免费一区二区三区高清不卡| 日本在线视频免费播放| 少妇裸体淫交视频免费看高清| 一进一出抽搐gif免费好疼| 国产av一区在线观看免费| 午夜亚洲福利在线播放| 久久热精品热| 久久人妻av系列| 亚洲国产精品成人综合色| 欧美又色又爽又黄视频| 亚洲欧美日韩高清在线视频| 日日干狠狠操夜夜爽| 午夜福利在线在线| 日韩视频在线欧美| 中国美女看黄片| 18禁黄网站禁片免费观看直播| 看黄色毛片网站| 精品欧美国产一区二区三| 搡女人真爽免费视频火全软件|