• <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)境
    如日韩欧美国产精品一区二区三区| 日本黄色视频三级网站网址| 18禁国产床啪视频网站| 一级片免费观看大全| 首页视频小说图片口味搜索| 免费一级毛片在线播放高清视频 | 一区二区日韩欧美中文字幕| 久久亚洲真实| 久久久国产一区二区| 午夜福利影视在线免费观看| 最好的美女福利视频网| videosex国产| av福利片在线| 91精品国产国语对白视频| 一级黄色大片毛片| 国产精品免费一区二区三区在线| 国产精品成人在线| 两个人看的免费小视频| 亚洲熟女毛片儿| 亚洲国产精品999在线| 亚洲av电影在线进入| 欧美人与性动交α欧美精品济南到| 最近最新中文字幕大全电影3 | 男女高潮啪啪啪动态图| 精品少妇一区二区三区视频日本电影| 午夜免费观看网址| 又大又爽又粗| 精品一区二区三区av网在线观看| 每晚都被弄得嗷嗷叫到高潮| 女性生殖器流出的白浆| 日本vs欧美在线观看视频| 国产1区2区3区精品| 亚洲一区二区三区色噜噜 | 欧美中文综合在线视频| 久久午夜综合久久蜜桃| 看黄色毛片网站| 999精品在线视频| 亚洲男人天堂网一区| 久久精品成人免费网站| 国产aⅴ精品一区二区三区波| 老熟妇乱子伦视频在线观看| av视频免费观看在线观看| 精品福利永久在线观看| 久久国产精品人妻蜜桃| 国产99白浆流出| 99精品欧美一区二区三区四区| 久久精品成人免费网站| 黄色怎么调成土黄色| 夜夜看夜夜爽夜夜摸 | 人妻久久中文字幕网| 不卡av一区二区三区| 国产亚洲精品久久久久久毛片| 丁香欧美五月| 久久亚洲真实| 天堂影院成人在线观看| 久久精品91蜜桃| 久久亚洲精品不卡| 这个男人来自地球电影免费观看| 亚洲情色 制服丝袜| 亚洲熟女毛片儿| 亚洲,欧美精品.| 99riav亚洲国产免费| 成人影院久久| 亚洲人成网站在线播放欧美日韩| 88av欧美| 午夜精品在线福利| 50天的宝宝边吃奶边哭怎么回事| 国产无遮挡羞羞视频在线观看| 午夜福利,免费看| 大香蕉久久成人网| 免费观看精品视频网站| 欧美+亚洲+日韩+国产| 久久久久久久午夜电影 | 99国产精品一区二区三区| 男人舔女人的私密视频| 大陆偷拍与自拍| 极品人妻少妇av视频| 国产熟女xx| 777久久人妻少妇嫩草av网站| 老司机在亚洲福利影院| 丝袜在线中文字幕| 国产男靠女视频免费网站| 国产高清激情床上av| 日韩三级视频一区二区三区| cao死你这个sao货| 很黄的视频免费| 国产高清国产精品国产三级| 亚洲av成人一区二区三| 侵犯人妻中文字幕一二三四区| 91字幕亚洲| 日韩欧美一区视频在线观看| 90打野战视频偷拍视频| 女警被强在线播放| 在线国产一区二区在线| 亚洲精品一区av在线观看| 国产av精品麻豆| 老司机午夜福利在线观看视频| 日韩国内少妇激情av| 欧美在线黄色| 亚洲aⅴ乱码一区二区在线播放 | 国产一区二区三区综合在线观看| 男人操女人黄网站| 脱女人内裤的视频| 精品电影一区二区在线| 久久久精品欧美日韩精品| av片东京热男人的天堂| 一本大道久久a久久精品| 最好的美女福利视频网| 日本a在线网址| 法律面前人人平等表现在哪些方面| 女同久久另类99精品国产91| 日本精品一区二区三区蜜桃| 99在线视频只有这里精品首页| 麻豆一二三区av精品| 精品人妻1区二区| 51午夜福利影视在线观看| 欧洲精品卡2卡3卡4卡5卡区| 女同久久另类99精品国产91| 中文字幕人妻丝袜一区二区| 久久香蕉激情| 国产成人精品久久二区二区免费| 日本五十路高清| 男女之事视频高清在线观看| 国产亚洲欧美精品永久| 国产精品免费视频内射| 日本精品一区二区三区蜜桃| 无人区码免费观看不卡| 女人爽到高潮嗷嗷叫在线视频| 天天影视国产精品| 1024视频免费在线观看| 久久国产精品人妻蜜桃| 国产高清视频在线播放一区| 成人亚洲精品一区在线观看| 久久久久久亚洲精品国产蜜桃av| 欧美精品啪啪一区二区三区| 国产av一区二区精品久久| 香蕉国产在线看| 美女福利国产在线| 日韩成人在线观看一区二区三区| 69av精品久久久久久| 91大片在线观看| 亚洲成人久久性| 一个人免费在线观看的高清视频| 国产区一区二久久| a级毛片黄视频| 岛国视频午夜一区免费看| 日韩有码中文字幕| 成人永久免费在线观看视频| 午夜福利欧美成人| 男人的好看免费观看在线视频 | 欧美成人性av电影在线观看| 国产精品永久免费网站| 色哟哟哟哟哟哟| 亚洲性夜色夜夜综合| 精品日产1卡2卡| 黄色成人免费大全| 99精品欧美一区二区三区四区| 午夜老司机福利片| 欧美 亚洲 国产 日韩一| 日韩欧美一区视频在线观看| 如日韩欧美国产精品一区二区三区| 国产成人精品久久二区二区免费| 免费少妇av软件| 日韩大码丰满熟妇| 真人做人爱边吃奶动态| 免费少妇av软件| 一区二区三区精品91| 久久久久久久久中文| 国产在线观看jvid| www.www免费av| 日日摸夜夜添夜夜添小说| 无遮挡黄片免费观看| 一级,二级,三级黄色视频| 亚洲精品久久午夜乱码| 国产野战对白在线观看| 一个人观看的视频www高清免费观看 | 亚洲色图 男人天堂 中文字幕| 亚洲午夜理论影院| 新久久久久国产一级毛片| 婷婷六月久久综合丁香| 欧美人与性动交α欧美软件| av网站在线播放免费| 91老司机精品| 精品卡一卡二卡四卡免费| 欧美日韩黄片免| 久久国产乱子伦精品免费另类| 亚洲久久久国产精品| 欧美日本中文国产一区发布| 在线看a的网站| a在线观看视频网站| 女人被狂操c到高潮| 人人澡人人妻人| 97人妻天天添夜夜摸| 在线视频色国产色| 日本精品一区二区三区蜜桃| 男女下面插进去视频免费观看| 日韩一卡2卡3卡4卡2021年| 日日爽夜夜爽网站| 色哟哟哟哟哟哟| 欧美黄色片欧美黄色片| 日韩高清综合在线| 男女高潮啪啪啪动态图| 黄色成人免费大全| 丰满的人妻完整版| 国产精品 国内视频| 午夜成年电影在线免费观看| 欧美激情久久久久久爽电影 | 亚洲熟女毛片儿| 亚洲男人天堂网一区| 超色免费av| 精品一区二区三区av网在线观看| 男女下面插进去视频免费观看| 国产欧美日韩一区二区精品| 亚洲精品国产色婷婷电影| 岛国在线观看网站| 成人三级黄色视频| 日韩视频一区二区在线观看| 国产成人av教育| 夜夜躁狠狠躁天天躁| 水蜜桃什么品种好| 国产蜜桃级精品一区二区三区| cao死你这个sao货| 校园春色视频在线观看| 手机成人av网站| 国产不卡一卡二| 日日干狠狠操夜夜爽| 久久人妻av系列| 波多野结衣av一区二区av| 香蕉国产在线看| 在线十欧美十亚洲十日本专区| 亚洲熟妇熟女久久| 国产成人啪精品午夜网站| 精品电影一区二区在线| 久久精品91无色码中文字幕| 男男h啪啪无遮挡| 国产高清激情床上av| 精品国产一区二区久久| av免费在线观看网站| 精品久久久久久久毛片微露脸| 美女扒开内裤让男人捅视频| 亚洲一区中文字幕在线| 少妇被粗大的猛进出69影院| 午夜精品在线福利| 岛国在线观看网站| 成人免费观看视频高清| 精品久久久久久成人av| 嫩草影视91久久| 在线观看午夜福利视频| 欧美乱码精品一区二区三区| 亚洲久久久国产精品| 熟女少妇亚洲综合色aaa.| 亚洲国产欧美日韩在线播放| 亚洲国产欧美一区二区综合| 女人爽到高潮嗷嗷叫在线视频| 99精品欧美一区二区三区四区| 国产成人欧美在线观看| 12—13女人毛片做爰片一| 日韩av在线大香蕉| 一区二区三区国产精品乱码| 久久天堂一区二区三区四区| 可以免费在线观看a视频的电影网站| 老汉色∧v一级毛片| 免费搜索国产男女视频| 精品国产国语对白av| 日本精品一区二区三区蜜桃| 亚洲avbb在线观看| 神马国产精品三级电影在线观看 | 亚洲欧美精品综合一区二区三区| 国产亚洲欧美在线一区二区| 久久久久精品国产欧美久久久| 亚洲成人精品中文字幕电影 | 黄网站色视频无遮挡免费观看| 国产有黄有色有爽视频| 亚洲精品av麻豆狂野| 亚洲视频免费观看视频| 可以免费在线观看a视频的电影网站| 亚洲国产欧美日韩在线播放| 一边摸一边抽搐一进一小说| 日韩免费高清中文字幕av| 免费高清视频大片| 法律面前人人平等表现在哪些方面| 波多野结衣高清无吗| 18禁裸乳无遮挡免费网站照片 | svipshipincom国产片| 久久久国产欧美日韩av| 69精品国产乱码久久久| 欧美精品亚洲一区二区| 精品国产乱码久久久久久男人| 可以在线观看毛片的网站| 亚洲片人在线观看| 亚洲成国产人片在线观看| 国产精品综合久久久久久久免费 | 99久久精品国产亚洲精品| 超碰97精品在线观看| 免费人成视频x8x8入口观看| 亚洲欧美一区二区三区久久| 亚洲精品国产精品久久久不卡| xxx96com| 国产精品国产av在线观看| 可以免费在线观看a视频的电影网站| 久久人人97超碰香蕉20202| 国产精品偷伦视频观看了| 国产精品99久久99久久久不卡| 叶爱在线成人免费视频播放| 久久国产亚洲av麻豆专区| 亚洲成av片中文字幕在线观看| 在线播放国产精品三级| 国产真人三级小视频在线观看| 精品国产一区二区久久| 欧美日韩亚洲高清精品| 嫁个100分男人电影在线观看| 久久久久国产精品人妻aⅴ院| 一个人免费在线观看的高清视频| 亚洲精品国产一区二区精华液| 亚洲精品国产区一区二| 一夜夜www| 在线国产一区二区在线| 国产午夜精品久久久久久| 最新美女视频免费是黄的| 久久精品aⅴ一区二区三区四区| 亚洲第一欧美日韩一区二区三区| 日本五十路高清| tocl精华| 9热在线视频观看99| 国产精品二区激情视频| 免费在线观看日本一区| 精品国产美女av久久久久小说| 老司机福利观看| 新久久久久国产一级毛片| 国产成人一区二区三区免费视频网站| 亚洲一区二区三区色噜噜 | 亚洲国产精品合色在线| 叶爱在线成人免费视频播放| 国产成人精品久久二区二区免费| 国产一区二区激情短视频| 久久久久久久久免费视频了| 少妇被粗大的猛进出69影院| 欧美中文综合在线视频| 男人舔女人下体高潮全视频| 精品午夜福利视频在线观看一区| 欧美日本亚洲视频在线播放| 久久久精品国产亚洲av高清涩受| 欧美亚洲日本最大视频资源| 青草久久国产| 国产亚洲av高清不卡| 亚洲精华国产精华精| 一边摸一边抽搐一进一出视频| 香蕉丝袜av| 色播在线永久视频| 美女高潮到喷水免费观看| 黄色怎么调成土黄色| 欧美黑人欧美精品刺激| 天堂影院成人在线观看| 亚洲视频免费观看视频| 高清在线国产一区| 丰满的人妻完整版| 在线播放国产精品三级| 另类亚洲欧美激情| 国产高清videossex| 亚洲少妇的诱惑av| 午夜精品在线福利| 大型av网站在线播放| 中文字幕色久视频| 黄网站色视频无遮挡免费观看| 久久精品人人爽人人爽视色| 高清在线国产一区| ponron亚洲| 黑丝袜美女国产一区| 夫妻午夜视频| 亚洲国产看品久久| 自拍欧美九色日韩亚洲蝌蚪91| 久久久国产欧美日韩av| 少妇的丰满在线观看| 欧美av亚洲av综合av国产av| 欧美中文日本在线观看视频| 日韩成人在线观看一区二区三区| 欧美中文日本在线观看视频| 搡老岳熟女国产| 精品日产1卡2卡| 亚洲国产看品久久| 多毛熟女@视频| 久久精品aⅴ一区二区三区四区| 午夜久久久在线观看| 两个人免费观看高清视频| 日韩欧美在线二视频| avwww免费| 国产精品永久免费网站| 亚洲激情在线av| 级片在线观看| 欧美黄色片欧美黄色片| 黄色怎么调成土黄色| 99在线人妻在线中文字幕| 人人澡人人妻人| 无遮挡黄片免费观看| av福利片在线| 精品国产亚洲在线| 91国产中文字幕| 久久婷婷成人综合色麻豆| 亚洲色图 男人天堂 中文字幕| 每晚都被弄得嗷嗷叫到高潮| 日本一区二区免费在线视频| av天堂在线播放| 欧美中文日本在线观看视频| 日韩 欧美 亚洲 中文字幕| 日韩欧美在线二视频| 他把我摸到了高潮在线观看| 欧美久久黑人一区二区| 亚洲熟妇熟女久久| 久久久久久久久中文| 国产精品av久久久久免费| 日韩免费高清中文字幕av| 久久精品国产亚洲av高清一级| 美女国产高潮福利片在线看| 麻豆av在线久日| 日本黄色视频三级网站网址| 亚洲一卡2卡3卡4卡5卡精品中文| 国产深夜福利视频在线观看| 一级作爱视频免费观看| 亚洲一区中文字幕在线| 亚洲性夜色夜夜综合| 国产成+人综合+亚洲专区| 满18在线观看网站| 黄色成人免费大全| 欧美精品一区二区免费开放| 老熟妇仑乱视频hdxx| 很黄的视频免费| 亚洲国产精品sss在线观看 | 麻豆成人av在线观看| 久久久久九九精品影院| 成人av一区二区三区在线看| 国产欧美日韩一区二区三| 亚洲欧洲精品一区二区精品久久久| 午夜免费观看网址| 中文字幕av电影在线播放| 黑丝袜美女国产一区| 老司机亚洲免费影院| 免费看十八禁软件| 国产欧美日韩一区二区三区在线| 女人被狂操c到高潮| 男人舔女人下体高潮全视频| 国产精品一区二区免费欧美| 90打野战视频偷拍视频| 国产精品国产高清国产av| 久久精品成人免费网站| 日本免费一区二区三区高清不卡 | 国产亚洲精品一区二区www| 亚洲一码二码三码区别大吗| 一边摸一边做爽爽视频免费| 黄色片一级片一级黄色片| 天天躁狠狠躁夜夜躁狠狠躁| 天堂影院成人在线观看| 女性生殖器流出的白浆| 在线十欧美十亚洲十日本专区| 在线观看午夜福利视频| 欧美日韩亚洲高清精品| 中文字幕色久视频| 欧美激情 高清一区二区三区| 日韩大尺度精品在线看网址 | 可以免费在线观看a视频的电影网站| 国产不卡一卡二| 丰满迷人的少妇在线观看| 国产三级黄色录像| 亚洲国产精品合色在线| 欧美成狂野欧美在线观看| www.自偷自拍.com| 国产精品久久电影中文字幕| 人人妻人人爽人人添夜夜欢视频| 欧美日韩国产mv在线观看视频| 亚洲七黄色美女视频| 亚洲中文字幕日韩| 日本vs欧美在线观看视频| avwww免费| 怎么达到女性高潮| 叶爱在线成人免费视频播放| 麻豆成人av在线观看| 欧美日韩一级在线毛片| 国产三级在线视频| 操美女的视频在线观看| 精品人妻1区二区| 国产97色在线日韩免费| 日韩欧美一区二区三区在线观看| 啦啦啦在线免费观看视频4| 日本免费a在线| 在线观看66精品国产| 精品熟女少妇八av免费久了| 巨乳人妻的诱惑在线观看| 久久热在线av| 丝袜美腿诱惑在线| 亚洲人成网站在线播放欧美日韩| 日本wwww免费看| 一区福利在线观看| 免费在线观看日本一区| 狠狠狠狠99中文字幕| 天堂俺去俺来也www色官网| x7x7x7水蜜桃| 日韩欧美一区二区三区在线观看| 麻豆一二三区av精品| 国产成人精品久久二区二区91| 妹子高潮喷水视频| 女人精品久久久久毛片| 宅男免费午夜| 亚洲av成人一区二区三| 1024视频免费在线观看| 日韩欧美国产一区二区入口| 久久伊人香网站| 亚洲精品美女久久av网站| 国产伦一二天堂av在线观看| 一区二区三区精品91| 免费看a级黄色片| 国产高清videossex| 午夜福利免费观看在线| 国产精品免费视频内射| 国产成人一区二区三区免费视频网站| 一区二区日韩欧美中文字幕| 在线观看免费高清a一片| 免费人成视频x8x8入口观看| 午夜精品国产一区二区电影| 国产精品久久电影中文字幕| 国产一区二区激情短视频| 日韩精品中文字幕看吧| 午夜精品久久久久久毛片777| 精品久久久久久,| xxx96com| 男女之事视频高清在线观看| 欧美亚洲日本最大视频资源| 纯流量卡能插随身wifi吗| 老司机午夜福利在线观看视频| 亚洲av美国av| 国产免费男女视频| 久久性视频一级片| 国产精品久久久av美女十八| 成人影院久久| 久久中文字幕一级| 国产99白浆流出| 老司机午夜福利在线观看视频| 天堂动漫精品| 色老头精品视频在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品999在线| 夜夜躁狠狠躁天天躁| 久久精品国产亚洲av高清一级| 老汉色∧v一级毛片| 精品日产1卡2卡| 在线天堂中文资源库| 一边摸一边抽搐一进一出视频| 精品熟女少妇八av免费久了| 一区二区日韩欧美中文字幕| 亚洲国产欧美网| 丝袜人妻中文字幕| 自线自在国产av| 伦理电影免费视频| 国产有黄有色有爽视频| 亚洲精品中文字幕在线视频| av网站在线播放免费| 亚洲中文av在线| 欧美日韩精品网址| 精品国产一区二区久久| av视频免费观看在线观看| 久久久久久久久免费视频了| 久久精品人人爽人人爽视色| 国产精品久久久久成人av| 18美女黄网站色大片免费观看| 国产深夜福利视频在线观看| 制服人妻中文乱码| 国产亚洲欧美在线一区二区| netflix在线观看网站| 天堂中文最新版在线下载| 老司机深夜福利视频在线观看| 这个男人来自地球电影免费观看| 久久人人精品亚洲av| 高潮久久久久久久久久久不卡| 亚洲欧美一区二区三区黑人| 女同久久另类99精品国产91| 一区在线观看完整版| 亚洲精品一区av在线观看| www.熟女人妻精品国产| 亚洲av美国av| 成人精品一区二区免费| 免费搜索国产男女视频| 亚洲伊人色综图| xxx96com| 69av精品久久久久久| 日韩欧美国产一区二区入口| 欧美久久黑人一区二区| 国产成人精品无人区| 波多野结衣高清无吗| 9色porny在线观看| cao死你这个sao货| 婷婷精品国产亚洲av在线| 91老司机精品| 夜夜爽天天搞| 国产亚洲精品第一综合不卡| 精品久久久久久久久久免费视频 | 中文字幕高清在线视频| 91国产中文字幕| 日韩精品免费视频一区二区三区| 精品国产国语对白av| 真人一进一出gif抽搐免费| 国产主播在线观看一区二区| 女人爽到高潮嗷嗷叫在线视频| 麻豆国产av国片精品| av视频免费观看在线观看| 成人永久免费在线观看视频| 欧美日韩一级在线毛片| 亚洲成人国产一区在线观看| 久久久久国产精品人妻aⅴ院| 亚洲人成伊人成综合网2020| 中文字幕人妻丝袜一区二区| 国产精品99久久99久久久不卡| 久久久久久亚洲精品国产蜜桃av| 啦啦啦 在线观看视频| 在线观看免费视频日本深夜| 欧美一级毛片孕妇| 国产成人精品在线电影| 老汉色∧v一级毛片|