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

    Kinetic Study of Esterification of Lactic Acid with Isobutanol and n-Butanol Catalyzed by Ion-exchange Resins*

    2009-05-15 06:17:36QuYixin屈一新PengShaojun彭少君WangShui王水ZhangZhiqiang張志強(qiáng)andWangJidong王際東
    關(guān)鍵詞:張志強(qiáng)王水瑪麗

    Qu Yixin (屈一新), Peng Shaojun (彭少君), Wang Shui (王水), Zhang Zhiqiang (張志強(qiáng)) and Wang Jidong (王際東)**

    ?

    Kinetic Study of Esterification of Lactic Acid with Isobutanol and-Butanol Catalyzed by Ion-exchange Resins*

    Qu Yixin (屈一新), Peng Shaojun (彭少君), Wang Shui (王水), Zhang Zhiqiang (張志強(qiáng)) and Wang Jidong (王際東)**

    Beijing Key Laboratory of Bioprocess, Department of Chemical Engineering, Beijing University of Chemical Technology, Beijing 100029, China

    The esterification reactions of lactic acid with isobutanol and-butanol have been studied in the presence of acid ion-exchange resin Weblyst D009. The influences of catalyst loading, stirrer speed, catalyst particle size, initial reactant molar ratio and temperature on the reaction rate have been examined. Experimental kinetic data were correlated by using the Pseudo-homogeneous, Langmuir-Hinshelwood and Eley-Rideal models. Nonideality of the liquid phase was taken into account by using activities instead of molar fractions. The activity coefficients were calculated according to the group contribution method UNIFAC. Provided that the nonideality of the liquid is taken into account, the esterification kinetics of lactic acid with isobutanol and-butanol catalyzed by the acid ion-exchange resin can be described using all three models with reasonable errors.

    kinetics, esterification, lactic acid, isobutanol,-butanol, ion exchange resin

    1 INTRODUCTION

    Lactic acid is used in wide applications in food industry for preservation and flavoring purposes, as well as in pharmaceutical and cosmetic industries. It is also used as an important raw material for the manufacture of biodegradable polymers [1, 2]. Its esters including isobutyl lactate and-butyl lactate are used as flavors and solvents with excellent properties [3]. The esterification reactions of lactic acid with alcohols are not only used in the production of the corresponding esters but also in the purification process of lactic acid. Therefore, kinetic data concerning esterification of lactic acid are necessary and important for design of these industrial processes.

    For the esterification of lactic acid with alcohols, both homogeneous and heterogeneous catalysts can be used. However, most of the recent studies use heterogeneous solid catalysts, like acid ion-exchange resins, which can avoid the drawbacks of homogeneous catalyst, such as equipment corrosion and side reactions [4-10]. Dassy. [11] studied the kinetics of liquid phase synthesis and hydrolysis of butyl lactate catalyzed by acid ion-exchange resin in dioxane and toluene in a batch reactor. They found that the reaction rate was first order with respect to the amount of catalyst and the concentration of acid without considering the nonideality of the liquid mixture. Kumar and Mahajani [12] described the esterification of lactic acid with-butanol by reactive distillation catalyzed by acid ion-exchange resin. They used pseudo-homogeneousmodel to simulate the kinetic data also without considering the nonideality of the liquid mixture. Sanz. [9] investigated the reaction of lactic acid with methanol using Amberlyst 15 as catalyst. Three models including the Pseudo-homogeneous (PH), the Langmuir-Hinshelwood (LH) and the Eley-Rideal (ER) were used to correlate the experimental data. The PH model was found to represent the experimental data fairly well. Zhang. [10] studied the esterification of lactic acid [80% (by mass)] with ethanol in the presence of five acid ion-exchange resins. The LH model based on the selective adsorption of water and ethanol on the catalyst was found to be a more appropriate model to describe the kinetic behavior of these systems.

    In this work, the esterification reactions of lactic acid with isobutanol and-butanol catalyzed by acid ion-exchange resins, Weblyst D009 and Weblyst D80 have been investigated. To provide a general kinetic model, experimental data were regressed by using different models based on homogeneous and heterogeneous approaches. To take into the nonideality of the liquid mixture, activities of the species in liquid phase, estimated by the group contribution method UNIFAC, were used instead of molar fractions.

    2 EXPERIMENTAL

    2.1 Chemical reagents

    Isobutanol [purity≥99.5% by mass)] and-butanol [purity≥99.7% (by mass)] were supplied by Hongyan Reagent Co., Tianjin. Lactic acid [purity≥85% (by mass)] was purchased from Xilong Chemical Co., Guangdong. All of them were used without further purification. Lactic acid has two functional groups, a carboxyl and a hydroxyl. When the esterification reaction is carried out in high concentration of lactic acid, polylactic acid will be formed due to self-esterification [13]. Therefore, in this study dilute aqueous solution of lactic acid [20% (by mass)] was used in order to avoid the formation of polymers. Dilute lactic acid was obtained by adding distilled water into the lactic acid [85% (by mass)]. After mixing, the solutions were heated at 353.15 K for 1 week to increase the rate of formation of various oligomers of lactic acid [14]. The amount of polymerized lactic acid in the dilute lactic acid was considered negligible after determined by back titration [14].

    2.2 Catalysts

    Acid ion-exchange resins used in this study were purchased from Xiqiao Organic Material Co., Beijing. The properties of the catalysts obtained from the company are in Table 1. Prior to using for esterification reactions, the catalysts were first washed with methanol, then with deionized water and finally dried at 373 K for 24 h to completely remove moisture.

    2.3 Apparatus and procedure

    Esterification reactions of lactic acid with isobutanol and with-butanol were carried out in a four-necked flask of 500 ml capacity fitted with a reflux condenser, a thermometer and a magnetic stirrer. The temperature was maintained within an accuracy of ±0.5 K by an electric-heated thermostatic oil bath. The aqueous lactic acid solution and the catalyst were first charged into the reactor and heated to the desired temperature. Then isobutanol (or-butanol) was fed into the reactor. This time was considered to be the starting point of the reaction. Samples were withdrawn at regular time intervals for analysis. The catalyst loading was expressed as the mass ratio of the dry catalyst to the pure lactic acid in this article.

    2.4 Analysis

    All samples were weighed accurately by using an electronic balance with an accuracy of 0.0001 g. The amount of lactic acid was determined by titration with a standard sodium hydroxide solution using phenolphthalein as an indicator. Parallel tests indicated that the average error of the titration method was less than 2%. The samples were also analyzed with a gas chromatograph (GC) equipped a SW-FFAP capillary column (30 m×0.32 mm×0.5 μm) and a flame ionization detector (FID). Nitrogen with a purity of 99.99% was used as the carrier gas. Temperature of the injector and the detector were kept at 473.15 K. Initial temperature of the column was 323.15 K. After 5 min, the column temperature was increased to 473.15 K with a rate of 15 K·min-1. Finally it was maintained at 473.15 K for 20 min. No significant differences were observed between the results obtained by GC and titration method.

    3 KINETIC MODELING

    Esterification reactions catalyzed by ion-exchange resins can be described using different kinetic models based on homogeneous and heterogeneous approaches. Although PH model does not take into account the adsorption effect of the species in the reactant medium, it has been successfully used in high polar reaction media [7, 9]. The LH and the ER models both include the adsorption effects of the species in the reactant medium. It is obvious that these two models are more complicated than the PH model. The basic assumption of LH model is that all reactants are adsorbed on the catalyst surface before chemical reactions occur. The ER model assumed that the reaction takes place between adsorbed and non-adsorbed reactants. A general kinetic expression for all the three models is written as

    To consider the nonideality of the liquid phase, activities of components were used instead of molar fractions. The activity coefficients were calculated by using group contribution method UNIFAC [15-20]. The equilibrium constants for the two esterification reactions were calculated from the component concentrations at the equilibrium through Eq. (3):

    Table 1 Physical properties of the ion-exchange resins

    wherex,γare the molar fraction and activity coefficient of theth component, respectively, at equilibrium.

    The kinetic parameters of the models were obtained by minimizing the sum of squared residuals between the experimental and calculated reaction rates as shown in Eq. (4) through the least square method:

    A fourth order Runge-Kutta method was used for integrating numerically the differential equations which describe the kinetic model with the previously determined parameters. The quality of the fit was estimated by the mean relative deviation of the molar fraction of lactic acid:

    The subscripts “exp” and “cal” in Eqs. (4) and (5) are referring to the experimental and the calculated results andfor the total number of samples.

    4 RESULTS AND DISCUSSION

    4.1 Catalyst performance

    At the beginning, two acid ion-exchange resins, Weblyst D009 and Weblyst D80, were selected as the candidate catalysts. Their catalytic activities were first tested. Fig. 1 shows the comparison of the conversion of lactic acid obtained with Weblyst D009 and Weblyst D80 in the esterification reactions of lactic acid with isobutanol and-butanol. Results presented in Fig. 1 indicate that both Weblyst D009 and Weblyst D80 are effective catalysts.

    isobutanol:□?Weblyst D009;○?Weblyst D80

    -butanol:■?Weblyst D009;●?Weblyst D80

    Although both catalysts have approximately the sameeqvalues, Weblyst D009 shows a noticeable higher activity than Weblyst D80 when the conversion of lactic acid is beyond 10%. During the experiments with Weblyst D80, it was noticed that in the reaction medium a larger part of Weblyst D80 cracked into small particles. This was not observed for Weblyst D009. The reasons for this difference are likely resulted from the differences in their compositions and/or structures. This means that besides the meq value, the composition and structure of an ion-exchange resin are also important factors in determining its catalytic effect. As a consequence, Weblyst D009 was chosen as the catalyst for further studies.

    4.2 Elimination of mass transfer resistance

    To evaluate the external mass-transfer resistance, the esterification reactions were carried out at different stirrer speeds of 200, 300, 400 r·min-1while keeping the other conditions unchanged. The results obtained were shown in Fig. 2. It can be seen that there was no significant difference in the conversion of lactic acid when the stirrer speeds are 300 and 400 r·min-1for both isobutanol and-butanol. This conclusion agrees with the previous works [21-23] which indicated that external diffusion does not control the overall reaction rate unless the stirrer speed is very low or the viscosity of reactant mixture is very high. To ensure the absence of external mass transfer resistance, the kinetic experiments in this work were all performed at a stirrer speed of 300 r·min-1.

    isobutanol:●?200 r·min-1;■?300 r·min-1;▲?400 r·min-1

    -butanol:○?200 r·min-1;□?300 r·min-1;△?400 r·min-1

    To evaluate the internal diffusion, the Weblyst D009 was screened into particles with three different size ranges, 0.8 mm-1.25 mm, 0.65 mm-0.8 mm, 0.3 mm-0.65 mm and esterification reactions with these three types of particles were carried out. No significant differences in the conversion of lactic acid were found for both esterification reactions with isobutanol and-butanol. This indicates that intra-particle diffusion resistance is negligible when the particle size is smaller than 1.25 mm for the esterification under the present reaction conditions. This is consistent with the previous studies [24, 25], which indicate that intraparticle diffusion resistances of the reactant in the ion exchange resin are not important.

    4.3 Effect of temperature

    In order to investigate the effect of temperature, both esterification reactions were carried out in the temperature range from 333.15 to 363.15 K. Typical results are shown in Figs. 3 and 4. It can be seen that the reaction rate increases substantially with the increasing temperature. However, the equilibrium conversion was nearly equal in the range of temperatures studied in this work. Esterification reactions of acids with alcohols occur generally with small absolute values of reaction enthalpies, so the equilibrium conversion of lactic acid is a week function of reaction temperature. The similar effect of reaction temperature on the equilibrium conversion of lactic acid has been observed by Delgado. [26] in the reaction of lactic acid with ethanol and Sanz. [9] in the reaction of lactic acid with methanol.

    ▼?333.15 K;▲?343.15 K;●?353.15 K;■?363.15 K

    ▼?333.15 K;▲?343.15 K;■?353.15 K;●?363.15 K

    From the equilibrium concentration of the reaction products, the equilibrium constants of the reactions in this work can be calculated according to

    4.4 Effect of catalyst loading

    The catalyst mass loading was varied from 1.5% to 6% for the both esterification reactions. The conversions of lactic acid as a function of time with different catalyst loadings are shown in Fig. 6 for isobutanol at 353.15 K and Fig. 7 for-butanol at 363.15 K. As it can be seen from these figures, with increasing catalyst loading the conversion of lactic acid increases. A higher loading of catalyst results in reduction of the time required to reach the reaction equilibrium. No significant change in the equilibrium conversion of lacticacid is observed when the catalyst loading is changed.

    Table 2 List of parameters of Eq. (1)

    ■-butanol;□?isobutanol

    catalyst mass loading: ▲?6%;■?3%;●?1.5%

    catalyst mass loading:▲?6%;■?3%;●?1.5%

    4.5 Effect of initial reactant molar ratio

    To investigate the effect of the initial molar ratio of alcohols to lactic acid on the conversion of lactic acid, the initial molar ratio was changed from 1 to 4 for isobutanol (ISOH/HL) and 1 to 3 for-butanol (BUOH/HL). The results are given in Figs. 8 and 9. It can be seen that the conversion of lactic acid increases with increasing the initial molar ratio of alcohols to lactic acid. For the esterification with isobutanol, the equilibrium conversion of lactic acid increases from 15% to 40% when the initial molar ratioISOH/HLis raised from 1 to 4. For the esterification with-butanol, the equilibrium conversion of lactic acid increases from 15% to 35% when the initial molar ratioBUOH/HLis raised from 1 to 3. The results observed for the effect of the initial molar ratio of alcohols to lactic acid indicated that the equilibrium conversion of lactic acid can be effectively enhanced by using a large excess of alcohols, which is consistent with the results of Delgado. [26] who studied the esterification reaction of lactic acid with ethanol.

    initial reactant molar ratioISOH/HL:▲?4︰1; ■?3︰1;●?2︰1;▼?1︰1

    initial reactant molar ratioBUOH/HL:▲?3︰1; ■?2︰1;●?1︰1

    4.6 Kinetic models

    From the comparison of the values ofRSandRD, it is noticed that the assumption of ideal behavior of the liquid phase results in large errors. This is expected since the reaction systems are high polarity medium. When the nonideality of the liquid phase is taken into account, the PH model gives comparable accuracy as the ER and the LH model. The ER model and LH model give very close values ofRSandRD, implying that the assumptions made for the two models are reasonable. Provided that the nonideality of liquid phase is taken into account, it is difficult to discriminate one model from the others solely based on the values ofRSandRDsince the differences between the values ofRSand ofRDfor different models are very small.

    The Arrhenius plots for the esterification of lactic acid with isobutanol and with-butanol are given in Fig. 10. The continuous line represents the results obtained with LH model for isobutanol and the dash line represents the results obtained with the PH model for-butanol. A good linear relation between lneand 1/is observed.

    然而早上7點(diǎn),瑪麗的父母卻發(fā)現(xiàn)孩子倒在衛(wèi)生間里,已經(jīng)停止了呼吸。醫(yī)生對悲傷的父母說,瑪麗可能是倒下時頭部受到撞擊,導(dǎo)致猝死。

    Figure 10 Arrhenius plot for the esterification of lactic acid with isobutanol and with-butanol (The continuous line represents the LH model, the dash line for the PH model.)

    ■?isobutanol;□-butanol

    Table 3 Parameters of the kinetic models for the ion-exchange catalyzed esterification of lactic acid with isobutanol and n-butanol

    Table 4 Comparison of the activation energies for the esterification of lactic acid with isobutanol and n-butanol

    A comparison of the activation energies obtained in this work with those found in references is given in Table 4. For the esterification of lactic acid with isobutanol catalyzed by ion-exchange resin, no value of the activation energy was found in the literature. The activation energy found in this work is about 58 kJ·mol-1for the ER and the LH model and 65.6 kJ·mol-1for the PH model. For the esterification of lactic acid with butanol, the activation energy found in this work is about 52 kJ·mol-1for the ER and the LH model and 54.3 kJ·mol-1for the PH model. These valuesare 4-6 kJ·mol-1higher than that of Dassy. [11]who used Amberlyst 15 as the catalyst and dioxane and toluene as the solvents.

    As compared to the previous work [11], the esterification of lactic acid with- and-butanol in aqueous solution was investigated using an acid ion-exchange resin Weblyst D009. Three kinetic models that take into account of the nonideality of the reaction medium were used to simulate the experimental data. When the nonideality of the reaction medium was considered, more accurate modeling results were obtained. Although all the three models that include the nonideality of the reaction medium give practically the same accuracy, it is prefer to use the PH model since this model has a simple form and is easy to use.

    5 CONCLUSIONS

    In this study, the kinetics of the esterification of lactic acid with isobutanol and-butanol catalyzed by acid ion-exchange resin Weblyst D009 was experimentally investigated. The reaction rate was found to increase with increasing reaction temperature and catalyst loading. The equilibrium conversion of lactic acid increased with increasing the molar ratio of the alcohols to lactic acid. The dependence of the equilibrium constant on the reaction temperature was in line with van’t Hoff equation. Three models, PH, LH and ER have been used to simulate the experimental data. To correct the nonideality of the liquid phase, activities by the group contribution method UNIFAC instead of molar fractions were used during the simulations. It was found that all these models were able to describe the kinetics of the esterification of lactic acid with butanol and isobutanol without high errors and gave comparable accuracy when the nonideality of the liquid phase was taken into account. The PH model is preferred due to its simple mathematical form.

    NOMENCLATURE

    activity

    A,eapparent activation energy, kJ·mol-1

    eqequilibrium constant

    adsorption coefficient

    eforward reaction rate constant, mol·g-1·min-1

    RDmean relative deviation

    number of experimental data

    catmass of catalyst, g

    nmole of component

    gas constant , kJ·mol-1·K-1

    reaction rate, mol·(g cat)-1·min-1

    RSminimum sum

    absolute temperature, K

    time, min

    molar fraction

    activity coefficient

    stoichiometric coefficient

    Subscripts

    BL butyl lactate

    BUOH-butanol

    calc calculated values

    e esterification reaction

    eq equilibrium

    exp experimental values

    HL lactic acid

    IL isobutyl lactate

    ISOH isobutanol

    ,components

    L isobutyl lactate or butyl lactate

    OH isobutanol or butanol

    W water

    1 Datta, R., Tsai, S.P., Bonsignore, P., Moon, S.H., Frank, J.R., “Technological and economic potential of poly (lactic acid) and lactic acid derivatives”,.., 16, 221-231 (1995).

    2 Karin, H., Bflrbel, H., “Factors affecting the fermentative lactic acid production from renewable resources”,.., 26, 87-107 (2000).

    3 Light Industry Research Institute in Jinan, Manual on Synthesis of Flavor, Light Industry Pub. Co., Beijing (1985). (in Chinese)

    4 Yadav, G.D., Thathagar, M.B., “Esterification of maleic acid with ethanol over cation-exchange resin catalysts”,..., 52, 99-110 (2002).

    5 Yadav, G.D., Thorat, T.S., “Kinetics of alkylation of-cresol with isobutylene catalyzed by sulfated zirconia”,...., 35, 721-732 (1996).

    6 Peters, T.A., Benes, N.E., Holmen, A., Keurentjes, J.T.F., “Comparison of commercial solid acid catalysts for the esterification of acetic acid with butanol”,..., 297, 182-188 (2006).

    7 Sanz, M.T., Murga, R., Beltran, S., Cabezas, J.L., Coca, J., “Autocatalyzed and ion-exchange-resin-catalyzed esterification kinetics of lactic acid with methanol”,...., 41, 512-517 (2002).

    8 Troupe, R.A., Dimilla, E., “Kinetics of the ethyl alcohol-lactic acid reaction”,..., 49, 847-855 (1957).

    9 Sanz, M.T., Murga, R., Beltran, S., Cabezas, J.L., Coca, J., “Kinetic study for the reactive system of lactic acid esterification with methanol: Methyl lactate hydrolysis reaction”,...., 43, 2049-2053 (2004).

    10 Zhang, Y., Ma, L., Yang, J., “Kinetics of esterification of lactic acid with ethanol catalyzed by cation-exchange resins”,..., 61, 101-114 (2004).

    11 Dassy, S., Wiame, H., Thyrion, F.C., “Kinetics of the liquid phase synthesis and hydrolysis of butyl lactate catalysed by cation-exchange resin”,...., 59, 149-156 (1994).

    12 Kumar, R., Mahajani, S.M., “Esterification of lactic acid with-butanol by reactive distillation”,...., 46, 6873-6882 (2007).

    13 Holten, C.H., Lactic Acid: Properties and Chemistry of Lactic Acid and Derivatives, Verlag Chemie, Copenhagen, Copenhagen (1971).

    14 Vu, D.T., Kolah, A.K., Asthana, N.S., “Oligomer distribution in concentrated lactic acid solutions”,., 236, 125-135 (2005).

    15 Rehfinger, A., Hoffman, U., “Kinetics of methyl tertiary butyl ether liquid phase synthesis catalyzed by ion exchange resin (I) Intrinsic rate expression in liquid phase activities”,..., 45, 1605-1607 (1990).

    16 Ladisch, M., Westgate, P., Hendrickson, R., Brewer, M., “Framework for correlating composition dependent equilibrium conversion in methyl-butyl ether formation by ion-exchange catalysts”,...., 34, 2811-2816 (1995).

    17 Rihko, L.K., Krause, A.O., “Kinetics of heterogeneously catalyzed-amyl methyl ether reactions in the liquid phase”,...., 34, 1172-1180 (1995).

    18 Rihko, L.K., Kiviranta-P??kk?nen, P.K., Krause, A.O., “Kinetic model for the etherification of isoamylenes with methanol”,...., 36, 614-621 (1997).

    19 Lee, M.J., Wu, H.T., Lin, H.M., “Kinetics of catalytic esterification of acetic acid and amyl alcohol over Dowex”,...., 39, 4094-4099 (2000).

    20 Reid, R.C., Prausnitz, J.M., Poling, B.E., The Properties of Gases and Liquids, McGraw-Hill, New York (1987).

    21 Chakrabarti, A., Sharma, M.M., “Cation ion exchange resin as catalyst”,.., 20, 1-45 (1993).

    22 Yu, W., Hidajat, K., Ray, A.K., “Determination of adsorption and kinetic parameters for methyl acetate esterification and hydrolysis reaction catalyzed by Amberlyst 15”,..., 260, 191-205 (2004).

    23 Darge, O., Thyrion, F.C., “Kinetics of the liquid phase esterification of acrylic acid with butanol catalysed by cation exchange resin”,...., 58, 351-355 (1993).

    24 Seo, Y., Hong, W.H., “Kinetics of esterification of lactic acid with methanol in the presence of cation exchange resin using a pseudo-homogeneous model”,..., 33, 128-133 (2000).

    25 Teo, H.T.R., Saha, B., “Heterogeneous catalysed esterification of acetic acid with isoamyl alcohol: Kinetic studies”,.., 228, 174-182 (2004).

    26 Delgado, P., Sanz, M.T., Beltran, S., “Kinetic study for esterification of lactic acid with ethanol and hydrolysis of ethyl lactate using an ion-exchange resin catalyst”,..., 126, 111-118 (2007).

    2009-03-19,

    2009-07-12.

    the National Basic Research Program of China (2007CB714300).

    ** To whom correspondence should be addressed. E-mail: jidongwang1963@yahoo.com.cn

    猜你喜歡
    張志強(qiáng)王水瑪麗
    Origin of tradeoff between movement velocity and attachment duration of kinesin motor on a microtubule
    集合中的易錯問題剖析
    難忘的八個字
    和我一起玩
    這才叫大藏家——張志強(qiáng)先生和他的精絕藏石
    寶藏(2018年12期)2019-01-29 01:50:50
    捉“妖”記(二)——讓“數(shù)列妖”無處藏身
    瑪麗的怪獸美容店
    冷暖
    冷暖
    短篇小說(2015年10期)2015-11-03 04:22:31
    志強(qiáng)的石
    中華奇石(2014年12期)2014-07-09 18:30:22
    av线在线观看网站| 狠狠婷婷综合久久久久久88av| 日本av手机在线免费观看| 午夜福利免费观看在线| 国产在线视频一区二区| 亚洲精品一二三| 欧美精品亚洲一区二区| 赤兔流量卡办理| 性色av一级| 在现免费观看毛片| 女性被躁到高潮视频| 亚洲成人手机| 免费av中文字幕在线| 91精品国产国语对白视频| 亚洲色图综合在线观看| av在线播放精品| 久久精品国产综合久久久| 久久久久国产精品人妻一区二区| 午夜91福利影院| 美国免费a级毛片| 亚洲,一卡二卡三卡| 色综合欧美亚洲国产小说| 日韩大片免费观看网站| 啦啦啦啦在线视频资源| 丰满乱子伦码专区| 在线看a的网站| 久久女婷五月综合色啪小说| 不卡视频在线观看欧美| 岛国毛片在线播放| 90打野战视频偷拍视频| 肉色欧美久久久久久久蜜桃| 天堂中文最新版在线下载| 日韩av不卡免费在线播放| 国产精品一国产av| 99热网站在线观看| 国产男人的电影天堂91| 国产人伦9x9x在线观看| 精品国产国语对白av| av国产精品久久久久影院| 一边亲一边摸免费视频| 日韩成人av中文字幕在线观看| 在线免费观看不下载黄p国产| 男人添女人高潮全过程视频| 国产精品 国内视频| 国产成人啪精品午夜网站| 免费在线观看视频国产中文字幕亚洲 | 在线观看人妻少妇| 免费观看人在逋| 十分钟在线观看高清视频www| 99久久99久久久精品蜜桃| 午夜精品国产一区二区电影| 亚洲成色77777| bbb黄色大片| 久久97久久精品| 丝袜在线中文字幕| 五月开心婷婷网| 美女主播在线视频| 精品卡一卡二卡四卡免费| 丝瓜视频免费看黄片| av有码第一页| 少妇人妻久久综合中文| 日韩免费高清中文字幕av| 女人高潮潮喷娇喘18禁视频| 51午夜福利影视在线观看| 一级毛片 在线播放| 亚洲欧美一区二区三区久久| 免费黄频网站在线观看国产| 纵有疾风起免费观看全集完整版| 男女高潮啪啪啪动态图| 久久久久久免费高清国产稀缺| 亚洲第一青青草原| 99香蕉大伊视频| 看免费av毛片| 日韩精品有码人妻一区| 亚洲欧洲国产日韩| 久久精品国产综合久久久| 涩涩av久久男人的天堂| 免费女性裸体啪啪无遮挡网站| 色吧在线观看| 午夜激情av网站| 午夜福利乱码中文字幕| 国产 一区精品| 国产99久久九九免费精品| 国产一级毛片在线| www.av在线官网国产| 精品国产乱码久久久久久小说| 久久97久久精品| 嫩草影视91久久| 免费高清在线观看日韩| 七月丁香在线播放| 欧美日韩av久久| 久热这里只有精品99| 最近最新中文字幕大全免费视频 | 99久久99久久久精品蜜桃| 男女午夜视频在线观看| 中文字幕人妻丝袜制服| 男女高潮啪啪啪动态图| 亚洲av综合色区一区| 老汉色av国产亚洲站长工具| bbb黄色大片| 人妻人人澡人人爽人人| 久久精品aⅴ一区二区三区四区| 日韩一本色道免费dvd| 欧美少妇被猛烈插入视频| 免费黄色在线免费观看| 黄色视频不卡| 国产成人欧美在线观看 | 一区福利在线观看| 黄片无遮挡物在线观看| 免费观看性生交大片5| 中文字幕亚洲精品专区| 日韩熟女老妇一区二区性免费视频| 日本vs欧美在线观看视频| 午夜福利影视在线免费观看| 考比视频在线观看| 美女主播在线视频| 在现免费观看毛片| 亚洲在久久综合| 别揉我奶头~嗯~啊~动态视频 | 人人澡人人妻人| 999精品在线视频| 成人亚洲精品一区在线观看| 国产亚洲欧美精品永久| 在线免费观看不下载黄p国产| 欧美黑人欧美精品刺激| 亚洲av成人精品一二三区| 中文字幕高清在线视频| 国产精品香港三级国产av潘金莲 | 国产片内射在线| 一边摸一边抽搐一进一出视频| 秋霞伦理黄片| 99精品久久久久人妻精品| 久久 成人 亚洲| 伊人久久国产一区二区| 久久久亚洲精品成人影院| 久久影院123| 丝袜在线中文字幕| 国产福利在线免费观看视频| 国产亚洲最大av| 蜜桃国产av成人99| 肉色欧美久久久久久久蜜桃| 中国国产av一级| 久久久久久人人人人人| 男女免费视频国产| 我的亚洲天堂| tube8黄色片| 亚洲国产欧美一区二区综合| 啦啦啦在线免费观看视频4| 如何舔出高潮| 国产男女超爽视频在线观看| 国产一区亚洲一区在线观看| 久久人人爽人人片av| 在线精品无人区一区二区三| 女人精品久久久久毛片| 侵犯人妻中文字幕一二三四区| 欧美黑人精品巨大| av福利片在线| 丝袜脚勾引网站| 尾随美女入室| 一级毛片 在线播放| 黄色视频在线播放观看不卡| 日韩,欧美,国产一区二区三区| 国产极品粉嫩免费观看在线| 午夜激情av网站| 捣出白浆h1v1| 99久久精品国产亚洲精品| 中文乱码字字幕精品一区二区三区| 国产av一区二区精品久久| 香蕉国产在线看| 亚洲婷婷狠狠爱综合网| 亚洲欧美成人精品一区二区| 女人精品久久久久毛片| 亚洲欧美激情在线| 性高湖久久久久久久久免费观看| 夫妻午夜视频| 亚洲精品国产av蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 国产精品女同一区二区软件| 涩涩av久久男人的天堂| 亚洲色图 男人天堂 中文字幕| 国产免费视频播放在线视频| 久久久精品免费免费高清| 最黄视频免费看| 在线 av 中文字幕| 国产成人精品在线电影| 黑人欧美特级aaaaaa片| 久久精品国产a三级三级三级| 大话2 男鬼变身卡| 亚洲精品av麻豆狂野| 少妇精品久久久久久久| 日本午夜av视频| 我要看黄色一级片免费的| 另类亚洲欧美激情| 午夜老司机福利片| 国产亚洲av高清不卡| 日韩熟女老妇一区二区性免费视频| 一级片'在线观看视频| 亚洲精品在线美女| 久久精品国产亚洲av高清一级| 精品一区二区三区av网在线观看 | 日韩不卡一区二区三区视频在线| 免费日韩欧美在线观看| 久久久欧美国产精品| 成年美女黄网站色视频大全免费| 99久久综合免费| 日韩中文字幕欧美一区二区 | 亚洲男人天堂网一区| 精品国产乱码久久久久久小说| www日本在线高清视频| 久久av网站| 亚洲伊人久久精品综合| 在线观看www视频免费| 看免费成人av毛片| 激情五月婷婷亚洲| 91成人精品电影| 汤姆久久久久久久影院中文字幕| 99久久99久久久精品蜜桃| 久久毛片免费看一区二区三区| 精品国产乱码久久久久久小说| 下体分泌物呈黄色| 亚洲中文av在线| 午夜久久久在线观看| 纯流量卡能插随身wifi吗| 国产成人精品久久二区二区91 | 宅男免费午夜| 免费不卡黄色视频| 一边摸一边做爽爽视频免费| 亚洲国产精品一区二区三区在线| 丁香六月欧美| 熟女av电影| 久久久久视频综合| 嫩草影视91久久| 汤姆久久久久久久影院中文字幕| av女优亚洲男人天堂| 美女高潮到喷水免费观看| 日本猛色少妇xxxxx猛交久久| 自线自在国产av| av网站免费在线观看视频| 久久97久久精品| av视频免费观看在线观看| 亚洲欧美成人综合另类久久久| 97人妻天天添夜夜摸| 亚洲av电影在线进入| 成人影院久久| 男人添女人高潮全过程视频| 亚洲一级一片aⅴ在线观看| 国产亚洲一区二区精品| 男的添女的下面高潮视频| 国产1区2区3区精品| 考比视频在线观看| 国产亚洲欧美精品永久| 日本91视频免费播放| 亚洲一卡2卡3卡4卡5卡精品中文| 麻豆精品久久久久久蜜桃| 亚洲欧美精品自产自拍| 在线观看免费高清a一片| 欧美变态另类bdsm刘玥| 亚洲国产欧美在线一区| 男女午夜视频在线观看| 国产一区二区在线观看av| 免费看av在线观看网站| 欧美日韩成人在线一区二区| 久久久久国产精品人妻一区二区| 精品福利永久在线观看| 欧美人与性动交α欧美软件| 国产免费现黄频在线看| 看免费av毛片| 久久精品国产亚洲av涩爱| 欧美国产精品一级二级三级| 一级爰片在线观看| 精品国产乱码久久久久久小说| 最近手机中文字幕大全| 午夜精品国产一区二区电影| 午夜91福利影院| 免费在线观看视频国产中文字幕亚洲 | 国产精品嫩草影院av在线观看| 三上悠亚av全集在线观看| 九草在线视频观看| 色综合欧美亚洲国产小说| 精品国产超薄肉色丝袜足j| 国产av精品麻豆| 久久久国产一区二区| 老司机在亚洲福利影院| 新久久久久国产一级毛片| 夫妻性生交免费视频一级片| 99热国产这里只有精品6| 99久久99久久久精品蜜桃| 国产探花极品一区二区| 免费av中文字幕在线| 国产精品麻豆人妻色哟哟久久| 欧美日韩一级在线毛片| 亚洲av综合色区一区| 视频在线观看一区二区三区| 久久久久精品性色| 国产免费福利视频在线观看| 99久国产av精品国产电影| 五月开心婷婷网| 日韩成人av中文字幕在线观看| 国产又色又爽无遮挡免| 久久久久精品国产欧美久久久 | www日本在线高清视频| 夫妻性生交免费视频一级片| 大片免费播放器 马上看| 亚洲精品美女久久av网站| 男女高潮啪啪啪动态图| 久热这里只有精品99| av又黄又爽大尺度在线免费看| 伊人久久大香线蕉亚洲五| 久久精品亚洲av国产电影网| 侵犯人妻中文字幕一二三四区| 国产福利在线免费观看视频| 又黄又粗又硬又大视频| 美女脱内裤让男人舔精品视频| 在线观看免费高清a一片| 亚洲精品美女久久久久99蜜臀 | 婷婷成人精品国产| 日韩人妻精品一区2区三区| 午夜91福利影院| 久久久久久久久久久免费av| 亚洲成人手机| 成人亚洲欧美一区二区av| 午夜福利视频精品| 人人妻人人澡人人爽人人夜夜| 色视频在线一区二区三区| a级片在线免费高清观看视频| 欧美中文综合在线视频| 国产淫语在线视频| 国产精品三级大全| h视频一区二区三区| 亚洲精品乱久久久久久| 日本欧美视频一区| 精品午夜福利在线看| 亚洲少妇的诱惑av| 国产精品亚洲av一区麻豆 | 丰满乱子伦码专区| 久久天堂一区二区三区四区| 男女边吃奶边做爰视频| 国产色婷婷99| 2018国产大陆天天弄谢| 无遮挡黄片免费观看| 久久这里只有精品19| 飞空精品影院首页| 免费在线观看黄色视频的| 亚洲国产欧美日韩在线播放| 波多野结衣一区麻豆| 久久久久久免费高清国产稀缺| 久久99精品国语久久久| 婷婷色av中文字幕| 香蕉国产在线看| 亚洲精品在线美女| 观看av在线不卡| 在线观看www视频免费| 久久99精品国语久久久| 亚洲精品国产区一区二| 看非洲黑人一级黄片| 中文字幕亚洲精品专区| 日本wwww免费看| 天天躁夜夜躁狠狠躁躁| 国产成人91sexporn| www.自偷自拍.com| 少妇人妻久久综合中文| 少妇精品久久久久久久| 久久久久精品久久久久真实原创| 国产黄色视频一区二区在线观看| 欧美成人精品欧美一级黄| 最近的中文字幕免费完整| 久久久欧美国产精品| 精品一区二区免费观看| 中文乱码字字幕精品一区二区三区| 亚洲第一青青草原| 国产精品免费大片| 免费看不卡的av| 男女高潮啪啪啪动态图| 人妻一区二区av| 久久热在线av| 精品少妇黑人巨大在线播放| 亚洲 欧美一区二区三区| 九草在线视频观看| 午夜91福利影院| 精品国产一区二区三区久久久樱花| 毛片一级片免费看久久久久| 精品酒店卫生间| 中文字幕高清在线视频| 国产有黄有色有爽视频| 老司机深夜福利视频在线观看 | 赤兔流量卡办理| 国产人伦9x9x在线观看| 十八禁网站网址无遮挡| 嫩草影院入口| 人体艺术视频欧美日本| av不卡在线播放| 国产av精品麻豆| 亚洲自偷自拍图片 自拍| 国产有黄有色有爽视频| 亚洲国产精品成人久久小说| 日韩大码丰满熟妇| 亚洲成人一二三区av| 香蕉国产在线看| 女人被躁到高潮嗷嗷叫费观| 少妇人妻 视频| 亚洲少妇的诱惑av| 又粗又硬又长又爽又黄的视频| 一级爰片在线观看| 精品一区在线观看国产| 免费在线观看完整版高清| 日本av手机在线免费观看| 国产成人精品福利久久| 9191精品国产免费久久| 肉色欧美久久久久久久蜜桃| 国产精品99久久99久久久不卡 | 乱人伦中国视频| 久久久精品国产亚洲av高清涩受| 亚洲欧美色中文字幕在线| 在线观看三级黄色| 亚洲视频免费观看视频| 观看美女的网站| 老汉色av国产亚洲站长工具| 9191精品国产免费久久| 天堂中文最新版在线下载| 亚洲欧美中文字幕日韩二区| 91国产中文字幕| 18禁动态无遮挡网站| 极品人妻少妇av视频| 国产成人系列免费观看| www.精华液| 亚洲第一青青草原| 在线观看免费高清a一片| 老汉色av国产亚洲站长工具| av卡一久久| 欧美日韩综合久久久久久| 卡戴珊不雅视频在线播放| 午夜福利在线免费观看网站| 免费在线观看视频国产中文字幕亚洲 | 国产亚洲av片在线观看秒播厂| 自拍欧美九色日韩亚洲蝌蚪91| av视频免费观看在线观看| 久久人人爽av亚洲精品天堂| av有码第一页| 国产精品熟女久久久久浪| 久久天躁狠狠躁夜夜2o2o | 午夜影院在线不卡| 国产高清国产精品国产三级| 日韩,欧美,国产一区二区三区| 国产1区2区3区精品| 叶爱在线成人免费视频播放| av女优亚洲男人天堂| 免费高清在线观看视频在线观看| 女的被弄到高潮叫床怎么办| av片东京热男人的天堂| 日韩欧美一区视频在线观看| 国产精品女同一区二区软件| 亚洲人成77777在线视频| www日本在线高清视频| 1024视频免费在线观看| 国产欧美日韩一区二区三区在线| 亚洲免费av在线视频| 午夜福利乱码中文字幕| 久久青草综合色| 国产成人精品在线电影| 亚洲欧美中文字幕日韩二区| 亚洲熟女精品中文字幕| 亚洲国产精品国产精品| 多毛熟女@视频| 不卡视频在线观看欧美| 欧美少妇被猛烈插入视频| 久久亚洲国产成人精品v| 国产成人免费无遮挡视频| 黄频高清免费视频| 韩国av在线不卡| 在线观看免费午夜福利视频| 久久久国产欧美日韩av| www.熟女人妻精品国产| 久久久久精品人妻al黑| 精品久久久精品久久久| 母亲3免费完整高清在线观看| 免费看av在线观看网站| 一级毛片 在线播放| 两个人免费观看高清视频| 久久99热这里只频精品6学生| 国产av精品麻豆| 国产精品一区二区精品视频观看| 美女主播在线视频| 波多野结衣av一区二区av| 国产精品三级大全| 精品人妻熟女毛片av久久网站| 大片免费播放器 马上看| 成人毛片60女人毛片免费| 亚洲精品成人av观看孕妇| 又大又黄又爽视频免费| 欧美在线黄色| 日本av手机在线免费观看| 精品国产露脸久久av麻豆| 9热在线视频观看99| 一区在线观看完整版| 日日啪夜夜爽| 欧美精品av麻豆av| 大香蕉久久网| 欧美人与性动交α欧美软件| 亚洲精品,欧美精品| 熟妇人妻不卡中文字幕| 欧美中文综合在线视频| 黄色一级大片看看| av不卡在线播放| 纵有疾风起免费观看全集完整版| 黄色 视频免费看| 极品人妻少妇av视频| 午夜福利一区二区在线看| xxxhd国产人妻xxx| 色网站视频免费| 欧美精品亚洲一区二区| 成人国产麻豆网| 午夜福利乱码中文字幕| 一二三四中文在线观看免费高清| 欧美精品一区二区大全| 久久久久久久国产电影| 午夜激情久久久久久久| 欧美日本中文国产一区发布| 在线 av 中文字幕| 免费观看a级毛片全部| 午夜影院在线不卡| 亚洲情色 制服丝袜| 免费高清在线观看日韩| 久久人人爽人人片av| 菩萨蛮人人尽说江南好唐韦庄| 国产乱人偷精品视频| 一级片'在线观看视频| 日韩,欧美,国产一区二区三区| 欧美国产精品va在线观看不卡| 亚洲精品视频女| 女人被躁到高潮嗷嗷叫费观| 亚洲精品国产一区二区精华液| 高清视频免费观看一区二区| 美国免费a级毛片| 在线观看免费日韩欧美大片| 免费黄频网站在线观看国产| 91老司机精品| 亚洲av日韩在线播放| 悠悠久久av| 久久久久国产一级毛片高清牌| 国产伦理片在线播放av一区| 亚洲国产最新在线播放| 又大又爽又粗| 老熟女久久久| av一本久久久久| 欧美日韩综合久久久久久| 亚洲熟女毛片儿| 日本vs欧美在线观看视频| 亚洲欧美成人精品一区二区| 亚洲欧美清纯卡通| 大香蕉久久成人网| 亚洲成人免费av在线播放| 搡老乐熟女国产| 国产高清不卡午夜福利| 99热全是精品| 一本—道久久a久久精品蜜桃钙片| 蜜桃在线观看..| 99九九在线精品视频| 只有这里有精品99| 一级爰片在线观看| 久久97久久精品| 久久免费观看电影| 女人高潮潮喷娇喘18禁视频| 人人澡人人妻人| 国产精品av久久久久免费| 国产日韩欧美视频二区| 亚洲国产av新网站| 久久av网站| 天天躁日日躁夜夜躁夜夜| 日韩av在线免费看完整版不卡| 中文字幕亚洲精品专区| 丰满乱子伦码专区| 国产有黄有色有爽视频| 国产一级毛片在线| 久久久国产一区二区| 久久久久久久久久久免费av| 狂野欧美激情性bbbbbb| 国语对白做爰xxxⅹ性视频网站| 欧美日韩成人在线一区二区| 欧美亚洲日本最大视频资源| 亚洲免费av在线视频| 1024香蕉在线观看| 不卡av一区二区三区| 免费在线观看黄色视频的| 久久毛片免费看一区二区三区| 一本久久精品| 亚洲色图 男人天堂 中文字幕| 久久婷婷青草| 久久精品国产综合久久久| 国产精品无大码| 成年动漫av网址| 咕卡用的链子| 国产精品成人在线| 国产免费又黄又爽又色| 毛片一级片免费看久久久久| 亚洲精品,欧美精品| 日韩 欧美 亚洲 中文字幕| 国产1区2区3区精品| 欧美亚洲 丝袜 人妻 在线| 97在线人人人人妻| 人人澡人人妻人| 男女下面插进去视频免费观看| av线在线观看网站| 久久久久精品国产欧美久久久 | 亚洲成人免费av在线播放| 精品少妇久久久久久888优播| 一级片免费观看大全| 97在线人人人人妻| 久久青草综合色| 这个男人来自地球电影免费观看 | 欧美日韩视频高清一区二区三区二| 亚洲七黄色美女视频| 国产欧美日韩一区二区三区在线| 久久久久人妻精品一区果冻| 最近中文字幕2019免费版| 亚洲欧美激情在线|