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

    羊舍自主巡檢機(jī)器人設(shè)計(jì)與試驗(yàn)

    2025-01-19 00:00:00趙寧磊程曼袁洪波王媛蔡振江

    摘要:

    針對規(guī)模化羊舍人工勞動強(qiáng)度大、作業(yè)效率低等問題,設(shè)計(jì)羊舍自主巡檢機(jī)器人系統(tǒng)。該系統(tǒng)基于2D激光雷達(dá)的即時定位與建圖(SLAM)技術(shù),自主導(dǎo)航系統(tǒng)使用基于擴(kuò)展卡爾曼濾波數(shù)據(jù)融合的AMCL定位算法實(shí)現(xiàn)機(jī)器人的定位,采用Gmapping算法構(gòu)建羊舍環(huán)境地圖,并采用A*和DWA算法規(guī)劃機(jī)器人巡檢路徑,最終實(shí)現(xiàn)舍內(nèi)地圖構(gòu)建、系統(tǒng)實(shí)時定位、多目標(biāo)點(diǎn)路徑規(guī)劃以及動態(tài)避障。機(jī)器人導(dǎo)航精度試驗(yàn)表明,巡檢機(jī)器人橫向偏差均值不超過8.3cm,縱向偏差均值不超過7.1cm;航向偏差的均值不超過6.1°。該系統(tǒng)構(gòu)建的環(huán)境地圖精度較高,能夠?qū)崿F(xiàn)對機(jī)器人的有效定位,并且能夠以實(shí)現(xiàn)羊舍的循環(huán)巡航和自動躲避周邊障礙物,導(dǎo)航精度滿足羊舍自主導(dǎo)航需求。

    關(guān)鍵詞:羊舍;巡檢機(jī)器人;多傳感器融合;路徑規(guī)劃;擴(kuò)展卡爾曼濾波

    中圖分類號:S24

    文獻(xiàn)標(biāo)識碼:A

    文章編號:2095-5553 (2025) 01-0078-07

    Design and experiment of autonomous inspection robot for sheep house

    Zhao Ninglei1, Cheng Man1, Yuan Hongbo1, Wang Yuan2, Cai Zhenjiang1

    (1. College of Mechanical and Electrical Engineering, Hebei Agricultural University, Baoding, 071001, China;

    2. Hengshui Zhihao Animal Husbandry Technology Co., Ltd., Hengshui, 053400, China)

    Abstract:

    In order to solve the problems of high labor intensity and low working efficiency of large-scale sheep house, an autonomous inspection robot system for sheep house is designed. The system is based on the real-time positioning and mapping (SLAM) technology of 2D LiDAR. The autonomous navigation system uses the AMCL positioning algorithm based on extended Kalman filter data fusion to realize the positioning of the robot, uses the Gmapping algorithm to construct the environment map of the sheep shed, and uses the A*"and DWA algorithms to plan the inspection path of the robot. Finally, it realizes the construction of the map in the house, real-time positioning of the system, multi-objective point path planning and dynamic obstacle avoidance. The robot navigation accuracy test shows that the average lateral deviation of the patrol robot does not exceed 8.3cm, the average longitudinal deviation does not exceed 7.1cm, and the average heading deviation does not exceed 6.1°. The environmental map constructed by this system has high accuracy, can effectively locate the robot, and can realize the cyclic cruising and cruising of the sheep house.

    Keywords:

    sheep house; inspection robot; multi-sensor fusion; path planning; extended kalman filtering

    0"引言

    近年來,中國畜禽養(yǎng)殖業(yè)規(guī)模和產(chǎn)值不斷增加,2022年總產(chǎn)值達(dá)到4.1萬億元,約占農(nóng)業(yè)總產(chǎn)值的30%[1],且畜禽產(chǎn)品種類、總量和產(chǎn)值等均距世界前列。然而,隨著畜禽養(yǎng)殖集約化程度不斷提高,高密度養(yǎng)殖對畜禽舍內(nèi)動物的健康狀況形成巨大的挑戰(zhàn)[2]。近年來,隨著養(yǎng)殖模式集約化和標(biāo)準(zhǔn)化程度的不斷提高,研發(fā)和應(yīng)用全自動移動式機(jī)器人,可有效促進(jìn)養(yǎng)殖舍內(nèi)智能化高效管理[3]。傳統(tǒng)的羊場巡檢方式中,人工巡檢方式不僅存在工作強(qiáng)度大、作業(yè)效率低等不足[4],且羊自身的傳染病還有很大幾率傳染給巡檢人員,存在較大的人身安全隱患問題[5]。傳統(tǒng)的視頻監(jiān)控系統(tǒng),由于存在較大的監(jiān)控盲區(qū)很難做到視頻監(jiān)控全方位覆蓋,同時由于整個系統(tǒng)使用攝像頭數(shù)量多、安裝布線復(fù)雜、綜合系統(tǒng)龐大,因此故障率較高,維護(hù)成本大[6]。使用智能移動機(jī)器人進(jìn)行羊場環(huán)境的巡檢,不僅可以為企業(yè)降低生產(chǎn)成本、提高生產(chǎn)效率、擴(kuò)大企業(yè)的行業(yè)競爭力[7],也可以有效避免人工巡檢及傳統(tǒng)視頻監(jiān)控過程中存在的缺陷和不足。

    針對羊場羊只巡檢的需要,設(shè)計(jì)一種上下位機(jī)結(jié)合的智能巡檢機(jī)器人,智能巡檢機(jī)器人采用PID控制的四輪差速移動底盤,保證其靈活性以及足夠的檢測范圍。應(yīng)用Gmapping算法并開發(fā)上位機(jī)界面,實(shí)現(xiàn)機(jī)器人的自主構(gòu)建地圖和實(shí)時獲取和存儲巡檢目標(biāo)點(diǎn)的功能。同時利用A*和DWA路徑規(guī)劃算法對每一個記錄的目標(biāo)點(diǎn)依次進(jìn)行路徑規(guī)劃和自主導(dǎo)航從而滿足羊場巡檢的需要。

    1"巡檢機(jī)器人系統(tǒng)設(shè)計(jì)

    1.1"巡檢機(jī)器人結(jié)構(gòu)設(shè)計(jì)

    機(jī)器人整體框架由40mm×40mm的歐標(biāo)鋁型材搭建,整車長為55cm、寬為55cm、高為65cm。機(jī)器人在室內(nèi)巡檢時的數(shù)據(jù)采集場景如圖1所示。整車采用四輪四驅(qū)獨(dú)立轉(zhuǎn)向機(jī)構(gòu),由4組直角直流減速電機(jī)、4組驅(qū)動器、4組AGV充氣輪組成,并通過4組長為255mm,寬為84mm,厚度為5mm的自制電機(jī)架將電機(jī)固定在鋁型材上,系統(tǒng)所選用的電機(jī)安裝于內(nèi)側(cè),電機(jī)與電機(jī)架的整體寬度約為95mm。機(jī)器人正前方中央放置一個RPLIDAR S1單線激光雷達(dá),可以獲取其探測范圍內(nèi)雷達(dá)平面范圍內(nèi)的二維點(diǎn)云數(shù)據(jù)。底盤中央安置了一組24V 40Ah的鋰電池,是所有傳感器、動力裝置的動力來源,能保證平臺在設(shè)備滿功率狀態(tài)下續(xù)航4h。主控器和從控制器單元放置在電源旁邊的盒子里。

    1.2"巡檢機(jī)器人硬件設(shè)計(jì)與選型

    機(jī)器人主要硬件組成如表1所示,該機(jī)器人硬件主要由傳感器模塊、控制模塊和驅(qū)動模塊組成。機(jī)器人系統(tǒng)的硬件框架如圖2所示。Jetson Nano主控制器通過串口驅(qū)動激光雷達(dá)獲取平面點(diǎn)云信息,STM32通過獲取編碼器信息,計(jì)算得出機(jī)器人的線速度與角速度,并通過PWM驅(qū)動電機(jī)控制板,進(jìn)而控制電機(jī)的速度,從而控制機(jī)器人的移動。此外,還通過模擬I2C協(xié)議獲取JY901的陀螺儀和加速度信息,將機(jī)器人的速度、陀螺儀和加速度等信息通過串口上傳到主控制器。主控制器通過獲取從控制器發(fā)過來的速度信息計(jì)算得出機(jī)器人的里程計(jì)、陀螺儀和加速度信息。

    2"巡檢功能設(shè)計(jì)

    2.1"基于擴(kuò)展卡爾曼濾波數(shù)據(jù)融合的AMCL定位算法

    機(jī)器人定位采用擴(kuò)展卡爾曼濾波數(shù)據(jù)融合的AMCL算法,AMCL定位算法基于粒子濾波,確定機(jī)器人當(dāng)前的位置[8]。首先,該算法將代表機(jī)器人位姿的粒子隨機(jī)分布于已知環(huán)境地圖中。其次,通過使用擴(kuò)展卡爾曼濾波融合IMU傳感器與里程計(jì)數(shù)據(jù),建立運(yùn)動模型預(yù)測所有隨機(jī)粒子的位姿。擴(kuò)展卡爾曼濾波算法通過系統(tǒng)的觀測數(shù)據(jù)及系統(tǒng)狀態(tài)對全局最優(yōu)狀態(tài)進(jìn)行推算,并通過上一狀態(tài)和下一狀態(tài)來預(yù)測當(dāng)前狀態(tài),其計(jì)算如式(1)所示。

    ut-=Gtut-1+Htut

    (1)

    式中:

    Gt——

    當(dāng)前時刻的狀態(tài)轉(zhuǎn)移矩陣(就是指從上一狀態(tài)轉(zhuǎn)變?yōu)橄乱粻顟B(tài)的關(guān)系矩陣);

    Ht——

    當(dāng)前時刻的控制矩陣(就是指影響控制量的控制矩陣);

    ut——

    當(dāng)前時刻的控制量;

    ut-1——

    上一時刻最優(yōu)估計(jì)值;

    ut-——當(dāng)前時刻估計(jì)值。

    預(yù)測均值的協(xié)方差來計(jì)算Kalman增益,其中協(xié)方差的計(jì)算如式(2)所示。

    t=Gtt-1GtT+Rt

    (2)

    式中:

    t-1——上一時刻預(yù)測值方差矩陣;

    t——不確定性的協(xié)方差矩陣;

    Rt——當(dāng)前時刻測量值噪聲矩陣。

    Kalman增益的計(jì)算如式(3)所示。

    Kt=tWtT(Wtt+Qt)-1

    (3)

    式中:

    Kt——卡爾曼增益(指的就是權(quán)值);

    Wt——測量方程;

    Qt——觀察量的協(xié)方差矩陣。

    通過Kalman增益來計(jì)算狀態(tài)估計(jì)值,其估計(jì)位姿計(jì)算如式(4)所示。

    pt=ut-+Kt[zt-h(ut-)]

    (4)

    式中:

    zt——傳感器在t時刻的測量值;

    pt——估計(jì)位姿。

    預(yù)測當(dāng)前狀態(tài)需要用到上一狀態(tài)的協(xié)方差,還需要計(jì)算當(dāng)前狀態(tài)的協(xié)方差用于下一次迭代。其中,當(dāng)前狀態(tài)的協(xié)方差矩陣計(jì)算如式(5)所示。

    t=(1-KtWt)t

    (5)

    由式(4)和式(5)得到了里程計(jì)與IMU數(shù)據(jù)融合后的最佳估計(jì)位姿ut,以及它的協(xié)方差矩陣t。最后根據(jù)激光雷達(dá)數(shù)據(jù)建立測量模型來更新所有隨機(jī)粒子的權(quán)重,粒子最終的收斂位置即為機(jī)器人在地圖中的實(shí)際位置。若機(jī)器人初始位姿未知或有誤,此時全局將注入更多隨機(jī)粒子,隨機(jī)器人的移動粒子再次更新,重新定位機(jī)器人在羊舍地圖中的位置。

    2.2"羊舍地圖構(gòu)建

    采用Gmapping算法對羊舍進(jìn)行地圖的構(gòu)建。Gmapping是一種基于Rao-Blackwellized粒子濾波的二維激光SLAM方法[9]。即將定位和建圖過程分離,先進(jìn)行定位再進(jìn)行建圖。為了解決RBPF頻繁重采樣的缺點(diǎn),Gmapping提出了選擇性重采樣,即減少重采樣的次數(shù),設(shè)定閾值;當(dāng)粒子權(quán)重大于所設(shè)閾值才執(zhí)行重采樣,這樣便減少了采樣的次數(shù),即減緩了粒子的退化[10]。里程計(jì)根據(jù)測量的速度可以推導(dǎo)出機(jī)器人的路程和方向,Gmapping算法依據(jù)里程計(jì)的測量數(shù)據(jù)估計(jì)機(jī)器人的運(yùn)動模型,高精度的里程計(jì)對于建圖和定位性能的提高十分重要。因此,首先對Gmapping功能包中的雷達(dá)最大檢測距離、粒子個數(shù)以及機(jī)器人激光匹配頻率進(jìn)行調(diào)試并取得最優(yōu)值,最后利用Gmapping功能包訂閱RPLIDAR S1發(fā)出的雷達(dá)話題信息和擴(kuò)展卡爾曼濾波融合的位姿信息進(jìn)行小車的精確定位和建圖。

    2.3"巡檢路線規(guī)劃

    2.3.1"巡檢目標(biāo)點(diǎn)獲取

    在巡檢機(jī)器人對羊舍環(huán)境進(jìn)行建模并得到二維柵格地圖后,需要獲取地圖上的目標(biāo)點(diǎn)坐標(biāo)信息從而實(shí)現(xiàn)巡檢的功能,為了獲取地圖上的目標(biāo)點(diǎn)坐標(biāo)信息,用Python QT構(gòu)建一個簡單的界面,界面如圖3所示。

    該界面通過點(diǎn)擊獲取節(jié)點(diǎn)的按鍵監(jiān)聽base_footprint與map坐標(biāo)的tf,從而得到機(jī)器人實(shí)時位置,通過點(diǎn)擊記錄位點(diǎn)的按鍵將獲取的機(jī)器人位置以json的數(shù)據(jù)格式存儲在當(dāng)前功能包的data目錄下??梢栽谏瞎?jié)進(jìn)行手動建圖時獲取并記錄各個目標(biāo)點(diǎn)的位置信息,以便在導(dǎo)航時可以直接讀取各個目標(biāo)點(diǎn)的位置信息進(jìn)行巡檢導(dǎo)航。

    2.3.2"巡檢路線生成

    根據(jù)已經(jīng)獲取的巡檢線路上的目標(biāo)點(diǎn)信息,生成一條巡檢路線,并發(fā)布到路徑規(guī)劃器中。開啟前需要將機(jī)器人放置到巡檢的固定起點(diǎn),打開定位節(jié)點(diǎn),并載入建好的羊舍地圖,機(jī)器人實(shí)際巡檢路線如圖4所示。

    如圖4所示,機(jī)器人將沿著A→B→C→D→E→D→C→B→A路線行駛,車體的正前方為x軸正方向,左方為y軸正方向。2.3.1節(jié)已經(jīng)通過QT界面獲取到各個目標(biāo)點(diǎn)相對于地圖的xy坐標(biāo)和朝向信息并存儲在指定文件夾下。用C++中的ifstream類依次讀取文件中的坐標(biāo),之后用ROS中的LINE_STRIP類型將每個點(diǎn)用作一組連接的線中的頂點(diǎn),其中點(diǎn)A連接到點(diǎn)B,點(diǎn)B連接到點(diǎn)C,點(diǎn)C連接到點(diǎn)D等,這樣就生成一條巡檢線路,最后以trajectory命名建立一個話題,把坐標(biāo)點(diǎn)位姿信息發(fā)布出去。在RVIZ可視化界面訂閱此話題,路線就會顯示在建好的地圖上。

    機(jī)器人的巡檢方式實(shí)際是依次遍歷巡檢線路上的每個目標(biāo)點(diǎn)位置信息從而實(shí)現(xiàn)循環(huán)往返,這是通過訂閱發(fā)布的trajectory話題將規(guī)劃線路傳入到自定義的全局規(guī)劃器。此過程首先通過監(jiān)聽base_footprint與map坐標(biāo)系的tf獲取機(jī)器人當(dāng)前位置,以當(dāng)前位置作為開始的規(guī)劃變量。然后,在for循環(huán)中依次遍歷插入巡檢路線上的各個目標(biāo)點(diǎn)的位置,同時對tf進(jìn)行監(jiān)聽獲取導(dǎo)航線上各個點(diǎn)的位置,獲取到所有目標(biāo)點(diǎn)位置后判斷機(jī)器人當(dāng)前位置與目標(biāo)點(diǎn)的相對距離,其中,目標(biāo)點(diǎn)與當(dāng)前位置相對距離的計(jì)算如式(6)所示。最后,根據(jù)式(6)將距離機(jī)器人當(dāng)前位置最近的目標(biāo)位置相對距離插入規(guī)劃。這里設(shè)置為首先規(guī)劃距離當(dāng)前位置最近的目標(biāo)點(diǎn)是為了防止機(jī)器人在巡檢過程中出現(xiàn)意外情況被迫停止巡檢時需要回到起點(diǎn)從頭開始巡檢的情況,這樣即使在巡檢過程中出現(xiàn)情況也可以在當(dāng)前位置直接規(guī)劃離自己最近的目標(biāo)點(diǎn)從而之后按照目標(biāo)點(diǎn)順序繼續(xù)巡檢。

    Distance=[p2(x)-p1(x)]2+[p2(y)-p1(y)]2

    (6)

    式中:

    Distance——目標(biāo)點(diǎn)與當(dāng)前位置相對距離;

    p2(x)——目標(biāo)點(diǎn)相對于地圖的x坐標(biāo);

    p1(x)——當(dāng)前位置點(diǎn)相對于地圖的x坐標(biāo);

    p2(y)——目標(biāo)點(diǎn)相對于地圖的y坐標(biāo);

    p1(y)——當(dāng)前位置點(diǎn)相對于地圖的y坐標(biāo)。

    除了發(fā)布相對距離之外還需要發(fā)表各個目標(biāo)點(diǎn)朝向來使巡檢過程更加順暢,機(jī)器人的姿態(tài)信息通過訂閱IMU話題信息得到,其姿態(tài)信息由四元數(shù)表示。由圖4的巡檢路線可知,A-B段路線上的目標(biāo)點(diǎn)朝向只需和機(jī)器人起點(diǎn)的姿態(tài)保持一致,而B-C段路線上的目標(biāo)點(diǎn)朝向相比起點(diǎn)姿態(tài)需要繞Z軸逆時針旋轉(zhuǎn)90°;C-D段路線上的目標(biāo)點(diǎn)朝向相比起點(diǎn)姿態(tài)同樣跟開機(jī)方向保持一致;D-E段路線上的目標(biāo)點(diǎn)朝向相比起點(diǎn)姿態(tài)需要繞Z軸順時針旋轉(zhuǎn)90°。給定一個歐拉旋轉(zhuǎn)(X,Y,Z)(即分別繞x軸、y軸和z軸旋轉(zhuǎn)X、Y、Z角度),歐拉角與四元數(shù)的轉(zhuǎn)換計(jì)算如式(7)所示。

    x=

    sinY2sinZ2cosX2+cosY2

    cosZ2sinX2

    y=

    sinY2cosZ2cosX2+cosY2

    sinZ2sinX2

    z=

    cosY2sinZ2cosX2-sinY2

    cosZ2sinX2

    w=

    cosY2sinZ2cosX2-sinY2

    cosZ2sinX2

    (7)

    機(jī)器人開機(jī)的朝向即為機(jī)器人正方向,如圖4所示機(jī)器人在A-B段路線的目標(biāo)點(diǎn)朝向與開機(jī)方向一致,所以A-B段路線上目標(biāo)點(diǎn)的歐拉旋轉(zhuǎn)為(0,0,0)對應(yīng)的四元數(shù)為:x=0,y=0,z=0,w=1。B-C段路線上的目標(biāo)點(diǎn)的歐拉旋轉(zhuǎn)為(0,0,-90),對應(yīng)的四元數(shù)為:x=0,y=0,z=-0.707,w=-0.707。C-D段路線上目標(biāo)點(diǎn)的歐拉旋轉(zhuǎn)為(0,0,0)對應(yīng)的四元數(shù)為:x=0,y=0,z=0,w=1。D-E段路線上目標(biāo)點(diǎn)的歐拉旋轉(zhuǎn)為(0,0,90),對應(yīng)的四元數(shù)為:x=0,y=0,z=0.707,w=0.707。

    最后,將各個目標(biāo)點(diǎn)之間的相對距離和朝向信息依次發(fā)送到move_base全局規(guī)劃模塊,返回原點(diǎn)的各段目標(biāo)點(diǎn)朝向可以根據(jù)上述公式類推得到。該模塊將通過ROS話題global_plan進(jìn)行發(fā)布進(jìn)行路徑規(guī)劃,然后將由本地規(guī)劃模塊接收并發(fā)布控制指令。這樣就實(shí)現(xiàn)了機(jī)器人巡檢線路生成以及路線上目標(biāo)點(diǎn)位姿的發(fā)布。

    2.4"路徑規(guī)劃算法

    2.4.1"A*全局路徑規(guī)劃

    在將目標(biāo)點(diǎn)路徑和朝向信息發(fā)布后,機(jī)器人按照規(guī)劃的路徑進(jìn)行移動。在全局路徑規(guī)劃中通過點(diǎn)對點(diǎn)路徑規(guī)劃可以實(shí)現(xiàn)巡檢機(jī)器人從起始點(diǎn)到目標(biāo)點(diǎn)的最優(yōu)路徑選擇,本文采用A*算法進(jìn)行全局路徑規(guī)劃。

    A*算法因?yàn)榭梢钥焖俚匾?guī)劃出最短路徑,所以多用在全局路徑規(guī)劃算法。A*算法將啟發(fā)式算法融入到評估函數(shù)里[11],其評估函數(shù)計(jì)算如式(8)所示。

    F(n)=G(n)+H(n)

    (8)

    式中:

    F(n)——綜合代價(jià)值,即評估函數(shù);

    G(n)——起點(diǎn)到當(dāng)前點(diǎn)的實(shí)際代價(jià)值;

    H(n)——當(dāng)前點(diǎn)到目標(biāo)點(diǎn)的估計(jì)代價(jià)值。

    A*算法中的估計(jì)代價(jià)值H(n)是優(yōu)于其他算法的核心,在路徑規(guī)劃過程中必須保證當(dāng)前點(diǎn)到目標(biāo)點(diǎn)的估計(jì)代價(jià)值低于實(shí)際距離。選用歐幾里得距離作為啟發(fā)函數(shù)以保證降低A*算法計(jì)算量的同時,進(jìn)一步優(yōu)化全局最優(yōu)的表達(dá)。

    2.4.2"DWA局部路徑規(guī)劃

    由于羊場環(huán)境的復(fù)雜性,若機(jī)器人在巡檢過程中突然出現(xiàn)一個動態(tài)障礙物,這時僅使用A*算法進(jìn)行路徑規(guī)劃是不能實(shí)現(xiàn)動態(tài)避障的。DWA算法是一種局部路徑規(guī)劃算法,適用于移動機(jī)器人等需要避開障礙物的場景。其主要流程包含:初始化、采樣速度樣本、樣本評分以及發(fā)布方案四部分[12]。首先,通過類加載模塊載入ROS局部路徑規(guī)劃包 Base Local Planner中子類的實(shí)例,并調(diào)用其初始化函數(shù)獲取機(jī)器人的初始狀態(tài)信息,進(jìn)一步得到巡檢機(jī)器人的軌跡運(yùn)動模型。基于軌跡運(yùn)動模型,機(jī)器人根據(jù)其運(yùn)動速度即可推算出其運(yùn)動軌跡。為了獲得機(jī)器人運(yùn)動速度樣本,需在一定時間間隔范圍內(nèi)通過傳感器對機(jī)器人各個維度的線速度和角速度進(jìn)行采集,并將其以結(jié)構(gòu)體的形式存儲到相應(yīng)的容器內(nèi)。在得到機(jī)器人速度樣本后,根據(jù)機(jī)器人采樣速度模擬推演出其相對應(yīng)的運(yùn)動軌跡,并通過軌跡評價(jià)函數(shù)對每一條軌跡進(jìn)行評價(jià)。評價(jià)函數(shù)的計(jì)算如式(9)所示。

    G(v,w)=

    max(φhead(v,w)+βdist(v,w))+δvelo(v,w))

    (9)

    式中:

    head(v,w)——速度軌跡與目標(biāo)點(diǎn)的接近程度;

    dist(v,w)——

    該采樣速度下運(yùn)動估計(jì)距離最近障礙物的距離;

    velo(v,w)——

    機(jī)器人在該速度組下的前進(jìn)效率;

    φ、β、δ——

    三個不同子項(xiàng)在評價(jià)函數(shù)中所占比重。

    利用式(9)評價(jià)所有軌跡,選擇評分最高的軌跡為局部規(guī)劃的結(jié)果輸出給move_base控制機(jī)器人移動。機(jī)器人在路徑規(guī)劃過程中通過RPLIDAR S1激光雷達(dá)檢測周圍環(huán)境障礙物信息,若機(jī)器人在導(dǎo)航過程中出現(xiàn)了動態(tài)障礙物,根據(jù)move_base中的柵格代價(jià)地圖(Cost map)模塊進(jìn)行檢測障礙物與機(jī)器人能否發(fā)生碰撞。代價(jià)地圖通常采用網(wǎng)格(grid)的形式,每個網(wǎng)格可以存儲0~255之間不同的數(shù)值,不同的數(shù)值代表該網(wǎng)格離障礙物的遠(yuǎn)近程度。根據(jù)數(shù)值所在區(qū)間判斷當(dāng)前屬于行駛自由區(qū)域還是受限區(qū)域。若為行駛自由區(qū)域,則不需要進(jìn)行其他操作,與之前軌跡保持一致即可;若為行駛受限區(qū)域,則進(jìn)一步獲取機(jī)器人模擬軌跡的角速度,當(dāng)角速度超過一定閾值(若以該閾值進(jìn)行移動,將進(jìn)入受限區(qū)域)后,對該軌跡的權(quán)重加大,并降低軌跡評分,進(jìn)而選擇更優(yōu)軌跡進(jìn)行運(yùn)動控制。最終實(shí)現(xiàn)機(jī)器人的局部路徑規(guī)劃與實(shí)時避障。

    3"現(xiàn)場試驗(yàn)與結(jié)果分析

    現(xiàn)場試驗(yàn)地點(diǎn)位于河北省衡水市某公司的8號和9號種羊舍。兩個羊舍分別長為60m,寬為8m;總計(jì)70個羊圈。在試驗(yàn)過程中需要將參數(shù)設(shè)置統(tǒng)一,線速度和角速度分別調(diào)整為0.30m/s和0.3rad/s。另外采用RVIZ進(jìn)行地圖可視化。RVIZ是ROS中常用的可視化工具,可以在建模地圖中顯示各傳感器的數(shù)據(jù)。保證PC端與機(jī)器人端連上統(tǒng)一的局域網(wǎng),然后通過終端控制使得在PC端輸入指令可以控制機(jī)器人。

    3.1"機(jī)器人建圖試驗(yàn)

    在實(shí)際測試環(huán)境中分別采用Hector SLAM、Cartographer SLAM及Gmapping SLAM算法構(gòu)建環(huán)境地圖,并分析相同環(huán)境下的地圖精度。3種SLAM在現(xiàn)場環(huán)境中構(gòu)建的地圖如圖5所示。3種算法都完成了羊舍地圖的構(gòu)建,但從圖5(a)中橢圓標(biāo)注的地方可看出,Hector SLAM的建圖結(jié)果產(chǎn)生了與實(shí)際環(huán)境不符的偏移。Cartographer SLAM算法雖采用回環(huán)檢測以校正誤差,但構(gòu)建的環(huán)境地圖出現(xiàn)了部分墻壁線條傾斜偏離實(shí)際環(huán)境的問題,如圖5(b)圓圈所示。相比而言,采用經(jīng)過擴(kuò)展卡爾曼濾波融合里程計(jì)和IMU信息后的Gmapping算法構(gòu)建的環(huán)境地圖邊緣更平整且整體方正,如圖5(c)所示,最符合羊舍環(huán)境的實(shí)際情況。

    3.2"機(jī)器人定位精度測試試驗(yàn)

    為驗(yàn)證巡檢機(jī)器人定位精度,試驗(yàn)將在羊舍場景下對機(jī)器人定位精度進(jìn)行測試。評價(jià)機(jī)器人導(dǎo)航精度指標(biāo)包括橫向偏差、縱向偏差、航向偏角[13]。設(shè)置機(jī)器人起點(diǎn)為坐標(biāo)原點(diǎn),隨機(jī)選取巡航路徑上的4個巡航點(diǎn)。將設(shè)置好的巡航點(diǎn)設(shè)置為(x0,y0),以機(jī)器人底盤中心為測量基準(zhǔn)點(diǎn),機(jī)器人行駛至目標(biāo)位置,記錄其基準(zhǔn)點(diǎn)坐標(biāo)位置為(x1,y1),巡檢機(jī)器人定位精度測試試驗(yàn)示意圖如圖6所示,其中橫向偏差dx和縱向偏差dy的計(jì)算如式(10)和式(11)所示。

    dx=x1-x0

    (10)

    dy=y1-y0

    (11)

    航向偏角為機(jī)器人正方向與y軸之間的夾角。現(xiàn)場測量采樣點(diǎn)與實(shí)際到達(dá)點(diǎn)的位置偏差如圖7所示。圖7(b)中箭頭為設(shè)置的實(shí)際目標(biāo)點(diǎn),兩條實(shí)線為到達(dá)指定位置時機(jī)器人的鋁型材橫杠相對于地面的投影,虛線為兩條實(shí)線中心點(diǎn)的連線,其用于確定機(jī)器人到達(dá)目標(biāo)點(diǎn)位置時機(jī)器人的中心,虛線中心的短橫線即為機(jī)器人到達(dá)指定位置時的底盤中心點(diǎn)。

    其中,該機(jī)器人的起點(diǎn)和目標(biāo)點(diǎn)在全局坐標(biāo)系中的位置和角度如表2所示。

    將各目標(biāo)點(diǎn)處10次試驗(yàn)的各指標(biāo)結(jié)果分別取平均值作為巡檢機(jī)器人在該目標(biāo)點(diǎn)處的定位偏差。機(jī)器人到達(dá)各個采樣點(diǎn)的實(shí)際位置與采樣點(diǎn)位置的橫向絕對偏差、縱向絕對偏差和航向角度絕對偏差如圖8所示。

    由圖8可知,巡檢機(jī)器人在0.3m/s的速度下,橫向平均絕對位置偏差小于7cm,縱向平均絕對位置偏差小于8cm,航向偏角平均絕對位置偏差小于6°。該精度滿足巡檢機(jī)器人在羊舍環(huán)境中的定位要求。

    3.3"機(jī)器人導(dǎo)航精度測試試驗(yàn)

    機(jī)器人在羊舍進(jìn)行巡檢作業(yè)時,需要沿著羊舍可通行區(qū)域進(jìn)行直線行駛。路徑規(guī)劃如圖4所示。機(jī)器人將沿著A→B→C→D→E→D→C→B→A路線進(jìn)行行駛。機(jī)器人在巡檢線路上設(shè)置了77個目標(biāo)點(diǎn),這些目標(biāo)點(diǎn)的位姿信息以json文件格式存儲在系統(tǒng)中,在導(dǎo)航時通過訂閱amcl_pose話題來實(shí)時獲取機(jī)器人在map坐標(biāo)系下位姿信息。通過比較設(shè)置目標(biāo)點(diǎn)位姿與導(dǎo)航時實(shí)時獲取到的目標(biāo)點(diǎn)位姿的誤差來評價(jià)其導(dǎo)航精度。將各目標(biāo)點(diǎn)處3次試驗(yàn)的各指標(biāo)結(jié)果分別取平均值作為巡檢機(jī)器人在該目標(biāo)點(diǎn)處的導(dǎo)航偏差。

    機(jī)器人到達(dá)各個目標(biāo)點(diǎn)時實(shí)時采集的實(shí)際位置與設(shè)置目標(biāo)點(diǎn)的橫向絕對誤差、縱向絕對誤差和航向偏角絕對誤差如圖9所示。

    總體而言,巡檢機(jī)器人在對設(shè)置目標(biāo)點(diǎn)進(jìn)行巡檢導(dǎo)航時橫向偏差最大值小于21cm,平均值不超過8.3cm,縱向偏差最大值小于19cm,平均值不超過7.1cm;航向偏差的平均值不超過6.1°。該精度可滿足巡檢機(jī)器人在羊舍環(huán)境下的導(dǎo)航精度需求。

    3.4"機(jī)器人避障試驗(yàn)

    為了測試機(jī)器人的路徑規(guī)劃和導(dǎo)航避障實(shí)際效果,進(jìn)行導(dǎo)航避障試驗(yàn),試驗(yàn)步驟與定位精度試驗(yàn)大體一致。不同之處在于路徑規(guī)劃之后,在其規(guī)劃的路徑上設(shè)置動態(tài)障礙物,當(dāng)RPLIDAR S1激光雷達(dá)檢測到障礙物并通過代價(jià)地圖模塊判斷到機(jī)器人即將與障礙物相撞時,機(jī)器人內(nèi)部程序啟動DWA算法進(jìn)行避障繞行,經(jīng)一系列計(jì)算后,機(jī)器人得到最優(yōu)軌跡作為避障路徑,這條路徑對應(yīng)的角速度和線速度組合則通過ROS系統(tǒng)傳遞給STM32控制板塊,機(jī)器人得到速度信息后進(jìn)行接下來的避讓行進(jìn),如圖10所示。

    由圖10可知,該機(jī)器人在巡檢過程中能夠檢測并識別出動態(tài)人形障礙物,并且可以順利避開動態(tài)障礙物并到達(dá)目標(biāo)點(diǎn),綜上所述,該機(jī)器人符合在羊舍中的避障要求。

    4"結(jié)論

    1) 設(shè)計(jì)一種基于2D激光雷達(dá)SLAM技術(shù)的羊舍巡檢機(jī)器人導(dǎo)航控制系統(tǒng),提出利用擴(kuò)展卡爾曼濾波融合里程計(jì)和IMU數(shù)據(jù)得到的最佳位姿作為預(yù)測所有隨機(jī)粒子的位姿,該方法提高激光雷達(dá)作為測量模型進(jìn)行更新的準(zhǔn)確性,實(shí)現(xiàn)羊舍的建圖與定位功能。

    2) 針對傳統(tǒng)多目標(biāo)點(diǎn)路徑規(guī)劃需要在每次巡檢之前在RVIZ界面上手動設(shè)置目標(biāo)點(diǎn)的問題,利用QT界面在手動控制機(jī)器人建圖時進(jìn)行目標(biāo)點(diǎn)的設(shè)置與保存,在之后的巡檢中直接依次遍歷保存的目標(biāo)點(diǎn)信息即可實(shí)現(xiàn)羊舍的往返巡檢。

    3) 將建圖定位與路徑規(guī)劃算法應(yīng)用到巡檢機(jī)器人上,結(jié)果表明該機(jī)器人可以完整地在羊舍實(shí)現(xiàn)往返巡檢并規(guī)避障礙物。巡檢機(jī)器人導(dǎo)航橫向平均偏差小于8.3cm,導(dǎo)航縱向平均偏差小于7.1cm,航向平均偏差小于 6.1°。通過試驗(yàn)證明,該方法能夠滿足羊舍的自主導(dǎo)航需求。

    參"考"文"獻(xiàn)

    [1]

    何沛桐, 張建華, 張凝, 等. 基于視覺感知的畜禽智慧養(yǎng)殖管理與疫病診斷研究進(jìn)展[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào), 2023, 28(10): 141-165.

    He Peitong, Zhang Jianhua, Zhang Ning, et al. Research progress in intelligent livestock and poultry breeding management and disease diagnosis based on visual perception [J]. Journal of China Agricultural University, 2023, 28(10): 141-165.

    [2]

    姚文, 鄧為, 許毅, 等. 家禽亞健康狀態(tài)監(jiān)測與健康預(yù)警展望[J]. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2023, 46(4): 635-644.

    Yao Wen, Deng Wei, Xu Yi, et al. Poultry sub-health status monitoring and health warning prospect [J]. Journal of Nanjing Agricultural University, 2023, 46(4): 635-644.

    [3]

    Hu D, Zhong H, Li S, et al. Segmenting areas of potential contamination for adaptive robotic disinfection in built environments [J]. Building and Environment, 2020, 184: 107226.

    [4]

    賀曉峰, 廖志偉, 肖華明, 等. 智能巡檢機(jī)器人傳感與控制系統(tǒng)研究與設(shè)計(jì)[J]. 中國煤炭, 2022, 48(S1): 1-5.

    He Xiaofeng, Liao Zhiwei, Xiao Huaming, et al. Research and design of intelligent sensing and control system of inspection robot [J]. China Coal, 2022, 48 (S1): 1-5.

    [5]

    汪陽, 劉文鍇, 張凌, 等. 羊巴氏桿菌病研究進(jìn)展[J]. 畜牧與獸醫(yī), 2023, 55(8): 126-132.

    Wang Yang, Liu Wenkai, Zhang Ling, et al. Research progress on pasteurellosis in sheep and goats [J]. Animal Husbandry Veterinary Medicine, 2023, 55(8): 126-132.

    [6]

    曾志雄, 董冰, 呂恩利, 等. 豬舍環(huán)境無線多點(diǎn)多源遠(yuǎn)程監(jiān)測系統(tǒng)設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2020, 51(2): 332-340,349.

    Zeng Zhixiong, Dong Bing, Lü Enli, et al. Design and experiment of wireless multi-point and multi-source remote monitoring system for pig house environment [J]. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(2): 332-340,349.

    [7]

    安石莉. 防疫消殺巡檢機(jī)器人定位導(dǎo)航系統(tǒng)設(shè)計(jì)[D].銀川: 寧夏大學(xué), 2023.

    An Shili. Design of positioning and navigation system for inspection robot epidemic prevention and elimination [D]. Yinchuan: Ningxia University, 2023.

    [8]

    宋懷波, 段援朝, 李嶸, 等.基于激光SLAM的牛場智能推翻草機(jī)器人自主導(dǎo)航系統(tǒng)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2023, 54(2): 293-301.

    Song Huaibo, Duan Yuanchao, Li Rong, et al. Autonomous navigation system for pasture intelligent overthrowing grass robot based on laser SLAM[J]. Transactions of the Chinese Society for Agricultural Machinery, 2023, 54(2): 293-301.

    [9]

    薛金林, 王培曉, 周俊, 等. 基于改進(jìn)Gmapping算法的果園二維環(huán)境地圖精準(zhǔn)構(gòu)建[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2023, 54(7): 26-34, 55.

    Xue Jinlin, Wang Peixiao, Zhou Jun, et al. Accurate construction of orchard two-dimensional environmental map based on improved Gmapping algorithm [J]. Transactions of the Chinese Society for Agricultural Machinery, 2023, 54(7): 26-34, 55.

    [10]

    Zhao J, Li J, Zhou J. Research ontwo-round self-balancing robot SLAM based on the Gmapping algorithm [J]. Sensors, 2023, 23(5): 2489.

    [11]

    溫淑慧, 問澤藤, 劉鑫, 等. 基于ROS的移動機(jī)器人自主建圖與路徑規(guī)劃[J]. 沈陽工業(yè)大學(xué)學(xué)報(bào), 2022, 44(1): 90-94.

    Wen Shuhui, Wen Zeteng, Liu Xin, et al. Autonomous mapping and path planning of mobile robot based on ROS [J]. Journal of Shenyang University of Technology, 2022, 44(1): 90-94.

    [12]

    王豪杰, 馬向華, 代婉玉, 等. 改進(jìn)DWA算法的移動機(jī)器人避障研究[J]. 計(jì)算機(jī)工程與應(yīng)用, 2023, 59(6): 326-332.

    Wang Haojie, Ma Xianghua, Dai Wanyu, et al. Research on obstacle avoidance of mobile robot with improved DWA algorithm [J]. Computer Engineering and Applications, 2023, 59(6): 326-332.

    [13]

    姜賽珂. 溫室巡檢機(jī)器人導(dǎo)航控制系統(tǒng)研究[D].鎮(zhèn)江:江蘇大學(xué), 2023.

    Jiang Saike. Research on navigation control system of greenhouse inspection robot [D]. Zhenjiang: Jiangsu University, 2023.

    国产成人精品福利久久| 天美传媒精品一区二区| 久久久久网色| 高清在线视频一区二区三区| 欧美区成人在线视频| 欧美日韩在线观看h| 夫妻午夜视频| a级毛色黄片| 亚洲丝袜综合中文字幕| 日韩av在线免费看完整版不卡| 大码成人一级视频| 黄色视频在线播放观看不卡| 男人狂女人下面高潮的视频| 国产成人freesex在线| 亚洲不卡免费看| 黑人巨大精品欧美一区二区蜜桃 | 如日韩欧美国产精品一区二区三区 | 久久国产亚洲av麻豆专区| 久久97久久精品| 能在线免费看毛片的网站| 欧美日韩av久久| 国产一级毛片在线| 国产日韩欧美亚洲二区| 午夜视频国产福利| 99国产精品免费福利视频| 99精国产麻豆久久婷婷| 最新的欧美精品一区二区| 亚洲欧洲国产日韩| 男女国产视频网站| 天堂8中文在线网| 欧美xxxx性猛交bbbb| 精品视频人人做人人爽| 国产极品天堂在线| 偷拍熟女少妇极品色| a级毛色黄片| 97超碰精品成人国产| 午夜日本视频在线| 国产精品一区二区在线不卡| 一级av片app| 在线精品无人区一区二区三| 3wmmmm亚洲av在线观看| 午夜91福利影院| 亚洲av在线观看美女高潮| 女人久久www免费人成看片| 久久精品国产自在天天线| 久久97久久精品| 亚洲国产成人一精品久久久| 欧美区成人在线视频| 美女国产视频在线观看| 春色校园在线视频观看| 国产精品麻豆人妻色哟哟久久| 美女脱内裤让男人舔精品视频| 久久久久久伊人网av| 美女内射精品一级片tv| 乱人伦中国视频| 久久久久久久久久人人人人人人| 国产亚洲91精品色在线| 天美传媒精品一区二区| 久久女婷五月综合色啪小说| 各种免费的搞黄视频| 国产高清三级在线| 青春草国产在线视频| 午夜福利视频精品| 一本久久精品| 亚洲成人一二三区av| 亚洲av成人精品一区久久| av视频免费观看在线观看| 亚洲精品亚洲一区二区| 综合色丁香网| 五月天丁香电影| 亚洲精品一区蜜桃| 日本黄色日本黄色录像| 中文天堂在线官网| 国产精品偷伦视频观看了| av福利片在线| 有码 亚洲区| 亚洲成人手机| 99热这里只有是精品在线观看| 老女人水多毛片| 色94色欧美一区二区| 3wmmmm亚洲av在线观看| 国产男女内射视频| 国产精品女同一区二区软件| 美女内射精品一级片tv| 麻豆成人午夜福利视频| 香蕉精品网在线| 精品一品国产午夜福利视频| 国产精品一区二区在线观看99| 日产精品乱码卡一卡2卡三| 女人精品久久久久毛片| 国产亚洲最大av| 99九九在线精品视频 | 黄色日韩在线| 一区二区三区精品91| 五月玫瑰六月丁香| 国模一区二区三区四区视频| 亚洲国产精品国产精品| 插逼视频在线观看| 国产精品一区二区在线不卡| 国产精品蜜桃在线观看| 成人亚洲精品一区在线观看| 一级a做视频免费观看| 国产中年淑女户外野战色| 国产伦在线观看视频一区| 69精品国产乱码久久久| 亚洲综合精品二区| 熟女av电影| 少妇人妻一区二区三区视频| 美女cb高潮喷水在线观看| 日韩精品免费视频一区二区三区 | 国产精品国产av在线观看| 国产精品福利在线免费观看| 毛片一级片免费看久久久久| 欧美97在线视频| 男女啪啪激烈高潮av片| 国产av国产精品国产| 亚洲精品乱码久久久v下载方式| h日本视频在线播放| 亚洲av男天堂| 欧美精品亚洲一区二区| 日韩,欧美,国产一区二区三区| 日本午夜av视频| 久久国产亚洲av麻豆专区| 久久久久久久久久久丰满| 色网站视频免费| 国产亚洲欧美精品永久| 欧美丝袜亚洲另类| 国产熟女欧美一区二区| 国产亚洲91精品色在线| 日韩av免费高清视频| 日本vs欧美在线观看视频 | 五月开心婷婷网| 性高湖久久久久久久久免费观看| 老司机影院毛片| 国产av码专区亚洲av| 热re99久久精品国产66热6| 婷婷色av中文字幕| 亚洲国产毛片av蜜桃av| 国产欧美日韩精品一区二区| 岛国毛片在线播放| 高清黄色对白视频在线免费看 | 亚洲综合精品二区| 欧美日韩亚洲高清精品| 熟女av电影| 国产精品一区二区性色av| 亚洲国产精品专区欧美| 欧美日韩视频高清一区二区三区二| 亚洲精品日本国产第一区| 97在线人人人人妻| 女性生殖器流出的白浆| 国产黄频视频在线观看| 免费久久久久久久精品成人欧美视频 | 国产国拍精品亚洲av在线观看| 精品少妇内射三级| 免费黄频网站在线观看国产| 天天躁夜夜躁狠狠久久av| 人人澡人人妻人| 99热这里只有精品一区| 中文天堂在线官网| 男男h啪啪无遮挡| 久久久久久久亚洲中文字幕| av播播在线观看一区| 成人影院久久| 国产高清有码在线观看视频| 成人特级av手机在线观看| 免费观看无遮挡的男女| 日本wwww免费看| 亚洲内射少妇av| 好男人视频免费观看在线| 日本黄色片子视频| 少妇被粗大的猛进出69影院 | 国产淫片久久久久久久久| 色婷婷久久久亚洲欧美| 国产淫片久久久久久久久| 乱人伦中国视频| 欧美老熟妇乱子伦牲交| 国产女主播在线喷水免费视频网站| 七月丁香在线播放| 一区二区三区免费毛片| 久久精品国产a三级三级三级| 久久精品国产亚洲av涩爱| 只有这里有精品99| 精品一区在线观看国产| a 毛片基地| 成年女人在线观看亚洲视频| 中文天堂在线官网| 国产成人精品婷婷| 男人和女人高潮做爰伦理| 免费黄频网站在线观看国产| 亚洲电影在线观看av| 七月丁香在线播放| 久久韩国三级中文字幕| 日韩成人伦理影院| 日韩av免费高清视频| 香蕉精品网在线| 男女免费视频国产| av国产精品久久久久影院| 久久国产精品男人的天堂亚洲 | 国产精品久久久久久久电影| 狂野欧美激情性xxxx在线观看| 日韩 亚洲 欧美在线| 亚洲三级黄色毛片| 少妇人妻久久综合中文| 国内精品宾馆在线| 国产在线男女| 亚洲精品乱久久久久久| 在线亚洲精品国产二区图片欧美 | 高清av免费在线| 女的被弄到高潮叫床怎么办| 亚洲精品乱码久久久久久按摩| 亚洲成人手机| 熟女av电影| 中国三级夫妇交换| 如何舔出高潮| 久久99精品国语久久久| 成年人午夜在线观看视频| 看免费成人av毛片| 丰满迷人的少妇在线观看| 日本91视频免费播放| 亚洲美女黄色视频免费看| videossex国产| 观看免费一级毛片| 亚洲精华国产精华液的使用体验| 插逼视频在线观看| 七月丁香在线播放| 色视频在线一区二区三区| av女优亚洲男人天堂| 乱码一卡2卡4卡精品| 熟女电影av网| 久久久久视频综合| 日本欧美国产在线视频| 成人美女网站在线观看视频| 国产成人精品无人区| 欧美另类一区| 久久99热6这里只有精品| 欧美激情国产日韩精品一区| 婷婷色综合大香蕉| 国产成人a∨麻豆精品| 极品少妇高潮喷水抽搐| 国产精品久久久久久久电影| 国产一区亚洲一区在线观看| 麻豆成人午夜福利视频| 午夜视频国产福利| av播播在线观看一区| 天美传媒精品一区二区| 国产成人精品一,二区| 哪个播放器可以免费观看大片| 国产精品久久久久久精品电影小说| 女人久久www免费人成看片| 男女国产视频网站| 日韩精品免费视频一区二区三区 | 久久精品熟女亚洲av麻豆精品| 少妇人妻 视频| 三级国产精品片| 寂寞人妻少妇视频99o| av在线播放精品| 80岁老熟妇乱子伦牲交| 大片电影免费在线观看免费| 亚洲国产毛片av蜜桃av| 久久久久久伊人网av| 一级av片app| 天天躁夜夜躁狠狠久久av| 最后的刺客免费高清国语| 在线精品无人区一区二区三| 亚洲av电影在线观看一区二区三区| 国产色爽女视频免费观看| av福利片在线| 免费看日本二区| 欧美精品人与动牲交sv欧美| 热99国产精品久久久久久7| 国产精品.久久久| 99九九线精品视频在线观看视频| 国产成人午夜福利电影在线观看| √禁漫天堂资源中文www| 日本91视频免费播放| 看非洲黑人一级黄片| 午夜激情福利司机影院| 午夜福利网站1000一区二区三区| 大香蕉97超碰在线| 夜夜看夜夜爽夜夜摸| 久久国产亚洲av麻豆专区| 日日摸夜夜添夜夜添av毛片| 激情五月婷婷亚洲| 成人午夜精彩视频在线观看| 精品久久国产蜜桃| 欧美97在线视频| 99热网站在线观看| 精品亚洲乱码少妇综合久久| 一区二区av电影网| 九九爱精品视频在线观看| 国产一区二区三区av在线| 久久久久国产精品人妻一区二区| 国产永久视频网站| 黄色视频在线播放观看不卡| 男的添女的下面高潮视频| 美女cb高潮喷水在线观看| 一个人看视频在线观看www免费| 亚洲欧洲国产日韩| 欧美人与善性xxx| 纵有疾风起免费观看全集完整版| 国产国拍精品亚洲av在线观看| 欧美丝袜亚洲另类| 午夜免费观看性视频| 国产av码专区亚洲av| 少妇的逼好多水| 日本欧美国产在线视频| 国产色爽女视频免费观看| 亚洲精品aⅴ在线观看| 五月玫瑰六月丁香| 22中文网久久字幕| 一级片'在线观看视频| 一本一本综合久久| 在线观看免费视频网站a站| 亚洲av中文av极速乱| 精品少妇久久久久久888优播| 十八禁高潮呻吟视频 | 国产日韩一区二区三区精品不卡 | 精品人妻熟女av久视频| 日日摸夜夜添夜夜添av毛片| 成人国产麻豆网| 久久久久精品性色| 伊人久久精品亚洲午夜| 最近最新中文字幕免费大全7| 国产片特级美女逼逼视频| 久久久久久久久久成人| 丰满人妻一区二区三区视频av| 欧美xxxx性猛交bbbb| 国产精品.久久久| 久久久国产一区二区| 色婷婷av一区二区三区视频| 亚洲内射少妇av| 好男人视频免费观看在线| 久久久精品免费免费高清| 国产 一区精品| 99久久精品国产国产毛片| 欧美人与善性xxx| 亚州av有码| 一区二区三区免费毛片| 观看美女的网站| 九九爱精品视频在线观看| 男女边摸边吃奶| 一级二级三级毛片免费看| av在线app专区| 欧美日韩一区二区视频在线观看视频在线| av不卡在线播放| av不卡在线播放| 国产免费又黄又爽又色| a级毛片在线看网站| 午夜久久久在线观看| 免费看av在线观看网站| 赤兔流量卡办理| 狠狠精品人妻久久久久久综合| 国产欧美亚洲国产| 草草在线视频免费看| 麻豆乱淫一区二区| 日本免费在线观看一区| 亚洲情色 制服丝袜| 美女视频免费永久观看网站| 久久久久久久久久人人人人人人| 亚洲美女黄色视频免费看| 欧美日韩精品成人综合77777| 2021少妇久久久久久久久久久| 国产伦在线观看视频一区| 中国三级夫妇交换| 久久精品国产亚洲网站| 插逼视频在线观看| 在线观看美女被高潮喷水网站| 少妇精品久久久久久久| 成人黄色视频免费在线看| 亚洲va在线va天堂va国产| 一本一本综合久久| kizo精华| 多毛熟女@视频| 精品一区二区三卡| 97超视频在线观看视频| 亚洲av不卡在线观看| 三级国产精品欧美在线观看| 大片免费播放器 马上看| 国产成人一区二区在线| 免费不卡的大黄色大毛片视频在线观看| 在线观看www视频免费| 高清黄色对白视频在线免费看 | 亚洲精品一二三| 久久久精品94久久精品| 日韩电影二区| 中文精品一卡2卡3卡4更新| 涩涩av久久男人的天堂| 又大又黄又爽视频免费| 婷婷色av中文字幕| 日产精品乱码卡一卡2卡三| 日韩亚洲欧美综合| 岛国毛片在线播放| 日韩伦理黄色片| av福利片在线| 国内少妇人妻偷人精品xxx网站| 中文字幕人妻丝袜制服| 99九九线精品视频在线观看视频| 国产精品久久久久久精品电影小说| 亚洲国产精品成人久久小说| 亚洲欧美日韩东京热| av一本久久久久| 99久久人妻综合| 22中文网久久字幕| 热re99久久国产66热| 国产成人freesex在线| 我要看日韩黄色一级片| 一个人免费看片子| 搡女人真爽免费视频火全软件| 秋霞在线观看毛片| 亚洲国产精品一区二区三区在线| 高清不卡的av网站| 另类亚洲欧美激情| 精品久久久久久久久av| 国产乱人偷精品视频| 精品国产乱码久久久久久小说| 99热6这里只有精品| 国产成人精品无人区| 人妻人人澡人人爽人人| 女人久久www免费人成看片| 我要看日韩黄色一级片| 久久精品熟女亚洲av麻豆精品| 中国三级夫妇交换| 国产欧美日韩一区二区三区在线 | 男的添女的下面高潮视频| 亚洲综合精品二区| 久久久国产精品麻豆| 观看美女的网站| 少妇被粗大的猛进出69影院 | 在线观看人妻少妇| 亚洲av成人精品一二三区| 欧美日韩亚洲高清精品| 中文欧美无线码| 欧美xxxx性猛交bbbb| 国产免费一级a男人的天堂| 亚洲婷婷狠狠爱综合网| www.av在线官网国产| 黑人猛操日本美女一级片| 黑人巨大精品欧美一区二区蜜桃 | 91精品国产九色| 亚洲av综合色区一区| 日日摸夜夜添夜夜添av毛片| 国产伦精品一区二区三区视频9| 欧美日韩av久久| 成人午夜精彩视频在线观看| 久久av网站| 亚洲怡红院男人天堂| 高清欧美精品videossex| 青春草国产在线视频| 高清午夜精品一区二区三区| av在线播放精品| 美女中出高潮动态图| www.色视频.com| 晚上一个人看的免费电影| 最新的欧美精品一区二区| 丝袜脚勾引网站| 色网站视频免费| 日本爱情动作片www.在线观看| 日韩欧美一区视频在线观看 | 亚洲va在线va天堂va国产| 精品少妇黑人巨大在线播放| 黑人猛操日本美女一级片| 高清在线视频一区二区三区| 一区二区av电影网| 国产在线一区二区三区精| 国产又色又爽无遮挡免| 精品一区二区免费观看| 午夜激情久久久久久久| 人妻夜夜爽99麻豆av| 久久青草综合色| 成人美女网站在线观看视频| 十八禁高潮呻吟视频 | 久久人人爽av亚洲精品天堂| 免费黄频网站在线观看国产| 日韩伦理黄色片| 韩国av在线不卡| 亚洲av日韩在线播放| 日韩免费高清中文字幕av| 精品卡一卡二卡四卡免费| 蜜桃久久精品国产亚洲av| 国产精品偷伦视频观看了| 国产亚洲精品久久久com| 欧美激情极品国产一区二区三区 | 国产伦理片在线播放av一区| 夫妻午夜视频| 精品视频人人做人人爽| 777米奇影视久久| 又黄又爽又刺激的免费视频.| 七月丁香在线播放| 色视频www国产| 一级爰片在线观看| 国产探花极品一区二区| 少妇猛男粗大的猛烈进出视频| 黄色日韩在线| 春色校园在线视频观看| 亚洲电影在线观看av| 简卡轻食公司| 久久久午夜欧美精品| 欧美3d第一页| 亚洲精品日韩av片在线观看| 女的被弄到高潮叫床怎么办| 精品久久久久久电影网| 精品人妻一区二区三区麻豆| 七月丁香在线播放| 一个人免费看片子| 国产精品免费大片| 亚洲精品一二三| 观看美女的网站| 看十八女毛片水多多多| 能在线免费看毛片的网站| 永久网站在线| 97精品久久久久久久久久精品| 街头女战士在线观看网站| 成年av动漫网址| 国产成人免费观看mmmm| 亚洲av成人精品一区久久| 色婷婷av一区二区三区视频| 欧美精品高潮呻吟av久久| 大片电影免费在线观看免费| av视频免费观看在线观看| 欧美性感艳星| 久久久久国产网址| 好男人视频免费观看在线| 欧美日韩一区二区视频在线观看视频在线| 在现免费观看毛片| 伊人亚洲综合成人网| 久久99热6这里只有精品| 亚洲熟女精品中文字幕| 亚洲精品乱码久久久v下载方式| 99久久中文字幕三级久久日本| 男女无遮挡免费网站观看| 亚洲欧洲精品一区二区精品久久久 | 免费观看无遮挡的男女| kizo精华| 91久久精品国产一区二区三区| 久久这里有精品视频免费| 精品一区在线观看国产| 中文资源天堂在线| 免费大片18禁| 久久久精品94久久精品| 国产精品麻豆人妻色哟哟久久| 日日撸夜夜添| 久久久久久久亚洲中文字幕| 国产成人精品婷婷| 午夜视频国产福利| 在线观看国产h片| 午夜福利,免费看| 人人妻人人添人人爽欧美一区卜| 国产成人午夜福利电影在线观看| 精品一区在线观看国产| 欧美变态另类bdsm刘玥| 噜噜噜噜噜久久久久久91| 麻豆成人av视频| 久久 成人 亚洲| 色网站视频免费| 在现免费观看毛片| 91午夜精品亚洲一区二区三区| 国产视频内射| 精品久久久久久久久亚洲| 久久久久久人妻| 日日撸夜夜添| 男人添女人高潮全过程视频| 激情五月婷婷亚洲| 高清视频免费观看一区二区| 高清不卡的av网站| 欧美区成人在线视频| 下体分泌物呈黄色| 只有这里有精品99| 狠狠精品人妻久久久久久综合| 色婷婷av一区二区三区视频| 亚洲第一av免费看| av一本久久久久| 国产成人freesex在线| 91久久精品国产一区二区成人| 精品国产一区二区三区久久久樱花| 色视频www国产| 亚洲av福利一区| 亚洲中文av在线| 亚洲精品成人av观看孕妇| 亚洲一级一片aⅴ在线观看| 亚洲在久久综合| 欧美日韩av久久| 亚洲怡红院男人天堂| 尾随美女入室| 久久久久网色| 99久久精品热视频| 麻豆成人av视频| 亚洲国产精品成人久久小说| 精品久久久久久电影网| 高清视频免费观看一区二区| 亚洲国产色片| 亚洲综合色惰| 久久精品国产亚洲av涩爱| 黄色视频在线播放观看不卡| 亚洲综合色惰| 国产在线视频一区二区| 51国产日韩欧美| 国产精品国产三级国产av玫瑰| 欧美日本中文国产一区发布| 97精品久久久久久久久久精品| 亚洲一级一片aⅴ在线观看| 成年人免费黄色播放视频 | 国产精品不卡视频一区二区| 亚洲四区av| 亚洲美女视频黄频| 国产成人一区二区在线| 国产精品久久久久久精品电影小说| 菩萨蛮人人尽说江南好唐韦庄| 欧美激情国产日韩精品一区| 水蜜桃什么品种好| 久久国内精品自在自线图片| 最新的欧美精品一区二区| 狠狠精品人妻久久久久久综合| 22中文网久久字幕| 成人特级av手机在线观看| 国产又色又爽无遮挡免|