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

    背水面齒角度對齒形迷宮流道灌水器水力特性及抗堵性能的影響

    2025-04-28 00:00:00鄭健陳婭趙雨露張國維王燕

    摘要: 為優(yōu)化齒形迷宮流道灌水器的流道結(jié)構(gòu),應(yīng)用CFD數(shù)值模擬方法,對背水面齒角度分別為0°,10°,19°,30°和45°的齒形迷宮流道灌水器進(jìn)行固液兩相流數(shù)值模擬計(jì)算,研究不同背水面齒角度流道內(nèi)流速、湍動能、渦旋分布規(guī)律及顆粒運(yùn)動軌跡等.結(jié)果表明:齒形迷宮流道灌水器流態(tài)指數(shù)和流量系數(shù)均與背水面齒角度呈正相關(guān)關(guān)系,背水面齒角度為0°時(shí),水力性能最佳;隨著背水面齒角度的增大,灌水器流道內(nèi)流速整體呈減小趨勢;流道內(nèi)湍動能隨背水面齒角度的增大而增強(qiáng);通過Ω渦識別法發(fā)現(xiàn)隨著背水面齒角度的增大,流道內(nèi)渦旋運(yùn)動區(qū)域面積減小且強(qiáng)度減弱;顆粒在5種背水面齒角度流道內(nèi)最長停留時(shí)間分別為1.92,3.77,4.37,0.87,0.82 s,背水面齒角度為45° 時(shí),灌水器抗堵塞性能相對較優(yōu).模擬結(jié)果證明適當(dāng)改變背水面齒角度,采用不對稱流道形式,有利于提高灌水器的性能.

    關(guān)鍵詞: 迷宮灌水器;背水面;齒角度;水力特性;抗堵性能

    中圖分類號: S275.6 文獻(xiàn)標(biāo)志碼: A 文章編號: 1674-8530(2025)04-0357-07

    DOI:10.3969/j.issn.1674-8530.23.0217

    鄭健, 陳婭, 趙雨露,等.背水面齒角度對齒形迷宮流道灌水器水力特性及抗堵性能的影響[J]. 排灌機(jī)械工程學(xué)報(bào),2025,43(4):357-363.

    ZHENG Jian, CHEN Ya, ZHAO Yulu, et al. Influence of tooth angle of back water surface on hydraulic characteristics and anti-clogging performance of dental labyrinth emitter[J]. Journal of drainage and irrigation machinery engineering (JDIME), 2025, 43(4): 357-363. (in Chinese)

    Influence of tooth angle of back water surface on

    hydraulic characteristics and anti-clogging

    performance of dental labyrinth emitter

    ZHENG Jian1,2*, CHEN Ya1,2, ZHAO Yulu1,2, ZHANG Guowei3," WANG Yan1

    (1. College of Energy and Power Engineering, Lanzhou University of Technology, Lanzhou, Gansu 730050, China;" 2. Key Laboratory of Multi-supply System with Solar Energy and Biomass, Lanzhou, Gansu 730050, China; 3. Dahua Water Saving Technology Co., Ltd., Gulang, Gansu 733101, China)

    Abstract: In order to optimize the flow channel structure of toothed labyrinth flow channel emitter, CFD numerical simulation method was used to simulate the solid-liquid two-phase flow numerical calculations on toothed labyrinth emitters with back surface tooth angles of 0°, 10°, 19°, 30° and 45°, and the differences of flow velocity, turbulent kinetic energy, vortex distribution and particle trajectory in the flow channel under different tooth angles of back water surface were studied. The results show that the flow index and the discharge coefficient of the labyrinth emitter are positively correlated with the tooth angle of the back water surface, the hydraulic performance is the best when the tooth angle of the back water surface is 0°. As the" tooth angle of the back water surface increases, the flow velocity in the emitter flow channel decreases as a whole, and the turbulent kinetic energy intensity increases. Through the Ω vortex identification, it is found that with the increase of the tooth angle of the back water surface, the area of the vortex motion area in the flow channel decreases and the intensity decreases. The longest particle residence time in the flow channel of the five kinds of back water surface tooth angles is 1.92, 3.77, 4.37, 0.87, 0.82 s, respectively. The emitter′s anti-clogging performance is optimal when the tooth angle of the back water surface is 45°. The simulation results show that the performance of the emitter can be improved by changing the tooth angle of the back water surface and adopting the form of an asymmetric flow channel.

    Key words: labyrinth emitter;back water surface;tooth angle;hydraulic characteristics;anti-clogging performance

    流道結(jié)構(gòu)是影響灌水器水力特性和抗堵性能的重要因素[1-2],合適的流道尺寸和形狀可降低其堵塞風(fēng)險(xiǎn).迷宮流道灌水器結(jié)構(gòu)相對簡單、水力性能高效、成本低,是目前使用最廣泛的灌水器[3].目前研究主要集中于迷宮流道灌水器結(jié)構(gòu)參數(shù)的優(yōu)化和組合,如喻黎明等[4]采用水沙運(yùn)動CFD-DEM 耦合數(shù)值模擬分析了齒角度、上底寬、齒高、齒尖參差量、流道寬等結(jié)構(gòu)參數(shù)對迷宮灌水器的影響,提出齒角度、流道寬是影響灌水器流道內(nèi)水流特性的主要流道結(jié)構(gòu)參數(shù);張傳杰等[5]發(fā)現(xiàn)齒形迷宮流道灌水器的低速區(qū)和旋渦區(qū)主要分布于齒尖的背水面,但未進(jìn)行深入研究.目前,根據(jù)流道內(nèi)低速區(qū)和旋渦區(qū)分布情況來優(yōu)化灌水器的水力特性和抗堵塞性能的研究較少.

    基于此,文中應(yīng)用CFD數(shù)值模擬方法,以達(dá)華節(jié)水科技股份有限公司生產(chǎn)的3 L/h的齒形迷宮流道灌水器E0作為基本流道形式,分析背水面齒角度為0°, 10°, 19° (E0),30° 和45° 時(shí)流道內(nèi)流速、湍動能、渦旋分布規(guī)律及顆粒運(yùn)動軌跡等,研究背水面齒角度對灌水器水力特性及抗堵塞性能的影響.以期為齒形迷宮流道灌水器的設(shè)計(jì)和優(yōu)化提供理論依據(jù).

    1 材料與方法

    1.1 物理模型

    采用達(dá)華節(jié)水科技股份有限公司生產(chǎn)的3 L/h的齒形迷宮流道灌水器E0作為基本流道,并以E0結(jié)構(gòu)形式為基礎(chǔ),僅將背水面齒角度α改變?yōu)?°,10°,30°和45°,形成其余4種齒形迷宮流道灌水器.圖1為灌水器結(jié)構(gòu)參數(shù)示意圖.流道截面為0.8 mm×1.0 mm,流道深度d為0.8 mm,流道高度H為 2.0 mm,齒高h(yuǎn)為1.0 mm,流道單元長度S為2.22 mm,齒角度θ為38°,背水面齒角度α為19°、出水口半徑r為0.9 mm,L1, L2表示流道深1/2位置(d=0.4 mm)截面上截取的流道內(nèi)縱向和橫向的2條直線,便于后期觀察流道中該處流速和湍動能隨背水面齒角度的變化規(guī)律.

    1.2 數(shù)值模擬方法

    1.2.1 數(shù)值模型與邊界條件

    應(yīng)用SolidWorks軟件對5種不同背水面齒角度的灌水器進(jìn)行建模.在ANSYS的Workbench平臺的Fluid Flow(Fluent)組件中進(jìn)行模擬,在Meshing中對三維物理模型進(jìn)行六面體非結(jié)構(gòu)網(wǎng)格劃分,采用標(biāo)準(zhǔn)k-ε湍流模型對不同網(wǎng)格數(shù)量的灌水器進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證.經(jīng)驗(yàn)證5種不同背水面齒角度灌水器的網(wǎng)格大小為0.06 mm時(shí),灌水器出水流量計(jì)算結(jié)果前后差異最大為0.3%,同時(shí)考慮到計(jì)算精度與計(jì)算效率等因素,最終確定5種灌水器網(wǎng)格大小為0.06 mm,網(wǎng)格數(shù)量范圍為43~46萬.

    用標(biāo)準(zhǔn)k-ε湍流模型對單相流進(jìn)行數(shù)值模擬,數(shù)值計(jì)算采用有限體積法離散控制方程,壓力項(xiàng)等參數(shù)采用二階迎風(fēng)格式,采用SIMPLE算法對速度和壓力的耦合求解,流道壁面采用標(biāo)準(zhǔn)壁面函數(shù)法,收斂精度設(shè)為10-4.將液相視作黏性不可壓縮流體,即密度為常數(shù),湍流模型結(jié)果驗(yàn)證時(shí)液相設(shè)為水,單相流模擬時(shí)液相設(shè)為水和沼液混合液體,其物理性質(zhì)參考文獻(xiàn)[6]中沼液與水的配比為1∶4時(shí)的參數(shù),密度為992.8 kg/m3,黏度為0.001 06 Pa·s.考慮重力作用,忽略表面張力影響.將連續(xù)相計(jì)算收斂后,采用離散相模型,加入顆粒相進(jìn)行模擬.考慮到顆粒受湍流影響的隨機(jī)性和流體湍流脈動會引起顆粒擴(kuò)散,采用隨機(jī)軌道模型和雙向耦合計(jì)算,將最大計(jì)算步數(shù)設(shè)置為50 000,保證顆粒全部逸出流道.顆粒入射方式設(shè)置為surface,入射顆粒流將從流道進(jìn)水口面以與進(jìn)口水流相同的速度垂直射入灌水器流道,顆粒體積分?jǐn)?shù)設(shè)為1%,壁面反彈系數(shù)設(shè)為0.8,其中,顆粒成分參考文獻(xiàn)[7]設(shè)為SiO2,密度設(shè)為2 200 kg/m3,粒徑為25 μm.

    1.2.2 湍流模型驗(yàn)證

    為驗(yàn)證文中所選湍流模型的合理性,對原型灌水器進(jìn)行物理水力性能試驗(yàn)測試.試驗(yàn)用水為自來水,每次測試時(shí)間為3 min,設(shè)置3條滴灌帶,每條隨機(jī)取3個(gè)灌水器測量流量,采用稱量法進(jìn)行流量測定,取平均值作為每個(gè)壓力條件下的流量.數(shù)值模擬采用標(biāo)準(zhǔn)k-ε湍流模型,進(jìn)口壓力設(shè)為0.03~0.15 MPa,出口壓力設(shè)為0,液相設(shè)為水.圖2為壓力在0.03~0.15 MPa條件下灌水器流量的試驗(yàn)和模擬結(jié)果,圖中Q為流量,p為壓力.試驗(yàn)與模擬相對應(yīng)的水力性能曲線變化一致且誤差為2%~7%,說明文中采用標(biāo)準(zhǔn)k-ε湍流模型進(jìn)行數(shù)值模擬是可行的.

    1.3 壓力-流量關(guān)系式

    灌水器流道的水力性能可由流態(tài)指數(shù)和流量系數(shù)反映,流態(tài)指數(shù)越小,相應(yīng)的灌水器水力性能越好.可根據(jù)不同壓力下得到的流量應(yīng)用公式(1)在SPSS中進(jìn)行冪函數(shù)擬合得到灌水器的流量系數(shù)與流態(tài)指數(shù).

    Q=ahbp,(1)

    式中:a為流量系數(shù);hp為壓力水頭;b為流態(tài)指數(shù).

    1.4 Ω渦識別法

    Ω渦識別法基于速度梯度張量V的分解,將其分解為2部分:對稱張量A和反對稱張量B.這2部分分別表示流場某一點(diǎn)的變形和旋轉(zhuǎn),使得該方法能夠同時(shí)辨識出流場中的強(qiáng)渦和弱渦,并通過Ω值反映剪切流與旋轉(zhuǎn)流的存在[8].Ω方法的物理意義是旋轉(zhuǎn)部分渦量大小占總渦量大小的比例,即

    Ω=‖B‖2F‖A‖2F+‖B‖2F+φ,(2)

    式中:φ為一個(gè)極小的正數(shù),旨在避免分母為0的情況.通過式(2),得知Ω的取值為0 ≤Ωlt;1.00.當(dāng)閾值超過0.50時(shí),說明反對稱張量B相對于對稱張量A有優(yōu)勢,所以可選擇略大于0.50的閾值作為渦識別的依據(jù)[9].LIU等[10]建議在實(shí)際應(yīng)用中選用Ω=0.52作為固定閾值,當(dāng)Ω低于0.52時(shí)表示流道內(nèi)運(yùn)動主要為剪切運(yùn)動,而當(dāng)Ω超過0.52時(shí)表示運(yùn)動以旋轉(zhuǎn)為主.

    1.5 數(shù)據(jù)處理方法

    文中采用CFD-Post,Tecplot 360,SPSS及Origin 2018軟件對模擬結(jié)果及數(shù)據(jù)進(jìn)行處理分析,通過F檢驗(yàn)法檢驗(yàn)回歸方程的顯著性,t檢驗(yàn)法檢驗(yàn)回歸系數(shù)的顯著性.

    2 結(jié)果與分析

    2.1 背水面齒角度變化對灌水器水力特性的影響

    根據(jù)5種流道形式在壓力為0.03~0.15 MPa條件下的流量運(yùn)用式(1)可擬合得出其流量系數(shù)和流態(tài)指數(shù),如表1所示.由表可以看出5種流道形式下的決定系數(shù)R2值均為1.00,回歸方程和回歸系數(shù)的顯著性均小于0.001,說明回歸方程擬合度高且擬合結(jié)果可靠.流態(tài)指數(shù)x與背水面齒角度成正相關(guān)關(guān)系,表明隨著背水面齒角度增大,流量對壓力的敏感程度增強(qiáng),水力特性逐漸下降.同時(shí)流量系數(shù)也隨背水面齒角度的增加而增大,且增幅比流態(tài)指數(shù)的明顯.

    2.2 背水面齒角度變化對流道內(nèi)流場的影響

    圖3為進(jìn)口壓力為0.08 MPa,流道深1/2位置截面的流速分布圖,v為流速.由圖可知,5種不同流道結(jié)構(gòu)形式灌水器的流道流場分布規(guī)律相似,齒尖迎水面流速較大,背水面流速較小,可分為以下4個(gè)區(qū)域:A高速區(qū)、B中速區(qū)、C渦旋區(qū)以及D低速區(qū),A高速區(qū)與B中速區(qū)組成了流道的主流區(qū)域,C渦旋區(qū)和D低速區(qū)為非主流區(qū)域.

    背水面齒角度α為0°,10°,19°,30°和45°時(shí)流道深1/2位置截面A高速區(qū)流速分別是2.60~4.25,2.57~4.17,2.55~4.14,2.56~3.70,2.58~3.94 m/s,B中速區(qū)流速分別是1.18~2.48,1.11~2.44,1.08~2.40,0.84~2.42,0.92~2.58 m/s,C渦旋區(qū)流速分別是0.10~0.60,0.16~0.91,0.12~0.61,0.14~0.55,0.15~0.65 m/s,D低速區(qū)流速分別是0.15~1.81,0.15~1.67,0.22~1.50,0.22~1.39,0.09~1.35 m/s.當(dāng)背水面齒角度增大時(shí),A高速區(qū)最大流速呈減小趨勢,最小流速變化很?。籅中速區(qū)最大流速隨著背水面齒角度增大呈先減小再增大的趨勢,最小流速呈減小趨勢;由于C渦旋區(qū)位于非主流區(qū)域中心位置,所以背水面齒角度的變化主要影響的是渦旋區(qū)的位置和大小,對渦旋區(qū)的速度影響不大;D低速區(qū)的最大和最小流速均隨背水面齒角度的增大呈減小趨勢.

    圖4為5種不同背水面齒角度流道,流道深1/2位置截面縱向L1和橫向L2這2條直線上速度隨背水面齒角度的變化,圖中y為流道深1/2位置截面上縱向L1直線從上到下y方向位置變化;x為流道深1/2位置截面上橫向L2直線從左到右x方向位置變化.可以看出流道橫向和縱向上速度隨背水面齒角度的增加而減小,反映了流道內(nèi)整體流速隨著背水面齒角度的增大而減小的現(xiàn)象,與上述結(jié)論相同.

    2.3 背水面齒角度變化對流道湍動能的影響

    湍動能是分析灌水器流道內(nèi)能量耗散機(jī)理的一個(gè)重要物理量.湍動能越大表明在該區(qū)域內(nèi)湍流從時(shí)均流提取的能量越多,即水流能量轉(zhuǎn)換成湍動能的比例越大,水流也越容易維持湍流的流動狀態(tài)[2],紊流強(qiáng)度就越大,同時(shí)湍流的強(qiáng)剪切力能促進(jìn)顆粒物質(zhì)團(tuán)簇,增加了顆粒物質(zhì)從流道流出的概率,從而降低了流道堵塞的風(fēng)險(xiǎn).

    圖5為進(jìn)口壓力為0.08 MPa時(shí),流道深1/2位置截面的湍動能云圖,圖中Ke為湍動能.由圖可以看出,5種流道形式的湍動能的分布規(guī)律都相同,并且有明顯的湍流梯度,湍動能較高的區(qū)域主要分布于齒尖的迎水面,即流道的主流區(qū)域的中速區(qū),而非主流區(qū)域湍動能較低.隨著背水面?zhèn)三X角度增大,高湍動能區(qū)域面積增加,湍流穩(wěn)定性也隨之變強(qiáng),特別是當(dāng)背水面齒角度為45°時(shí)低湍動能區(qū)域幾乎消失.

    圖6為5種不同背水面齒角度流道深1/2位置截面縱向L1和橫向L2這2條直線上湍動能隨背水面齒角度的變化,由圖可以看出隨著齒角度的增大, L1和L2的峰值增加,說明湍動能強(qiáng)度與背水面齒角度成正比.

    2.4 背水面齒角度變化對流道內(nèi)渦旋的影響

    渦流特性對灌水器的水力特性和抗堵塞性能有著關(guān)鍵作用[11].為進(jìn)一步研究不同背水面齒角度對流道內(nèi)渦流特性的影響,采用Ω渦識別方法對流道內(nèi)旋渦結(jié)構(gòu)進(jìn)行捕捉,以更準(zhǔn)確地探究背水面齒角度對渦旋區(qū)域的影響.

    圖7為5種不同背水面齒角度在工作壓力為0.08 MPa時(shí)流道深1/2位置截面處Ω≥0.52渦旋強(qiáng)度云圖,同時(shí)也展現(xiàn)了速度矢量.由圖7可知流道內(nèi)Ω≥0.52時(shí)為渦旋運(yùn)動區(qū),否則為剪切運(yùn)動區(qū).渦旋運(yùn)動區(qū)主要有2個(gè)區(qū)域:一是Ω值較高的強(qiáng)渦區(qū),主要分布于C渦旋區(qū)與B中速區(qū)、D低速區(qū)相鄰的回流區(qū)域;另一個(gè)是Ω值較低的弱渦區(qū).隨著背水面齒角度的增大,非主流區(qū)域面積減少,流道內(nèi)渦旋運(yùn)動區(qū)域面積減小且強(qiáng)度減弱.

    2.5 離散相顆粒運(yùn)動軌跡分布情況

    借助 CFD-Post 查看不同形式灌水器流道內(nèi)固體顆粒濃度的分布情況,圖8為顆粒折減系數(shù)為3.5,體積分?jǐn)?shù)為1%時(shí),固體顆粒在5種形式灌水器內(nèi)的運(yùn)動軌跡,圖中T為顆粒停留的時(shí)間.

    如圖8所示,大部分固體顆??梢噪S灌水器流道主流區(qū)水流流出,少數(shù)固體顆粒會在齒尖背水面區(qū)域做無規(guī)則旋轉(zhuǎn)運(yùn)動,齒尖背水面區(qū)域?yàn)闇u旋區(qū),該區(qū)域流速較低,存在渦旋,導(dǎo)致顆粒在此區(qū)域做旋轉(zhuǎn)運(yùn)動不易脫離,增加了顆粒在流道內(nèi)的運(yùn)移時(shí)間,也增大了灌水器發(fā)生堵塞的風(fēng)險(xiǎn).顆粒在背水面齒角度為0°,10°,19°,30°和45°流道內(nèi)的最長停留時(shí)間分別為1.92,3.77,4.37,0.87和0.82 s.當(dāng)α=45°時(shí),顆粒在灌水器內(nèi)停留時(shí)間最短,灌水器抗堵塞性能相對較優(yōu),與E0(α=19°)相比,增大或減小背水面齒角度均有利于降低灌水器堵塞概率.

    3 討 論

    文中原型灌水器E0為對稱流道結(jié)構(gòu)形式,其余4種為不對稱流道形式.當(dāng)背水面齒角度減小為10°和0°時(shí),流態(tài)指數(shù)減小,水力性能提升;流道內(nèi)最大流速增大,渦旋運(yùn)動區(qū)域增大,與E0(α=19°)相比,顆粒在流道內(nèi)滯留的最長時(shí)間由原來的4.37 s變?yōu)榱?.77 s和1.92 s.這主要是因?yàn)楸乘纨X角度減小時(shí),流道內(nèi)背水面面積增加,擴(kuò)展了非主流區(qū)域水流流動范圍,使得渦旋區(qū)發(fā)展充分,提升了灌水器流道的水力性能和自清洗能力[12-13].當(dāng)背水面齒角度增大至30° 和45° 時(shí),流道非主流區(qū)域面積占比減小,且湍動能增強(qiáng),顆粒在流道內(nèi)滯留的最長時(shí)間分別是0.87 s和0.82 s,較E0明顯下降,灌水器抗堵塞性能大大提升.這可能是由于背水面齒角度的增大,減小了低速區(qū)與渦旋區(qū)范圍,湍動能也隨背水面齒角度的增加而增強(qiáng),在較大的湍流剪切力下,黏附在壁上的堵塞物質(zhì)更加無序,增加了堵塞物質(zhì)流出流道的可能性,同時(shí)渦旋流隨背水面齒角度的增大而只留下主流區(qū)部分,主流區(qū)速度較大,顆粒會被主流區(qū)水流帶出流道[14].

    綜上分析,可以得出當(dāng)背水面齒角度減小至0°時(shí),渦旋區(qū)發(fā)展最為充分,有助于提升灌水器流道的自清洗能力和增強(qiáng)流道內(nèi)的能量消耗,水力性能最佳;當(dāng)背水面齒角度增大至45°時(shí),流道內(nèi)湍動能增強(qiáng),顆粒在流道內(nèi)滯留的最長時(shí)間最短,抗堵塞性能最佳;通過適當(dāng)增大或減小背水面齒角度,采用不對稱流道形式,有利于提高灌水器性能.

    4 結(jié) 論

    1) 齒形迷宮流道灌水器流態(tài)指數(shù)和流量系數(shù)均隨背水面齒角度增大而增大,但流態(tài)指數(shù)變化幅度較小,背水面齒角度越小,水力性能越好.

    2) 5種不同流道結(jié)構(gòu)形式灌水器的流道流場分布規(guī)律相似,齒尖迎水面流速較大,背水面流速較?。浑S著背水面齒角度的增大,流道內(nèi)流場流速整體呈減小趨勢,湍動能呈增大趨勢;通過Ω渦識別法發(fā)現(xiàn)隨著背水面齒角度的增大,流道內(nèi)渦旋運(yùn)動區(qū)域面積減小且渦旋運(yùn)動強(qiáng)度降低.

    3) 適當(dāng)增大或減小背水面齒角度,采用不對稱流道形式,有利于提高灌水器性能.

    參考文獻(xiàn)(References)

    [1] 王惠,王世雷,王文娥,等.邊縫式迷宮滴灌帶灌水器流道內(nèi)泥沙沉積影響因素分析[J].農(nóng)業(yè)工程學(xué)報(bào),2023,39(24):92-99.

    WANG Hui, WANG Shilei, WANG Wen′e, et al. In-fluence of the emitter with tooth-shape labyrinth flow channel on sediment deposition[J]. Transactions of the CSAE, 2023, 39(24): 92-99.(in Chinese)

    [2] YANG B, WANG J D, ZHANG Y Q, et al. Anti-clogging performance optimization for dentiform labyrinth emitters[J]. Irrigation science, 2020, 38: 275-285.

    [3] YU L M, LI N, LIU Q L, et al. Influence of dentation angle of labyrinth channel of drip emitters on hydraulic and anti-clogging performance[J]. Irrigation and drai-nage, 2019, 68(2): 256-267.

    [4] 喻黎明,徐霞,楊啟良,等.滴灌灌水器迷宮流道結(jié)構(gòu)對泥沙運(yùn)動的影響[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2017, 48(2): 255-261.

    YU Liming, XU Xia, YANG Qiliang, et al. Influence of geometrical parameters of labyrinth passage of drip irrigation emitter on sand movement[J]. Transactions of the CSAM, 2017, 48(2): 255-261. (in Chinese)

    [5] 張傳杰,牛勇,劉逍航.齒間角對迷宮灌水器水力特性影響的數(shù)值模擬[J]. 排灌機(jī)械工程學(xué)報(bào), 2022, 40(7): 751-756.

    ZHANG Chuanjie, NIU Yong, LIU Xiaohang. Numeri-cal simulation of influence of tooth angle on hydraulic characteristics of labyrinth emitter[J]. Journal of drainage and irrigation machinery engineering, 2022, 40(7): 751-756. (in Chinese)

    [6] 張恩繼.水—沼液一體穴灌土壤水分及銨態(tài)氮遷移特征試驗(yàn)研究[D].蘭州: 蘭州理工大學(xué), 2016.

    [7] 鄭杰, 陳紅, 孟令劍,等.不同類型灌水器滴頭對沼液抗堵塞性能的影響[J].中國農(nóng)業(yè)科技導(dǎo)報(bào), 2022, 24(7): 108-116.

    ZHENG Jie, CHEN Hong, MENG Lingjian, et al. Effect of different emitter types on anti-clogging perfor-mance of biogas slurry concentrations[J]. Journal of agricultural science and technology, 2022, 24(7): 108-116. (in Chinese)

    [8] 閆龍龍,高波,張寧,等.基于Omega渦識別方法的誘導(dǎo)輪內(nèi)非定??栈鲃臃治觯跩].工程熱物理學(xué)報(bào),2023,44(3):654-659.

    YAN Longlong, GAO Bo, ZHANG Ning, et al. Analysis of unsteady cavitating flow in an inducer based on Omega vortex identification[J]. Journal of engineering thermophysics, 2023,44(3):654-659. (in Chinese)

    [9] 趙建濤,袁建平,裴吉,等. Ω渦識別方法在雙進(jìn)口兩級中開泵中的應(yīng)用[J]. 排灌機(jī)械工程學(xué)報(bào), 2022, 40(5): 447-453.

    ZHAO Jiantao, YUAN Jianping, PEI Ji, et al. Application of Ω vortex identification method in two-stage double-suction centrifugal pump[J]. Journal of drainage and irrigation machinery engineering, 2022, 40(5): 447-453. (in Chinese)

    [10] LIU C Q, WANG Y Q, YANG Y, et al. New Omega vortex identification method[J]. Science China (physics, mechanics amp; astronomy), 2016, 59(8): 62-70.

    [11] WANG C C, LI Z Q, MA J J. Influence of emitter structure on its hydraulic performance based on the vortex[J]. Agriculture, 2021, 11(6): 508-520.

    [12] LI Y F, FENG X Y, HAN X C, et al. Analysis of internal flow characteristics and determination of the core structural parameters of a stellate labyrinth channel irrigation emitter[J]. Biosystems engineering,2023,231:1-19.

    [13] QIU X F, CHEN G F, WANG H T, et al. Vertical optimisation of tooth shape to improve the anti-clogging performance of emitters in drip irrigation systems[J]. Biosystems engineering, 2023, 233:193-203.

    [14] XU T Y, ZHANG L X. Influence and analysis of structure design and optimization on the performance of a pit drip irrigation emitter[J]. Irrigation and drainage, 2020, 69(4): 633-645.

    (責(zé)任編輯 朱漪云)

    收稿日期: 2023-10-26; 修回日期: 2024-03-10; 網(wǎng)絡(luò)出版時(shí)間: 2025-04-01

    網(wǎng)絡(luò)出版地址: https://link.cnki.net/urlid/32.1814.TH.20250331.1720.002

    基金項(xiàng)目: 國家自然科學(xué)基金資助項(xiàng)目(51969012);國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2022YFD1900503);甘肅省教育廳產(chǎn)業(yè)支撐計(jì)劃項(xiàng)目(2021CYZC-27,2021CYZC-33);濟(jì)南市水務(wù)科技項(xiàng)目(JNSWKJ202206);2023年高校教師創(chuàng)新基金資助項(xiàng)目(2023B-431)

    第一作者簡介: 鄭健(1981—),男,甘肅會寧人,教授,博士生導(dǎo)師(通信作者,zhj16822@126.com),主要從事農(nóng)業(yè)水土工程研究.

    第二作者簡介: 陳婭(1998—),女,貴州息烽人,碩士研究生(c18394139271@163.com),主要從事農(nóng)業(yè)水土工程研究.

    午夜老司机福利剧场| 一级作爱视频免费观看| 男人狂女人下面高潮的视频| 少妇的逼水好多| 两性午夜刺激爽爽歪歪视频在线观看| 99国产极品粉嫩在线观看| 男女下面进入的视频免费午夜| 亚洲av日韩精品久久久久久密| 久久天躁狠狠躁夜夜2o2o| 中出人妻视频一区二区| 91狼人影院| 国产一级毛片七仙女欲春2| 久久精品夜夜夜夜夜久久蜜豆| 18禁在线播放成人免费| 18禁黄网站禁片午夜丰满| 国产 一区 欧美 日韩| 两性午夜刺激爽爽歪歪视频在线观看| 18禁黄网站禁片免费观看直播| 欧美性感艳星| xxxwww97欧美| 国产亚洲精品综合一区在线观看| 国产野战对白在线观看| 免费高清视频大片| 精品人妻熟女av久视频| 色综合婷婷激情| ponron亚洲| 香蕉av资源在线| 成人午夜高清在线视频| 亚洲欧美精品综合久久99| 国产伦在线观看视频一区| 久久精品人妻少妇| 免费大片18禁| 亚洲av不卡在线观看| 特级一级黄色大片| 国内揄拍国产精品人妻在线| 日韩欧美精品免费久久 | 欧美黄色片欧美黄色片| 色播亚洲综合网| 欧美xxxx性猛交bbbb| 亚洲在线自拍视频| 欧美三级亚洲精品| www.999成人在线观看| 久久亚洲真实| 99久久成人亚洲精品观看| h日本视频在线播放| 少妇被粗大猛烈的视频| 在线观看午夜福利视频| 国产大屁股一区二区在线视频| av福利片在线观看| 国产精品av视频在线免费观看| 女人十人毛片免费观看3o分钟| 老司机午夜福利在线观看视频| 在线国产一区二区在线| 在线观看av片永久免费下载| 国产精品自产拍在线观看55亚洲| 免费看日本二区| 欧美日韩国产亚洲二区| 亚洲精品成人久久久久久| 黄色视频,在线免费观看| 亚洲人与动物交配视频| 内射极品少妇av片p| 久久国产乱子伦精品免费另类| 99国产综合亚洲精品| 国产亚洲欧美98| 亚洲,欧美,日韩| 精品人妻视频免费看| 免费大片18禁| 在线播放国产精品三级| 999久久久精品免费观看国产| 国产91精品成人一区二区三区| 精品久久久久久久末码| 午夜精品一区二区三区免费看| 亚洲一区二区三区色噜噜| 日本精品一区二区三区蜜桃| 精品无人区乱码1区二区| 欧美高清成人免费视频www| av中文乱码字幕在线| 国产精品精品国产色婷婷| 日本黄色片子视频| 国产精品1区2区在线观看.| 国产蜜桃级精品一区二区三区| 一区二区三区免费毛片| 人妻丰满熟妇av一区二区三区| 国内精品美女久久久久久| 国产av不卡久久| 亚洲人成电影免费在线| 亚洲,欧美,日韩| eeuss影院久久| 俄罗斯特黄特色一大片| 精品99又大又爽又粗少妇毛片 | 97碰自拍视频| 听说在线观看完整版免费高清| 美女高潮的动态| 国产成+人综合+亚洲专区| 久久人人精品亚洲av| 我要搜黄色片| 亚洲第一电影网av| 精品一区二区三区人妻视频| 综合色av麻豆| 99热这里只有精品一区| 永久网站在线| 色综合站精品国产| 国产精品久久电影中文字幕| 亚洲中文字幕日韩| 90打野战视频偷拍视频| 一卡2卡三卡四卡精品乱码亚洲| 国产精品乱码一区二三区的特点| 色吧在线观看| 午夜激情欧美在线| 老司机深夜福利视频在线观看| 亚洲成人中文字幕在线播放| 久久午夜福利片| 观看免费一级毛片| 人人妻,人人澡人人爽秒播| 欧美一区二区亚洲| www.熟女人妻精品国产| 成人国产一区最新在线观看| 岛国在线免费视频观看| 亚洲av成人精品一区久久| 99久久久亚洲精品蜜臀av| 欧美日韩中文字幕国产精品一区二区三区| 91麻豆av在线| 日本成人三级电影网站| 久久久久久久午夜电影| 免费黄网站久久成人精品 | 1024手机看黄色片| 真人做人爱边吃奶动态| 色在线成人网| 午夜老司机福利剧场| 淫妇啪啪啪对白视频| 嫩草影视91久久| 搡老妇女老女人老熟妇| 亚洲人成网站高清观看| 成人一区二区视频在线观看| 此物有八面人人有两片| 色精品久久人妻99蜜桃| 99精品在免费线老司机午夜| 色噜噜av男人的天堂激情| 十八禁人妻一区二区| 中文在线观看免费www的网站| 我的老师免费观看完整版| 亚洲av电影在线进入| 99久国产av精品| 变态另类成人亚洲欧美熟女| 午夜福利免费观看在线| 国产av一区在线观看免费| 亚洲一区二区三区色噜噜| or卡值多少钱| 一本精品99久久精品77| 亚洲精品一区av在线观看| 亚洲欧美日韩卡通动漫| 国产精品野战在线观看| 最新中文字幕久久久久| 淫秽高清视频在线观看| 亚洲av电影在线进入| 欧美日韩福利视频一区二区| 欧美日韩瑟瑟在线播放| АⅤ资源中文在线天堂| 国产精品久久久久久久久免 | 国内毛片毛片毛片毛片毛片| 直男gayav资源| 我的老师免费观看完整版| 极品教师在线免费播放| 国内久久婷婷六月综合欲色啪| 国产一级毛片七仙女欲春2| 嫩草影院精品99| 亚洲男人的天堂狠狠| 国产精品久久久久久久电影| 搡老妇女老女人老熟妇| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲国产精品成人综合色| 久久人妻av系列| 91午夜精品亚洲一区二区三区 | 蜜桃久久精品国产亚洲av| 午夜福利在线观看免费完整高清在 | 狂野欧美白嫩少妇大欣赏| 人人妻人人看人人澡| 日韩中字成人| 狂野欧美白嫩少妇大欣赏| 欧美激情久久久久久爽电影| 精品午夜福利在线看| 97超级碰碰碰精品色视频在线观看| 一个人看的www免费观看视频| 欧美精品国产亚洲| 国产亚洲欧美98| 在线观看免费视频日本深夜| 琪琪午夜伦伦电影理论片6080| 国产精品久久久久久久久免 | 中亚洲国语对白在线视频| 中文字幕久久专区| 草草在线视频免费看| 成人性生交大片免费视频hd| 99久久九九国产精品国产免费| 长腿黑丝高跟| 欧美日韩亚洲国产一区二区在线观看| 亚洲专区中文字幕在线| 亚洲欧美清纯卡通| 婷婷亚洲欧美| 亚洲国产欧美人成| 国产aⅴ精品一区二区三区波| 男人舔奶头视频| 日本免费a在线| 一边摸一边抽搐一进一小说| 99热这里只有是精品50| 一级a爱片免费观看的视频| 国产亚洲av嫩草精品影院| 国产伦在线观看视频一区| 久久精品国产清高在天天线| 如何舔出高潮| 好男人电影高清在线观看| 亚洲精品影视一区二区三区av| 国产综合懂色| 日本与韩国留学比较| 久久久久九九精品影院| 欧美不卡视频在线免费观看| 久久久精品欧美日韩精品| 天堂影院成人在线观看| 亚洲成人中文字幕在线播放| 美女 人体艺术 gogo| 精品福利观看| 听说在线观看完整版免费高清| 最近最新中文字幕大全电影3| 精品人妻偷拍中文字幕| 国产91精品成人一区二区三区| 亚洲精品久久国产高清桃花| 亚洲狠狠婷婷综合久久图片| 国产一区二区三区视频了| 精品人妻偷拍中文字幕| 97人妻精品一区二区三区麻豆| 亚洲av一区综合| 国产中年淑女户外野战色| 级片在线观看| 全区人妻精品视频| 男人舔女人下体高潮全视频| 亚洲精华国产精华精| 男女那种视频在线观看| 国产一级毛片七仙女欲春2| 亚洲综合色惰| 久久午夜福利片| 国产 一区 欧美 日韩| 小说图片视频综合网站| avwww免费| 中亚洲国语对白在线视频| 国产在视频线在精品| ponron亚洲| 久久久久久久午夜电影| 久9热在线精品视频| 精品人妻一区二区三区麻豆 | 亚洲人成网站在线播放欧美日韩| 久久99热这里只有精品18| 一区二区三区免费毛片| 又黄又爽又刺激的免费视频.| 国产精品一区二区性色av| 国产一区二区亚洲精品在线观看| 51国产日韩欧美| 亚洲成av人片在线播放无| 村上凉子中文字幕在线| 男女之事视频高清在线观看| 欧美日韩亚洲国产一区二区在线观看| av福利片在线观看| 精品乱码久久久久久99久播| 简卡轻食公司| 成人永久免费在线观看视频| 精品久久久久久久久亚洲 | 午夜精品久久久久久毛片777| 久久久久久久久中文| 自拍偷自拍亚洲精品老妇| 一个人观看的视频www高清免费观看| 精品午夜福利视频在线观看一区| 亚洲va日本ⅴa欧美va伊人久久| 赤兔流量卡办理| 欧美一区二区精品小视频在线| 一个人免费在线观看的高清视频| 特大巨黑吊av在线直播| 国产精品国产高清国产av| 色综合亚洲欧美另类图片| 成人午夜高清在线视频| 国产视频内射| 午夜视频国产福利| 自拍偷自拍亚洲精品老妇| 69av精品久久久久久| 最近最新中文字幕大全电影3| 色哟哟哟哟哟哟| 色综合站精品国产| 成人鲁丝片一二三区免费| 成人国产综合亚洲| 我的女老师完整版在线观看| 最近在线观看免费完整版| 国产亚洲精品av在线| 国产精品一区二区性色av| 丰满的人妻完整版| 成人无遮挡网站| 亚洲一区高清亚洲精品| 国内精品一区二区在线观看| 日本精品一区二区三区蜜桃| 青草久久国产| 国产亚洲精品av在线| 午夜福利欧美成人| 国产成人影院久久av| 久久国产精品人妻蜜桃| 动漫黄色视频在线观看| 久久久久久国产a免费观看| 欧美黄色淫秽网站| 少妇高潮的动态图| 日韩欧美三级三区| 日韩亚洲欧美综合| 国产精品人妻久久久久久| 国产一区二区在线观看日韩| 在线观看午夜福利视频| 免费人成在线观看视频色| 日本与韩国留学比较| 久久午夜福利片| 窝窝影院91人妻| 亚洲精品日韩av片在线观看| 99riav亚洲国产免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产高清激情床上av| 亚洲欧美日韩高清专用| 国产伦一二天堂av在线观看| 一进一出好大好爽视频| 成年免费大片在线观看| 一级黄色大片毛片| 国产成年人精品一区二区| 桃色一区二区三区在线观看| 757午夜福利合集在线观看| 久久久久久久久大av| 欧美色欧美亚洲另类二区| 国产爱豆传媒在线观看| 久99久视频精品免费| 亚洲电影在线观看av| 一边摸一边抽搐一进一小说| 久久香蕉精品热| 久久亚洲真实| av天堂中文字幕网| 一边摸一边抽搐一进一小说| 中文亚洲av片在线观看爽| 久久亚洲真实| 久久国产精品影院| 久久久久九九精品影院| 亚洲第一区二区三区不卡| www.www免费av| 免费无遮挡裸体视频| а√天堂www在线а√下载| 在线观看舔阴道视频| a级毛片a级免费在线| av专区在线播放| 哪里可以看免费的av片| 看片在线看免费视频| 一二三四社区在线视频社区8| 中文在线观看免费www的网站| 精品一区二区免费观看| 中文在线观看免费www的网站| 国产精品久久久久久精品电影| 校园春色视频在线观看| 最新中文字幕久久久久| 直男gayav资源| 国产伦精品一区二区三区视频9| 亚洲男人的天堂狠狠| 国产不卡一卡二| 欧美成人免费av一区二区三区| 无遮挡黄片免费观看| 久久久国产成人精品二区| 最好的美女福利视频网| 首页视频小说图片口味搜索| eeuss影院久久| 老司机福利观看| 成人三级黄色视频| 精品日产1卡2卡| 搡老岳熟女国产| 极品教师在线视频| 草草在线视频免费看| 日本三级黄在线观看| 久久精品国产亚洲av香蕉五月| 久久精品国产亚洲av天美| 九九久久精品国产亚洲av麻豆| 蜜桃久久精品国产亚洲av| 国产主播在线观看一区二区| 级片在线观看| 一二三四社区在线视频社区8| 亚洲三级黄色毛片| 久久久久久大精品| 亚洲,欧美,日韩| 最近中文字幕高清免费大全6 | 熟女电影av网| 亚洲最大成人手机在线| 国产爱豆传媒在线观看| 老司机午夜十八禁免费视频| 搡老岳熟女国产| 成人亚洲精品av一区二区| 少妇丰满av| 亚洲人成网站在线播| 最新中文字幕久久久久| 国产 一区 欧美 日韩| 亚洲精品色激情综合| 怎么达到女性高潮| 综合色av麻豆| 日韩国内少妇激情av| 久久久久久久精品吃奶| 在线观看免费视频日本深夜| 精品欧美国产一区二区三| 不卡一级毛片| 在线免费观看不下载黄p国产 | 日本免费a在线| 少妇熟女aⅴ在线视频| 亚洲av熟女| 看十八女毛片水多多多| 久久性视频一级片| 91字幕亚洲| 非洲黑人性xxxx精品又粗又长| 自拍偷自拍亚洲精品老妇| 国产单亲对白刺激| 黄色女人牲交| 日韩大尺度精品在线看网址| 一个人免费在线观看的高清视频| 免费看光身美女| 男女视频在线观看网站免费| 午夜两性在线视频| 日韩有码中文字幕| 精品人妻1区二区| 国产视频内射| 国产高清激情床上av| 国产大屁股一区二区在线视频| 亚洲18禁久久av| 可以在线观看毛片的网站| 亚洲av电影在线进入| 搞女人的毛片| 久久天躁狠狠躁夜夜2o2o| 中文亚洲av片在线观看爽| 麻豆国产av国片精品| 日韩欧美国产一区二区入口| 精品日产1卡2卡| 国产欧美日韩一区二区三| 男女那种视频在线观看| 久久久久久大精品| 日韩精品中文字幕看吧| 男插女下体视频免费在线播放| 国产高潮美女av| 免费在线观看亚洲国产| 最新中文字幕久久久久| 久久久精品大字幕| 亚洲最大成人av| 欧美精品啪啪一区二区三区| 成人无遮挡网站| 悠悠久久av| 在线看三级毛片| 欧美bdsm另类| 成人美女网站在线观看视频| or卡值多少钱| 国产白丝娇喘喷水9色精品| av在线观看视频网站免费| 性色avwww在线观看| 不卡一级毛片| 99热这里只有是精品在线观看 | 天堂网av新在线| 国产亚洲欧美在线一区二区| 日本熟妇午夜| 一区二区三区免费毛片| ponron亚洲| 性欧美人与动物交配| 国内揄拍国产精品人妻在线| 人妻丰满熟妇av一区二区三区| 噜噜噜噜噜久久久久久91| 国内精品一区二区在线观看| 毛片女人毛片| 最好的美女福利视频网| 在线播放国产精品三级| 女生性感内裤真人,穿戴方法视频| 中文字幕高清在线视频| 国产精品日韩av在线免费观看| 国内精品久久久久精免费| 亚洲av.av天堂| 99久国产av精品| 日韩精品青青久久久久久| 一本一本综合久久| 欧美日韩综合久久久久久 | 黄片小视频在线播放| 最近最新中文字幕大全电影3| 国产蜜桃级精品一区二区三区| 人妻久久中文字幕网| 少妇人妻一区二区三区视频| 国内精品久久久久久久电影| 嫩草影视91久久| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼| 最近在线观看免费完整版| 色哟哟哟哟哟哟| 欧美成狂野欧美在线观看| 日韩中文字幕欧美一区二区| 最后的刺客免费高清国语| 亚洲精品一卡2卡三卡4卡5卡| 日韩av在线大香蕉| 婷婷亚洲欧美| 欧美黑人巨大hd| 色哟哟哟哟哟哟| 日本黄色片子视频| 国产91精品成人一区二区三区| 波多野结衣巨乳人妻| 欧美zozozo另类| 永久网站在线| 国产成人欧美在线观看| 国产私拍福利视频在线观看| 亚洲成av人片在线播放无| 天天躁日日操中文字幕| 国产乱人视频| 亚洲精品乱码久久久v下载方式| АⅤ资源中文在线天堂| 亚洲欧美日韩东京热| 给我免费播放毛片高清在线观看| 午夜免费成人在线视频| 级片在线观看| 亚洲性夜色夜夜综合| 在线a可以看的网站| 无人区码免费观看不卡| 草草在线视频免费看| av专区在线播放| 欧美激情在线99| 人妻夜夜爽99麻豆av| 久久草成人影院| 亚洲18禁久久av| 搞女人的毛片| 真人一进一出gif抽搐免费| 亚洲国产精品sss在线观看| 国产色爽女视频免费观看| 听说在线观看完整版免费高清| 一级黄色大片毛片| 免费在线观看成人毛片| 久久久成人免费电影| 亚洲中文日韩欧美视频| 一区二区三区免费毛片| 亚洲精品在线观看二区| 国产av一区在线观看免费| a级毛片a级免费在线| 变态另类成人亚洲欧美熟女| 免费在线观看亚洲国产| 宅男免费午夜| 一本综合久久免费| 青草久久国产| 美女免费视频网站| 免费人成在线观看视频色| 狠狠狠狠99中文字幕| 我的老师免费观看完整版| 婷婷丁香在线五月| 99在线人妻在线中文字幕| 成人国产综合亚洲| 亚洲人与动物交配视频| 色综合亚洲欧美另类图片| 91在线观看av| 禁无遮挡网站| av国产免费在线观看| 亚洲三级黄色毛片| 免费人成视频x8x8入口观看| 又粗又爽又猛毛片免费看| 久久伊人香网站| 亚洲av电影在线进入| 亚洲人成网站在线播放欧美日韩| 最新在线观看一区二区三区| 亚洲精品色激情综合| 日本熟妇午夜| 亚洲成人久久爱视频| 国产免费av片在线观看野外av| 99久久99久久久精品蜜桃| 俺也久久电影网| 少妇人妻一区二区三区视频| 99热精品在线国产| 欧美成人a在线观看| 99国产极品粉嫩在线观看| 九色国产91popny在线| 精品免费久久久久久久清纯| 12—13女人毛片做爰片一| 欧美高清性xxxxhd video| 国产一区二区激情短视频| 欧美激情久久久久久爽电影| 狂野欧美白嫩少妇大欣赏| 香蕉av资源在线| 国产精品99久久久久久久久| 真人一进一出gif抽搐免费| 久久精品久久久久久噜噜老黄 | 国产真实伦视频高清在线观看 | 国产三级中文精品| 日韩欧美国产在线观看| 国产黄色小视频在线观看| 日韩欧美 国产精品| 国产精品永久免费网站| 在线观看美女被高潮喷水网站 | 日韩高清综合在线| 午夜影院日韩av| 精品国内亚洲2022精品成人| 中亚洲国语对白在线视频| 亚洲美女黄片视频| 成人一区二区视频在线观看| 欧美成人一区二区免费高清观看| 亚洲成人精品中文字幕电影| 免费观看精品视频网站| 色精品久久人妻99蜜桃| 亚洲成人精品中文字幕电影| 窝窝影院91人妻| 18禁裸乳无遮挡免费网站照片| 丁香欧美五月| 小蜜桃在线观看免费完整版高清| 男人的好看免费观看在线视频| 18禁黄网站禁片免费观看直播| 中文字幕高清在线视频| 免费观看的影片在线观看| 男人舔女人下体高潮全视频| 9191精品国产免费久久| 真人一进一出gif抽搐免费| 国产伦精品一区二区三区视频9| 国产伦一二天堂av在线观看| 精品熟女少妇八av免费久了| 久久精品国产清高在天天线| 人妻制服诱惑在线中文字幕| 国产麻豆成人av免费视频| 亚州av有码| 婷婷精品国产亚洲av| 久久精品国产99精品国产亚洲性色|