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

    基于PBD 算法的海上補給裝備柔性體仿真技術(shù)

    2023-04-02 10:56:36李楠任鴻翔
    兵工學(xué)報 2023年3期
    關(guān)鍵詞:索具附著物視景

    李楠, 任鴻翔

    (1. 海軍勤務(wù)學(xué)院,天津 300450; 2. 大連海事大學(xué),遼寧 大連 116085)

    0 引言

    海上補給是提升海軍艦艇的續(xù)航力和作戰(zhàn)半徑的有效手段[1]。海上補給裝備系統(tǒng)復(fù)雜,技術(shù)要求高,訓(xùn)練難度大。開展仿真訓(xùn)練成本較低、效率較高。目前開發(fā)的海上補給仿真訓(xùn)練系統(tǒng)在訓(xùn)練架設(shè)索道的環(huán)節(jié)中,一般采用兩端掛點實時距離計算并進行直線連點的方式[2],與補給實際過程嚴重不符,無法模擬由于橫風(fēng)載荷和船體搖擺等帶來的索具振動,導(dǎo)致官兵在訓(xùn)練時,無法開展索具振動控制技術(shù)訓(xùn)練。

    索具等柔性物體建模的研究始于20 世紀初,在當(dāng)時的計算條件和數(shù)值仿真技術(shù)條件下,對索具的研究僅限于穩(wěn)態(tài)特性。隨著計算機計算性能的大幅提升,對索具的研究也有了深入提高。目前索具類柔性體主要使用的動力學(xué)模型包括:靜力學(xué)平面非線性模型、質(zhì)量彈簧系統(tǒng)模型、有限元模型和多體動力學(xué)模型等[3]。

    基于位置動力學(xué)(PBD)算法具有穩(wěn)定性強、允許較大的時間步長以及較好的收斂性等優(yōu)點。本文從多體動力學(xué)模型出發(fā),建立補給裝備索具柔性體模型,并運用PBD 算法開展海上補給裝備索具視景仿真研究,旨在提高海上補給仿真訓(xùn)練系統(tǒng)的真實感和訓(xùn)練效能。

    1 海上補給裝備柔性體仿真建模

    1.1 索具多體動力學(xué)模型

    索具多體動力學(xué)模型是將繩索離散成若干索具段,每個索具段視為剛體,各索具段之間由連接件連接便可建立一個多剛體系統(tǒng)[4]。

    索具多剛體動力學(xué)模型的基礎(chǔ)上,在索具單元之間的連接件上添加適當(dāng)?shù)淖枘岷蛷椥孕巫儯缙揭谱枘?、平移剛度、旋轉(zhuǎn)阻尼、旋轉(zhuǎn)剛度陣等幾種矩陣,使剛體單元之間的聯(lián)接變成柔性連接,從而將索具的多剛體動力學(xué)模型轉(zhuǎn)化為索具剛?cè)峄旌蟿恿W(xué)模型[5]。

    在此基礎(chǔ)上,進一步將各索具單元進行柔性處理成為具備6 階模態(tài)向量的柔性簡支梁,剛?cè)峄旌象w系統(tǒng)就轉(zhuǎn)換成了多柔體系統(tǒng)[6-7]。

    海上補給裝備柔性體包括高架索、牽引索、輸油軟管等,動力學(xué)特性均與繩索相近。因此,本文采用索具多體動力學(xué)模型建模方法進行海上補給裝備柔性體視景仿真建模。

    1.2 海上補給裝備柔性體視景仿真建模

    海上補給裝備中柔性體包括輸油軟管、高架索、牽索、引纜等。其運動形式屬于三維動態(tài)運動,為準確反映其運動,在建模時將索具視為剛?cè)峄旌象w系統(tǒng)。如圖2 所示,將索具在空間上離散為N段單元索具段,索具段之間通過無質(zhì)量的質(zhì)點連接,將所有外載荷(阻力、重力、浮力等)的作用集中到索具段上,并計入慣性質(zhì)量的影響,建立索具段點的運動方程組,通過相應(yīng)的數(shù)值計算方法求[8-10]。

    應(yīng)用牛頓第二定律,得到索具節(jié)點的基本運動控制方程:

    式中:Mi為索具節(jié)點的質(zhì)量矩陣,包括索具節(jié)點自身的慣性質(zhì)量mi和其在流體中的附加質(zhì)量Mai;為第i個節(jié)點的加速度;Fi為作用于節(jié)點i上的外力,包括索具節(jié)點的張力ΔTi、內(nèi)部阻尼力ΔCi、浮力Bi和重力Gi等。式(1)中的mi、可以表示為

    式中:Pi為3×3 階單位矩陣,ωi為繩索單元轉(zhuǎn)動慣量;ln為節(jié)點間未拉伸長度;ρ為繩索材料密度;Ca為附加質(zhì)量系數(shù);si為為微元段的橫截面積;Pi為微元段在浮動坐標系上的投影,下標i-1/2和i+1/2 分別表示節(jié)點i-1 到i、節(jié)點i到i+1 之間的物理量,如li+1/2=si+1-si表示節(jié)點i到節(jié)點i+1 間的長度;θ、φ分別為浮動坐標系與地球坐標系之間的夾角。

    索具動力學(xué)模型解算如圖1 所示。

    圖1 柔性索具動力學(xué)模型解算流程Fig. 1 Computation flow of the dynamic model of the flexible rigging

    2 海上補給裝備柔性體視景仿真

    利用多體動力學(xué)建立的海上補給裝備索具等柔性體模型,通過計算,可以很好地反映索具的拉伸、彎曲、扭轉(zhuǎn)等變形效應(yīng),但在三維可視化仿真中,實時解算仍存在迭代次數(shù)多、收斂難、失效性較差等問題,而在使用PBD 方法模擬索具,需要多次迭代才能使索具的距離約束條件收斂,容易引起模擬的索具產(chǎn)生彈性失真。本文在普通索具模擬的約束條件中增加長距離約束、引腳約束等約束,從而解決了海上補給裝備柔性體在視景仿真中失真的問題。

    2.1 PBD 算法概述

    在PBD 方法中,索具用約束耦合的粒子表示,在每個時間步長中迭代地求解這些約束條件[11]。對索具粒子的模擬迭代計算分為3 個步驟:1)根據(jù)速度和外力移動索具粒子;2)使索具粒子滿足約束條件;3)基于給定數(shù)據(jù)和時間步長Δt,執(zhí)行時間 積分。

    索具粒子由N粒子和M個約束組成。粒子i∈[1, ···,N]的質(zhì)量為mi、位置為xi和速度為xi,約束j∈[1, ···,M]有如下5 種特性:

    1)約束基數(shù)nj,即第j個約束影響的粒子數(shù)目為nj。

    3){i1,… ,inj},ik∈ {1,… ,N}。

    4)約束對應(yīng)的剛度參數(shù)kj∈[ 0,1]。

    5)約束類型分為單邊約束Cj(xi1, …,xinj)≥0與雙邊約束。

    在模擬步驟中初始化索具粒子的位置,速度和質(zhì)量倒數(shù)。把所有不能轉(zhuǎn)換為位置約束的力(如重力)記為fext(xi),并根據(jù)fext(xi)的值進行一次數(shù)值積分預(yù)測速度vi。然后直接對速度進行衰減操作模擬阻尼。計算完速度后再計算位置的預(yù)測值pi。在迭代計算的每個時間步長中生成附加的碰撞約束Mcoll,碰撞頂點的數(shù)量決定了碰撞約束Mcoll的數(shù)量。得到M個內(nèi)部約束和Mcoll個外部約束之后,對這些約束進行迭代求解,得到校正后的粒子位置。最后根據(jù)摩擦和恢復(fù)系數(shù)修改速度vi。

    2.2 海上補給裝備柔性體的約束設(shè)置與計算

    海上補給裝備柔性體在橫風(fēng)、船體搖擺、懸掛重物移動等外部載荷作用下,會出現(xiàn)拉升、彎曲、附著、碰撞等變形。能否在視景中逼真模擬以上各種柔性體變形是衡量海上補給仿真訓(xùn)練系統(tǒng)逼真度高低的重要指標[12]。

    PBD 方法通過修改索具粒子的位置關(guān)系,可以模擬出索具等柔性體拉升、彎曲、附著、碰撞的變形。在模擬柔性體拉升、彎曲時,本文使用了拉伸約束、彎曲約束,避免了索具的過度拉伸和彎曲失真;在模擬柔性體和附著物連接時,使用了長距離附著物約束和單邊引腳約束,避免了附著物的拉伸變形失真,附著物與柔性體連接點的兩個粒子合并為一個粒子進行計算[13-14]。

    2.2.1 海上補給裝備柔性體拉伸約束的設(shè)置與計算

    海上補給裝備中的高架索是實施航行橫向補給時作為承載質(zhì)量的索具。在補給過程中,始終呈現(xiàn)緊張狀態(tài)。在其上要懸掛輸油軟管鞍座或干貨行車,并承受橫風(fēng)、船搖等載荷,因此,在仿真過程中,需要實時計算其張力,并保證其柔性拉伸狀態(tài)隨載荷變化而產(chǎn)生拉伸變形。本文使用拉伸約束,考慮高架索為鋼索,抗拉伸能力強,因此設(shè)置拉伸約束系數(shù)為1,仿真效果如圖2 所示。其中圖2(a)顯示的是高架索的組成粒子及受力情況(綠色矢量為高架索垂直方向上的力,黃色矢量為高架索水平方向上的力),圖2(b)為高架索組成粒子的蒙皮截圖,圖2(c)為高架索的仿真效果。

    圖2 海上補給裝備高架索拉升效應(yīng)視景仿真效果Fig. 2 Visual simulation results of the hoisting effect of the viaduct of the replenishment equipment

    2.2.2 海上補給裝備柔性體彎曲約束的設(shè)置與計算

    除了拉伸效應(yīng)外,在海上補給裝備柔性體仿真中也要考慮彎曲效應(yīng)。柔性體模型上3 個相鄰粒子位置可以決定柔性體的彎曲效果。p1、p2、p3為柔性體3 個相鄰粒子,以p1、p3連線為基準線,通過計算p2距基準線的距離來確定柔性體的彎矩程度。由于粒子數(shù)量和解算率的不同,可能會出現(xiàn)柔性體彎曲效果失真。為此需要增加彎曲約束來保證仿真效果。

    輸油(水)軟管是海上補給裝備,在橫向補給過程中用于輸送油(水)液貨的柔性體設(shè)備。其具有高抗彎曲度但低拉伸剛度的物理特性,因此在添加拉伸約束和彎曲約束時,取拉伸約束值較低,約為0.5,取彎曲約束值較大,約為0.8,仿真效果如 圖3所示。其中圖3(a)顯示的是輸油軟管的組成粒子及受力情況(綠色矢量為高架索垂直方向上的力,黃色矢量為高架索水平方向上的力),圖3(b)為輸油軟管組成粒子的蒙皮截圖,圖3(c)為輸油軟管的仿真效果。

    圖3 海上補給裝備輸油軟管彎曲效應(yīng)視景仿真效果Fig. 3 Visual simulation results of the bending effect of the oil hose of the replenishment equipment

    2.2.3 海上補給裝備柔性體長距離附著約束的設(shè)置與計算

    海上補給裝備柔性體通常存在自由端,且自由端會附著輸油軟管、貨物、掛鉤等設(shè)備,此時柔性體會隨著附著物的運動而出現(xiàn)擺動。在使用PBD方法對上述存在附著物時柔性體的模擬中發(fā)現(xiàn),若將附著物看作柔性體的末端節(jié)點進行迭代計算時,則容易產(chǎn)生附著物被一同拉升而出現(xiàn)彈性失真。長距離附著約束只計算從附著物與柔性體接觸點粒子到柔性體自由粒子的整個路徑的拉伸,不對附著物產(chǎn)生影響[15]。

    海上干貨補給是在高架索上設(shè)置行車,行車下部通過4 根吊帶懸掛貨物小車。4 根吊帶屬于柔性體,且作業(yè)過程中呈拉伸狀態(tài),而下部懸吊的貨物小車屬于剛體,在行車運行時,貨物小車會出現(xiàn)擺動。若將貨物小車與吊帶看作一個整體進行柔性體仿真,則施加拉伸約束后,貨物小車也會出現(xiàn)拉伸變形。因此,需要在貨物小車與吊帶之間添加長距離附著約束,從而實現(xiàn)具有附著物的吊帶自由端在貨物小車重力作用下自由擺動時不出現(xiàn)貨物小車的變形失真。仿真效果如圖4 所示。

    圖4 吊掛索具與吊筐之間的擺動效應(yīng)的視景 仿真效果Fig. 4 Visual simulation results of the swing effect between the hanging rigging and the hanging basket

    2.2.4 海上補給裝備柔性體與附著物相互約束的設(shè)置與計算

    海上補給裝備柔性體端部一般都會連接附著物,如發(fā)射引纜。發(fā)射引纜是海上航行橫向補給作業(yè)的第一步,目的是便于接收艦通過引纜接收高架索等索具。引纜是一根纖維軟繩,前段掛有引纜球,通過撇纜槍發(fā)射到接收艦。發(fā)射引纜時,引纜隨撇纜球到達接收艦的過程中,柔性體和附著物都可以進行3 自由度運動,且相互影響。為模擬柔性體與附著物相互作用,本文通過銷子、銷接頭等引腳約束連接多個物體,引腳約束可以確保這些物體在每一個時間步長更新的位置與所連接物體的附著點位置一致。使用單邊引腳約束,以確保物體頂點位于動態(tài)物體的外部,從而避免了物體之間穿透現(xiàn)象的產(chǎn)生。仿真效果如圖5 所示。

    圖5 引纜與撇纜球之間的相互作用效應(yīng)的視景 仿真效果Fig. 5 Visual simulation results of the interaction effect between the lead cable and the mooring ball

    2.3 海上補給仿真訓(xùn)練時效性計算

    在確定PBD 算法后,還需確定仿真迭代計算方法,用于確保在較少的迭代步數(shù)下,模型計算盡快收斂。PBD 算法通常采用非線性高斯-賽德爾或牛頓迭代法兩種迭代算法進行計算。本文通過以下算例比較非線性高斯-賽德爾或牛頓迭代法兩種迭代算法,確定海上補給柔性體仿真中迭代算法。

    海上補給仿真訓(xùn)練系統(tǒng)將柔性體分段長度設(shè)為0.2 m,重力加速度為9.8 m/s2,比較非線性高斯-賽德爾迭代和牛頓迭代梁總算法的相對誤差與迭代次數(shù)[16],發(fā)現(xiàn)迭代次數(shù)相同的情況下,高斯-賽德爾迭代法在收斂速度上優(yōu)于牛頓法,如表1 所示,采用兩種不同的方法計算得到不同的迭代次數(shù)下所用的時間和對應(yīng)的相對誤差。從表1 中可以看出,在 5 次迭代計算時,PBD 方法得到的相對誤差為9.02×10-3、耗時為28 ms,而牛頓迭代法得到的相對誤差為4.5×10-2、耗時為422 ms。兩種方法在5 次迭代內(nèi)均已經(jīng)能夠滿足仿真精度的需求,牛頓迭代算法更加耗時。因此,本文采用的高斯-賽德爾迭代方法在較低的迭代次數(shù)時能較好地滿足仿真系統(tǒng)的要求,且表現(xiàn)出更好的穩(wěn)定性能。從仿真實時性的角度考慮,此方法在交互應(yīng)用中更具有優(yōu)勢。

    表1 PBD 方法的性能結(jié)果Table 1 Performance results of the PBD method

    2.4 海上補給仿真訓(xùn)練系統(tǒng)

    本文從多體動力學(xué)模型出發(fā),建立補給裝備索具柔性體模型,并運用PBD 算法開展海上補給裝備索具視景仿真研究,使用了拉伸約束、彎曲約束、長距離附著物約束和單邊引腳約束,避免了變形失真。通過提高索具、輸油軟管等柔性體的動態(tài)逼真度,復(fù)現(xiàn)海上補給實際狀態(tài),并以此為基礎(chǔ)研發(fā)了海上補給仿真訓(xùn)練系統(tǒng),用于開展海上航行橫向干、液貨補給仿真訓(xùn)練[17-20]。

    該系統(tǒng)為國內(nèi)首套用于海上補給訓(xùn)練的仿真系統(tǒng),可開展海上航行橫向補給全流程訓(xùn)練。通過部隊試用,反映出本文系統(tǒng)索具、輸油軟管等柔性體的動態(tài)逼真度極高,能夠很好地復(fù)現(xiàn)海上補給實際狀態(tài),訓(xùn)練效果基本與實裝訓(xùn)練一致,可有效提升海軍部隊海上補給的訓(xùn)練水平。

    3 結(jié)論

    本文針對海上補給仿真訓(xùn)練系統(tǒng)尚存在柔性體變形失真的問題,建立海上補給裝備索具柔性體多體動力學(xué)仿真模型,通過增加拉伸約束、彎曲約束、長距離約束和附著物相互約束實現(xiàn)了柔性體的動態(tài)變形三維模擬。最后通過比較高斯-賽德爾與牛頓算法的實時計算效果,確定本文采用的高斯-賽德爾迭代方法來提高視景仿真的效果,避免因為計算時間長導(dǎo)致畫面失真。最終開發(fā)出了海上補給仿真訓(xùn)練系統(tǒng),極大地提高了系統(tǒng)的操作真實感,為海上補給專業(yè)官兵盡快形成戰(zhàn)斗力提供了有效途徑。

    猜你喜歡
    索具附著物視景
    大型結(jié)構(gòu)物吊裝索具計算分析
    膨脹彎海上吊裝設(shè)計理論方法研究
    中國水運(2021年5期)2021-07-01 19:10:08
    兵馬俑修復(fù):為你,千千萬萬遍
    船底涂層有什么作用?
    vega及其在戰(zhàn)場環(huán)境視景仿真系統(tǒng)中的應(yīng)用
    電子測試(2018年4期)2018-05-09 07:28:19
    輸電線路附著物測算系統(tǒng)測算功能模塊的研究
    杜亞中:欲正人者先正己
    鍛造與沖壓(2015年7期)2015-06-22 14:38:15
    全回轉(zhuǎn)拖輪運動建模與視景仿真
    基于Vortex與Vega Prime的車輛模擬器動力學(xué)與視景仿真
    菹草(Potamogeton crispus)附著物對水體氮、磷負荷的響應(yīng)*
    街头女战士在线观看网站| 极品人妻少妇av视频| 国产精品三级大全| 国产淫片久久久久久久久| 亚洲国产成人一精品久久久| 中国美白少妇内射xxxbb| 全区人妻精品视频| 女人精品久久久久毛片| 久久久久久久精品精品| 哪个播放器可以免费观看大片| 欧美3d第一页| 三级国产精品片| 日本91视频免费播放| 大香蕉久久网| 哪个播放器可以免费观看大片| 9色porny在线观看| 超碰97精品在线观看| 日韩精品免费视频一区二区三区 | 视频中文字幕在线观看| 免费看日本二区| 日韩在线高清观看一区二区三区| 亚洲精品,欧美精品| 亚洲怡红院男人天堂| 91aial.com中文字幕在线观看| 精品久久久久久久久亚洲| 伊人亚洲综合成人网| 国产精品一区二区在线观看99| .国产精品久久| 日韩av在线免费看完整版不卡| xxx大片免费视频| 国模一区二区三区四区视频| 国产乱人偷精品视频| 岛国毛片在线播放| 亚洲综合色惰| 婷婷色综合www| 九九爱精品视频在线观看| 内射极品少妇av片p| 男的添女的下面高潮视频| 亚洲自偷自拍三级| 亚洲综合精品二区| 青春草视频在线免费观看| 国产成人freesex在线| 在线观看免费日韩欧美大片 | 国产精品久久久久久精品电影小说| 国产亚洲av片在线观看秒播厂| 成人亚洲精品一区在线观看| 国产伦精品一区二区三区视频9| 久久6这里有精品| 国产精品久久久久久精品古装| 国产伦理片在线播放av一区| 国产午夜精品久久久久久一区二区三区| 狂野欧美激情性xxxx在线观看| 水蜜桃什么品种好| 国产毛片在线视频| 久久久久精品久久久久真实原创| videossex国产| 亚洲不卡免费看| 国产精品.久久久| av福利片在线| 丰满少妇做爰视频| 久久精品国产鲁丝片午夜精品| 看十八女毛片水多多多| 又黄又爽又刺激的免费视频.| 国产成人精品福利久久| 一本色道久久久久久精品综合| 精品国产国语对白av| 久久精品国产亚洲av天美| 国产成人精品福利久久| 亚洲怡红院男人天堂| 美女内射精品一级片tv| 在线观看一区二区三区激情| 午夜福利在线观看免费完整高清在| 亚洲精品色激情综合| 能在线免费看毛片的网站| 日日啪夜夜爽| 国产成人精品无人区| 久久人人爽人人爽人人片va| 黄色视频在线播放观看不卡| 日韩成人伦理影院| 亚洲国产av新网站| 在线 av 中文字幕| 久久亚洲国产成人精品v| 成人毛片a级毛片在线播放| 欧美亚洲 丝袜 人妻 在线| 国产精品成人在线| 久久久久精品久久久久真实原创| 免费播放大片免费观看视频在线观看| 天美传媒精品一区二区| 亚洲国产成人一精品久久久| 插逼视频在线观看| 一级片'在线观看视频| 亚洲欧美中文字幕日韩二区| 欧美老熟妇乱子伦牲交| 日本av免费视频播放| 久久99热这里只频精品6学生| 丝袜在线中文字幕| 伦理电影大哥的女人| www.色视频.com| 能在线免费看毛片的网站| 亚洲国产毛片av蜜桃av| 亚洲精品乱码久久久久久按摩| 日韩强制内射视频| 欧美 日韩 精品 国产| 九草在线视频观看| 国产精品人妻久久久影院| 中文乱码字字幕精品一区二区三区| 国产欧美日韩精品一区二区| 如日韩欧美国产精品一区二区三区 | 久久久亚洲精品成人影院| 国产成人精品婷婷| 亚洲激情五月婷婷啪啪| 18禁动态无遮挡网站| 中文字幕免费在线视频6| av线在线观看网站| 国产午夜精品一二区理论片| 成人毛片a级毛片在线播放| av免费在线看不卡| 国产成人免费观看mmmm| 狠狠精品人妻久久久久久综合| 久久久久久人妻| 香蕉精品网在线| av免费观看日本| 天天躁夜夜躁狠狠久久av| av网站免费在线观看视频| 街头女战士在线观看网站| 草草在线视频免费看| 日本色播在线视频| 国产成人精品婷婷| 国产成人a∨麻豆精品| 国产美女午夜福利| 女的被弄到高潮叫床怎么办| 亚洲av成人精品一区久久| 九草在线视频观看| 日韩人妻高清精品专区| 国产成人免费观看mmmm| 午夜免费鲁丝| 精品人妻熟女av久视频| 国产日韩欧美亚洲二区| 人妻制服诱惑在线中文字幕| 交换朋友夫妻互换小说| 色视频www国产| tube8黄色片| 日韩,欧美,国产一区二区三区| 久久人妻熟女aⅴ| 久久久久久久久久久久大奶| 观看美女的网站| 丝瓜视频免费看黄片| 欧美日韩一区二区视频在线观看视频在线| 国产男人的电影天堂91| 免费看日本二区| 最近中文字幕2019免费版| 国产国拍精品亚洲av在线观看| 国产美女午夜福利| 国产精品一区二区在线观看99| 啦啦啦在线观看免费高清www| 免费黄色在线免费观看| 女性生殖器流出的白浆| 久久久精品94久久精品| 国产乱人偷精品视频| 永久免费av网站大全| 国产日韩一区二区三区精品不卡 | 国产中年淑女户外野战色| 不卡视频在线观看欧美| 啦啦啦在线观看免费高清www| 久久99热6这里只有精品| 三级经典国产精品| 精品久久国产蜜桃| 黑人巨大精品欧美一区二区蜜桃 | 少妇高潮的动态图| 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 午夜福利在线观看免费完整高清在| av天堂久久9| 纯流量卡能插随身wifi吗| 欧美三级亚洲精品| 久久午夜综合久久蜜桃| 极品少妇高潮喷水抽搐| 国产 一区精品| 能在线免费看毛片的网站| 午夜视频国产福利| 看非洲黑人一级黄片| 国产精品嫩草影院av在线观看| 午夜福利,免费看| 亚洲综合色惰| 男人狂女人下面高潮的视频| 国产免费一级a男人的天堂| 久久久久国产精品人妻一区二区| 涩涩av久久男人的天堂| 亚洲熟女精品中文字幕| 老司机亚洲免费影院| 国产精品不卡视频一区二区| 亚洲美女视频黄频| 国产精品国产三级国产专区5o| 女性生殖器流出的白浆| 久久久久人妻精品一区果冻| 亚洲美女搞黄在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲综合精品二区| av有码第一页| 一级二级三级毛片免费看| 新久久久久国产一级毛片| 91午夜精品亚洲一区二区三区| 国产免费视频播放在线视频| 亚洲欧美成人精品一区二区| 久久人人爽人人爽人人片va| av线在线观看网站| 一区二区三区精品91| 激情五月婷婷亚洲| av不卡在线播放| 97精品久久久久久久久久精品| 免费少妇av软件| 精品少妇黑人巨大在线播放| 日产精品乱码卡一卡2卡三| 一区二区av电影网| 国产免费视频播放在线视频| 精品99又大又爽又粗少妇毛片| 日韩,欧美,国产一区二区三区| 看非洲黑人一级黄片| 卡戴珊不雅视频在线播放| 中文字幕免费在线视频6| 精品酒店卫生间| 亚洲美女视频黄频| 日韩免费高清中文字幕av| 精品亚洲成国产av| 你懂的网址亚洲精品在线观看| 日韩中文字幕视频在线看片| 中文乱码字字幕精品一区二区三区| 午夜免费观看性视频| 男女啪啪激烈高潮av片| 丰满迷人的少妇在线观看| 亚洲精品乱码久久久v下载方式| 亚洲激情五月婷婷啪啪| 国产老妇伦熟女老妇高清| 国产成人精品福利久久| a级毛片免费高清观看在线播放| 日日啪夜夜爽| 69精品国产乱码久久久| 一级黄片播放器| 久久久久久久久久人人人人人人| av在线播放精品| 精品视频人人做人人爽| 日日撸夜夜添| 纯流量卡能插随身wifi吗| 一级,二级,三级黄色视频| 热re99久久国产66热| 久久亚洲国产成人精品v| 亚洲精品色激情综合| 中文字幕亚洲精品专区| 91aial.com中文字幕在线观看| 女性被躁到高潮视频| 亚洲欧美日韩东京热| 男女边吃奶边做爰视频| 成人综合一区亚洲| 最新中文字幕久久久久| 一级av片app| a级毛片在线看网站| 精品卡一卡二卡四卡免费| 亚洲无线观看免费| 97在线视频观看| 97在线人人人人妻| 久久精品久久久久久噜噜老黄| 欧美精品亚洲一区二区| 亚洲欧美一区二区三区黑人 | 欧美三级亚洲精品| 91精品伊人久久大香线蕉| 国产精品国产av在线观看| 国产成人免费观看mmmm| av线在线观看网站| 视频区图区小说| 久久人人爽人人片av| 夜夜看夜夜爽夜夜摸| 亚洲,一卡二卡三卡| 两个人免费观看高清视频 | kizo精华| 亚洲精品国产成人久久av| 又黄又爽又刺激的免费视频.| 日韩伦理黄色片| 午夜激情福利司机影院| a级毛片在线看网站| 高清欧美精品videossex| 国产欧美亚洲国产| 久久青草综合色| 亚洲第一av免费看| 国产免费福利视频在线观看| 又大又黄又爽视频免费| 国产精品久久久久久av不卡| 精品99又大又爽又粗少妇毛片| 精品熟女少妇av免费看| 亚洲精品国产色婷婷电影| 国产一区有黄有色的免费视频| 视频区图区小说| 性色av一级| 美女脱内裤让男人舔精品视频| 最新中文字幕久久久久| 校园人妻丝袜中文字幕| 老司机影院毛片| 大陆偷拍与自拍| 一边亲一边摸免费视频| 天堂俺去俺来也www色官网| 日韩视频在线欧美| 三级国产精品欧美在线观看| 国产成人免费观看mmmm| 简卡轻食公司| 亚洲国产日韩一区二区| 十分钟在线观看高清视频www | 亚洲国产精品国产精品| 三级国产精品欧美在线观看| 成人18禁高潮啪啪吃奶动态图 | 欧美精品亚洲一区二区| 免费看光身美女| 久久久久久久久大av| 欧美xxxx性猛交bbbb| 精品少妇黑人巨大在线播放| 最近最新中文字幕免费大全7| 国产精品一区www在线观看| 日本黄色日本黄色录像| 精品午夜福利在线看| 国产在线一区二区三区精| 美女主播在线视频| 美女视频免费永久观看网站| 亚洲欧美精品专区久久| 久久午夜福利片| 国产一级毛片在线| 久久精品久久精品一区二区三区| 国产精品国产三级国产专区5o| 日韩成人伦理影院| 久久99蜜桃精品久久| 成人无遮挡网站| 99久久精品一区二区三区| 国产真实伦视频高清在线观看| 18禁在线无遮挡免费观看视频| 99久久精品一区二区三区| 国产乱来视频区| 不卡视频在线观看欧美| 日韩三级伦理在线观看| 久久人人爽人人片av| 久久久久久久久久久免费av| 免费不卡的大黄色大毛片视频在线观看| 欧美xxxx性猛交bbbb| 精品人妻熟女毛片av久久网站| 亚洲av男天堂| 乱人伦中国视频| 日韩中字成人| 18禁动态无遮挡网站| 亚洲熟女精品中文字幕| 晚上一个人看的免费电影| 一区二区三区精品91| 久久精品国产亚洲网站| 国产精品福利在线免费观看| 国产黄色视频一区二区在线观看| 国产av一区二区精品久久| 久久精品久久久久久噜噜老黄| 一本一本综合久久| 久久久精品免费免费高清| 久久精品国产自在天天线| 91久久精品国产一区二区三区| 国产日韩欧美视频二区| 大香蕉久久网| 国产精品免费大片| 纯流量卡能插随身wifi吗| 男女免费视频国产| 天美传媒精品一区二区| 日韩视频在线欧美| 国产熟女午夜一区二区三区 | av免费观看日本| 精品一区二区免费观看| 搡老乐熟女国产| 高清毛片免费看| 欧美精品一区二区大全| 亚洲,一卡二卡三卡| 国产午夜精品一二区理论片| 亚洲电影在线观看av| 多毛熟女@视频| 自线自在国产av| 国产黄频视频在线观看| 能在线免费看毛片的网站| 久久久久久久久久久免费av| 一边亲一边摸免费视频| 欧美成人精品欧美一级黄| 久久99一区二区三区| 久久 成人 亚洲| 精品一品国产午夜福利视频| 亚洲av.av天堂| 纯流量卡能插随身wifi吗| a级毛片在线看网站| 亚洲高清免费不卡视频| 免费少妇av软件| 日韩中字成人| 国产老妇伦熟女老妇高清| h日本视频在线播放| 一级毛片久久久久久久久女| 97在线人人人人妻| 街头女战士在线观看网站| 久久久久网色| 色网站视频免费| 天堂8中文在线网| 中文在线观看免费www的网站| 久久精品久久久久久噜噜老黄| 欧美 亚洲 国产 日韩一| 日韩在线高清观看一区二区三区| 亚洲怡红院男人天堂| 免费看日本二区| 亚洲欧美精品自产自拍| 91在线精品国自产拍蜜月| 如何舔出高潮| 街头女战士在线观看网站| 国产亚洲精品久久久com| 亚洲欧美日韩东京热| 美女国产视频在线观看| 亚洲av不卡在线观看| 日韩人妻高清精品专区| 亚洲国产最新在线播放| 成人美女网站在线观看视频| 日日撸夜夜添| 男人狂女人下面高潮的视频| a级毛片在线看网站| 99久国产av精品国产电影| 人妻制服诱惑在线中文字幕| 免费av中文字幕在线| a级一级毛片免费在线观看| 女人精品久久久久毛片| 嘟嘟电影网在线观看| 欧美精品人与动牲交sv欧美| 国产免费一区二区三区四区乱码| 另类亚洲欧美激情| 在线观看三级黄色| 欧美成人午夜免费资源| 中文欧美无线码| 一本大道久久a久久精品| 国产伦在线观看视频一区| av一本久久久久| 亚洲第一av免费看| 五月玫瑰六月丁香| 欧美3d第一页| 亚洲精品久久久久久婷婷小说| 人妻人人澡人人爽人人| 99re6热这里在线精品视频| 久久久久人妻精品一区果冻| 夜夜爽夜夜爽视频| 国产一区二区三区av在线| 国产亚洲精品久久久com| 久久久国产精品麻豆| 国产综合精华液| 亚洲欧美日韩东京热| 国产欧美另类精品又又久久亚洲欧美| 国产精品嫩草影院av在线观看| 蜜桃在线观看..| 亚州av有码| 免费黄色在线免费观看| 国产精品久久久久久久电影| 黄色毛片三级朝国网站 | 一区二区av电影网| 曰老女人黄片| 人妻制服诱惑在线中文字幕| 国国产精品蜜臀av免费| 国产91av在线免费观看| a级毛片在线看网站| 高清欧美精品videossex| 十八禁网站网址无遮挡 | 青春草视频在线免费观看| 91aial.com中文字幕在线观看| 下体分泌物呈黄色| 少妇高潮的动态图| 一级毛片久久久久久久久女| 国产又色又爽无遮挡免| 日本午夜av视频| 日本黄色片子视频| 纵有疾风起免费观看全集完整版| 偷拍熟女少妇极品色| 日本欧美视频一区| 蜜臀久久99精品久久宅男| 少妇猛男粗大的猛烈进出视频| 亚洲,欧美,日韩| 国产男人的电影天堂91| 国产精品蜜桃在线观看| 日韩欧美 国产精品| 亚洲va在线va天堂va国产| 久久女婷五月综合色啪小说| 亚洲av电影在线观看一区二区三区| 麻豆成人午夜福利视频| 成人免费观看视频高清| 国产黄片视频在线免费观看| av视频免费观看在线观看| 国产亚洲欧美精品永久| 99热这里只有精品一区| 欧美精品一区二区大全| 久久韩国三级中文字幕| 99久久精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看| a级毛色黄片| 日本午夜av视频| 国产爽快片一区二区三区| 精品国产一区二区久久| 国产视频内射| 国产一区二区三区av在线| 免费在线观看成人毛片| 国产永久视频网站| 下体分泌物呈黄色| 久久久久久久久久久久大奶| 久久精品夜色国产| av免费观看日本| 国产一区有黄有色的免费视频| av天堂久久9| .国产精品久久| 国产精品99久久99久久久不卡 | 人人澡人人妻人| 女人久久www免费人成看片| 99久久精品国产国产毛片| 欧美亚洲 丝袜 人妻 在线| 久久精品国产亚洲av天美| 国产精品三级大全| 五月开心婷婷网| 国产精品熟女久久久久浪| 国产精品成人在线| 国产精品麻豆人妻色哟哟久久| 一级a做视频免费观看| 看免费成人av毛片| av不卡在线播放| 欧美精品一区二区免费开放| 国产高清三级在线| 综合色丁香网| 寂寞人妻少妇视频99o| 超碰97精品在线观看| 97超视频在线观看视频| 欧美+日韩+精品| 久久国内精品自在自线图片| 欧美日韩在线观看h| 99热这里只有是精品50| 在线观看美女被高潮喷水网站| 欧美日韩一区二区视频在线观看视频在线| 国产亚洲一区二区精品| 亚洲精品乱久久久久久| 尾随美女入室| 国产亚洲精品久久久com| 欧美精品一区二区免费开放| 国产日韩欧美在线精品| 综合色丁香网| 我的老师免费观看完整版| 国产无遮挡羞羞视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产精品999| 能在线免费看毛片的网站| 蜜臀久久99精品久久宅男| 国产日韩欧美亚洲二区| 99九九线精品视频在线观看视频| 中文字幕人妻丝袜制服| 嘟嘟电影网在线观看| 国产成人精品婷婷| 女性生殖器流出的白浆| 少妇人妻精品综合一区二区| 国产一区二区在线观看av| 亚洲精华国产精华液的使用体验| 久久 成人 亚洲| 日韩av不卡免费在线播放| 久久精品国产鲁丝片午夜精品| av天堂久久9| 在线天堂最新版资源| 超碰97精品在线观看| 中文精品一卡2卡3卡4更新| 九草在线视频观看| 51国产日韩欧美| 亚洲怡红院男人天堂| 免费看不卡的av| 亚洲婷婷狠狠爱综合网| 只有这里有精品99| 日本欧美视频一区| 国产成人免费无遮挡视频| 成人毛片a级毛片在线播放| 亚洲精品乱码久久久v下载方式| 国产视频首页在线观看| 肉色欧美久久久久久久蜜桃| 免费看av在线观看网站| 高清av免费在线| 久久女婷五月综合色啪小说| 天堂俺去俺来也www色官网| 大码成人一级视频| 人妻制服诱惑在线中文字幕| 国产黄色视频一区二区在线观看| 色婷婷av一区二区三区视频| 人人妻人人爽人人添夜夜欢视频 | 久久国产精品大桥未久av | 99久久中文字幕三级久久日本| 国产精品99久久久久久久久| 国产精品久久久久久精品电影小说| 一区二区av电影网| 亚洲无线观看免费| 精品一区二区三区视频在线| 亚洲怡红院男人天堂| 色5月婷婷丁香| 中文乱码字字幕精品一区二区三区| 少妇人妻 视频| 免费看不卡的av| 国内精品宾馆在线| 亚洲欧洲国产日韩| 色5月婷婷丁香| 麻豆成人午夜福利视频| 日本wwww免费看| 国产成人一区二区在线| 亚洲国产精品成人久久小说| 亚洲伊人久久精品综合| 国产探花极品一区二区| 亚洲欧美一区二区三区国产| 中国三级夫妇交换| 黄色一级大片看看| 最近的中文字幕免费完整| 18禁裸乳无遮挡动漫免费视频| 黄色一级大片看看| 国产精品国产三级国产专区5o| 国产成人精品婷婷| 成人亚洲欧美一区二区av| 日本黄色片子视频| 日本猛色少妇xxxxx猛交久久| 精品国产露脸久久av麻豆| 18禁动态无遮挡网站| 在线观看国产h片|