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

    基于Gazebo 的深海網(wǎng)箱洗網(wǎng)機(jī)器人動(dòng)力學(xué)仿真平臺(tái)搭建

    2024-03-06 01:45:32黃小華龐國(guó)良袁太平
    南方水產(chǎn)科學(xué) 2024年1期
    關(guān)鍵詞:網(wǎng)衣網(wǎng)箱浮力

    邢 旺,黃小華 ,李 根 ,龐國(guó)良 ,袁太平

    1. 浙江海洋大學(xué),浙江 舟山 316022

    2. 中國(guó)水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所/農(nóng)業(yè)農(nóng)村部外海漁業(yè)可持續(xù)利用重點(diǎn)實(shí)驗(yàn)室/廣東省網(wǎng)箱工程技術(shù)研究中心,廣東 廣州 510300

    3. 三亞熱帶水產(chǎn)研究院,海南 三亞 572024

    我國(guó)是海洋大國(guó),發(fā)展海水養(yǎng)殖能充分發(fā)揮海域廣闊的優(yōu)勢(shì),滿足人們對(duì)于優(yōu)質(zhì)水產(chǎn)動(dòng)物蛋白的需求[1-2]。由于近海養(yǎng)殖日趨飽和,養(yǎng)殖壓力不斷攀升,資源和空間的問(wèn)題日益加劇,因此人們將目光投向了具有水域面積遼闊、水質(zhì)條件優(yōu)良、水體交換速度快等優(yōu)點(diǎn)[3]的深遠(yuǎn)海區(qū)域。

    網(wǎng)箱是深海養(yǎng)殖中不可或缺的組成,大型深海網(wǎng)箱具有一般網(wǎng)箱沒(méi)有的大容量、抗風(fēng)浪和高效益等優(yōu)點(diǎn),同時(shí)大型深海網(wǎng)箱往往配備成套的智能系統(tǒng),易實(shí)現(xiàn)無(wú)人值守和遠(yuǎn)程管控[4]。網(wǎng)衣清洗是深海網(wǎng)箱養(yǎng)殖必不可少的環(huán)節(jié)。網(wǎng)衣得不到及時(shí)清潔會(huì)影響魚類生長(zhǎng),同時(shí)污損生物會(huì)增加網(wǎng)衣的重量,導(dǎo)致網(wǎng)衣、網(wǎng)箱的結(jié)構(gòu)破壞風(fēng)險(xiǎn)增加[5-6]。

    目前,網(wǎng)衣清洗可分為:人工清洗法、機(jī)械沖洗法和洗網(wǎng)機(jī)器人清洗法。人工清洗和機(jī)械清洗均費(fèi)時(shí)費(fèi)力,且會(huì)造成危險(xiǎn)[7]。洗網(wǎng)機(jī)器人具有操作簡(jiǎn)單、安全便攜和自動(dòng)化程度高等優(yōu)點(diǎn)。其主要由控制單元、運(yùn)動(dòng)單元和清洗單元等構(gòu)成[8]??刂萍斑\(yùn)動(dòng)單元解決機(jī)器人貼網(wǎng)移動(dòng)、水下穿梭等問(wèn)題,而清洗單元?jiǎng)t可以通過(guò)多種不同的清洗工具達(dá)到對(duì)網(wǎng)箱網(wǎng)衣的清理。

    總體上,在結(jié)構(gòu)設(shè)計(jì)方面,研究者探索了各種不同的洗網(wǎng)方式,例如高壓旋轉(zhuǎn)水流、毛刷、空化射流等,旨在實(shí)現(xiàn)高效、便捷的網(wǎng)箱清洗。可分為清洗工具、清洗水流亦或兩者綜合使用。在清洗工具方面:通過(guò)清洗盤清洗噴嘴[9]或提供多種類型噴嘴的清洗單元來(lái)進(jìn)行洗網(wǎng)作業(yè)[10],以及采用多個(gè)旋轉(zhuǎn)刀片構(gòu)成洗網(wǎng)清潔頭實(shí)現(xiàn)高效清洗[11]等。在清洗水流方面:通過(guò)高壓旋轉(zhuǎn)水流[12]和通過(guò)空化射流進(jìn)行清洗作業(yè)[13]。綜合兩者還有:采用水流和毛刷作為清洗方式[14]和通過(guò)清洗模塊的形式進(jìn)行清洗作業(yè)[15]等組合模式。

    在運(yùn)動(dòng)控制方面,采用自動(dòng)導(dǎo)航[16-17]、自適應(yīng)抓鉤控制[18]等控制策略,以實(shí)現(xiàn)洗網(wǎng)機(jī)器人的自主貼網(wǎng)移動(dòng)。為應(yīng)對(duì)深海環(huán)境中的復(fù)雜動(dòng)力學(xué),建立動(dòng)力學(xué)模型、開(kāi)展仿真研究也成為熱點(diǎn)?;谝曈X(jué)、聲吶等各種傳感器的控制方案也受到關(guān)注。

    好的洗網(wǎng)機(jī)器人應(yīng)具備水下抗擾動(dòng)、低能耗的特點(diǎn)。然而洗網(wǎng)機(jī)器人在貼網(wǎng)移動(dòng)、水下穿梭的過(guò)程中,受力復(fù)雜、非線性程度高、狀態(tài)難以觀測(cè),要求大量、全面的實(shí)驗(yàn)數(shù)據(jù)來(lái)支撐控制器的優(yōu)化,該方式成本高、耗時(shí)長(zhǎng)。針對(duì)此問(wèn)題,本文在物理仿真環(huán)境Gazebo 的基礎(chǔ)上,設(shè)計(jì)并開(kāi)發(fā)了一套洗網(wǎng)機(jī)器人動(dòng)力學(xué)及運(yùn)動(dòng)學(xué)仿真平臺(tái),通過(guò)仿真來(lái)獲得洗網(wǎng)機(jī)器人水下工作的完備數(shù)據(jù),開(kāi)展不同控制器的全面性能驗(yàn)證,助力洗網(wǎng)機(jī)器人的控制器開(kāi)發(fā)及優(yōu)化。

    1 仿真平臺(tái)架構(gòu)

    針對(duì)機(jī)械臂、無(wú)人駕駛等機(jī)器人仿真問(wèn)題,研究者開(kāi)發(fā)了多款機(jī)器人仿真軟件,包括Webots、V-REP (現(xiàn)為CoppeliaSim)、MuJoCo 和Gazebo 等(圖1),目前暫無(wú)針對(duì)水下機(jī)器人的仿真平臺(tái)[19]。

    本文選擇Gazebo 作為基礎(chǔ)仿真平臺(tái),主要原因是:1) 機(jī)器人操作系統(tǒng)ROS[20]提供了豐富的Gazebo 接口,因此Gazebo 方便聯(lián)合ROS 開(kāi)展機(jī)器人控制器開(kāi)發(fā);2) Gazebo 支持多種物理引擎,如ODE 和Bullet,且有直觀、易操作的三維仿真界面 (圖1-d);3) Gazebo 社區(qū)開(kāi)源了UUV Simulator 項(xiàng)目,該項(xiàng)目可以支持搭建和測(cè)試水下機(jī)器人,仿真水動(dòng)力對(duì)機(jī)器人的作用為本文的開(kāi)發(fā)提供了參考[21]。

    在Gazebo 和ROS 的基礎(chǔ)上,本文設(shè)計(jì)了圖2 所示的仿真系統(tǒng)架構(gòu)。其中綠色部分為基礎(chǔ)平臺(tái)已經(jīng)具備的功能,橙色部分為需要開(kāi)發(fā)的功能。ROS 提供了良好的進(jìn)程間通訊機(jī)制,針對(duì)機(jī)器人控制策略設(shè)計(jì)問(wèn)題,需要開(kāi)發(fā)狀態(tài)估計(jì)模塊、控制模塊、日志記錄模塊等;Gazebo 提供了機(jī)器人建模語(yǔ)言及編譯系統(tǒng)、3D 渲染功能及物理引擎,在此基礎(chǔ)上,開(kāi)發(fā)洗網(wǎng)機(jī)器人本體,并開(kāi)發(fā)水動(dòng)力、推進(jìn)器、傳感器等插件。本文針對(duì)的是動(dòng)力學(xué)仿真問(wèn)題,因此,主要介紹水動(dòng)力插件和推進(jìn)器插件相關(guān)內(nèi)容。通過(guò)ROS-Gazebo 通訊接口,可實(shí)現(xiàn)兩平臺(tái)的信息交換[22]。

    圖2 仿真系統(tǒng)架構(gòu)Fig. 2 Framework of simulation system

    2 洗網(wǎng)機(jī)器人結(jié)構(gòu)建模

    結(jié)構(gòu)對(duì)機(jī)器人動(dòng)力學(xué)和控制性能有非常重要的影響。結(jié)構(gòu)尺寸、外表面影響碰撞、摩擦、附加質(zhì)量[23]等參數(shù),推進(jìn)器布置影響機(jī)器人的姿態(tài)控制,重心浮心位置影響機(jī)器人的穩(wěn)定性,質(zhì)量和慣量決定機(jī)器人的響應(yīng)速度,上述因素綜合影響著機(jī)器人的運(yùn)動(dòng)性能和操作能力。

    構(gòu)建洗網(wǎng)機(jī)器人本體結(jié)構(gòu)步驟如下:1) 先在SolidWorks 軟件中設(shè)計(jì)出機(jī)器人各零部件模型,包括機(jī)器人本體、推進(jìn)器、洗盤3 個(gè)部分,并對(duì)零部件進(jìn)行裝配得到洗網(wǎng)機(jī)器人三維模型 (圖3);2) 利用SolidWorks to URDF Exporter 插件設(shè)置機(jī)體、推進(jìn)器和洗盤的坐標(biāo)系、旋轉(zhuǎn)軸以及父子關(guān)系;3) 獲取機(jī)器人中各部件的相對(duì)位置關(guān)系和關(guān)節(jié)之間的關(guān)系,利用這些信息編寫機(jī)器人模型xacro 文件;4)將各零部件導(dǎo)出為STL 格式用于描述機(jī)器人模型在Gazebo 中的外觀元素部分;5) 設(shè)置機(jī)器人水動(dòng)力模型參數(shù),并將其數(shù)據(jù)編寫軟件包。

    圖3 洗網(wǎng)機(jī)器人本體結(jié)構(gòu)圖Fig. 3 Structure diagram of washing robot body

    3 洗網(wǎng)機(jī)器人動(dòng)力學(xué)建模

    一旦洗網(wǎng)機(jī)器人的受力已知,就可模擬機(jī)器人在水下無(wú)輸入情況下的運(yùn)動(dòng)軌跡,還可模擬機(jī)器人在控制作用下的運(yùn)動(dòng)軌跡,有助于深入理解洗網(wǎng)機(jī)器人的運(yùn)動(dòng)行為,設(shè)計(jì)及驗(yàn)證控制器性能[24]。

    洗網(wǎng)機(jī)器人在工作過(guò)程中受力復(fù)雜,本文僅考慮影響較大的力,包括:浮力、重力、水阻尼力、附加質(zhì)量效應(yīng)、水流作用、推進(jìn)器及洗盤作用。忽略了臍帶纜、湍流、渦流等的影響[25],主要原因?yàn)椋?) 考慮到機(jī)器人移動(dòng)速度慢,設(shè)計(jì)時(shí)選用相對(duì)大功率的推進(jìn)器,因此湍流、渦流對(duì)機(jī)器人的控制影響較?。?) 臍帶纜在設(shè)計(jì)時(shí)會(huì)選用零浮力纜,在大部分場(chǎng)景下對(duì)機(jī)器人運(yùn)動(dòng)無(wú)影響 (少數(shù)情況的影響我們將在后續(xù)深入研究)。

    本文依據(jù)Fossen 手冊(cè)[26]海洋運(yùn)載器的六自由度運(yùn)動(dòng)學(xué)方程建立洗網(wǎng)機(jī)器人動(dòng)力學(xué)模型:

    式中:v為洗網(wǎng)機(jī)器人速度向量,即機(jī)器人本體坐標(biāo)系在世界坐標(biāo)系下的速度在本體坐標(biāo)系下的表示;v˙為洗網(wǎng)機(jī)器人加速度向量,即機(jī)器人本體坐標(biāo)系在世界坐標(biāo)系下的加速度在本體坐標(biāo)系下的表示;η6×1為洗網(wǎng)機(jī)器人在世界坐標(biāo)系下的姿態(tài)向量;M6×6=MRB+MA,為系統(tǒng)的慣性矩陣,MRB和MA分別為剛體和附加質(zhì)量引起的慣性矩陣;C6×6=CRB(vr)+CA(vr),為系統(tǒng)的科里奧利向心矩陣 (下文簡(jiǎn)寫為科氏力),CRB(vr)和CA(vr) 分別為剛體和附加質(zhì)量導(dǎo)致的科里奧利向心矩陣;D6×6=DL+DN,為系統(tǒng)的阻尼矩陣,DL和DN分別為線性和非線性阻尼矩陣;g(η)6×1為系統(tǒng)的重力、浮力的合力矩矢量;g0為用于預(yù)調(diào)平 (壓載控制) 的向量;τ∈R6×1,為推進(jìn)器或其他動(dòng)力元件產(chǎn)生的力和力矩;τb∈R6×1,為風(fēng)、浪和流環(huán)境等不確定因素產(chǎn)生的干擾力。

    下文將分別對(duì)公式 (1) 中出現(xiàn)的附加質(zhì)量、科氏力、阻尼等作用力對(duì)機(jī)器人運(yùn)動(dòng)的影響進(jìn)行仿真。需要注意的是,由于洗網(wǎng)機(jī)器人移動(dòng)范圍小且速度慢,因此地球旋轉(zhuǎn)導(dǎo)致的科氏力可忽略,機(jī)器人旋轉(zhuǎn)運(yùn)動(dòng)導(dǎo)致的科氏力由物理引擎自動(dòng)處理,因此,下文不專門探討科氏力對(duì)機(jī)器人的影響[27]。

    3.1 推進(jìn)器和洗盤

    推進(jìn)器通過(guò)螺旋槳旋轉(zhuǎn)推水給機(jī)器人施加推力和轉(zhuǎn)矩 (圖4),推力和轉(zhuǎn)矩與螺旋槳的轉(zhuǎn)速、水流密度推力系數(shù)等相關(guān),計(jì)算公式[28]為:

    圖4 螺旋槳洗盤對(duì)機(jī)器人施加的推力和轉(zhuǎn)矩Fig. 4 Thrust and torque applied by propeller to robot

    式中:τ為螺旋槳的推力;Q為螺旋槳的轉(zhuǎn)矩;n為螺旋槳的轉(zhuǎn)速;D為螺旋槳的直徑;ρ為流體密度;KTp為螺旋槳的推力系數(shù);KQ為轉(zhuǎn)矩系數(shù)。

    洗盤和推進(jìn)器對(duì)機(jī)器人的作用類似,即提供一個(gè)推力和轉(zhuǎn)矩,因此,采用和推進(jìn)器一樣的公式(2) 和(3),調(diào)整系數(shù)來(lái)獲得洗盤對(duì)機(jī)器人的作用。只是動(dòng)態(tài)特性時(shí)間常數(shù)需要減小,旋轉(zhuǎn)系數(shù)需要增大。

    3.2 重力和浮力

    靜止在水中的機(jī)器人,會(huì)受到重力和浮力作用。如果重心和浮心不重合,重力和浮力將產(chǎn)生恢復(fù)力矩,使得機(jī)器人姿態(tài)在平衡狀態(tài)附近震蕩[29]。水下機(jī)器人受到的浮力和重力如圖5 所示,在NED 坐標(biāo)系中,作用于浮心的浮力,作用于重心上的重力,重力與浮力組合作用可表示為:

    圖5 重力、浮力共同作用示意圖Fig. 5 Schematic diagram of combined action of gravity and buoyancy

    式中:xg、yg、zg、xb、yb、zb分別表示重心和浮心在其坐標(biāo)軸上相應(yīng)的分量;B和W分別代表在z軸方向上的重力和浮力。

    由于恢復(fù)效應(yīng)為重心和浮心不重合所致,故z軸方向的旋轉(zhuǎn)不受恢復(fù)效應(yīng)影響,因此只考慮x和y軸方向波動(dòng)所導(dǎo)致的轉(zhuǎn)矩對(duì)機(jī)器人的影響。

    在僅考慮重力和浮力的影響下,設(shè)置機(jī)器人重心 (0, 0, 0),浮心 (0, 0,zb),分別對(duì)機(jī)器人的x、y和z軸方向施加100 Nm 的沖擊轉(zhuǎn)矩,仿真機(jī)器人的運(yùn)動(dòng)軌跡 (圖6)。

    圖6 Gazebo 中對(duì)機(jī)器人施加轉(zhuǎn)矩 (以x 軸為例)Fig. 6 Applying torque to robots in Gazebo(Taking x-axis as an example)

    由圖7 可見(jiàn),由于重力和浮力的恢復(fù)力距作用,機(jī)器人的x、y軸角速度作周期性的變化,而z軸角速度則保持在一個(gè)定值,這是因?yàn)楦×椭亓Ψ较蚓鶠樨Q直方向,無(wú)法在z軸方向上產(chǎn)生恢復(fù)力矩。

    圖7 重力、浮力作用下的角速度沖擊響應(yīng)Fig. 7 Angular velocity shock response under action of gravity and buoyancy

    3.3 流體阻尼作用

    機(jī)器人在水下運(yùn)動(dòng)會(huì)受到水流的碰撞和摩擦,導(dǎo)致阻尼作用,即速度越大,受到速度反方向的力就越大。阻尼作用可以建模為線性阻尼矩陣DL和非線性阻尼矩陣DN。

    DL用于描述水動(dòng)力引起的線性阻尼效應(yīng)。該阻尼力矩陣是一個(gè)方陣,其元素表示不同自由度上的阻尼系數(shù)。它被用于分析水下機(jī)器人在運(yùn)動(dòng)過(guò)程中受到速度相關(guān)阻尼力的影響。線性阻尼力矩陣的規(guī)模和結(jié)構(gòu)受水下機(jī)器人的幾何形狀、運(yùn)動(dòng)速度以及周圍水流特性的影響 (圖8)。其形式為:

    圖8 阻尼效應(yīng)示意圖Fig. 8 Schematic diagram of damping effect

    式中:Xu、Yv、Zw、Kp、Mq和Nr為一階阻尼系數(shù)。

    DN用于描述水動(dòng)力引起的非線性阻尼效應(yīng)。該矩陣包含了一些非線性阻尼項(xiàng),例如速度的高階項(xiàng)和非線性耦合項(xiàng)。這些項(xiàng)反映了水下機(jī)器人在運(yùn)動(dòng)過(guò)程中非線性阻尼力的貢獻(xiàn),通常是由于機(jī)器人的非線性形變或流體動(dòng)力學(xué)效應(yīng)所引起的 (圖8)。其形式為:

    式中:Xu|u|、Yv|v|、Zw|w|、Kp|p|、Mq|q|、Nr|r|為二階水動(dòng)力系數(shù)。

    兩者結(jié)合即為阻尼矩陣D,表示為:

    參考Berg[29]線性阻尼矩陣DL以及非線性阻尼矩陣DN可設(shè)置為:

    組合兩項(xiàng)作為其阻尼矩陣D為:

    阻尼對(duì)機(jī)器人運(yùn)動(dòng)狀態(tài)影響的仿真情況如下:

    1) 阻尼對(duì)x、y、z軸方向平動(dòng)線速度的影響。對(duì)x、y、z軸方向分別在0.5、0.2 及0.7 s 同時(shí)施加100 N 的單位沖擊力,并繪制在同一幅圖表上(圖9-a)??梢杂^察到由于平動(dòng)不涉及恢復(fù)效應(yīng),3 個(gè)方向均在單位沖擊時(shí)達(dá)到最大線速度,并由于阻尼效應(yīng)隨時(shí)間減?。黄渲衳、y軸方向的速度最終趨近于0,而z軸的速度并不為0,這是由于仿真過(guò)程重力略大于浮力,導(dǎo)致其在初始狀態(tài)有一個(gè)初速度,并且最終在阻尼效應(yīng)下速度回歸于該初速度。

    圖9 阻尼效應(yīng)對(duì)機(jī)器人平動(dòng)和轉(zhuǎn)動(dòng)的影響Fig. 9 Effect of damping on translation and rotation of robots

    2) 阻尼對(duì)x、y、z軸角速度的影響。對(duì)x和z軸在第0.8 秒時(shí)以及y軸在第1.25 秒時(shí)施加100 Nm 的單位沖擊扭矩,并同時(shí)繪制在圖9-b 上。可以觀察到x和y軸的角速度變化規(guī)律完全相同。由于恢復(fù)效應(yīng)的影響,它們均作周期性的阻尼震蕩運(yùn)動(dòng)。且隨著時(shí)間的增加,周期性幅度逐漸減小。而z軸方向由于不受恢復(fù)效應(yīng)的影響,不會(huì)產(chǎn)生周期性的震蕩,而是沿著光滑的凹曲線緩慢下降。

    3.4 慣量及附加質(zhì)量

    附加質(zhì)量的產(chǎn)生是由于機(jī)器人的移動(dòng)帶動(dòng)周圍流體運(yùn)動(dòng)狀態(tài)變化而對(duì)機(jī)器人產(chǎn)生的反作用,用符號(hào)MA表示,其為6×6 的矩陣,包含36 個(gè)參數(shù)[30],由于機(jī)器人本體的對(duì)稱性,本文將附加質(zhì)量矩陣簡(jiǎn)化為對(duì)角矩陣,即附加質(zhì)量不會(huì)改變機(jī)器人速度方向,只改變大小,設(shè)置:

    根據(jù)歐拉力學(xué)[30],剛體質(zhì)量和慣量對(duì)運(yùn)動(dòng)的影響可建模為MRB,即:

    因此,慣量和質(zhì)量共同作用下的影響矩陣為:M6×6=MRB+MA。討論在慣量及附加質(zhì)量共同影響下的仿真情況如下:

    1) 無(wú)附加質(zhì)量下x、y、z軸線性加速度的變化。為探究無(wú)附加質(zhì)量對(duì)線性方向加速度的影響情況,分別對(duì)x、y、z軸施加100 N 的單位沖擊力(圖10-a)。除因仿真中誤差導(dǎo)致的重力略大于浮力使得z軸的線性加速度略小于x和y軸外,其余加速度影響完全相同。

    圖10 單位沖擊對(duì)有無(wú)附加質(zhì)量線性加速度的影響Fig. 10 Influence of unit impact on linear acceleration with or without added mass

    2) 有無(wú)附加質(zhì)量對(duì)線性加速度的影響。經(jīng)過(guò)分析仿真結(jié)果,在有附加質(zhì)量時(shí)x、y、z軸的線性加速度變化完全相同,故只考慮有無(wú)附加質(zhì)量對(duì)線性加速度的影響。以x軸作為觀察附加質(zhì)量對(duì)線性加速度影響的參考。對(duì)機(jī)器人施加公式 (11) 的附加質(zhì)量矩陣,同時(shí)對(duì)其施加100 N 的沖擊力。根據(jù)理論,在施加沖擊力F不變的情況下,增加附加質(zhì)量會(huì)導(dǎo)致其真實(shí)質(zhì)量m增大,機(jī)器人的加速度a應(yīng)該減小。但由于模擬仿真環(huán)境不像真實(shí)世界中時(shí)間為連續(xù)的變量,而是離散變量,故圖10-b 中的反向加速度在真實(shí)時(shí)間中應(yīng)和正向加速度同時(shí)發(fā)生,疊加之后的加速度a就比無(wú)附加質(zhì)量時(shí)小,故其符合理論情況。

    3.5 海流場(chǎng)力作用

    海流是由海洋水體的流動(dòng)引起的,具有一定的速度和方向。因此,機(jī)器人在海流中運(yùn)動(dòng)時(shí),會(huì)受海流作用,產(chǎn)生阻力和推力[31]。考慮到水流對(duì)機(jī)器人的影響非常復(fù)雜,本文對(duì)該作用進(jìn)行簡(jiǎn)化,即水流場(chǎng)的作用等效為在機(jī)器人質(zhì)心處施加一個(gè)速度矢量,具體如圖11 所示。

    圖11 流場(chǎng)圖Fig. 11 Flow field diagram

    其中與機(jī)器人前進(jìn)方向x軸夾角也為水平夾角設(shè)為α,與xOy平面的夾角也為豎直夾角設(shè)為β。且其中的速度、水平夾角和豎直夾角如上圖的參數(shù)。所以其各速度分量可以表示為:

    式中:currentVel 為當(dāng)前水流速度矢量,隨著回調(diào)函數(shù)實(shí)時(shí)更新。

    為增加仿真的真實(shí)性,引入噪聲來(lái)模擬水下環(huán)境中的擾動(dòng)。采用高斯馬爾可夫過(guò)程 (Gauss-Markov process) 來(lái)模擬噪聲的影響。通過(guò)引入噪聲,能夠更準(zhǔn)確地模擬真實(shí)環(huán)境中的隨機(jī)變動(dòng),從而提高仿真的真實(shí)性。高斯馬爾可夫過(guò)程是一種常用于噪聲建模和處理的數(shù)學(xué)模型[32]。它通常用來(lái)描述隨機(jī)信號(hào)或噪聲的變化過(guò)程,具有隨機(jī)性、平穩(wěn)性等特點(diǎn)。本模擬仿真使用下述公式來(lái)引入高斯噪聲:

    式中:var 表示當(dāng)前時(shí)刻的過(guò)程狀態(tài)或變量的值;step 表示時(shí)間步長(zhǎng);mu 表示過(guò)程的阻尼系數(shù);noiseAmp 表示噪聲幅值或強(qiáng)度;random 表示服從均勻分布的隨機(jī)數(shù),用來(lái)引入隨機(jī)性或噪聲成分,模擬真實(shí)世界中不確定的隨機(jī)變動(dòng)。

    在不引入噪聲的情況下,將速度矢量設(shè)為1 m·s?1,且與水平面和豎直面的夾角均為。觀察線速度的變化曲線,海流場(chǎng)作用的仿真情況見(jiàn)圖12。在3 s之前的圖像滿足模擬仿真結(jié)果。但3 s 之后出現(xiàn)了較大波動(dòng),這是因?yàn)闄C(jī)器人設(shè)置的深度為5 m,在波動(dòng)時(shí)刻機(jī)器人沖出水面,然后在阻尼等的作用下出現(xiàn)波動(dòng)并最終趨于穩(wěn)定。

    圖12 無(wú)高斯噪聲流場(chǎng)下的線速度零輸入零狀態(tài)響應(yīng)Fig. 12 Linear velocity zero input zero state response in no gaussian noise flow field

    將mu 設(shè)置為0.3,noiseAmp 設(shè)置為0.5,繼續(xù)觀察對(duì)機(jī)器人速度的影響情況。由圖13 可見(jiàn),噪聲對(duì)線速度和角速度的影響非常大。如何通過(guò)控制策略來(lái)消除噪聲對(duì)機(jī)器人平穩(wěn)運(yùn)行的影響是下一步要研究的問(wèn)題,本文只涉及機(jī)器人模型的建立,未給出解決方案。

    圖13 流場(chǎng)高斯噪聲下的線速度和角速度零輸入零狀態(tài)響應(yīng)Fig. 13 Zero input and zero state response of linear velocity and angular velocity with Gaussian noise of flow field

    3.6 網(wǎng)衣作用

    洗網(wǎng)機(jī)器人在貼網(wǎng)移動(dòng)時(shí),會(huì)受到網(wǎng)衣的作用。網(wǎng)衣對(duì)機(jī)器人施加摩擦力、支持力,影響機(jī)器人在網(wǎng)面和垂直網(wǎng)面的運(yùn)動(dòng)。同時(shí),網(wǎng)衣將網(wǎng)箱的移動(dòng)作用到機(jī)器人上,改變機(jī)器人的狀態(tài);此外,網(wǎng)衣可能改變水動(dòng)力特性,進(jìn)而改變機(jī)器人的狀態(tài)。在本文模擬仿真中,假設(shè)網(wǎng)箱固定,重點(diǎn)關(guān)注摩擦力和支持力的影響 (圖14),忽略其他受力方面,這種簡(jiǎn)化有助于減少計(jì)算的復(fù)雜度。

    圖14 機(jī)器人在網(wǎng)箱表面的受力情況和貼網(wǎng)示意圖Fig. 14 Force situation of robot on surface of cage and schematic diagram of sticking net

    4 結(jié)論

    本文基于Gazebo 仿真平臺(tái),搭建了洗網(wǎng)機(jī)器人仿真環(huán)境,實(shí)現(xiàn)了洗網(wǎng)機(jī)器人三維建模和動(dòng)力學(xué)仿真。通過(guò)單元仿真的方式,分別建?;蚍抡媪送七M(jìn)器及洗盤、流體阻尼、慣量及附加質(zhì)量、海流場(chǎng)以及網(wǎng)衣對(duì)洗網(wǎng)機(jī)器人運(yùn)動(dòng)的影響。驗(yàn)證了仿真平臺(tái)的有效性,為后續(xù)設(shè)計(jì)抗干擾能力強(qiáng)、能耗低的控制器奠定了基礎(chǔ)。但目前的仿真平臺(tái)中洗網(wǎng)機(jī)器人在水下的真實(shí)運(yùn)動(dòng)情況仍有提升的空間,體現(xiàn)在:1) 附加質(zhì)量未考慮機(jī)器人的結(jié)構(gòu);2) 流場(chǎng)模型和現(xiàn)實(shí)有差別;3) 網(wǎng)衣作用缺乏有效模型;4) 未考慮推進(jìn)器時(shí)滯效應(yīng);5) 未基于真實(shí)實(shí)驗(yàn)數(shù)據(jù)辨識(shí)機(jī)器人的受力情況。在后續(xù)的工作中將針對(duì)上述問(wèn)題進(jìn)行優(yōu)化,盡可能地模擬洗網(wǎng)機(jī)器人在真實(shí)工況下的行為,為控制器設(shè)計(jì)提供數(shù)據(jù)支撐。

    猜你喜歡
    網(wǎng)衣網(wǎng)箱浮力
    遮蔽效應(yīng)下雙平面網(wǎng)衣和圓形網(wǎng)衣的水動(dòng)力性能
    基于CFD的剛性養(yǎng)殖網(wǎng)衣流場(chǎng)數(shù)值模擬及不確定度分析
    “浮力”知識(shí)鞏固
    我們一起來(lái)“制服”浮力
    浮力大小由誰(shuí)定
    遼墓出土網(wǎng)衣編法三例
    10kV配網(wǎng)箱變局放檢測(cè)案例分析
    電子制作(2017年24期)2017-02-02 07:14:45
    全國(guó)首例網(wǎng)箱豢養(yǎng)江豚繁育成功
    池塘小網(wǎng)箱養(yǎng)殖黃鱔
    遼墓出土銅絲網(wǎng)衣修復(fù)與復(fù)原報(bào)告
    99久久精品一区二区三区| 搞女人的毛片| 午夜福利18| 十八禁网站免费在线| 免费看十八禁软件| 欧美另类亚洲清纯唯美| 特级一级黄色大片| 国产成人系列免费观看| 中文字幕熟女人妻在线| 亚洲人与动物交配视频| 免费观看精品视频网站| 三级国产精品欧美在线观看 | 亚洲国产色片| 国产在线精品亚洲第一网站| 精品一区二区三区视频在线观看免费| 国产精品av视频在线免费观看| 欧美日韩乱码在线| av片东京热男人的天堂| 国产精品99久久久久久久久| 无遮挡黄片免费观看| 亚洲av日韩精品久久久久久密| 久久午夜亚洲精品久久| 日日夜夜操网爽| 精品国产超薄肉色丝袜足j| 国产免费av片在线观看野外av| 国产成人精品无人区| 女人被狂操c到高潮| 日韩欧美在线乱码| 久久人人精品亚洲av| 九九热线精品视视频播放| 操出白浆在线播放| 中文字幕av在线有码专区| 悠悠久久av| 88av欧美| 日韩有码中文字幕| 丁香六月欧美| 999久久久国产精品视频| 亚洲精品美女久久久久99蜜臀| 高清毛片免费观看视频网站| 天堂网av新在线| 午夜两性在线视频| 欧美黑人欧美精品刺激| 欧美日韩一级在线毛片| 国产成人精品久久二区二区免费| 久久人人精品亚洲av| 偷拍熟女少妇极品色| x7x7x7水蜜桃| 亚洲自偷自拍图片 自拍| 人人妻人人澡欧美一区二区| 一个人免费在线观看电影 | 国产真人三级小视频在线观看| 蜜桃久久精品国产亚洲av| 国内精品一区二区在线观看| 国产免费av片在线观看野外av| 成人午夜高清在线视频| 国内精品久久久久久久电影| 黑人操中国人逼视频| 男人舔奶头视频| 久久精品综合一区二区三区| aaaaa片日本免费| 精品久久久久久久人妻蜜臀av| 999久久久国产精品视频| 女生性感内裤真人,穿戴方法视频| 两性夫妻黄色片| 久久久久精品国产欧美久久久| 老司机深夜福利视频在线观看| 色在线成人网| 欧美黑人巨大hd| 精品99又大又爽又粗少妇毛片 | 最近最新中文字幕大全免费视频| 午夜成年电影在线免费观看| 网址你懂的国产日韩在线| 国内揄拍国产精品人妻在线| 亚洲成av人片免费观看| 国产精品野战在线观看| 嫩草影视91久久| 午夜福利免费观看在线| 99久久综合精品五月天人人| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产精品sss在线观看| 亚洲av电影在线进入| 三级毛片av免费| 看黄色毛片网站| 精品不卡国产一区二区三区| 精品人妻1区二区| 亚洲成人精品中文字幕电影| 亚洲国产精品合色在线| 三级男女做爰猛烈吃奶摸视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲欧美在线一区二区| 制服丝袜大香蕉在线| 精品熟女少妇八av免费久了| 国产久久久一区二区三区| 色综合婷婷激情| 国产午夜福利久久久久久| 18禁黄网站禁片免费观看直播| 亚洲国产欧美人成| 熟女人妻精品中文字幕| 国内精品久久久久精免费| 国产亚洲精品久久久久久毛片| 亚洲国产精品sss在线观看| 亚洲精品国产精品久久久不卡| 亚洲成av人片免费观看| 两个人看的免费小视频| 日韩精品中文字幕看吧| 亚洲美女黄片视频| 欧美一区二区国产精品久久精品| 久久久久九九精品影院| 久久伊人香网站| 9191精品国产免费久久| 97碰自拍视频| 国产真人三级小视频在线观看| 最近视频中文字幕2019在线8| a级毛片在线看网站| 国内精品久久久久久久电影| 伦理电影免费视频| 国产精品电影一区二区三区| 日韩国内少妇激情av| 一a级毛片在线观看| 精品久久久久久久末码| 国产成人av激情在线播放| 日韩人妻高清精品专区| 最近视频中文字幕2019在线8| 狂野欧美激情性xxxx| 国产毛片a区久久久久| 国产又色又爽无遮挡免费看| 亚洲欧美日韩高清专用| 一个人看视频在线观看www免费 | 长腿黑丝高跟| 美女 人体艺术 gogo| 亚洲无线观看免费| 精品国产美女av久久久久小说| 欧美成狂野欧美在线观看| 国产精品野战在线观看| 国产男靠女视频免费网站| 免费看美女性在线毛片视频| 国产又黄又爽又无遮挡在线| 日本五十路高清| 欧美激情在线99| 在线播放国产精品三级| www.自偷自拍.com| 51午夜福利影视在线观看| 欧洲精品卡2卡3卡4卡5卡区| 免费在线观看影片大全网站| 欧美zozozo另类| 国产高清videossex| 国产极品精品免费视频能看的| 亚洲精品粉嫩美女一区| 国产激情久久老熟女| 夜夜看夜夜爽夜夜摸| 麻豆国产97在线/欧美| 首页视频小说图片口味搜索| 身体一侧抽搐| 亚洲欧美日韩卡通动漫| 亚洲人成网站高清观看| 婷婷精品国产亚洲av在线| 欧美成人性av电影在线观看| 免费观看精品视频网站| 日韩 欧美 亚洲 中文字幕| 午夜免费激情av| 99久久精品一区二区三区| 成人三级黄色视频| 亚洲国产色片| 亚洲精品一卡2卡三卡4卡5卡| 深夜精品福利| 中出人妻视频一区二区| 一区二区三区激情视频| 国产伦人伦偷精品视频| 国内精品久久久久精免费| 日本与韩国留学比较| 男女视频在线观看网站免费| 亚洲中文日韩欧美视频| 在线免费观看的www视频| 可以在线观看的亚洲视频| 亚洲自偷自拍图片 自拍| 午夜免费观看网址| 两人在一起打扑克的视频| 网址你懂的国产日韩在线| 亚洲精品美女久久av网站| 桃色一区二区三区在线观看| 老汉色∧v一级毛片| 亚洲九九香蕉| 午夜福利成人在线免费观看| 亚洲 国产 在线| av视频在线观看入口| 丰满的人妻完整版| 一级作爱视频免费观看| 又紧又爽又黄一区二区| 偷拍熟女少妇极品色| 午夜久久久久精精品| 国内少妇人妻偷人精品xxx网站 | 国产午夜福利久久久久久| av在线天堂中文字幕| 日本在线视频免费播放| 日韩欧美精品v在线| 日韩高清综合在线| 日韩国内少妇激情av| 国产成年人精品一区二区| 久久久久性生活片| 每晚都被弄得嗷嗷叫到高潮| 国模一区二区三区四区视频 | 成人午夜高清在线视频| 日韩欧美在线乱码| 国产欧美日韩精品亚洲av| 免费看日本二区| 日本黄大片高清| 免费观看人在逋| 国产人伦9x9x在线观看| 在线视频色国产色| 美女cb高潮喷水在线观看 | 岛国在线免费视频观看| 国产亚洲精品综合一区在线观看| 波多野结衣高清作品| 欧美黄色淫秽网站| 国产乱人视频| 天天躁日日操中文字幕| 国产免费男女视频| 深夜精品福利| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人不卡在线观看播放网| 久久亚洲真实| 国产成人精品久久二区二区免费| 亚洲美女黄片视频| 香蕉国产在线看| 国产成+人综合+亚洲专区| 国产精品久久久久久精品电影| 狂野欧美白嫩少妇大欣赏| h日本视频在线播放| 一个人免费在线观看的高清视频| 亚洲第一电影网av| 国产三级中文精品| 一本一本综合久久| 亚洲专区国产一区二区| 亚洲av熟女| 成人一区二区视频在线观看| 久久亚洲精品不卡| 在线看三级毛片| 亚洲电影在线观看av| 国产伦人伦偷精品视频| 成人欧美大片| 日本 av在线| 18禁黄网站禁片午夜丰满| 国语自产精品视频在线第100页| 亚洲无线在线观看| 欧美最黄视频在线播放免费| 19禁男女啪啪无遮挡网站| 少妇熟女aⅴ在线视频| 久久精品影院6| 午夜精品在线福利| 欧美日韩综合久久久久久 | 亚洲人成电影免费在线| 熟女电影av网| 国产精品久久久久久人妻精品电影| 国产1区2区3区精品| 国产精品av视频在线免费观看| 中文字幕精品亚洲无线码一区| 麻豆成人av在线观看| 我要搜黄色片| 色综合婷婷激情| 色吧在线观看| 国产成人欧美在线观看| 波多野结衣高清无吗| 俺也久久电影网| 日本熟妇午夜| 国内毛片毛片毛片毛片毛片| 亚洲av片天天在线观看| 亚洲精品国产精品久久久不卡| 国内久久婷婷六月综合欲色啪| 亚洲真实伦在线观看| 一本一本综合久久| 国语自产精品视频在线第100页| 欧美日韩福利视频一区二区| 岛国在线观看网站| 国产黄色小视频在线观看| 日日摸夜夜添夜夜添小说| 婷婷六月久久综合丁香| 欧美日韩亚洲国产一区二区在线观看| 在线视频色国产色| 最新在线观看一区二区三区| 欧美日韩综合久久久久久 | 亚洲 欧美一区二区三区| 国产伦精品一区二区三区视频9 | 91老司机精品| 久久99热这里只有精品18| 国产亚洲欧美98| 色哟哟哟哟哟哟| 亚洲国产精品成人综合色| 国产熟女xx| 长腿黑丝高跟| 观看免费一级毛片| 国产91精品成人一区二区三区| 久久草成人影院| 久久欧美精品欧美久久欧美| 亚洲人成电影免费在线| 成人av一区二区三区在线看| 麻豆国产av国片精品| 变态另类成人亚洲欧美熟女| 69av精品久久久久久| 午夜日韩欧美国产| 视频区欧美日本亚洲| 男女那种视频在线观看| 在线免费观看的www视频| 成人亚洲精品av一区二区| 精品人妻1区二区| 精品一区二区三区视频在线 | 亚洲av免费在线观看| 亚洲男人的天堂狠狠| 神马国产精品三级电影在线观看| 欧美日韩福利视频一区二区| 老熟妇仑乱视频hdxx| 国产精华一区二区三区| 久久香蕉精品热| 久久亚洲真实| 国产精品99久久99久久久不卡| 极品教师在线免费播放| 国产亚洲精品av在线| 久久久久久国产a免费观看| 国产在线精品亚洲第一网站| 岛国视频午夜一区免费看| 日本在线视频免费播放| 国产免费av片在线观看野外av| 欧美一级毛片孕妇| 国产成人系列免费观看| 亚洲欧美精品综合久久99| 神马国产精品三级电影在线观看| 亚洲人成伊人成综合网2020| 亚洲精品在线观看二区| 欧美色欧美亚洲另类二区| 一区二区三区高清视频在线| 欧美成狂野欧美在线观看| 九九在线视频观看精品| 这个男人来自地球电影免费观看| 国产精品亚洲一级av第二区| 18禁裸乳无遮挡免费网站照片| 激情在线观看视频在线高清| 欧美zozozo另类| 又黄又粗又硬又大视频| 制服丝袜大香蕉在线| 中文字幕av在线有码专区| 久久香蕉国产精品| 国产激情偷乱视频一区二区| 国产一区二区激情短视频| 国产精品一区二区免费欧美| 欧美黄色片欧美黄色片| 亚洲在线观看片| 岛国视频午夜一区免费看| 制服人妻中文乱码| 国产精品免费一区二区三区在线| 国产aⅴ精品一区二区三区波| 他把我摸到了高潮在线观看| 久久人妻av系列| 免费av毛片视频| 日韩欧美在线二视频| 国产精品香港三级国产av潘金莲| 日韩欧美在线二视频| 一区福利在线观看| 在线播放国产精品三级| 精品乱码久久久久久99久播| 精品国产乱子伦一区二区三区| 五月玫瑰六月丁香| 国内精品久久久久精免费| 男女床上黄色一级片免费看| 亚洲欧美日韩高清在线视频| 亚洲av成人av| 美女被艹到高潮喷水动态| 久久精品影院6| 中文字幕久久专区| 久久久久久国产a免费观看| 日日摸夜夜添夜夜添小说| 床上黄色一级片| 国产激情欧美一区二区| 久久久久久国产a免费观看| 国产伦精品一区二区三区四那| 国产av不卡久久| 精品国产乱子伦一区二区三区| 男女床上黄色一级片免费看| 国产av麻豆久久久久久久| 久久香蕉国产精品| 亚洲国产日韩欧美精品在线观看 | 亚洲欧美激情综合另类| 亚洲欧美日韩无卡精品| 欧美另类亚洲清纯唯美| 他把我摸到了高潮在线观看| 精品久久久久久久久久免费视频| 长腿黑丝高跟| 日韩欧美在线乱码| 一个人免费在线观看电影 | 毛片女人毛片| ponron亚洲| 一进一出抽搐动态| 午夜福利18| 欧美3d第一页| 欧美在线黄色| 日本三级黄在线观看| 精品国产美女av久久久久小说| 一本精品99久久精品77| 五月玫瑰六月丁香| 一区二区三区国产精品乱码| 久久中文字幕一级| 成人高潮视频无遮挡免费网站| 欧美成人性av电影在线观看| 一级黄色大片毛片| 后天国语完整版免费观看| 国产一区二区三区视频了| 久久久色成人| 亚洲av成人不卡在线观看播放网| 国产单亲对白刺激| 久久亚洲精品不卡| 国产在线精品亚洲第一网站| 国产亚洲欧美在线一区二区| 国产伦在线观看视频一区| 国产黄色小视频在线观看| 日韩精品中文字幕看吧| 十八禁人妻一区二区| 国产一区二区三区视频了| 欧美日韩福利视频一区二区| 欧美日本亚洲视频在线播放| 精品99又大又爽又粗少妇毛片 | 美女被艹到高潮喷水动态| 淫秽高清视频在线观看| 超碰成人久久| 夜夜爽天天搞| 精品国产超薄肉色丝袜足j| 亚洲国产色片| 午夜影院日韩av| 久久草成人影院| 国产高清激情床上av| 亚洲欧洲精品一区二区精品久久久| 18禁国产床啪视频网站| 免费一级毛片在线播放高清视频| 黄色视频,在线免费观看| 在线十欧美十亚洲十日本专区| av视频在线观看入口| 国产精品1区2区在线观看.| www.www免费av| 999精品在线视频| 国产成人av教育| 久久人人精品亚洲av| 欧美极品一区二区三区四区| 午夜免费观看网址| 国产精品野战在线观看| 一个人免费在线观看电影 | 精品熟女少妇八av免费久了| 国产日本99.免费观看| 久久欧美精品欧美久久欧美| 精品日产1卡2卡| 国产成人系列免费观看| 级片在线观看| 久久久久久国产a免费观看| 搡老熟女国产l中国老女人| 亚洲av熟女| 黑人巨大精品欧美一区二区mp4| 一本一本综合久久| 大型黄色视频在线免费观看| 国产97色在线日韩免费| 国产免费av片在线观看野外av| 久久久精品大字幕| 俄罗斯特黄特色一大片| 亚洲精品一卡2卡三卡4卡5卡| 成年人黄色毛片网站| 男人舔女人的私密视频| 久久久久国产一级毛片高清牌| 在线永久观看黄色视频| 免费一级毛片在线播放高清视频| 精品免费久久久久久久清纯| 变态另类成人亚洲欧美熟女| 欧美日韩福利视频一区二区| 91av网一区二区| 亚洲欧洲精品一区二区精品久久久| 天堂√8在线中文| 国产熟女xx| 日本 欧美在线| 麻豆久久精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| 99re在线观看精品视频| 久久热在线av| 天天添夜夜摸| 亚洲专区国产一区二区| 一进一出抽搐动态| 色噜噜av男人的天堂激情| 久久人妻av系列| 亚洲无线在线观看| 国产乱人伦免费视频| tocl精华| 亚洲欧美日韩无卡精品| 观看免费一级毛片| 国产精品精品国产色婷婷| h日本视频在线播放| 最近在线观看免费完整版| 女人高潮潮喷娇喘18禁视频| 蜜桃久久精品国产亚洲av| 亚洲熟妇中文字幕五十中出| 欧美乱码精品一区二区三区| 老鸭窝网址在线观看| 亚洲中文字幕一区二区三区有码在线看 | 国产av一区在线观看免费| 哪里可以看免费的av片| 又黄又粗又硬又大视频| 嫁个100分男人电影在线观看| 日韩精品青青久久久久久| 精品午夜福利视频在线观看一区| 在线免费观看的www视频| 亚洲午夜理论影院| 国产精品久久视频播放| 日韩精品青青久久久久久| 国产激情久久老熟女| 99在线人妻在线中文字幕| 亚洲自偷自拍图片 自拍| 高清毛片免费观看视频网站| 国产精品影院久久| 日韩欧美在线二视频| 亚洲精品国产精品久久久不卡| 少妇熟女aⅴ在线视频| 欧美成狂野欧美在线观看| 村上凉子中文字幕在线| 久久热在线av| 日韩欧美免费精品| 动漫黄色视频在线观看| 国产单亲对白刺激| 国产高清视频在线观看网站| 一本综合久久免费| 禁无遮挡网站| 日韩欧美精品v在线| xxx96com| 91av网站免费观看| 999久久久精品免费观看国产| 99久久精品国产亚洲精品| 日韩欧美一区二区三区在线观看| 亚洲国产欧美网| 最近最新免费中文字幕在线| 男人舔奶头视频| 国内久久婷婷六月综合欲色啪| 欧美日韩中文字幕国产精品一区二区三区| 久久久久国产精品人妻aⅴ院| 国产欧美日韩一区二区精品| 中文字幕人成人乱码亚洲影| 精品国产乱子伦一区二区三区| 性色av乱码一区二区三区2| 精品一区二区三区视频在线观看免费| 欧美中文日本在线观看视频| 久久久久久久久久黄片| 国产精品亚洲美女久久久| 国产成人影院久久av| 2021天堂中文幕一二区在线观| 中文在线观看免费www的网站| a级毛片a级免费在线| 在线观看66精品国产| 国产麻豆成人av免费视频| 久久精品夜夜夜夜夜久久蜜豆| 看片在线看免费视频| 91久久精品国产一区二区成人 | 五月伊人婷婷丁香| 怎么达到女性高潮| 两个人视频免费观看高清| 久久久久久人人人人人| 国产久久久一区二区三区| 欧美性猛交黑人性爽| 又紧又爽又黄一区二区| 日本黄色片子视频| 亚洲国产欧洲综合997久久,| 国产毛片a区久久久久| 女警被强在线播放| 国产蜜桃级精品一区二区三区| 9191精品国产免费久久| 黄片大片在线免费观看| 久久人人精品亚洲av| 亚洲在线观看片| 欧美激情在线99| 欧美+亚洲+日韩+国产| 国产av一区在线观看免费| 久久中文字幕人妻熟女| 午夜福利在线观看免费完整高清在 | 露出奶头的视频| 99视频精品全部免费 在线 | 亚洲激情在线av| 国产精品亚洲美女久久久| 真人一进一出gif抽搐免费| 国产一区二区在线观看日韩 | 不卡av一区二区三区| 亚洲人成电影免费在线| 色综合亚洲欧美另类图片| 香蕉久久夜色| 国产欧美日韩一区二区精品| 亚洲美女视频黄频| 亚洲国产欧美网| 国产精华一区二区三区| 中文字幕人成人乱码亚洲影| 一区二区三区国产精品乱码| 一本久久中文字幕| 搡老岳熟女国产| 天堂av国产一区二区熟女人妻| 这个男人来自地球电影免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 成人欧美大片| www.www免费av| 搡老岳熟女国产| 一级毛片高清免费大全| 丝袜人妻中文字幕| 搡老岳熟女国产| 一区二区三区国产精品乱码| 男插女下体视频免费在线播放| 欧美在线黄色| 色综合婷婷激情| 成年人黄色毛片网站| 久久亚洲精品不卡| 一本综合久久免费| 欧美日韩福利视频一区二区| 91麻豆精品激情在线观看国产| 老司机深夜福利视频在线观看| 国产精品美女特级片免费视频播放器 | 两人在一起打扑克的视频| 亚洲av成人av| 在线a可以看的网站| 中出人妻视频一区二区| 久久久久久大精品|