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

    180°彎曲水槽內(nèi)床面沖淤演化過程數(shù)值模擬

    2014-08-18 07:09:28,,
    長江科學(xué)院院報 2014年9期
    關(guān)鍵詞:凹岸演化過程輸沙

    , ,

    (1.天津大學(xué) 仁愛學(xué)院建筑工程系,天津 301636;

    2.大連理工大學(xué) 水利工程學(xué)院,大連 116024;

    3.長江科學(xué)院 水力學(xué)研究所,武漢 430010)

    1 研究背景

    河道內(nèi)的水流和泥沙輸運以及河床沖淤變化是水利工程中主要的研究對象,而蜿蜒曲折又是天然河道最為常見的河流形態(tài),彎曲河道內(nèi)的水流由于受離心力和水面橫比降造成的壓力梯度共同作用的影響,形成彎曲河道典型的二次流(橫向環(huán)流)特征,其與主流向流動的結(jié)合形成復(fù)雜的三維螺旋流動,具有顯著的三維特性[1-4]。

    由于二次流在河道斷面上表現(xiàn)出表層向凹岸流動,底層流向凸岸的環(huán)流特性,從而造成近河床處泥沙不斷從凹岸向凸岸輸移,加之河床沖淤過程中的橫向底坡使輸沙方向偏離水流方向,并且隨著橫向底坡的越來越大,這種影響愈發(fā)明顯,尤其對于二次流強(qiáng)度顯著和橫向坡度較大的強(qiáng)彎曲河道[5],以往二維數(shù)值模擬中半經(jīng)驗性的二次流校正[6]和考慮橫向底坡造成輸沙偏離流動方向的模型[6-8]難以適用這些問題,為此本文采用Engelund[9-10]橫向輸沙模型,并考慮河床變形所引起的臨界起動切應(yīng)力變化[11],采用流動和地形沖淤相同時間步的方法模擬180°彎曲水槽內(nèi)床面非恒定的地形沖淤演化過程。

    2 數(shù)值模型

    2.1 水動力模型

    關(guān)于水動力過程的求解,本文基于變量交錯定義的三角形網(wǎng)格,采用有限體積分層積分方法對RANS方程進(jìn)行離散[12],標(biāo)準(zhǔn)k-ε紊流模型和顯式2步投影方法求解RANS方程。首先求解包含對流項和擴(kuò)散項的動量方程,得到中間時刻的速度場;其次求解由離散連續(xù)性方程和動量方程得到的關(guān)于壓力的Poisson方程,得到壓力校正項;最后將其帶入動量方程對中間時刻速度場進(jìn)行修正,從而得到新時刻的速度場和水位[8]。

    2.2 泥沙輸運模型

    計算地形演化過程的泥沙連續(xù)性方程可以表示為

    (1)

    式中:p′表示孔隙率;qsx和qbx分別表示x方向上的懸移質(zhì)和推移質(zhì)輸沙率;qsy和qby分別表示y方向上的懸移質(zhì)和推移質(zhì)輸沙率。

    選用Van Rijn推移質(zhì)輸沙率公式[13]計算推移質(zhì)輸沙率

    圖4顯示非恒定的地形沖淤演變過程,分別顯示50,100,200,480 min 4個時刻的地形沖淤厚度等值線。從圖4中可以看出彎曲水槽的沖淤變化基本都發(fā)生在彎段內(nèi),沿彎段凸岸一側(cè)發(fā)生淤積,凹岸側(cè)發(fā)生沖刷。隨著水流的不斷作用,凸岸側(cè)不論是淤積量還是淤積范圍都在增長,相應(yīng)地凹岸一側(cè)的沖刷深度和沖刷范圍也在不斷增加,受彎道二次流橫向輸沙的影響,凹岸一側(cè)泥沙不斷地被水流搬運至凸岸側(cè),并且隨著時間推移局部淤積和沖刷量不斷增加,變化越來越緩慢,最終地形適應(yīng)水流的變化開始穩(wěn)定下來,達(dá)到平衡沖淤狀態(tài)。

    (2)

    (3)

    隨著數(shù)字網(wǎng)絡(luò)化技術(shù)的發(fā)展,圖書館的管理模式、服務(wù)功能也發(fā)生了巨大變化,從單一的靜默閱讀服務(wù),發(fā)展成為團(tuán)隊研討,為雙創(chuàng)發(fā)展提供共享學(xué)習(xí)空間。為此,圖書館辦館方針任務(wù)也要隨之調(diào)整,要打造特色空間、智慧空間,整合文獻(xiàn)資源,建設(shè)特色資源,提高館藏利用率,處理好傳統(tǒng)文獻(xiàn)存儲空間與特色共享學(xué)習(xí)空間的相互依存、相互促進(jìn)的關(guān)系。這也是我們特色建館的關(guān)鍵。

    (4)

    (5)

    Z100 3.165 mmol Zn(NO3)2·6H2O溶解于40 mL MeOH/DI混合溶液,25.32 mmol HmIM溶解于20 mL MeOH,前者迅速倒入后者,并持續(xù)攪拌120 min。離心操作后于60 ℃真空烘箱中干燥24 h得白色粉末[13]。

    彎段相對長度Lc/B=7.85,符合Leopold[16]定義的相對長度<10的自由蜿蜒河道特征,寬度水深之比B/h=8.5,屬于相對窄深一類彎道,曲率半徑寬度之比Rc/B=2.5,而Zeng[5]認(rèn)為(Rc/B)<3屬于強(qiáng)彎曲渠道,并且參照Vasquez[6]對二次流強(qiáng)度(離心加速度和重力加速度之比)的定義ac/g=0.006,可以認(rèn)為該彎曲水槽具有比較強(qiáng)的二次流效應(yīng),從而對其中的流動和輸沙計算帶來一定的難度。

    (6)

    式中β2表示橫向底坡。

    評定兩組患者的疼痛情況,采用VRS 法對患者的疼痛情況進(jìn)行全面性的評估:Ⅰ級表示為輕度疼痛,對患者的睡眠和生活的不會造成嚴(yán)重性的影響,患者也不需要服用止痛藥物;Ⅱ級表示為中度疼痛,患者的疼痛程度比較明顯,無法忍受,對患者的生活和睡眠造成了影響,需要服用止痛藥物;Ⅲ級表示為重度疼痛,疼痛程度非常劇烈,對患者的生活和睡眠造成了嚴(yán)重的影響,并且還伴隨被動體位、自主神經(jīng)紊亂等現(xiàn)象,需要給予藥物止痛。

    對于彎曲水槽二次流輸沙所形成的橫向底坡,其不僅造成輸沙方向的偏移,而且還是河床平衡二次流逆底坡輸沙,從而達(dá)到?jīng)_淤平衡的關(guān)鍵因素,本文采用方便實用的Engelund方法[9-10]考慮橫向底坡對輸沙方向的影響。

    (7)

    式中:φ表示橫向底坡造成輸沙方向偏移的角度;?zb/?n表示橫向底坡;Θ表示Shields數(shù)。

    3 彎曲水槽試驗?zāi)M

    3.1 試驗布置

    圖1 LFM 180°彎曲水槽

    應(yīng)用模型對代爾伏特理工大學(xué)(DUT)流體力學(xué)實驗室(LFM)180°彎曲水槽試驗進(jìn)行模擬[15],通過試驗測量與模擬結(jié)果的對比驗證進(jìn)一步分析地形平衡沖淤的演化過程。該水槽試驗的平面布置如圖1所示,試驗參數(shù)見表1。

    表1 LFM 180°彎曲水槽試驗參數(shù)

    在彎曲角度為180°的矩形彎段上游和下游分別連接有長度為10 m的矩形直段,彎段的曲率半徑Rc=4.25 m,寬度B=1.7 m,入口處水深h=0.2 m,流速v=0.5 m/s,上游恒定入流量Q=0.17 m3/s,形成的水面坡度為1.8‰,謝才系數(shù)C=26.4,弗勞德數(shù)Fr=0.36,彎段長度Lc=13.35 m。推移質(zhì)輸沙率qb=0.03 kg/(m·s),試驗以推移質(zhì)輸沙為主,沙粒中值粒徑d50=0.78 mm,標(biāo)準(zhǔn)偏移σg=1.15,故認(rèn)為是均勻顆粒泥沙,因此不考慮泥沙的粗化現(xiàn)象,只模擬推移質(zhì)輸沙過程。

    會計作為一門應(yīng)用型學(xué)科,其內(nèi)容從來不是一成不變的。全球經(jīng)濟(jì)形勢的不斷發(fā)展,要求會計、審計規(guī)則和方法與之相適應(yīng),其提供和反映的信息必須更可靠、更合理。人區(qū)別于機(jī)器和其他生物的最大特征之一在于人類具有優(yōu)秀的學(xué)習(xí)實踐和創(chuàng)新能力。而會計領(lǐng)域十分強(qiáng)調(diào)學(xué)習(xí)能力,繼續(xù)教育種類和有關(guān)認(rèn)證考試繁多。只有不斷學(xué)習(xí)掌握高精尖的會計領(lǐng)域?qū)I(yè)知識,具備探索和創(chuàng)新精神,才能在人工智能的浪潮中站穩(wěn)腳跟、不被淘汰。

    國外對于步進(jìn)式直線壓電驅(qū)動器的研究開始于上世紀(jì)六十年代,國內(nèi)的相關(guān)研究略晚于國外,起步于九十年代。多年來,有越來越多的大學(xué)和研究機(jī)構(gòu)深入研究了步進(jìn)式直線壓電驅(qū)動器,取得了豐富的研究成果。

    在彎曲河道沖淤過程中,二次流和河床橫向底坡的變化對推移質(zhì)輸沙量和輸沙方向起到?jīng)Q定性作用。由彎道彎曲效應(yīng)所形成的二次流在近床面處將沙粒從凹岸向凸岸輸移,從而使凹岸側(cè)地形不斷降低而凸岸測地形不斷抬升,形成橫向底坡,處于底坡上的沙粒又受到自身重力的作用有從凸岸向凹岸滑動的趨勢,如此,在二次流逆底坡輸沙和沙粒順底坡重力分量的共同作用下,沙粒在河床上發(fā)生運動,直至二者適應(yīng)當(dāng)前的水流狀態(tài)達(dá)到平衡為止,此時,地形沖淤即認(rèn)為達(dá)到?jīng)_淤平衡狀態(tài)。

    3.2 模擬效果分析

    圖2 平衡沖淤狀態(tài)下水位分布等值線

    圖2為沖淤平衡狀態(tài)下的水位等值線圖,從圖2中可以看出在水槽彎曲段凹岸和凸岸有明顯的水位差,水流由上游直段進(jìn)入彎段之后受彎段彎曲效應(yīng)的影響,在整個彎段內(nèi)開始出現(xiàn)水面橫比降,表現(xiàn)為凹岸水位抬升,凸岸水位降低,并且在彎頂前后水位差達(dá)到最大,直至出彎段進(jìn)入下游直段后,彎段彎曲效應(yīng)消失,凹岸和凸岸水位迅速恢復(fù)至順直河段水平。

    圖3為分別距離凸岸和凹岸0.34 m所選取的2個縱向剖面在給定的水流狀態(tài)下達(dá)到?jīng)_淤平衡時,計算與測量的相對水深沿程變化對比情況。圖中橫縱坐標(biāo)均進(jìn)行了無量綱化處理。橫軸表示彎段相對長度Lc/B,如圖1所示,(0~7.85)B為水槽彎曲段,之后為下游直段;縱軸表示相對水深h/h0,h0=0.2 m,相對水深<1表示河床發(fā)生淤積,相對水深>1則表示河床發(fā)生沖刷,相對水深=1則河床為不沖不淤狀態(tài)。從圖3中可以看出凸岸一側(cè)發(fā)生淤積,而凹岸一側(cè)則發(fā)生沖刷,其中最大淤積和最大沖刷位置對稱分布在彎段(2~3)B之間,此外在彎段與下游直段的交界位置8B前后也有局部的沖淤極值對稱分布。數(shù)值模擬分別采用Struiksma和Engelund兩種考慮推移質(zhì)導(dǎo)致輸沙方向偏移方法,除Struiksma方法有較大差異之外,Engelund方法結(jié)果總體上模擬出了剖面沿程的沖淤變化情況,并且除模擬的最大淤積位置稍微靠下游之外,不論是最大沖刷位置還是局部沖淤極值位置都有不錯的模擬效果,而且沖淤變化量也基本與試驗測量值保持接近。只是在彎段下半段6B前后凸岸側(cè)淤積量和凹岸側(cè)沖刷量均達(dá)到最小,模擬結(jié)果雖然也體現(xiàn)出了這種趨勢,但與測量結(jié)果仍然存在一定的差距,尤其以凹岸側(cè)更甚。在出彎段進(jìn)入下游直段之后(8B之后),凸岸淤積和凹岸沖刷量陡然減少,甚至在接近10B位置凹岸凸岸沖淤態(tài)勢發(fā)生扭轉(zhuǎn),表現(xiàn)為凸岸沖刷而凹岸淤積,顯然,凸岸測的模擬結(jié)果要優(yōu)于凹岸一側(cè)。

    式中:?表示沙粒休止角;β1表示縱向底坡。

    考慮地形沖淤所造成的底坡變化對平底河床臨界起動切應(yīng)力的影響,采用Van Rijn的修正公式[11, 14]對平底河床臨界起動切應(yīng)力進(jìn)行修正

    圖3 平衡沖淤狀態(tài)下縱向剖面相對水深計算與測量結(jié)果對比

    3.3 沖淤演化過程

    模型采用非結(jié)構(gòu)化三角形網(wǎng)格,計算區(qū)域共劃分7 020個三角形單元,3 756個節(jié)點。模擬參照表1所示流動和泥沙參數(shù)給定相應(yīng)的初始和邊界條件,固壁邊界采用滑移不穿透邊界條件,為了減少邊界條件對流動過程的影響上下游直段均延長至10 m,在上游入口處給定流量邊界,下游出口給定常數(shù)水位邊界,進(jìn)出口保持試驗平衡輸沙率0.034 kg/(m·s),糙率n=0.025,計算時間步長Δt=0.01 s,直至480 min后沖淤達(dá)到平衡為止,即(?zb/?t→ 0)。

    由于流動變化的時間尺度小于地形沖淤變化的時間尺度,采用流動和地形沖淤不同時間步的模型不能真實地模擬出河床沖淤演化過程[14],本文流動和沖淤模型采用相同時間步(即在每一個時間步內(nèi),計算出流場、懸移質(zhì)和推移質(zhì)輸沙的變化后,即對地形進(jìn)行修改,而不是流場和推移質(zhì)輸沙計算若干步之后再對地形進(jìn)行修改)模擬出平衡沖淤前的河床沖淤演化過程。

    圖4 地形沖淤演化過程等值線

    圖5 河床切應(yīng)力演化過程等值線

    不比不知道,一比嚇一跳。在茄子移栽后34天,洋豐百倍邦套餐肥展現(xiàn)了強(qiáng)勁效果。觀摩會當(dāng)天,眾人冒著高溫,在新洋豐工作人員帶領(lǐng)下來到茄子示范田。剛進(jìn)入茄子地,眾人就被左右兩塊茄子地長勢差距震驚了,示范田茄子植株明顯比對照田要高出一大截。在眾人驚奇之余,盧瓊珍高興地為觀摩人員介紹示范田和對照田的施肥方案。他一共租80畝地,其中一半轉(zhuǎn)租出去。因盧瓊珍是新洋豐肥料零售商,剩余40畝地在新洋豐技術(shù)員指導(dǎo)下,施用了百倍邦套餐肥。

    圖5顯示非恒定的河床切應(yīng)力演化過程,分別包含50,100,200,480 min 4個時刻的床面切應(yīng)力等值線,表示近河床處的水流強(qiáng)度情況。最大和最小切應(yīng)力始終對稱分布在彎段出口位置,且凹岸一側(cè)大于凸岸一側(cè),除了彎段內(nèi)切應(yīng)力為3等值線不斷從靠近凸岸側(cè)偏向凹岸一側(cè)之外,其余等值線隨時間變化不十分明顯,切應(yīng)力為3的等值線的遷移過程與彎段內(nèi)主流從凸岸向凹岸遷移的規(guī)律一致。

    圖6選取試驗最大橫向坡度發(fā)生57°斷面的橫向地形沖淤演化過程,圖中橫軸-0.85≤ξ≤0.85,指向凹岸為正,無量化處理表示從凸岸到凹岸的相對距離。從最開始第10 min緊靠凸岸岸壁位置少量的淤積至50 min凸岸淤積量和淤積范圍的急劇增長,再到100,200,480 min淤積量和淤積范圍變小并且趨于穩(wěn)定,不沖不淤位置從最初的-0.4遷移至-0.1,且不斷向凹岸側(cè)遷移,最終穩(wěn)定在軸線位置。從圖6中同樣可以看出,該斷面的橫向坡度隨時間推移越來越大,因為平衡沖淤是在二次流逆底坡輸沙和沙粒順底坡重力分量共同作用下達(dá)成,橫向坡度的不斷增大,代表沙粒重力分量平衡二次流輸沙的能力不斷增強(qiáng),最終二次流的逆底坡輸沙被沙粒重力分量順底坡輸移作用完全抵消,即達(dá)到?jīng)_淤平衡狀態(tài)。

    就像文章開頭里面提到的一樣,過去的一年里,Sarah憑借著超人的毅力邁出了事業(yè)家庭雙豐收。與此同時,女性葡萄酒專家的身影在葡萄酒行業(yè)上也越來越常見。特別是目前在香港常駐的三位MW均為女性,這會是一種巧合嗎?

    圖6 57°斷面不同時刻橫向相對水深變化過程

    4 結(jié) 論

    文章采用三維水動力和泥沙輸運模型模擬彎曲水槽內(nèi)的三維流動、泥沙輸運和內(nèi)床面沖淤演化問題,直接采用三維模型模擬流動的二次流特性,避免二維模型中對二次流的半經(jīng)驗校正,對于二次流強(qiáng)度較強(qiáng),橫向底坡坡度較大的問題采用Engelund橫向輸沙模型模擬輸沙方向的偏移角度,采用Van Rijn對河床臨界起動切應(yīng)力進(jìn)行修正,流動和沖淤同步的模擬方法模擬非恒定的地形沖淤演化過程。得到的主要結(jié)論有:

    (1) 模型成功模擬了180°彎曲水槽的主要沖淤特性,基于非結(jié)構(gòu)化網(wǎng)格的三維地形沖淤模型在實際蜿蜒河道的數(shù)值模擬中有廣闊的應(yīng)用前景。

    (2) 對于強(qiáng)彎曲、橫向底坡較大的水槽沖淤問題,不需人為指定參數(shù),方便實用的Engelund橫向輸沙模型以及考慮對河床臨界起動切應(yīng)力進(jìn)行修正方法的應(yīng)用,能取得不錯的模擬效果。

    禮貌和不禮貌的研究傾向于忽略面子威脅目標(biāo)對象的回應(yīng)。話語分析者進(jìn)行語篇分析時,揭示了受話人對于不禮貌的回應(yīng)關(guān)系到對話的進(jìn)展。因此,Culpeper et al.(2003)提出了語篇分析中不禮貌話語回應(yīng)策略的框架。

    (3) 三維流動和地形沖淤相同時間步計算的模型能有效地模擬出非恒定的地形沖淤過程,有助于認(rèn)識沖淤平衡前的地形演化過程。

    參考文獻(xiàn):

    [1] 吳持恭,水力學(xué)[M].北京:高等教育出版社, 2003.(WU Chi-gong.Hydraulics[M].Beijing: Higher Education Press, 2003.(in Chinese))

    [2] 張武紅, 呂 昕.彎道水力學(xué)[M].北京:水利電力出版社, 1993.(ZHANG Wu-hong, LV Xin.Hydraulics of Curved Channel[M].Beijing: Water Resources and Electric Power Press, 1993.(in Chinese))

    [3] 董耀華.彎道水流的基本特性及數(shù)值模擬[J].長江科學(xué)院院報,1996,13(1): 25-29.(DONG Yao-hua.Essential Features and Numerical Modeling of Bend Flow[J].Journal of Yangtze River Scientific Research Institute, 1996, 13(1): 25-29.(in Chinese))

    [4] 趙根生,盧金友,汪 鵬.連續(xù)彎道演變機(jī)理研究綜述[J].長江科學(xué)院院報,2009,26(6): 1-3.(ZHAO Gen-sheng, LU Jin-you, WANG Peng.Review on Research for Evolution of Continuously Curved Channels[J].Journal of Yangtze River Scientific Research Institute, 2009, 26(6): 1-3.(in Chinese))

    [5] ZENG J, CONSTANTINESCU G, WEBER L.A 3D Non-hydrostatic Model to Predict Flow and Sediment Transport in Loose-bed Channel Bends[J].Journal of Hydraulic Research, 2008,46(3): 356-372.

    [6] VASQUEZ J A, MILLAR R G, STEFFLER P M.River2D Morphology, Part II: Curved Alluvial Channels[C]∥Canadian Society for Civil Engineering(CSCE).17th Canadian Hydrotechnical Conference, Edmonton, Alberta, August 17-19, 2005: 259-268.

    [7] ZENG J, CONSTANTINESCU G, WEBER L.Prediction of Equilibrium Flow and Bed Load Transport in a Curved Bend[C]∥American Society of Civil Engineers(ASCE) World Environmental and Water Resource Congress Omaha, Nebraska, United States[J].May 21-25, 2006: 1-11.

    [8] 邢 巖, 金 生, 艾叢芳, 等.明渠彎道三維沖淤數(shù)值模擬[J].泥沙研究,2013,(6): 1-5.(XING Yan, JIN Sheng, AI Cong-fang,etal.Three-dimensional Numerical Modeling of Bed Evolution in Open Channel Bends[J].Journal of Sediment Research, 2013, (6): 1-5.(in Chinese))

    [9] ENGELUND F.Flow and Bed Topography in Channel Bends[J].Journal of the Hydraulics Division,1974,100(11): 1631-1648.

    [10]RUTHER N, OLSEN N R.Three-dimensional Modeling of Sediment Transport in a Narrow 90 Channel Bend[J].Journal of Hydraulic Engineering, 2005,131(10): 917-920.

    [11]DUAN J G, JULIEN P Y.Numerical Simulation of the Inception of Channel Meandering[J].Earth Surface Processes and Landforms, 2005, 30(9): 1093-1110.

    [12]邢 巖, 艾叢芳, 金 生.強(qiáng)彎渠道三維水流數(shù)值模擬.水道港口[J].2013,34(4): 335-343.(XING Yan, AI Cong-fang, JIN Sheng.Three-dimensional Modeling of Flow in a Sharp Open Channel Bend[J].Journal of Waterway and Harbor, 2013, 34(4): 335-343.(in Chinese))

    [13]WU W.Computational River Dynamics[M].London: Taylor & Francis, 2008.

    [14]KHOSRONEJAD A, RENNIE C D, SALEHI N S,etal.3D Numerical Modeling of Flow and Sediment Transport in Laboratory Channel Bends[J].Journal of Hydraulic Engineering, 2007,133(10): 1123-1134.

    [15]STRUIKSMA N, OLESEN K W, FLOKSTRA C,etal.Bed Deformation in Curved Alluvial Channels[J].Journal of Hydraulic Research, 1985,23(1): 57-79.

    [16]KASSEM A A, CHAUDHRY M H.Numerical Modeling of Bed Evolution in Channel Bends[J].Journal of Hydraulic Engineering, 2002,128(5): 507-514.

    猜你喜歡
    凹岸演化過程輸沙
    生命的演化過程
    模因論視角下韓語“??”表“喝”動作演化過程研究
    時間非齊次二態(tài)量子游蕩的演化過程分析
    重慶萬盛石林的形成時代及發(fā)育演化過程
    河流凹凸岸芻議
    例談河流凹凸岸的判斷方法
    例談河流凹凸岸的判斷方法
    湖南省四水流域森林、徑流、輸沙變化特性
    基于分布式水文模型的流域輸沙過程模擬
    塔克拉瑪干沙漠腹地塔中地區(qū)的風(fēng)動力環(huán)境
    成年av动漫网址| 伦理电影大哥的女人| 亚洲欧美清纯卡通| 日本wwww免费看| 国产 一区 欧美 日韩| 乱系列少妇在线播放| 精品久久久久久久久亚洲| 亚洲三级黄色毛片| 有码 亚洲区| 中文字幕av成人在线电影| 伊人久久精品亚洲午夜| 18禁动态无遮挡网站| 国产精品精品国产色婷婷| 伦理电影大哥的女人| 成年av动漫网址| 在现免费观看毛片| 亚洲精品乱码久久久久久按摩| 高清在线视频一区二区三区| 六月丁香七月| 色吧在线观看| 九九在线视频观看精品| 又黄又爽又刺激的免费视频.| 日韩欧美精品免费久久| 国产乱人视频| 中国美白少妇内射xxxbb| 成人二区视频| 日韩欧美精品免费久久| 综合色av麻豆| 久久99热这里只有精品18| av在线观看视频网站免费| 免费看光身美女| 国产爽快片一区二区三区| 国产淫片久久久久久久久| 在线观看美女被高潮喷水网站| 一级片'在线观看视频| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久久久久精品电影| 岛国毛片在线播放| 国产一区有黄有色的免费视频| 国产精品一二三区在线看| 偷拍熟女少妇极品色| 少妇人妻久久综合中文| 我的老师免费观看完整版| 国产成人午夜福利电影在线观看| 蜜桃亚洲精品一区二区三区| 免费观看的影片在线观看| 成人国产av品久久久| 天天躁夜夜躁狠狠久久av| h日本视频在线播放| 亚洲国产精品999| 日韩视频在线欧美| 秋霞在线观看毛片| 下体分泌物呈黄色| 毛片女人毛片| 国模一区二区三区四区视频| 免费观看无遮挡的男女| 深爱激情五月婷婷| 成人漫画全彩无遮挡| 国产亚洲最大av| 一级a做视频免费观看| 国产精品一区二区性色av| 特级一级黄色大片| av线在线观看网站| 国产老妇伦熟女老妇高清| 亚洲国产av新网站| 人人妻人人澡人人爽人人夜夜| 午夜免费观看性视频| 免费观看无遮挡的男女| 视频中文字幕在线观看| 最新中文字幕久久久久| 久久久亚洲精品成人影院| 成年版毛片免费区| 国产淫语在线视频| 久热久热在线精品观看| 美女cb高潮喷水在线观看| 综合色av麻豆| 免费播放大片免费观看视频在线观看| 精品熟女少妇av免费看| 成人综合一区亚洲| 国产精品久久久久久精品古装| 欧美成人一区二区免费高清观看| 欧美国产精品一级二级三级 | 我要看日韩黄色一级片| 午夜日本视频在线| 日韩亚洲欧美综合| 成人亚洲精品av一区二区| 丰满少妇做爰视频| 亚洲内射少妇av| 亚洲精品,欧美精品| 一区二区三区精品91| 2021天堂中文幕一二区在线观| 国产黄色免费在线视频| 久久久久久久久久久免费av| 国产淫语在线视频| 一级爰片在线观看| 国产 一区 欧美 日韩| 亚洲精品一区蜜桃| 少妇熟女欧美另类| 亚洲综合精品二区| 婷婷色综合www| 日本一二三区视频观看| 大香蕉97超碰在线| 久久国内精品自在自线图片| 丝袜美腿在线中文| 日本爱情动作片www.在线观看| 色播亚洲综合网| 一本色道久久久久久精品综合| 亚洲自拍偷在线| 永久免费av网站大全| 国产精品.久久久| 欧美高清成人免费视频www| 亚洲天堂国产精品一区在线| 亚洲精品国产av蜜桃| 秋霞在线观看毛片| 午夜爱爱视频在线播放| 好男人视频免费观看在线| 亚洲综合精品二区| 国产av国产精品国产| 亚洲av一区综合| 国产精品人妻久久久影院| 免费av毛片视频| 只有这里有精品99| 国产精品.久久久| 精品视频人人做人人爽| 少妇人妻精品综合一区二区| 午夜福利视频1000在线观看| 国模一区二区三区四区视频| 校园人妻丝袜中文字幕| 国产视频首页在线观看| 午夜激情久久久久久久| 18禁裸乳无遮挡动漫免费视频 | 国产精品av视频在线免费观看| 亚洲精品,欧美精品| 欧美精品一区二区大全| 精品国产露脸久久av麻豆| 久久这里有精品视频免费| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久精品古装| 日韩一区二区三区影片| 亚洲最大成人手机在线| av国产精品久久久久影院| a级毛片免费高清观看在线播放| 欧美性猛交╳xxx乱大交人| 91在线精品国自产拍蜜月| 尾随美女入室| 97在线人人人人妻| 国产亚洲av片在线观看秒播厂| 成人综合一区亚洲| 小蜜桃在线观看免费完整版高清| 国产又色又爽无遮挡免| 少妇的逼水好多| 免费看不卡的av| 尾随美女入室| 日韩国内少妇激情av| 大陆偷拍与自拍| 91久久精品电影网| av播播在线观看一区| 最新中文字幕久久久久| 日韩大片免费观看网站| 国产高清国产精品国产三级 | 亚洲欧洲日产国产| 久久久精品欧美日韩精品| 午夜亚洲福利在线播放| 性插视频无遮挡在线免费观看| 乱系列少妇在线播放| 欧美 日韩 精品 国产| av黄色大香蕉| 看黄色毛片网站| 看免费成人av毛片| 欧美区成人在线视频| 国产成人a∨麻豆精品| 欧美成人a在线观看| 国产精品久久久久久精品古装| 久久久精品免费免费高清| 日韩国内少妇激情av| av在线观看视频网站免费| 九色成人免费人妻av| 欧美精品国产亚洲| 久久精品综合一区二区三区| 91在线精品国自产拍蜜月| 久久久久网色| 在线免费十八禁| 亚洲欧美日韩另类电影网站 | 日本wwww免费看| xxx大片免费视频| 嫩草影院精品99| 亚洲欧美日韩卡通动漫| 欧美成人一区二区免费高清观看| 国内精品宾馆在线| 免费黄色在线免费观看| 成年免费大片在线观看| 日韩一区二区三区影片| 精品99又大又爽又粗少妇毛片| 国产亚洲av片在线观看秒播厂| 黄色欧美视频在线观看| 18+在线观看网站| 成人高潮视频无遮挡免费网站| 人体艺术视频欧美日本| 男的添女的下面高潮视频| 亚洲天堂av无毛| 国产黄a三级三级三级人| 热99国产精品久久久久久7| 亚洲人与动物交配视频| 亚洲久久久久久中文字幕| 日韩一区二区视频免费看| 99精国产麻豆久久婷婷| 日本色播在线视频| 国产国拍精品亚洲av在线观看| 一区二区三区免费毛片| 建设人人有责人人尽责人人享有的 | 亚洲激情五月婷婷啪啪| 夜夜看夜夜爽夜夜摸| 国产成人精品久久久久久| 哪个播放器可以免费观看大片| 中文欧美无线码| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产欧美人成| 99九九线精品视频在线观看视频| 一本久久精品| 国产中年淑女户外野战色| 亚洲天堂国产精品一区在线| 在线播放无遮挡| h日本视频在线播放| 在线看a的网站| 欧美日韩视频精品一区| 一本色道久久久久久精品综合| 亚洲av电影在线观看一区二区三区 | 午夜爱爱视频在线播放| 国产高清三级在线| 亚洲无线观看免费| 国产毛片在线视频| 免费看av在线观看网站| 又爽又黄a免费视频| 亚洲精华国产精华液的使用体验| 人体艺术视频欧美日本| 只有这里有精品99| 一本色道久久久久久精品综合| 免费大片18禁| 亚洲欧美中文字幕日韩二区| 成人一区二区视频在线观看| 岛国毛片在线播放| 只有这里有精品99| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av免费在线观看| 美女主播在线视频| 精华霜和精华液先用哪个| 久久久亚洲精品成人影院| eeuss影院久久| 久久ye,这里只有精品| 亚洲国产最新在线播放| 如何舔出高潮| 国产午夜精品一二区理论片| 亚洲av日韩在线播放| 最近最新中文字幕免费大全7| 午夜福利视频精品| 久久精品久久久久久久性| 成人午夜精彩视频在线观看| 大话2 男鬼变身卡| 国产精品人妻久久久影院| 色视频在线一区二区三区| 国产美女午夜福利| 免费在线观看成人毛片| av国产精品久久久久影院| 亚洲av免费高清在线观看| 亚洲精品中文字幕在线视频 | 欧美高清成人免费视频www| 91狼人影院| 麻豆久久精品国产亚洲av| 久久99精品国语久久久| 高清av免费在线| 18禁裸乳无遮挡免费网站照片| 如何舔出高潮| 青春草视频在线免费观看| 在线天堂最新版资源| 久久99热这里只有精品18| 啦啦啦中文免费视频观看日本| 国产一级毛片在线| 伦精品一区二区三区| 午夜福利在线在线| 国产亚洲91精品色在线| 精品久久久精品久久久| 亚洲不卡免费看| www.av在线官网国产| 亚洲欧美日韩无卡精品| 激情五月婷婷亚洲| 亚洲最大成人手机在线| 熟妇人妻不卡中文字幕| 91狼人影院| 亚洲精品aⅴ在线观看| 黑人高潮一二区| 久久6这里有精品| 国产精品99久久99久久久不卡 | 免费不卡的大黄色大毛片视频在线观看| 男插女下体视频免费在线播放| 欧美极品一区二区三区四区| 亚洲国产精品国产精品| 内射极品少妇av片p| 国产免费福利视频在线观看| 亚洲美女视频黄频| 欧美少妇被猛烈插入视频| 亚洲精品国产av成人精品| 内地一区二区视频在线| 国产精品人妻久久久久久| 欧美三级亚洲精品| 日韩伦理黄色片| 美女xxoo啪啪120秒动态图| 亚洲最大成人手机在线| 久久精品国产亚洲av天美| 97超碰精品成人国产| 麻豆久久精品国产亚洲av| 99久久精品热视频| 观看美女的网站| 在线免费观看不下载黄p国产| 久久久久久久久久久丰满| 亚洲性久久影院| 国产黄色免费在线视频| 中文精品一卡2卡3卡4更新| 亚洲欧美清纯卡通| 亚洲国产精品国产精品| 欧美+日韩+精品| 欧美精品人与动牲交sv欧美| 五月玫瑰六月丁香| 久久精品夜色国产| 亚洲,欧美,日韩| 哪个播放器可以免费观看大片| 久久精品久久精品一区二区三区| 国产精品蜜桃在线观看| 免费av不卡在线播放| 99热网站在线观看| 禁无遮挡网站| 中文字幕久久专区| 亚洲av中文av极速乱| 一级二级三级毛片免费看| 一本一本综合久久| 涩涩av久久男人的天堂| 免费看光身美女| 国产极品天堂在线| 日韩强制内射视频| 99热6这里只有精品| 少妇的逼水好多| 日产精品乱码卡一卡2卡三| 国产精品爽爽va在线观看网站| 51国产日韩欧美| 国产欧美日韩精品一区二区| 欧美精品人与动牲交sv欧美| 成人特级av手机在线观看| 99热全是精品| 亚洲性久久影院| 在线观看一区二区三区| 久久6这里有精品| 亚洲人成网站在线播| 伦理电影大哥的女人| 精品久久久久久电影网| 插阴视频在线观看视频| 国产一区有黄有色的免费视频| 国产午夜福利久久久久久| 国产男人的电影天堂91| 免费观看的影片在线观看| 亚洲三级黄色毛片| 91久久精品电影网| 又粗又硬又长又爽又黄的视频| 在现免费观看毛片| 黄色欧美视频在线观看| 中文天堂在线官网| 欧美zozozo另类| 亚洲精品一区蜜桃| 肉色欧美久久久久久久蜜桃 | 久久精品久久久久久久性| 久久久久国产网址| 九九久久精品国产亚洲av麻豆| 一区二区三区乱码不卡18| 3wmmmm亚洲av在线观看| 国产精品精品国产色婷婷| 美女视频免费永久观看网站| 99久久九九国产精品国产免费| 国产伦精品一区二区三区视频9| 日韩欧美精品免费久久| 18+在线观看网站| 建设人人有责人人尽责人人享有的 | 麻豆成人av视频| 黄片wwwwww| 日韩欧美精品免费久久| 欧美一级a爱片免费观看看| 日韩欧美 国产精品| 中国国产av一级| 久久久久久久大尺度免费视频| 看十八女毛片水多多多| 久久久久久久久大av| 色婷婷久久久亚洲欧美| 国产精品人妻久久久久久| 久久精品国产自在天天线| 91精品伊人久久大香线蕉| 亚洲天堂国产精品一区在线| 午夜福利网站1000一区二区三区| 国产视频首页在线观看| 性色av一级| 色综合色国产| 爱豆传媒免费全集在线观看| 亚洲人与动物交配视频| 搞女人的毛片| 欧美日韩精品成人综合77777| 精品久久久久久久末码| 午夜免费观看性视频| 蜜臀久久99精品久久宅男| 一区二区三区四区激情视频| 国产综合懂色| 国产永久视频网站| 一级毛片久久久久久久久女| 精品亚洲乱码少妇综合久久| 国产高潮美女av| 啦啦啦中文免费视频观看日本| 日日啪夜夜爽| 免费高清在线观看视频在线观看| 久久韩国三级中文字幕| 波多野结衣巨乳人妻| av国产免费在线观看| 色婷婷久久久亚洲欧美| 伊人久久精品亚洲午夜| 九色成人免费人妻av| 亚洲人成网站在线播| 中文天堂在线官网| 久久韩国三级中文字幕| 秋霞伦理黄片| 日韩中字成人| 国产片特级美女逼逼视频| 高清视频免费观看一区二区| 精品一区二区免费观看| 人妻系列 视频| 国产av码专区亚洲av| 久久久精品94久久精品| 男人狂女人下面高潮的视频| 三级男女做爰猛烈吃奶摸视频| 国产精品爽爽va在线观看网站| 久久精品国产亚洲av天美| 中文字幕av成人在线电影| 卡戴珊不雅视频在线播放| 午夜免费观看性视频| 国产伦精品一区二区三区视频9| 18禁裸乳无遮挡免费网站照片| 国产精品国产av在线观看| 国内揄拍国产精品人妻在线| 99久国产av精品国产电影| 日日摸夜夜添夜夜添av毛片| 日韩不卡一区二区三区视频在线| 亚洲欧美一区二区三区国产| h日本视频在线播放| 欧美另类一区| 直男gayav资源| av在线观看视频网站免费| 18禁裸乳无遮挡免费网站照片| 伊人久久精品亚洲午夜| 久久精品综合一区二区三区| 九九爱精品视频在线观看| 国产女主播在线喷水免费视频网站| 偷拍熟女少妇极品色| 亚洲人成网站高清观看| 日本与韩国留学比较| 午夜老司机福利剧场| 亚洲国产日韩一区二区| 有码 亚洲区| 内地一区二区视频在线| 精品人妻偷拍中文字幕| av国产久精品久网站免费入址| kizo精华| 成人国产av品久久久| 免费在线观看成人毛片| 十八禁网站网址无遮挡 | 国产一区有黄有色的免费视频| 97超碰精品成人国产| 插逼视频在线观看| 亚洲精品久久午夜乱码| 国产精品麻豆人妻色哟哟久久| 精品久久久久久久久亚洲| 最近2019中文字幕mv第一页| 久久久久久久久久人人人人人人| 在线观看美女被高潮喷水网站| 国产探花在线观看一区二区| 七月丁香在线播放| 男女下面进入的视频免费午夜| 欧美成人a在线观看| 免费av观看视频| xxx大片免费视频| 高清在线视频一区二区三区| 亚洲欧美日韩无卡精品| 国产美女午夜福利| 成人亚洲欧美一区二区av| 视频中文字幕在线观看| 97超碰精品成人国产| 精品久久国产蜜桃| 久久久久久久久久久丰满| 久久精品国产自在天天线| 久久久久国产精品人妻一区二区| 欧美日韩综合久久久久久| 久久99热这里只频精品6学生| 欧美潮喷喷水| 国产精品一及| 中文字幕制服av| 欧美日韩在线观看h| 国产精品国产三级国产专区5o| 国产乱人视频| 久久久久网色| 99精国产麻豆久久婷婷| 欧美性猛交╳xxx乱大交人| 夜夜爽夜夜爽视频| 日韩不卡一区二区三区视频在线| 三级国产精品片| 寂寞人妻少妇视频99o| 国产免费一区二区三区四区乱码| 特级一级黄色大片| 一级毛片电影观看| 男女那种视频在线观看| 最新中文字幕久久久久| 青春草国产在线视频| 男的添女的下面高潮视频| 国产在视频线精品| 欧美激情久久久久久爽电影| 五月玫瑰六月丁香| 伊人久久国产一区二区| 日韩伦理黄色片| 国产69精品久久久久777片| 亚洲av日韩在线播放| 日韩欧美精品免费久久| 亚洲最大成人中文| 成人特级av手机在线观看| 日韩视频在线欧美| 七月丁香在线播放| 日本三级黄在线观看| eeuss影院久久| 一级av片app| 街头女战士在线观看网站| 国产成人aa在线观看| 纵有疾风起免费观看全集完整版| 久久久久九九精品影院| 欧美区成人在线视频| 久久韩国三级中文字幕| 精品少妇黑人巨大在线播放| 欧美性感艳星| 国产精品国产av在线观看| 少妇的逼好多水| 精品久久久噜噜| 丰满人妻一区二区三区视频av| 亚洲精品一区蜜桃| 久久精品久久久久久噜噜老黄| 91精品国产九色| 亚洲自拍偷在线| 啦啦啦啦在线视频资源| 国产精品福利在线免费观看| 狂野欧美激情性bbbbbb| 91久久精品电影网| 一本一本综合久久| 久久久久久久久久人人人人人人| 国产精品久久久久久av不卡| 黄色配什么色好看| 国产黄片视频在线免费观看| 99re6热这里在线精品视频| 黄色怎么调成土黄色| 一个人看视频在线观看www免费| 国产av国产精品国产| 国产一区二区三区综合在线观看 | a级一级毛片免费在线观看| 99久久中文字幕三级久久日本| 尤物成人国产欧美一区二区三区| 久久久久精品久久久久真实原创| 午夜精品国产一区二区电影 | 色视频在线一区二区三区| 亚洲性久久影院| 午夜免费鲁丝| 特级一级黄色大片| 嫩草影院新地址| 国产欧美亚洲国产| 丝袜美腿在线中文| 老司机影院毛片| 国产精品久久久久久av不卡| 黄色日韩在线| 一个人观看的视频www高清免费观看| 女人久久www免费人成看片| 久久韩国三级中文字幕| 三级经典国产精品| 亚洲高清免费不卡视频| 久久久久久久久久人人人人人人| 亚洲aⅴ乱码一区二区在线播放| 欧美激情久久久久久爽电影| 日日摸夜夜添夜夜爱| 十八禁网站网址无遮挡 | 国产精品不卡视频一区二区| 久久久久久久午夜电影| 赤兔流量卡办理| 亚洲最大成人av| 亚洲性久久影院| 亚洲国产精品成人综合色| 国产黄片美女视频| 天天躁日日操中文字幕| 国产老妇女一区| 色婷婷久久久亚洲欧美| 22中文网久久字幕| 午夜激情福利司机影院| 国产午夜精品久久久久久一区二区三区| 香蕉精品网在线| 亚洲av成人精品一二三区| 老师上课跳d突然被开到最大视频| 又大又黄又爽视频免费| 九九爱精品视频在线观看| 草草在线视频免费看| 久久精品国产鲁丝片午夜精品| 一级毛片 在线播放| 真实男女啪啪啪动态图| 亚洲av福利一区| 亚洲丝袜综合中文字幕| a级毛片免费高清观看在线播放| 在线精品无人区一区二区三 | 亚洲欧洲日产国产| 欧美高清成人免费视频www| 少妇裸体淫交视频免费看高清| 国产成人精品福利久久| 国语对白做爰xxxⅹ性视频网站|