• <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.

    久久亚洲精品不卡| 久99久视频精品免费| 中出人妻视频一区二区| 日日夜夜操网爽| 在线观看免费视频日本深夜| 18禁黄网站禁片免费观看直播| 女人高潮潮喷娇喘18禁视频| 久久中文看片网| 国产黄色小视频在线观看| 制服诱惑二区| 日本黄色视频三级网站网址| 国产精品香港三级国产av潘金莲| 国产精品九九99| 亚洲成av人片免费观看| 成熟少妇高潮喷水视频| 熟妇人妻久久中文字幕3abv| 国产av在哪里看| 女性生殖器流出的白浆| 久久久国产成人免费| 国产在线精品亚洲第一网站| 中文亚洲av片在线观看爽| 巨乳人妻的诱惑在线观看| 亚洲熟妇熟女久久| 精品欧美国产一区二区三| 欧美大码av| 亚洲精品一卡2卡三卡4卡5卡| 国产av在哪里看| 淫秽高清视频在线观看| 人人妻人人澡人人看| 亚洲在线自拍视频| 欧美三级亚洲精品| 国产精品免费一区二区三区在线| 亚洲av熟女| 国产单亲对白刺激| 免费在线观看视频国产中文字幕亚洲| 色哟哟哟哟哟哟| 精品乱码久久久久久99久播| 丝袜人妻中文字幕| www.熟女人妻精品国产| av免费在线观看网站| 男人舔女人下体高潮全视频| 日韩av在线大香蕉| 99热这里只有精品一区 | 国产成人影院久久av| 校园春色视频在线观看| 激情在线观看视频在线高清| 中文字幕人成人乱码亚洲影| 日日干狠狠操夜夜爽| 国产成人欧美在线观看| 欧美乱色亚洲激情| 村上凉子中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 美女 人体艺术 gogo| 国产精品野战在线观看| 久久久久久人人人人人| 亚洲精品色激情综合| 十八禁人妻一区二区| 高清在线国产一区| 麻豆成人av在线观看| 国内揄拍国产精品人妻在线 | 国产成人欧美在线观看| 非洲黑人性xxxx精品又粗又长| 黄色丝袜av网址大全| 亚洲欧美精品综合一区二区三区| 美女免费视频网站| 搡老妇女老女人老熟妇| 久久久久国产一级毛片高清牌| 人人妻人人看人人澡| 欧美又色又爽又黄视频| 嫩草影院精品99| 午夜老司机福利片| 精品卡一卡二卡四卡免费| 亚洲av成人不卡在线观看播放网| 国产一区二区在线av高清观看| 给我免费播放毛片高清在线观看| 国产又黄又爽又无遮挡在线| 麻豆久久精品国产亚洲av| 亚洲熟妇中文字幕五十中出| 可以在线观看的亚洲视频| 麻豆一二三区av精品| 亚洲专区中文字幕在线| videosex国产| aaaaa片日本免费| 一本综合久久免费| 国产又爽黄色视频| 成人午夜高清在线视频 | 国产一卡二卡三卡精品| 成年人黄色毛片网站| 精品第一国产精品| 亚洲精品国产区一区二| 琪琪午夜伦伦电影理论片6080| 99在线视频只有这里精品首页| 国产单亲对白刺激| 亚洲国产精品合色在线| 久久狼人影院| 一级片免费观看大全| 成人18禁在线播放| 国产99白浆流出| 亚洲av成人av| 中文字幕精品亚洲无线码一区 | 国产精品久久久av美女十八| av视频在线观看入口| 亚洲男人的天堂狠狠| 一级a爱视频在线免费观看| 国产精品久久久人人做人人爽| 69av精品久久久久久| 亚洲国产精品久久男人天堂| 黄网站色视频无遮挡免费观看| 欧美在线黄色| 人成视频在线观看免费观看| 中文在线观看免费www的网站 | 亚洲成人国产一区在线观看| 国产精华一区二区三区| 成人欧美大片| 亚洲真实伦在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美色欧美亚洲另类二区| 超碰成人久久| 国产免费av片在线观看野外av| 亚洲黑人精品在线| 黑人欧美特级aaaaaa片| 18禁观看日本| 国产黄色小视频在线观看| 日韩欧美在线二视频| 欧美成人午夜精品| 国产精品日韩av在线免费观看| 亚洲成人久久爱视频| 91在线观看av| 国产又色又爽无遮挡免费看| 欧美黑人巨大hd| 观看免费一级毛片| 久久久久久人人人人人| 国产精品av久久久久免费| 人妻久久中文字幕网| 一二三四在线观看免费中文在| 国产野战对白在线观看| 亚洲国产精品成人综合色| 中文在线观看免费www的网站 | 国产成人精品久久二区二区91| 精品午夜福利视频在线观看一区| 美女国产高潮福利片在线看| 精品乱码久久久久久99久播| 亚洲人成伊人成综合网2020| 老司机在亚洲福利影院| 男女床上黄色一级片免费看| 一进一出抽搐动态| 女性被躁到高潮视频| 精品国产美女av久久久久小说| 欧美黄色片欧美黄色片| 免费在线观看成人毛片| 老熟妇乱子伦视频在线观看| 欧美午夜高清在线| 别揉我奶头~嗯~啊~动态视频| 国产97色在线日韩免费| 一级毛片高清免费大全| 一区二区三区激情视频| 欧美黄色片欧美黄色片| 亚洲免费av在线视频| 成年版毛片免费区| www.999成人在线观看| 国内精品久久久久久久电影| 老熟妇乱子伦视频在线观看| 在线天堂中文资源库| 十分钟在线观看高清视频www| 香蕉久久夜色| 一进一出好大好爽视频| 久久久精品国产亚洲av高清涩受| 在线观看舔阴道视频| 日韩精品免费视频一区二区三区| 男人的好看免费观看在线视频 | 亚洲 欧美 日韩 在线 免费| 视频区欧美日本亚洲| 亚洲第一欧美日韩一区二区三区| 欧美精品啪啪一区二区三区| 丰满的人妻完整版| 女警被强在线播放| 成年版毛片免费区| 女生性感内裤真人,穿戴方法视频| 99精品久久久久人妻精品| 国产免费男女视频| 亚洲第一青青草原| 日韩一卡2卡3卡4卡2021年| 日本熟妇午夜| 欧美黄色片欧美黄色片| aaaaa片日本免费| av免费在线观看网站| 久久香蕉国产精品| 亚洲一区二区三区色噜噜| 叶爱在线成人免费视频播放| 精品第一国产精品| 亚洲国产毛片av蜜桃av| 后天国语完整版免费观看| 国产伦在线观看视频一区| 国产欧美日韩一区二区三| 99久久国产精品久久久| 少妇 在线观看| 视频区欧美日本亚洲| 国产激情久久老熟女| 少妇粗大呻吟视频| 免费人成视频x8x8入口观看| 中文字幕av电影在线播放| 久久狼人影院| bbb黄色大片| 成人国产综合亚洲| 欧美日韩中文字幕国产精品一区二区三区| 亚洲人成伊人成综合网2020| 亚洲人成网站在线播放欧美日韩| 成人免费观看视频高清| 日韩欧美在线二视频| 男人舔女人的私密视频| 1024视频免费在线观看| netflix在线观看网站| 久热这里只有精品99| 女警被强在线播放| 国产亚洲av高清不卡| 亚洲av电影在线进入| svipshipincom国产片| 日韩免费av在线播放| 国产成人精品久久二区二区91| 最新美女视频免费是黄的| 亚洲专区中文字幕在线| 亚洲精品久久国产高清桃花| 中文字幕高清在线视频| 国产精品,欧美在线| 日韩一卡2卡3卡4卡2021年| 黑人巨大精品欧美一区二区mp4| 99在线人妻在线中文字幕| 精品日产1卡2卡| 一区二区三区高清视频在线| 中国美女看黄片| www.熟女人妻精品国产| 一级黄色大片毛片| 丝袜美腿诱惑在线| 亚洲 欧美一区二区三区| 99re在线观看精品视频| bbb黄色大片| 精品少妇一区二区三区视频日本电影| 成人免费观看视频高清| 久久热在线av| 亚洲第一av免费看| 国产麻豆成人av免费视频| 婷婷精品国产亚洲av| 亚洲av电影在线进入| 九色国产91popny在线| 亚洲自拍偷在线| 麻豆国产av国片精品| 又黄又爽又免费观看的视频| 欧美性猛交黑人性爽| 久久午夜综合久久蜜桃| 国产精品久久久av美女十八| 无人区码免费观看不卡| 夜夜看夜夜爽夜夜摸| 男女之事视频高清在线观看| 久久国产精品人妻蜜桃| 一进一出抽搐gif免费好疼| 久久国产亚洲av麻豆专区| 在线观看66精品国产| 午夜a级毛片| 亚洲五月婷婷丁香| 成人一区二区视频在线观看| 神马国产精品三级电影在线观看 | 亚洲精品一区av在线观看| 国产一区二区激情短视频| 国内精品久久久久久久电影| 在线十欧美十亚洲十日本专区| 日韩欧美一区视频在线观看| 一级片免费观看大全| 看免费av毛片| 神马国产精品三级电影在线观看 | 制服丝袜大香蕉在线| www.熟女人妻精品国产| 一进一出抽搐gif免费好疼| 日韩视频一区二区在线观看| 女人被狂操c到高潮| 美女国产高潮福利片在线看| 99热只有精品国产| 亚洲欧美激情综合另类| 好看av亚洲va欧美ⅴa在| 人人妻人人澡欧美一区二区| 久久中文字幕人妻熟女| 三级毛片av免费| 精品国产国语对白av| 久久国产精品男人的天堂亚洲| 亚洲九九香蕉| 哪里可以看免费的av片| 美女 人体艺术 gogo| 两性午夜刺激爽爽歪歪视频在线观看 | 99久久综合精品五月天人人| 欧美不卡视频在线免费观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 老熟妇乱子伦视频在线观看| 日本a在线网址| 波多野结衣巨乳人妻| 欧美激情 高清一区二区三区| 国产成人欧美在线观看| 国产精品一区二区精品视频观看| 亚洲片人在线观看| 黄色片一级片一级黄色片| 中文亚洲av片在线观看爽| 日韩av在线大香蕉| 在线十欧美十亚洲十日本专区| 欧美激情极品国产一区二区三区| 国产精品永久免费网站| 国产主播在线观看一区二区| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品在线观看二区| 精品熟女少妇八av免费久了| 国产99白浆流出| 1024视频免费在线观看| 老司机在亚洲福利影院| 欧美在线一区亚洲| 亚洲va日本ⅴa欧美va伊人久久| 19禁男女啪啪无遮挡网站| xxx96com| 亚洲,欧美精品.| 校园春色视频在线观看| 一级a爱片免费观看的视频| 国产日本99.免费观看| 一区二区三区高清视频在线| 国产熟女午夜一区二区三区| 亚洲精品中文字幕在线视频| 午夜精品在线福利| 夜夜看夜夜爽夜夜摸| 久久精品91无色码中文字幕| 国产黄片美女视频| 老鸭窝网址在线观看| www日本在线高清视频| 香蕉丝袜av| 国产成人精品久久二区二区免费| 天堂√8在线中文| 一本精品99久久精品77| 日韩欧美一区二区三区在线观看| 人妻久久中文字幕网| 国产在线观看jvid| 国产日本99.免费观看| 成人国产一区最新在线观看| 热99re8久久精品国产| 久久婷婷成人综合色麻豆| 国产亚洲av嫩草精品影院| 久热爱精品视频在线9| 亚洲av熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 一个人免费在线观看的高清视频| 99精品在免费线老司机午夜| 哪里可以看免费的av片| 国产99久久九九免费精品| 男女床上黄色一级片免费看| 人妻丰满熟妇av一区二区三区| 国产精品久久久av美女十八| 狂野欧美激情性xxxx| 欧美日本亚洲视频在线播放| 欧美丝袜亚洲另类 | 亚洲一区中文字幕在线| 69av精品久久久久久| 欧美日韩瑟瑟在线播放| 欧美在线一区亚洲| 婷婷精品国产亚洲av| 99久久无色码亚洲精品果冻| 白带黄色成豆腐渣| 成人午夜高清在线视频 | 日韩高清综合在线| 国产精品亚洲一级av第二区| 中国美女看黄片| 在线观看66精品国产| 午夜福利在线观看吧| 黄色成人免费大全| 桃红色精品国产亚洲av| 亚洲成人免费电影在线观看| 成年人黄色毛片网站| 国产一级毛片七仙女欲春2 | 国产男靠女视频免费网站| 三级毛片av免费| 韩国精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 一进一出抽搐gif免费好疼| 午夜精品久久久久久毛片777| 国产亚洲av高清不卡| 又大又爽又粗| 欧美日韩黄片免| 亚洲成av人片免费观看| 国产高清视频在线播放一区| x7x7x7水蜜桃| 波多野结衣av一区二区av| 在线av久久热| 久久国产精品男人的天堂亚洲| 欧美黄色淫秽网站| 一a级毛片在线观看| 老司机福利观看| av欧美777| 欧美另类亚洲清纯唯美| 香蕉av资源在线| 最新美女视频免费是黄的| 在线免费观看的www视频| 最近最新中文字幕大全电影3 | 中出人妻视频一区二区| 欧美黑人精品巨大| 日韩欧美国产在线观看| 99精品欧美一区二区三区四区| 国产片内射在线| 日韩欧美一区二区三区在线观看| 亚洲天堂国产精品一区在线| 日本精品一区二区三区蜜桃| 麻豆一二三区av精品| 亚洲av电影在线进入| 国产熟女xx| 搡老妇女老女人老熟妇| 亚洲午夜精品一区,二区,三区| 久久中文字幕人妻熟女| 成人国产综合亚洲| a在线观看视频网站| 动漫黄色视频在线观看| 99精品久久久久人妻精品| 自线自在国产av| 精品一区二区三区四区五区乱码| 18禁观看日本| 国产麻豆成人av免费视频| 久久久久亚洲av毛片大全| 日本撒尿小便嘘嘘汇集6| 在线视频色国产色| √禁漫天堂资源中文www| 国产91精品成人一区二区三区| 日本在线视频免费播放| 国产精品一区二区三区四区久久 | 男女下面进入的视频免费午夜 | 国产亚洲精品第一综合不卡| 午夜福利成人在线免费观看| 精品国产美女av久久久久小说| 国产片内射在线| 亚洲av美国av| 一级毛片高清免费大全| 亚洲av第一区精品v没综合| 亚洲一区中文字幕在线| 亚洲五月婷婷丁香| 国产高清视频在线播放一区| 久久亚洲精品不卡| 亚洲成av片中文字幕在线观看| 精品欧美一区二区三区在线| 久热这里只有精品99| 精品第一国产精品| 亚洲精品粉嫩美女一区| 一本综合久久免费| 最近最新免费中文字幕在线| 亚洲精品国产区一区二| 日本撒尿小便嘘嘘汇集6| 美女午夜性视频免费| 日本三级黄在线观看| 国产精品久久久久久精品电影 | 久久精品aⅴ一区二区三区四区| √禁漫天堂资源中文www| 国产精品免费一区二区三区在线| 午夜免费鲁丝| 亚洲九九香蕉| 少妇 在线观看| 性欧美人与动物交配| 丁香欧美五月| 天堂影院成人在线观看| 在线天堂中文资源库| 黑人操中国人逼视频| 日本a在线网址| 精品少妇一区二区三区视频日本电影| 久久午夜亚洲精品久久| 国产精品久久久久久亚洲av鲁大| 桃色一区二区三区在线观看| 久久亚洲精品不卡| 50天的宝宝边吃奶边哭怎么回事| 色综合站精品国产| 不卡av一区二区三区| 母亲3免费完整高清在线观看| 国产精品1区2区在线观看.| 男人的好看免费观看在线视频 | 国产成年人精品一区二区| 中文字幕av电影在线播放| 久久人妻福利社区极品人妻图片| 成人国产综合亚洲| 国产av又大| 午夜免费成人在线视频| 久久伊人香网站| 精品乱码久久久久久99久播| 精品国产国语对白av| 久久久久久久精品吃奶| 精品一区二区三区视频在线观看免费| 日韩有码中文字幕| 美女午夜性视频免费| 精品第一国产精品| 亚洲性夜色夜夜综合| 老熟妇乱子伦视频在线观看| 国产熟女午夜一区二区三区| 国产在线观看jvid| 久久久久久久久免费视频了| 欧美性猛交╳xxx乱大交人| 欧美激情久久久久久爽电影| 亚洲色图 男人天堂 中文字幕| 好男人电影高清在线观看| 91麻豆av在线| 国产成年人精品一区二区| 精品国产国语对白av| 国产精品久久久久久人妻精品电影| 午夜激情av网站| 欧美色视频一区免费| 村上凉子中文字幕在线| 制服人妻中文乱码| 国产一区二区三区在线臀色熟女| 老司机深夜福利视频在线观看| 欧美黑人巨大hd| ponron亚洲| 精品日产1卡2卡| 中文字幕人妻丝袜一区二区| 久久久久免费精品人妻一区二区 | 亚洲av片天天在线观看| 啪啪无遮挡十八禁网站| 国产精品综合久久久久久久免费| 97超级碰碰碰精品色视频在线观看| 国产精品久久久av美女十八| 国产在线精品亚洲第一网站| 香蕉丝袜av| 黄网站色视频无遮挡免费观看| 99精品久久久久人妻精品| www日本黄色视频网| 欧美久久黑人一区二区| 免费观看人在逋| 亚洲中文字幕日韩| 成年免费大片在线观看| 日韩欧美三级三区| 欧美乱色亚洲激情| 国产色视频综合| 久久国产乱子伦精品免费另类| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区激情短视频| 亚洲午夜精品一区,二区,三区| 久久性视频一级片| 国产野战对白在线观看| 九色国产91popny在线| 俄罗斯特黄特色一大片| 国产精品亚洲一级av第二区| 亚洲成a人片在线一区二区| 1024香蕉在线观看| 亚洲国产日韩欧美精品在线观看 | 老熟妇仑乱视频hdxx| 亚洲av电影不卡..在线观看| 一二三四社区在线视频社区8| 中文资源天堂在线| 99riav亚洲国产免费| 2021天堂中文幕一二区在线观 | 一二三四社区在线视频社区8| 欧美日本视频| 91国产中文字幕| 国产亚洲精品第一综合不卡| 人人澡人人妻人| 亚洲欧美一区二区三区黑人| 精品国产超薄肉色丝袜足j| 狂野欧美激情性xxxx| 亚洲av成人av| 老熟妇乱子伦视频在线观看| 国产成人一区二区三区免费视频网站| 亚洲国产精品sss在线观看| 亚洲人成网站高清观看| 精品免费久久久久久久清纯| 999精品在线视频| 狂野欧美激情性xxxx| 亚洲精品美女久久av网站| 日韩精品中文字幕看吧| 国产麻豆成人av免费视频| 国产精品久久久久久精品电影 | cao死你这个sao货| 色哟哟哟哟哟哟| 校园春色视频在线观看| 久久久久久大精品| 熟女电影av网| www国产在线视频色| 在线天堂中文资源库| 国产乱人伦免费视频| 国内精品久久久久久久电影| 黄片小视频在线播放| 亚洲va日本ⅴa欧美va伊人久久| 午夜精品在线福利| 久久国产精品男人的天堂亚洲| 亚洲在线自拍视频| 欧美在线一区亚洲| 亚洲精品在线观看二区| 国产精品一区二区免费欧美| 大型av网站在线播放| 日韩欧美国产一区二区入口| 男女那种视频在线观看| 大型av网站在线播放| 久久欧美精品欧美久久欧美| 亚洲天堂国产精品一区在线| 男女之事视频高清在线观看| 日本免费一区二区三区高清不卡| 国产激情久久老熟女| 亚洲人成77777在线视频| 国产精品免费视频内射| 十分钟在线观看高清视频www| 欧美黄色片欧美黄色片| 又黄又爽又免费观看的视频| 非洲黑人性xxxx精品又粗又长| 色综合亚洲欧美另类图片| 一本大道久久a久久精品| 国产av又大| 日韩欧美一区视频在线观看| 韩国av一区二区三区四区| 好男人在线观看高清免费视频 | 啦啦啦观看免费观看视频高清| 亚洲av第一区精品v没综合| 12—13女人毛片做爰片一| 国产精品 欧美亚洲| 久久中文字幕一级| 99久久久亚洲精品蜜臀av| √禁漫天堂资源中文www| 免费看美女性在线毛片视频| 国产高清视频在线播放一区| 天堂影院成人在线观看| 国产高清有码在线观看视频 | 99re在线观看精品视频| 色播在线永久视频| 一本一本综合久久|