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

    A Fast Charging–Cooling Coupled Scheduling Method for a Liquid Cooling-Based Thermal Management System for Lithium-Ion Batteries

    2022-01-26 06:39:22SiqiChenNengshengBoAkhilGrgXionginPengLingGo
    Engineering 2021年8期

    Siqi Chen,Nengsheng Bo,Akhil Grg,Xiongin Peng,Ling Go,*

    a State Key Laboratory of Digital Manufacturing Equipment and Technology,School of Mechanical Science and Engineering,Huazhong University of Science and Technology,Wuhan 430074,China

    b Intelligent Manufacturing Key Laboratory of Ministry of Education,Shantou University,Shantou 515063,China

    ABSTRACT Efficient fast-charging technology is necessary for the extension of the driving range of electric vehicles.However,lithium-ion cells generate immense heat at high-current charging rates.In order to address this problem,an efficient fast charging–cooling scheduling method is urgently needed.In this study,a liquid cooling-based thermal management system equipped with mini-channels was designed for the fastcharging process of a lithium-ion battery module.A neural network-based regression model was proposed based on 81 sets of experimental data,which consisted of three sub-models and considered three outputs:maximum temperature,temperature standard deviation,and energy consumption.Each sub-model had a desirable testing accuracy(99.353%,97.332%,and 98.381%)after training.The regression model was employed to predict all three outputs among a full dataset,which combined different charging current rates(0.5C,1C,1.5C,2C,and 2.5C(1C=5 A))at three different charging stages,and a range of coolant rates(0.0006,0.0012,and 0.0018 kg·s-1).An optimal charging–cooling schedule was selected from the predicted dataset and was validated by the experiments.The results indicated that the battery module’s state of charge value increased by 0.5 after 15 min,with an energy consumption lower than 0.02 J.The maximum temperature and temperature standard deviation could be controlled within 33.35 and 0.8 °C,respectively.The approach described herein can be used by the electric vehicles industry in real fast-charging conditions.Moreover,optimal fast charging–cooling schedule can be predicted based on the experimental data obtained,that in turn,can significantly improve the efficiency of the charging process design as well as control energy consumption during cooling.

    Keywords:Lithium-ion battery module Fast-charging Neural network regression Scheduling State of charge Energy consumption

    1.Introduction

    Owing to the significant challenges of fossil fuel shortages and greenhouse gas emissions,the global transportation industry is replacing traditional vehicles with clean energy automobiles[1].Rechargeable batteries are the central technology of electric vehicles(EVs),and lithium-ion battery cells are mainly used in energy storage devices for their many advantages:①high energy density,②high power density,③long cycle life,and ④low self-discharge rate.Lithium-ion battery cells are also widely used in various kinds of energy storage devices[2].However,despite its many advantages,a lithium-ion cell’s performance is sensitive to the ambient temperature.Previous research has shown that temperature strongly affects a lithium-ion battery cell’s charging/discharging efficiency,and this hinders the development of fast-charging technology.Severe temperatures accelerate the aging process,shorten the battery cells’lifetime,and can even result in thermal runaway.It is important to control the maximum temperature during the fast-charging process[2–4].Moreover,the uniformity of a battery module’s temperature distribution affects inconsistencies among the battery cells,which leads to the uneven rates of aging and a decay of the battery module’s lifespan.Considering the severe heat generated during the fast-charging process,it is vital to maintain the maximum temperature of a lithium-ion battery module within the range of 25–40 °C.Also,its temperature standard deviation(TSD)should be controlled at below 5 °C[5,6].

    A battery thermal management system(BTMS)consists of a cooling function and a heating function.This study focuses on the cooling technology for the fast-charging process.Several kinds of cooling strategies have been developed for battery thermal management,such as air cooling,liquid cooling,a cooling method based on the phase change principle,and combinations of these.Air cooling,with its advantages of lower consumption and a simple structure,has become a commonly used method in EVs.Various studies have investigated this method before,through theoretical analysis and experimental validation,which were performed to enhance the cooling performance.These studies include flow channel structural design[7,8],battery cell layout design[9–11],and certain controlling methods of the BTMS[12].However,air cooling cannot provide sufficient cooling performance in several applications,especially under extreme high ambient temperatures or high charging/discharging current rates with insufficient heat conductivity of the coolant(air)[13].

    Phase change materials(PCMs)are currently employed in cooling systems because they improve the uniformity of the battery module’s temperature distribution[14].However,PCMs cannot be widely employed in EVs because of their inadequate longterm cooling performance,which results from their intrinsic lower heat conductivity[15].Some recent studies were proposed to investigate PCM cooling efficiency under fast-charging.To enhance the heat conductivity,composite PCMs were studied by combining paraffin wax with aluminum,copper,graphite,graphene,or hexagonal boron nitride.However,the latent heat capacity of the PCMs decreased because of these combinations.Moreover,there is a significant volume consumption when PCM-based thermal management systems are used under fast-charging conditions.Under an 8C(1C=5 A)charging current rate,the volume consumption increases to 1.50–1.62 times that of traditional utilization[4,16,17].

    Heat pipes have also been used in BTMS for their gas–liquid phase-changing characteristics[18].Rao et al.[19]studied the cooling performance of heat pipes in experiments.Their results indicate that the maximum temperature can be controlled within 50°C if a single cell’s heat-generating power is less than 50 W.Wu et al.[20]conducted a comparison between natural convection,forced convection,and heat pipe-cooling,and concluded that heat pipe-cooling is efficient in maintaining temperature rise.However,the effective thermal conduction area of traditional tubular heat pipes has the feature of ‘‘line contact” between the heat pipes’ evaporators and battery cells’ surface,which results in insufficient cooling performance and considerable TSD of the battery module[21].Zhao et al.[22]designed a flat plate-based heat pipe with mini-channels inside,which contacted well with prismatic battery cells,and it has become a candidate choice for cooling structures.

    The advantages and disadvantages of various battery thermal management methods are compared in Table 1.Compared with the other cooling strategies,liquid cooling is a more efficient approach owing to its higher specific heat capacity and a uniform temperature distribution.However,liquid cooling-based thermal management systems can be error prone,and the assembly process and maintenance are much more complicated and expensive because the BTMS includes some additional equipment(pipes,pumps,or heat exchangers)[23].Moreover,the pressure drop between the inlet and outlet leads to unavoidable energy consumption[23].Furthermore,the battery cells of a battery module should be arranged under adiabatic conditions to avoid thermal runaway.Mei et al.[24]proposed three three-dimensional layered electrochemical–thermal models with various unit numbers.These were validated by experiments at various discharging current rates(0.5C,1C,1.5C,and 2C).The results indicated that a one-unit model needed a lower computational load,and this was helpful for the adiabatic studies.Previous studies on liquid cooling were conducted mainly on the structural design of flow channels[25–27]and on the enhancement of heat conduction by the coolingmedium[28–30].Chen et al.[31]proposed a comprehensive optimization design for a liquid cooling-based battery module and validated the design by experiments.The results indicated that the battery module’s thermal performance(maximum temperature and TSD)was significantly enhanced after the multi-objective structural optimization design.A liquid cooling-based BTMS was also verified to be effective in the discharging process of a lithium-ion battery module under a wide range of current rates.Panchal et al.[32]designed a liquid cooling plate for a lithium-ion battery module discharged at a 4C current rate.A cascade-based liquid cooling structure has been shown to have the ability to reduce a lithiumion battery module’s TSD from 7 to 2°C[33].

    Table 1 Advantages and disadvantages of various battery thermal management methods.

    Previous studies have revealed the superiority of liquid coolingbased battery modules.It has been shown that a battery module’s maximum temperature and TSD can be maintained within an appropriate range.However,these studies focused on the structural design of the thermal management system and on the thermal performance of liquid cooling-based battery modules under various working conditions.There is a lack of studies investigating a fast charging–cooling coupled scheduling which is necessary in order to enhance the driving range and charging efficiency.Moreover,in order to ensure driving safety and avoid overheating or an uneven temperature distribution,thermal management during the fast-charging process is also significant.There is a limited number of studies on combining fast-charging process scheduling with thermal management.The schedule selection reaches a reasonable tradeoff between the charging speed,cooling efficiency,and energy consumption that is applicable for EVs.

    This study addressed the above issues through

    (1)Designing a liquid cooling-based thermal management system equipped with mini-channels for the fast-charging process of a lithium-ion battery module;

    (2)Proposing a neural network-based regression model in order to design an efficient fast charging–cooling coupled scheduling method;

    (3)Employing the proposed regression model to predict three results from a full dataset,which combines a range of charging current rates(0.5C,1C,1.5C,2C,and 2.5C)at three different charging stages,and a range of coolant rates(0.0006,0.0012,and 0.0018 kg·s-1);

    (4)Selecting an optimal fast charging–cooling coupled schedule from the predicted dataset by setting constraints;

    (5)Validating the optimal schedule through the experiments.

    2.Design of liquid cooling-based battery module

    2.1.Cooling structure design for fast-charging

    A liquid cooling-based battery module is shown in Fig.1.A kind of 5 A·h lithium-ion cell was selected,with its working voltage ranging from 3.20 to 3.65 V.Eight prismatic battery cells(considering the battery testing machine’s voltage limitation)were arranged on the top surface of a cooling plate.The heat generated from the battery cells was conducted and dissipated to the coolant flowing through the mini-channels,which resulted in a lower volume cost and provided sufficient thermal conductivity.The side surfaces between the single cells were separated by heatinsulating cotton designed to avoid the interaction between single cells under the effect of overheating.This structural design is accessible for EVs in real applications.

    Aluminum alloy 7075 was selected for manufacturing of the cooling plate by milling,considering the high thermal conductivity and low cost of this material.Antifreeze was selected as the coolant,which consisted of ethylene glycol(50%)and water(50%).Antifreeze is commonly used in EVs and traditional vehicles for its lower ice point,especially in high-latitude areas.

    2.2.Experimental configuration design

    The experimental setup is shown in Fig.2.Eight single cells were connected in series,and the positive and negative electrodes of the module were connected to a battery testing machine(BTS-4000,Neware Electronics Co.,Ltd.,China).Four K-type thermocouples were attached to the same sidewall of each cell to measure the real-time temperature during the 15 min fast-charging process.A pressure sensor was attached to the inlet of the cooling plate;the real-time inlet pressure was also measured.The battery module’s temperature and pressure were recorded by two data loggers(THM001,Penghe Electronics Co.,Ltd.,China),with the values shown on the screens.A motor peristaltic tubing pump(Kamoer,China)was used to supply the coolant flow for the cooling system.The ambient temperature in the laboratory was set at 25 °C.The initial state of charge(SOC)value of the battery module was discharged to 0.2 before each set of experiments.

    Fig.1.Schematic of the liquid cooling-based lithium-ion battery module.(a)Diagram of lithium-ion battery module;(b)diagram of mini-channel-based cooling plate.

    Fig.2.Image of the experimental setup.

    3.Artificial neural network regression model

    3.1.Artificial neural network regression

    An artificial neural network is a kind of machine learning model employed for data classification or data prediction.The model structure is constructed based on data and learning rules[34].A neural network regression model is trained with data,based on a training algorithm,to predict a subsequent set of data.

    As shown in Fig.3,an artificial neural network model consists of some nodes/neurons,which are set in multiple layers:an input layer,one or more hidden layers,and an output layer[34–36].Each node/neuron has an activation function,which calculates how much neuron is ‘‘stimulated.” At each layer,the collections of nodes/neurons transform the input parameters;these parameters are distributed to the next layer,which is described by Eqs.(1)–(3):

    Fig.3.Diagram of neural network. x represents the input to the first layer; z represents the first layer’s output; y represents the output; i and j represent the neural network node index; xi represents the ith input to the input layer; zj represents the jth output of the first layer; wji(l) represents the weight between the jth node in the first layer and the ith node in the input layer.

    wherexrepresents the input to the first layer;zrepresents the first layer’s output;iandjrepresent the neural network node index;nrepresents the index of the layer;represents the weight between thejth node in the first layer and theith node in the input layer;arepresents the input of the output layer;represents the output value ofith node in(n+ 1)th layer after being activated by the activation function;andwandw0represent the weight and bias between the neurons,respectively,which measures the significance of the data passed along the link(synapse).F(a)represents the activation function,which employs the hidden layer’s aggregated output to calculate the outputy[35,37].

    The initial weights and biases are assigned randomly,and the training process continues until the desired output is obtained,which is evaluated by the cost function Eq.(4):

    whereorepresents the desired output;E(w)represents the cost function evaluating the training process;wrepresents the weight;andkrepresents the index of the cost function calculation.The Levenberg–Marquardt(LM)algorithm is used in the neural network training process,which is a variation of gradient descent.The weight and bias of the neural network model are changed during the training process to minimize the error,which is described by Eq.(5):

    where J=?E/?wrepresents the full-scale Jacobian matrix related tow,I represents the identity matrix,μ is combination coefficient,anderepresents the prediction error.

    The LM algorithm starts with a forward computation by Eqs.(1)–(3).The prediction errors of the output layer and hidden layer are calculated by Eqs.(6)–(8):

    where δ represents the prediction error of the output layer.

    As shown in Eqs.(9)and(10),the Jacobian is calculated by a backpropagation process:

    3.2.Configuration of neural network model

    As shown in Fig.4,the entire charging process is divided into three stages;each stage lasts 5 min.I1,I2,andI3represent the charging current value at stages 1,2,and 3,respectively.Qrepresents the coolant flow rate of the fast-charging process.I1,I2,I3,andQare set as input parameters of the neural network models.

    This study had three objectives:①maximum temperature(Tmax),② TSD,and ③ energy consumption(W)of the fastcharging process.These three parameters were set as output parameters of three neural network models(NN1,NN2,and NN3).

    The value of the liquid cooling system’s energy consumption is calculated by Eq.(11):

    Fig.4.Schematic of the three stages during the fast-charging process.

    wherePrepresents the power consumption of the cooling system;trepresents time.

    The TSD is calculated based on all of the measured temperature values from the temperature sensors and is described by Eq.(12):

    whereTrrepresents the temperature value of each temperature sensor,represents the average temperature value of all the temperature sensors,andrrepresents the battery cell index.

    All of the training data are from experiments and three charging current rates were selected for the three stages of the experimental charging process:①0.5C,②1.5C,and ③2.5C.Three coolant flow rates were also selected for the charging process:①36 mL·min-1,②72 mL·min-1,and ③108 mL·min-1.As shown in Fig.5,a total of 81 sets of experimental data were obtained for the training process.

    The theoretical capacity is defined based on the assumption that all of the active substances in an electrode participate in the battery cells’ electrochemical reaction,which is treated as the highest theoretical value calculated according to Faraday’s law.In this study,the theoretical capacity of the battery cell is 5 A·h.The SOC value of the battery module is employed to evaluate the charging effect,which is calculated by Eq.(13):

    whereIrepresents the charging current.

    Specific information about the training dataset and corresponding SOC value is listed in Table 2.

    Table 2 (continued)

    Table 2 Training dataset of the neural network model.

    Fig.5.Composition of the experiment-based training dataset.

    3.3.Structure of the neural network regression model

    Unlike some battery cells’ behaviour-based governing equations,a neural network regression model is a ‘‘black box” model,which does not provide insight into how the hidden model’s structure actually works.Moreover,there are two problems for a nonlinear neural network regression model:① underfitting or overfitting,and ②the possible stagnation of the training process at the local optimal stage.

    Some solutions are proposed to obtain a better regression performance:①input parameter preprocessing,so that representative ‘‘features” can be extracted from the data;②increasing the overall percentage of training datasets to cover the anticipated operational range;and ③ comprehensive cross-validation to enhance training performance when encountering an unknown usage profile[34].

    The configurations of three neural network regression models are shown in Table 3.Each model consists of three layers,and these models are designed with the same activation function(a feedforward LM backpropagation approach)in the training process.As Fig.6 shows,the traditional neural network model is split into three sub-models,one for each objective(output parameter).Each input layer consists of the four input parameters(I1,I2,I3,andQ)and is followed by a hidden layer.The output parameter estimation responds to each of the input layers.The number of nodes in the input layer shows the exact number of input parameters used in the model.To obtain a continuous result,the hidden layers feed into another aggregation layer,which combines the results from the previous layers and feeds them to the output layer.To avoid overfitting and underfitting,two nodes/neurons in the hidden layer were implemented to get the desirable regression performance.Among the 81 sets of experimental data,a number of 73(90%)sets of experimental data were employed for training;four(5%)sets of experimental data were employed for validation;four(5%)sets of experimental data were employed for testing.One thousand epochs were set for the neural network regression training process.

    Table 3 Configurations of three neural network regression models.

    3.4.Regression performance evaluation of neural network model

    As mentioned in Section 3.3,LM backpropagation was employed in order to enhance the model training rate related to pure error backpropagation or steepest descent,and this algorithm maintained the accuracy of the trained model.

    The neural network regression model was trained with the designed model structure,input parameters,and number of nodes.The accuracy of both the training and the prediction/estimation was evaluated by the mean absolute error(MSE),which is expressed as Eq.(14).

    Fig.6.Diagram of the proposed neural network model.

    whereeavgrepresents the average absolute error,Mrepresents the number of data points,Outputmrepresents themth estimated output parameter,and Outputr,mrepresents themth reference output parameter.The minimum error for the neural network estimation was found with the number of nodes per output parameter estimation layer set as two.

    The training process of these three neural network models is shown in Fig.7.With the above neural network configuration,the estimation performance was improved,and the three regression models achieved a desirable speed of error convergence after initialization,during the training process.

    As shown in Fig.8,high fitting accuracy rates were obtained after the training process.For the training process,the fitting accuracy rates of NN1,NN2,and NN3were 95.123%,82.866%,and 98.142%,respectively.For the testing process,the fitting accuracy rates of NN1,NN2,and NN3were 99.353%,97.332%,and 98.381%,respectively.The overall fitting accuracy rates of NN1,NN2,and NN3were 95.478%,83.496%,and 98.049%,respectively.All of the three neural network models achieved high accuracy rates,especially in testing;and the predicted values of the output parameters were close to the reference values.There was a slight difference between the measured and predicted energy consumption values,but the difference was within the tolerance limits.The results indicate that the relationship between the input parameters and output parameters was accurately modelled.With good regression performance of these three models,further prediction/estimation should be reliable.

    Fig.7.Training process of three neural network models.(a)Training process of NN1;(b)training process of NN2;(c)training process of NN3.Mu:reduction factor;Gamk:display frequency;ssX:inflation factor;Val fail:validation failure time.

    Fig.8.Regression performance of the three neural network models.(a)Regression performance of NN1;(b)regression performance of NN2;(c)regression performance of NN3.Y = T:trend line in ideal condition with 100% fitting rate;Y:output;T:target.

    The error histogram plots of these three neural network models are shown in Fig.9.For the three neural network models,a number of 20 bins were distributed following the Gaussian distribution rules,which proved the representativeness of the training data.

    4.Optimal charging–cooling schedule selection and experimental validation

    4.1.Estimation of fast charging–cooling schedules based on the trained regression model

    The trained neural network regression model was employed to estimate the results of some given input combinations.Some estimation results were compared with reference to some known experimental data to determine the accuracy and robustness of the trained model.

    According to the current rates of fast-charging technology features,a total of five charging current rates were selected for the three stages of the charging process:①0.5C,②1C,③1.5C,④2C,and ⑤2.5C.Three coolant flow rates were also selected for the charging process:①36 mL·min-1,②72 mL·min-1,and③ 108 mL·min-1.Fig.10 shows the composition of the fast-charging designs.

    There were 375 sets of fast charging–cooling schedules in a full dataset;and the maximum temperature,TSD,and energy consumption of these charging–cooling schedules were predicted using the proposed three neural network regression models(NN1,NN2,and NN3).Specific information and the SOC values of these fast charging–cooling schedules are listed as in Appendix A Table S1.

    Fig.9.Error histogram of three neural network models.(a)Error histogram of NN1;(b)error histogram of NN2;(c)error histogram of NN3.Errors = Targets–Outputs.

    4.2.Optimal charging–cooling schedule selection

    To obtain an efficient fast-charging schedule,the thermal performance of the battery module should be controlled(maximum temperature and temperature distribution uniformity),the SOC value of the battery module must be charged to a high level within a short time,and the energy consumption of the cooling system should also be controlled.

    In this study,the specific requirements for the optimal fast charging–cooling schedule are listed as follows:

    (1)The SOC value of the battery module increases by 0.5 after the fast-charging process(15 min);

    Fig.10.Composition of the fast-charging designs.

    (2)The maximum temperature of the battery module should be lower than 306.5 K(33.35 °C);

    (3)The TSD should be controlled within 0.8 °C;

    (4)The energy consumption of the cooling system during the fast-charging process should be less than 0.02 J.

    Addressing all of these constraints is essential for a comprehensive and reliable charging–cooling schedule design.The cooling effect,uniform temperature distribution,charging speed,and energy consumption should all be considered as having the same degree of importance.These selection standards are described by Eq.(15):

    Schedule 50 was selected as the optimal fast charging–cooling schedule,and it was validated in the experiments using the above-described experimental setup.

    4.3.Experimental validation

    Figs.11 and 12 show the temperature and pressure curves of the fast charging–cooling design Schedule 50 in the experiments.

    A comparison between the selected optimal fast charging–cooling schedule and the experimental data is presented in Table 4.

    Table 4 Comparison between the selected optimal fast charging–cooling schedule and the experimental data.

    As shown in the results,the differences between the predicted and experimental data are within the limits of tolerance.The deviation of maximum temperature,TSD,and the energy consumption of the cooling system are 0.468 °C,0.045327 °C,and 0.001694 J,respectively.This proves the prediction accuracy of the three proposed neural network regression models.Consequently,these models can be efficiently employed for fast charging–cooling coupled scheduling in real applications,with much lower costs.

    Fig.11.Temperature curve of the fast charging–cooling design Schedule 50.

    Fig.12.Pressure curve of the fast charging–cooling design Schedule 50.

    5.Conclusions and further work

    This study proposed a neural network-based regression model for fast charging–cooling coupled scheduling,which significantly saves time and cost during the fast charging-and-cooling design process.The proposed neural network model configurations achieve the desirable regression performance,and the designed configuration makes the estimation more realistic.The main conclusions of this study are listed below.

    (1)The proposed model was trained based on 81 sets of experimental data,which consisted of three sub-models considering three targets:maximum temperature,TSD,and energy consumption of the cooling system.The regression results achieved high regression performance in tests,with 99.353%,97.332%,and 98.381% accuracy,respectively.

    (2)The regression model was utilized to predict three target values for all of the combinations among a wide range of charging current rates(0.5C,1C,1.5C,2C,and 2.5C),at three different charging stages and a range of coolant rates(0.0006,0.0012,and 0.0018 kg·s-1).Moreover,an optimal fast charging–cooling coupled schedule was selected from the predicted dataset based on the constraints of the three targets.

    (3)The optimal schedule was validated by experiments,and the results show that the battery module’s SOC value increased by 0.5 after 15 min of the fast-charging process.The energy consumption of this process was controlled within 0.02 J,and the maximum temperature and TSD were lower than 33.35 and 0.8°C,respectively.

    The proposed method can be employed by the EVs industry in real fast-charging conditions,and the experimental test data can be used for regression training to predict the optimal fast charging–cooling schedule,which significantly promotes the efficiency of the charging process design and controls the energy consumption during the cooling process.

    Further work will aim to develop and train neural network models to cover wider applications,and the effects of temperature changes in the environment and battery degradation may also be included.

    Acknowledgements

    This work was supported by the Program for Huazhong University of Science and Technology(HUST)Academic Frontier Youth Team(2017QYTD04),the Program for HUST Graduate Innovation and Entrepreneurship Fund(2019YGSCXCY037).Authors acknowledge Grant DMETKF2018019 by State Key Laboratory of Digital Manufacturing Equipment and Technology,Huazhong University of Science and Technology.This study was also financially supported by the Guangdong Science and Technology Project(2016B020240001),and the Guangdong Natural Science Foundation(2018A030310150).

    Compliance with ethics guidelines

    Siqi Chen,Nengsheng Bao,Akhil Garg,Xiongbin Peng,and Liang Gao declare that they have no conflict of interest or financial conflicts to disclose.

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.ultsonch.2019.104640.

    免费观看人在逋| 国产aⅴ精品一区二区三区波| 国产片内射在线| 久久九九热精品免费| 欧美日韩中文字幕国产精品一区二区三区| 精品福利观看| 女同久久另类99精品国产91| 91九色精品人成在线观看| 午夜久久久在线观看| 中文字幕最新亚洲高清| 黄色丝袜av网址大全| 一级毛片女人18水好多| 日韩国内少妇激情av| 精品国产亚洲在线| 人人妻,人人澡人人爽秒播| 少妇被粗大的猛进出69影院| 99热只有精品国产| 在线观看一区二区三区| 十八禁网站免费在线| 99在线人妻在线中文字幕| 亚洲国产精品成人综合色| 欧美成狂野欧美在线观看| 丝袜人妻中文字幕| 久热这里只有精品99| 国产一区二区在线av高清观看| or卡值多少钱| 欧美日韩乱码在线| 中亚洲国语对白在线视频| 一本久久中文字幕| 91国产中文字幕| 欧美色视频一区免费| 久久精品成人免费网站| 一进一出抽搐动态| 国产区一区二久久| 精品无人区乱码1区二区| 首页视频小说图片口味搜索| 首页视频小说图片口味搜索| www国产在线视频色| 日韩视频一区二区在线观看| 久久中文字幕人妻熟女| 欧美性猛交╳xxx乱大交人| 国产高清有码在线观看视频 | 777久久人妻少妇嫩草av网站| 欧美激情久久久久久爽电影| 日日摸夜夜添夜夜添小说| 黄色a级毛片大全视频| 久久草成人影院| 午夜老司机福利片| 人成视频在线观看免费观看| 俺也久久电影网| 青草久久国产| 国产区一区二久久| 国产97色在线日韩免费| 久久久久精品国产欧美久久久| 国产伦一二天堂av在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲中文字幕一区二区三区有码在线看 | 夜夜爽天天搞| 久久香蕉激情| 精品久久久久久成人av| 丁香六月欧美| 欧美日韩精品网址| 黑人巨大精品欧美一区二区mp4| 女人高潮潮喷娇喘18禁视频| 99精品欧美一区二区三区四区| 亚洲久久久国产精品| 免费女性裸体啪啪无遮挡网站| 1024视频免费在线观看| 最近在线观看免费完整版| 女警被强在线播放| 亚洲色图 男人天堂 中文字幕| 免费人成视频x8x8入口观看| 久久天躁狠狠躁夜夜2o2o| 久久九九热精品免费| 老司机午夜福利在线观看视频| 国产人伦9x9x在线观看| 91在线观看av| 国产成人啪精品午夜网站| 欧美性猛交黑人性爽| 国产三级黄色录像| 亚洲aⅴ乱码一区二区在线播放 | 一夜夜www| 色婷婷久久久亚洲欧美| 欧美精品亚洲一区二区| 波多野结衣高清无吗| 1024香蕉在线观看| 波多野结衣高清作品| 俄罗斯特黄特色一大片| 51午夜福利影视在线观看| 中国美女看黄片| 婷婷丁香在线五月| 黑人操中国人逼视频| 每晚都被弄得嗷嗷叫到高潮| 久热这里只有精品99| 国产精品日韩av在线免费观看| 大型av网站在线播放| 嫁个100分男人电影在线观看| 国产真实乱freesex| 久久伊人香网站| 精品久久久久久久人妻蜜臀av| 国产熟女xx| 亚洲精品国产区一区二| 一区福利在线观看| 黄色视频,在线免费观看| 中文字幕人妻熟女乱码| 婷婷丁香在线五月| 99精品在免费线老司机午夜| 国产亚洲精品一区二区www| 欧美日韩黄片免| 在线永久观看黄色视频| 久久久精品欧美日韩精品| 一夜夜www| 免费在线观看视频国产中文字幕亚洲| 国产激情偷乱视频一区二区| 色婷婷久久久亚洲欧美| 亚洲精品在线观看二区| 波多野结衣高清作品| 国产精品 国内视频| 精品熟女少妇八av免费久了| 中文字幕最新亚洲高清| 欧美不卡视频在线免费观看 | 亚洲av成人不卡在线观看播放网| 欧美激情 高清一区二区三区| 亚洲性夜色夜夜综合| 亚洲欧美一区二区三区黑人| 99热6这里只有精品| 人人澡人人妻人| 一本综合久久免费| 精品国产超薄肉色丝袜足j| 国产又色又爽无遮挡免费看| 丝袜美腿诱惑在线| 一进一出抽搐动态| 亚洲aⅴ乱码一区二区在线播放 | 深夜精品福利| 99热只有精品国产| 一个人免费在线观看的高清视频| 亚洲欧美一区二区三区黑人| 国产av又大| 午夜成年电影在线免费观看| 欧美一级a爱片免费观看看 | 国产1区2区3区精品| 婷婷亚洲欧美| 免费看a级黄色片| 久久久国产精品麻豆| 母亲3免费完整高清在线观看| 又大又爽又粗| 亚洲一码二码三码区别大吗| 一区二区三区精品91| 久久香蕉激情| 99久久精品国产亚洲精品| 黄色毛片三级朝国网站| 国产高清视频在线播放一区| 国产精品久久久久久精品电影 | 成人精品一区二区免费| 国产午夜福利久久久久久| 午夜福利一区二区在线看| 一个人免费在线观看的高清视频| 无人区码免费观看不卡| 校园春色视频在线观看| 精品高清国产在线一区| 国语自产精品视频在线第100页| 麻豆久久精品国产亚洲av| 人人妻,人人澡人人爽秒播| 亚洲成av片中文字幕在线观看| 黄色视频,在线免费观看| 大香蕉久久成人网| av片东京热男人的天堂| 欧美大码av| 国产成人精品久久二区二区免费| 国产精品1区2区在线观看.| 波多野结衣av一区二区av| 国产一卡二卡三卡精品| 日本a在线网址| 国产麻豆成人av免费视频| 国产av一区在线观看免费| 欧美精品啪啪一区二区三区| 精品久久久久久,| 日韩欧美三级三区| 淫秽高清视频在线观看| 99国产精品一区二区蜜桃av| 久久香蕉激情| 最新美女视频免费是黄的| www日本在线高清视频| 又紧又爽又黄一区二区| 日本 av在线| 一二三四社区在线视频社区8| 一个人观看的视频www高清免费观看 | 一本久久中文字幕| 国产国语露脸激情在线看| 人人澡人人妻人| 曰老女人黄片| 91成人精品电影| 亚洲人成77777在线视频| 免费在线观看成人毛片| 黄片播放在线免费| 国语自产精品视频在线第100页| 在线视频色国产色| 亚洲精品美女久久av网站| 国产国语露脸激情在线看| 午夜日韩欧美国产| 亚洲电影在线观看av| 亚洲中文字幕一区二区三区有码在线看 | 国产精品久久视频播放| 亚洲欧美精品综合久久99| 九色国产91popny在线| 日韩精品中文字幕看吧| 97碰自拍视频| 黄片小视频在线播放| 1024视频免费在线观看| 在线av久久热| 一区二区三区国产精品乱码| 国语自产精品视频在线第100页| 亚洲欧美日韩无卡精品| 日韩欧美免费精品| 成年人黄色毛片网站| 国产精品久久电影中文字幕| 国产亚洲精品综合一区在线观看 | 美女午夜性视频免费| 麻豆av在线久日| 欧美中文综合在线视频| 曰老女人黄片| 丁香欧美五月| 女生性感内裤真人,穿戴方法视频| 满18在线观看网站| 午夜老司机福利片| 精品久久久久久久久久免费视频| 亚洲成人久久性| 久久久精品欧美日韩精品| 国产区一区二久久| 国产精品99久久99久久久不卡| 免费看a级黄色片| 久久久水蜜桃国产精品网| 欧美色视频一区免费| 黄色毛片三级朝国网站| av天堂在线播放| 亚洲熟妇中文字幕五十中出| 啪啪无遮挡十八禁网站| 亚洲avbb在线观看| 欧美国产精品va在线观看不卡| 两个人免费观看高清视频| 中文字幕另类日韩欧美亚洲嫩草| 99在线视频只有这里精品首页| 久久国产精品人妻蜜桃| 美女国产高潮福利片在线看| 男人舔女人下体高潮全视频| 免费人成视频x8x8入口观看| 亚洲avbb在线观看| 久热爱精品视频在线9| 一a级毛片在线观看| 中国美女看黄片| 高清在线国产一区| 曰老女人黄片| 啦啦啦免费观看视频1| 制服丝袜大香蕉在线| 国产精品亚洲av一区麻豆| 人人澡人人妻人| 777久久人妻少妇嫩草av网站| 国产在线观看jvid| 岛国视频午夜一区免费看| 日本一区二区免费在线视频| 久久久久久久久中文| 老司机靠b影院| 日本 av在线| 国产熟女午夜一区二区三区| 午夜日韩欧美国产| 免费看a级黄色片| 俄罗斯特黄特色一大片| 午夜福利成人在线免费观看| 动漫黄色视频在线观看| 午夜福利免费观看在线| 国内精品久久久久久久电影| 一本综合久久免费| 嫩草影视91久久| 亚洲av第一区精品v没综合| 可以在线观看的亚洲视频| 91在线观看av| 亚洲真实伦在线观看| or卡值多少钱| 欧美丝袜亚洲另类 | 久热这里只有精品99| 99riav亚洲国产免费| 91老司机精品| 亚洲三区欧美一区| 波多野结衣高清作品| 成年人黄色毛片网站| 波多野结衣巨乳人妻| 丁香欧美五月| 免费搜索国产男女视频| 高潮久久久久久久久久久不卡| 一进一出抽搐gif免费好疼| 亚洲熟妇熟女久久| 女人被狂操c到高潮| 亚洲五月色婷婷综合| 日韩 欧美 亚洲 中文字幕| 两个人视频免费观看高清| 午夜久久久在线观看| 变态另类成人亚洲欧美熟女| 久久狼人影院| 视频在线观看一区二区三区| 中文字幕最新亚洲高清| 老熟妇仑乱视频hdxx| 免费人成视频x8x8入口观看| 成人18禁在线播放| 欧美三级亚洲精品| 亚洲av电影在线进入| 精品久久蜜臀av无| 国产精品1区2区在线观看.| 日韩欧美一区视频在线观看| 国产精品久久久久久人妻精品电影| 99热这里只有精品一区 | 可以免费在线观看a视频的电影网站| 岛国在线观看网站| 久久婷婷人人爽人人干人人爱| 欧美成人性av电影在线观看| 色尼玛亚洲综合影院| 成年版毛片免费区| 草草在线视频免费看| 国产色视频综合| 色播在线永久视频| 一级a爱片免费观看的视频| 给我免费播放毛片高清在线观看| av电影中文网址| 欧美日本亚洲视频在线播放| 精品国产乱码久久久久久男人| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲精品久久久久久毛片| 日韩欧美一区二区三区在线观看| 久久久久九九精品影院| 国产成人精品久久二区二区免费| 中文字幕人成人乱码亚洲影| 一区二区三区精品91| 国产激情偷乱视频一区二区| www.自偷自拍.com| 亚洲七黄色美女视频| 成人18禁高潮啪啪吃奶动态图| 久久婷婷成人综合色麻豆| 精品久久久久久久毛片微露脸| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 久久99热这里只有精品18| 男人的好看免费观看在线视频 | 又黄又爽又免费观看的视频| 亚洲久久久国产精品| 国产区一区二久久| 精品福利观看| 亚洲第一青青草原| 99久久精品国产亚洲精品| 午夜两性在线视频| 免费高清在线观看日韩| √禁漫天堂资源中文www| 老汉色∧v一级毛片| 国产一区二区三区在线臀色熟女| 非洲黑人性xxxx精品又粗又长| 99久久99久久久精品蜜桃| 亚洲第一av免费看| 可以在线观看的亚洲视频| 夜夜看夜夜爽夜夜摸| 久久久久久久精品吃奶| 成人精品一区二区免费| 久久天堂一区二区三区四区| bbb黄色大片| 女人爽到高潮嗷嗷叫在线视频| 国产一区在线观看成人免费| 这个男人来自地球电影免费观看| 一区二区三区国产精品乱码| 成人三级黄色视频| 久久性视频一级片| 国产精品电影一区二区三区| avwww免费| 无人区码免费观看不卡| 亚洲欧美精品综合一区二区三区| 国产精品爽爽va在线观看网站 | 亚洲专区国产一区二区| 99在线人妻在线中文字幕| 老司机福利观看| 亚洲精品一区av在线观看| 午夜精品久久久久久毛片777| 精品高清国产在线一区| 日韩欧美一区视频在线观看| 老司机午夜十八禁免费视频| 婷婷精品国产亚洲av| 亚洲国产欧洲综合997久久, | 搡老岳熟女国产| 久久久国产成人免费| 老熟妇乱子伦视频在线观看| 国产精品九九99| 校园春色视频在线观看| 日韩免费av在线播放| 亚洲精华国产精华精| 亚洲午夜理论影院| www.999成人在线观看| 免费看日本二区| 身体一侧抽搐| 18禁美女被吸乳视频| 一级毛片女人18水好多| 午夜免费观看网址| 最新在线观看一区二区三区| 亚洲中文字幕日韩| 欧美性猛交╳xxx乱大交人| 中文字幕最新亚洲高清| 巨乳人妻的诱惑在线观看| 亚洲人成伊人成综合网2020| 国产精品1区2区在线观看.| 亚洲中文av在线| www.熟女人妻精品国产| 在线视频色国产色| 可以在线观看的亚洲视频| 黄网站色视频无遮挡免费观看| 欧美日韩亚洲综合一区二区三区_| 欧美大码av| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品 国内视频| 一本大道久久a久久精品| 亚洲久久久国产精品| 校园春色视频在线观看| 亚洲 欧美 日韩 在线 免费| 精品免费久久久久久久清纯| 一区二区三区国产精品乱码| 午夜免费观看网址| 99热6这里只有精品| 久久国产精品影院| 国产亚洲精品久久久久5区| 啦啦啦观看免费观看视频高清| 国产不卡一卡二| 一进一出抽搐动态| 女同久久另类99精品国产91| 亚洲狠狠婷婷综合久久图片| 50天的宝宝边吃奶边哭怎么回事| www日本黄色视频网| 亚洲人成伊人成综合网2020| 悠悠久久av| 午夜福利视频1000在线观看| 亚洲三区欧美一区| 男人的好看免费观看在线视频 | 国产视频一区二区在线看| 久久久水蜜桃国产精品网| 老司机福利观看| 日本 av在线| 九色国产91popny在线| 日韩三级视频一区二区三区| 在线国产一区二区在线| 欧美成人性av电影在线观看| 黄网站色视频无遮挡免费观看| 别揉我奶头~嗯~啊~动态视频| 亚洲精品在线观看二区| 精品久久久久久成人av| 亚洲成av片中文字幕在线观看| 一级毛片高清免费大全| 999久久久国产精品视频| 亚洲色图 男人天堂 中文字幕| 麻豆久久精品国产亚洲av| 亚洲成人久久爱视频| 亚洲成人精品中文字幕电影| 特大巨黑吊av在线直播 | 啦啦啦韩国在线观看视频| 国产三级黄色录像| 美国免费a级毛片| 国产成人精品久久二区二区91| 麻豆成人av在线观看| 给我免费播放毛片高清在线观看| 欧美又色又爽又黄视频| 性欧美人与动物交配| 久久99热这里只有精品18| 看免费av毛片| 国产伦人伦偷精品视频| 亚洲成人免费电影在线观看| 99热6这里只有精品| 亚洲国产看品久久| 九色国产91popny在线| 久久久久久久午夜电影| 日本熟妇午夜| 午夜福利视频1000在线观看| 午夜福利在线在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av美国av| 午夜日韩欧美国产| a级毛片在线看网站| 免费在线观看影片大全网站| www国产在线视频色| 1024香蕉在线观看| 精品卡一卡二卡四卡免费| 欧美成人性av电影在线观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲久久久国产精品| 9191精品国产免费久久| 色精品久久人妻99蜜桃| 亚洲全国av大片| 特大巨黑吊av在线直播 | 香蕉久久夜色| 日本 欧美在线| 亚洲中文字幕日韩| xxxwww97欧美| 午夜免费观看网址| 91在线观看av| 日韩高清综合在线| 自线自在国产av| 精品久久久久久成人av| 国产黄a三级三级三级人| 亚洲va日本ⅴa欧美va伊人久久| 国产精品 国内视频| 午夜视频精品福利| 久久精品国产99精品国产亚洲性色| 亚洲精品久久成人aⅴ小说| 国产av一区在线观看免费| 女性生殖器流出的白浆| 国产黄片美女视频| bbb黄色大片| 男女下面进入的视频免费午夜 | 欧美午夜高清在线| 国产亚洲精品综合一区在线观看 | 国产人伦9x9x在线观看| 免费高清在线观看日韩| 侵犯人妻中文字幕一二三四区| 久9热在线精品视频| 午夜福利18| 精品第一国产精品| 听说在线观看完整版免费高清| 亚洲专区国产一区二区| x7x7x7水蜜桃| 亚洲激情在线av| 亚洲中文日韩欧美视频| 夜夜爽天天搞| 国产爱豆传媒在线观看 | 久热爱精品视频在线9| 人人妻人人澡人人看| 久久精品aⅴ一区二区三区四区| 可以在线观看毛片的网站| 日本免费a在线| 久久久久久免费高清国产稀缺| 在线av久久热| 中文字幕人妻丝袜一区二区| 亚洲va日本ⅴa欧美va伊人久久| 成年人黄色毛片网站| 亚洲国产精品sss在线观看| 淫妇啪啪啪对白视频| 黄频高清免费视频| 老司机靠b影院| 精品日产1卡2卡| 19禁男女啪啪无遮挡网站| 老熟妇乱子伦视频在线观看| 国产国语露脸激情在线看| 久热这里只有精品99| 久久久久久久久中文| 麻豆av在线久日| 国产av在哪里看| 九色国产91popny在线| 欧美日韩乱码在线| 悠悠久久av| 熟女少妇亚洲综合色aaa.| 国产亚洲欧美在线一区二区| 久久久久国产精品人妻aⅴ院| 黄片大片在线免费观看| 亚洲国产精品999在线| 一级毛片女人18水好多| 国产精品久久久久久人妻精品电影| 亚洲专区中文字幕在线| 亚洲男人的天堂狠狠| 久久久国产精品麻豆| 日本成人三级电影网站| 欧美黄色淫秽网站| 国产精品久久久人人做人人爽| 久久精品国产综合久久久| 这个男人来自地球电影免费观看| 国产三级在线视频| 国内精品久久久久精免费| 欧美zozozo另类| 国产高清视频在线播放一区| 日本三级黄在线观看| 国产黄a三级三级三级人| 婷婷亚洲欧美| 热99re8久久精品国产| 宅男免费午夜| 欧美国产精品va在线观看不卡| av福利片在线| 亚洲国产看品久久| 亚洲av片天天在线观看| 国产亚洲av嫩草精品影院| 国产真人三级小视频在线观看| 精品第一国产精品| 欧美成人性av电影在线观看| 日韩免费av在线播放| 亚洲欧洲精品一区二区精品久久久| 琪琪午夜伦伦电影理论片6080| 色婷婷久久久亚洲欧美| 无人区码免费观看不卡| 日韩欧美三级三区| 日韩精品中文字幕看吧| 三级毛片av免费| 精品乱码久久久久久99久播| 搡老熟女国产l中国老女人| 一本久久中文字幕| 美女午夜性视频免费| 色综合婷婷激情| 给我免费播放毛片高清在线观看| 国产在线观看jvid| 精品久久久久久久人妻蜜臀av| 久久久久九九精品影院| 欧美在线黄色| 人人妻人人看人人澡| 黄片大片在线免费观看| 一a级毛片在线观看| 9191精品国产免费久久| 欧美成狂野欧美在线观看| 亚洲国产精品999在线| 麻豆一二三区av精品| 亚洲精品久久国产高清桃花| 搡老岳熟女国产| 久久 成人 亚洲| 国产精华一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 日日夜夜操网爽| www.www免费av| 精品国产一区二区三区四区第35| 中文字幕av电影在线播放|