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

    變周期減搖水艙控制與仿真

    2012-07-07 14:19:04孫寶國謝廣東
    中國艦船研究 2012年5期
    關(guān)鍵詞:水艙被動式氣閥

    孫寶國 謝廣東

    中國艦船研究設(shè)計(jì)中心,上海201102

    變周期減搖水艙控制與仿真

    孫寶國 謝廣東

    中國艦船研究設(shè)計(jì)中心,上海201102

    減搖水艙的設(shè)計(jì)首先要考慮與船體的適配,同時在設(shè)計(jì)中要重點(diǎn)考慮水艙的阻尼與周期特性,其中周期特性是水艙最基本的特性參數(shù),適用周期范圍越寬,減搖效果越理想。針對變周期減搖水艙的相當(dāng)長度與一般U型水艙相當(dāng)長度計(jì)算的不同之處——其連通道寬度可調(diào)且邊艙面積隨邊艙外壁角變化而變化,通過積分方法推導(dǎo)出了可變周期減搖水艙相當(dāng)長度的計(jì)算公式,提出了一種適用于變周期減搖水艙的估算方法。所設(shè)計(jì)的變周期減搖水艙采用連通道擋板控制改變周期,避免了氣閥控制的時間延遲問題,還可根據(jù)船體運(yùn)動規(guī)律調(diào)整連通道阻尼擋板以提高減搖效果。以某型船為例,建立了船舶減搖水艙仿真模型。經(jīng)減搖水艙自由衰減振蕩數(shù)值仿真方法驗(yàn)證,所提出的周期估算公式可靠。采用該模型研究了變周期減搖水艙的控制方法,結(jié)果表明,在其工作范圍內(nèi),可使橫搖角始終保持在8°以內(nèi),可以在5.5~18 s的周期范圍內(nèi)有效減搖,大大超過了被動式減搖水艙的工作頻帶。

    變周期減搖水艙;固有周期;自由振蕩;減搖效果

    0 引 言

    減搖水艙是目前船舶最常用的減搖裝置之一,依靠橫搖運(yùn)動所產(chǎn)生的能量使水艙內(nèi)的水流動,從而產(chǎn)生減搖力矩以減少船舶的橫搖。它具有結(jié)構(gòu)簡單,造價低廉,便于維護(hù)與保養(yǎng)等優(yōu)點(diǎn),目前普遍應(yīng)用于RORO船、科學(xué)考察船、運(yùn)輸船等。

    減搖水艙依據(jù)控制方式可分為被動式減搖水艙、可控被動式減搖水艙和主動式減搖水艙。被動式減搖水艙由于結(jié)構(gòu)形式固定,減搖頻帶較窄,減搖能力相對較差;主動式減搖水艙由于消耗能量較多,目前應(yīng)用很少;可控被動式減搖水艙通過一定的控制,拓寬了減搖的頻帶,同時耗能極小,因此被廣泛應(yīng)用。國際上,可控被動式減搖水艙主要有兩種形式,一種是INTERING氣閥式可控被動式水艙,采用氣閥開關(guān)控制空氣連通道,延長水艙的固有周期;另一種為JFE和Stabilio的可變周期減搖水艙,通過擋板控制水連通道的連通面積,拓展減搖周期的范圍,通過控制不同位置和數(shù)量的擋板開關(guān)來調(diào)節(jié)減搖水艙的工作狀態(tài),使其適應(yīng)不同的航行環(huán)境,達(dá)到最佳的減搖效果[1-2]。

    本文通過建立“船舶—減搖水艙”數(shù)學(xué)仿真模型,首先利用計(jì)算流體力學(xué)的方法,計(jì)算不同狀態(tài)下減搖水艙的固有周期,然后根據(jù)船舶的運(yùn)動規(guī)律設(shè)計(jì)控制策略,基于某船完成變周期減搖水艙減搖效果計(jì)算,進(jìn)行減搖能力分析。

    1 數(shù)學(xué)模型

    1.1 海浪仿真模型

    海浪運(yùn)動屬于隨機(jī)過程,多以功率譜的形式表示頻率特性。至今所提出的大多數(shù)譜都具有如下結(jié)構(gòu)形式:

    式中,A及B中包含風(fēng)要素(風(fēng)速、風(fēng)時、風(fēng)距)或波要素(波高、周期)。文中仿真采用P-M譜,式中,g為重力加速度,H1/3為有義波高。

    海浪對船舶的橫搖干擾力矩主要與波傾角有關(guān),在研究船舶的橫搖運(yùn)動及控制系統(tǒng)時,需要考慮海浪波傾角αe對船舶及其控制系統(tǒng)的影響,海浪波傾角可視為零均值的平穩(wěn)隨機(jī)過程。

    波傾角譜和波高頻譜的關(guān)系為:

    式中,K?為有效波傾系數(shù)。

    當(dāng)船以一定航速和浪向航行時,波浪實(shí)際作用于船體上的頻率已不同于自然頻率,而表現(xiàn)為遭遇頻率ωe。根據(jù)能量等效原則,遭遇頻率能量譜密度與自然頻率能量譜密度函數(shù)之間有如下關(guān)系:

    式中,V為航速,kn;β為遭遇角,(°)。

    根據(jù)能量守恒,波傾角遭遇譜如下所示:

    時域內(nèi)的波傾角可以通過線性疊加法,采用下式進(jìn)行仿真:

    1.2 船舶水艙系統(tǒng)數(shù)學(xué)模型

    船舶橫搖運(yùn)動數(shù)學(xué)模型為:

    式中,φ為船舶橫搖角,(°);I1+Jt為船舶自身的質(zhì)量和附加質(zhì)量對x軸的慣性矩,(t·m2);D為船舶的排水量,t;為初穩(wěn)心高,m;c為邊艙軸線到船舶縱中剖面的距離,m;2Nφ為船舶的橫搖阻尼力矩比例系數(shù);S0為單邊艙的自由液面面積,m2;z值表示水艙右邊艙中的液位向上移動的位移,m;b2為水艙軸線對于船舶橫搖軸的靜矩,具有面積因次,可以表示為:式(6)中Jt可以表示為:

    水艙內(nèi)液體振蕩的微分方程式為:

    式中,2Nt為減搖水艙的橫搖阻尼力矩比例系數(shù);為水艙內(nèi)水柱的相當(dāng)長度,m。

    考慮到船舶在不規(guī)則波中的橫搖情況,減搖裝置的效果可以用減搖前和減搖后的三一橫搖幅值的比值來表示[3-4],稱為減搖比K,即

    式中,m0為橫搖能譜的0次矩;m2為橫搖能譜的2次矩;m4為橫搖能譜的4次矩。

    2 可變周期減搖水艙工作原理

    可變周期減搖水艙基本結(jié)構(gòu)如圖1所示。水艙邊艙用擋板隔開,分成Tank A和Tank B兩個部分,對應(yīng)設(shè)置兩個空氣連通道,用兩個閥F1和F2控制其開關(guān)。在底部連通道設(shè)置3個旋轉(zhuǎn)擋板D1,D2和D3。

    圖1 可變周期減搖水艙實(shí)物簡圖Fig.1 Physical diagram of variable-cycle anti-rolling tanks

    當(dāng)船舶橫搖頻率超出水艙的工作范圍時,水艙減搖效果不理想,甚至?xí)霈F(xiàn)增搖,為船舶的安全性帶來不利影響,需要通過氣閥F1和F2關(guān)閉減搖水艙。當(dāng)船舶橫搖頻率在水艙的工作范圍內(nèi)時,可以通過控制F1和F2兩個氣閥對A,B艙進(jìn)行切換。進(jìn)一步通過對水艙擋板D1,D2和D3的轉(zhuǎn)動進(jìn)行連續(xù)控制,改變水艙的固有周期。擋板葉片可以在0°~90°內(nèi)任意角度旋轉(zhuǎn),根據(jù)船舶的橫搖角、橫搖角速度等信息以及水艙內(nèi)液體的流動狀態(tài)預(yù)測系統(tǒng)運(yùn)動參數(shù)的變化,以此為反饋控制信號,利用葉片的角度調(diào)整水艙內(nèi)液體振蕩參數(shù),使液體的流動適應(yīng)船舶橫搖周期的變化,達(dá)到有效減搖的目的[5]。

    可變周期減搖水艙可以工作在表1所示的不同狀態(tài)下。

    表1 可變周期減搖水艙的4種工作狀態(tài)Tab.1 Four working states of variable-cycle anti-rolling tanks

    3 可變周期減搖水艙周期估算方法

    通過變周期減搖水艙的工作原理可知水艙固有周期的計(jì)算尤為重要。精確的周期估算對于研究水艙自身的特性,以及水艙的控制都很有必要。

    可變周期減搖水艙的周期估算公式為:

    式中,g為重力加速度,g=9.81 m/s2。

    式中,A0為邊艙截面積,m2。

    式中,Bc為水艙邊艙內(nèi)壁間距,m;Dc為連通道高度,m;Lc為連通道寬度,m;Lt為兩邊艙長度,m;Bt為水艙總寬度,m;h為靜水面深度,m;θ為邊艙傾斜角度,(°)。具體如圖2所示。

    圖2 可變周期減搖水艙結(jié)構(gòu)簡圖Fig.2 Configuration diagram of variable-cycle anti-rolling tank

    可變周期減搖水艙周期估算公式中水艙的相當(dāng)長度不同于一般U型水艙相當(dāng)長度的計(jì)算。其連通道寬度可調(diào),邊艙面積隨邊艙外壁角度的變化而變化??紤]到以上兩個方面的不同,通過積分的方法推導(dǎo)出了可變周期減搖水艙相當(dāng)長度的計(jì)算公式。

    可變周期減搖水艙也不同于一般氣閥控制的減搖水艙,通過水通道擋板控制避免了氣閥控制的時間延遲問題,使得控制更加精確、有效。為了便于分析,只考慮水艙工作在B,C,D狀態(tài)下,把水艙模型簡化成只有擋板控制的模型。

    為了更加精確地驗(yàn)證水艙的周期特性,采用VOF模型對減搖水艙進(jìn)行仿真分析。該問題的氣液兩相流連續(xù)方程為[6-7]:

    式中,υ為速度矢量,α表示水或空氣在單元體中的體積分?jǐn)?shù),它介于0和1之間。同時滿足:

    式中,下標(biāo)w和a分別表示水和空氣。

    動量方程為:

    式中,ρ和 μ表示單元的流體密度和動力學(xué)粘性系數(shù);p為壓力,MPa;g為重力加速度;F為體積力,kN。

    ρ由下式確定(μ的確定方法類似):

    4 仿真分析

    以某船為例,基本的船型數(shù)據(jù)如下。

    表2 基本的船型數(shù)據(jù)Tab.2 Basic ship dimensions

    根據(jù)優(yōu)化設(shè)計(jì)方法,設(shè)計(jì)了水艙的結(jié)構(gòu)以及在船體上的布置,如圖3所示。

    圖3 水艙在船體上的布置Fig.3 Configuration of tanks on the ship

    4.1 周期仿真驗(yàn)證

    本節(jié)利用減搖水艙自由衰減振蕩數(shù)值仿真的方法對水艙的周期計(jì)算公式進(jìn)行驗(yàn)證。分別針對上述3種工作狀態(tài)下的水艙進(jìn)行了自由衰減振蕩仿真研究,水艙自由振蕩初始條件為左邊艙液位高度高于右邊艙,并通過在邊艙內(nèi)設(shè)置波高儀的方法對邊艙液位進(jìn)行實(shí)時監(jiān)控。圖4給出了水艙狀態(tài)C時自由振蕩數(shù)值仿真氣液分布云圖。

    圖4 狀態(tài)C自由振蕩過程水艙氣液分布云圖Fig.4 The gas-liquid contours of tanks in the process of free oscillation at state C

    經(jīng)過仿真計(jì)算,可以得出如圖5所示的右邊艙液位隨時間變化曲線圖,從而可以計(jì)算出水艙周期。經(jīng)仿真與公式計(jì)算得出的水艙周期如表3所示。

    圖5 右邊艙液位隨時間變化曲線Fig.5 Time history of the liquid level in right side tank

    4.2 減搖效果仿真驗(yàn)證

    根據(jù)仿真結(jié)果,可以得出水艙有效減搖周期范圍:狀態(tài)B,7~11 s;狀態(tài)C,9~15 s;狀態(tài)D,11~17 s。

    通過監(jiān)測船舶的橫搖角,可以計(jì)算出船舶的平均橫搖周期。如果平均橫搖周期不在減搖系統(tǒng)設(shè)計(jì)的工作范圍內(nèi),則水艙關(guān)閉;如果平均橫搖周期在水艙工作范圍內(nèi),則根據(jù)船舶橫搖周期確定水艙的控制方案。每隔一段固定的時間統(tǒng)計(jì)一次船舶平均橫搖周期,不斷選擇和調(diào)整可變周期減搖水艙的工作狀態(tài),使其始終發(fā)揮良好的減搖作用。

    根據(jù)上述方案,建立“船舶—減搖水艙”仿真模型,分別對安裝可變周期減搖水艙與被動式減搖水艙的橫搖角與橫搖減搖效果進(jìn)行仿真[8-10],得到的仿真結(jié)果曲線如圖6和圖7所示。

    5 結(jié) 論

    本文研究了變周期減搖水艙的工作原理,推導(dǎo)出了變周期減搖水艙的周期估算公式,并應(yīng)用數(shù)值仿真的方法對這一公式的可靠性進(jìn)行了驗(yàn)證。從仿真結(jié)果可知:

    圖6 安裝可變周期減搖水艙和被動式減搖水艙時的橫搖響應(yīng)對比圖Fig.6 The rolling response comparison between variable-cycle anti-rolling tanks and passive anti-rolling tanks

    圖7 安裝可變周期減搖水艙和被動式減搖水艙時的減搖效果對比圖Fig.7 Stabilized effect comparison between variable-cycle anti-rolling tanks and passive anti-rolling tanks

    1)通過水艙自由振蕩仿真計(jì)算的結(jié)果與周期估算公式計(jì)算的結(jié)果相差非常小,僅為3.7%左右。因此,利用周期估算公式對水艙的周期進(jìn)行計(jì)算是可行的。

    2)從橫搖響應(yīng)上看,可變周期減搖水艙在其工作范圍內(nèi)可使橫搖角始終保持在8°以內(nèi),具有良好的減搖效果。

    3)從減搖范圍上看,可變周期減搖水艙可以在5.5~18 s的周期范圍內(nèi)有效減搖,大大超出了被動式減搖水艙的工作頻帶,而且不產(chǎn)生增搖效果。因此,合理控制的可變周期減搖水艙性能穩(wěn)定,安全實(shí)用。

    [1] STIGTER C.The performance of U-tanks as a passive anti-rolling device[J].International Shipbuilding Prog? ress,1966,13:249-275.

    [2] MOALEJI R,GREIG A R.On the development of ship anti-roll tanks[J].Ocean Engineering,2007,34(1):103-121.

    [3] 賈正余,黃宇,曾啟盛,等.可控式被動減搖水艙氣道的阻尼[J].船舶工程,2011,33(3):41-44.

    JIA Z Y,HUANG Y,ZENG Q S,et al.Damping of air duct of passive controlled anti-rolling tanks[J].Ship Engineering,2011,33(3):41-44.

    [4] 梁利華,劉志坤,張松濤.減搖水艙的固有頻率和阻尼研究與試驗(yàn)分析[J].船舶工程,2010,32(3):35-38.

    LIANG L H,LIU Z K,ZHANG S T.Study and experi?ment analysis on the natural frequency and damping of anti-rolling tank[J].Ship Engineering,2010,32(3):35-38.

    [5] PHAIROH T,HANG J K.Adaptive ship roll mitigation by using a U-tube tank[J].Ocea Engineering,2007,34(3/4):403-415.

    [6] 張虹,姜德龍,梁利華,等.減搖水艙試驗(yàn)裝置研究[J].機(jī)床與液壓,2009(12):77-79.

    ZHANG H,JIANG D L,LIANG L H,et al.The re?search on test equipment of anti-rolling tank[J].Ma?chine Tool and Hydraulics,2009(12):77-79.

    [7] 于立君,金鴻章,王輝.減搖水艙半物理仿真裝置的試驗(yàn)研究[J].哈爾濱工程大學(xué)學(xué)報(bào),2008,29(6):635-640.

    YU L J,JIN H Z,WANG H.Research on experiment of anti-rolling tank semi-physical simulation equipment[J].Journal of Harbin Engineering University,2008,29(6):635-640.

    [8] 蔣頡,丁勇.減搖水艙擋板開度對穩(wěn)定矩影響的研究[J].中國艦船研究,2011,6(3):36-39.

    JIANG J,DING Y.Stabilizing moment influenced by angle of damping plate of anti-rolling tank[J].Chinese Journal of Ship Research,2011,6(3):36-39.

    [9] 楊民.航行狀態(tài)下減搖水艙減搖能力仿真研究[J].自動化技術(shù)與應(yīng)用,2008,27(9):1-5.

    YANG M.Simulation study of the efficiency of the an?ti-rolling tank[J].Techniques of Automation and Ap?plications,2008,27(9):1-5.

    [10] 金鴻章,張宏瀚.基于被動式PD控制策略的“雙水艙”減搖系統(tǒng)研究[J].船舶力學(xué),2011,15(10):1082-1089.

    JIN H Z,ZHANG H H.Research on double-tank sta?bilizers based on passive PD control strategy[J].Jour?nal of Ship Mechanics,2011,15(10):1082-1089.

    [責(zé)任編輯:喻 菁]

    Control and Simulation of Variable-Cycle Anti-Rolling Tanks

    SUN Bao-guo XIE Guang-dong
    Shanghai Division,China Ship Development and Design Center,Shanghai 201102,China

    During the design of marine anti-rolling tanks,the suitability of the tank to the ship should be a priority,and the damping and cycle characteristics should be considered at the same time.Particularly,the tank cycle is a fundamental characteristic parameter:the wider range of application it provides,the su?perior stabilized effect is obtained.However,the calculation of the considerable length for variable-cycle anti-rolling tanks is distinct from that for general U type tanks,where the width of water channels for vari?able-cycle anti-rolling tanks can be adjusted optionally,and the area of side tanks changes along with the outer corner.In view of the problem,this paper proposed a calculation formula for the considerable length of variable-cycle anti-rolling tanks through the integral approach.The target tank adopted water channel baffle control to change the cycle,which avoids time delay caused by the air valve control,and the optimal stabilized effect could be thus realized by adjusting the damping baffle according to ship motion.To vali?date the algorithm,a certain ship model was established with free oscillation numerical simulations con?ducted.The results show that the roll angle remains within a range of±8°,and the roll stabilization is effec?tive with a period range of 5.5 s to 18 s,which significantly exceeds the working band of passive anti-roll?ing tanks.

    variable-cycle anti-rolling tank;natural period;free oscillation;stabilized effect

    U662.2

    A

    1673-3185(2012)05-22-05

    10.3969/j.issn.1673-3185.2012.05.004

    2012-02-03

    孫寶國(1970-),男,高級工程師。研究方向:艦船總體研究與設(shè)計(jì)。E?mail:sun1bao1guo@163.com

    謝廣東(1985-),男,碩士,助理工程師。研究方向:艦船總體研究與設(shè)計(jì)。E?mail:xieguangdong0@sina.com

    謝廣東。

    猜你喜歡
    水艙被動式氣閥
    阻尼被動減搖水艙流體動力特性分析及效果評估
    U 型減搖水艙的流體動力特性分析及控制系統(tǒng)設(shè)計(jì)
    可變周期被動式減搖水艙設(shè)計(jì)與仿真研究
    神奇的潛水艇
    被動式建筑節(jié)能技術(shù)
    被動式建筑節(jié)能技術(shù)
    簡化供暖在被動式公寓的應(yīng)用
    被動式建筑節(jié)能技術(shù)
    往復(fù)式壓縮機(jī)氣閥泄漏故障診斷處理一例
    循環(huán)氫壓縮機(jī)氣閥故障分析
    黄色一级大片看看| www.色视频.com| 国产精品精品国产色婷婷| 夜夜骑夜夜射夜夜干| 晚上一个人看的免费电影| 日韩电影二区| 九九在线视频观看精品| 97超视频在线观看视频| 亚洲国产最新在线播放| 午夜精品国产一区二区电影| 岛国毛片在线播放| 在线观看一区二区三区| 国产精品女同一区二区软件| 欧美3d第一页| 精品一区二区三卡| 久久鲁丝午夜福利片| 日韩一区二区视频免费看| 成人漫画全彩无遮挡| 在线观看免费高清a一片| 国产精品久久久久久精品电影小说 | 国产亚洲午夜精品一区二区久久| 国产男女内射视频| 亚洲国产欧美在线一区| 九色成人免费人妻av| av一本久久久久| 亚洲精品久久午夜乱码| 三级国产精品片| 久久毛片免费看一区二区三区| 1000部很黄的大片| 成人特级av手机在线观看| 成年人午夜在线观看视频| 免费少妇av软件| 日本免费在线观看一区| 久久人人爽人人爽人人片va| 狂野欧美白嫩少妇大欣赏| 久久婷婷青草| 国产乱人视频| 一级二级三级毛片免费看| 亚洲精品第二区| 国产一区有黄有色的免费视频| 高清午夜精品一区二区三区| 日本vs欧美在线观看视频 | 成年av动漫网址| 国产一区二区三区综合在线观看 | 在现免费观看毛片| 成人国产麻豆网| 国产一区二区三区av在线| 草草在线视频免费看| 又粗又硬又长又爽又黄的视频| 久久青草综合色| 国产精品嫩草影院av在线观看| 久久久久久久国产电影| 国产免费又黄又爽又色| 日韩一区二区视频免费看| 男女无遮挡免费网站观看| kizo精华| 国产精品一区二区三区四区免费观看| 99精国产麻豆久久婷婷| 国产真实伦视频高清在线观看| 久久久成人免费电影| 搡老乐熟女国产| 在线观看免费高清a一片| 麻豆成人午夜福利视频| 超碰97精品在线观看| 蜜桃久久精品国产亚洲av| 日韩电影二区| 国产极品天堂在线| av国产精品久久久久影院| 国产v大片淫在线免费观看| 国产在线视频一区二区| 久久99热这里只有精品18| av卡一久久| 内地一区二区视频在线| 日本-黄色视频高清免费观看| 18禁裸乳无遮挡动漫免费视频| 国产亚洲精品久久久com| 又黄又爽又刺激的免费视频.| 高清午夜精品一区二区三区| 国产亚洲欧美精品永久| 久久精品国产亚洲网站| 看十八女毛片水多多多| 日本vs欧美在线观看视频 | 久久久久性生活片| 国产在线免费精品| 成人国产av品久久久| 免费观看性生交大片5| 卡戴珊不雅视频在线播放| 久久久久网色| 男人狂女人下面高潮的视频| 夜夜爽夜夜爽视频| 久久人妻熟女aⅴ| 2018国产大陆天天弄谢| 成年免费大片在线观看| 性高湖久久久久久久久免费观看| 日韩中文字幕视频在线看片 | 一区二区三区精品91| 一级毛片久久久久久久久女| 亚洲国产色片| 妹子高潮喷水视频| 精品一品国产午夜福利视频| 国产精品偷伦视频观看了| 亚洲av电影在线观看一区二区三区| 国产亚洲一区二区精品| 久久女婷五月综合色啪小说| 一边亲一边摸免费视频| 天堂俺去俺来也www色官网| 麻豆成人av视频| 观看av在线不卡| 久久韩国三级中文字幕| 国产日韩欧美在线精品| 日产精品乱码卡一卡2卡三| 视频中文字幕在线观看| 欧美 日韩 精品 国产| 女性生殖器流出的白浆| 七月丁香在线播放| 少妇人妻精品综合一区二区| 国产欧美亚洲国产| 日日撸夜夜添| 少妇 在线观看| 麻豆成人av视频| 成人国产麻豆网| 中文字幕精品免费在线观看视频 | 免费大片18禁| 丝袜喷水一区| 自拍偷自拍亚洲精品老妇| 国产黄片视频在线免费观看| 99热这里只有是精品50| 精品人妻视频免费看| 久久6这里有精品| 国产男人的电影天堂91| 最近的中文字幕免费完整| 久久久久久久久久成人| 精品国产乱码久久久久久小说| videossex国产| 国产成人精品一,二区| av网站免费在线观看视频| 搡老乐熟女国产| 晚上一个人看的免费电影| 国产高清不卡午夜福利| 噜噜噜噜噜久久久久久91| 免费看日本二区| 女人久久www免费人成看片| 日本vs欧美在线观看视频 | 免费人成在线观看视频色| 亚洲精品日韩在线中文字幕| 亚洲成人手机| 赤兔流量卡办理| 美女主播在线视频| 久久久久网色| 看十八女毛片水多多多| 国产精品欧美亚洲77777| 国产亚洲一区二区精品| 国产成人精品久久久久久| 亚洲精品国产成人久久av| 蜜桃亚洲精品一区二区三区| 亚洲av日韩在线播放| 国产精品爽爽va在线观看网站| 嘟嘟电影网在线观看| 99精国产麻豆久久婷婷| 国产精品女同一区二区软件| 精品亚洲成国产av| 在线观看免费日韩欧美大片 | 免费大片18禁| 日韩视频在线欧美| 久久国产乱子免费精品| 99久久精品国产国产毛片| 小蜜桃在线观看免费完整版高清| 久久99热这里只有精品18| 男男h啪啪无遮挡| 一级av片app| 国产伦理片在线播放av一区| 免费人妻精品一区二区三区视频| 亚洲国产精品一区三区| 成年免费大片在线观看| 中国美白少妇内射xxxbb| 亚洲欧美精品自产自拍| 国产伦精品一区二区三区视频9| 观看美女的网站| 久久久久国产精品人妻一区二区| 青春草视频在线免费观看| 夜夜看夜夜爽夜夜摸| 免费av中文字幕在线| 成年免费大片在线观看| 欧美人与善性xxx| 国产真实伦视频高清在线观看| 久久久精品94久久精品| 麻豆成人午夜福利视频| videos熟女内射| 成人亚洲欧美一区二区av| 久久精品久久久久久噜噜老黄| 一区二区三区四区激情视频| 国产69精品久久久久777片| 久热久热在线精品观看| 97在线人人人人妻| 国产深夜福利视频在线观看| 国产成人精品福利久久| 久久99精品国语久久久| 欧美另类一区| 国产在视频线精品| 国产久久久一区二区三区| 91狼人影院| 久久99精品国语久久久| 18禁在线播放成人免费| 亚洲av成人精品一区久久| 国产极品天堂在线| 久久久久网色| 欧美日韩综合久久久久久| 久久久久人妻精品一区果冻| 观看免费一级毛片| 久久精品国产自在天天线| 国产淫片久久久久久久久| 内射极品少妇av片p| 爱豆传媒免费全集在线观看| 99热这里只有是精品在线观看| 免费av不卡在线播放| 国产视频首页在线观看| 老师上课跳d突然被开到最大视频| 国产极品天堂在线| 国产片特级美女逼逼视频| 少妇的逼水好多| 蜜臀久久99精品久久宅男| 日韩中文字幕视频在线看片 | 午夜免费男女啪啪视频观看| 午夜福利高清视频| 多毛熟女@视频| 国产中年淑女户外野战色| 少妇的逼水好多| 成人国产麻豆网| 狂野欧美白嫩少妇大欣赏| 免费av不卡在线播放| 噜噜噜噜噜久久久久久91| 老师上课跳d突然被开到最大视频| 亚洲精品,欧美精品| 婷婷色综合www| 精品少妇黑人巨大在线播放| 国产男人的电影天堂91| 内地一区二区视频在线| 亚洲天堂av无毛| 亚洲精品日本国产第一区| 少妇高潮的动态图| 欧美性感艳星| 小蜜桃在线观看免费完整版高清| 777米奇影视久久| 午夜精品国产一区二区电影| 黄片无遮挡物在线观看| 日日摸夜夜添夜夜爱| 国产在线一区二区三区精| 蜜臀久久99精品久久宅男| 高清欧美精品videossex| 不卡视频在线观看欧美| 中国国产av一级| 国产成人精品一,二区| 久久精品国产自在天天线| 少妇人妻精品综合一区二区| 国产女主播在线喷水免费视频网站| 国产免费视频播放在线视频| 欧美变态另类bdsm刘玥| 中文精品一卡2卡3卡4更新| 18禁动态无遮挡网站| 91午夜精品亚洲一区二区三区| 好男人视频免费观看在线| 色5月婷婷丁香| 中文精品一卡2卡3卡4更新| 欧美高清性xxxxhd video| 午夜福利在线在线| 国产 精品1| 最近中文字幕2019免费版| 综合色丁香网| 蜜桃在线观看..| 一本色道久久久久久精品综合| 久热久热在线精品观看| 中文欧美无线码| 国产 精品1| 你懂的网址亚洲精品在线观看| 天美传媒精品一区二区| av专区在线播放| 七月丁香在线播放| 亚洲国产精品国产精品| 亚洲激情五月婷婷啪啪| 人妻少妇偷人精品九色| 国产精品一区二区性色av| 国产精品久久久久久av不卡| 日本黄色日本黄色录像| 天堂8中文在线网| 精品视频人人做人人爽| 国产成人精品婷婷| 国产在线一区二区三区精| 欧美日韩视频精品一区| 人妻系列 视频| 国产高潮美女av| 欧美成人一区二区免费高清观看| 伦理电影大哥的女人| 激情 狠狠 欧美| 免费人成在线观看视频色| av免费在线看不卡| 亚洲国产成人一精品久久久| 午夜福利在线在线| 国产中年淑女户外野战色| av在线蜜桃| 欧美成人一区二区免费高清观看| av福利片在线观看| 亚洲欧美日韩卡通动漫| 久久毛片免费看一区二区三区| xxx大片免费视频| 久久精品国产a三级三级三级| 国产 一区精品| 韩国av在线不卡| 成人毛片60女人毛片免费| 国模一区二区三区四区视频| 日本黄色日本黄色录像| 各种免费的搞黄视频| 欧美丝袜亚洲另类| 美女脱内裤让男人舔精品视频| 国产精品免费大片| 人妻系列 视频| 欧美激情极品国产一区二区三区 | 在线观看免费高清a一片| 一级爰片在线观看| 男女边摸边吃奶| 日日摸夜夜添夜夜爱| 黄色欧美视频在线观看| 亚洲精品视频女| 亚洲色图综合在线观看| 国产色婷婷99| 免费观看a级毛片全部| 18禁裸乳无遮挡免费网站照片| 免费人成在线观看视频色| 2022亚洲国产成人精品| 国精品久久久久久国模美| 99热这里只有是精品50| av播播在线观看一区| 亚洲综合精品二区| 亚洲,欧美,日韩| 黑人高潮一二区| 建设人人有责人人尽责人人享有的 | 欧美变态另类bdsm刘玥| 搡女人真爽免费视频火全软件| 我要看黄色一级片免费的| 久热这里只有精品99| 婷婷色麻豆天堂久久| 欧美成人a在线观看| 日韩一区二区视频免费看| 国产一区二区在线观看日韩| 国内精品宾馆在线| 麻豆成人av视频| 国产欧美另类精品又又久久亚洲欧美| 最近最新中文字幕大全电影3| 亚洲国产精品专区欧美| 久久精品国产a三级三级三级| 男女下面进入的视频免费午夜| 蜜桃亚洲精品一区二区三区| 日韩伦理黄色片| 嘟嘟电影网在线观看| 欧美日韩在线观看h| 国产精品一二三区在线看| 这个男人来自地球电影免费观看 | 国产成人a∨麻豆精品| 又爽又黄a免费视频| 国产亚洲5aaaaa淫片| 成人无遮挡网站| 日韩一区二区视频免费看| 国产真实伦视频高清在线观看| 最后的刺客免费高清国语| av视频免费观看在线观看| 精品久久久久久久久亚洲| 色婷婷久久久亚洲欧美| 欧美激情国产日韩精品一区| 噜噜噜噜噜久久久久久91| 另类亚洲欧美激情| 久久久久人妻精品一区果冻| 亚洲久久久国产精品| 天美传媒精品一区二区| 亚洲中文av在线| 少妇裸体淫交视频免费看高清| 国产精品麻豆人妻色哟哟久久| 国产精品人妻久久久久久| 国产成人91sexporn| 最近2019中文字幕mv第一页| 国产在线视频一区二区| 有码 亚洲区| 蜜桃久久精品国产亚洲av| 这个男人来自地球电影免费观看 | 伊人久久国产一区二区| 亚洲精品自拍成人| 美女福利国产在线 | 久久国产乱子免费精品| 在线免费十八禁| 亚洲第一区二区三区不卡| 久久国内精品自在自线图片| 国产免费福利视频在线观看| 欧美精品亚洲一区二区| 欧美亚洲 丝袜 人妻 在线| 成年免费大片在线观看| 看非洲黑人一级黄片| 91精品国产九色| 少妇被粗大猛烈的视频| 99久久人妻综合| 99久久中文字幕三级久久日本| 精品亚洲成a人片在线观看 | 久久久久久久久久成人| 91精品国产国语对白视频| .国产精品久久| 国产在视频线精品| 欧美精品一区二区大全| 少妇猛男粗大的猛烈进出视频| 午夜免费男女啪啪视频观看| 免费不卡的大黄色大毛片视频在线观看| 啦啦啦啦在线视频资源| 久久国产精品大桥未久av | 国产精品久久久久久av不卡| 另类亚洲欧美激情| 丰满人妻一区二区三区视频av| av在线app专区| 亚洲精品色激情综合| 亚洲三级黄色毛片| av一本久久久久| 高清视频免费观看一区二区| 国产免费视频播放在线视频| 亚洲电影在线观看av| 国内精品宾馆在线| 男人舔奶头视频| 妹子高潮喷水视频| 亚洲美女视频黄频| 一级爰片在线观看| 国产精品爽爽va在线观看网站| 国产伦在线观看视频一区| 久久青草综合色| av福利片在线观看| 男人爽女人下面视频在线观看| 欧美精品人与动牲交sv欧美| 国产无遮挡羞羞视频在线观看| 久久精品国产自在天天线| 国产乱人偷精品视频| 亚洲不卡免费看| 久久久成人免费电影| 少妇精品久久久久久久| 男女边吃奶边做爰视频| 99国产精品免费福利视频| 久久 成人 亚洲| 最近中文字幕2019免费版| 2021少妇久久久久久久久久久| 22中文网久久字幕| 欧美3d第一页| 日本vs欧美在线观看视频 | 边亲边吃奶的免费视频| 久久精品久久久久久久性| 99热国产这里只有精品6| 极品少妇高潮喷水抽搐| 精品久久国产蜜桃| 不卡视频在线观看欧美| 国模一区二区三区四区视频| 国产爱豆传媒在线观看| 美女视频免费永久观看网站| 我的女老师完整版在线观看| 日韩一区二区视频免费看| 天堂8中文在线网| 日本黄色日本黄色录像| 婷婷色综合www| 中文字幕av成人在线电影| 韩国高清视频一区二区三区| 国内揄拍国产精品人妻在线| 欧美97在线视频| 国产亚洲最大av| 免费看光身美女| 欧美成人一区二区免费高清观看| 交换朋友夫妻互换小说| 国产精品一二三区在线看| 免费播放大片免费观看视频在线观看| 国产免费又黄又爽又色| 久久人人爽av亚洲精品天堂 | 国产片特级美女逼逼视频| 久久99精品国语久久久| 亚洲精品456在线播放app| 精品人妻视频免费看| 久久99热6这里只有精品| 亚洲av中文av极速乱| 人妻少妇偷人精品九色| 免费观看av网站的网址| 成人影院久久| 免费人成在线观看视频色| 国产精品国产三级国产av玫瑰| av在线app专区| 美女国产视频在线观看| 欧美高清性xxxxhd video| 久久精品国产自在天天线| 国产精品久久久久久久久免| 久久国产乱子免费精品| 国产视频内射| 亚洲一区二区三区欧美精品| 下体分泌物呈黄色| 你懂的网址亚洲精品在线观看| 国产免费一区二区三区四区乱码| 国产在线免费精品| 久久精品久久精品一区二区三区| 欧美人与善性xxx| 久久影院123| 久久毛片免费看一区二区三区| 97在线视频观看| 久久久精品免费免费高清| 精品一区在线观看国产| 亚洲欧美成人精品一区二区| 舔av片在线| 中文资源天堂在线| 丰满迷人的少妇在线观看| 欧美老熟妇乱子伦牲交| 伦理电影免费视频| 80岁老熟妇乱子伦牲交| 久久人人爽人人爽人人片va| 亚洲av成人精品一区久久| 国产精品国产三级国产av玫瑰| 99久久人妻综合| 亚洲怡红院男人天堂| 亚洲不卡免费看| 视频区图区小说| 日本黄色片子视频| 国产有黄有色有爽视频| 久久99热这里只频精品6学生| 久久亚洲国产成人精品v| 五月伊人婷婷丁香| 亚洲精品日本国产第一区| 三级国产精品欧美在线观看| 大香蕉97超碰在线| 欧美三级亚洲精品| 国内精品宾馆在线| 亚洲真实伦在线观看| 又爽又黄a免费视频| 国产精品爽爽va在线观看网站| 黄色日韩在线| 一级黄片播放器| 国产伦精品一区二区三区视频9| 99久久精品一区二区三区| 免费观看无遮挡的男女| 国产永久视频网站| 久久精品国产亚洲av天美| 有码 亚洲区| 性高湖久久久久久久久免费观看| 国产精品99久久99久久久不卡 | 国产乱人视频| 直男gayav资源| 国产视频首页在线观看| 人体艺术视频欧美日本| 国产精品一二三区在线看| 99久久中文字幕三级久久日本| 日日摸夜夜添夜夜添av毛片| 中文字幕精品免费在线观看视频 | 涩涩av久久男人的天堂| 美女xxoo啪啪120秒动态图| 欧美日韩一区二区视频在线观看视频在线| 久久精品久久久久久久性| 亚州av有码| 伊人久久精品亚洲午夜| 国产精品久久久久久精品电影小说 | 在线 av 中文字幕| 久久精品国产自在天天线| h视频一区二区三区| 国产精品久久久久久av不卡| 国产黄色免费在线视频| 极品少妇高潮喷水抽搐| 天堂8中文在线网| 男女免费视频国产| 在线亚洲精品国产二区图片欧美 | 亚洲一级一片aⅴ在线观看| 亚洲,欧美,日韩| 在线 av 中文字幕| 免费播放大片免费观看视频在线观看| 少妇的逼水好多| 一本久久精品| 亚洲欧美日韩无卡精品| 精品久久久噜噜| 91在线精品国自产拍蜜月| 日韩欧美 国产精品| 中文字幕免费在线视频6| 一级毛片电影观看| 熟女电影av网| 国产v大片淫在线免费观看| 有码 亚洲区| 久久国产精品大桥未久av | 一个人看视频在线观看www免费| av不卡在线播放| 少妇丰满av| 欧美高清成人免费视频www| 日本一二三区视频观看| 欧美日韩精品成人综合77777| 久热这里只有精品99| 最后的刺客免费高清国语| 国产真实伦视频高清在线观看| 亚洲三级黄色毛片| 99久久综合免费| av播播在线观看一区| 联通29元200g的流量卡| 精品久久久久久久末码| 欧美精品一区二区免费开放| 亚洲三级黄色毛片| 好男人视频免费观看在线| 黄片无遮挡物在线观看| 成人一区二区视频在线观看| videossex国产| 国产真实伦视频高清在线观看| 一级毛片aaaaaa免费看小| 欧美日韩国产mv在线观看视频 | 免费黄频网站在线观看国产| av在线播放精品| 亚洲精品自拍成人| 国产成人a区在线观看| 国产av国产精品国产| 久久久午夜欧美精品| 免费少妇av软件| 18禁在线无遮挡免费观看视频| 欧美精品一区二区免费开放| h视频一区二区三区| 亚洲va在线va天堂va国产| 国产一区二区在线观看日韩| 蜜桃在线观看..| 国产老妇伦熟女老妇高清| 亚洲美女视频黄频|