亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        水下仿生側(cè)線感知研究進(jìn)展

        2023-03-18 01:51:10翟宇凡熊明磊謝廣明
        水下無人系統(tǒng)學(xué)報 2023年1期
        關(guān)鍵詞:側(cè)線游動機器

        翟宇凡 ,熊明磊 ,王 晨 ,謝廣明 ,3*

        (1.北京大學(xué) 工學(xué)院 智能仿生設(shè)計實驗室,北京,100871;2.北京大學(xué) 軟件工程國家工程研究中心,北京,100871;3.北京大學(xué) 海洋研究院,北京,100871)

        0 引言

        作為海洋大國,我國擁有豐富的海洋資源。自從黨的十八大報告中提出“建設(shè)海洋強國”戰(zhàn)略需求以來,國內(nèi)海洋裝備發(fā)展獲得前所未有的重視,特別是針對水下機器人的研究變得非?;钴S。

        隨著先進(jìn)制造與人工智能技術(shù)的飛速發(fā)展,機器人可以完成越來越多的陸地工程任務(wù),例如路徑規(guī)劃、自主巡航、資源探測與開采、群體協(xié)作等。相比之下,水下機器人面臨著極端復(fù)雜的工作環(huán)境,如復(fù)雜的水流運動與水底結(jié)構(gòu)、高壓環(huán)境、低能見度、海水的腐蝕性以及對儀器設(shè)備的干擾等,因此,發(fā)展水下機器人的難度遠(yuǎn)高于陸地機器人。一方面,水下機器人無法憑借簡單的探測工具獲取周圍水環(huán)境信息;另一方面,對機器人的運動控制也會受到水流的巨大干擾。這些都是在海洋探索與開發(fā)過程中亟待解決的問題。

        仿生機器人學(xué)是近年來的熱點研究問題之一。自然界中動植物的生理結(jié)構(gòu)或行為特征能夠為研究人員設(shè)計機器人提供靈感。魚類在海洋中生活了數(shù)十億年,他們的身體結(jié)構(gòu)不斷進(jìn)化,行為不斷地適應(yīng)著復(fù)雜的水下環(huán)境,最終呈現(xiàn)出如今的形態(tài)。在運動方面,魚類的游動具有機動性高、推進(jìn)效率高、能量效率高等特征。在感知方面,水生生物特有的一種感知器官——側(cè)線能夠幫助它們獲取水下的壓強與流速信息,進(jìn)而由神經(jīng)系統(tǒng)進(jìn)行行為決策。

        以墨西哥盲魚為例(見圖1),它們長期生活于黑暗的水下洞穴環(huán)境中,視覺系統(tǒng)退化,但是仍然能夠憑借側(cè)線以及其他感官系統(tǒng)實現(xiàn)捕食等行為[1]。受此啟發(fā),為了應(yīng)對水下環(huán)境中聲吶傳感器的散射和多路徑傳播、光學(xué)傳感器在黑暗渾濁的海洋環(huán)境中的失明等問題[2],研究者們開發(fā)了由一系列傳感器組成的人工側(cè)線系統(tǒng)作為新型傳感系統(tǒng),用于水下探測。相比于感知遠(yuǎn)場聲波的聲吶系統(tǒng),人工側(cè)線系統(tǒng)主要用于感知近場產(chǎn)生的局部流動,能夠彌補聲吶在近距離感知上的盲區(qū)。此外,隨著隱身技術(shù)的不斷發(fā)展,常規(guī)的傳感系統(tǒng)存在失效問題[3],而水下物體運動產(chǎn)生的流場變化幾乎是不可消除的,人工側(cè)線系統(tǒng)對于反隱身技術(shù)的發(fā)展有著重要意義。為了進(jìn)一步提升人工側(cè)線系統(tǒng)的性能,完成各種各樣的水下任務(wù),研究人員以仿生機器魚作為人工側(cè)線系統(tǒng)的載體,進(jìn)行了一系列感知、控制實驗。文中以水下仿生側(cè)線感知為主題,綜述國內(nèi)外研究現(xiàn)狀,并介紹筆者團隊在此領(lǐng)域的研究成果。

        圖1 墨西哥盲魚Fig.1 Mexican blind fish

        1 真實魚類側(cè)線感官

        側(cè)線是水生生物皮膚感官中最高度分化的組織結(jié)構(gòu),能夠幫助生物在水下黑暗環(huán)境中感知壓強與速度信息,獲取周圍環(huán)境特征,實現(xiàn)捕食、集群等行為[4]。側(cè)線主要通過神經(jīng)丘感受外界刺激,具體而言可以分為2 類: 表面神經(jīng)丘和管道神經(jīng)丘[5]。這2 種神經(jīng)丘有著不同的生理結(jié)構(gòu),且經(jīng)過億萬年的進(jìn)化,在魚體表面呈現(xiàn)出各自特定的數(shù)目與分布特征(見圖2[6]),能夠獲取不同類型的環(huán)境信息,幫助生物體達(dá)到最好的感知效果[7]。

        圖2 魚類側(cè)線神經(jīng)丘結(jié)構(gòu)與分布Fig.2 Structure and distribution of the neuromasts of the real fish lateral line

        表面神經(jīng)丘位于皮膚表面,具有一定的凸起結(jié)構(gòu),數(shù)目較多且較為分散,能夠幫助魚類感知水流流向與流速[8]。在物理模型中,表面神經(jīng)丘可以表示為2 段首尾相連且彎曲剛度不同的梁(見圖3[9])。當(dāng)生物體游動時,表面神經(jīng)丘的梁結(jié)構(gòu)在流速為U∞的水流作用下,沿z軸發(fā)生彎曲,撓度為v(z),位于皮膚下方的神經(jīng)丘產(chǎn)生相應(yīng)的神經(jīng)沖動,進(jìn)而傳遞到中樞神經(jīng)系統(tǒng)進(jìn)行處理,得到關(guān)于水流速度與方向的信息。表面神經(jīng)丘靈敏度由撓度v(z)與水流速度U∞的相對大小關(guān)系決定[9]。

        圖3 表面神經(jīng)丘簡化物理模型Fig.3 Simplified physical model of superficial neuromast

        管道神經(jīng)丘在體表一般呈直線分布,位于表皮下的側(cè)線管中,側(cè)線管內(nèi)充滿著粘液,并且與外界水環(huán)境通過小孔相通,能夠幫助魚類感知水下的壓強梯度[10]。具體的感知機理為: 外部流場的壓強差p1?p2通過小孔傳遞到側(cè)線管中,引起側(cè)線管中的局部流動,從而觸發(fā)管道神經(jīng)丘產(chǎn)生神經(jīng)沖動[11](見圖4[2])。

        圖4 管道神經(jīng)丘簡化物理模型Fig.4 Simplified physical model of canal neuromast

        在這2 類神經(jīng)丘的幫助下,魚類能夠充分地獲取外界流場的速度信息與壓強信息,再經(jīng)由大腦實現(xiàn)對于這2 類信號的處理,即使處于黑暗環(huán)境,也能夠較好地感知水流信息,獲取周圍障礙、生物的位置,在水下更好地生存[12]。

        2 基于不同感知原理的人工側(cè)線傳感器

        為了克服聲學(xué)和光學(xué)原理在水下復(fù)雜環(huán)境中的失效問題,研究者們以真實側(cè)線為靈感,基于不同的傳感原理制作了多種用于水下檢測的人工側(cè)線傳感器,以此為單元設(shè)計出人工側(cè)線傳感器陣列,開展多種水下檢測研究,并與仿生機器魚結(jié)合進(jìn)行感知與控制研究[13-14]。

        現(xiàn)有的人工側(cè)線傳感器主要采用以下傳感方式模擬魚類的側(cè)線感知功能: 壓阻效應(yīng)、壓電效應(yīng)、電容效應(yīng)、光學(xué)原理以及熱線式風(fēng)速儀原理等,部分人工側(cè)線傳感器如圖5 所示。

        圖5 基于不同傳感原理的人工側(cè)線傳感器Fig.5 Several artificial lateral line sensors based on different sensing principles

        2.1 壓阻效應(yīng)

        壓阻式傳感器主要基于半導(dǎo)體材料的壓阻效應(yīng)感知外部信息。半導(dǎo)體材料受到外力后發(fā)生形變,電阻也會產(chǎn)生變化。通過設(shè)計合適的測量電路,讀取有關(guān)電學(xué)信息,能夠反推出壓強、張力等待測的力學(xué)量,進(jìn)而得到水流的速度、壓強等信息。真實的側(cè)線由2 類神經(jīng)丘組成,以此為靈感設(shè)計的壓阻式人工側(cè)線傳感器也可分為2 類,分別能夠測量水下的流速與壓強。

        Fan 等[15]制造出第1 個壓阻式人工側(cè)線流速傳感器。該傳感器(見圖5(a))主要由1 根垂直的人造纖毛、1 個固定式的半導(dǎo)體自由懸臂和1 個應(yīng)變儀組成。類似于魚類側(cè)線的感知機理,外部水流流經(jīng)時,連接在懸臂自由端的纖毛產(chǎn)生彎曲,導(dǎo)致懸臂發(fā)生形變,位于懸臂梁底部的應(yīng)變儀測量出這一形變大小,反映水流信息。該傳感器能夠測量流速范圍為0.1~1 m/s 的層流。后續(xù)的研究者將這一傳感器與熱線式風(fēng)速儀進(jìn)行對比,省略了懸臂設(shè)計,將纖毛直接連接在基板上,更加貼近真實側(cè)線的結(jié)構(gòu)(見圖5(b))。這一改進(jìn)設(shè)計使得傳感器能夠用于不同溫度條件與多種復(fù)雜流動下的流速檢測,提高了傳感器的適用范圍與穩(wěn)定性,并且能夠大規(guī)模集成,用于分布式流量傳感[16]。此外,研究者們還從傳感器結(jié)構(gòu)、半導(dǎo)體材料選擇、制作工藝等角度進(jìn)行優(yōu)化(見圖5(d)(e)),提高了人工側(cè)線流速傳感器的靈敏度、測量范圍及魯棒性等性能[17-26]。

        人工側(cè)線流速傳感器主要基于懸臂或纖毛結(jié)構(gòu),而另一類用于測量壓強的人工側(cè)線傳感器主要采用平面結(jié)構(gòu)設(shè)計,將壓阻器件直接安裝在基板上,檢測水下壓強場的分布與變化。Vicente 等[27]提出一種用于檢測方形和圓柱形障礙物的人工側(cè)線壓強傳感器陣列(見圖5(c))。該陣列由數(shù)百個微型壓強傳感器組成,這些傳感器都被集成在半導(dǎo)體材料上。該陣列的壓強檢測極限為1 Pa。除了集成現(xiàn)有的壓強傳感器外,研究者們還自主設(shè)計了多種用于人工側(cè)線系統(tǒng)的壓強傳感器,并將其通過一定的排列方式組合成傳感陣列(見圖5(f)),用于檢測水下的壓強信息以及物體信息。通過使用液晶聚合物等新型材料,傳感器的化學(xué)穩(wěn)定性、密封性以及強度得到了大幅提升,能夠用于極其惡劣的深海環(huán)境[28-31]。

        2.2 壓電效應(yīng)

        壓電效應(yīng)是指特定材料在受到外力作用時表面上產(chǎn)生電荷,破壞原有電中性的現(xiàn)象。這種現(xiàn)象也為人工側(cè)線傳感器的研發(fā)提供了新的靈感?;诓牧蠅弘娦?yīng)的傳感器能夠通過收集電信息來感知來自外部環(huán)境的力學(xué)信息。Asadnia 等[32-33]制作了一款基于Pb(Zr0.52Ti0.48)O3薄膜材料的壓電式壓強傳感器(見圖5(g))。其采用立體光刻工藝,將聚合物纖毛連接在薄膜上,并且通過浮動電極來測量材料產(chǎn)生的壓電效應(yīng)大小,進(jìn)而得到纖毛感受的外部環(huán)境信息。并基于纖毛束結(jié)構(gòu)設(shè)計了一款微型流量傳感器(見圖5(h)),其主要部件為用于感知外部信息的聚合物細(xì)絲與壓電納米纖維。該團隊使用精密滴鑄和膨脹工藝,在纖毛束周圍形成了一個圓頂狀水凝膠殼斗,提升了傳感器的穩(wěn)定性。該傳感器的檢測閾值為8 μm/s[34]。此外,離子聚合物-金屬復(fù)合材料(IPMC)、聚偏二氟乙烯(PVDF)等壓電智能材料也能夠用于制作壓電式人工側(cè)線傳感器[35-36]。相比于壓阻式傳感器而言,壓電式傳感器元件能夠在外力作用下自主產(chǎn)生可測的電信號,不需要添加外部測量電路,因此體積較小,易于集成。

        2.3 電容效應(yīng)

        電容器具有高靈敏度和低功耗的特點,因此也被應(yīng)用到許多不同類型的傳感器中。電容式傳感器中最關(guān)鍵的組件是電容示值器,它能夠?qū)⑼獠看碳さ淖饔眯ЧD(zhuǎn)換為電容的變化,從而提供了一種檢測水下壓強和流速的有效方法。類似于壓阻式傳感器,附著在基板上的纖毛在局部流動的作用下發(fā)生形變,進(jìn)而改變基板電極間的距離,根據(jù)電容的變化信息就能夠得到外部流場的信息。Krijnen 等[37]受蟋蟀感知啟發(fā),將纖毛連接到氮化硅薄膜上,薄膜間的鉻電極形成了可變的電容器(見圖5(i))。該傳感器的結(jié)構(gòu)還可以進(jìn)一步優(yōu)化[38-39]。Stocking 等[40]受海豹觸須啟發(fā)制作了一種電容式晶須傳感器(見圖5(j))。將剛性人造晶須連接到錐型平行板電容器基座上,流場變化引起晶須形變,進(jìn)而導(dǎo)致電容的變化。電容式傳感器能夠用于測量流場信息,但是嚴(yán)格來說,其結(jié)構(gòu)和感知原理與魚類的側(cè)線并不相同。目前對于電容式傳感器的研究僅局限于傳感器的制作以及校準(zhǔn),并未有大量的實際應(yīng)用研究。

        2.4 光學(xué)原理與熱線式風(fēng)速儀原理

        光學(xué)原理與熱線式風(fēng)速儀原理也在一定程度上推動了人工側(cè)線傳感器的發(fā)展?;诠鈱W(xué)原理的傳感器具有類似于側(cè)線的纖毛結(jié)構(gòu),纖毛在水流中的彎曲可以通過紅外發(fā)光二極管等光學(xué)方法進(jìn)行測量,進(jìn)而推算水流速度(見圖5(k))[41-44]。熱線式風(fēng)速儀是一種常用的測量儀器,主要結(jié)構(gòu)為放置在流體中且被加熱過的金屬線,當(dāng)有流體流經(jīng)時,熱量的損失導(dǎo)致溫度和電阻的變化,進(jìn)而可以通過電信號測量流體速度。研究者們基于這一原理制作了用于測量流場速度的傳感器陣列(見圖5(l))[45-48]。這2 類傳感器雖然能夠測量流場速度,用于水下檢測研究,但是與魚類的側(cè)線結(jié)構(gòu)有著很大的差異,難以安裝到仿生機器魚載體上進(jìn)行感知實驗探究。

        2.5 小結(jié)

        目前大部分的研究主要關(guān)注壓阻式傳感器的設(shè)計與優(yōu)化。相比于其他傳感器,壓阻式傳感器在感知原理上最接近魚類的側(cè)線神經(jīng)丘,制作工藝更加成熟,商業(yè)產(chǎn)品化程度更高,穩(wěn)定性更好,因此,研究者們更多地會選擇基于壓阻效應(yīng)的壓強傳感器與流速傳感器進(jìn)行水流檢測研究,或者與仿生機器魚相結(jié)合?,F(xiàn)有的人工側(cè)線傳感器已經(jīng)能夠用于初步的水下探測研究,但是這些傳感器大多只是簡單模仿了魚類側(cè)線神經(jīng)丘的結(jié)構(gòu),并不能深度還原具體的感知機理,在信息處理性能方面還有很大的差距。

        為了更好地模擬真實的魚類側(cè)線,除了對于單個傳感器的設(shè)計之外,研究者們還對陣列中的傳感器數(shù)目、位置分布優(yōu)化進(jìn)行了深入的研究。針對各種流動問題(如均勻流、偶極子振蕩源產(chǎn)生的流動),基于勢流理論與保角變換計算特定翼型周圍的速度場與壓強場,從場的數(shù)據(jù)分布特征、系統(tǒng)經(jīng)驗可觀測性等角度進(jìn)行分析,提出一些描述傳感器性能的指標(biāo),如最優(yōu)測度、對流速以及目標(biāo)物體信息的估計誤差等,然后利用優(yōu)化方法得到最優(yōu)的傳感器數(shù)目以及位置分布,并且與真實魚類的側(cè)線分布進(jìn)行了對比[49-53]。這一類研究探究的問題較為簡單,首先,勢流理論模型沒有考慮真實流動中的粘性效應(yīng);其次,保角變換只能探究特定外形周圍的流動,較難進(jìn)行推廣;最后,對于一些較為復(fù)雜的非定常流動問題(例如魚體擺動前進(jìn)時的流場感知),這種方法將極為復(fù)雜。為了更好地模擬真實流場,研究者們通過計算流體動力學(xué)(computational fluid dynamics,CFD)仿真以及實驗獲取流動中的流場數(shù)據(jù),提出描述傳感器數(shù)據(jù)有效性、差異性和冗余性的指標(biāo),通過貝葉斯估計、神經(jīng)網(wǎng)絡(luò)等方法求解傳感器的最優(yōu)分布,得到更為可靠的結(jié)果[54-58](見圖6[56],圖中序號1~20 為傳感器最優(yōu)位置順序)。但是總體而言,對于這一類問題的探究,目前大多只停留在工程優(yōu)化的層面,不同的問題有著不同的優(yōu)化目標(biāo)以及不同的側(cè)線分布特征,尚未有一個從真實魚類側(cè)線生物機理與進(jìn)化角度出發(fā),能夠真正解釋側(cè)線分布、指導(dǎo)傳感器陣列設(shè)計的結(jié)論。

        圖6 壓強傳感器與剪應(yīng)力傳感器分布優(yōu)化結(jié)果Fig.6 Optimization results for the distribution of pressure sensors and shear stress sensors

        3 基于人工側(cè)線系統(tǒng)的水流檢測研究

        自然界中的水流十分復(fù)雜,并且上文提到現(xiàn)有的水下探測方法具有一定的局限性,隨著各種人工側(cè)線傳感器的出現(xiàn),由這些傳感器組成的人工側(cè)線系統(tǒng)逐漸被用于感知流體環(huán)境信息。目前這類研究主要集中在以下幾個方面: 流場分類、流速與流向檢測、渦街檢測和振蕩源檢測等。部分人工側(cè)線系統(tǒng)及其載體如圖7所示。

        1) 在流場分類方面,Salum?e 等[59]使用1 條安裝有5 個壓強傳感器(S1~S5)的機器魚(見圖7(a))開展流動分類實驗,發(fā)現(xiàn)在不同的流動環(huán)境下(例如均勻流動和周期性湍流),壓強傳感器數(shù)據(jù)呈現(xiàn)不同的特征,這一點能夠幫助機器魚識別所處的流場環(huán)境;此后又通過增加傳感器的數(shù)目(ps0~ps10)提升了人工側(cè)線系統(tǒng)的性能(圖7(b)分別為軸線AA'上的主視圖與俯視圖),新的人工側(cè)線系統(tǒng)能夠幫助機器魚基于壓強波動分布區(qū)分層流與湍流[60]。除了對于壓強數(shù)據(jù)的基本分析外,學(xué)習(xí)方法也能夠幫助人工側(cè)線系統(tǒng)識別流場環(huán)境。Liu 等[61]使用由23 個壓強傳感器組成人工側(cè)線系統(tǒng)(見圖7(c)),基于可視化的壓差矩陣與一個4 層卷積神經(jīng)網(wǎng)絡(luò)完成了流場分類任務(wù)。

        2) 水流最重要的參數(shù)為流速與流向,對這2個參數(shù)的準(zhǔn)確估計是精準(zhǔn)、穩(wěn)定控制機器魚運動的基礎(chǔ)。Salum?e 等[59]基于伯努利方程,擬合出機器魚(見圖7(a))兩側(cè)的平均壓強與水流速度的關(guān)系,同時基于機器魚受到水流沖擊的一側(cè)壓強更大,提出根據(jù)機器魚兩側(cè)壓強差識別流向的方法。此外,核嶺回歸以及神經(jīng)網(wǎng)絡(luò)等學(xué)習(xí)方法也被用來分析人工側(cè)線數(shù)據(jù)與水流流速、流向的關(guān)系[62-65]。

        圖7 多種人工側(cè)線系統(tǒng)及載體Fig.7 Several artificial lateral line systems and the carriers

        3) 渦街在自然界中十分常見,例如水流流經(jīng)巖石障礙物產(chǎn)生的尾流、魚類擺尾產(chǎn)生的尾跡等。對渦街特征的估計能夠幫助機器魚識別上游物體的尺寸、位置以及周圍魚類的相關(guān)信息,這對機器魚在自然水流環(huán)境下的感知與運動具有重要意義[66]。Yang 等[47]首先使用人工側(cè)線系統(tǒng)分析圓柱體產(chǎn)生的Kármán 渦街的空間速度分布。Ren等[67]從理論角度研究了真魚如何利用側(cè)線感知渦街,基于勢流理論,建立了魚體附近壓強場分布與渦街尺寸、傳播速度和方向等參數(shù)的聯(lián)系[41,59,68]。此外,渦街頻率、強度和距離等參數(shù)也被證明能夠通過人工側(cè)線數(shù)據(jù)來估計[69]。在此基礎(chǔ)上,Free等[70-71]通過經(jīng)驗觀測確定了游經(jīng)渦街的最優(yōu)路徑,采用閉環(huán)控制策略,實現(xiàn)了機器魚在渦街中的導(dǎo)航與巡游。

        4) 振蕩源也是一種模擬魚類尾跡的方法。利用人工側(cè)線系統(tǒng)獲取的速度與壓強信息,估計振蕩源的位置與運動方向,對于進(jìn)一步實現(xiàn)水下機器魚的最優(yōu)軌跡控制與跟蹤游動有著重要意義?;趧萘骼碚?Yen 等[72]將機器魚的尾鰭視為振蕩源(見圖7(d)),建立了流場變化的數(shù)學(xué)模型,描述了尾鰭擺動參數(shù)與人工側(cè)線壓強變化的關(guān)系,設(shè)計了讓機器魚根據(jù)人工側(cè)線信息調(diào)整擺動幅度、頻率和偏航角的跟蹤控制方法。除了流體力學(xué)理論模型外,研究者們還使用CFD 仿真和真實實驗的方式,基于其他人工側(cè)線模型(圖7(e)和(f))獲取相關(guān)數(shù)據(jù),圖7(e)中字母所示為傳感器位置。使用最小二乘法、波束成形算法以及神經(jīng)網(wǎng)絡(luò)等方法,較好地估計了振蕩源的參數(shù)[3,44,51,73-88],實現(xiàn)了對振蕩源運動軌跡的追蹤[89-90]。

        4 基于人工側(cè)線系統(tǒng)的仿生機器魚感知與控制研究

        水流檢測主要基于靜態(tài)的機器魚模型或者靜態(tài)的人工側(cè)線陣列,而對于自由運動的機器魚而言,使用運動的側(cè)線感知非定常流動時,難度將大大增加。然而,Chambers 等[91]的研究表明,運動狀態(tài)下的人工側(cè)線感知性能優(yōu)于靜態(tài)側(cè)線,這也與自然界中魚類的感知方式一致。因此,搭載有人工側(cè)線的、運動的機器魚能夠提供更有價值的實驗數(shù)據(jù),可以嘗試?yán)眠@些信息,讓機器魚實現(xiàn)類似于真實魚類的感知功能與自主決策行為。由于完全自由游動情形較為復(fù)雜,研究者們通常首先探究具有一定運動約束的側(cè)線感知問題。目前,相關(guān)領(lǐng)域已有較多的研究成果,例如運動模式識別、運動參數(shù)(速度、攻角等)估計與控制、障礙物識別與避障、鄰近機器魚的感知、定位以及能量節(jié)省等。這也是筆者團隊的主攻方向。部分仿生機器魚與人工側(cè)線系統(tǒng)如圖8所示。

        在運動模式識別方面,Liu 等[92]發(fā)現(xiàn)機器魚以不同的模態(tài)(直游、轉(zhuǎn)彎、上升和下潛等)運動時,壓強數(shù)據(jù)呈現(xiàn)出不同的特征。他們從9 個壓強傳感器數(shù)據(jù)中提取特征點,基于減法聚類算法識別出機器魚的運動模態(tài)。該方法有助于從側(cè)線數(shù)據(jù)中快速獲得機器魚的運動狀態(tài),有利于后續(xù)迅速、高效施加控制決策。

        在運動參數(shù)估計以及機器魚控制方面,Kruusmaa 等[93]基于壓強傳感器(見圖7(a))測得的數(shù)據(jù)給出速度估算公式,并且提出一種線性控制律來調(diào)節(jié)機器魚尾鰭擺動頻率,從而保持位置的穩(wěn)定。此外,還基于魚頭兩側(cè)的壓差估計魚頭與水流的夾角,并且設(shè)計控制方法讓機器魚根據(jù)壓差數(shù)據(jù)調(diào)整游動參數(shù),最終實現(xiàn)了機器魚在穩(wěn)態(tài)水流以及障礙物后方的位置估計以及位置穩(wěn)定[94]。美國馬里蘭大學(xué)研究團隊[95-97]設(shè)計了一種類似于Joukowski 翼型的柔性機器魚(見圖8(a),圖中物理量詳見文獻(xiàn)[95]),并在體表安裝分布式壓強傳感器,同時基于貝葉斯濾波器提出一種檢測流速和攻角的方法,實現(xiàn)了對機器魚游速的反饋控制。此外,相關(guān)研究表明,基于神經(jīng)網(wǎng)絡(luò)等訓(xùn)練方法,壓傳數(shù)據(jù)還能夠反映機器魚的俯仰角、擾動力和擾動力矩等參數(shù),有助于提升對機器魚的運動控制精度[98-100]。

        避障功能也是機器魚水下游動過程中必須具備的功能之一,這一方面也有著許多研究成果。Martiny 等[101]使用一款裝配有4 個熱線式人工側(cè)線傳感器的自主水下航行器(見圖8(b)),基于測得的航行器周圍的流速信息,通過理論分析與實驗方法估計出了障礙物與航行器的距離。DeVries等[96]對上游有障礙物時機器魚(見圖8(c))周圍的流場進(jìn)行分析,從理論的角度得到了估計障礙物相對位置的非線性模型,并且使用遞歸貝葉斯濾波實現(xiàn)了機器魚在障礙物后方的位置穩(wěn)定。Yen等[102]將機器魚(見圖8(d),圖中A,B,C 分別為傳感器位置)的尾鰭等效為一個振蕩偶極子源,基于3 個壓強傳感器測得的數(shù)據(jù)提出了一種控制機器魚沿著直線壁面游動并且保持一定距離的方法。此外,Yang 等[103]使用壓差傳感器與流速傳感器融合感知,使得機器魚(見圖8(e),圖中P1~P3為壓強傳感器,V4~V6為速度傳感器)在不同的攻角下都能夠估計出與壁面的距離?,F(xiàn)有的避障研究都基于最簡單的情形,目前尚未有能夠控制機器魚在復(fù)雜障礙物群之間自由游動的研究結(jié)果,距離真正的避障需求差距較大。

        圖8 多種仿生機器魚及人工側(cè)線系統(tǒng)Fig.8 Several bio-inspired robotic fish and artificial lateral line systems

        實際應(yīng)用中,希望能夠控制機器魚群完成單條機器魚無法完成的復(fù)雜任務(wù)。因此,除了對于單魚的感知、控制研究外,多個機器魚相互間的感知與編隊控制也是一個熱點問題。Li 等[104]基于仿真數(shù)據(jù)訓(xùn)練了人工側(cè)線壓強數(shù)據(jù)與雙魚距離的關(guān)系。Yen 等[105]用一個振蕩源模型模擬前魚的魚尾,基于壓傳數(shù)據(jù)估計出前魚的位置,并且設(shè)計了反饋跟蹤策略。Li 等[106]基于CFD 仿真與快速傅里葉變換(fast Fourier transform,FFT)方法,分析了雙魚前后距離、左右距離以及相位差等參數(shù)對于壓強信號的影響。筆者團隊[107-110]在這一方面也開展了大量研究,發(fā)現(xiàn)基于人工側(cè)線測得的壓強信息,后方的機器魚能夠感知到前魚的尾鰭擺動信息(例如頻率、幅度和偏航角)以及兩魚的相對位姿信息(例如相對垂直距離、相對偏航角、相對俯仰角和相對側(cè)傾角等)。為了更好地進(jìn)行感知,基于伯努利原理提出了雙魚游動時體表動壓變化的理論模型,并且采用可視化實驗和CFD 仿真對集群游動現(xiàn)象進(jìn)行了一定的分析。目前基于側(cè)線的集群游動研究只局限于數(shù)目較少的群體間的感知,尚未有基于側(cè)線數(shù)據(jù)的集群編隊控制方面的研究以及大規(guī)模群體的研究。

        除了上述研究主題外,人工側(cè)線也有助于實現(xiàn)機器魚的定位。Muhammad 等[111]提出了一種水下地理感知技術(shù),使得機器魚能夠根據(jù)側(cè)線數(shù)據(jù)在半自然或自然環(huán)境中識別出曾經(jīng)到達(dá)過的位置。Francisco 等[112]提出一種基于地圖的定位技術(shù),使用CFD 方法計算環(huán)境的流場速度分布圖,再根據(jù)人工側(cè)線系統(tǒng)獲取的信息估計所處的位置。此外,對于魚類、機器魚的研究還有推進(jìn)、節(jié)能等主題[59],由于上述研究中并不涉及人工側(cè)線感知,所以不再贅述。但是在未來,也可以從側(cè)線感知的角度探討這些問題。

        5 基于仿箱鲀機器魚的仿生側(cè)線感知研究進(jìn)展

        筆者團隊以自然界中的箱鲀(見圖9)為仿生對象,設(shè)計了2 款機器魚(見圖10)。箱鲀外形獨特,像一個盒子,盒子邊緣有獨特的脊骨結(jié)構(gòu),且皮膚表面呈現(xiàn)凹凸的特征。在游動時,這種結(jié)構(gòu)能幫助箱鲀產(chǎn)生有利于姿態(tài)穩(wěn)定的渦[113-115]。此外,盒子狀的機器魚能夠在體內(nèi)搭載更多用于感知、控制的硬件設(shè)備,有利于研究的開展。

        圖9 真實的箱鲀魚Fig.9 Real boxfish

        圖10 2款仿箱鲀機器魚Fig.10 Two box-fish-inspired robots

        5.1 2 款仿生機器魚及其動力學(xué)建模

        對于2 款仿生機器魚,其中一款外形接近箱鲀的特征,由多鰭肢驅(qū)動(尾鰭、胸鰭等)(見圖10(a))。在另一款仿生機器魚的設(shè)計過程中,考慮到上述提及的箱鲀游動時產(chǎn)生的渦結(jié)構(gòu)以及胸鰭的擺動會影響體側(cè)側(cè)線的感知效果,將機器魚殼體簡化為規(guī)則盒狀,并且去除了其他鰭肢,采用單尾鰭驅(qū)動(見圖10(b))。2 款機器魚分別搭載有9 個和11 個壓強傳感器作為人工側(cè)線系統(tǒng),以及慣性導(dǎo)航模塊(inertial measurement unit,IMU)測量三維空間內(nèi)的加速度與姿態(tài)角信息。2 款機器魚都基于中樞模式發(fā)生器(central pattern generator,CPG)來控制鰭肢的擺動,可以調(diào)整鰭肢的擺動頻率、幅度和偏置,實現(xiàn)直游、轉(zhuǎn)彎、上升、下降和橫滾等運動。

        在開展感知與控制研究之前,首先需要對于機器魚在水下的運動特性進(jìn)行分析。針對2 款機器魚建立了自由游動的三維動力學(xué)模型,該模型能夠在不同的運動方式及大范圍的運動參數(shù)下有效估計機器魚在水下的速度、姿態(tài)和軌跡等信息(見圖11[116-117])。

        圖11 盤旋運動下基于動力學(xué)模型的軌跡估計結(jié)果Fig.11 Trajectory estimation results based on dynamic model under spiral motion

        5.2 基于人工側(cè)線系統(tǒng)的運動狀態(tài)估計

        運動狀態(tài)估計是仿生機器魚研究中的基本問題之一,只有實現(xiàn)了準(zhǔn)確的運動速度、軌跡估計,才能更好地對仿生機器魚進(jìn)行控制,開展下一步研究。筆者團隊基于2 款仿生機器魚,在運動參數(shù)估計方法上取得了一定的進(jìn)展。

        機器魚的游動速度是最重要的運動參數(shù)之一,現(xiàn)有的速度估計方法大多基于伯努利方程,機器魚頭部的靜態(tài)壓強、動態(tài)壓強以及游動速度有如下關(guān)系

        式中:pt為運動時頭部的駐點壓強;ps為靜止時的壓強;v為運動速度;ρ為水流密度。在此基礎(chǔ)上,發(fā)現(xiàn)機器魚游動時身體會進(jìn)行節(jié)律性擺動導(dǎo)致駐點位置的偏移,因此在速度估計模型中引入了角速度ω,有

        式中:p為魚頭處3 個壓強傳感器減去靜水壓強后的均值;a,b,c,d為待定系數(shù)。根據(jù)回歸結(jié)果能夠基于壓傳數(shù)據(jù)與角速度反解游動速度,準(zhǔn)確、實時地估計機器魚的游動速度[118](見圖12(a)[119])。

        此外,除了游動速度,還需要考慮機器魚的轉(zhuǎn)彎角速度、上升/下降速度,才能實現(xiàn)對機器魚軌跡的準(zhǔn)確估計?;贚ighthill 的理論分析結(jié)果與非定常伯努利方程,首次推導(dǎo)出機器魚體表動態(tài)壓強變化量?pdynamic(絕對壓強減去靜水壓強)與主要運動參數(shù)的關(guān)系模型

        式中:U表示游動速度;ω表示偏航角速度;θ表示俯仰角速度;每一項前的系數(shù)都是只依賴于機器魚幾何外形的參數(shù)?;趯嶒灉y得的數(shù)據(jù),采用最小二乘法標(biāo)定出適用于不同運動狀態(tài)(直線、轉(zhuǎn)彎、上升和盤旋)下的模型參數(shù)。標(biāo)定完成后,在后續(xù)的實驗中,基于側(cè)線壓傳測得的壓強數(shù)據(jù),反解出機器魚的運動參數(shù),最終實現(xiàn)對軌跡的準(zhǔn)確估計(見圖12(b)[120])。

        圖12 運動狀態(tài)估計結(jié)果Fig.12 Estimation results of the motion states

        該工作首次實現(xiàn)了對于自由游動機器魚的軌跡估計,且方法能夠推廣至所有擁有剛性殼體的機器魚,也證明了人工側(cè)線系統(tǒng)在運動參數(shù)、軌跡估計方面的有效性,并具備成為水下機器人通用傳感系統(tǒng)基本元件的潛力。在未來,會關(guān)注壓強傳感器在野外環(huán)境中的標(biāo)定,并且嘗試將人工側(cè)線系統(tǒng)應(yīng)用于自然環(huán)境中的軌跡估計研究。此外,還會將人工側(cè)線系統(tǒng)得到的軌跡信息作為閉環(huán)控制回路中的反饋信息,實現(xiàn)對水下機器魚的軌跡跟蹤控制。

        5.3 基于人工側(cè)線系統(tǒng)的姿態(tài)保持控制

        攻角即機器魚朝向與來流的夾角。為了能夠控制機器魚像真實魚類一樣在水下自由游動,除了需要控制機器魚的游動速度外,還需要控制機器魚在來流中的攻角,實現(xiàn)姿態(tài)的穩(wěn)定。在此基礎(chǔ)上,從傳感器融合感知的角度進(jìn)一步研究這一問題。對比基于3 種反饋信息(人工側(cè)線與神經(jīng)網(wǎng)路估計的攻角數(shù)據(jù)、IMU 測得的偏航角數(shù)據(jù)以及二者融合數(shù)據(jù))下攻角控制策略的效果,可以發(fā)現(xiàn),雙傳感器融合信息方法優(yōu)于單傳感器方法(見圖13[121])。

        圖13 單尾鰭驅(qū)動機器魚基于雙傳感器融合方法的攻角控制結(jié)果Fig.13 Control results of angle of attack of single-tail driven robotic fish based on dual-sensor fusion method

        為了更好地解決這一問題,首次引入強化學(xué)習(xí)框架。由于流體的強非線性,機器魚在游動時周圍的流場變化無法用簡單的理論模型表達(dá)。為此,采用深度神經(jīng)網(wǎng)絡(luò)模擬機器魚游動時的流固耦合環(huán)境?;谶@一虛擬仿真環(huán)境不斷生成數(shù)據(jù),使用確定性策略梯度(deep deterministic policy gradient,DDPG)算法訓(xùn)練機器魚利用壓強傳感器數(shù)據(jù)與IMU 數(shù)據(jù)實現(xiàn)姿態(tài)穩(wěn)定的控制策略。最終,強化學(xué)習(xí)框架能夠直接將人工側(cè)線系統(tǒng)測得的傳感器原始數(shù)據(jù)以端到端的方式直接映射到連續(xù)控制信號,實現(xiàn)對攻角的控制(見圖14[122])。由于虛擬環(huán)境的引入,避免了大量獲取實驗數(shù)據(jù)的過程。相比于傳統(tǒng)的模型預(yù)測控制(model predictive control,MPC)方法而言,該方法能夠與MPC 方法相媲美甚至效果更好。此外,MPC 方法依賴大量的計算資源,需要在電腦上運行部署,而基于學(xué)習(xí)的方法能夠直接在機器魚搭載的樹莓派上運行,這對于自主式機器魚的感知與控制有著重要意義[122]。

        圖14 單尾鰭驅(qū)動機器魚基于深度強化學(xué)習(xí)框架的攻角控制結(jié)果Fig.14 Control results of angle of attack of single-tail driven robotic fish based on deep reinforcement learning framework

        這個方法能夠很好地用于實驗室環(huán)境下的姿態(tài)保持控制,但是自然環(huán)境下的實驗具有更多的不確定性。當(dāng)機器魚能夠自由游動時,控制其相對于水流的攻角也會更加復(fù)雜。未來,將嘗試結(jié)合速度、攻角的估計與控制研究,在自然環(huán)境中讓機器魚像真實魚類那樣游動。

        5.4 基于人工側(cè)線系統(tǒng)的機器魚鄰近感知

        目前關(guān)于仿生機器魚群的感知與控制研究較少,筆者團隊較早地基于人工側(cè)線系統(tǒng)開展了對于多魚鄰近感知的探究分析。一方面,通過可視化與CFD 仿真技術(shù)分析了機器魚擺尾產(chǎn)生的渦街特性(見圖15[107])。另一方面,通過采集尾鰭擺動引起的流場壓強變化數(shù)據(jù),發(fā)現(xiàn)人工側(cè)線壓強傳感器能夠檢測到前方尾鰭的擺動幅度、頻率、偏置角以及相對垂直高度、相對偏航角、相對俯仰角以及相對橫滾角等相對位姿信息[107]。在此基礎(chǔ)上,基于線性回歸、隨機森林和神經(jīng)網(wǎng)絡(luò)等方法通過壓強傳感器數(shù)據(jù)對上述信息進(jìn)行估計,得到了較好的估計效果(見圖16[108])。

        圖15 尾鰭后方瞬時渦環(huán)圖Fig.15 Instantaneous vortex core region behind an individual oscillating caudal fin

        圖16 基于隨機森林方法的前方尾鰭擺動幅度估計結(jié)果Fig.16 Estimation results of oscillating amplitude of the front caudal fin based on random forest method

        此外,還分析了領(lǐng)航者-跟隨者的雙魚隊形下,人工側(cè)線壓強數(shù)據(jù)與前后距離的關(guān)系[109],如圖17[110]所示。圖中: Δd1和Δd2為器件尺寸;Slogitudinal和Slateral為機器魚前后、左右距離;Vu和Vd為兩魚速度;ωu和ωd為兩魚角速度。進(jìn)一步地,基于勢流理論與數(shù)據(jù)驅(qū)動方法,提出了雙魚隊形下的壓強理論模型,發(fā)現(xiàn)運動兩魚相對一側(cè)的傳感器會受到兩魚運動的影響,而兩魚相背一側(cè)的傳感器只會受到自身運動的影響。實驗結(jié)果也驗證了這一模型的準(zhǔn)確性[110]。

        圖17 雙魚領(lǐng)航者-跟隨者隊形示意圖Fig.17 Leader-follower formation diagram of two robotic fishes

        以上進(jìn)行的機器魚鄰近感知研究大部分都是簡化的問題,機器魚的部分自由度被約束。目前尚未有關(guān)于完全自由的集群游動時機器魚對周圍個體的感知研究。進(jìn)行此項研究的基礎(chǔ)是首先需要基于一些方法讓機器魚以特定的隊形穩(wěn)定地進(jìn)行游動,目前尚未有基于側(cè)線反饋信息實現(xiàn)隊形控制的研究,而能夠?qū)崿F(xiàn)隊形保持的機器魚模型也并未搭載人工側(cè)線。未來可以將這兩方面研究相結(jié)合,將集群控制與鄰近感知研究拓展到完全自由的情形,讓機器魚群能夠像魚群一樣進(jìn)行游動與感知。

        6 總結(jié)與展望

        由于水下探測的需要以及現(xiàn)有水下探測方法的局限性,研究人員以魚類的側(cè)線感官為靈感,開發(fā)了由一系列傳感器組成的人工側(cè)線系統(tǒng),并與水下機器人技術(shù)相結(jié)合,實現(xiàn)了多項感知以及控制任務(wù)。人工側(cè)線系統(tǒng)在水下探測領(lǐng)域展現(xiàn)出巨大的潛力。

        文中概述了基于不同原理的人工側(cè)線傳感器及其用于水流檢測、機器魚感知與控制的各項研究。關(guān)于人工側(cè)線系統(tǒng)的研究已經(jīng)有了許多成果,但是現(xiàn)有的人工側(cè)線系統(tǒng)性能仍然與真實的魚類側(cè)線相去甚遠(yuǎn)。

        從傳感器性能的角度分析,目前用于人工側(cè)線系統(tǒng)的傳感器陣列從靈敏度、穩(wěn)定性、協(xié)調(diào)性和信息處理能力等許多方面來看,都無法與真實魚類的側(cè)線相比。為了提升單個傳感器的靈敏度與穩(wěn)定性,可以從大量新型材料中選取具有傳感能力、力學(xué)性能優(yōu)良、化學(xué)性質(zhì)穩(wěn)定的材料(例如石墨烯等)作為基礎(chǔ),設(shè)計相應(yīng)的敏感元件、轉(zhuǎn)換元件與處理電路,并作防水處理,使得傳感器能夠在惡劣的溫度和壓力條件下保持足夠的精準(zhǔn)度。除了提升單個傳感器的性能外,傳感器陣列的協(xié)調(diào)作用也十分重要。現(xiàn)有的人工側(cè)線系統(tǒng)一般只由一種傳感器(壓強傳感器或流速傳感器)組成,并且呈直線或方形排列。相比之下,真實的側(cè)線由表面神經(jīng)丘和管道神經(jīng)丘組成,并且具有復(fù)雜的分布方式。在后續(xù)的研究中,可以在機器魚表面同時安裝壓強傳感器和流速傳感器,進(jìn)行流場感知與機器魚控制實驗,與現(xiàn)有的傳感器陣列性能進(jìn)行對比分析。目前用來優(yōu)化人工側(cè)線分布的指標(biāo)大多只能針對某一特定的目標(biāo)問題,具有一定的局限性,并不能解釋魚類側(cè)線的分布規(guī)律,也無法推廣到所有的感知任務(wù)中。因此,需要研究魚類側(cè)線分布的生物學(xué)原理,提出具有普適意義的傳感器性能指標(biāo),從而為不同任務(wù)場景下人工側(cè)線系統(tǒng)的分布設(shè)計提供理論指導(dǎo)。在對側(cè)線信息的處理方面,魚類能夠迅速對感受到的信號做出反應(yīng)。隨著側(cè)線硬件的進(jìn)步,也需要開發(fā)更高效的濾波、感知以及行為決策方法,從算法的角度提升人工側(cè)線的性能。除了傳統(tǒng)的算法外,機器學(xué)習(xí)方法與人工智能的發(fā)展為這一問題提供了新的方向。

        在水流檢測方面,研究者們針對不同的問題提出了多種感知方法。但是,目前的研究大多在實驗室環(huán)境下進(jìn)行,人工側(cè)線載體運動簡單(例如靜止或直線運動),水環(huán)境穩(wěn)定。而在自然環(huán)境中,真實的水流與魚類的游動都要復(fù)雜的多,這將大大增加感知實驗的難度。隨著人工側(cè)線性能的提升,未來需要更多地開展自然環(huán)境中的實驗,例如使用人工側(cè)線系統(tǒng)在真實的河流中實時估計水流速度、河流中的障礙物位置以及真魚的位置與擺動狀態(tài)等,進(jìn)一步挖掘人工側(cè)線的感知潛力,使其越來越貼近真實的魚類側(cè)線功能。

        水下機器魚的流場輔助控制是一個前景十分可觀的研究方向,未來在海洋勘探中可以得到廣泛的應(yīng)用。作為一切控制任務(wù)與感知任務(wù)的前提,機器魚首先需要能夠準(zhǔn)確估計自身的速度、位置和姿態(tài)等信息。因此,需要基于流體力學(xué)理論或數(shù)據(jù)驅(qū)動方法建立水下機器魚的運動學(xué)與動力學(xué)模型,幫助機器魚在無法使用視覺的情況下完成運動模式的自我識別與軌跡估計。在此基礎(chǔ)上,才能夠設(shè)計感知與控制方法,豐富機器魚在水下的功能。避障行為是實現(xiàn)機器魚自主運動必須要解決的難題之一。機器魚需要根據(jù)人工側(cè)線獲取的流場信息,識別障礙物的大小和位置,重新構(gòu)建出周圍的地圖環(huán)境并且規(guī)劃最佳導(dǎo)航路徑。在實驗室環(huán)境中,可以探究機器魚對于簡單障礙物單元(例如墻壁、圓形障礙和三角障礙)以及障礙物群的感知,設(shè)計相應(yīng)的避障控制策略。在具備一定的感知能力后,在自然河流實驗環(huán)境中進(jìn)一步實現(xiàn)機器魚的避障游動行為。除了對于單條機器魚的控制研究外,基于人工側(cè)線系統(tǒng),還可以研究機器魚群的相互感知與集群控制,例如在三角隊形或方形隊形下,探究機器魚利用人工側(cè)線感知周圍個體的位置信息,并且基于相對位置信息設(shè)計反饋控制策略,實現(xiàn)隊形的穩(wěn)定。在未來的水下任務(wù)中,相比于單魚而言,使用魚群能夠極大地提高水下探索的效率與水下任務(wù)的成功率。

        7 結(jié)束語

        水下仿生側(cè)線感知技術(shù)涉及生物學(xué)、機器人學(xué)、流體力學(xué)以及智能科學(xué)等多個學(xué)科與領(lǐng)域,雖然研究難度較大,但具有較高的研究價值與廣闊的應(yīng)用前景。人工側(cè)線系統(tǒng)的發(fā)展為海洋勘探與開發(fā)提供了強大的工具,尤其是在惡劣環(huán)境下聲吶、紅外等傳感方式失效時。文中主要介紹了水下仿生側(cè)線感知的研究現(xiàn)狀以及筆者團隊的有關(guān)工作。目前現(xiàn)有的研究主要集中于人工側(cè)線傳感器的研發(fā)、基于人工側(cè)線系統(tǒng)的水流檢測、基于人工側(cè)線系統(tǒng)的仿生機器魚感知與控制等方面。但是目前的大部分研究對問題作了簡化,與真實魚類在水中的游動以及感知還有一定差距,要實現(xiàn)與真實魚類相媲美的感知能力與自主控制決策能力,還有很長的路要走。在未來,隨著硬件與感知決策算法的進(jìn)一步發(fā)展,人工側(cè)線系統(tǒng)有望進(jìn)一步提升水下機器魚的感知能力以及其他綜合能力,幫助更好地探索海洋、保護海洋,建設(shè)海洋強國。

        猜你喜歡
        側(cè)線游動機器
        永不停歇的魚
        機器狗
        機器狗
        球軸承用浪型保持架徑向游動量的測量
        哈爾濱軸承(2021年1期)2021-07-21 05:43:16
        把手放進(jìn)袋子里
        未來機器城
        電影(2018年8期)2018-09-21 08:00:06
        關(guān)于側(cè)線股道LKJ數(shù)據(jù)的編制
        捉迷藏
        無敵機器蛛
        父親
        天津詩人(2014年4期)2014-11-14 19:05:52
        久久亚洲精品国产亚洲老地址| 久久久精品免费观看国产| 天天躁日日操狠狠操欧美老妇| 加勒比日本东京热1区| 国产麻豆极品高清另类| 亚洲国产精品无码一线岛国| 亚洲精品午夜无码电影网 | 正在播放亚洲一区| 亚洲中文有码字幕青青| 久久精品国产72国产精福利| 人妻少妇被猛烈进入中文| 夜夜躁日日躁狠狠久久av| 一本无码人妻在中文字幕免费| 亚洲国产精品自拍一区| 精品国产亚洲av成人一区| 亚洲精品一区二区高清| 久久精品国产免费观看| 国产欧美日韩久久久久| 久久天堂av色综合| 久久精品国产亚洲av豆腐| 97人伦影院a级毛片| 一道久在线无码加勒比| 岛国视频在线无码| 亚洲最大在线视频一区二区| 亚洲熟女乱综合一区二区| 亚洲成a人片在线网站| 国产精品天干天干在线观蜜臀| 精品久久人妻一区二区| 国产高清在线精品一区app| av蓝导航精品导航| 亚洲红杏AV无码专区首页| 午夜一区二区视频在线观看| 一夲道无码人妻精品一区二区| 中文乱码人妻系列一区二区| 亚洲中文字幕无线乱码va| 穿着白丝啪啪的av网站| 国产成人无码a区在线观看视频 | 亚洲免费观看在线视频| 中文字幕亚洲精品人妻| 亚洲一区二区三区少妇| 中国人妻被两个老外三p|