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

    Relationship on the Head Deviation Rate and Laying Length of Porous Tubes

    2019-11-05 00:59:42YamengWangandHuanfangLiu

    Yameng Wang and Huanfang Liu

    (College of Water Conservancy and Architecture Engineering, Shihezi University, Shihezi 832000, Xinjiang, China)

    Abstract: To improve the utilization efficiency of water, the laying length of porous tubes was analyzed. Considering long porous tubes, the irrigation uniformity of porous tubes was studied by analyzing the deviation rate of the working head of porous tubes and through an impact analysis of the deviation ratio of flow rate based on laying length, slope, pressure head, and perforation spacing. Results showed that the deviation ratio of the flow rate could be reduced by decreasing the length of tubes and increasing the slope. The deviation ratio of flow rate was less influenced by pressure head and perforation spacing. A best-laid slope resulted in a longer laying length of the porous tubes, and they were not constant.

    Keywords: flow deviation rate; porous tubes; deviation rate of working head; laying length; head loss

    1 Introduction

    Porous tubes have been employed in many areas, such as petrochemicals, agriculture irrigation, water supply, and drainage[1-2]. The uniformity of fluid distribution is an important indicator of engineering quality[3-4], which is used to determine the pressure head and the length in many irrigation systems, and the design requirements of uniformity should be achieved to ensure quality and efficiency. In actual engineering systems, uniformity can be affected by many factors, including working pressure changes of emitter, manufacturing deviation of emitter, congestion, water temperature changes, and micro topography changes[5-10]. To date, some results for calculating flow uniformity have been achieved by domestic scholars. Examples for the method of hydraulics design, as studied by Kang[11], is in accordance with both average flow rate and emitter uniformity[11-12]. The methods for calculating comprehensive flow and ultimate flow deviation rates were proposed based on the analysis of the existing rates of hydraulic deviation and flow deviation[3,6-9]. Scaloppi and Allen[13]applied a differential approach to multiple outlet pipes with constant and continuously variable outflows and simulated the pressure distributions along uniform sprinkle systems, trickle irrigation laterals, manifolds, and gated pipes, considering the effect of ground slope and velocity head on the pipeline hydraulics. Hathoot et al.[14-15]investigated the problem of a lateral pipe with equally spaced emitters and a uniform slope, and estimated the head loss between emitters using the Darcy-Weisbach formula with variation in the Reynolds number, different zones on the Moody diagram, and a friction coefficient formula corresponding to each zone. Many results were achieved by predecessors scholars, but few have studied the interaction relation between flow uniformity and laying length, pressure head, laying slope, as well as perforation spacing[16-21], which has effects on practical engineering applications.

    Another problem with the momentum method is that it ignores the influence of the slope on the flow in the process of establishing a variable mass momentum equation, which is only regarded as a horizontal slope. According to mechanics, the momentum method only considers the influence of hydrodynamic pressure, friction, and momentum changes on the flow, but neglects the effect of gravity. For most perforated pipes used in micro-irrigation, slope is not horizontal and gravity might have a significant influence on the flow. Therefore, gravity should be introduced into the momentum equation of the variable mass to investigate the influence of the slope.

    Flow uniformity in porous tubes is significantly affected by the hole system. Gravitational pressure and forced pressure have different effects on the porous tubes distribution systems. Being the primary problem of flow uniformity, free porous flow relies on gravity and can be affected by many factors, so it is necessary to be further studied. This paper analyzes the flow uniformity in porous tubes, which is affected by laying length, pressure head, laying slope, and perforation spacing. The results of this study are helpful to improve the flow uniformity of porous tubes, reduce investment, and guide the construction and operation of porous tubes. The objectives are to develop a momentum equation of variable mass flow using the momentum conservation method under the condition of a certain slope, as well as solve the flow deviation rate and the deviation rate of working head.

    2 Expression of the Flow Deviation Rate

    While designing water distribution projects, engineers only consider the influence of flow uniformity affected by the manufacturing deviation of the emitter and the hydraulic friction loss. Of the two factors, manufacturing deviation is determined by the production technology level, and if it is high, flow uniformity is low, despite of hydraulic friction loss. Although manufacturers have improved the level of technology and produced emitters with smaller manufacturing deviations, the flow uniformity of the system remains poor without a reasonable hydraulic design. In the present study, manufacturing deviation was ignored and the factors of hydraulic design were considered to analyze data. The flow uniformity can be expressed in many ways, and the following computational formula[22]was adopted in this study:

    1)The flow deviation rate:

    (1)

    whereqvrepresents the flow deviation rate of emitter;qmax(L/h) represents the maximum flow of emitter;qmin(L/h) represents the minimum flow of emitter; andqa(L/h) represents the flow of design.

    When the flow index is 0.5, the flow formula of the pipe hole of porous tubes[8]is

    (2)

    whereqi(L/h) represents the flow of perforation;h(m) represents the working pressure head of perforation; andkprepresents the coefficient.

    If the flow of pipe hole is regarded as the design flow, the formula can be expressed as

    (3)

    Combining with the formulas in Eqs.(1)-(3), the following formula is obtained:

    (4)

    wherehmax(m) represents the maximum of working head, andhmin(m) represents the minimum of working head.

    2)Deviation rate of the working head.

    In irrigation systems, the pressure head varies according to the position of the emitter, thus causing a change in flow, which leads to flow deviation. The difference between the working pressure and the design of the pressure heads is called head deviation. The difference between the maximum and the minimum deviations in the entire system to the design pressure head is the deviation rate of the working head. Therefore, it can be written as

    (5)

    whereha(m) represents the design of working head;Δhmax(m) represents the maximum of head deviation; andΔhmin(m) represents the minimum of head deviation.

    In light of previous research results, the formula of the relationship between the deviation of working head and the flow deviation rate is obtained by using the binomial expansion method[23]as

    qv=khv

    (6)

    wherekrepresents the state index of emitter.

    3 Analysis of the Deviation Rate in the Working Head

    Since mass decreases along the flow in the porous tube flow or the embranchment flow, it is characterized as variable mass flow. Pressure distribution is affected by momentum exchange, head loss, and slope. The velocity of the flow decreases with changes in mass because of the holes in the sidewall, so the momentum changes along the flow and a certain amount of kinetic energy is converted into the pressure head, which causes the pressure to increase along the flow. Friction produces pressure losses that lead to the decrease of the pressure, and at the same time, the flow might result in rough waves at the outlet of the hole, which can increase energy losses. The head change induced by the slope of the pipe has a significant effect on the pressure distribution. Therefore, it is necessary to simultaneously take these influence factors into consideration if the pressure change is maintained in a certain range.

    Fig.1 shows the schematic of the pipeline with multiple outlets. It can be observed that the perforated tube is closed at the end of one side and the cross-sectional area is the same along the flow, which has one row of radial holes. In practical engineering, factors such as laying length, laying slope, pressure head, and perforation spacing simultaneously exert effects, which influence the deviation rate of the working head of the system concurrently. The average pressure head should be considered as the design pressure head during the design process to make the average pressure consistent with the design pressure.

    Fig.1 Schematic of pipeline with multiple holes

    Then, Eq.(5) can be rewritten as

    (7)

    The average, maximum, and minimum pressure heads with different laying lengths, pressure heads, and slopes can be computed by the calculation formula for the flow, as well as processing the head loss and the pressure head along the pipe[24-25]. In practical engineering, long porous tubes are more commonly used than other porous tubes. Therefore, the following formula investigates the influence of the deviation rate of the working head as a factor that affects flow uniformity in long porous tubes. According to the general principle, the conditions of 2 000

    The different formula for resistance coefficient along the way was adopted in the range of different Reynolds numbers. Different results of the calculating head loss were analyzed, and the formula for calculating the head loss along the pipe is as follows[24]:

    (8)

    whereχrepresents correction coefficient;frepresents friction coefficient;D(m) represents the diameter of pipe;Q0(m3/s) represents the flow of porous tubes inlet; andLrepresents the length of pipe. This formula is closer to the true head loss in the physical model than those in previous studies.

    3.1 Analysis of the Deviation Rate of the Working Head on a Flat Slope

    In the case of a flat slope, the pressure curve was in decline, the pressure head of the first hole was the maximum, and that of the last hole was the minimum. Therefore,

    hmax=h1,hmin=hz=h1-Hf

    whereh1represents the pressure head of the first hole, andhzrepresents the pressure head of the last hole.

    Combined with Eq.(7), then

    (9)

    It can be inferred from Eq.(9) that the deviation rate of the working head as well as the difference between the first pressure head and the frictional head loss were inversely related, provided that there was a flat slope for the long porous tubes without a velocity head. The greater the increase in the head loss was, the greater the deviation in the increase in the working head was. The deviation rate of the working head decreased with an increase in the first pressure head, which showed that the flow uniformity was high when the deviation rate of the working head was small. Thus, the first pressure head could be increased to improve the flow uniformity and reduce the deviation rate of the working head.

    The frictional head loss along the pipe was related to the total flow and the laying length of the porous tubes. If the length of the pipe was longer, the head loss was greater, the deviation rate of working head was bigger, and the flow uniformity was smaller. The head loss is controlled by the first pressure head and the total flow of the pipe, so a steady flow occurred when they were balanced. When the first pressure head increased, the deviation rate of the working head also decreased. Moreover, the total flow increased with the increase in the first pressure head, while it caused a deviation in the rate of increase in the working head. Therefore, the deviation rate of the working head decreased with an increase in the first pressure head, though the influence was not obvious, which was also observed for flow uniformity.

    3.2 Analysis of the Deviation Rate of the Working Head on an Adverse Slope

    In adverse slope (I<0) conditions, the pressure curve declined. The pressure head of the first hole was the maximum, and that of the last hole was the minimum. Combining the maximum and the minimum pressure heads, as well as the average pressure head and Eq. (7), we can obtain

    (10)

    whereLrepresents the length of pipe, andIrepresents the slope of pipe.

    Compared with Eqs.(9) and (10), the numerator obtained in Eq.(10) was greater than that in Eq.(9), and the denominator in Eq.(10) was smaller than that in Eq.(9). Therefore, in the same conditions, the deviation rate of the working head of the adverse slope was greater than that of the flat slope, because under the condition of an adverse slope, the first pressure head not only compensated for the head loss caused by the adverse slope but also overcame the frictional head loss. From Eq.(10), it can also be observed that the influence was similar to the laying length and the first pressure head for the deviation rate of the working head of the adverse slope. Thus, it can improve flow uniformity and decrease the deviation rate of the working head by reducing the laying length or increasing the first pressure head.

    3.3 Analysis of the Deviation Rate of the Working Head on a Positive Slope

    In positive slope (I>0) conditions, the head change along the porous tubes was complicated. The maximum pressure head appeared on the first hole section or the end section, and the minimum pressure head was random. Therefore, the analysis of the deviation rate of the working head should depend on the location of the minimum pressure head. When the minimum pressure head is in the middle of the pipe, the formula for the minimum pressure head is

    (11)

    hz=h0+IL-Hf

    (12)

    (13)

    The first condition is that the pressure head of the end section equals to the pressure head of the first section, and the maximum pressure head appears in the first section.IL=Hfandhmax=h1. According to Eq. (7), the deviation rate of the working head can be calculated as

    (14)

    A comparison of Eqs.(14) and (9) showed that the structures of the two formulas were similar. The deviation rate of the working head decreased as the first pressure head increased. The frictional head loss was larger, which could lead to an increase in the deviation rate of the working head. Under the same conditions, the numerator in Eq.(14) was smaller than that in Eq.(9), and the denominator in Eq.(14) was greater than that in Eq.(9). The result illustrated that the deviation rate of the working head of the positive slope was smaller than that of the flat slope in the same conditions. Thus, the flow uniformity of the flat slope was higher than that of the flat slope in the same condition. Besides, it could improve flow uniformity and decrease the deviation rate of the working head by reducing the laying length or increasing the first pressure head.

    The second condition is that the pressure head of the end section is smaller than that of the first section, and the maximum pressure head appears in the first section.hmax=h1and 0

    (15)

    In this formula, the numerator in Eq.(15) was smaller than that in Eq.(9), and the denominator in Eq.(14) was larger than that in Eq.(9). Thus, the deviation rate of the working head of the positive slope was smaller than that of the flat slope in the same conditions. The flow uniformity of the flat slope was bigger than that of the flat slope in the same condition. This is because the head loss was compensated by the head increment on the positive slope. Hence, improved flow uniformity could be improved by increasing the slope appropriately. In construction design, the slope of the terrain can be adopted to improve the flow uniformity.

    The third condition is that the pressure head of the end section is bigger than that of the first section, and the maximum pressure head occurs in the end section. 0

    (16)

    For Eq.(16), under the same conditions, the slope increased, and both the denominator and the numerator increased. However, the numerator increased faster than the denominator, indicating that the deviation rate of the working head of the positive slope was bigger than that of the flat slope in this condition. Thus, increasing only the slope would decrease the flow uniformity.

    According to the above analysis, whenI>0, the deviation rate of the working head could be decreased by reducing the laying length and increasing the first pressure head, which could also improve the flow uniformity. When 0

    The fourth condition is that the pressure head on the end section is the maximum, and the first section is the minimum.

    IL/Hf>2.75,hmin=h1,hmax=h1+IL-Hf

    According to Eq.(7), the deviation rate of the working head is as follows:

    (17)

    For Eq.(17), the increase of the first pressure head could reduce the deviation rate of the working head. Moreover, if (IL-Hf) is increased, the deviation rate of the working head is increased. Thus, the flow uniformity can be improved by increasing the first pressure head or reducing (IL-Hf).

    4 Analysis of the Laying Length

    The laying length of the porous tubes was dependent on the flow uniformity. Thus, the utilization efficiency of water was the best when the flow uniformity was higher and the flow deviation rate was lower. However, it would increase the investment and the operating cost. Therefore, flow deviation rate should be determined on the basis of the requirement with a suitable laying length. Currently, the flow deviation rate of less than 20% of the system is considered good. Although the chemical industry has a higher demand for flow deviation rates, most rates fall below 10%. Thus, the total flow computational formula, the head loss computational formula, and the pressure head computational formula can be utilized in different conditions, such as different slopes and pressure heads, to calculate the average pressure head, the maximum pressure head, and the minimum pressure head. Hence, Eq.(4) was used to calculate the flow deviation rate in this study.

    4.1 Effect of the Slope on Laying Length

    When the pressure head of inletH=1 m, the perforation spacingS=0.2 m, the diameter of pipeD=40 mm, and the diameter of perforationd=2 mm. The porous tubes are made of steel. The calculated results are listed in Table 1.

    It can be seen from Fig.2 that the flow deviation rates and the laying lengths of the porous pipe were different when the slope was fixed. With the decrease of the flow deviation rate, the laying length of the porous tube was gradually reduced. This is because the longer the laying length is, the greaterHfis, and the greater the flow deviation rate is, the smaller the outflow uniformity is. Under the premise of keeping the flow deviation rate unchanged, with the laying slope changed from the inverse slope (I<0) to the positive slope (0≤I≤0.1), the laying length first increased and then decreased. Table 1 shows different flow deviation rates for the best-laid slope in different conditions. For example, whenqv=15%, the longest laying length wasI=0.003, and the pipe could be laid on this slope with the length of 38.2 m. With the pressure head being constant, an appropriate increase in the slope could increase the laying length. The reason is that the head loss is less than the slope head increment caused by the slope. However, when the slope increased to a certain extent and continued to increase asI>0.003, the laying length would become shorter, because the friction head loss was less than the enhanced head due to the slope. It indicated that there is an optimal slope that can ensure the pipe has the longest length. From the table, it can be inferred that when the slope was constant, the laying length was longer, the flow deviation rate was bigger, and the flow uniformity was lower. The cause is that the increase of the length of the laying pipe accelerates the head loss, leading to the porous tube hole flow unstable and the difference between the maximum flow and the minimum flow greater, thus resulting in larger flow deviation rate. In summary, the maximum laying length of the porous pipe can be obtained at the slope of 0.002-0.003.

    Table 1 Laying length L(m) under H=1 m, S=0.2 m, D=40 mm (d=2 mm)

    Fig.2 Scatter plot of Table 1

    4.2 Effect of the Slope on Pressure Head

    To explore the relationship between the pressure head and the laying length, the laying length in a different slope condition was discussed when the perforation spacingS=0.2 m, the diameter of pipeD=40 mm, the diameter of perforationd=2 mm, and the flow deviation rateqv=15%, as shown in Table 2.

    It can be seen from Table 2 that when the slope is on the inverse slope (I<0) and remains unchanged, the laying length of the porous tube increased with the increase of pressure head. When the slope 00.004 and is unchanged, the laying length increased first and then decreased with the increase of pressure head. Under the premise ofqv=15%, where the best laying slope was different under different pressure heads, indicating that the optimum laying length was different. It can be seen from Table 2 that the laying length increased with the increase in the pressure head when the pressure headH=0.4 m and the slopeI≤0.003. IfI≥0.004, for a constant slope, increasing the pressure head properly could increase the laying length under the premise of ensuring the flow deviation rate. However, if the pressure head is too large, the laying length should be reduced to ensure the flow deviation rate.

    4.3 Effect of Slope on Perforation Spacing

    Table 3 shows the laying length of porous tubes with different perforation spacing and slopes, when the diameter of pipeD=40 mm, the diameter of perforationd=2 mm, and the flow deviation rateqv=15%.

    Table 2 Laying length L(m) under S=0.2 m, qv=15%, D=40 mm (d=2 mm)

    Table 3 Laying length L(m) under qv=15%, D=40 mm (d=2 mm)

    It can be seen in Table 3 that different perforation spacing corresponded to an optimal slope. In this slope, the porous tubes were longer than the lengths of other slopes. The optimal laying slope corresponding to different perforation spacing was varied. For example, whenS=0.2 m, the optimal slopeI=0.003, and the laying length was 38.2 m. WhenS=0.3 m, the optimal slopeI=0.007, and the laying length was 53.1 m. When the slope was not changed, the laying length increased with the increase in perforation spacing. This occurred because with a larger perforation spacing, the total flow of the porous tubes was smaller, which could reduce the head loss.

    In summary, whenI>0 and the pressure head is not changed, there is always a best slope for laying the longest porous tubes. In this case, the head loss is exactly equal to the incremental head obtained by the slope, and the maximum length can thus be achieved. Notably, influencing factors of flow uniformity are laying length, pressure head, pipe diameter, and perforation spacing. As these factors change, the total flow, first pressure head, and total head loss along the pipe change. Therefore, the best laying length is variable. The above test was only performed in a certain pipe diameter and a certain perforation diameter to get the maximum laying length and the best slope, while the maximum laying length would be different with different pipe diameters or perforation spacing. Therefore, the best slope can only be obtained in a positive slope, because only in the positive slope will an incremental head be produced to compensate for the head loss, and the slope can reach the best only when the two are balanced. In the flat slope and the adverse slope, the weight of the component was inversely related to the flow, thus making the pressure distribution gap larger and the flow uniformity smaller, which however, did not appear to be the best slope. The above analysis was consistent with the law of the deviation rate of the working head analysis. Considering the problem of ground flatness in practical engineering, the head loss would increase slightly, and the optimum laying length would be reduced accordingly.

    5 Conclusions

    For flat slope, the longer the length was, the larger the deviation rate of the working head was. The flow uniformity decreased withthe increase in laying length. The change in the first pressure head had only a slight effect on the deviation rate of the working head, and the increase in pressure head was not significant for increasing the flow uniformity. Under the same condition, for the adverse slope, the deviation rate of the working head was larger than that of the flat slope of the pipe. Thus, the flow uniformity was less than that of the flat slope. The effects of the laying length and the pressure head on the deviation rate of the working head of the adverse slope were similar to those of the flat slope porous tubes.

    The positive slope based on the ratio of IL andHfwas used to analyze the effect of the slope, the pressure head, and other factors on the deviation rate of the working head and the flow uniformity of the impact. For maintaining a certain flow deviation rate, there was an optimal slope that could make the pipe obtain the longest laying length, which was greater than the laying length in the other slope. Moreover, this optimal slope must be in the positive slope. Increasing the pressure head appropriately would increase the length of the pipe. However, if the pressure head was too large to ensure the flow deviation rate, the length of the pipe would decrease. For the same slope, increasing the perforation spacing would increase the length of pipe.

    国产三级黄色录像| 在线视频色国产色| 亚洲电影在线观看av| 精品一区二区三区视频在线 | 成年免费大片在线观看| 日本三级黄在线观看| 国产欧美日韩精品一区二区| 国产精品乱码一区二三区的特点| 黄色日韩在线| 欧美三级亚洲精品| 日本免费一区二区三区高清不卡| 一a级毛片在线观看| 国产欧美日韩一区二区三| 亚洲精品国产精品久久久不卡| 久久天躁狠狠躁夜夜2o2o| 在线看三级毛片| www.熟女人妻精品国产| 在线播放无遮挡| 老熟妇乱子伦视频在线观看| 精品一区二区三区视频在线观看免费| 亚洲精品亚洲一区二区| 国产蜜桃级精品一区二区三区| 床上黄色一级片| 成人午夜高清在线视频| 狂野欧美白嫩少妇大欣赏| 一卡2卡三卡四卡精品乱码亚洲| 熟女少妇亚洲综合色aaa.| 国产成人欧美在线观看| 尤物成人国产欧美一区二区三区| 欧美大码av| 国产又黄又爽又无遮挡在线| svipshipincom国产片| 制服人妻中文乱码| 精品一区二区三区av网在线观看| 美女被艹到高潮喷水动态| 午夜日韩欧美国产| 亚洲国产精品合色在线| 国产精品久久久久久久电影 | 亚洲欧美一区二区三区黑人| 91久久精品国产一区二区成人 | 男女那种视频在线观看| a级一级毛片免费在线观看| 国产一区二区三区视频了| 精品乱码久久久久久99久播| 在线观看日韩欧美| 久久久久久久久中文| 亚洲专区中文字幕在线| 我要搜黄色片| 特大巨黑吊av在线直播| 国产免费男女视频| 日韩欧美在线乱码| 欧美成人a在线观看| 噜噜噜噜噜久久久久久91| 亚洲国产欧美人成| 香蕉丝袜av| bbb黄色大片| 怎么达到女性高潮| 国产熟女xx| 在线国产一区二区在线| 亚洲人成网站在线播放欧美日韩| 1000部很黄的大片| 亚洲片人在线观看| 叶爱在线成人免费视频播放| tocl精华| 欧美另类亚洲清纯唯美| 午夜视频国产福利| 国产一级毛片七仙女欲春2| 午夜福利欧美成人| 亚洲av成人av| 深爱激情五月婷婷| 最近最新免费中文字幕在线| 1024手机看黄色片| 久久久国产成人免费| 国产69精品久久久久777片| 国产精品永久免费网站| 亚洲在线自拍视频| 在线播放无遮挡| 在线观看66精品国产| 日韩精品青青久久久久久| 搡老岳熟女国产| 国产成年人精品一区二区| 露出奶头的视频| 色播亚洲综合网| 一级黄色大片毛片| 日韩欧美国产一区二区入口| av天堂在线播放| 搡女人真爽免费视频火全软件 | 亚洲第一电影网av| 日韩欧美 国产精品| 久久久久久久午夜电影| av片东京热男人的天堂| 又黄又粗又硬又大视频| 久久久久久久久大av| 亚洲av免费高清在线观看| 99精品在免费线老司机午夜| 国产野战对白在线观看| 一区二区三区免费毛片| 可以在线观看毛片的网站| 嫩草影视91久久| 叶爱在线成人免费视频播放| 国产高清videossex| 亚洲 欧美 日韩 在线 免费| 亚洲无线在线观看| 熟女少妇亚洲综合色aaa.| 在线a可以看的网站| 国产私拍福利视频在线观看| 一级黄片播放器| 好看av亚洲va欧美ⅴa在| 一级a爱片免费观看的视频| 夜夜爽天天搞| 国产综合懂色| 成人国产一区最新在线观看| 男女之事视频高清在线观看| 欧美中文综合在线视频| 日韩免费av在线播放| 欧美一区二区亚洲| 色综合欧美亚洲国产小说| 精品福利观看| 亚洲av美国av| 色视频www国产| 成人高潮视频无遮挡免费网站| 一区福利在线观看| 国产美女午夜福利| 国产精品99久久久久久久久| 黄片小视频在线播放| 一本一本综合久久| 两性午夜刺激爽爽歪歪视频在线观看| 岛国视频午夜一区免费看| 麻豆成人午夜福利视频| 岛国在线观看网站| 88av欧美| 国产成人av教育| 内射极品少妇av片p| 搡女人真爽免费视频火全软件 | 欧美日韩精品网址| 亚洲精品国产精品久久久不卡| 夜夜夜夜夜久久久久| 亚洲最大成人中文| 岛国视频午夜一区免费看| 热99re8久久精品国产| 天天一区二区日本电影三级| 老汉色av国产亚洲站长工具| 精品午夜福利视频在线观看一区| 黄色成人免费大全| 国产国拍精品亚洲av在线观看 | 悠悠久久av| 亚洲性夜色夜夜综合| 婷婷丁香在线五月| 一进一出抽搐gif免费好疼| 国产精品自产拍在线观看55亚洲| 亚洲久久久久久中文字幕| 精品国产亚洲在线| 一夜夜www| 亚洲第一电影网av| 国产三级在线视频| 亚洲av免费高清在线观看| 9191精品国产免费久久| 色综合亚洲欧美另类图片| tocl精华| 少妇熟女aⅴ在线视频| 黄色女人牲交| x7x7x7水蜜桃| 亚洲人成网站在线播放欧美日韩| 免费大片18禁| 久久久精品大字幕| 国产精品香港三级国产av潘金莲| 成人性生交大片免费视频hd| 在线十欧美十亚洲十日本专区| 黄片小视频在线播放| 中文字幕精品亚洲无线码一区| 91九色精品人成在线观看| 婷婷六月久久综合丁香| 无人区码免费观看不卡| 精品不卡国产一区二区三区| 国产三级黄色录像| 欧洲精品卡2卡3卡4卡5卡区| 欧美+亚洲+日韩+国产| 91在线精品国自产拍蜜月 | av在线天堂中文字幕| 蜜桃久久精品国产亚洲av| 欧美三级亚洲精品| 欧美zozozo另类| 精品一区二区三区av网在线观看| 午夜福利欧美成人| 亚洲人成伊人成综合网2020| 精品国产亚洲在线| 亚洲av一区综合| 制服人妻中文乱码| 亚洲成人中文字幕在线播放| 亚洲黑人精品在线| 成熟少妇高潮喷水视频| 亚洲中文日韩欧美视频| 亚洲中文字幕一区二区三区有码在线看| 又黄又粗又硬又大视频| netflix在线观看网站| 午夜两性在线视频| 少妇丰满av| 国产一区二区三区在线臀色熟女| 国产欧美日韩一区二区精品| 日韩欧美国产在线观看| 亚洲色图av天堂| 日本与韩国留学比较| 久久久国产精品麻豆| 99精品久久久久人妻精品| 国产高清三级在线| 三级男女做爰猛烈吃奶摸视频| 一a级毛片在线观看| 女人十人毛片免费观看3o分钟| 色av中文字幕| 不卡一级毛片| 91在线精品国自产拍蜜月 | 极品教师在线免费播放| 欧美国产日韩亚洲一区| 久久人人精品亚洲av| www.www免费av| 少妇人妻一区二区三区视频| 日韩亚洲欧美综合| 黄色日韩在线| 午夜亚洲福利在线播放| 一区二区三区国产精品乱码| 88av欧美| 99久久成人亚洲精品观看| 欧美+亚洲+日韩+国产| 中国美女看黄片| 欧美日本视频| 亚洲内射少妇av| 亚洲av电影不卡..在线观看| 精品人妻1区二区| 午夜老司机福利剧场| 亚洲午夜理论影院| 岛国视频午夜一区免费看| 免费看日本二区| 欧美日韩黄片免| 岛国在线免费视频观看| 久久精品国产综合久久久| 宅男免费午夜| 国产精品久久电影中文字幕| 午夜免费成人在线视频| 久99久视频精品免费| 亚洲av五月六月丁香网| 俄罗斯特黄特色一大片| 看片在线看免费视频| 婷婷六月久久综合丁香| 欧美一级毛片孕妇| 最好的美女福利视频网| bbb黄色大片| 欧美日韩综合久久久久久 | 国产午夜精品久久久久久一区二区三区 | 国产一区二区在线av高清观看| 法律面前人人平等表现在哪些方面| 波多野结衣高清无吗| 午夜老司机福利剧场| 久久精品国产亚洲av涩爱 | 在线观看免费午夜福利视频| 999久久久精品免费观看国产| 两个人视频免费观看高清| 国产欧美日韩精品一区二区| 国产高清视频在线观看网站| 午夜久久久久精精品| 欧美av亚洲av综合av国产av| 97超级碰碰碰精品色视频在线观看| 国产极品精品免费视频能看的| 狠狠狠狠99中文字幕| 99热精品在线国产| 国产精品免费一区二区三区在线| 特大巨黑吊av在线直播| 在线视频色国产色| 国产亚洲欧美98| 两个人的视频大全免费| а√天堂www在线а√下载| 天堂av国产一区二区熟女人妻| 国产精品一及| 精品久久久久久久久久久久久| 最好的美女福利视频网| 国产中年淑女户外野战色| 国产极品精品免费视频能看的| 99久久无色码亚洲精品果冻| 91九色精品人成在线观看| 国产av不卡久久| av黄色大香蕉| 久久这里只有精品中国| 中出人妻视频一区二区| 村上凉子中文字幕在线| 老司机午夜十八禁免费视频| ponron亚洲| 十八禁人妻一区二区| 日韩精品中文字幕看吧| 最近最新中文字幕大全免费视频| 一本综合久久免费| 亚洲狠狠婷婷综合久久图片| 小蜜桃在线观看免费完整版高清| 变态另类丝袜制服| 免费人成在线观看视频色| 18禁国产床啪视频网站| 99久久成人亚洲精品观看| 丁香欧美五月| 99久久综合精品五月天人人| 91av网一区二区| 美女 人体艺术 gogo| av在线蜜桃| www.www免费av| 国产精品一区二区三区四区免费观看 | 最新美女视频免费是黄的| 色噜噜av男人的天堂激情| 欧美色欧美亚洲另类二区| 一进一出抽搐gif免费好疼| 高清毛片免费观看视频网站| 久久国产精品影院| 色精品久久人妻99蜜桃| 国产黄片美女视频| 51午夜福利影视在线观看| 看黄色毛片网站| 国产av不卡久久| 中文在线观看免费www的网站| 欧美bdsm另类| 精品久久久久久成人av| 精品国产三级普通话版| 日韩欧美在线二视频| 亚洲精品日韩av片在线观看 | 国产免费男女视频| 18禁国产床啪视频网站| 欧美色视频一区免费| 国产亚洲精品综合一区在线观看| 亚洲 国产 在线| 嫩草影院入口| av专区在线播放| 在线视频色国产色| 91久久精品国产一区二区成人 | 岛国在线观看网站| 色噜噜av男人的天堂激情| 日韩高清综合在线| 免费在线观看日本一区| 精品久久久久久,| 欧美日韩福利视频一区二区| 国产精品乱码一区二三区的特点| 美女黄网站色视频| 久久精品国产99精品国产亚洲性色| 国产一区二区三区视频了| 免费人成视频x8x8入口观看| 国产精品野战在线观看| 男人舔奶头视频| 久久久久九九精品影院| 国产探花极品一区二区| 久久久久久久午夜电影| 国产色婷婷99| 国产精品香港三级国产av潘金莲| 欧美日韩福利视频一区二区| 国产毛片a区久久久久| 成年女人永久免费观看视频| 欧美另类亚洲清纯唯美| 91字幕亚洲| 真人一进一出gif抽搐免费| 香蕉久久夜色| 深夜精品福利| 18禁黄网站禁片午夜丰满| 又黄又粗又硬又大视频| 91麻豆精品激情在线观看国产| 两人在一起打扑克的视频| 小说图片视频综合网站| 狠狠狠狠99中文字幕| 日日摸夜夜添夜夜添小说| 悠悠久久av| 一个人免费在线观看电影| 国产高清videossex| 757午夜福利合集在线观看| 99精品欧美一区二区三区四区| 国产麻豆成人av免费视频| 在线观看免费视频日本深夜| 国内毛片毛片毛片毛片毛片| 免费人成在线观看视频色| 国产高清激情床上av| 国产日本99.免费观看| 欧美日韩精品网址| 88av欧美| 51国产日韩欧美| 一级黄色大片毛片| 久久久久久久午夜电影| 亚洲天堂国产精品一区在线| 一进一出抽搐gif免费好疼| 天美传媒精品一区二区| 亚洲av中文字字幕乱码综合| 亚洲国产精品sss在线观看| 久久九九热精品免费| 激情在线观看视频在线高清| 88av欧美| 国产成人欧美在线观看| av在线蜜桃| 久久久久久久精品吃奶| 国产色爽女视频免费观看| 人人妻,人人澡人人爽秒播| 日日干狠狠操夜夜爽| 在线观看舔阴道视频| 免费av毛片视频| 97超级碰碰碰精品色视频在线观看| 精品不卡国产一区二区三区| 色综合站精品国产| 色av中文字幕| 亚洲专区中文字幕在线| 国产精品 国内视频| 免费一级毛片在线播放高清视频| 一级a爱片免费观看的视频| 久久久久国内视频| 真人做人爱边吃奶动态| 久久久成人免费电影| 首页视频小说图片口味搜索| 两性午夜刺激爽爽歪歪视频在线观看| 黄片小视频在线播放| 99久久精品国产亚洲精品| 免费无遮挡裸体视频| 成人特级av手机在线观看| 亚洲天堂国产精品一区在线| 亚洲美女黄片视频| 麻豆成人午夜福利视频| 亚洲最大成人中文| 亚洲av一区综合| 成人精品一区二区免费| 国产91精品成人一区二区三区| 又爽又黄无遮挡网站| 国产极品精品免费视频能看的| 免费人成视频x8x8入口观看| 又黄又粗又硬又大视频| 欧美一区二区亚洲| 一区二区三区国产精品乱码| 五月玫瑰六月丁香| 中出人妻视频一区二区| 亚洲乱码一区二区免费版| 非洲黑人性xxxx精品又粗又长| 亚洲av二区三区四区| 乱人视频在线观看| 丁香六月欧美| 久久久久久久久久黄片| 成熟少妇高潮喷水视频| 国产男靠女视频免费网站| 精品无人区乱码1区二区| 人妻夜夜爽99麻豆av| 欧美成人性av电影在线观看| 久久人妻av系列| 91字幕亚洲| 免费看十八禁软件| 成人国产综合亚洲| 制服丝袜大香蕉在线| 中文资源天堂在线| or卡值多少钱| 午夜视频国产福利| 久久精品影院6| 久久久久精品国产欧美久久久| 久久亚洲精品不卡| 色视频www国产| 欧美三级亚洲精品| 老鸭窝网址在线观看| 国产熟女xx| 欧美黑人欧美精品刺激| 在线看三级毛片| 精品不卡国产一区二区三区| 搡老妇女老女人老熟妇| 中文字幕久久专区| 18禁黄网站禁片免费观看直播| 久久午夜亚洲精品久久| 在线观看免费视频日本深夜| 日韩欧美精品免费久久 | 中文字幕久久专区| 色综合站精品国产| 午夜日韩欧美国产| av在线天堂中文字幕| 国产亚洲精品av在线| 超碰av人人做人人爽久久 | 欧美乱妇无乱码| 欧美极品一区二区三区四区| 麻豆成人av在线观看| 男女之事视频高清在线观看| 最新美女视频免费是黄的| 欧美日本亚洲视频在线播放| 国产亚洲av嫩草精品影院| 欧美精品啪啪一区二区三区| 国产精品精品国产色婷婷| 夜夜躁狠狠躁天天躁| 国产主播在线观看一区二区| 成年女人毛片免费观看观看9| 亚洲成人精品中文字幕电影| 99久久综合精品五月天人人| 国产一区二区在线观看日韩 | 欧美日韩瑟瑟在线播放| а√天堂www在线а√下载| 91九色精品人成在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲国产日韩欧美精品在线观看 | 国产aⅴ精品一区二区三区波| 18禁国产床啪视频网站| 国产精品一区二区免费欧美| 热99在线观看视频| 久久久国产成人精品二区| 国产三级在线视频| 成人鲁丝片一二三区免费| 少妇人妻一区二区三区视频| 成年女人毛片免费观看观看9| 美女免费视频网站| 狂野欧美激情性xxxx| 成人国产一区最新在线观看| 国产日本99.免费观看| 亚洲国产精品999在线| 久久精品国产99精品国产亚洲性色| 精品久久久久久,| 午夜福利成人在线免费观看| av天堂在线播放| 一级作爱视频免费观看| 中文字幕av成人在线电影| 国产中年淑女户外野战色| 亚洲欧美日韩无卡精品| 又爽又黄无遮挡网站| 亚洲精品日韩av片在线观看 | 国产一区二区三区视频了| 国产真实乱freesex| 嫩草影视91久久| 舔av片在线| 99久久99久久久精品蜜桃| 我要搜黄色片| 国产精品国产高清国产av| 色播亚洲综合网| 国产精品98久久久久久宅男小说| 中亚洲国语对白在线视频| 日韩中文字幕欧美一区二区| 禁无遮挡网站| av在线蜜桃| 熟女人妻精品中文字幕| 精品人妻1区二区| 99久国产av精品| 热99在线观看视频| 天堂影院成人在线观看| 国产欧美日韩精品亚洲av| 免费一级毛片在线播放高清视频| 最近视频中文字幕2019在线8| 国产亚洲精品av在线| 69人妻影院| 精品福利观看| 美女cb高潮喷水在线观看| 日韩欧美 国产精品| xxxwww97欧美| 日韩 欧美 亚洲 中文字幕| 亚洲天堂国产精品一区在线| 99精品欧美一区二区三区四区| 国产激情欧美一区二区| 久久久久久大精品| 9191精品国产免费久久| 国产黄色小视频在线观看| www日本在线高清视频| 国产精品久久久久久人妻精品电影| 婷婷精品国产亚洲av| 一区二区三区免费毛片| 日本黄色视频三级网站网址| 在线观看日韩欧美| 免费看a级黄色片| 校园春色视频在线观看| 高清日韩中文字幕在线| 99热这里只有是精品50| 丰满乱子伦码专区| 国内精品久久久久精免费| 19禁男女啪啪无遮挡网站| 午夜日韩欧美国产| 舔av片在线| 欧美+日韩+精品| 制服丝袜大香蕉在线| 欧美不卡视频在线免费观看| 久久午夜福利片| 国产一级毛片七仙女欲春2| 亚洲成人久久爱视频| 男女边摸边吃奶| 亚洲一区高清亚洲精品| 日日啪夜夜撸| 国产视频内射| 自拍偷自拍亚洲精品老妇| 亚洲久久久久久中文字幕| 精品午夜福利在线看| 亚洲精品一二三| 最近中文字幕高清免费大全6| 内射极品少妇av片p| 91久久精品国产一区二区成人| 99视频精品全部免费 在线| 99久国产av精品国产电影| 亚洲精品亚洲一区二区| 99久国产av精品| 国产精品一及| 成人一区二区视频在线观看| 国产精品日韩av在线免费观看| 欧美高清成人免费视频www| 日日摸夜夜添夜夜爱| 免费观看av网站的网址| 亚洲人与动物交配视频| 精华霜和精华液先用哪个| 中文字幕免费在线视频6| 哪个播放器可以免费观看大片| 中文字幕人妻熟人妻熟丝袜美| 床上黄色一级片| 国产成人一区二区在线| 卡戴珊不雅视频在线播放| 亚洲国产高清在线一区二区三| 好男人在线观看高清免费视频| 亚洲激情五月婷婷啪啪| 特大巨黑吊av在线直播| 亚洲av男天堂| 成人一区二区视频在线观看| 18禁在线播放成人免费| 国产白丝娇喘喷水9色精品| 国内精品宾馆在线| 日日啪夜夜撸| 99热这里只有精品一区| 国产男人的电影天堂91| 看黄色毛片网站| 久久精品久久精品一区二区三区| 美女脱内裤让男人舔精品视频| 精品人妻熟女av久视频| 亚洲精品久久久久久婷婷小说| 久久久久久久大尺度免费视频| 日日撸夜夜添| 波多野结衣巨乳人妻| 国精品久久久久久国模美|