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

    水下航行器對接過程下的動力學建模與運動仿真研究

    2023-09-22 01:48:32李德軍楊申申王小慶
    船舶力學 2023年9期
    關鍵詞:航速攻角航行

    周 超,張 偉,李德軍,王 璇,楊申申,王小慶

    (1.中國人民解放軍92228部隊,北京 100072;2.中國船舶科學研究中心,深海載人裝備國家重點實驗室,江蘇 無錫 214082;3.深海技術科學太湖實驗室,江蘇 無錫 214082)

    0 引 言

    水下航行器作為一種海洋探測裝備,其體積小、隱身性好、機動性強、成本低、可快速組網,探測能力強,可廣泛用于海底科學考察、水域勘察與識別、海洋生物群落調查、水下線纜維修和海底打撈等。水下航行器在長時間航行作業(yè)時,續(xù)航能力直接影響水下航行器探測能力,需要搭載在母艇上,以補充能量或者交換數據[1]。

    世界范圍內開展水下無人航行器研發(fā)的國家主要為美、俄、歐亞等國,其中,美國在水下無人航行器技術領域始終處于領先地位。目前,許多強國重視無人航行器與水下平臺的搭載。例如,瑞典的SAABSUBROV 航行器,可由潛艇進行操控;美國的“曼塔”(Manta)UUV、俄羅斯的“波塞冬”UUV 都計劃搭載于潛艇,輔助潛艇執(zhí)行水下作戰(zhàn)任務;美國新建第三批“弗吉尼亞”級核潛艇也采用了新型艇首結構,為多個無人潛水器的搭載提供了平臺基礎[2]。

    在近距離對接過程中,水下航行器和母平臺間存在復雜的水動力干擾。由于水下航行器相對于平臺尺寸較小,其航行姿態(tài)受影響較大,給對接回收帶來了困難。在對接過程中,需要對水下航行器進行不斷控制,以抵消水動力干擾,保證對接的精準完成。為了有效控制水下航行器,需要弄清對接過程中的水動力干擾作用機理與規(guī)律,為實現對接提供技術支撐。在水下航行器對接水動力干擾研究方面,文獻[3]基于CFD技術,研究了潛艇布放回收UUV過程中的相關參數對流場和水動力性能的影響。文獻[4]采用CFD技術,計算了大型水下搭載平臺與微小型水下機器人間的水動力干擾性能,分析了在不同來流速度和不同間距比時的微小型水下機器人的阻力、升力和傾覆力矩性能,此外,對微小型水下機器人的中縱剖面的壓力和速度分布進行了細致的分析。文獻[5]采用CFD 技術,采用有限體積法和非結構動網格技術,對AUV 對接過程在不同雷諾數下和不同攻角下的非定場流場進行了模擬研究。

    水下航行器對接過程主要分為兩步:當水下航行器距離目標平臺距離較遠時,根據低頻聲納的信號,確定目標平臺的相對位置和方位,通過遠程導引,航行靠近目標平臺,并到達平臺的垂直面后上方;當水下航行器距離目標平臺距離較近時,根據高頻聲納或者光學導引方式,開啟自動控制功能,完成水下航行器與平臺的對接。其中,近距離雙體干擾下的對接建模與控制是研究的難點。

    本文采用CFD 數值計算,利用RANS 方法求解湍流方程,研究干擾水動力特性及規(guī)律,完成水下航行器和母平臺不同距離下的水動力干擾計算,進而建立水下航行器對接時的動力學模型。在此基礎上,針對不同航速和不同距離下的對接工況,開展干擾水動力對姿態(tài)和速度影響的仿真研究;選擇干擾影響最大的工況,開展對接運動控制研究,基于動力學模型,采用PID 和滑??刂品椒?,完成對接控制器設計,開展對接運動仿真。

    1 水下航行器運動學

    1.1 運動坐標系

    為研究小型水下航行器對接時的運動特性,首先應定義O-XYZ為大地坐標系,其固定于地球;定義o1-x1y1z1為水下航行器坐標系,坐標原點為水下航行器的重心;定義o2-x2y2z2為平臺坐標系,坐標原點為母平臺的重心(如圖1所示)。水下航行器相對母平臺尺寸較小,其運動受干擾水動力影響較大,而母平臺受干擾水動力影響可以忽略,因此本文主要研究對接過程下的水下航行器運動建模問題。

    圖1 對接空間坐標系Fig.1 Spatial coordinate system under docking

    1.2 運動學模型

    水下航行器隨體坐標系與慣性系轉移矩陣[6]為

    水下航行器速度關系為

    水下航行器角速度關系為

    其中,

    上式中,X、Y、Z為水下航行器位置,u、v、w為線速度,p、q、r為角速度,φ、θ、ψ為姿態(tài)角。

    2 水下航行器對接動力學模型

    2.1 對接干擾力模型

    當水下航行器間進行對接作業(yè)時,兩者存在相互干擾的水動力。本文采用商業(yè)軟件STARCCM+,選用SSTk-ω模型進行數值計算分析,該模型的適用性在文獻[7]中得到了驗證,表明其在近壁面流動分離以及尾流場特性的數值預報上具有較好的適用性。

    目前通常采用三維不可壓縮粘性Navier-Stokes(N-S)方程作為流體運動的控制方程,求解水下運載平臺周圍的流場流動,其表達式為

    式中,u是流體速度矢量,μ是流體動力粘性系數,p是壓力,ρ為流體密度。

    首先,開展數值方法的驗證工作,選擇風洞模型試驗結果與計算結果進行比對,不同攻角的水下航行器流體動力系數計算結果和模型試驗結果如圖2所示。

    圖2 數值方法驗證結果Fig.2 Validation results of numerical method

    從圖中可以看出,本文計算結果與模型試驗結果吻合較好,相關水動力系數值最大偏差為5.65%,證明了計算方法的有效性。

    本文以Suboff為對接母平臺,選取長方體作為計算域(見圖3)。計算域的長度為Suboff 艇長的6倍,其中前方的面距離Suboff首部1.5倍艇長,后方的面距離潛艇尾部3.5 倍艇長,側向的四個面距離Suboff中心均2倍艇長,假定Suboff航速為3.5 kn。

    圖3 計算域邊界條件設定Fig.3 Boundary condition setting in calculation domain

    對于整個計算域采用多面體網格、棱柱層網格以及表面重構這三者相結合的方法進行網格的生成。Lx為水下航行器與Suboff 圍殼后緣之間的水平距離,h為水下航行器與Suboff 甲板之間的垂向距離,α表示速度攻角,Lw表示平行中體長度,選取對接縱向位置Lx=Lw。分別對水下航行器相距母平臺不同距離、不同攻角下的水動力特性進行了預報(典型流動圖像見圖4)。

    圖4 水下航行器對接CFD計算流場圖Fig.4 CFD calculation flow field diagram of underwater vehicle docking

    水下航行器的近距離對接過程在母平臺的垂直面上方完成,干擾水動力系數和對接距離的關系如圖5~8所示,通過分段函數擬合可得到對接干擾力模型。其中,X’、Z’和M’分別為水下航行器縱向力、垂向力和俯仰力矩系數。表1是對接距離為3.5L和無界下的縱向水動力系數。由表1和圖6~8可得到結論:在h≥3.5L時,水下航行器的水動力系數和無平臺干擾下的水動力系數幾乎一致,水動力幾乎沒有干擾影響,可視為無界環(huán)境。

    表1 距離h=3.5L和無界下的X’水動力系數(×10-3)Tab.1 Hydrodynamic coefficients of X’at a distance of h=3.5L and unbounded(×10-3)

    圖5 無界下的水動力系數(沒有干擾下)Fig.5 Hydrodynamic coefficients under unbounded conditions

    圖6 干擾水動力系數dX’和對接距離的關系曲線Fig.6 Relationship between interference hydrodynamic coefficient dX’and docking distance

    圖7 干擾水動力系數dZ’和對接距離的關系曲線Fig.7 Relationship between interference hydrodynamic coefficient dZ’and docking distance

    圖8 干擾水動力系數dM’和對接距離的關系曲線Fig.8 Relationship between interference hydrodynamic coefficient dM’and docking distance

    計算結果表明:(1)在h≥3.5L時,可視為無界環(huán)境,在無界環(huán)境下,隨著攻角增大,縱向力幾乎無變化,垂向力和攻角呈負相關關系,縱向力矩和攻角呈正相關關系;(2)在對接過程中,水下航行器的縱向力隨著對接距離的減小呈現逐漸增加的趨勢,攻角在α≤4°范圍內,對接距離對縱向力系數影響不大,且在0°攻角時達到最小值,也就是說,在0°攻角時,水下航行器的阻力最小;(3)在3.5 kn航速下的對接過程中,垂向力隨著攻角的增大呈現逐漸減小的趨勢,在0°攻角時垂向力最??;垂向干擾水動力系數隨著對接距離的減小先增大后減小,呈現非線性干擾特性;(4)在3.5 kn航速下的對接過程中,俯仰力矩隨著攻角的增大呈現逐漸增加的趨勢,在0°攻角時俯仰力矩最?。桓┭龈蓴_水動力系數隨著距離的減小先增大后減小,呈現非線性干擾特性。

    2.2 舵/推進器模型

    無人水下航行器所受到的控制力和力矩主要由舵機和推進器產生,舵包括艏升降舵、艉升降舵和艉方向舵,舵力表達式如式(6)所示。

    式中,δb為艏升降舵角,δs為艉升降舵舵角,δr為艉方向舵舵角,其他項無量綱水動力系數可參考文獻[8]。

    在與母艇航行對接過程中,通過主推保持水下航行器航行狀態(tài),艏舵控制航行深度變化,艉舵控制航行縱傾變化;在母艇懸停靜態(tài)對接過程中,通過垂向推進器定深和縱傾控制完成水下對接過程。

    2.3 對接動力學模型

    水下航行器對接過程的運動方程如下:

    式中:水下航行器η1=[X1,Y1,Z1,φ1,θ1,ψ1]T,v1=[u1,v1,w1,p1,q1,r1]T;母平臺η2=[X2,Y2,Z2,φ2,θ2,ψ2]T,v2=[u2,v2,w2,p2,q2,r2]T;M1為航行器的慣性矩陣,C1為航行器的科氏力和向心力矩陣,D1為航行器的水動力矩陣,g1為航行器的靜力矩陣,τT1為航行器的推進器推力,τδ1為航行器的舵力,F(η1,η2)為航行器與母平臺間的干擾水動力。

    采用美國海軍艦船研究與發(fā)展中心(NSRDC)報告中的標準方程,以水下航行器六自由度運動數學模型為基礎,考慮航行器和平臺間的動力學干擾,得到水下航行器對接過程的動力學模型方程組[9-10]。水下航行器與母平臺的對接假設在母平臺中縱垂面上進行,縱向、垂向和縱傾方程如下:

    縱向方程為

    垂向方程式

    縱傾方程為

    式中:Fi(x,y,z)為對接干擾水動力,與兩者距離h有關,例如式Fi(x,y,z)中,擬合項干擾dZ'ww如式(11)所示;水下航行器其它參數可參考文獻[11-12]。

    3 水下航行器對接仿真研究

    3.1 水下航行器對接干擾仿真研究

    (1)為研究不同對接距離下干擾水動力對姿態(tài)的影響,選取母平臺以3.5 kn 航速進行對接,此時水下航行器完成縱向相對位置的跟蹤,其航速保持在3.5 kn,當初始對接距離h0=0.25L和h0=0.4L時對應的垂向干擾水動力和縱傾水動力矩較大,對垂向速度和縱傾角影響如圖9和圖10所示。

    圖9 不同對接距離下的垂向速度變化Fig.9 Vertical velocity change under different docking distances

    圖10 不同對接距離下的縱傾角變化Fig.10 Pitch change under different docking distances

    仿真結果表明,h0≥3.5L時,水動力干擾影響可忽略,由于上下不對稱,導致產生微小垂向速度和縱傾;h0=0.25L時,干擾水動力對垂向速度和縱傾角影響最大,垂向速度達到0.02 m/s,縱傾值將近4°,尤其在航速3.5 kn時對動態(tài)對接帶來了難度,垂向速度和縱傾變化和干擾水動力相關,在實際對接時需要引起重視。

    (2)為研究不同對接航速下干擾水動力對姿態(tài)影響,選取初始對接距離為h0=0.25L,水下航行器分別以2.5 kn、3 kn、3.5 kn 航速進行對接時,垂向速度和縱傾角隨時間的變化規(guī)律如圖11 和圖12所示。

    圖11 不同航速下的垂向速度變化Fig.11 Vertical speed change at different speeds

    圖12 不同航速下的縱傾角變化Fig.12 Pitch change at different speeds

    從仿真結果看出,隨著航速的增大,干擾水動力對縱傾和垂向速度的影響增大,其中,在3.5 kn航速下,干擾產生的縱傾角小于4°,該縱傾角并不大,反映該水下航行器外形設計和水動力性能優(yōu)良。在對接中,可進一步通過對接控制實現精準對接。

    3.2 水下航行器對接控制仿真研究

    對接時,水下航行器航行至對接母平臺上方,首先通過縱向位置控制器,完成水下航行器與母平臺縱向相對位置的控制,保證水下航行器和母平臺的對接縱向相對位置偏差在0.2 m 范圍內,然后通過定深和縱傾調節(jié)完成與母平臺的對接。選擇航速3.5 kn時開啟自動對接,定深控制采用艏舵,縱傾控制采用艉舵,縱向位置跟蹤和縱傾控制采用PID 控制方法,深度控制采用滑??刂品椒?。對接控制方案如圖13所示。

    圖13 對接控制結構框圖Fig.13 Block diagram of docking control structure

    將深度誤差取為e(t),切換面選擇為s(t),定義滑模函數為

    在式(12)中,滿足Hurwitz條件,取c>0,并選擇指數逼近定律,如式(13)所示。

    在對接過程中,縱傾值不大,取Z¨ ≈w˙,θ˙≈q,由此,可得到深度控制器設計滑模控制律,如式(14)所示,為消除抖震影響,在控制中,采用飽和函數sat(s)代替切換函數sgn(s)。

    式中各參數表達可參考文獻[13]。其中,F為

    飽和函數表達式為

    從仿真結果(圖14~16)可以看出:深度控制方面,水下航行器在28 s時達到了期望深度,在50 s后穩(wěn)定,滑模控制對水動力干擾的影響具有一定適應性;縱傾控制方面,結合3.1 節(jié)仿真結果,考慮縱傾值不大和系統的穩(wěn)定性兩方面因素,采用PID 控制,由于干擾水動力影響,控制過程縱傾角有一定波動,但縱傾角不大,整體變化在θ≤2°范圍,最終控制精度達0.5°。

    圖14 對接控制縱向相對位置值變化Fig.14 Change in longitudinal relative position of docking control

    圖15 對接控制深度值變化Fig.15 Change of depth of docking control

    圖16 對接控制縱傾值變化Fig.16 Change of trim value of docking control

    4 結 論

    本文研究了不同距離下的雙體水動力干擾規(guī)律,建立了小型水下航行器對接過程的運動模型,研究了干擾水動力對姿態(tài)和速度的影響。在此基礎上,采用滑模控制方法,完成定深和縱傾控制,并完成了相應的控制仿真計算。研究結果表明:(1)在h≥3.5L時,可視為無界環(huán)境,在無界環(huán)境下,隨著攻角增大,縱向力幾乎無變化,垂向力和攻角呈負相關關系,縱向力矩和攻角呈正相關關系;(2)在對接過程中,水下航行器的縱向力隨著對接距離的減小呈現逐漸增加的趨勢,垂向干擾水動力系數隨著距離的減小先增大后減小,俯仰干擾水動力系數隨著距離的減小先增大后減小,呈現非線性干擾特性;(3)仿真結果表明,對接時的水動力干擾主要影響航行器的縱傾和垂向速度,對接運動模型具有很強的非線性特點,滑??刂茖λ聦舆\動控制具有很好的適用性。

    猜你喜歡
    航速攻角航行
    VLCC在波浪中的航速優(yōu)化與能效優(yōu)化分析
    提升全回轉港作拖輪航速的有效途徑
    水上消防(2022年1期)2022-06-16 08:06:56
    到慧骃國的航行
    風標式攻角傳感器在超聲速飛行運載火箭中的應用研究
    低速水面目標航速精度分析及精確解算
    大攻角狀態(tài)壓氣機分離流及葉片動力響應特性
    小舟在河上航行
    中學生英語(2017年6期)2017-07-31 21:28:55
    航行
    青年歌聲(2017年6期)2017-03-13 00:57:56
    基于CFD的波浪滑翔機航速預測
    附加攻角效應對顫振穩(wěn)定性能影響
    振動與沖擊(2015年2期)2015-05-16 05:37:34
    亚洲人成电影观看| 国产精品二区激情视频| 色婷婷久久久亚洲欧美| 久久国产亚洲av麻豆专区| 国产高清国产精品国产三级| 成熟少妇高潮喷水视频| 国产亚洲欧美98| 99久久精品国产亚洲精品| 热99re8久久精品国产| 一区在线观看完整版| 亚洲色图 男人天堂 中文字幕| 国产av又大| 国产精华一区二区三区| 两性夫妻黄色片| 午夜精品在线福利| 丝袜人妻中文字幕| 日韩人妻精品一区2区三区| 夜夜爽天天搞| 啪啪无遮挡十八禁网站| 亚洲aⅴ乱码一区二区在线播放 | 一进一出好大好爽视频| 国产成人系列免费观看| 精品一区二区三卡| av中文乱码字幕在线| 亚洲精品在线观看二区| 久久人人精品亚洲av| 国产野战对白在线观看| 美女大奶头视频| 国产又色又爽无遮挡免费看| 欧美黑人精品巨大| 国产av一区二区精品久久| 精品人妻在线不人妻| 国产精品免费一区二区三区在线| 国产一卡二卡三卡精品| 日韩欧美三级三区| 91麻豆精品激情在线观看国产 | 日韩欧美在线二视频| 两人在一起打扑克的视频| 国产成人欧美在线观看| 久久久精品国产亚洲av高清涩受| 国产亚洲欧美98| 亚洲成a人片在线一区二区| 69精品国产乱码久久久| 热99re8久久精品国产| 久久久久久亚洲精品国产蜜桃av| aaaaa片日本免费| 亚洲欧美激情在线| 久久青草综合色| 国产在线精品亚洲第一网站| 在线av久久热| 亚洲一区高清亚洲精品| 久久久久久亚洲精品国产蜜桃av| 国产成人欧美在线观看| 大码成人一级视频| 午夜久久久在线观看| 两人在一起打扑克的视频| 亚洲国产中文字幕在线视频| 级片在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲色图综合在线观看| 可以在线观看毛片的网站| 麻豆一二三区av精品| 啦啦啦在线免费观看视频4| 国产精品电影一区二区三区| 中文字幕av电影在线播放| 丝袜美足系列| 不卡av一区二区三区| 又黄又爽又免费观看的视频| 黄色成人免费大全| 精品一区二区三区视频在线观看免费 | 在线av久久热| 中文字幕色久视频| 欧美中文综合在线视频| 精品久久久久久成人av| 黄色丝袜av网址大全| 精品第一国产精品| 可以在线观看毛片的网站| 一区二区三区国产精品乱码| 亚洲自偷自拍图片 自拍| 757午夜福利合集在线观看| 国产乱人伦免费视频| 国产成人欧美| 日韩免费av在线播放| 12—13女人毛片做爰片一| 亚洲伊人色综图| 久热这里只有精品99| 亚洲激情在线av| 久久天堂一区二区三区四区| 午夜福利,免费看| 成人免费观看视频高清| 亚洲熟妇中文字幕五十中出 | 中文字幕av电影在线播放| 看片在线看免费视频| 高清欧美精品videossex| 精品福利观看| 中文亚洲av片在线观看爽| 亚洲成人久久性| 日韩精品青青久久久久久| 婷婷精品国产亚洲av在线| 午夜a级毛片| 国产免费av片在线观看野外av| 亚洲精品一二三| 欧美日本中文国产一区发布| 美女高潮到喷水免费观看| 亚洲精品国产精品久久久不卡| 国产一区二区三区综合在线观看| 国产精品自产拍在线观看55亚洲| 极品教师在线免费播放| 1024视频免费在线观看| 如日韩欧美国产精品一区二区三区| 男人的好看免费观看在线视频 | 91国产中文字幕| 99热国产这里只有精品6| 国产成人欧美在线观看| 国产精品久久电影中文字幕| 精品第一国产精品| 国产亚洲精品综合一区在线观看 | 又黄又爽又免费观看的视频| 热re99久久精品国产66热6| 男女午夜视频在线观看| 欧美黑人精品巨大| 亚洲,欧美精品.| 成人亚洲精品一区在线观看| 日本免费a在线| 久久99一区二区三区| 一区二区三区精品91| 三级毛片av免费| 法律面前人人平等表现在哪些方面| 精品第一国产精品| 91九色精品人成在线观看| 欧美精品亚洲一区二区| 身体一侧抽搐| 国产野战对白在线观看| 另类亚洲欧美激情| 久久久国产精品麻豆| 亚洲精品中文字幕一二三四区| 一区二区三区精品91| 夜夜夜夜夜久久久久| 法律面前人人平等表现在哪些方面| 一级毛片女人18水好多| 热99国产精品久久久久久7| 国产高清视频在线播放一区| 欧美av亚洲av综合av国产av| 亚洲成av片中文字幕在线观看| 亚洲精品中文字幕一二三四区| 免费在线观看黄色视频的| 91麻豆av在线| 欧美在线一区亚洲| av在线天堂中文字幕 | 国产欧美日韩综合在线一区二区| 99久久综合精品五月天人人| 日韩大尺度精品在线看网址 | 精品高清国产在线一区| 久久久久国产精品人妻aⅴ院| 黄网站色视频无遮挡免费观看| 欧美日本亚洲视频在线播放| 亚洲,欧美精品.| 欧美日韩乱码在线| 后天国语完整版免费观看| 欧美一级毛片孕妇| 精品国内亚洲2022精品成人| 9热在线视频观看99| 好男人电影高清在线观看| 在线永久观看黄色视频| 99国产综合亚洲精品| 亚洲五月婷婷丁香| 国产精品亚洲av一区麻豆| 亚洲精品久久午夜乱码| 老熟妇乱子伦视频在线观看| ponron亚洲| 日韩欧美三级三区| 99国产精品一区二区蜜桃av| 亚洲色图 男人天堂 中文字幕| 婷婷丁香在线五月| 亚洲午夜理论影院| 九色亚洲精品在线播放| 男人舔女人的私密视频| av视频免费观看在线观看| 国产精品自产拍在线观看55亚洲| 国产精品久久视频播放| 97超级碰碰碰精品色视频在线观看| 制服人妻中文乱码| 老司机深夜福利视频在线观看| 国产深夜福利视频在线观看| 亚洲激情在线av| 狠狠狠狠99中文字幕| 18禁黄网站禁片午夜丰满| 黑人猛操日本美女一级片| 黑人操中国人逼视频| 999久久久精品免费观看国产| 男人舔女人下体高潮全视频| 国产av精品麻豆| 欧美日韩一级在线毛片| 精品第一国产精品| 欧美日韩瑟瑟在线播放| 波多野结衣高清无吗| 国产一区二区在线av高清观看| 青草久久国产| √禁漫天堂资源中文www| 欧美国产精品va在线观看不卡| ponron亚洲| 亚洲视频免费观看视频| 亚洲aⅴ乱码一区二区在线播放 | 可以免费在线观看a视频的电影网站| 我的亚洲天堂| 亚洲 欧美 日韩 在线 免费| 在线观看免费视频网站a站| 国产亚洲精品久久久久5区| 黄色片一级片一级黄色片| 午夜免费鲁丝| 国产精品亚洲一级av第二区| av欧美777| 黄频高清免费视频| 亚洲人成伊人成综合网2020| 久久青草综合色| ponron亚洲| 亚洲熟妇中文字幕五十中出 | 激情在线观看视频在线高清| 欧美+亚洲+日韩+国产| 国产一区二区三区综合在线观看| 69av精品久久久久久| 无人区码免费观看不卡| 国产成人av激情在线播放| 变态另类成人亚洲欧美熟女 | 一区二区日韩欧美中文字幕| 男女之事视频高清在线观看| 国产极品粉嫩免费观看在线| 国产在线精品亚洲第一网站| 九色亚洲精品在线播放| 国产精品久久久av美女十八| 琪琪午夜伦伦电影理论片6080| tocl精华| 最新美女视频免费是黄的| 老司机深夜福利视频在线观看| 神马国产精品三级电影在线观看 | 啪啪无遮挡十八禁网站| 成年女人毛片免费观看观看9| 国产熟女xx| 老熟妇仑乱视频hdxx| 国产欧美日韩精品亚洲av| 午夜免费激情av| 亚洲avbb在线观看| 国产成+人综合+亚洲专区| 身体一侧抽搐| 精品乱码久久久久久99久播| 一夜夜www| 首页视频小说图片口味搜索| 日日爽夜夜爽网站| 精品一区二区三区视频在线观看免费 | 男女午夜视频在线观看| 大陆偷拍与自拍| 日本wwww免费看| 成人精品一区二区免费| 久久久久久人人人人人| 91精品三级在线观看| 香蕉国产在线看| 国产日韩一区二区三区精品不卡| 美女国产高潮福利片在线看| 涩涩av久久男人的天堂| 国产精品国产av在线观看| 夜夜爽天天搞| 成人18禁高潮啪啪吃奶动态图| 欧美成人午夜精品| 免费高清视频大片| 在线观看免费午夜福利视频| 午夜91福利影院| 亚洲国产精品sss在线观看 | 国产av一区在线观看免费| 亚洲成人免费电影在线观看| 国产一区二区三区综合在线观看| 18禁黄网站禁片午夜丰满| 黑人操中国人逼视频| 在线观看一区二区三区激情| 国产精品 国内视频| 天堂影院成人在线观看| 国产精品综合久久久久久久免费 | 欧美黑人欧美精品刺激| av天堂久久9| 欧美在线黄色| 亚洲欧美激情综合另类| 欧美一区二区精品小视频在线| 日韩高清综合在线| 亚洲精品一卡2卡三卡4卡5卡| 亚洲九九香蕉| 亚洲国产欧美网| 91精品三级在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩无卡精品| 最新在线观看一区二区三区| 好男人电影高清在线观看| 国产色视频综合| av免费在线观看网站| 日韩欧美免费精品| 欧美日韩瑟瑟在线播放| 久久人妻福利社区极品人妻图片| 亚洲午夜理论影院| 国产精品爽爽va在线观看网站 | 亚洲一卡2卡3卡4卡5卡精品中文| 精品一区二区三区视频在线观看免费 | 91麻豆av在线| 在线观看一区二区三区激情| 狂野欧美激情性xxxx| 嫩草影视91久久| 桃红色精品国产亚洲av| 亚洲熟女毛片儿| 一个人观看的视频www高清免费观看 | 看片在线看免费视频| 88av欧美| 人人妻人人添人人爽欧美一区卜| 少妇被粗大的猛进出69影院| 夫妻午夜视频| 80岁老熟妇乱子伦牲交| www.熟女人妻精品国产| 一区福利在线观看| 美女扒开内裤让男人捅视频| 无遮挡黄片免费观看| 美女 人体艺术 gogo| 亚洲精品国产区一区二| 午夜免费激情av| 国产精品98久久久久久宅男小说| 久久久久久久久久久久大奶| 9热在线视频观看99| 51午夜福利影视在线观看| 在线观看舔阴道视频| 女性生殖器流出的白浆| 国产1区2区3区精品| 亚洲熟妇中文字幕五十中出 | 成在线人永久免费视频| 国产亚洲精品综合一区在线观看 | 美女国产高潮福利片在线看| 久久久久久久精品吃奶| 亚洲av五月六月丁香网| 电影成人av| 久久精品亚洲熟妇少妇任你| а√天堂www在线а√下载| 欧美日韩国产mv在线观看视频| 久久久久久大精品| 黑人操中国人逼视频| 亚洲av日韩精品久久久久久密| 亚洲精品久久午夜乱码| 国产午夜精品久久久久久| 人人妻人人添人人爽欧美一区卜| 国产成年人精品一区二区 | 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品国产一区二区精华液| 亚洲欧美日韩无卡精品| 99久久综合精品五月天人人| 男人舔女人下体高潮全视频| 黑人猛操日本美女一级片| 日韩精品免费视频一区二区三区| 夫妻午夜视频| 中文字幕人妻熟女乱码| 欧美一区二区精品小视频在线| 久久久国产成人免费| 亚洲一区二区三区色噜噜 | 午夜激情av网站| 欧美在线一区亚洲| 欧美 亚洲 国产 日韩一| 亚洲av美国av| 首页视频小说图片口味搜索| 变态另类成人亚洲欧美熟女 | 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产亚洲av香蕉五月| av网站在线播放免费| 美女午夜性视频免费| 欧美成人免费av一区二区三区| 无遮挡黄片免费观看| 国产精品亚洲av一区麻豆| 久久欧美精品欧美久久欧美| 国产成人影院久久av| 国产免费男女视频| 久久99一区二区三区| 欧美日韩瑟瑟在线播放| 最近最新免费中文字幕在线| 国产视频一区二区在线看| 久久精品91无色码中文字幕| 天堂俺去俺来也www色官网| 久久久久九九精品影院| 亚洲精品国产区一区二| 国产国语露脸激情在线看| 亚洲全国av大片| 操出白浆在线播放| 亚洲精品一区av在线观看| 99久久精品国产亚洲精品| 新久久久久国产一级毛片| 最新在线观看一区二区三区| 午夜福利欧美成人| 在线播放国产精品三级| 国产成人av教育| 老司机靠b影院| 亚洲人成电影免费在线| 午夜免费成人在线视频| 亚洲第一欧美日韩一区二区三区| 一区二区日韩欧美中文字幕| 两个人看的免费小视频| 日本黄色日本黄色录像| 99久久99久久久精品蜜桃| 叶爱在线成人免费视频播放| 日韩精品中文字幕看吧| 亚洲一区二区三区欧美精品| 久久精品91蜜桃| 大香蕉久久成人网| ponron亚洲| 精品一区二区三区四区五区乱码| 国产在线观看jvid| 日韩欧美在线二视频| 亚洲视频免费观看视频| 久久青草综合色| 琪琪午夜伦伦电影理论片6080| 亚洲av片天天在线观看| 51午夜福利影视在线观看| 老熟妇乱子伦视频在线观看| 黄片大片在线免费观看| 女人被躁到高潮嗷嗷叫费观| 免费不卡黄色视频| 精品国产一区二区三区四区第35| 99精品在免费线老司机午夜| 91成年电影在线观看| 91在线观看av| 日本a在线网址| 激情在线观看视频在线高清| 大香蕉久久成人网| 亚洲全国av大片| 一进一出抽搐动态| 韩国精品一区二区三区| www国产在线视频色| 欧美在线黄色| 最近最新中文字幕大全电影3 | 麻豆久久精品国产亚洲av | 99re在线观看精品视频| 国产精品一区二区在线不卡| 黄色视频不卡| av电影中文网址| 69av精品久久久久久| 国产亚洲av高清不卡| 香蕉国产在线看| 大型黄色视频在线免费观看| 亚洲成人国产一区在线观看| 久久青草综合色| 久久久国产精品麻豆| 成人永久免费在线观看视频| 母亲3免费完整高清在线观看| 亚洲成人免费电影在线观看| 女警被强在线播放| 久久久国产欧美日韩av| 国产精品永久免费网站| 亚洲熟妇熟女久久| 午夜福利影视在线免费观看| 色精品久久人妻99蜜桃| 成人av一区二区三区在线看| 成人国语在线视频| 丰满的人妻完整版| 中文字幕av电影在线播放| 成在线人永久免费视频| 午夜激情av网站| 欧美成狂野欧美在线观看| 亚洲专区中文字幕在线| 亚洲第一欧美日韩一区二区三区| 国产一卡二卡三卡精品| 成人精品一区二区免费| 高清黄色对白视频在线免费看| 女同久久另类99精品国产91| bbb黄色大片| 丁香欧美五月| 制服人妻中文乱码| 90打野战视频偷拍视频| 欧美激情久久久久久爽电影 | 国产在线精品亚洲第一网站| 免费高清在线观看日韩| 欧美不卡视频在线免费观看 | 97碰自拍视频| 校园春色视频在线观看| 欧美黄色片欧美黄色片| 欧美日韩国产mv在线观看视频| 日韩av在线大香蕉| 麻豆久久精品国产亚洲av | 丰满的人妻完整版| 一级毛片精品| 精品免费久久久久久久清纯| 一区二区三区激情视频| 免费不卡黄色视频| 日韩欧美一区二区三区在线观看| 国产精品野战在线观看 | 成人特级黄色片久久久久久久| 人人妻人人澡人人看| 丁香欧美五月| 人人澡人人妻人| 亚洲av电影在线进入| 手机成人av网站| tocl精华| 亚洲国产精品999在线| 级片在线观看| 免费高清在线观看日韩| 久热爱精品视频在线9| 欧美日韩一级在线毛片| 国产又色又爽无遮挡免费看| 一级毛片女人18水好多| 亚洲av第一区精品v没综合| 亚洲中文字幕日韩| e午夜精品久久久久久久| 少妇被粗大的猛进出69影院| 美女福利国产在线| 老鸭窝网址在线观看| 亚洲一区二区三区色噜噜 | 天天躁夜夜躁狠狠躁躁| 久久久国产一区二区| 亚洲熟妇中文字幕五十中出 | 99riav亚洲国产免费| 大型黄色视频在线免费观看| 欧美中文综合在线视频| 五月开心婷婷网| 国产精品国产高清国产av| 亚洲色图 男人天堂 中文字幕| 亚洲美女黄片视频| 精品第一国产精品| 香蕉丝袜av| 亚洲人成网站在线播放欧美日韩| 色婷婷av一区二区三区视频| 国产一卡二卡三卡精品| 国产熟女xx| 日韩大码丰满熟妇| 午夜免费成人在线视频| 搡老乐熟女国产| 亚洲av电影在线进入| 一级a爱视频在线免费观看| 国产欧美日韩综合在线一区二区| 在线观看免费视频网站a站| 国产激情欧美一区二区| 人成视频在线观看免费观看| 999久久久国产精品视频| 熟女少妇亚洲综合色aaa.| 国产黄a三级三级三级人| 色婷婷久久久亚洲欧美| 99精品久久久久人妻精品| 午夜福利欧美成人| 一级毛片精品| 老鸭窝网址在线观看| 亚洲 国产 在线| 欧美中文综合在线视频| 91成年电影在线观看| 一二三四在线观看免费中文在| 精品国产一区二区久久| 欧美午夜高清在线| 在线观看舔阴道视频| 嫁个100分男人电影在线观看| 国产三级在线视频| xxx96com| 天堂中文最新版在线下载| 麻豆久久精品国产亚洲av | 国产精品日韩av在线免费观看 | 男人舔女人下体高潮全视频| 国产精品野战在线观看 | 纯流量卡能插随身wifi吗| 熟女少妇亚洲综合色aaa.| 亚洲中文av在线| 一二三四在线观看免费中文在| 很黄的视频免费| 18禁国产床啪视频网站| 高清欧美精品videossex| 波多野结衣一区麻豆| 国产精品香港三级国产av潘金莲| 中文字幕av电影在线播放| 亚洲五月天丁香| 人人妻,人人澡人人爽秒播| 成人18禁高潮啪啪吃奶动态图| 91av网站免费观看| 国产成人欧美| 欧美中文综合在线视频| 久久久精品欧美日韩精品| 日韩中文字幕欧美一区二区| 欧美中文日本在线观看视频| 亚洲久久久国产精品| 欧美日韩瑟瑟在线播放| 国产成人av教育| 午夜免费观看网址| 国产精品久久久久成人av| 如日韩欧美国产精品一区二区三区| 精品一区二区三卡| 精品高清国产在线一区| 午夜免费激情av| 丁香欧美五月| 久久国产亚洲av麻豆专区| 这个男人来自地球电影免费观看| 日韩成人在线观看一区二区三区| xxxhd国产人妻xxx| 母亲3免费完整高清在线观看| av天堂在线播放| 夜夜躁狠狠躁天天躁| 免费搜索国产男女视频| 精品人妻1区二区| 高潮久久久久久久久久久不卡| 免费人成视频x8x8入口观看| 亚洲精品在线观看二区| 人人澡人人妻人| 精品福利永久在线观看| 国产深夜福利视频在线观看| 黑人欧美特级aaaaaa片| 日日干狠狠操夜夜爽| www.自偷自拍.com| 两个人看的免费小视频| 亚洲aⅴ乱码一区二区在线播放 | 欧洲精品卡2卡3卡4卡5卡区| 极品教师在线免费播放| 悠悠久久av| 亚洲精品国产色婷婷电影| 精品免费久久久久久久清纯| 91字幕亚洲| 久久人妻av系列| 久久久久久久精品吃奶| av天堂在线播放| 怎么达到女性高潮| 久久天躁狠狠躁夜夜2o2o| 亚洲自拍偷在线| 麻豆一二三区av精品| 久久精品亚洲av国产电影网|