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

    基于風雪兩相流的高寒動車組轉向架防冰雪擾流技術研究*

    2016-10-31 07:46:31何德華王剛義陳厚嫦李超輝
    鐵道機車車輛 2016年4期
    關鍵詞:擾流擾流板轉向架

    何德華, 王剛義, 陳厚嫦, 李超輝

    (1 中國鐵道科學研究院 機車車輛研究所, 北京 100081;2 華東交通大學 機電工程學院, 江西南昌 330013)

    ?

    基于風雪兩相流的高寒動車組轉向架防冰雪擾流技術研究*

    何德華1, 王剛義2, 陳厚嫦1, 李超輝2

    (1中國鐵道科學研究院機車車輛研究所, 北京 100081;2華東交通大學機電工程學院, 江西南昌 330013)

    針對高寒動車組冬季運行時轉向架處存在冰雪堆積的問題,運用商用軟件STAR-CCM+和拉格朗日顆粒模型,采用更為直觀的風雪兩相流技術對擾流設計后轉向架區(qū)域冰雪附著情況進行模擬分析,并以壁面上的雪顆粒入射質量通量(IMF)來評價轉向架區(qū)域的擾流板防冰雪效果。計算結果表明:雪強、雪顆粒附著特性、車速、擾流板高度對轉向架冰雪附著率均有影響;擾流效果比較好的2種方案是,①改進的流線型擾流板方案可使轉向架冰雪附著速率降至無擾流板車型的54%,同時自身承受24%的附著量;②高200 mm的擾流板可使轉向架冰雪附著速率降至無擾流板車型的70%,同時自身承受9%的附著量。并且擾流板上的附著量不會結冰,會隨著氣流流動的方向流向地面。

    高寒動車組; 轉向架; 擾流板; 風雪兩相流;

    高寒動車組轉向架區(qū)域冰雪堆積是現(xiàn)階段轉向架防冰雪亟待解決的問題[1]。在國外,北山茂等的論文《基于轉向架區(qū)域氣動特性改進的研究》通過相關數(shù)值方法討論了在700系列車轉向架上增加擾流板的方法防冰雪侵害,北山茂等建立了轉向架區(qū)域的局部CFD模型,采用單相氣流的模擬輸入方法,以壓力系數(shù)為冰雪附著的表征變量,開展了幾組數(shù)值模擬分析轉向架擾流板的防冰雪效果[2];而我國現(xiàn)階段李俊民,單永林,林鵬的論文《高速動車組轉向架防冰雪導流罩的空氣動力學性能分析》通過改進車下的設備艙導流罩設計為傾斜結構,應用CFD對3編組動車組模型進行防風雪單向流分析[3]。相較于單向流模擬,兩相流有許多優(yōu)點,可更好的模擬現(xiàn)實環(huán)境,直觀的反映附著情況,這對于理解轉向架冰雪附著機制、改進轉向架擾流板防冰雪的設計有重要意義。目前,兩相流技術在汽車、飛機和建筑上的應用已經很廣泛[4-6],但應用于高寒動車組轉向架區(qū)域擾流技術,進行冰雪附著情況的研究尚未見報道。基于風雪兩相流計算模型,運用STAR-CCM+模擬轉向架區(qū)域對不同高度的擾流板的冰雪附著情況,以期為防冰雪擾流技術提供性能優(yōu)良的擾流板。

    1 控制方程和計算方法

    1.1控制方程

    描述流場流動除采用基本控制方程N-S方程[7],還要添加兩相流的控制方程。跟單相流相比,兩相流的計算不但涉及兩相間相互作用,其中的顆粒相與邊界的相互作用還涉及到運動分析,其復雜性和計算量高于單相流的計算。目前可用于風雪兩相流計算的多相流模型主要有兩種:(1) 歐拉兩相流混合模型;(2)拉格朗日顆粒模型。

    對于本文的數(shù)值計算,由于歐拉兩相流模型將第二相視為流體,從而較難描述壁面附著,且第二相顆粒尺寸受到網格分辨率的影響,難以用于合理地研究冰雪附著問題。因此使用拉格朗日顆粒模型,以突顯干濕雪顆粒在與壁面碰撞時的特點差異[8-9]。

    拉格朗日兩相流模型分為歐拉流體部分和拉格朗日顆粒部分,其中歐拉流體部分與常規(guī)流動方程類似,而拉格朗日顆粒部分的基本方程如下:

    (1)

    (2)

    (3)

    顆粒在遇到壁面后可能發(fā)生彈開和附著,相應的方程如下:

    (4)

    (5)

    1.2計算方法

    本文研究的問題可視為定常不可壓縮問題,其計算方法簡單描述為:采用基于RANS(雷諾平均法)的SST(剪切應力輸運)兩方程k-ω模型描述湍流行為;采用All y+方式描述壁面;采用有限體積法離散質量和動量控制方程,但模型中聯(lián)立能量控制方程,將過程近似視為等溫;采用二階迎風格式離散對流項和擴散項,避免數(shù)值假擴散帶來的誤差;采用分離式SIMPLE求解器處理離散所得方程組[10]。

    2 模型的建立

    2.1計算域的確定

    列車運行時的外流場是無限大的,但在數(shù)值模擬計算時,一般采用有限計算區(qū)域來代替無限計算區(qū)域。計算區(qū)域的選取很重要,要保證流場發(fā)展充分,應取得足夠大;而計算區(qū)域過大,網格數(shù)和計算時間將也會增大,因此要有效地準確地進行數(shù)值模擬,需選擇合適的計算區(qū)域。綜合這兩方面的考慮,采用的高寒動車組計算模型及計算區(qū)域如圖1所示,計算區(qū)域為長方體,計算區(qū)域的選為:車頭前100 m,車后250 m,車身上方100 m,車身左右各100 m。

    圖1 計算模型

    2.2計算模型及網格劃分

    考慮到計算工作量大,采用某型三編組的局部計算模型,如圖1所示。使用CFD軟件STAR-CCM+[11]進行模擬計算。根據試算,在本項目的計算中,采用Trim + Prism Layer的體網格類型進行網格劃分。對于網格劃分,在整個空間計算域采用較大尺寸的網格,在流場變化大的區(qū)域以及本例所關心的區(qū)域進行網格細化,在擾流板附近、中間車的轉向架及其周邊區(qū)域進行特別細化,從細網格到粗網格采用逐層過渡的方案,這樣就能保證計算要求也能提高計算效率。

    治理和防范高寒列車走行部冰雪附著問題有不少技術方案,最為根本和有效的是改變走行部的氣動條件,減少冰雪附著量,避免關鍵部位出現(xiàn)高強度附著。在不改變現(xiàn)有列車的基本結構條件下,在適當區(qū)域添加擾流裝置,影響走行部氣動條件,是一種簡單而高效的方法。所以,針對研究的某型動車組底部三維幾何結構,設計了與之相匹配的擾流結構,其幾何結構如圖2。通過試算,長度a和寬度h的長度取為定值2 000 mm和330 mm,而高b的值取0,100,200,300 mm。

    圖2 擾流板幾何結構

    2.3邊界條件的確定

    (1)入口邊界條件:由于計算量限制,僅能選取局部車體進行兩相流計算。通過分析可以發(fā)現(xiàn),頭車中部區(qū)域速度場結構較為穩(wěn)定,表現(xiàn)為外部流速較均勻,接近于自由來流流速,而車底流速則明顯較小,但在垂直于流向的平面上整體上較均勻。因此設置內外兩個列車風的速度入口條件,外部為高速區(qū),內部(主要是底部)為低速區(qū),經過試算發(fā)現(xiàn)這種設置所得下游流場與整車計算很相近。因此在每個入口區(qū)內,假設入口邊界來流的三維速度分布沒有受到模型的擾動,除運動方向外,另外兩個方向的速度分量為零,沿運動方向的速度在另外兩個方向上均勻分布,并且雪顆粒與周圍氣流具有相同流速。

    (2)出口邊界條件:取一個標準大氣壓強;雪顆??梢宰杂闪鞒?。

    (3)地面邊界條件:采用移動地板的做法,以消除地面效應的影響,假定高寒動車組風條件下地面雪顆粒揚起和沉積達到動態(tài)平衡,在地面處設置為顆粒和壁面完全彈性碰撞。

    (4)列車表面邊界條件:由于在列車表面存在邊界層效應的影響,故列車表面設定為有摩擦的固壁邊界(無滑移邊界),接近于實際情況,可較為精確的計算出列車表面的摩擦阻力,即黏性剪切力,以及表面的壓力分布等參數(shù);由于本研究需要準確計算轉向架流動,因此車輪的轉動不可忽略,在車輪和車軸上設置了轉動條件,轉動軸為相應車軸中心,轉動角速度由運行速度求出。

    3 轉向架風雪兩相流的模擬分析

    3.1局部模型流場分析

    在目前的計算條件下開展三編組整車兩相流計算是非常困難的,并且對于研究轉向架擾流結構這一局部問題也沒有必要。針對走行部流動特點和兩相流計算的實際條件,以三編組頭車后半段和中間車前半段的下半部分作為計算對象。為了檢測采用局部模型是否有效,采用流線型擾流板(長2 000 mm,高300 mm,寬330 mm)、高寒動車組車速為300 km/h、降雪強度為50 mm/24 h、濕雪的條件為例進行模擬分析。首先給出在車速300 km/h車體和轉向架表面壓力分布云圖和各個特征部位的壓力分布云圖,如下圖3。

    圖3 局部模型表面壓力分布云圖(300 km/h)

    由圖3可知,前轉向架區(qū)域的壓力峰值點出現(xiàn)在前輪緣中下方、制動盤中下部以及設備吊艙前緣等迎風面上,低壓區(qū)則分布在高壓區(qū)周圍;后輪壓力峰值點位置類似,但壓力值較小,輪軸附近存在一個低壓區(qū),由于擾流板和轉向架圍擋的遮蔽作用,轉向架上部表面并未出現(xiàn)明顯的高、低壓區(qū);從壓力云圖可以看出,擾流板和斜面圍擋在其身后形成了一個滯流區(qū),由于斜面的存在使得擾流板不能緊靠轉向架,氣流在擾流板后逐漸向上恢復,從而影響了滯流區(qū)的大小,同時在撞擊中軸后進一步向上,在轉向架上部區(qū)域產生了較大的流速;滯流區(qū)內存在一個局部環(huán)流結構,整個轉向架區(qū)域上部并未形成一個整體的環(huán)流,而是分散為多個環(huán)流。

    而后轉向架區(qū)域起到遮蔽效應的主要是直立的轉向架圍擋,且其效果比較明顯,后轉向架區(qū)域的壓力峰值主要出現(xiàn)在后輪輪緣、制動盤以及側方結構的迎風面處,低壓區(qū)出現(xiàn)在其附近,此時整個轉向架區(qū)域上部形成了一個較大的局部環(huán)流,雖然其流速值相對較小,但環(huán)流下部距離高速氣流很近,可能吸入含雪氣流,進而產生局部冰雪附著??梢钥闯?,斜面圍擋對后方車廂連接處的導流作用明顯,而擾流板由于尺寸和距離原因很難發(fā)揮作用,且其將會產生附帶的增阻影響,其增阻影響與其所阻滯的氣流速度分布有關。

    通過比較發(fā)現(xiàn),局部模型計算結果與整車模型計算結果很相近,這說明采用的局部模型符合要求并且設置的各種條件都滿足實際情況。

    3.2擾流板高度對流場的顆粒分布的模擬分析

    為了模擬出風雪兩相流的顆粒分布,通過壁面上的雪顆粒入射質量通量(Incident Mass Flux Snow,IMF)來評價轉向架區(qū)域冰雪的附著情況,在設置條件下該值等同于壁面處的顆粒附著速率。在300 km/h車速、50 mm/24 h雪強、濕雪條件下,圖4給出了不同擾流板高度條件下的顆粒附著速率分布。

    圖4 不同擾流板高度條件下的顆粒附著速率分布

    可以看到其與流動結構有明確的對應,前轉向架不受影響,完全相同。而隨著擾流板高度的增加,后轉向架上的附著速率呈減小趨勢,剎車片、吊艙等部件上的附著區(qū)域明顯減小,而與此同時,擾流板上的附著速率則明顯增加,這說明擾流板在替后方轉向架遮擋高速氣流的同時,本身也攔截了部分冰雪附著量,因此其作用是同時通過遮蔽(即導流)和代附著兩種機制發(fā)揮作用的。通過模擬計算分析,流線型擾流板(長2 000 mm,寬330 mm,高300 mm,高200 mm)方案為效果比較好,因為該擾流板方案可使轉向架冰雪附著速率降至無擾流板車型的54%,同時自身承受24%的附著量;而高200 mm的擾流板可使轉向架冰雪附著速率降至無擾流板車型的70%,同時自身承受9%的附著量。并且擾流板上的附著量不會結冰,會隨著氣流流向地面。

    3.3計算結果分析

    擾流板高度是影響附著率的一個關鍵因素,高度較大的后流線型擾流板能有效降低轉向架動壓,進而可能減少冰雪附著。圖5的結果驗證了這一結論,隨著擾流板高度的增加,后轉向架上的附著速率出現(xiàn)了明顯的下降,這是由于擾流板的增高加強了其遮蔽效應。由于局部來流速度更低,300 mm高度時后轉向架的冰雪附著甚至優(yōu)于較為理想的前轉向架;而若將后轉向架和擾流板綜合考慮,附著率降低效果在擾流板200 mm高時就已經穩(wěn)定,在100 mm以下則無明顯效果。需要指出,這里討論基于最不利工況,即默認參數(shù)+濕雪條件。

    圖5 前后轉向架附著速率與擾流板高度的關系

    4 結 論

    基于轉向架擾流板的風雪兩相流的初步研究,通過理論分析和運用STAR-CCM+模擬在風雪兩相流輸入

    情況下高寒動車組轉向架擾流板高度的防冰雪效果。結果表明:

    (1)拉格朗日兩相流模型計算高寒動車組冰雪附著問題是可靠而實用,能更直觀、有效的體現(xiàn)高寒動車組冰雪實際附著的情況,對改進擾流板更有效果;

    (2)在進行模擬分析時,發(fā)現(xiàn)雪強、雪顆粒附著特性、車速、擾流板高度對轉向架冰雪附著率均有影響;

    (3)通過模擬計算分析,流線型擾流板方案可使轉向架冰雪附著速率降至無擾流板車型的54%,同時自身承受24%的附著量;高200 mm的擾流板可使轉向架冰雪附著速率降至無擾流板車型的70%,同時自身承受9%的附著量。并且擾流板上的附著量不會結冰,而是隨著氣流流向地面。

    [1]別玉華. 高寒動車組車體保溫性能研究與優(yōu)化[D]. 青島:青島理工大學,2014.

    [2]北山茂,等. 臺車周りの空力特性改善に関する研究[C].日本機械學會第14回交通物流部門大會演講論文集(2005.12.7-9川崎).

    [3]李俊民,單永林,林鵬. 高速動車組轉向架防冰雪導流罩的空氣動力學性能分析[J].計算機輔助工程,2013,22(2):20-26,80.

    [4]李夢. 車身外部擾流器對汽車氣動特性的影響研究[D].吉林:吉林大學,2013.

    [5]盛強,刑玉明,何超. 基于CFD的機翼結冰過程分析[J]. 航空計算技術. 2009,39(2):37-40.

    [6]劉曉述. 風雪流運動對屋面積雪荷載分布影響數(shù)值模擬分析[D]. 沈陽:沈陽建筑大學,2013.

    [7]鐘錫昌,劉學宗. 計算流體力學[M]. 北京:科學出版社,1983.

    [8]黃學輝. 轉向架對高速列車氣動特性影響研究[D].成都:西南交通大學,2008.

    [9]Evert Andersson, HIGH-SPEED TRAIN OPERATION IN WINTER CLIMATEA STUDY ON WINTER RELATED PROBLEMS AND SOLUTIONS APPLIED IN SWEDEN[R], NORWAY AND FINLAND, 2006, ETH, BVF5 Winter R1.1 060703.

    [10]謝飛. 低馬赫數(shù)、低雷諾數(shù)翼型非定常流數(shù)值模擬[D].西安:西北工業(yè)大學,2003.

    [11]李明,李明高. STAR-CCM+與流場計算[M].北京:機械工業(yè)出版社,2011.

    Analysis on Two-phase Flowing Anti-ice/Snow of Low-temperature EMU Bogie

    HEDehua1,WANGGangyi2,CHENHouchang1,LIChaohui2

    (1Locomotive&Car Research Institute,China Academy of Railway Sciences, Beijing 100081, China; 2Electromechanical Engineering Institute,East China Jiaotong University, Nanchang 330013 Jiangxi, China)

    Facing the snow and ice accumulated problem of low-temperature EMU bogie in winter, we are using more intuitive technology of two-phase flowing to simulate and analyze the snow and ice accumulated situation of high-cold EMU bogie after designing of turbulence and evaluating the anti-ice/snow effect of bogie by the wall of snow particles incident on mass flux (IMF), which using commercial software STAR-CCM+ and Lagrangian particle model. The results show that: The snow and ice accumulated bogie are affected by strength of snowing, the attached characteristics of snow particles, speed of high-cold EMU and the height of spoiler; The improved streamlined spoilers(long 2000, high 300, wide 330) can cut the rate of snow and ice accumulated bogie by 54% than the train of non-spoilers while the spoilers bear 24% of attached volume; The spoilers(long 2000, high 200, wide 330) can cut the rate of snow and ice accumulated bogie by 70% than the train of non-spoilers while the spoilers bear 9% of attached volume. The volume which is attached on the spoiler will not freeze and flow to the surface with air.

    high-cold EMU; bogie; spoiler; two-phase flowing

    1008-7842 (2016) 04-0038-05

    ??)男,助理研究員(

    2016-03-08)

    U271.91

    Adoi:10.3969/j.issn.1008-7842.2016.04.09

    *中國鐵路總公司科技項目(2013YJ001)

    猜你喜歡
    擾流擾流板轉向架
    給飛機剎車的擾流板
    汽車擾流板的結構和試驗設計
    20t軸重米軌轉向架裝用搖枕的研制
    基于SPS模式的轉向架軸箱裝配線仿真研究
    淺談飛機擾流板靈活性試驗驗證方法
    工程與試驗(2019年4期)2019-03-27 03:33:06
    超聲流量計探頭擾流影響的流聲耦合研究
    測控技術(2018年5期)2018-12-09 09:04:32
    基于結構應力的轉向架構架焊縫疲勞強度研究
    連續(xù)擾流促進日光溫室番茄生長
    農家顧問(2016年11期)2017-01-06 12:58:52
    209 P型轉向架軸箱彈簧斷裂的原因及預防
    管內擾流元件的強化傳熱原理與性能指標研究進展
    化工進展(2014年1期)2014-08-08 09:52:04
    亚洲精品一卡2卡三卡4卡5卡| 国产精品野战在线观看| 日本撒尿小便嘘嘘汇集6| 国产人妻一区二区三区在| 别揉我奶头~嗯~啊~动态视频| 真人一进一出gif抽搐免费| 亚洲人成伊人成综合网2020| 国产综合懂色| 天堂影院成人在线观看| 看片在线看免费视频| 国产激情偷乱视频一区二区| 在线观看舔阴道视频| 精品福利观看| 九九在线视频观看精品| 欧美日韩黄片免| 欧美又色又爽又黄视频| eeuss影院久久| 黄片小视频在线播放| 九九热线精品视视频播放| 床上黄色一级片| 一个人免费在线观看电影| 中国美女看黄片| 日韩有码中文字幕| 99久久精品国产亚洲精品| 亚洲av电影不卡..在线观看| 亚洲人成网站高清观看| 亚洲午夜理论影院| 少妇熟女aⅴ在线视频| 午夜两性在线视频| 99国产极品粉嫩在线观看| 美女免费视频网站| 99热只有精品国产| 亚洲性夜色夜夜综合| 嫩草影院新地址| 嫩草影院入口| 久99久视频精品免费| 一级av片app| 淫妇啪啪啪对白视频| 国产精品人妻久久久久久| 激情在线观看视频在线高清| 身体一侧抽搐| 少妇被粗大猛烈的视频| 88av欧美| 免费看美女性在线毛片视频| 国产成人影院久久av| 三级毛片av免费| 国产精品久久久久久久久免 | 亚洲久久久久久中文字幕| 午夜精品一区二区三区免费看| 国产伦一二天堂av在线观看| 欧美区成人在线视频| 国产三级中文精品| 久久这里只有精品中国| 床上黄色一级片| 色播亚洲综合网| 日韩亚洲欧美综合| 色视频www国产| 国产aⅴ精品一区二区三区波| 日韩大尺度精品在线看网址| 欧美日韩瑟瑟在线播放| 午夜两性在线视频| 好男人电影高清在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产欧美日韩精品亚洲av| 精品免费久久久久久久清纯| 在现免费观看毛片| 国产单亲对白刺激| 怎么达到女性高潮| 国产精品国产高清国产av| 一进一出抽搐gif免费好疼| 亚洲精品日韩av片在线观看| 男女床上黄色一级片免费看| 制服丝袜大香蕉在线| 99在线视频只有这里精品首页| 欧美成人免费av一区二区三区| 午夜久久久久精精品| 国产亚洲精品av在线| 99视频精品全部免费 在线| 亚洲真实伦在线观看| 宅男免费午夜| 性色av乱码一区二区三区2| 18禁裸乳无遮挡免费网站照片| 免费在线观看成人毛片| 男女那种视频在线观看| 久久久久久久久久成人| 色哟哟哟哟哟哟| 久久精品夜夜夜夜夜久久蜜豆| 丰满乱子伦码专区| 99riav亚洲国产免费| 欧美3d第一页| 99热这里只有是精品50| 少妇熟女aⅴ在线视频| 又粗又爽又猛毛片免费看| 久久这里只有精品中国| 深爱激情五月婷婷| 亚洲五月婷婷丁香| 久久草成人影院| 亚洲中文日韩欧美视频| 亚洲专区国产一区二区| 亚洲成人精品中文字幕电影| 99热只有精品国产| 欧美激情在线99| 少妇人妻精品综合一区二区 | 久久久久精品国产欧美久久久| 国产精品人妻久久久久久| 欧美三级亚洲精品| 国产一区二区在线观看日韩| 欧美一区二区国产精品久久精品| 国产成人影院久久av| 国产v大片淫在线免费观看| 欧美潮喷喷水| 女同久久另类99精品国产91| 精品久久久久久,| 国产黄a三级三级三级人| 精品久久久久久久末码| 久久久国产成人精品二区| 亚洲内射少妇av| 最新中文字幕久久久久| 女同久久另类99精品国产91| 一级a爱片免费观看的视频| 久久人人爽人人爽人人片va | 在线播放国产精品三级| 国产午夜精品论理片| 91久久精品电影网| 一进一出好大好爽视频| 搡老妇女老女人老熟妇| 国产成人av教育| av在线老鸭窝| 国产精品一区二区三区四区免费观看 | 国产午夜福利久久久久久| 日本与韩国留学比较| 亚洲熟妇中文字幕五十中出| 乱人视频在线观看| 国产黄色小视频在线观看| 51国产日韩欧美| 国产高清三级在线| 亚洲经典国产精华液单 | 日本撒尿小便嘘嘘汇集6| 中文字幕久久专区| 午夜两性在线视频| 一区福利在线观看| 免费观看人在逋| 欧美色欧美亚洲另类二区| 成人美女网站在线观看视频| 中文亚洲av片在线观看爽| 少妇丰满av| 深爱激情五月婷婷| 热99re8久久精品国产| 日韩大尺度精品在线看网址| 婷婷亚洲欧美| 色在线成人网| 国产一区二区亚洲精品在线观看| 青草久久国产| 欧美国产日韩亚洲一区| 日韩精品中文字幕看吧| 日韩欧美精品v在线| 国产精品爽爽va在线观看网站| 别揉我奶头~嗯~啊~动态视频| 久久人妻av系列| 在线观看66精品国产| av欧美777| 搡老妇女老女人老熟妇| 成年免费大片在线观看| 亚洲中文日韩欧美视频| 国产精品综合久久久久久久免费| www.色视频.com| 丝袜美腿在线中文| 午夜福利成人在线免费观看| 好看av亚洲va欧美ⅴa在| av天堂在线播放| 成年女人永久免费观看视频| 亚洲成av人片免费观看| 精品人妻一区二区三区麻豆 | 欧美在线一区亚洲| 成年女人毛片免费观看观看9| 搡老熟女国产l中国老女人| 免费在线观看影片大全网站| 深夜精品福利| 男插女下体视频免费在线播放| 99久久成人亚洲精品观看| 伦理电影大哥的女人| av国产免费在线观看| 一个人免费在线观看电影| 麻豆av噜噜一区二区三区| 在线播放无遮挡| 神马国产精品三级电影在线观看| 女生性感内裤真人,穿戴方法视频| 国产成人aa在线观看| 看免费av毛片| 动漫黄色视频在线观看| 91麻豆精品激情在线观看国产| 欧美成人a在线观看| 久久草成人影院| 国产精品国产高清国产av| 欧美在线黄色| 亚洲欧美日韩卡通动漫| 精品午夜福利视频在线观看一区| 日韩国内少妇激情av| 亚洲乱码一区二区免费版| 噜噜噜噜噜久久久久久91| 美女 人体艺术 gogo| 日日摸夜夜添夜夜添小说| 日韩欧美 国产精品| 十八禁网站免费在线| 禁无遮挡网站| 中文字幕久久专区| 老熟妇仑乱视频hdxx| 国产av不卡久久| 99在线视频只有这里精品首页| 国产精品电影一区二区三区| 免费在线观看影片大全网站| 哪里可以看免费的av片| 国产真实乱freesex| 99久久精品一区二区三区| 欧美一区二区精品小视频在线| 国产精品三级大全| 51国产日韩欧美| netflix在线观看网站| 日本与韩国留学比较| 日本a在线网址| 国产精品亚洲av一区麻豆| 好看av亚洲va欧美ⅴa在| 村上凉子中文字幕在线| 一夜夜www| 国产精品久久电影中文字幕| 老熟妇乱子伦视频在线观看| 国产主播在线观看一区二区| 中亚洲国语对白在线视频| 级片在线观看| 国产亚洲精品久久久久久毛片| 国产欧美日韩一区二区三| 久久久久九九精品影院| 亚洲av成人精品一区久久| 悠悠久久av| 丰满乱子伦码专区| 亚洲av熟女| 国产日本99.免费观看| 一本综合久久免费| 亚洲国产精品sss在线观看| 久久久久久久精品吃奶| 亚洲第一电影网av| 99视频精品全部免费 在线| 亚洲自偷自拍三级| 18禁黄网站禁片午夜丰满| 18+在线观看网站| 日本与韩国留学比较| 一区二区三区四区激情视频 | 在线天堂最新版资源| www.999成人在线观看| 亚洲成人久久爱视频| 99久久99久久久精品蜜桃| 亚洲精品一区av在线观看| 欧美成人性av电影在线观看| 亚洲人成网站在线播| 又紧又爽又黄一区二区| 国产91精品成人一区二区三区| 午夜精品在线福利| 国产精品人妻久久久久久| 国产午夜精品久久久久久一区二区三区 | 精品无人区乱码1区二区| 91在线观看av| 美女免费视频网站| 午夜福利在线观看吧| 亚洲在线自拍视频| 精品人妻偷拍中文字幕| 我要搜黄色片| 欧美另类亚洲清纯唯美| 在线观看66精品国产| 亚洲av成人av| 在线免费观看不下载黄p国产 | 国产亚洲精品久久久久久毛片| 久久精品影院6| 国产精品伦人一区二区| 亚洲色图av天堂| 日韩欧美三级三区| 精品人妻熟女av久视频| 国产亚洲精品av在线| 在线观看一区二区三区| 国产成人影院久久av| 午夜视频国产福利| 亚洲美女黄片视频| 亚洲国产精品成人综合色| 一本精品99久久精品77| 麻豆久久精品国产亚洲av| 午夜福利免费观看在线| 国产v大片淫在线免费观看| 日韩中文字幕欧美一区二区| 久久人人爽人人爽人人片va | 99久久九九国产精品国产免费| 色吧在线观看| 国产麻豆成人av免费视频| 欧美乱色亚洲激情| 久久99热这里只有精品18| 日韩亚洲欧美综合| 韩国av一区二区三区四区| 色噜噜av男人的天堂激情| 国产成人影院久久av| 国产伦精品一区二区三区四那| 3wmmmm亚洲av在线观看| 男插女下体视频免费在线播放| 丁香六月欧美| 精华霜和精华液先用哪个| 一本久久中文字幕| 日韩欧美在线乱码| 男女下面进入的视频免费午夜| 给我免费播放毛片高清在线观看| 成人欧美大片| 日韩中字成人| 18禁裸乳无遮挡免费网站照片| 看片在线看免费视频| 欧美在线黄色| 久久久国产成人免费| 亚洲成av人片在线播放无| 少妇被粗大猛烈的视频| 久久午夜亚洲精品久久| 国产午夜精品论理片| 国产主播在线观看一区二区| 精品一区二区三区视频在线| 免费看光身美女| 国产真实伦视频高清在线观看 | 别揉我奶头~嗯~啊~动态视频| 欧美+日韩+精品| 久久久久久久久大av| 美女 人体艺术 gogo| 国产白丝娇喘喷水9色精品| 精品午夜福利在线看| 亚洲一区高清亚洲精品| .国产精品久久| 久久精品综合一区二区三区| 色av中文字幕| 亚洲一区二区三区色噜噜| 国内毛片毛片毛片毛片毛片| 午夜久久久久精精品| 99国产极品粉嫩在线观看| 午夜福利在线观看吧| 99久久99久久久精品蜜桃| 亚洲综合色惰| 在线免费观看不下载黄p国产 | 日韩免费av在线播放| 69av精品久久久久久| 国产av麻豆久久久久久久| 最新中文字幕久久久久| 精品熟女少妇八av免费久了| 成人高潮视频无遮挡免费网站| 亚洲精品色激情综合| 两人在一起打扑克的视频| 在线观看免费视频日本深夜| 亚洲性夜色夜夜综合| 乱码一卡2卡4卡精品| 免费看a级黄色片| 一本久久中文字幕| 成人永久免费在线观看视频| 亚洲激情在线av| 免费在线观看亚洲国产| 国产黄片美女视频| 精品久久国产蜜桃| 欧美高清性xxxxhd video| 亚洲 国产 在线| 18美女黄网站色大片免费观看| 麻豆一二三区av精品| 99久久精品热视频| 99久久精品国产亚洲精品| 国产精品久久久久久久电影| 超碰av人人做人人爽久久| 九色国产91popny在线| 人人妻人人澡欧美一区二区| 国产精华一区二区三区| 国产免费男女视频| 国产国拍精品亚洲av在线观看| 少妇人妻精品综合一区二区 | 熟妇人妻久久中文字幕3abv| www.999成人在线观看| 国产v大片淫在线免费观看| 国产高潮美女av| 亚洲电影在线观看av| a级毛片免费高清观看在线播放| 校园春色视频在线观看| 国产爱豆传媒在线观看| 色精品久久人妻99蜜桃| 国产精品不卡视频一区二区 | 男女那种视频在线观看| 搡老岳熟女国产| 国产美女午夜福利| av天堂中文字幕网| 国产精品自产拍在线观看55亚洲| 国产 一区 欧美 日韩| 无遮挡黄片免费观看| 亚洲人与动物交配视频| av天堂中文字幕网| 一级毛片久久久久久久久女| 啦啦啦观看免费观看视频高清| 国产主播在线观看一区二区| 久久99热6这里只有精品| 午夜福利在线观看免费完整高清在 | 嫩草影院新地址| 波多野结衣巨乳人妻| 一进一出好大好爽视频| 国产精品久久久久久久久免 | 少妇丰满av| 日韩欧美精品免费久久 | 精华霜和精华液先用哪个| 亚洲五月婷婷丁香| 精品福利观看| 欧美性感艳星| 国产探花在线观看一区二区| 国产高清激情床上av| 免费在线观看成人毛片| 人妻夜夜爽99麻豆av| 欧美xxxx黑人xx丫x性爽| 亚洲国产日韩欧美精品在线观看| 99热只有精品国产| 亚洲avbb在线观看| 国产av在哪里看| 一级作爱视频免费观看| av福利片在线观看| 亚洲av一区综合| 亚洲av中文字字幕乱码综合| 欧美黄色片欧美黄色片| 免费黄网站久久成人精品 | 国产单亲对白刺激| 在线观看av片永久免费下载| 97热精品久久久久久| 内射极品少妇av片p| 欧美精品国产亚洲| 黄色丝袜av网址大全| 亚洲真实伦在线观看| 最近在线观看免费完整版| 亚洲色图av天堂| 亚洲av一区综合| a级一级毛片免费在线观看| 香蕉av资源在线| 国产v大片淫在线免费观看| 大型黄色视频在线免费观看| 美女cb高潮喷水在线观看| 99久久九九国产精品国产免费| 久久久精品欧美日韩精品| 免费人成在线观看视频色| avwww免费| 一级av片app| 久久精品人妻少妇| 亚洲成av人片免费观看| 婷婷精品国产亚洲av| 69人妻影院| 好看av亚洲va欧美ⅴa在| 国产一区二区激情短视频| 国产成人a区在线观看| 久久久成人免费电影| 中出人妻视频一区二区| 此物有八面人人有两片| 国产欧美日韩一区二区精品| 波野结衣二区三区在线| 日日干狠狠操夜夜爽| 国产亚洲精品av在线| 变态另类丝袜制服| 男女那种视频在线观看| 国产精品av视频在线免费观看| 欧美成人免费av一区二区三区| 午夜福利在线在线| 狂野欧美白嫩少妇大欣赏| 国产一级毛片七仙女欲春2| 性色avwww在线观看| 国产精品乱码一区二三区的特点| 亚洲无线在线观看| 久久久精品大字幕| 亚洲,欧美精品.| 真实男女啪啪啪动态图| 国内揄拍国产精品人妻在线| 国产精品久久久久久久久免 | 一个人免费在线观看电影| 欧美高清成人免费视频www| 久久久久久久午夜电影| 国产精品电影一区二区三区| www.999成人在线观看| 老女人水多毛片| 99国产极品粉嫩在线观看| 悠悠久久av| 久久九九热精品免费| 全区人妻精品视频| 久久天躁狠狠躁夜夜2o2o| 99久久九九国产精品国产免费| 男女之事视频高清在线观看| 国产亚洲精品av在线| 床上黄色一级片| 人妻丰满熟妇av一区二区三区| 97超视频在线观看视频| 亚洲精品456在线播放app | eeuss影院久久| 久久99热6这里只有精品| 久久精品人妻少妇| 成人av在线播放网站| 国产精品一区二区性色av| 十八禁人妻一区二区| 波野结衣二区三区在线| 欧美zozozo另类| 国产精品99久久久久久久久| 在线观看美女被高潮喷水网站 | 搡老岳熟女国产| 首页视频小说图片口味搜索| 精品人妻1区二区| 美女 人体艺术 gogo| 成人毛片a级毛片在线播放| 给我免费播放毛片高清在线观看| 女同久久另类99精品国产91| 亚洲精品影视一区二区三区av| 久久久久国产精品人妻aⅴ院| 嫩草影院新地址| 欧美绝顶高潮抽搐喷水| 天美传媒精品一区二区| 老司机深夜福利视频在线观看| 我要看日韩黄色一级片| 赤兔流量卡办理| 国产一区二区亚洲精品在线观看| 免费观看精品视频网站| 97超视频在线观看视频| 中文字幕免费在线视频6| 精华霜和精华液先用哪个| 国产单亲对白刺激| 亚洲国产精品sss在线观看| 男人狂女人下面高潮的视频| av在线观看视频网站免费| 中文字幕精品亚洲无线码一区| 此物有八面人人有两片| 一级黄片播放器| 亚洲精品在线观看二区| 99热这里只有是精品在线观看 | 久久香蕉精品热| www.www免费av| 身体一侧抽搐| 午夜精品久久久久久毛片777| 精品人妻偷拍中文字幕| 久久午夜福利片| 亚洲国产高清在线一区二区三| 欧美成狂野欧美在线观看| 一进一出抽搐gif免费好疼| 一进一出抽搐动态| 日韩欧美国产一区二区入口| xxxwww97欧美| 亚洲成人中文字幕在线播放| 免费黄网站久久成人精品 | 欧美另类亚洲清纯唯美| 丰满的人妻完整版| 淫妇啪啪啪对白视频| 有码 亚洲区| 91九色精品人成在线观看| 国产精品女同一区二区软件 | 九九在线视频观看精品| 毛片女人毛片| 成人毛片a级毛片在线播放| 国产中年淑女户外野战色| 日本五十路高清| 国产v大片淫在线免费观看| 亚洲精品成人久久久久久| 久久伊人香网站| 欧美黑人欧美精品刺激| 成人一区二区视频在线观看| 丰满乱子伦码专区| 国产伦一二天堂av在线观看| 两个人视频免费观看高清| 在线a可以看的网站| 国产精品嫩草影院av在线观看 | 国内精品一区二区在线观看| 性插视频无遮挡在线免费观看| 全区人妻精品视频| 简卡轻食公司| 2021天堂中文幕一二区在线观| 如何舔出高潮| 色吧在线观看| 国产av麻豆久久久久久久| 国内毛片毛片毛片毛片毛片| 在线十欧美十亚洲十日本专区| 一个人观看的视频www高清免费观看| 在线观看免费视频日本深夜| 精品一区二区三区视频在线| 精品人妻1区二区| 国产精品久久久久久久久免 | 日本一二三区视频观看| 欧美一区二区精品小视频在线| 亚洲成a人片在线一区二区| 婷婷色综合大香蕉| 一本精品99久久精品77| 国产三级黄色录像| 国产伦一二天堂av在线观看| 制服丝袜大香蕉在线| 黄色日韩在线| 亚洲最大成人手机在线| 男人舔奶头视频| 老司机福利观看| 在线免费观看的www视频| 窝窝影院91人妻| 亚洲,欧美精品.| 波多野结衣巨乳人妻| 狠狠狠狠99中文字幕| 熟女电影av网| 欧美xxxx性猛交bbbb| 日韩欧美国产一区二区入口| 老司机福利观看| 97超视频在线观看视频| 日韩欧美在线乱码| 亚洲aⅴ乱码一区二区在线播放| 国产视频一区二区在线看| 亚洲成人久久爱视频| 亚洲国产精品999在线| 中文在线观看免费www的网站| 有码 亚洲区| 欧美绝顶高潮抽搐喷水| www.www免费av| 日韩精品青青久久久久久| 亚洲av不卡在线观看| 久久久久久久久久黄片| 日韩欧美在线乱码| 一个人观看的视频www高清免费观看| 亚洲国产高清在线一区二区三| 中文字幕精品亚洲无线码一区| 最新中文字幕久久久久| 啦啦啦观看免费观看视频高清|