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

    大跨徑懸索橋施工貓道氣動力系數(shù)數(shù)值分析

    2015-03-25 07:12:52李勝利胡亞楠杜道傾王東煒
    關(guān)鍵詞:氣動力雷諾數(shù)攻角

    李勝利,胡亞楠,杜道傾,王東煒

    (鄭州大學(xué) 土木工程學(xué)院,河南 鄭州450001)

    0 引言

    貓道是大跨度柔性透風(fēng)索結(jié)構(gòu)且剛度小,對風(fēng)敏感.貓道在風(fēng)荷載作用下,若產(chǎn)生較大變位將會影響主纜的施工安全甚至成橋后主纜的受力[1-3].為了加快施工進度,降低費用,保證正常的通航要求,國內(nèi)外均開始使用不設(shè)抗風(fēng)纜的貓道,這就對該類貓道的抗風(fēng)安全提出了更高的要求[4-5].為了確保無抗風(fēng)纜貓道安全,有必要對其氣動力系數(shù)進行分析. 在過去相當(dāng)長的一段時間里,風(fēng)工程中常假定具有尖銳邊角斷面的實體(即鈍體)周圍流場不依賴于雷諾數(shù)的變化[6].實踐表明,在鈍體斷面的橋梁結(jié)構(gòu)流場中存在雷諾數(shù)效應(yīng)[7]隨著懸索橋跨度越來越大,貓道也向大跨、輕質(zhì)方向發(fā)展,雷諾數(shù)對貓道斷面氣動力系數(shù)的影響日漸突出,忽略其影響會給貓道抗風(fēng)設(shè)計帶來危險[8-10].因此,研究雷諾數(shù)對貓道斷面氣動力系數(shù)的影響顯得非常必要,隨著計算機技術(shù)和計算流體力學(xué)方法的迅速發(fā)展,CFD 數(shù)值模擬方法已經(jīng)成為研究貓道斷面氣動力系數(shù)的新手段,該方法不僅可以按照實際尺寸進行計算,而且費用少、速度快、資料完備,同時,還具有模擬真實和理想條件的能力[11].筆者以某大跨徑懸索橋貓道斷面氣動力系數(shù)為研究對象,利用計算流體力學(xué)軟件Fluent,通過采用3 種不同湍流模型,即標(biāo)準(zhǔn)k-ε 模型、RNG k-ε 模型和SST 模型,對貓道斷面的氣動力系數(shù)在不同風(fēng)速、不同攻角和不同雷諾數(shù)情況下進行CFD 方法識別,與節(jié)段模型風(fēng)洞試驗的試驗值進行對比分析.

    1 貓道氣動力系數(shù)模擬參數(shù)確定

    1.1 工程概況

    圖1 ~3 為某大跨徑懸索橋貓道的平面圖、橫斷面圖和側(cè)面圖. 西南交通大學(xué)風(fēng)工程試驗研究中心對該貓道進行了節(jié)段模型靜力試驗,試驗風(fēng)速為:v=10,20,30 m/s;試驗攻角為:α = -5° ~+5°,節(jié)段模型采用1∶6 的幾何縮尺比,其模型長L=2.100 m,寬B=0.692 7 m,高H=0.241 7 m,長寬比L/B=3.03[12].本研究采用數(shù)值模擬的參數(shù)設(shè)置方法,CFD 模型所采用的尺寸、比例、風(fēng)速、攻角和風(fēng)洞試驗相同,計算采用風(fēng)軸坐標(biāo)系,其定義如圖4 所示.

    1.2 計算參數(shù)的設(shè)置

    本研究采用Gambit 軟件建立數(shù)值模型.由于風(fēng)洞試驗一般要求測出氣動力系數(shù)在-5° ~5°區(qū)間的變化曲線,因此通過數(shù)值模擬得到貓道斷面風(fēng)攻角在-5° ~5°范圍內(nèi)的氣動力系數(shù)是合理的、可行的.貓道斷面用橫向單排直徑為4 mm 和豎向單排直徑為2 mm 的圓來模擬,模擬透風(fēng)率分別為84.09%和85.6%,貓道實際透風(fēng)率分別為84.2%和85.5%,模擬值和實際值差別很?。P涂s尺比取1/6,計算域為24 m ×8 m 的矩形區(qū)域(如圖5 所示),入口距迎風(fēng)面的距離取8 m,計算域上下邊界距模型的距離分別為3.75 m 和4 m,背風(fēng)面距出口的距離取15.3 m,這樣可以盡量避免物體后部卷起的分離漩渦打到外邊界又反射回來,同時也使外邊界附近的流場參數(shù)分布較好地與所假定的邊界條件相容,求解可以達到很好的收斂性[13]. 邊界條件:迎風(fēng)面邊界條件定義為速度進口邊界條件,空氣密度取1.225 kg/m3;背風(fēng)面邊界條件定義為壓力出口邊界條件,相對壓力值取為零;上下邊界均為無滑移的墻面;貓道斷面外壁定義為無滑移的墻面(如圖5 所示). 由于貓道斷面沿縱向基本不變,因此貓道計算模型可近似按二維模型處理,橫風(fēng)向空氣的流動處理為二維定常不可壓,在計算分析中不考慮溫度變化和能量交換,靜力計算按定常來處理,流場計算方法選擇SIMPLE 算法,離散格式選擇二階迎風(fēng)格式.同時采用非結(jié)構(gòu)化網(wǎng)格對貓道斷面計算域進行網(wǎng)格劃分,并對貓道斷面壁面處采取加密網(wǎng)格,使數(shù)值模擬更接近于實際情況,迭代殘差值控制為10-6.因為網(wǎng)格劃分的質(zhì)量好壞直接影響模擬結(jié)果的精度,所以在進行計算分析之前,要對網(wǎng)格劃分的質(zhì)量進行評估,對于網(wǎng)格劃分質(zhì)量較差的要進行優(yōu)化[14].筆者在劃分網(wǎng)格時通過數(shù)次優(yōu)化調(diào)節(jié),最終得到較合理的網(wǎng)格劃分形式,如圖6 所示.

    圖5 貓道斷面計算域Fig.5 Computational domain of catwalk cross section

    圖6 貓道斷面CFD 計算網(wǎng)格Fig.6 The catwalk cross section of CFD computational grid

    1.3 湍流模型的選擇

    目前,在實際工程中被廣泛采用的鈍體的繞流湍流模式是渦黏性模式即兩方程模型里的標(biāo)準(zhǔn)k-ε 模型、RNG k-ε 模型和SST 模型.

    標(biāo)準(zhǔn)k-ε 模型是應(yīng)用最廣泛同時也是經(jīng)驗積累最多的湍流模型.其方程形式簡單,計算量比較小,也能很好地模擬實際工程的湍流運動,并能取得較好的模擬結(jié)果,因此是模擬湍流流動的首選模型[15].由于標(biāo)準(zhǔn)k-ε 模型假定動黏性系數(shù)是各向同性的標(biāo)量并且在彎曲流線的情況下,湍流流動是各向異性的,從而在模擬強旋流、彎曲流線或彎曲壁面時出現(xiàn)一定的失真. 為了彌補標(biāo)準(zhǔn)k-ε模型的缺陷,許多研究者對標(biāo)準(zhǔn)k-ε 模型進行了修正,RNG k-ε 模型就是其中的一種修正方案.RNG k-ε 模型是標(biāo)準(zhǔn)k-ε 模型的改進方案,與標(biāo)準(zhǔn)k-ε 模型相比,RNG k-ε 模型修正了湍流的黏度和在ε 方程中增加了反映流體的時均應(yīng)變率項.為了彌補標(biāo)準(zhǔn)k-ε 模型對湍流模擬的缺陷,Mente 通過混合函數(shù)將二者結(jié)合起來而產(chǎn)生的SST 模型,既可以在靠近壁面的區(qū)域采用k-ω 湍流模型,又可以在邊界層邊緣以及自由剪切層采用k-ω 形式的k-ε 模型并且過渡區(qū)域可以通過混合函數(shù)來過渡[16],所以SST 模型集結(jié)了標(biāo)準(zhǔn)k-ε湍流模型和k-ω 湍流模型的優(yōu)點[17].為了比較選擇合適的貓道數(shù)值模擬湍流模型,筆者進行了3種湍流模型的計算.

    1.4 氣動力系數(shù)的定義

    對于阻力系數(shù)來說,由圖7 可見,在-5° ~0°攻角范圍內(nèi),SST 湍流模型的模擬值和試驗值非常接近,但試驗值比模擬值略大;在0° ~5°攻角范圍內(nèi),標(biāo)準(zhǔn)k-ε 湍流模型的模擬值均比試驗值大;-5° ~5°攻角范圍內(nèi),SST 湍流模型的模擬值變化趨勢總體較其他計算模型更符合試驗值變化趨勢.可見,不同風(fēng)攻角時,理論模擬值和試驗值差別大小不同,但整體上SST 湍流模型的模擬值更接近試驗值.

    圖7 貓道斷面阻力系數(shù)與風(fēng)洞試驗值對比Fig.7 Compared drag coefficients of the catwalk cross section with results of wind tunnel test

    作用在貓道上的靜風(fēng)荷載可分為橫向風(fēng)載(阻力)D,豎向風(fēng)載(升力)L. 在靜風(fēng)作用下,結(jié)構(gòu)的姿態(tài)發(fā)生改變,靜風(fēng)與貓道截面的相對攻角隨之改變,氣動力系數(shù)是隨攻角變化而改變的,同時平均風(fēng)速也隨著貓道的變形而改變,因此引入有效風(fēng)攻角和有效風(fēng)速的概念,將靜風(fēng)荷載表示為有效風(fēng)速、氣動力系數(shù)和有效風(fēng)攻角的函數(shù),具體氣動力系數(shù)的定義如下

    對于升力系數(shù)來說,由圖8 可見,在所有風(fēng)攻角范圍內(nèi),3 種湍流模型的模擬值差別微?。?在-5° ~-3°攻角范圍內(nèi),標(biāo)準(zhǔn)k-ε 湍流模型和試驗值有較大差別;在-5° ~-3°攻角范圍內(nèi),試驗值均比RNG k-ε 和SST 湍流模型模擬值大;在-3° ~5°攻角范圍內(nèi),試驗值與模擬值變化趨勢基本相同;整體來說,RNG k-ε 和SST 湍流模型的模擬值與風(fēng)洞試驗值更為接近.

    圖8 貓道斷面升力系數(shù)與風(fēng)洞試驗值對比Fig.8 Compared life coefficients of the catwalk cross section with results of wind tunnel test

    式中:Cd、Cl分別為阻力系數(shù)、升力系數(shù);D,L 分別為風(fēng)軸坐標(biāo)系下貓道所受到的氣動阻力、升力;ρ為空氣密度,取1.225 kg/m3;U 為離斷面足夠遠處來流的平均風(fēng)速;H 表示節(jié)段模型的高度;l 表示節(jié)段模型的長度;B 為節(jié)段模型寬度.

    2 數(shù)值模擬與試驗結(jié)果對比分析

    3 種不同湍流模型的數(shù)值計算結(jié)果與試驗結(jié)果如圖7 和8 所示.

    由上述分析可知,在-5° ~5°風(fēng)攻角范圍內(nèi),CFD 計算的阻力和升力系數(shù)與試驗值變化趨勢基本一致,并且在阻力系數(shù)曲線吻合程度上,SST模型曲線的變化趨勢與試驗結(jié)果趨勢基本一致,能較好地反映試驗值的變化規(guī)律. 雖然3 種模型計算結(jié)果與試驗值均有一定差距,但是考慮到絕對誤差值均較小,因此阻力和升力系數(shù)的模擬值和試驗值的差別是可以接受的,同時也說明本研究所用貓道透風(fēng)率的二維模擬方法是可行的.

    3 貓道氣動力系數(shù)的雷諾數(shù)效應(yīng)

    由上節(jié)的計算結(jié)果可知,SST 模型為貓道斷面氣動力系數(shù)CFD 方法識別的最合理的模型.因此筆者采用SST 模型,研究雷諾數(shù)Re 變化對不同風(fēng)攻角下貓道斷面氣動力系數(shù)的影響,雷諾數(shù)的設(shè)置方法如表1 所示. CFD 模擬計算不同雷諾數(shù)對貓道氣動力系數(shù)的影響結(jié)果如圖9 和10 所示.

    表1 貓道試驗斷面風(fēng)速和雷諾數(shù)表Tab.1 The table of catwalk test cross section wind speed and reynolds numbers

    由圖9 可見,當(dāng)雷諾數(shù)為104~106數(shù)量級時,0°攻角附近,不同雷諾數(shù)阻力系數(shù)相差很小;在-5° ~0°攻角范圍內(nèi),雷諾數(shù)較小時,阻力系數(shù)較大,隨著雷諾數(shù)的增大,阻力系數(shù)減小,雷諾數(shù)達到105數(shù)量級時,阻力系數(shù)非常接近;在0° ~5°攻角范圍內(nèi),隨著雷諾數(shù)的增大,阻力系數(shù)減小,且雷諾數(shù)較大時,阻力系數(shù)相差微?。w來看,隨雷諾數(shù)的增大,阻力系數(shù)CFD 模擬的相對誤差最大為15.4%,模擬值和試驗值[17]相對誤差最大為13.8%,誤差不大.

    由圖10 可見,當(dāng)雷諾數(shù)為104~106數(shù)量級時,在-5° ~0°攻角范圍內(nèi),隨著雷諾數(shù)的增大,升力系數(shù)增大,當(dāng)雷諾數(shù)較大時,升力系數(shù)幾乎相同;在0° ~5°攻角范圍內(nèi),雷諾數(shù)的大小和升力系數(shù)沒有明顯變化規(guī)律,但不同雷諾數(shù)之間升力系數(shù)差別很?。w來看,隨著雷諾數(shù)的增大,個別升力系數(shù)CFD 模擬的相對誤差最大為41.3%,模擬值和試驗值相對誤差最大為41.6%,但從圖上總體來看,模擬值與試驗值變化趨勢基本吻合,模擬結(jié)果是可以接受的.

    另外,由圖9 和圖10 可知,當(dāng)雷諾數(shù)為103數(shù)量級時,阻力系數(shù)和升力系數(shù)均和風(fēng)洞試驗數(shù)據(jù)有相當(dāng)大差距.阻力系數(shù)普遍偏大,相對誤差都在20%以上;個別升力系數(shù)差值甚至達到4 倍.阻力系數(shù)與升力系數(shù)變化無規(guī)律可循,對其模擬數(shù)據(jù)分析意義不大.

    由上述分析可知,當(dāng)雷諾數(shù)大小數(shù)量級與風(fēng)洞試驗一致,或相差不大時,雷諾數(shù)的變化對阻力系數(shù)影響較小,但對個別升力系數(shù)影響較大,雷諾數(shù)越接近風(fēng)洞試驗,雷諾數(shù)的數(shù)量級誤差越小,因此在進行貓道的數(shù)值模擬時,選擇與風(fēng)洞試驗一致的雷諾數(shù)很有必要.

    4 結(jié)論

    (1)當(dāng)數(shù)值模擬時設(shè)置的雷諾數(shù)大小數(shù)量級與試驗工況本身實際的雷諾數(shù)一致,或與試驗工況本身雷諾數(shù)本身相差不大時,雷諾數(shù)的變化對阻力系數(shù)相對誤差影響較小,但對個別升力系數(shù)影響相對較大,所以在進行貓道的數(shù)值模擬時,選擇接近與試驗工況本身實際的雷諾數(shù)很有必要;

    (2)在-5° ~5°風(fēng)攻角范圍內(nèi),SST 湍流模型為貓道斷面阻力和升力系數(shù)CFD 方法識別的較合理的模型.

    [1] 楊學(xué)祥,張明卓,朱艷,等.大跨度懸索橋施工貓道若干問題綜述[J].世界橋梁,2007(2):61 -64.

    [2] 吳勝東,馮兆祥,蔣波.特大跨徑懸索橋上部結(jié)構(gòu)施工關(guān)鍵技術(shù)研究[J]. 土木工程學(xué)報,2007,40(4):32 -37.

    [3] 賈寧,劉健新,劉萬鋒.懸索橋施工貓道靜風(fēng)失穩(wěn)機理分析[J].公路交通科技,2008,25(3):99 -102.

    [4] 李勝利. 大跨徑懸索橋施工期暫態(tài)結(jié)構(gòu)抗風(fēng)性能及控制[D]. 哈爾濱:哈爾濱工業(yè)大學(xué)土木工程學(xué)院,2010.

    [5] 李勝利,胡亞楠,王東煒.臺風(fēng)-波浪耦合作用下鋼管樁施工平臺力學(xué)行為分析[J].鄭州大學(xué)學(xué)報:工學(xué)版,2014,35(6):82 -85.

    [6] 何晗欣,劉健新. 大跨度懸索橋抗風(fēng)措施及其優(yōu)化研究[J]. 公 路 交 通 科 技,2011,28 (3):93 -97,102.

    [7] 趙凱,李永樂,歐陽韋,等.大跨度斜拉橋的靜風(fēng)失穩(wěn)特性[J].公路交通科技,2011,28(8):67 -72.

    [8] SELVAM R P,BOSCH H,JOSHI R. Comparison of 2D and 3D CFD Modeling of Bridge Aerodynamics[C]∥The 5th International Symposium On Computational Wind Engineering (CWE2010). Chapel Hill:North Carolina,2010:23 -27.

    [9] 李加武,崔欣,張宏杰,等. 粗糙度對雷諾數(shù)效應(yīng)的影響[J].長安大學(xué)學(xué)報:自然科學(xué)版,2009,29(2):56 - 59.

    [10]李薇,胡兆同,李加武. CFD 方法研究橋梁斷面三分力系數(shù)的雷諾數(shù)效應(yīng)[J]. 長安大學(xué)學(xué)報:自然科學(xué)版,2010,30(6):44 -49.

    [11]譚紅霞,陳振清. CFD 在橋梁斷面靜力三分力系數(shù)計算中的應(yīng)用[J]. 工程力學(xué),2009 ,26(11 ):68 -72.

    [12]廖海黎.舟山西堠門大橋貓道抗風(fēng)穩(wěn)定性風(fēng)洞模型試驗[R].成都:西南交通大學(xué)風(fēng)工程試驗研究中心,2005.

    [13]瞿偉廉,劉琳娜.基于CFD 的橋梁三分力系數(shù)識別的數(shù)值研究[J]. 武漢理工大學(xué)學(xué)報,2007,29(7):85 -88.

    [14]林斌,孫曉穎,武岳. 大慶石油學(xué)院體育館屋面風(fēng)荷載的風(fēng)洞試驗及CFD 數(shù)值模擬[J]. 沈陽建筑大學(xué)學(xué)報:自然科學(xué)版,2006,22(3):357 -361.

    [15] WATANABE S,F(xiàn)UMOTO K. Aerodynamic study of slotted box girder using computational fluid dynamics[J]. Wind Engineering and IndustrialAerodynamics,2008,96:1885 -1894.

    [16] SIRIBOONLUCKUL N,JUNTASARO V. Turbulence modeling for wall-bounded particle-laden flow with separation[J]. International Communications in Heat and Mass Transfer,2007,34(3):331 -338.

    [17] BAGLIETTO E,NINOKATA H. A turbulence model study for simulating flow inside tight lattice rod bundles[J]. Nuclear Engineering and Design,2005,235(7):773 -784.

    猜你喜歡
    氣動力雷諾數(shù)攻角
    飛行載荷外部氣動力的二次規(guī)劃等效映射方法
    風(fēng)標(biāo)式攻角傳感器在超聲速飛行運載火箭中的應(yīng)用研究
    大攻角狀態(tài)壓氣機分離流及葉片動力響應(yīng)特性
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    側(cè)風(fēng)對拍動翅氣動力的影響
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計研究
    民機高速風(fēng)洞試驗的阻力雷諾數(shù)效應(yīng)修正
    附加攻角效應(yīng)對顫振穩(wěn)定性能影響
    振動與沖擊(2015年2期)2015-05-16 05:37:34
    民用飛機攻角傳感器安裝定位研究
    欧美av亚洲av综合av国产av| 99在线视频只有这里精品首页| 91麻豆精品激情在线观看国产| 又紧又爽又黄一区二区| 亚洲国产中文字幕在线视频| 久久精品成人免费网站| avwww免费| 亚洲午夜精品一区,二区,三区| 在线播放国产精品三级| 免费搜索国产男女视频| 亚洲片人在线观看| 精品一区二区三区视频在线观看免费| √禁漫天堂资源中文www| tocl精华| 欧美性猛交黑人性爽| 国产三级在线视频| 波多野结衣巨乳人妻| 欧美丝袜亚洲另类 | av电影中文网址| 欧美丝袜亚洲另类 | 2021天堂中文幕一二区在线观 | 亚洲激情在线av| 国产一区二区三区视频了| 欧美乱码精品一区二区三区| 欧美色欧美亚洲另类二区| 欧美性猛交╳xxx乱大交人| 精品久久久久久久毛片微露脸| 99在线人妻在线中文字幕| 成人特级黄色片久久久久久久| 成人三级做爰电影| 久久久久久久午夜电影| 欧美三级亚洲精品| 国产人伦9x9x在线观看| 麻豆国产av国片精品| 19禁男女啪啪无遮挡网站| 精品人妻1区二区| 亚洲专区字幕在线| 国产v大片淫在线免费观看| 久久中文字幕一级| 国产精品香港三级国产av潘金莲| 婷婷亚洲欧美| 久久青草综合色| 很黄的视频免费| 亚洲黑人精品在线| 免费av毛片视频| 国产99白浆流出| 免费在线观看日本一区| 国产成人欧美在线观看| 日韩高清综合在线| 亚洲狠狠婷婷综合久久图片| 他把我摸到了高潮在线观看| 日韩免费av在线播放| 日韩中文字幕欧美一区二区| 久久精品国产亚洲av香蕉五月| 国产乱人伦免费视频| 国产精品野战在线观看| 亚洲国产精品久久男人天堂| 99精品久久久久人妻精品| 亚洲av熟女| 变态另类丝袜制服| 一a级毛片在线观看| 亚洲成人久久性| 欧美精品亚洲一区二区| 看黄色毛片网站| 亚洲国产精品久久男人天堂| 在线观看免费午夜福利视频| 老鸭窝网址在线观看| 亚洲精品粉嫩美女一区| 国内揄拍国产精品人妻在线 | 午夜福利18| 亚洲一区二区三区色噜噜| 日韩精品中文字幕看吧| 黑人欧美特级aaaaaa片| 成年女人毛片免费观看观看9| 国产真实乱freesex| 亚洲精品久久国产高清桃花| videosex国产| 免费看a级黄色片| 精品国产亚洲在线| 欧美日韩瑟瑟在线播放| 国产精品一区二区免费欧美| 精品久久蜜臀av无| 国产精品 国内视频| 香蕉丝袜av| 久久久久久久久久黄片| 午夜福利18| 国产一区二区三区视频了| 久热爱精品视频在线9| 国产亚洲精品av在线| 婷婷丁香在线五月| 麻豆国产av国片精品| 中文字幕另类日韩欧美亚洲嫩草| 日韩三级视频一区二区三区| 久久中文看片网| 国产亚洲精品第一综合不卡| 国产一区二区三区视频了| 又黄又爽又免费观看的视频| 久久久久久久久中文| 免费在线观看成人毛片| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美中文日本在线观看视频| 一级a爱片免费观看的视频| 亚洲欧美精品综合一区二区三区| 国产激情偷乱视频一区二区| 国产片内射在线| 亚洲第一欧美日韩一区二区三区| tocl精华| 国产高清有码在线观看视频 | www.自偷自拍.com| netflix在线观看网站| 久久久久久久久免费视频了| svipshipincom国产片| 国产成人一区二区三区免费视频网站| a在线观看视频网站| 18禁美女被吸乳视频| 成在线人永久免费视频| 日本精品一区二区三区蜜桃| 在线免费观看的www视频| av在线播放免费不卡| 97人妻精品一区二区三区麻豆 | 天天躁狠狠躁夜夜躁狠狠躁| 不卡一级毛片| 欧美三级亚洲精品| 婷婷亚洲欧美| 欧美日韩精品网址| 欧美精品啪啪一区二区三区| 久久精品亚洲精品国产色婷小说| 91在线观看av| 少妇熟女aⅴ在线视频| 欧美 亚洲 国产 日韩一| 亚洲国产日韩欧美精品在线观看 | 亚洲国产欧美网| 成人精品一区二区免费| 日韩大尺度精品在线看网址| 精品国内亚洲2022精品成人| 成人一区二区视频在线观看| 麻豆一二三区av精品| 日本五十路高清| 999精品在线视频| 国产一级毛片七仙女欲春2 | 亚洲黑人精品在线| 夜夜爽天天搞| 中国美女看黄片| 99riav亚洲国产免费| 成人手机av| 在线观看免费午夜福利视频| 午夜老司机福利片| 一进一出抽搐动态| 可以在线观看的亚洲视频| 中文字幕人妻熟女乱码| 亚洲人成电影免费在线| 欧美人与性动交α欧美精品济南到| 两性夫妻黄色片| 香蕉国产在线看| 国产真实乱freesex| 精品一区二区三区四区五区乱码| 无人区码免费观看不卡| 欧美日韩一级在线毛片| 别揉我奶头~嗯~啊~动态视频| 黄色片一级片一级黄色片| 99久久精品国产亚洲精品| 亚洲欧美精品综合久久99| 精品日产1卡2卡| 国产精品永久免费网站| 午夜久久久久精精品| 精品第一国产精品| 午夜免费鲁丝| 亚洲中文日韩欧美视频| 搡老熟女国产l中国老女人| 欧美乱色亚洲激情| 色老头精品视频在线观看| 777久久人妻少妇嫩草av网站| 国产黄片美女视频| 麻豆久久精品国产亚洲av| 欧美又色又爽又黄视频| 久久人人精品亚洲av| 一级a爱片免费观看的视频| 色综合亚洲欧美另类图片| 美女高潮喷水抽搐中文字幕| 国产精品一区二区三区四区久久 | 一本久久中文字幕| 好男人在线观看高清免费视频 | 精品国产亚洲在线| 一二三四在线观看免费中文在| 欧美精品啪啪一区二区三区| 正在播放国产对白刺激| 亚洲男人的天堂狠狠| 国产片内射在线| 一级毛片女人18水好多| 人人妻人人看人人澡| av电影中文网址| 又紧又爽又黄一区二区| 亚洲精品国产区一区二| 999久久久国产精品视频| 男女午夜视频在线观看| 久久亚洲精品不卡| 亚洲成av片中文字幕在线观看| 极品教师在线免费播放| av视频在线观看入口| 国产精华一区二区三区| 桃红色精品国产亚洲av| 无人区码免费观看不卡| 国产单亲对白刺激| 观看免费一级毛片| 欧美日韩精品网址| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品aⅴ一区二区三区四区| 国产成年人精品一区二区| 女生性感内裤真人,穿戴方法视频| 亚洲无线在线观看| 亚洲欧美一区二区三区黑人| 麻豆av在线久日| 久9热在线精品视频| 婷婷丁香在线五月| 我的亚洲天堂| 亚洲国产日韩欧美精品在线观看 | 亚洲国产精品合色在线| 国产精品亚洲一级av第二区| 久久久久久人人人人人| 国产野战对白在线观看| 国产主播在线观看一区二区| 日韩欧美三级三区| videosex国产| 免费搜索国产男女视频| 好男人在线观看高清免费视频 | 久久精品人妻少妇| 中文字幕人成人乱码亚洲影| 精品乱码久久久久久99久播| 亚洲av成人av| 性色av乱码一区二区三区2| 久久精品成人免费网站| av福利片在线| 色av中文字幕| 成人18禁高潮啪啪吃奶动态图| 免费在线观看成人毛片| 欧美成人免费av一区二区三区| 两性夫妻黄色片| 人人妻人人澡欧美一区二区| or卡值多少钱| av视频在线观看入口| 亚洲色图 男人天堂 中文字幕| 成年免费大片在线观看| 国产精品 欧美亚洲| 精品高清国产在线一区| 国产成人欧美| 美女大奶头视频| 免费电影在线观看免费观看| 日韩三级视频一区二区三区| 成人午夜高清在线视频 | 女性被躁到高潮视频| 动漫黄色视频在线观看| 大型黄色视频在线免费观看| 搡老岳熟女国产| 国产欧美日韩精品亚洲av| 搡老妇女老女人老熟妇| 91国产中文字幕| 老司机午夜十八禁免费视频| 美女扒开内裤让男人捅视频| 国产精品电影一区二区三区| 亚洲一区二区三区不卡视频| 久热爱精品视频在线9| 亚洲精品在线美女| 99国产极品粉嫩在线观看| 国产私拍福利视频在线观看| 国产一区二区三区在线臀色熟女| 又大又爽又粗| 亚洲一区高清亚洲精品| 九色国产91popny在线| av电影中文网址| 999久久久精品免费观看国产| 亚洲国产精品999在线| 免费电影在线观看免费观看| 香蕉丝袜av| 亚洲成av人片免费观看| 精品久久蜜臀av无| 国产成人精品久久二区二区免费| 亚洲精品在线美女| 99re在线观看精品视频| 亚洲狠狠婷婷综合久久图片| 亚洲成av人片免费观看| 法律面前人人平等表现在哪些方面| 亚洲国产日韩欧美精品在线观看 | 在线观看免费视频日本深夜| 人人妻,人人澡人人爽秒播| 亚洲,欧美精品.| cao死你这个sao货| 国产成人一区二区三区免费视频网站| 亚洲天堂国产精品一区在线| 少妇熟女aⅴ在线视频| www.自偷自拍.com| 欧美激情 高清一区二区三区| 亚洲欧美一区二区三区黑人| 美女 人体艺术 gogo| 久久人妻av系列| 不卡av一区二区三区| 国产爱豆传媒在线观看 | 午夜福利欧美成人| 午夜福利免费观看在线| 两性夫妻黄色片| www.熟女人妻精品国产| www日本在线高清视频| 男人舔奶头视频| 天堂√8在线中文| a级毛片在线看网站| 可以免费在线观看a视频的电影网站| 精品高清国产在线一区| 成人av一区二区三区在线看| 久久久水蜜桃国产精品网| 久久精品成人免费网站| 欧美一级毛片孕妇| 精品国产国语对白av| 黄色 视频免费看| 午夜免费观看网址| 国产成人一区二区三区免费视频网站| 人妻丰满熟妇av一区二区三区| 久久久久国产精品人妻aⅴ院| 人人妻,人人澡人人爽秒播| 日韩欧美三级三区| 男女视频在线观看网站免费 | 桃红色精品国产亚洲av| 成人亚洲精品av一区二区| 亚洲午夜理论影院| 精品电影一区二区在线| 男女床上黄色一级片免费看| 国产精品久久视频播放| 成人亚洲精品av一区二区| av电影中文网址| 嫩草影视91久久| 国产av不卡久久| 久久久久九九精品影院| 日本在线视频免费播放| 亚洲最大成人中文| 精品久久久久久,| 黑丝袜美女国产一区| 真人一进一出gif抽搐免费| 亚洲天堂国产精品一区在线| 国产精品一区二区精品视频观看| 男女午夜视频在线观看| 国产精品爽爽va在线观看网站 | www.熟女人妻精品国产| 久久精品亚洲精品国产色婷小说| 免费在线观看视频国产中文字幕亚洲| 国产伦在线观看视频一区| 久久国产精品人妻蜜桃| www.精华液| a在线观看视频网站| www.999成人在线观看| 精品免费久久久久久久清纯| 男人操女人黄网站| 精品一区二区三区av网在线观看| 91字幕亚洲| 亚洲精品美女久久av网站| 国产精华一区二区三区| 黄片播放在线免费| 亚洲成人免费电影在线观看| 久久精品aⅴ一区二区三区四区| 亚洲七黄色美女视频| 首页视频小说图片口味搜索| 国产伦一二天堂av在线观看| 波多野结衣高清无吗| 欧美日韩中文字幕国产精品一区二区三区| 人成视频在线观看免费观看| 久久精品国产99精品国产亚洲性色| 中文资源天堂在线| 97碰自拍视频| 国产一区二区三区在线臀色熟女| 久99久视频精品免费| 免费一级毛片在线播放高清视频| 老司机靠b影院| 亚洲精品美女久久av网站| 久9热在线精品视频| 精品国内亚洲2022精品成人| 12—13女人毛片做爰片一| 欧美黑人欧美精品刺激| 久久久久九九精品影院| 精品一区二区三区av网在线观看| 夜夜看夜夜爽夜夜摸| 亚洲精品色激情综合| 麻豆国产av国片精品| 日日夜夜操网爽| 亚洲成人精品中文字幕电影| 99国产精品一区二区蜜桃av| 国内久久婷婷六月综合欲色啪| 一边摸一边做爽爽视频免费| 十八禁网站免费在线| 麻豆成人午夜福利视频| 波多野结衣巨乳人妻| 无遮挡黄片免费观看| av视频在线观看入口| 伦理电影免费视频| 亚洲精品国产区一区二| 久久香蕉国产精品| 午夜老司机福利片| 国产伦在线观看视频一区| 国产精品久久久av美女十八| 人成视频在线观看免费观看| 亚洲天堂国产精品一区在线| 少妇粗大呻吟视频| 国产精品自产拍在线观看55亚洲| 欧美成人免费av一区二区三区| 久久国产精品影院| 十八禁网站免费在线| 亚洲人成伊人成综合网2020| 欧美性长视频在线观看| 香蕉av资源在线| 在线观看日韩欧美| 国产伦一二天堂av在线观看| 国产高清激情床上av| 久久九九热精品免费| 波多野结衣巨乳人妻| 欧美在线黄色| 欧美乱妇无乱码| 亚洲国产日韩欧美精品在线观看 | 99久久久亚洲精品蜜臀av| 日本 av在线| a级毛片在线看网站| svipshipincom国产片| netflix在线观看网站| 国产成人系列免费观看| 国产精品 国内视频| 我的亚洲天堂| 波多野结衣高清无吗| 日韩av在线大香蕉| 欧美一级a爱片免费观看看 | 精品国产超薄肉色丝袜足j| 亚洲精品国产精品久久久不卡| 日本撒尿小便嘘嘘汇集6| 亚洲精品一卡2卡三卡4卡5卡| 国产av不卡久久| 女人被狂操c到高潮| 中文亚洲av片在线观看爽| 露出奶头的视频| 成人18禁在线播放| 好男人电影高清在线观看| 日本免费a在线| 黄色丝袜av网址大全| 欧美一区二区精品小视频在线| 日韩欧美免费精品| 正在播放国产对白刺激| 伊人久久大香线蕉亚洲五| 午夜福利在线观看吧| 国产成人av激情在线播放| 2021天堂中文幕一二区在线观 | 久久精品国产综合久久久| 日韩 欧美 亚洲 中文字幕| 国产精华一区二区三区| 亚洲av片天天在线观看| 国内揄拍国产精品人妻在线 | 中文字幕人妻丝袜一区二区| 久9热在线精品视频| 少妇熟女aⅴ在线视频| 亚洲久久久国产精品| 十分钟在线观看高清视频www| 婷婷亚洲欧美| 淫秽高清视频在线观看| 国内久久婷婷六月综合欲色啪| 精品电影一区二区在线| 久久精品91无色码中文字幕| 久久 成人 亚洲| 久久中文字幕人妻熟女| 中亚洲国语对白在线视频| 99精品在免费线老司机午夜| 无遮挡黄片免费观看| 国产色视频综合| 一本大道久久a久久精品| 免费电影在线观看免费观看| av超薄肉色丝袜交足视频| 一区二区三区国产精品乱码| 在线永久观看黄色视频| 两人在一起打扑克的视频| 国产视频一区二区在线看| 欧美性猛交╳xxx乱大交人| 中文字幕最新亚洲高清| 亚洲精品在线观看二区| 麻豆av在线久日| 9191精品国产免费久久| 91麻豆精品激情在线观看国产| 久久这里只有精品19| 激情在线观看视频在线高清| 国产一区二区在线av高清观看| 亚洲国产中文字幕在线视频| 97碰自拍视频| 十八禁人妻一区二区| 日本熟妇午夜| 欧美日韩中文字幕国产精品一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲精品第一综合不卡| 男人操女人黄网站| 亚洲av五月六月丁香网| 视频区欧美日本亚洲| 国产人伦9x9x在线观看| 无遮挡黄片免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国产综合亚洲| 精品久久久久久久毛片微露脸| 久久午夜综合久久蜜桃| 一级a爱片免费观看的视频| 在线观看免费午夜福利视频| 精品国产超薄肉色丝袜足j| 精品久久久久久久末码| 国产精品99久久99久久久不卡| 午夜亚洲福利在线播放| 久久久久免费精品人妻一区二区 | 少妇被粗大的猛进出69影院| xxxwww97欧美| 每晚都被弄得嗷嗷叫到高潮| 麻豆av在线久日| www日本在线高清视频| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩高清在线视频| 少妇被粗大的猛进出69影院| 可以免费在线观看a视频的电影网站| 老熟妇仑乱视频hdxx| 麻豆av在线久日| 婷婷亚洲欧美| 欧美激情高清一区二区三区| 日日爽夜夜爽网站| 88av欧美| 高清在线国产一区| 午夜a级毛片| 亚洲第一欧美日韩一区二区三区| 波多野结衣高清无吗| 午夜精品久久久久久毛片777| 中文亚洲av片在线观看爽| 免费搜索国产男女视频| 国产伦在线观看视频一区| 国产一区二区在线av高清观看| 精品久久久久久久久久免费视频| 日韩欧美一区二区三区在线观看| 中出人妻视频一区二区| 男女午夜视频在线观看| 人人妻人人澡人人看| 每晚都被弄得嗷嗷叫到高潮| 桃红色精品国产亚洲av| 黄色成人免费大全| av视频在线观看入口| 亚洲va日本ⅴa欧美va伊人久久| 极品教师在线免费播放| 99国产精品99久久久久| 亚洲熟妇熟女久久| 中文字幕久久专区| 亚洲,欧美精品.| 午夜久久久久精精品| 久久婷婷人人爽人人干人人爱| 成熟少妇高潮喷水视频| 亚洲美女黄片视频| av片东京热男人的天堂| 久久久久亚洲av毛片大全| 国产极品粉嫩免费观看在线| 老司机靠b影院| 后天国语完整版免费观看| 日本黄色视频三级网站网址| 人成视频在线观看免费观看| 亚洲国产精品合色在线| 18禁美女被吸乳视频| 久久 成人 亚洲| 19禁男女啪啪无遮挡网站| 亚洲欧美日韩高清在线视频| 亚洲,欧美精品.| 成年免费大片在线观看| 亚洲精品在线观看二区| 午夜激情福利司机影院| 国产高清视频在线播放一区| 麻豆国产av国片精品| 在线观看www视频免费| 亚洲五月婷婷丁香| 国产精品亚洲一级av第二区| 99久久无色码亚洲精品果冻| 久久久久国产精品人妻aⅴ院| 51午夜福利影视在线观看| 国产精品亚洲美女久久久| 欧美性猛交╳xxx乱大交人| 真人做人爱边吃奶动态| 国产三级黄色录像| 久久 成人 亚洲| 在线天堂中文资源库| 日日爽夜夜爽网站| 亚洲精品国产精品久久久不卡| 久久青草综合色| 中文字幕人妻丝袜一区二区| 黑人欧美特级aaaaaa片| 亚洲第一青青草原| 久久人妻av系列| 亚洲精品色激情综合| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品色激情综合| 亚洲在线自拍视频| 男女床上黄色一级片免费看| 亚洲在线自拍视频| 亚洲熟妇熟女久久| 99国产精品99久久久久| 19禁男女啪啪无遮挡网站| 国产99久久九九免费精品| 国产av不卡久久| 黄色毛片三级朝国网站| 视频区欧美日本亚洲| 看黄色毛片网站| 黑丝袜美女国产一区| 99久久99久久久精品蜜桃| 黄色视频,在线免费观看| 一本大道久久a久久精品| 18禁美女被吸乳视频| 老熟妇乱子伦视频在线观看| 狠狠狠狠99中文字幕| 欧美激情久久久久久爽电影| а√天堂www在线а√下载| 日韩精品中文字幕看吧| 啪啪无遮挡十八禁网站| av免费在线观看网站| 1024视频免费在线观看| 一二三四在线观看免费中文在| 成人午夜高清在线视频 | 亚洲精品色激情综合| 最近最新中文字幕大全电影3 | АⅤ资源中文在线天堂|