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

    Specific Enthalpy Based Heat Stress Index for Indoor Environments without Radiation Effect

    2022-05-09 07:00:28LIUShuxiao劉書曉LIUYunZOUYue

    LIU Shuxiao(劉書曉), LIU Yun(劉 赟), ZOU Yue(鄒 鉞)*

    1 College of Environmental Science and Engineering, Donghua University, Shanghai 201620, China2 Dianwei Science Co., Ltd., Shanghai 201620, China

    Abstract: Over 100 human thermal indices have been developed to predict the combined thermal impact on the body. In principle, these indices based on energy thermal budget equations should not only be the most complex but also be the most accurate. However, the simple indices based on algebraic or statistical models [e.g., the wet-bulb globe temperature (WBGT)] continue to be the most popular. A new heat stress index, the enthalpy dry-bulb temperature (EnD) for indoor environments is developed and validated in this study. The EnD index is unique in that it uses the air specific enthalpy, not the wet-bulb temperature, to measure the latent heat transfer from the skin to the surrounding environment. Theoretically, the EnD index can be treated as the equivalent temperature based on the convective heat transfer coefficient hc. Comparison is made between the EnD index and the widely used WBGT index based on the experimental data taken from three independent studies available in the scientific literature. The results show that the EnD index can reduce the overestimation of the dry-bulb air temperature and thus reduce heat stress in most cases, especially for hot and humid environments. It can be concluded that the EnD index has the potential to replace the WBGT index as the standard heat stress index in the future.

    Key words: heat stress; wet-bulb globe temperature(WBGT); specific enthalpy; indoor environment; maximum allowable duration of heat exposure

    Introduction

    According to the World Health Organization (WHO), occupational heat stress is becoming a significant problem in several regions, and global climate change will make the situation worse. In addition to the significant public health threat, excessive heat exposure affects economic productivity and reduces the hourly work capacity[1-2].

    The human heat strain response is the result of the interplay among many factors, which include the thermal state of the surrounding environment (e.g., air temperature, air velocity, humidity, and solar radiation), human physical and behavioral conditions (e.g., activity level, clothing, and posture), and individual characteristics. A large number of human thermal climate indices have been developed to predict the combined thermal impact on the body or the corresponding effects on physiological variables that differ according to the number of variables, sophistication level of the employed human body heat balance theory, and applications[3-5]. De Freitas and Grigorieva[4]sorted a total of 162 indices from the year of 1905 to the year of 2012 into eight primary classification classes: simulation devices for integrated measurement such as thermal manikins[6](group A: 11 indices); single-sensor indices such as the first thermal index (i.e., the wet-bulb temperature, which can be traced to Haldane[7]) (groups B: 6 indices); indices based on algebraic or statistical models such as the wet-bulb globe temperature (WBGT)[8](group C: 40 indices); two proxy thermal indices (strain and stress) such as the predicted four-hour sweat rate (P4SR)[9]and effective temperature (ET)[10](groups D and E: 28 indices); two energy thermal budget equation indices (strain and stress) such as the body-atmosphere energy exchange index (BIODEX)[11], and predicted heat strain (PHS)[12](groups F and G: 67 indices); special-purpose indices such as the draught risk index (DRI)[13](group H: 10 indices).

    Groups F and G are the largest groups of indices, and they are based on human body heat balance equations. Normally, these kinds of indices apply all climatic and personal condition parameters as the input data for simulating the heat exchange between the human body and the ambient environment, and output the detailed information on the corresponding physiological reaction parameters. In principle, these indices are not the most complex accurate[14]. Current computational power allows meteorological data for even complex models to be handled on a handheld device. However, the simple indices in group C based on algebraic or statistical models (e.g., WBGT) continue to be the most popular indices[5, 15].

    Most of the algebraic or statistical thermal climate indices in group C[4]include two or three of the following variables related to the different heat exchange processes between the human body and the ambient environment: dry-bulb temperaturetaand air velocity for convective heat transferVa; wet-bulb temperaturetwb, relative humidity (RH)RH, dewpoint temperaturetdew, and water vapor pressurePafor latent heat transfer; globe temperaturetgand solar radiationSfor radiation heat transfer; metabolic rateMand clothing insulationClofor describing the human body condition. Table 1 summarizes all parameters involved in 38 of 40 total indices; two Russian indices (i.e., the severity rating, and the thermal insulation characteristics of clothing) are not included because of the missing relevant literatures. In addition to above six fundamental parameters, the environmental air enthalpyHahas also been applied once as an index[16].

    Table 1 Occurrence frequency of variables in 38 statistical indices of group C

    Table 1 indicates that the different temperatures (ta,twb,tdew, andtg) are the most common variables applied in statistical indices. Consequently, many corresponding indices, such as WBGT[9], Oxford index (WD)[17], discomfort index (DI)[18], and wet-bulb dry temperature (WBDT)[19], are expressed in the format of a temperature-equivalent value or the ostensive temperature value. Most people can easily understand the thermal effect of any given temperature based on their experience. This utility may be a big reason that simple statistical indices are still the most popular. However, according to the classical heat and mass transfer theory[20], the enthalpy difference is the driving potential for the simultaneous transfer of heat and mass in an air-water/vapor mixture,i.e. sweat.

    The objective of this study is to develop and validate a new index in order to assess heat stress in hot indoor environments. The index is developed based on the heat balance equation between a nude individual and the environment in which the effects of radiation can be neglected. We will show in this paper that this new heat stress index is capable of predicting the maximum allowable duration of heat exposure and physiological variables under moderate and light workloads.

    1 Development of Heat Stress Index

    We assumes that the nude individual works in a hot indoor environment, and the individual has a total metabolic rate. The internal wall surface temperature is close to the skin surface temperature, and the radiative heat transfer between the individual and surroundings is negligible[21]. According to the ASHRAE handbook[20], the heat balance equation of the human body per skin surface area can be written as

    M-W=qsk+qres+S,

    (1)

    whereMis the total metabolic rate of metabolic heat production (W/m2),Wis the rate of mechanical work (W/m2),qskis the total rate of heat loss from the skin (W/m2), andqresis the total rate of heat loss through respiration (W/m2).

    For most human activities,Wis very small compared with the human metabolic rate. According to the ASHRAE handbook[20],Wcan be regarded as zero.qresis also quite limited in a hot indoor environment. Therefore,qskis the key parameter that indicates the extent to which the human body removes metabolic heat and keeps heat storage below hazardous values.

    For a nude individual with negligible radiative heat transfer, the heat lost from the skin dissipates to the immediate surroundings by the convection and evaporation of sweat. Hence,

    qsk=C+Esk=hc(tsk-ta)+whc(Psk, S-Pa) ,

    (2)

    whereCis the rate of convective heat transfer from the skin (W/m2),Eskis the rate of evaporative heat loss from the skin (W/m2),hcis the convective heat transfer coefficient [W/(m2·K)],tskis the skin surface temperature (K),wis the skin wettedness dimensionless,heis the evaporative heat transfer coefficient (analogous tohc) [W/(m2·kPa)],Psk, Sis the water vapor pressure at the skin surface which is assumed to be the saturated water vapor attsk(kPa).

    The Lewis ratioLcan be used to relate the surface convective and evaporative heat transfer coefficients at a low mass transfer rate. According to ASHRAE handbook[20], the Lewis ratio is about 16.5 K/kPa for typical indoor conditions.

    L=he/hc=16.5.

    (3)

    Substituting Eq. (3) into Eq. (2), and rearranging the parameters based on the skin surface and ambient air yield

    qsk=hc(tsk-ta)+16.5whc(Psk, S-Pa)=

    hc[(tsk+16.5wPsk, S)-(ta+16.5wPa)].

    (4)

    The water vapor pressurePaof the moist air can be expressed in terms of the total or barometric pressurePtand the humidity ratio of moist aird.

    (5)

    wheredrefers to the humidity ratio of moist air (kg/kg).

    Considering that the total pressurePtis about 101 kPa and the water vapor pressurePais always within a range of 1-5 kPa under normal indoor conditions, Eq. (5) can be simplified as

    (6)

    Substitute Eq. (6) into Eq. (4):

    qsk=hc[(tsk+26.53wdsk, S)-(ta+26.53wda)] ,

    (7)

    wheredsk, Sis the humidity ratio of air at the skin surface, which is assumed as the humidity ratio of saturated moist air attsk(kg/kg);dais the humidity ratio of ambient air (kg/kg).

    Enthalpy is one of the most important measures of energy in a thermodynamic system. For moist air, considering that the water vapor content is very small, the specific enthalpy of the moist air can be expressed as

    H=1.006t+2 500d,

    (8)

    whereHrefers to the specific enthalpy of the moist air (kJ/kg).

    Comparing Eq. (8) with Eq. (7) shows that the difference in the specific enthalpy plays an important role in dissipating heat from the skin surface to the surrounding air. The larger the skin wettedness is, the larger the specific enthalpy difference is. Skin wettedness is a good measure of heat stress and evaporative potential of the environment. The maximum theoretical skin moisture is 1.0, but most indoor environments need to consider the maximum safe working hours[20], and the dimensionless parameter is usually in the range of 0.5-0.8. Alber-Wallerstr?m and Holmér[22]also studied the efficiency of sweat evaporation of six men working on a bicycle ergometer at different workloads and humidity levels, and the average skin wettedness was about 0.7. If we takew= 0.66 and insert Eq. (8) into Eq. (7),qskcan be rewritten as

    qsk=hc[(tsk+1 751dsk, S)-(ta+1 751da)]≈
    hc[(0.3tsk+0.7(1.006tsk+2 500dsk, S))-
    (0.3ta+0.7(1.006tsk+2 500da))]=
    hc[(0.3tsk+0.7Hsk)-(0.3ta+0.7Ha)],

    (9)

    whereHskis the specific enthalpy of air at the skin surface, which is assumed to be the specific enthalpy of saturated moist air attsk(kJ/kg).

    If the individual is working in a hot indoor environment, the mean skin temperature is relatively stable, where the variation in skin temperature is always less than 2 ℃[24]. The maximum allowable duration of heat exposure is then determined based on the environmental parameters or combinations of the environmental parameters such as the WBGT index. According to Eq. (9), in an indoor environment where radiative heat transfer is negligible, the new heat stress index can be treated as the equivalent temperature based on the convective heat transfer coefficienthc. We name this new index enthalpy dry-bulb temperature (EnD) asEnD,because it is a combination of the specific enthalpy and dry-bulb temperature of the surrounding air.

    EnD=0.3ta+0.7Ha.

    (10)

    We have devised a new rational heat stress index based on the heat balance equation between a nude individual and the environment without radiation. The mathematical expression of the EnD index is quite similar to that of the WBGT index for indoor environments. However, the difference here is that the wet-bulb temperature is replaced with the specific enthalpy. Since the EnD index is developed based on the heat balance equation between a nude individual and a typical indoor environment, this condition can be difficult to fulfil in the most situations. However, we can validate the EnD index by carrying out thermal stress tests on individuals dressed in clothing with 0.6-1.1 clo, which are quite common for indoor environments such as deep mining or submarine rooms in developing countries.

    2 Validation of EnD Index

    The main objective of a heat stress index is to accurately predict the maximum allowable duration of heat exposure. Many standards have recommended that the individual should be removed from the heat exposure if he/she is unable to continue working. The individual should also be removed if his/her heart rate exceeds 180 beats/min for certain time or if his/her rectal temperature reaches 39.0 ℃[24]. For these reasons, the EnD index should correlate well with key physiological variables such as the rectal temperature and the heart rate.

    To validate the new EnD index, we compared the WBGT and EnD indices based on the experimental data taken from three independent studies in the scientific literature. All of these data are available publicly and used for different purposes. To eliminate variations (which were common in field experiments), all measurement data were collected in well-controlled climate chambers. It was important to note that all tests were carried out with the internal wall surface temperature close to the air temperature, not human skin temperature. However, radiative heat transfer coefficienthris nearly constant as 4.7 W/(m2·K) which is less thanhcfor the most indoor calculations[20]. According to Eq. (3), total error of ignoring radiation heat transfer can be estimated less than 6%.

    2.1 Experiments based on Lind and Hellon[17]

    The work of Lind and Hellon[17]was to investigate the tolerance time and physiological responses of low labor subjects in the high temperature and humidity environment. Ten healthy, and acclimatized young men were exposed to eight different indoor environments, where four were exposed to humid climates of which the dry-bulb temperature and the relative humidity were 37.0 ℃/85%, 38.0 ℃/85%, 39.0 ℃/85%, and 40.0 ℃/85%, and the remaining four were exposed to dry climates of which the dry-bulb temperature and the relative humidity were 44.7℃/45%, 46.1℃/45%, 48.9℃/45%, and 50.3℃/45%. All of the tests were carried out under conditions from which the effective temperature scale was originally developed. Therefore, the tests were carried out under the following conditions: sedentary activity (66 W/m2), light clothing (0.5-0.7 clo), and low wind velocity. The maximum duration of heat exposure was 180 min for all tests.

    Figure 1 shows the comparison between WBGT and EnD indices based on the tolerance time and two physiological variables (rectal temperature in 60 min and heartbeat in 60 min). All of the experimental data were reevaluated from the work of Lind and Hellon[17]. Since not all of the subjects were able to complete the maximum duration of exposure of 180 min, and the physiological variables of each subject might be different before the subjects entered the test chamber, only the variations of the physiological variables after the first 60 min of heat exposure were compared.

    Fig. 1 Comparison between WBGT and EnD indices(low labor): (a)WBGT index based on tolerance time; (b)EnD index based on tolerance time; (c)WBGT index based on rectal temperature in 60 min; (d)EnD index based on rectal temperature in 60 min; (e)WBGT index based on heartbeat in 60 min; (f)EnD index based on heartbeat in 60 min

    For the WBGT index, all of the test results fall into two distinct series based on the relative humidity (RH ≈ 85%, and RH ≈ 45%) at the start of the test. Based on the results, it is evident that the WBGT index is not an adequate measure of heat stress or physiological stress within the range of environmental conditions investigated in this study. Although there are deviations at higher EnD index values, the results obtained from the two series of climates (humid and dry climates) begin to merge into a single curve. The results indicate that it is possible to accurately assess the severity of heat stress based on the EnD index within the range of environmental conditions considered in this study. However, the deviations at higher EnD values imply that the new heat stress index still overestimates the stress imposed by the high dry-bulb temperature in severe heat. The excessive response of the WBGT index towards the dry-bulb temperature could result in significantly underestimation of heat stress in environments where the evaporation of sweat is restricted by high humidity levels or poor ventilation[29].

    2.2 Experiments based on Pimental and Avellini[23]

    One objective of the work of Pimental and Avellini[23]was to investigate the tolerance time and physiological responses of moderate labor subjects in an equivalent WBGT environment. Eight healthy, and well-acclimated young male subjects attempted to complete 4 h of heat exposure in five different environments. The dry-bulb temperature and the relative humidity of the five climates were 35.0 ℃/80%, 49.0 ℃/25%, 43.0 ℃/60%, 49.0 ℃/35%, and 49.0 ℃ /39%. All of the subjects wore the U.S. Navy Utility Uniform with a thermal insulation value of 1.1 clo and the wind velocity was 1 m/s. For the test, the subjects were required to walk on a level treadmill at a pace of 1.3 m/s for 25 min and then sat for 5 min every 30 min. The time-weighted metabolic rate was about 270 W(i.e., moderate labor).

    Figure 2 shows the comparison between WBGT and EnD indices based on the tolerance time, rectal temperature in 60 min, and heartbeat in 60 min. All of the experimental data were reevaluated from the work of Pimental and Avellin[23]. It can be seen that the WBGT index is not a suitable predictor of heat stress for the tolerance time and physiological variables. In contrast, the EnD index shows a stronger correlation with all of variables in this case, especially for the tolerance time. Interestingly, there is a striking similarity in the patterns for rectal temperature and heartbeat between Figs. 1 and 2. We perceive that perhaps only one physiological parameter needs to be measured in the future.

    Fig. 2 Comparison between WBGT and EnD indices (moderate labor): (a)WBGT index based on tolerance time; (b)EnD index based on tolerance time; (c)WBGT index based on rectal temperature in 60 min; (d)EnD index based on rectal temperature in 60 min; (e)WBGT index based on heartbeat in 60 min; (f)EnD index based on heartbeat in 60 min

    2.3 Experiments based on Lü[25]

    In order to develop a physiological index for heat tolerance limits in China, an artificial climatic chamber was built in Tianjin University to study the effect of heat exposure on the physiological variables in different working environments. A part of the research was well-documented in Ref.[25]. Even though the thesis was written in Chinese, the relevant data were used in the works of Refs. [26-28], which were all published in a renowned international journal.

    In Lü’s study, 148 healthy young male students were tested under 45 experimental combinations. Each student only took one test to ensure that none of the subjects was acclimatized to the hot environment. All of the subjects were dressed in the same type of uniform with a thermal insulation value of 0.6 clo and the wind velocity was about 0.2 m/s. Each subject was randomly assigned to perform three levels of task: light, moderate, and heavy tasks. In order to ensure that the number of subjects was adequate as the previous two studies for statistical analysis, we only reanalyzed the results obtained for more than five subjects in one test condition. All of the cases that fulfilled this criterion could be classified as moderate workloads, where the time-weight metabolic rate was about 300 W(i.e., medium high labor).

    Since the rectal temperature was not available, we only presented the comparison between WBGT and EnD indices for the tolerance time and the heartbeat in 60 min, as shown in Fig. 3. Based on the results, it is difficult to identify which index is more suitable to predict the heat stress. We perceive that the tolerance time versus the EnD index may result in better curve fitting compared with the tolerance time versus the WBGT index. However, there is a large scatter in the data points for the heartbeat for both heat stress indices. It is likely that all of the subjects are not acclimated to the heat during the tests, which leads to the large scatter of the data points.

    Fig. 3 Comparison between WBGT and EnD indices(medium high labor): (a)WBGT index based on tolerance time; (b)EnD index based on tolerance time; (c)WBGT index based on heartbeat in 60 min; (d)EnD index based on heartbeat in 60 min

    Epstein and Moran[3]established a compromise approach to describe the effects of environmental and behavioral factors on human thermal stress: a change of 1 clo was equivalent to a change in 10 ℃ inTawhile exercising, and a change of 0.1 m/s in air movement could be offset by 0.5 ℃ inTa. Although Lü[25], and Pimental and Avellini[23]did tests with different thermal insulation values and air speeds, the combination effect of these two factors on human stress was quite close. Considering that the results of these two studies were also obtained under a similar time-weight metabolic rate, we presented the comparison between WBGT and EnD indices for these results in Fig. 4. The most intriguing finding here is that most of the data points for tolerance time obtained from the two independent sources can be fitted into a single curve for the EnD index. This indicates the validity of the EnD index for similar work intensities.

    Fig. 4 Comparison between WBGT and EnD indices: (a)WBGT index based on tolerance time; (b)EnD index based on tolerance time; (c)WBGT index based on heartbeat in 60 min; (d)EnD index based on heartbeat in 60 min

    3 Conclusions

    In this study, the EnD index is proposed and validated against experimental data available in the scientific literature, which is derived from the heat balance equation between a nude individual and the environment in which the effects of radiation are negligible. The Lewis ratio is used to integrate evaporation of sweat with the heat lost from the skin surface to the surrounding air by convection. The EnD index can be treated as the equivalent temperature based on the convective heat transfer coefficienthc, and the EnD index is expressed as the combination of specific enthalpy and dry-bulb temperature of the surrounding air. The ratio of specific enthalpy to dry-bulb temperature is dependent on the skin wettedness.

    The present EnD index is quite similar to the WBGT index for indoor environments, and the primary difference is that the wet-bulb temperature is replaced with the specific enthalpy. With this approach, the effect of dry-bulb temperature is reduced dramatically because the specific enthalpy is significantly larger than the corresponding wet-bulb temperature. The higher the water vapor content is in the air, the lower the effect of dry-bulb temperature on the EnD index is. The excessive response towards the dry-bulb temperature is indeed one of the disadvantages of the WBGT index. Theoretically, the EnD index can overcome this limitation over a wider range.

    We also reevaluated the experimental data taken from three independent studies in the scientific literature in order to validate the EnD index. In all cases, it is found that the EnD index shows better correlation with the tolerance time, rectal temperature variation, and heartbeat variation after the first 60 min of the test compared with the WBGT index. In addition, the EnD index has a much wider range (50-120) compared with the WBGT index (25-40), indicating that the EnD index is more sensitive to slight changes in the environmental conditions. The EnD index is also easy to use because of its simplicity.

    国产成人免费无遮挡视频| 99九九在线精品视频| 亚洲熟女毛片儿| 国产亚洲欧美在线一区二区| 国产成人av激情在线播放| 国产xxxxx性猛交| 国产麻豆69| 亚洲国产看品久久| 无限看片的www在线观看| 三级毛片av免费| 国产欧美日韩精品亚洲av| 午夜福利视频精品| 亚洲免费av在线视频| 18禁裸乳无遮挡动漫免费视频| 亚洲七黄色美女视频| 久久毛片免费看一区二区三区| 曰老女人黄片| 人妻久久中文字幕网| 大香蕉久久网| 精品久久久精品久久久| 久久天躁狠狠躁夜夜2o2o| 妹子高潮喷水视频| 亚洲国产av新网站| a级片在线免费高清观看视频| 亚洲av成人一区二区三| 一区在线观看完整版| 午夜福利欧美成人| 一级黄色大片毛片| 免费久久久久久久精品成人欧美视频| 亚洲久久久国产精品| 热99国产精品久久久久久7| 国产成人av激情在线播放| 人人妻人人澡人人看| 乱人伦中国视频| 国产亚洲午夜精品一区二区久久| 精品国产乱码久久久久久男人| 日日夜夜操网爽| 国产成人精品无人区| 丝瓜视频免费看黄片| 国产人伦9x9x在线观看| 热99久久久久精品小说推荐| av国产精品久久久久影院| 精品久久久久久久毛片微露脸| 在线播放国产精品三级| 免费黄频网站在线观看国产| 亚洲成人手机| 熟女少妇亚洲综合色aaa.| a级毛片在线看网站| 久久狼人影院| 久久国产亚洲av麻豆专区| 亚洲av欧美aⅴ国产| 亚洲专区中文字幕在线| 国产在线精品亚洲第一网站| 国产一区二区三区在线臀色熟女 | 老鸭窝网址在线观看| 一区二区三区乱码不卡18| 久久中文字幕一级| 淫妇啪啪啪对白视频| 精品乱码久久久久久99久播| 热re99久久国产66热| 精品国产一区二区三区久久久樱花| 国产无遮挡羞羞视频在线观看| 可以免费在线观看a视频的电影网站| 精品国产一区二区三区久久久樱花| 国产精品麻豆人妻色哟哟久久| 成人永久免费在线观看视频 | 国产成人免费无遮挡视频| 欧美亚洲 丝袜 人妻 在线| 天天躁狠狠躁夜夜躁狠狠躁| 免费观看人在逋| av天堂在线播放| 国产欧美日韩一区二区精品| 午夜福利,免费看| 一级黄色大片毛片| 久久这里只有精品19| 黄色视频,在线免费观看| 人妻 亚洲 视频| 97人妻天天添夜夜摸| 丰满人妻熟妇乱又伦精品不卡| 亚洲五月婷婷丁香| 天堂中文最新版在线下载| 精品视频人人做人人爽| 亚洲欧美日韩另类电影网站| 午夜免费鲁丝| 精品熟女少妇八av免费久了| 性色av乱码一区二区三区2| 国产亚洲精品第一综合不卡| 老司机深夜福利视频在线观看| 桃红色精品国产亚洲av| 亚洲第一青青草原| 国产在线观看jvid| 99国产精品免费福利视频| 老司机影院毛片| 啦啦啦免费观看视频1| 久久久久久久久免费视频了| 国产欧美日韩精品亚洲av| 午夜91福利影院| 精品国产一区二区久久| 99国产精品免费福利视频| 国产精品久久久久久精品电影小说| 亚洲色图 男人天堂 中文字幕| 久久精品aⅴ一区二区三区四区| 亚洲国产欧美网| 在线av久久热| 午夜福利在线免费观看网站| 国产精品影院久久| 欧美性长视频在线观看| 9色porny在线观看| 久久久国产成人免费| 大陆偷拍与自拍| 久久精品aⅴ一区二区三区四区| 一区二区三区激情视频| 9191精品国产免费久久| 日韩一卡2卡3卡4卡2021年| 午夜福利在线免费观看网站| 一区二区日韩欧美中文字幕| 久久婷婷成人综合色麻豆| 丰满少妇做爰视频| 国产片内射在线| 国产精品 国内视频| 99国产精品一区二区蜜桃av | 亚洲色图综合在线观看| 久久ye,这里只有精品| 视频区图区小说| 男女免费视频国产| 日韩欧美三级三区| 乱人伦中国视频| 亚洲avbb在线观看| 国产欧美日韩综合在线一区二区| 亚洲av电影在线进入| 黑人巨大精品欧美一区二区蜜桃| 人人澡人人妻人| 51午夜福利影视在线观看| 新久久久久国产一级毛片| 午夜福利视频精品| 麻豆国产av国片精品| 国产精品 国内视频| 国产一区二区激情短视频| 99国产极品粉嫩在线观看| 夜夜骑夜夜射夜夜干| 最新在线观看一区二区三区| 91精品三级在线观看| 久久久久久久国产电影| 岛国在线观看网站| 欧美成人免费av一区二区三区 | 我的亚洲天堂| 人人妻人人澡人人看| 18禁美女被吸乳视频| 久久精品aⅴ一区二区三区四区| 9191精品国产免费久久| 久久精品国产亚洲av高清一级| 考比视频在线观看| 成年版毛片免费区| 中文字幕最新亚洲高清| 午夜精品久久久久久毛片777| 精品久久久久久久毛片微露脸| 精品少妇黑人巨大在线播放| 巨乳人妻的诱惑在线观看| 在线亚洲精品国产二区图片欧美| 99国产精品一区二区三区| 一级毛片女人18水好多| 一本色道久久久久久精品综合| bbb黄色大片| 精品国产乱码久久久久久男人| 一个人免费看片子| 老司机亚洲免费影院| 人人妻人人添人人爽欧美一区卜| 高清av免费在线| 久久人妻熟女aⅴ| 在线观看一区二区三区激情| 又紧又爽又黄一区二区| 人妻久久中文字幕网| 大香蕉久久成人网| 午夜福利在线免费观看网站| 男女床上黄色一级片免费看| 国产高清视频在线播放一区| 午夜两性在线视频| 中文亚洲av片在线观看爽 | 一边摸一边抽搐一进一小说 | avwww免费| 嫩草影视91久久| 日本黄色视频三级网站网址 | 国产97色在线日韩免费| 国产熟女午夜一区二区三区| 久久亚洲真实| 久久中文字幕一级| 日韩视频一区二区在线观看| av视频免费观看在线观看| 免费看a级黄色片| 三上悠亚av全集在线观看| 国产精品久久久久久精品电影小说| 午夜91福利影院| 欧美乱码精品一区二区三区| 色尼玛亚洲综合影院| 国产亚洲欧美精品永久| 视频在线观看一区二区三区| 19禁男女啪啪无遮挡网站| 91字幕亚洲| 麻豆成人av在线观看| 久久人妻av系列| 精品欧美一区二区三区在线| 精品国内亚洲2022精品成人 | 午夜福利欧美成人| 999精品在线视频| 天堂俺去俺来也www色官网| 女人爽到高潮嗷嗷叫在线视频| 一级毛片精品| 成年女人毛片免费观看观看9 | 欧美午夜高清在线| 亚洲少妇的诱惑av| 国产在线一区二区三区精| 成人特级黄色片久久久久久久 | 国产精品久久久久久精品古装| 曰老女人黄片| 亚洲国产欧美网| 亚洲精品av麻豆狂野| 脱女人内裤的视频| 久久ye,这里只有精品| 中文字幕最新亚洲高清| av天堂久久9| 国产亚洲欧美精品永久| 国产成人啪精品午夜网站| 亚洲免费av在线视频| 国产熟女午夜一区二区三区| 久久精品成人免费网站| 一夜夜www| 狠狠狠狠99中文字幕| 亚洲久久久国产精品| av欧美777| 国产欧美日韩一区二区三| 国产在线观看jvid| 十八禁高潮呻吟视频| 999精品在线视频| 动漫黄色视频在线观看| 正在播放国产对白刺激| e午夜精品久久久久久久| 中亚洲国语对白在线视频| 欧美黑人欧美精品刺激| 午夜福利乱码中文字幕| 亚洲欧美色中文字幕在线| 丁香六月天网| 日本欧美视频一区| 成年动漫av网址| 欧美老熟妇乱子伦牲交| 18禁国产床啪视频网站| 精品久久久久久久毛片微露脸| 免费人妻精品一区二区三区视频| 交换朋友夫妻互换小说| 在线看a的网站| 男女免费视频国产| 亚洲七黄色美女视频| 久久久国产精品麻豆| 俄罗斯特黄特色一大片| 啦啦啦 在线观看视频| 久久久欧美国产精品| 亚洲精品美女久久久久99蜜臀| 热99久久久久精品小说推荐| 丁香欧美五月| 精品一品国产午夜福利视频| 国产高清激情床上av| 高清视频免费观看一区二区| 国产在线观看jvid| 欧美精品高潮呻吟av久久| 欧美精品高潮呻吟av久久| 精品一区二区三卡| tube8黄色片| 欧美黑人精品巨大| 日韩欧美一区二区三区在线观看 | 亚洲精品久久午夜乱码| 国产精品国产高清国产av | 国产高清国产精品国产三级| 久久久国产一区二区| 国产成人一区二区三区免费视频网站| 美女主播在线视频| 一级,二级,三级黄色视频| 亚洲精华国产精华精| 色精品久久人妻99蜜桃| 麻豆av在线久日| 在线观看免费日韩欧美大片| 中文欧美无线码| 欧美精品一区二区免费开放| bbb黄色大片| 日韩大片免费观看网站| 久久久精品94久久精品| 久久久久精品国产欧美久久久| 欧美大码av| 免费一级毛片在线播放高清视频 | 日韩大片免费观看网站| 国产一区二区在线观看av| www.精华液| 亚洲av成人一区二区三| 精品国产亚洲在线| 他把我摸到了高潮在线观看 | 国产麻豆69| 久久久久国产一级毛片高清牌| 搡老乐熟女国产| 精品国产乱子伦一区二区三区| 搡老熟女国产l中国老女人| 久久精品亚洲熟妇少妇任你| 757午夜福利合集在线观看| 国产精品亚洲一级av第二区| 黑人巨大精品欧美一区二区蜜桃| 久久人人97超碰香蕉20202| 一个人免费看片子| 91麻豆精品激情在线观看国产 | 国产熟女午夜一区二区三区| av福利片在线| 麻豆国产av国片精品| 亚洲色图 男人天堂 中文字幕| 国产色视频综合| 黄片大片在线免费观看| 日本av免费视频播放| 亚洲天堂av无毛| 国产成人av教育| 国产成人免费观看mmmm| 99热国产这里只有精品6| 首页视频小说图片口味搜索| av天堂久久9| 国产亚洲精品一区二区www | 欧美国产精品一级二级三级| 亚洲色图综合在线观看| 国产精品1区2区在线观看. | 国产成人一区二区三区免费视频网站| 欧美中文综合在线视频| av国产精品久久久久影院| 又黄又粗又硬又大视频| 建设人人有责人人尽责人人享有的| 蜜桃在线观看..| 岛国毛片在线播放| 免费看十八禁软件| bbb黄色大片| 两个人看的免费小视频| 国产精品久久久久久精品电影小说| 动漫黄色视频在线观看| 午夜福利欧美成人| 国产精品.久久久| 老鸭窝网址在线观看| 国产成人av教育| 搡老乐熟女国产| 又黄又粗又硬又大视频| 黑人巨大精品欧美一区二区蜜桃| 午夜福利一区二区在线看| 久久精品人人爽人人爽视色| 欧美在线一区亚洲| 国产成人av教育| 我的亚洲天堂| 亚洲中文字幕日韩| 亚洲欧美色中文字幕在线| 精品少妇黑人巨大在线播放| 久久精品国产亚洲av高清一级| av福利片在线| av天堂久久9| 亚洲国产欧美日韩在线播放| 亚洲av欧美aⅴ国产| 精品卡一卡二卡四卡免费| 亚洲成人国产一区在线观看| 国产熟女午夜一区二区三区| 啦啦啦在线免费观看视频4| 日本黄色视频三级网站网址 | 两个人看的免费小视频| 亚洲欧美一区二区三区久久| 又黄又粗又硬又大视频| 亚洲午夜精品一区,二区,三区| 久久人人97超碰香蕉20202| 99精品欧美一区二区三区四区| 久久亚洲真实| 欧美日韩视频精品一区| 国产人伦9x9x在线观看| 欧美日韩av久久| videos熟女内射| 国产欧美亚洲国产| 亚洲五月婷婷丁香| 亚洲avbb在线观看| 欧美激情久久久久久爽电影 | 欧美乱码精品一区二区三区| 国产成人欧美| 岛国毛片在线播放| 国产男女内射视频| 麻豆乱淫一区二区| 欧美激情 高清一区二区三区| 中国美女看黄片| 黑人猛操日本美女一级片| 女性被躁到高潮视频| 纯流量卡能插随身wifi吗| 日韩欧美三级三区| 久久亚洲真实| 欧美国产精品一级二级三级| 老司机亚洲免费影院| 欧美av亚洲av综合av国产av| 人人妻人人澡人人爽人人夜夜| 91国产中文字幕| 亚洲色图综合在线观看| 18禁观看日本| 国产一区二区三区在线臀色熟女 | 麻豆成人av在线观看| 久久久欧美国产精品| 免费看十八禁软件| 国产成人av激情在线播放| 日韩欧美国产一区二区入口| 天天操日日干夜夜撸| 少妇的丰满在线观看| 色婷婷久久久亚洲欧美| 国产精品久久久人人做人人爽| 国产精品成人在线| 啦啦啦视频在线资源免费观看| 国产在线一区二区三区精| 成年人免费黄色播放视频| 国产熟女午夜一区二区三区| 久久国产精品人妻蜜桃| 国产视频一区二区在线看| 国产精品免费一区二区三区在线 | 久热爱精品视频在线9| 国产成人系列免费观看| 国产成人精品在线电影| 久久九九热精品免费| 少妇的丰满在线观看| 精品久久久精品久久久| 亚洲全国av大片| 视频区欧美日本亚洲| 欧美人与性动交α欧美精品济南到| 国产精品电影一区二区三区 | 国产免费视频播放在线视频| av线在线观看网站| 成在线人永久免费视频| 一区二区av电影网| 一区二区三区国产精品乱码| 丁香六月天网| 亚洲一码二码三码区别大吗| 亚洲avbb在线观看| 女同久久另类99精品国产91| 十八禁网站网址无遮挡| 黄色片一级片一级黄色片| 午夜91福利影院| 国产亚洲精品第一综合不卡| 黄频高清免费视频| 日韩中文字幕视频在线看片| 免费日韩欧美在线观看| 成年动漫av网址| 一级,二级,三级黄色视频| 国产日韩一区二区三区精品不卡| 欧美 日韩 精品 国产| 91成年电影在线观看| 亚洲人成伊人成综合网2020| 亚洲午夜理论影院| 日韩成人在线观看一区二区三区| 国产精品久久久久久人妻精品电影 | 久久国产精品影院| 国产高清激情床上av| 亚洲精品自拍成人| 亚洲色图综合在线观看| 久久婷婷成人综合色麻豆| 国产成人精品在线电影| 男女高潮啪啪啪动态图| 国产精品99久久99久久久不卡| 欧美黄色片欧美黄色片| 极品人妻少妇av视频| av免费在线观看网站| 久久久久精品人妻al黑| 亚洲成国产人片在线观看| 黑人巨大精品欧美一区二区mp4| 高清视频免费观看一区二区| 成人永久免费在线观看视频 | 国产精品久久久av美女十八| 日韩视频一区二区在线观看| 国产精品欧美亚洲77777| 欧美日韩黄片免| 国产福利在线免费观看视频| 国产免费视频播放在线视频| 精品人妻在线不人妻| 久久香蕉激情| 国产精品99久久99久久久不卡| 在线天堂中文资源库| 成人国产av品久久久| 99久久人妻综合| 天堂中文最新版在线下载| 最近最新中文字幕大全电影3 | 18禁国产床啪视频网站| 热99re8久久精品国产| 国产精品免费视频内射| 视频区欧美日本亚洲| 国产精品免费大片| 亚洲国产欧美日韩在线播放| 俄罗斯特黄特色一大片| 亚洲欧洲日产国产| 成人av一区二区三区在线看| 满18在线观看网站| 黄色怎么调成土黄色| 巨乳人妻的诱惑在线观看| 欧美人与性动交α欧美精品济南到| 高清欧美精品videossex| 人人妻人人爽人人添夜夜欢视频| 亚洲色图 男人天堂 中文字幕| 成年人午夜在线观看视频| 欧美激情极品国产一区二区三区| 可以免费在线观看a视频的电影网站| 建设人人有责人人尽责人人享有的| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品国产综合久久久| 精品人妻在线不人妻| 99热国产这里只有精品6| 97人妻天天添夜夜摸| 成在线人永久免费视频| 91老司机精品| 久久久精品免费免费高清| 午夜老司机福利片| 蜜桃国产av成人99| 亚洲黑人精品在线| 亚洲精品av麻豆狂野| 一区福利在线观看| 18禁裸乳无遮挡动漫免费视频| 精品亚洲乱码少妇综合久久| 老司机午夜十八禁免费视频| 老司机靠b影院| 视频在线观看一区二区三区| 99久久99久久久精品蜜桃| 免费在线观看完整版高清| av又黄又爽大尺度在线免费看| 成人av一区二区三区在线看| 99国产精品一区二区蜜桃av | 久久久久视频综合| 国产精品久久久久久精品古装| 免费观看人在逋| 国产精品电影一区二区三区 | 极品少妇高潮喷水抽搐| 国产色视频综合| 欧美激情久久久久久爽电影 | 欧美黄色片欧美黄色片| 欧美成人免费av一区二区三区 | 啦啦啦中文免费视频观看日本| 美女午夜性视频免费| 精品一区二区三区四区五区乱码| 动漫黄色视频在线观看| 亚洲av日韩在线播放| 少妇猛男粗大的猛烈进出视频| 亚洲第一青青草原| 777米奇影视久久| av线在线观看网站| 69精品国产乱码久久久| 亚洲专区国产一区二区| 性色av乱码一区二区三区2| 中文字幕高清在线视频| a级片在线免费高清观看视频| 午夜激情久久久久久久| 波多野结衣av一区二区av| tocl精华| 亚洲成人免费电影在线观看| 国产伦人伦偷精品视频| 欧美日韩国产mv在线观看视频| 亚洲精品国产精品久久久不卡| 国产aⅴ精品一区二区三区波| 美女午夜性视频免费| aaaaa片日本免费| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区在线观看av| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 国产亚洲av高清不卡| 性色av乱码一区二区三区2| 一区二区三区国产精品乱码| 一本一本久久a久久精品综合妖精| 欧美另类亚洲清纯唯美| av天堂在线播放| 69精品国产乱码久久久| 精品福利观看| 国内毛片毛片毛片毛片毛片| 国产av国产精品国产| 亚洲第一青青草原| 国产激情久久老熟女| 高清视频免费观看一区二区| 三级毛片av免费| 亚洲精品中文字幕在线视频| 欧美乱妇无乱码| 日韩一区二区三区影片| 又大又爽又粗| 久久国产精品影院| 色精品久久人妻99蜜桃| 国产老妇伦熟女老妇高清| 国内毛片毛片毛片毛片毛片| 真人做人爱边吃奶动态| 亚洲精品国产精品久久久不卡| 真人做人爱边吃奶动态| 亚洲av日韩精品久久久久久密| 日日爽夜夜爽网站| 亚洲精品国产精品久久久不卡| 一级毛片女人18水好多| 亚洲精品美女久久av网站| 怎么达到女性高潮| 欧美黄色片欧美黄色片| 正在播放国产对白刺激| 国产91精品成人一区二区三区 | 亚洲成人免费电影在线观看| 国产黄频视频在线观看| 成人av一区二区三区在线看| 欧美日韩视频精品一区| 黑人巨大精品欧美一区二区蜜桃| 日韩大片免费观看网站| 久热这里只有精品99| www.熟女人妻精品国产| 国产精品久久久人人做人人爽| 黑人猛操日本美女一级片| 日韩欧美一区视频在线观看| 国产亚洲欧美精品永久| 久久国产精品男人的天堂亚洲| 日本撒尿小便嘘嘘汇集6| 午夜精品国产一区二区电影| 国产精品成人在线| 亚洲中文字幕日韩| 午夜福利在线观看吧| 老司机在亚洲福利影院| cao死你这个sao货| 欧美国产精品va在线观看不卡| 欧美激情高清一区二区三区| 中文字幕高清在线视频| 国产在线免费精品| 精品少妇久久久久久888优播|