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

    Coordinating and Evaluating of Multiple Key Performance Indicators for Manufacturing Equipment:Case Study of Distillation Column☆

    2014-07-17 09:10:24LiZhuHongyeSuShanLuYueWangQuanlingZhang

    LiZhu,Hongye Su*,Shan Lu,YueWang,Quan ling Zhang

    State Key Laboratory of Industrial Control Technology,Institute of Cyber-Sy stems and Control,Zhejiang University,Hangzhou 310027,China

    Coordinating and Evaluating of Multiple Key Performance Indicators for Manufacturing Equipment:Case Study of Distillation Column☆

    LiZhu,Hongye Su*,Shan Lu,YueWang,Quan ling Zhang

    State Key Laboratory of Industrial Control Technology,Institute of Cyber-Sy stems and Control,Zhejiang University,Hangzhou 310027,China

    A R T I C L E I N F O

    Article history:

    Received 1May 2013

    Received in revised form 24 January 2014 Accepted 3March 2014

    Available on line 20 June 2014

    Key performance indicator

    Coordination

    Manufacturing equipment

    Evaluation

    Manufacturing equipment takes the task of operation and directly effects on the manufacturing process.One single Key Performance Indicator(KPI)is mainly employed to evaluate equipment inmost studies,neither integrating the KPIs into a completed evaluation system nor considering the impact and conflict among KPIs.In this paper,a KPI evaluation architecture is presented to define and analyze KPIs,and then a common structure for KPI to obtain the KPI set of manufacturing equipment is introduced.An available multi-KPI coor dination model is proposed to discern and balance the relationship among multi-KPI.Finally,a case study is introduced to illustrate the applicability of the coor dination model by using multi-objective optimization strategy and an efficient solution is obtained.

    ?2014 Chemical Industry and Engineering Society of China,and Chemical Industry Press.All rights reserved.

    1.Introduction

    Manufacturing equipment plays a principal role in the manufacturing process and directly influences the yields and profits of manufacturing companies.There are manufacturing equipment design,operation and maintaining all the way through the entire manufacturing process.To ensure the stability and efficiency of production process,it is necessary to evaluate the performance of manufacturing equipments.

    Many investigations have been reported in the literatures about equipment evaluations.Ljung berg[1]introduced overall equipment effectiveness in the formulation and execution of a total productivity maintenance strategy to measure the performance of equipment as a metric.Yacoub and MacGregor[2]presented a data analysis method to assess the equipment performance in manufacturing processes.Chen et al.[3]proposed a performance indicator of process manufacturing time,which can be transform ed to production capacity and to evaluate equipment performance.Geng et al.[4]developed a fuzzy analytic hierarchy process method to get energy efficiency ind ices to assess energy utilization states of different equipments.Shen et al.[5]made a determination on the benchmark of manufacturing equipment and took production ratio as an instance.Younes et al.[6]presented the application of parameters design to improve both the product quality and equipment performance in a hot sheet rolling plant.Irshad et al.[7]provided total quality model and applied productivity and quality indicators to evaluate performance.Garza-Reyes et al.[8]investigated the relationship between overall equipment effectiveness and process capability for measuring the performance.

    However,most studies pay attention to use a single performance indicator to measure the manufacturing equipment without integrating the key performance indicators in to a completed evaluation system. And qualitative analysis methods of equipment performance are in the majority,which lacks systematicness,generality and accuracy.In addition,there are conflicts and inconsistence in the relationship of multi-KPI,thus it is helpful and essential to find out the trade off compromises to balance the multi-KPI,which is useful to narrow down the choices in decision-making in the manufacturing process[9].

    Considering these aspects for evaluating the performance of manufacturing equipment,evaluation architecture of key performance indicators for manufacturing executive system is introduced to define and analyze KPIs in this study.Based on the common structure of KPIs proposed by the international standard ISO 22400[10],a set of KPIs to measure manufacturing equipment performance is established and described.To discern and coordinate the relationships among the multi-KPI in this KPI set,a multi-KPI coordination model is proposed, which takes in to account coordinating objectives selected,mass balance,energy balance,quality and safety constraints,etc.The trade off among the multi-KPIis explored by collaboration model,which is not limited to evaluate the special manufacturing equipment,but also accounts for the evaluation of different levels in manufacturing process.To demonstrate the effectiveness of the presented model,themethanol-water packed distillation column is taken as a real case study. The trade off compromise between production rate and unit energy consumption is obtained by using multi-objective optimization strategy.

    2.KPI Evaluation Set of Manufacturing Equipment

    Key performance indicators(KPIs)are defined as quantifiable and strategic measurements that reflect the critical success factors in the manufacturing process.KPIs are very important for understanding, benchmarking and improving the performance of manufacturing executive system from both the manufacturing process perspective of eliminating waste and the corporate perspective of achieving strategic goals[10].Zhu,etal.[11]proposed KPI evaluation architecture and divided KPI analysis in to two parts including KPI definition and KPI utility, as shown in Fig.1.The International Organization for Standardization (ISO)22400 Automation systems and integration—key performance indicators for manufacturing operations manage men t[10]established a set of general evaluation system to express the objectives and critical factors,and proposed a common structure for standardizing KPI definition.The performance of manufacturing equipment directly influences production security,quality and efficiency,so it is necessary to evaluate manufacturing equipment from many aspects.There are different key performance indicators from different evaluation aspects,such as production ratio,unit energy consumption,quality ratio,equipment load rate,and overall equipment effectiveness,which form the KPI evaluation set of manufacturing equipment.According to KPI evaluation architecture and KPI common structure,Tables 1 and 2 give the definition and description of production ratio,and unit energy consumption.

    Production rate indicator is defined as output-input ratio and higher ratio means the efficiency is better.However,it is complicated to compute the production rate because of the interaction among the quantity of different products in the multi-input and multi-output process.For the manufacturing equipment,the ratio of target product and produced quantity is used to evaluate the operating state and production efficiency.

    Unit energy consumption indicator is used to evaluate energy consumed by equipment for energy savings,environmental protection and cost reduction.Though energy can be considered as a form of raw material,it helps to evaluate the consumption of energy using distinct indicators.

    The relationship between the good quantity and the produced quantity is denoted as quality ratio.The higher the quality ratio,the more products meet the quality requirement.The information about the ratio of produced quantity in relation to the maximum equipment production capacity is provided by equipment load rate indicator,which reflects production efficiency,technical performance,equipment production state and equipment utilization by researching the usage of equipment.The value of equipment load rate impacts the production cost and profit.The optimal results of manufacturing process in the nondisturbance condition assess production loss and improve product quality represented by over all equipment effectiveness,which could be applied to evaluate either single production equipment or production unit which consists of multiple production equipments.These indicators could be defined and described on the basis of KPI common structure, and form the KPI set for manufacturing equipment with other indicators.

    3.Multi-KPI Coordination Model of Manufacturing Equipment

    According to the KPI evaluation architecture and common structure, one single KPI reflects the effectiveness of the certain concern point in the equipment operation state.However,due to managing the manufacturing equipment more effectively and efficiently,it is necessary to analyze and coordinate the relationships among multi-KPI in the KPI evaluation set,as shown in Fig.2.

    The relationships of KPIs are complicated and they interact and impact each other in the KPI set,which affects decision making.For instance,production rate may have a direct effect on the overall equipment effectiveness of the manufacturing equipment[8].Raising production rate while simultaneously reducing unit energy consumption is often conflicted and inconsistent.Therefore,finding out the trade off compromises to keep the multi-KPI balance is crucial.According to the multi-objective optimization methodology,the multi-KPI coordination model of manufacturing equipment contains two parts:coordinate objectives and constraint condition.

    3.1.Coordinate objectives

    Minimizing or maximizing the focused KPIs is the coordinate objectives.The optimal overall objective is obtained by balancing the relationship among KPIs.

    Fig.1.Evaluation architecture of key performance indicators.

    Table 1Production ratio indicator

    whereχ1,χ2…,χmand y1,y2,…ynare variables for calculating each KPI, such as production yield for calculating production rate indicator,and fault time for computing mean operating time between failure indicator.

    3.2.Constraint condition

    To discern the relation among KPIs,it is required to analyze the effect of variables on KPIs in the constraint condition,which can be considered from the aspects of mass balance,energy balance,quality and safety constraints,etc.Mass balance is closely related to equipment performance and production process.In broad term s,the input material is equal to the output material for the certain equipment.Meanwhile, energy transfer takes p lace in the production process,which includes complicated physical and chemical changes.Hence,it is necessary to keep balance between produced and consumed energy.The quality and safety require men ts would be taken in to consideration in the process to ensure economical and social benefits for equipment.

    where z1,z2,…,zkare variables for the special equipment.For instance, material quantity,purity,energy medium quantity,temperature, pressure,etc.

    Table 2Unit energy consumption indicator

    4.Multi-KPI Coordinate Model Application for Distillation Column

    Distillation process is widely applied in the production of metallurgy and chemical enterprise and it is a unit operation with high energy consumption.It is of great practical significance for gaining anticipative economical and social benefits by keeping the balance between production yield and energy consumption.For the typical distillation equipment, production rate and unit energy consumption are the key performance indicators for evaluating the production process of distillation column. In the case of methanol-water packed distillation column,coordinating the relation between production rate and unit energy consumption is a typical way to evaluate the equipment on the condition of quality standard.

    4.1.Distillation process description

    In the distillation process,the input materials are separated to satisfy the desired purity through repeated vaporization and condensation by pre-heater,distillation column,condenser,re-boiler,re flux drum,etc. [12].For themethanol-water packed distillation column,the raw material flows in to distillation column after being p reheated by pre-heater. The rising steam is condensed in to liquid by overhead condenser and divided in to product and re flux.The production methanol is kept in the storage after cooling.During this process,the steam and water are taken as heating and cooling energy medium to meet the energy demand by energy transfer,respectively.Fig.3 shows the distillation process.

    F is the quantity of input material which contains methanol and water,B is the production quantity of column bottom,and D is the quantity of high purity methanol product.χF,χB,χDare the purity of input material,production of column and methanol production,respectively.L means the re flux quantity and V stands for the steam quantity in the column.

    Thevalues and hard limits of themethanol-water packed distillation column[13,14]are given in Tables 3 and 4.

    4.2.Formulation of the KPI coordination model for distillation column

    4.2.1.Coordinate objectives

    Two key performance indicators of production rate and unit energy consumption are considered to evaluate the distillation column in this case study.Themethanol-water packed distillation column is used for producing the methanol which satisfies the quantity requirement, while in this distillation process,the energy is consumed,such as the heating for re-boiler and the cooling for condenser.In order to raise the production rate,the unit energy consumption is increased.Thus,it is necessary to regulate these KPIs to achieve the balance.

    The production rate indicator of distillation is described as

    where D and B are the production of column top and bottom, respectively.

    The unit energy consumption indicator is expressed as

    where Wsteamis the heating steam quantity used in the re-boiler and Wwateris the cooling water quantity consumed in the condenser.

    Given the profit of the production and the cost of energy consumption,there are two in dependent objective functions of the coordinated model in this paper,namely,maximizing the production rate andminimizing the unit energy consumption at the same time,which is form ulated m athem atically as follows,

    Fig.2.Relationship among multi-KPI in the KPIset.

    Fig.3.Methanol-water packed distillation column.

    4.2.2.Mass balance constraints

    When the distillation column reaches the steady state,the input material and the outputmaterial achieve the balance.Thus,them ass balance constraint of the entire column is

    Them ass balance constraint of the product methanol is considered as

    The rising steam is congealed totally by the overhead condenser and divided into top column product and re flux.Hence,mass balance of the overhead re flux d rum is

    4.2.3.Energy balance constraints

    There are conversion phenomenon between vapor and liquid in the distillation.Because of the change of re flux ratio R=l/d,energy balance constrain t of gas and liquid is

    The energy balance constraint of the whole reactor-distillation is expressed as follows in the steady state,

    Table 3Nominal values for the variables of the reaction-distillation column

    Table 4Hard lim its on the flow rates of the reactor-distillation column

    4.2.4.Quality constraints

    The product methanol quality is the focus of the distillation process. According to the quality requirement,the quality constraints are

    4.2.5.Safety constraints

    In order to ensure safety in production,it is necessary to keep the equipment meet capacity requirem en ts,

    4.3.Solution strategy and result discussion

    In order to raise the production rate and reduce the unit energy consumption,the variables that affect the relationship between the two KPIs are analyzed using mechanism method,including methanol production quantity from column top(D),byproduct production quantity from column bottom(B),steam quantity in the column(V),re flux ratio(R),methanol production pu rity(χD),byproduct purity(χB), heating steam quantity(Wsteam)and cooling water quantity(Wwater). According to the characteristics of the coordinate model and the requirement of methanol production quantity,methanol production quantity from column top(D)and methanol production purity(χD) are considered as the main decision variables for optimization in this study.

    Due to the relationship between production rate and energy consumption,the appropriate set point value of one variable for maximizing the production rate may not be suitable for minimizing the unit energy consumption.Hence,the applicable trade-off solutions for both two KPI objectives should be considered.

    It is appropriate to solve the multi-KPI coordinate problem by using multi-objective optimization strategy,which aims at searching for one or more satisfying solutions in Pareto optimal set.Multi-objective optimization strategy contains two main algorithm s:one is converting the multi-objective problems to a single objective problem using some methods,and another is selecting the satisfying solutions using some trade-off criterion in the Pareto optimal set.To get the Pareto optimal set,mathematical programming approach and genetic algorithm(GA) are applicable[15,16].Following the vector evaluated genetic algorithm (VEGA)first proposed by S chaffer in 1985[17],many multi-objective evolutionary algorithms are studied,such as multiple objectives with particle swarm optimization algorithm(MOPSO)[18],the nondominate sorting genetic algorithm(NSGA)and the non-dominatesorting genetic algorithm II(NSGA-II)[19].In this study,εrestriction method is employed to acquire and evaluate the Pareto solution set due to containing strong nonlinear constraints[20].

    Table 5Distillation process parameters of Pareto optimal solution

    Fig.4.Pareto carve obtained by ε restriction method for the coordinated model.

    For this study,the production rate is considered to the critical object because it directly influences the enterprise's income,while unit energy consumption is taken as one of constraints,less than ε.The range of energy consumption is[UECmin,UECmax],where UECminand UECmaxare the minimum and maximum unit energy consumption respectively. In this paper,the minimum and maximum unit energy consumption is the solution of optimization model under the present constraints.The solution byε restriction method is given in Table 5.

    An efficient Pareto curve is generated by the optimal solutions of optimization problem.The conflict between the effects of the decision variables on the two KPI coordinated objective functions results in the optimum being a Pareto optimal set rather than a unique solution.The characteristic of the Pareto set for this studied case is that the two objects could not achieve the satisfied results.When one point in the Pareto optimal set is moved to another,one objective function is improved but the other becomes worse.Therefore,there is no dominated solution within the Pareto optimal set.The additional in formation like KPI priority is required to choose the p referred solution from the whole Pareto optimal set for operation[21-23].The obtained Pareto carve is shown in Fig.4.

    Fig.4 describes the relationship between the two KPIs,production rate and unit energy consumption,which provides a better perspective for operator to make an effective decision.The dotted points in Fig.4 give the maximum and minimum production rate with the change of unit energy consumption within the present constraints.It is shown that the rising of the production rate causes the increasing of unit energy consumption correspondingly.However,the change of gradient is not as enormous as it is in the beginning.Due to the equipment capacity bound and quality requirement,production rate remains the same with the growing of unit energy consumption.

    Each dotted point on the Pareto carve is an optimal solution under the present constraints,which is related with a set of decision variables. And according to the present constraints and requirements,three key points from the corresponding objective values of the solutions are listed in detail in Table 6,respectively.

    From this table,it is observed that there is a trade off set between maximizing production rate and minimizing unit energy consumption. Under the quality requirements,the energy consumption raises simultaneously with the growing of production yields,while the methanol production purity decreases.There are different emphasis points and analysis view s for these three cases.Case A produces high purity product preferred by the market and low unit energy consumption cost, but low production rate leads to less total benefit.On the contrast, Case C produces more low purity product with abundant heating steam and cooling water,and also receives little total benefit.Case B falls in between Case A and Case C,and though production rate does not reach the maximum and unit energy consumption is high,more benefit is achieved.Considering the actual constraints and require men ts,decision-makers balance the relation between the two KPIs and select one operating point to get the maximum benefit.

    5.Conclusions

    Manufacturing equipment plays an important role in the manufactuing process.It is useful and helpful to evaluate manufacturing equipment for improving its capability.The paper addresses the challenge of key performance indicator evaluation.Th rough describing and analyzing each KPI,the KPI evaluation set for manufacturing equipmentbased on KPIs evaluation architecture and common structure is obtained. To discern the relationship among KPIs,an available coordination model is presented and a case study about methanol-water packed distillation column is used to demonstrate the effectiveness.With coordinating the relation between production ratio and unit energy consumption based on multi-objective optimization strategy,an efficient solution set denoted as Pareto set is gained.The relation between these KPIs is explored and used for choosing the appropriate solution from the Pareto set. Furthermore,the coordination model and strategy are not only limited to evaluate the special manufacturing equipment,but also accounts for the evaluation of different levels in manufacturing process.

    Table 6Comparison between coordinated objectives and decision variables from three selected cases

    Nomenclature

    B production quantity of column bottom,kg·h?1

    b the molar mass of production from the column bottom, kmol·h?1

    Cpfmeththe specific heat of methanol,kJ·kg?1·°C?1

    Cpfwatthe specific heat of water,kJ·kg?1·°C?1

    D the quantity of high purity methanol product,kg·h?1

    d the molar mass of high purity methanol product,kmol·h?1

    E energy consumed in the manufacturing process

    Ef plate efficiency

    F the quantity of input material,t·a?1

    f the molar mass of input material,km ol·h?1

    ΔHmethanolevaporation latent heat of methanol,k J·kg?1

    ΔHsteamevaporation latent heat of steam,k J·kg?1

    ΔHwaterevaporation latent heat of water,kJ·kg?1

    L the re flux quantity,kg·h?1

    l the molar mass of re flux,kg·h?1

    M material used in the manufacturing process

    Mbaverage relative molecular mass of the production from the column bottom,g·m ol?1

    Mdaverage relative molecular mass of the high purity methanol product,g·m ol?1

    Mfaverage relative molecular mass of input material,g·m ol?1

    m ax maximum

    m in minimum

    n number of stages

    PR production rate

    Q quality requirement

    QBheat load of reboiler,k J·h?1

    QCheat load of condenser,kJ·h?1

    QC1heat load of chiller,kJ·h?1

    QDheat load of column top,kJ·h?1

    QFheat load of p reheater,kJ·h?1

    Qsheat load of steam,k J·h?1

    R re flux ratio,%

    S safety requirement

    t temperature,°C

    UEC unit energy consumption

    V the steam quantity in the column,kg·h?1

    v the molar mass of the steam in the column,km ol·h?1

    Wsteamthe heating steam quantity used in the re-boiler,kg·h?1

    Wwaterthe cooling water quantity consumed in the condenser,kg·h?1

    χBthe purity of production in the column bottom,%

    χbthe mole fraction of production in the column bottom,%

    χDthe purity of methanol production,%

    χdthe mole fraction of methanol production,%

    χFthe purity of input material,%

    χfthe mole fraction of input material,%

    χidata and variables for calculating each KPI

    yidata and variables for calculating each KPI

    zivariables or data for the special equipment

    αrelative volatility

    ηheat loss ratio in the column bottom,%

    [1]O.Ljungberg,Measurement of overall equipment effectiveness as a basis for TPM activities,Int.J.Oper.Prod.Manag.18(5)(1998)495-507.

    [2]F.Yacoub,J.F.MacGregor,Analysis of equipment performance in processes with complex serial-parallel structures,Ind.Eng.Chem.Res.48(2009)7181-7185.

    [3]K.S.Chen,M.L.Huang,P.L.Chang,Performance evaluation on manufacturing times, Int.J.Adv.Manuf.Technol.31(2006)335-341.

    [4]Z.Q.Geng,Y.M.Han,X.B.Gu,Q.X.Zhu,Energy efficiency estimation based on data fusion strategy:case study of ethy lene product industry,Ind.Eng.Chem.Res.51 (25)(2012)8526-8534.

    [5]Q.H.Shen,H.Y.Su,L.Zhu,S.Lu,Benchmark determination for key performance indicators ofmanufacturing equipment,CIESC J.63(9)(2012)2958-2964.

    [6]M.A.Younes,M.Shahtout,M.N.Damir,A parameters design approach to improve product quality and equipment performance in hot rolling,J.Mater.Process.Technol. 171(2006)83-92.

    [7]A.S.Irshad,Y.Jam il,K.M.Rau f,A.M.Syed,Evaluation of performance in manufacturing organization through productivity and quality,Afr.J.Bus.Manag.5(6)(2011) 2211-2219.

    [8]J.A.Garza-Reyes,S.Eld ridge,K.D.Barber,H.Soriano-Meier,Overall equipment effectiveness(OEE)and process capability(PC)measures,Int.J.Qual.Reliab.Manag.27 (1)(2010)48-62.

    [9]C.C.Yuen,S.K.Aatmeeyata Gupta,A.K.Ray,,Multi-objective optimization of membrane separation modules using genetic algorithm,J.Mem br.Sci.176(2000) 177-196.

    [10]ISO/DIS 22400-2.“Automation systems and integration—key performance indicators for manufacturing operations management-Part2:Definitions and descriptions”.

    [11]L.Zhu,H.Y.Su,Q.H.Shen,Evaluation architecture of manufacturing execution system based on key performance indicators in process industry,Com put.Integr. Manuf.Syst.18(12)(2012)2643-2649.

    [12]S.Q.Wang,L.K.Dai,L.Yu,Process Control Engineering,Chemical Industry Press, Beijing,2008.(in Chinese).

    [13]T.Tosukhowong,Dynamic Real-time Optimization and Control of an Integrated Plant,Ph.D.Thesis,Georgia Institute of Technology,Notrh Avenue,2006.

    [14]Y.Sh i,Chemical Engineering Manual,Chemical Industry Press,Beijing,1996. (in Chinese).

    [15]S.J.Li,H.Wang,F.Qian,Multi-objective genetic algorithm and its applications in chemical engineering,Com put.Appl.Chem.20(6)(2003)755-760(in Chinese).

    [16]F.Xu,X.L.Luo,Margin analysis and control design of FCCU regenerator based on dynamic optimization(II)Solution and resu ltanalysis,CIESCJ.60(3)(2009)683-690.

    [17]J.D.Schaffer,Multi-objective optimization with vector evaluated genetic algorithm s, Proceedings of the International Conference on Genetic Algorithm s and their Application,1985,pp.93-100,(Pittsburgh).

    [18]C.A.Coello,G.T.Pulido,M.S.Lechuga,Handing multiple objectives with particle swarm optimization,IEEETrans.Evol.Com put.8(3)(2004)256-279.

    [19]K.Deb,A.Pratap,S.Agrawal,et al.,A fast and elitist multi-objective genetic algorithm:NSGA-II,IEEETrans.Evol.Com put.6(2)(2002)182-197.

    [20]C.Schweiger,Process Synthesis,Design,and Control:Optimization with Dynamic Models and Discrete Decisions,Princeton University,Princeton,1999.

    [21]S.J.Mu,H.Y.Su,Y.Gu,J.Chu,Multi-objective optimization of industrial purified terephthalic acid oxidation process,Chin.J.Chem.Eng.11(5)(2003)536-541.

    [22]W.F.Hou,H.Y.Su,S.J.Mu,J.Chu,Multi-objective optimization of the industrial naphtha catalytic reforming process,Chin.J.Chem.Eng.15(1)(2007)75-80.

    [23]Y.Q.Jiao,H.Y.Su,Z.W.Liao,W.F.Hou,Modeling and multi-objective optimization of refinery hydrogen network,Chin.J.Chem.Eng.19(6)(2011)990-998.

    ☆Supported by the National Natural Science Foundation of China(61134007, 61320106009).

    *Corresponding author.

    E-mailaddress:hysu@iipc.zju.edu.cn(H.Su).

    青春草视频在线免费观看| 不卡视频在线观看欧美| 人妻系列 视频| 久久久久精品性色| 91久久精品国产一区二区三区| 成年女人毛片免费观看观看9 | 可以免费在线观看a视频的电影网站 | 久久精品夜色国产| 亚洲欧洲国产日韩| 成人手机av| 精品国产一区二区三区久久久樱花| 久久女婷五月综合色啪小说| 亚洲欧美精品自产自拍| 日本-黄色视频高清免费观看| 国产免费现黄频在线看| 久久女婷五月综合色啪小说| 黑丝袜美女国产一区| 欧美日本中文国产一区发布| 久久久久久久久久久免费av| 午夜福利在线免费观看网站| 熟女电影av网| 欧美精品人与动牲交sv欧美| 亚洲av国产av综合av卡| 国产精品无大码| 一区二区av电影网| 国产国语露脸激情在线看| 18禁裸乳无遮挡动漫免费视频| 亚洲内射少妇av| 久久99蜜桃精品久久| 精品国产一区二区久久| 国产一区亚洲一区在线观看| 国产亚洲最大av| 久久97久久精品| www日本在线高清视频| 最近的中文字幕免费完整| 在线观看www视频免费| 三上悠亚av全集在线观看| 老女人水多毛片| 考比视频在线观看| 80岁老熟妇乱子伦牲交| 99re6热这里在线精品视频| 黑丝袜美女国产一区| 成年美女黄网站色视频大全免费| 亚洲激情五月婷婷啪啪| 91午夜精品亚洲一区二区三区| 看非洲黑人一级黄片| 免费久久久久久久精品成人欧美视频| 国产亚洲最大av| 综合色丁香网| 色播在线永久视频| 欧美 亚洲 国产 日韩一| 侵犯人妻中文字幕一二三四区| 亚洲伊人色综图| 精品一品国产午夜福利视频| 国产成人精品福利久久| 男女下面插进去视频免费观看| 亚洲精品日韩在线中文字幕| 精品一区二区免费观看| 爱豆传媒免费全集在线观看| 久久久久人妻精品一区果冻| 中文字幕av电影在线播放| 久久国产亚洲av麻豆专区| 午夜老司机福利剧场| 亚洲人成电影观看| 国产精品不卡视频一区二区| 十八禁网站网址无遮挡| 美女脱内裤让男人舔精品视频| 久久女婷五月综合色啪小说| 一级毛片黄色毛片免费观看视频| 交换朋友夫妻互换小说| 欧美日韩一级在线毛片| 一级片'在线观看视频| 欧美精品av麻豆av| 亚洲色图综合在线观看| 亚洲精品日本国产第一区| 丰满饥渴人妻一区二区三| 乱人伦中国视频| 男女无遮挡免费网站观看| 欧美+日韩+精品| 国产成人精品在线电影| 老鸭窝网址在线观看| 最黄视频免费看| 可以免费在线观看a视频的电影网站 | 成年人免费黄色播放视频| 99国产综合亚洲精品| 1024视频免费在线观看| 日韩伦理黄色片| 亚洲av福利一区| 天堂8中文在线网| 丝袜美足系列| 两个人看的免费小视频| 9色porny在线观看| 91精品国产国语对白视频| 亚洲一区中文字幕在线| 如何舔出高潮| 美女午夜性视频免费| 亚洲精品美女久久久久99蜜臀 | av国产精品久久久久影院| 久久午夜福利片| 少妇 在线观看| 日日啪夜夜爽| 亚洲,欧美,日韩| 天天操日日干夜夜撸| 最黄视频免费看| 欧美另类一区| 精品国产国语对白av| 国产精品一二三区在线看| 久久国产精品男人的天堂亚洲| 老汉色av国产亚洲站长工具| 国产福利在线免费观看视频| 99热国产这里只有精品6| 国产亚洲最大av| 2021少妇久久久久久久久久久| 热re99久久精品国产66热6| 国产又爽黄色视频| 免费大片黄手机在线观看| 亚洲精品aⅴ在线观看| 精品一区二区三区四区五区乱码 | 久久毛片免费看一区二区三区| 色网站视频免费| 亚洲美女黄色视频免费看| 女人精品久久久久毛片| 黄片播放在线免费| 免费日韩欧美在线观看| 18禁动态无遮挡网站| 欧美在线黄色| 成人免费观看视频高清| freevideosex欧美| 一本色道久久久久久精品综合| 久久热在线av| 亚洲国产色片| 久久精品国产a三级三级三级| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日本vs欧美在线观看视频| 天堂8中文在线网| 亚洲欧美一区二区三区国产| 丁香六月天网| www.av在线官网国产| 日韩一本色道免费dvd| 各种免费的搞黄视频| 日日爽夜夜爽网站| 丰满饥渴人妻一区二区三| 午夜福利乱码中文字幕| 免费观看在线日韩| 亚洲欧洲国产日韩| 国产精品女同一区二区软件| 久久久久网色| 亚洲少妇的诱惑av| 欧美人与善性xxx| 永久免费av网站大全| 国产精品秋霞免费鲁丝片| 91久久精品国产一区二区三区| 女性被躁到高潮视频| 久久久久久久国产电影| 亚洲第一青青草原| av网站免费在线观看视频| a级毛片在线看网站| 肉色欧美久久久久久久蜜桃| 久久这里只有精品19| 国产精品女同一区二区软件| 97人妻天天添夜夜摸| 久久久久久人妻| 欧美激情 高清一区二区三区| 日韩一本色道免费dvd| 一级片'在线观看视频| 亚洲国产欧美网| 妹子高潮喷水视频| 蜜桃在线观看..| 午夜福利在线观看免费完整高清在| 日韩精品有码人妻一区| 午夜免费鲁丝| 最近2019中文字幕mv第一页| 国产精品麻豆人妻色哟哟久久| 女性生殖器流出的白浆| 两性夫妻黄色片| 午夜久久久在线观看| 亚洲国产精品成人久久小说| 波多野结衣av一区二区av| 伦理电影大哥的女人| a级毛片在线看网站| 亚洲av福利一区| 亚洲综合精品二区| 999久久久国产精品视频| 亚洲久久久国产精品| 久久99精品国语久久久| 夫妻午夜视频| 国产老妇伦熟女老妇高清| 新久久久久国产一级毛片| 丝袜美腿诱惑在线| 26uuu在线亚洲综合色| 青青草视频在线视频观看| 黄网站色视频无遮挡免费观看| 欧美变态另类bdsm刘玥| 精品福利永久在线观看| 成人国语在线视频| 中文字幕精品免费在线观看视频| 午夜福利视频在线观看免费| 美女高潮到喷水免费观看| 五月伊人婷婷丁香| 免费观看性生交大片5| 丰满饥渴人妻一区二区三| 尾随美女入室| 亚洲av国产av综合av卡| 午夜激情久久久久久久| 国产在视频线精品| 九色亚洲精品在线播放| 国产精品久久久av美女十八| tube8黄色片| 看十八女毛片水多多多| 国产精品无大码| 国产av国产精品国产| 日韩中文字幕欧美一区二区 | 视频区图区小说| 黑人欧美特级aaaaaa片| 免费看av在线观看网站| 国产探花极品一区二区| 香蕉国产在线看| 久久热在线av| 另类精品久久| 亚洲在久久综合| 少妇被粗大猛烈的视频| 91精品国产国语对白视频| 男人舔女人的私密视频| 在线精品无人区一区二区三| 欧美日韩一级在线毛片| 久久韩国三级中文字幕| 亚洲美女搞黄在线观看| 波野结衣二区三区在线| 日韩一卡2卡3卡4卡2021年| 婷婷色综合www| 波多野结衣av一区二区av| 国产熟女欧美一区二区| 色吧在线观看| 国产老妇伦熟女老妇高清| 青春草视频在线免费观看| 亚洲av欧美aⅴ国产| √禁漫天堂资源中文www| 一级a爱视频在线免费观看| 国产精品国产av在线观看| 亚洲激情五月婷婷啪啪| 激情五月婷婷亚洲| 一级黄片播放器| 人人妻人人爽人人添夜夜欢视频| 亚洲色图 男人天堂 中文字幕| 国精品久久久久久国模美| 中文字幕亚洲精品专区| 成人亚洲欧美一区二区av| 十分钟在线观看高清视频www| 天天躁日日躁夜夜躁夜夜| 巨乳人妻的诱惑在线观看| 交换朋友夫妻互换小说| 成人影院久久| 国产亚洲一区二区精品| 99国产精品免费福利视频| 汤姆久久久久久久影院中文字幕| 一级毛片电影观看| 欧美精品亚洲一区二区| 欧美老熟妇乱子伦牲交| av在线观看视频网站免费| 欧美国产精品一级二级三级| 9热在线视频观看99| 精品国产一区二区三区四区第35| 男女边摸边吃奶| 天美传媒精品一区二区| 亚洲av综合色区一区| 亚洲中文av在线| 亚洲av男天堂| 水蜜桃什么品种好| 中文字幕人妻丝袜一区二区 | 99热国产这里只有精品6| 亚洲欧美一区二区三区久久| 国产又爽黄色视频| 亚洲欧美精品自产自拍| 男女边吃奶边做爰视频| 黄网站色视频无遮挡免费观看| 不卡视频在线观看欧美| 欧美人与性动交α欧美精品济南到 | 国产片内射在线| 99热全是精品| 在线天堂中文资源库| 亚洲综合色惰| 久久久久网色| 夜夜骑夜夜射夜夜干| 亚洲图色成人| 中文字幕人妻丝袜制服| av网站在线播放免费| 麻豆精品久久久久久蜜桃| 欧美日本中文国产一区发布| 日韩制服骚丝袜av| 免费在线观看视频国产中文字幕亚洲 | 男人舔女人的私密视频| 丝袜美腿诱惑在线| 日韩中文字幕视频在线看片| 亚洲国产色片| 亚洲男人天堂网一区| 免费在线观看黄色视频的| av电影中文网址| 秋霞伦理黄片| 午夜日本视频在线| 国产av国产精品国产| 涩涩av久久男人的天堂| 亚洲国产精品一区二区三区在线| 亚洲国产精品一区三区| 国产深夜福利视频在线观看| 欧美激情高清一区二区三区 | 亚洲色图综合在线观看| 一个人免费看片子| 亚洲精品久久久久久婷婷小说| 免费观看无遮挡的男女| 丁香六月天网| 亚洲一码二码三码区别大吗| 男女高潮啪啪啪动态图| 国产一级毛片在线| 国产精品一区二区在线观看99| 亚洲av.av天堂| 亚洲精品第二区| 成年人免费黄色播放视频| 99久久中文字幕三级久久日本| a 毛片基地| 欧美+日韩+精品| 国产成人91sexporn| 国产精品久久久久成人av| av视频免费观看在线观看| 国产精品偷伦视频观看了| 久久 成人 亚洲| 久久这里有精品视频免费| 午夜av观看不卡| av国产久精品久网站免费入址| 两个人免费观看高清视频| 91aial.com中文字幕在线观看| 免费观看无遮挡的男女| 国产亚洲精品第一综合不卡| 国产免费现黄频在线看| 一区在线观看完整版| 久久精品国产综合久久久| 国产免费又黄又爽又色| 久久久亚洲精品成人影院| 精品一区二区三区四区五区乱码 | 亚洲国产成人一精品久久久| 国产有黄有色有爽视频| 一级毛片电影观看| 自线自在国产av| 伊人亚洲综合成人网| 亚洲国产精品一区三区| 三级国产精品片| 亚洲,欧美,日韩| 色婷婷av一区二区三区视频| 午夜激情av网站| 久久久久精品人妻al黑| 啦啦啦在线观看免费高清www| 久久久精品区二区三区| 久久久久久久久久久免费av| 777米奇影视久久| 欧美成人精品欧美一级黄| 亚洲精品美女久久av网站| 交换朋友夫妻互换小说| 男女啪啪激烈高潮av片| 少妇人妻 视频| 欧美少妇被猛烈插入视频| 中文乱码字字幕精品一区二区三区| 国产精品久久久久久精品古装| 亚洲人成网站在线观看播放| 久久精品人人爽人人爽视色| 高清黄色对白视频在线免费看| 日本wwww免费看| 久久久精品国产亚洲av高清涩受| 久久国产亚洲av麻豆专区| 国产免费又黄又爽又色| 天天躁夜夜躁狠狠久久av| 亚洲,一卡二卡三卡| 久久狼人影院| 日韩中字成人| 色婷婷久久久亚洲欧美| 看免费成人av毛片| av网站在线播放免费| 精品99又大又爽又粗少妇毛片| 欧美精品人与动牲交sv欧美| 国产av国产精品国产| 国产xxxxx性猛交| 国产精品无大码| 欧美人与性动交α欧美软件| 欧美黄色片欧美黄色片| 日韩一卡2卡3卡4卡2021年| 日韩中文字幕欧美一区二区 | 丰满饥渴人妻一区二区三| 夫妻性生交免费视频一级片| 大片电影免费在线观看免费| 自线自在国产av| 日产精品乱码卡一卡2卡三| 99久久综合免费| 人妻一区二区av| av又黄又爽大尺度在线免费看| 在线观看www视频免费| 欧美精品人与动牲交sv欧美| 99九九在线精品视频| 高清视频免费观看一区二区| 可以免费在线观看a视频的电影网站 | 中文精品一卡2卡3卡4更新| 性色av一级| 高清黄色对白视频在线免费看| 91在线精品国自产拍蜜月| 在线观看免费视频网站a站| 久久久精品国产亚洲av高清涩受| 麻豆精品久久久久久蜜桃| 欧美激情高清一区二区三区 | 日韩免费高清中文字幕av| 18禁国产床啪视频网站| 成人国产av品久久久| 女人精品久久久久毛片| 少妇熟女欧美另类| 90打野战视频偷拍视频| 一边亲一边摸免费视频| 久久久久久久国产电影| 一区二区三区精品91| 午夜福利乱码中文字幕| 五月伊人婷婷丁香| 成年人午夜在线观看视频| 午夜福利在线免费观看网站| 成人18禁高潮啪啪吃奶动态图| 少妇精品久久久久久久| 欧美日韩精品成人综合77777| 国产精品免费视频内射| 免费观看在线日韩| 日本爱情动作片www.在线观看| 80岁老熟妇乱子伦牲交| 国产97色在线日韩免费| 亚洲国产欧美日韩在线播放| a级毛片在线看网站| 日本欧美国产在线视频| 波野结衣二区三区在线| 日本av手机在线免费观看| 久久这里有精品视频免费| 色婷婷av一区二区三区视频| 欧美成人精品欧美一级黄| 久久久亚洲精品成人影院| 免费在线观看视频国产中文字幕亚洲 | 精品一品国产午夜福利视频| 国产成人91sexporn| 另类精品久久| a级毛片在线看网站| 一本大道久久a久久精品| 岛国毛片在线播放| av网站在线播放免费| 精品亚洲成a人片在线观看| 亚洲美女黄色视频免费看| 国产老妇伦熟女老妇高清| 黑丝袜美女国产一区| 日日爽夜夜爽网站| 亚洲男人天堂网一区| 国产亚洲午夜精品一区二区久久| 丁香六月天网| 免费播放大片免费观看视频在线观看| 韩国av在线不卡| 亚洲欧美精品自产自拍| 久久韩国三级中文字幕| 欧美黄色片欧美黄色片| av免费观看日本| 黄网站色视频无遮挡免费观看| 边亲边吃奶的免费视频| 热re99久久国产66热| 亚洲第一青青草原| 色哟哟·www| 青春草视频在线免费观看| 极品人妻少妇av视频| 建设人人有责人人尽责人人享有的| 午夜福利网站1000一区二区三区| 精品少妇内射三级| 久久久久国产网址| 9色porny在线观看| 日韩在线高清观看一区二区三区| 日韩av免费高清视频| 国产野战对白在线观看| 日本wwww免费看| 欧美成人精品欧美一级黄| 如日韩欧美国产精品一区二区三区| 亚洲欧美精品综合一区二区三区 | 婷婷色综合大香蕉| 国产成人aa在线观看| 欧美精品一区二区大全| 咕卡用的链子| 免费在线观看完整版高清| 国产人伦9x9x在线观看 | 丰满乱子伦码专区| 精品一区二区三卡| 大码成人一级视频| 欧美精品高潮呻吟av久久| 国产成人免费观看mmmm| 久久人人97超碰香蕉20202| 久久精品国产自在天天线| 国产一级毛片在线| 只有这里有精品99| 麻豆乱淫一区二区| 欧美日韩成人在线一区二区| 最近2019中文字幕mv第一页| 男女啪啪激烈高潮av片| 在线免费观看不下载黄p国产| 亚洲一区二区三区欧美精品| 一级a爱视频在线免费观看| 视频在线观看一区二区三区| 亚洲精品第二区| 亚洲激情五月婷婷啪啪| 九色亚洲精品在线播放| 日韩一本色道免费dvd| 侵犯人妻中文字幕一二三四区| 精品人妻偷拍中文字幕| 黄色视频在线播放观看不卡| 国产精品 国内视频| 人妻少妇偷人精品九色| 老鸭窝网址在线观看| 人妻少妇偷人精品九色| 精品人妻偷拍中文字幕| 欧美成人午夜精品| 国产精品成人在线| 国产精品一二三区在线看| 久久精品亚洲av国产电影网| 另类精品久久| 国产熟女午夜一区二区三区| 欧美 日韩 精品 国产| 观看av在线不卡| 九九爱精品视频在线观看| 国产成人91sexporn| 波多野结衣av一区二区av| 日韩一卡2卡3卡4卡2021年| 精品亚洲成a人片在线观看| 中文字幕色久视频| 亚洲三级黄色毛片| 多毛熟女@视频| 国产精品国产av在线观看| 欧美国产精品一级二级三级| a级片在线免费高清观看视频| 亚洲av电影在线观看一区二区三区| a级片在线免费高清观看视频| 亚洲情色 制服丝袜| 老汉色av国产亚洲站长工具| 啦啦啦视频在线资源免费观看| 最近中文字幕2019免费版| 亚洲欧美一区二区三区久久| 寂寞人妻少妇视频99o| 狂野欧美激情性bbbbbb| 欧美日韩成人在线一区二区| 老汉色∧v一级毛片| 国产极品粉嫩免费观看在线| 国产乱人偷精品视频| 免费不卡的大黄色大毛片视频在线观看| 久久久a久久爽久久v久久| 少妇的丰满在线观看| 成年动漫av网址| 中文字幕最新亚洲高清| 国产精品二区激情视频| 国产av国产精品国产| 国产一区二区激情短视频 | 久久精品久久久久久久性| 美女主播在线视频| 99久久人妻综合| 中文字幕色久视频| 日韩电影二区| 国产片内射在线| 99国产精品免费福利视频| 日日摸夜夜添夜夜爱| 国产成人精品在线电影| 一本大道久久a久久精品| 欧美xxⅹ黑人| 欧美精品亚洲一区二区| 丝袜人妻中文字幕| 亚洲成人一二三区av| 日韩精品有码人妻一区| 国产野战对白在线观看| 欧美精品国产亚洲| 久久这里有精品视频免费| av天堂久久9| 中文字幕人妻丝袜一区二区 | 不卡av一区二区三区| 成人毛片a级毛片在线播放| 国产精品久久久久久精品电影小说| 一级,二级,三级黄色视频| 超碰97精品在线观看| 欧美bdsm另类| 丁香六月天网| 777久久人妻少妇嫩草av网站| 汤姆久久久久久久影院中文字幕| 一区二区三区激情视频| av福利片在线| 丁香六月天网| 熟女少妇亚洲综合色aaa.| 亚洲欧美成人精品一区二区| 看免费av毛片| 国产免费现黄频在线看| 久久国内精品自在自线图片| 在线看a的网站| 亚洲精品第二区| 欧美日韩一级在线毛片| 十八禁高潮呻吟视频| 日韩av免费高清视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩中文字幕欧美一区二区 | 国产在线一区二区三区精| 亚洲精品乱久久久久久| 日韩一区二区视频免费看| 黄片播放在线免费| 久久久国产精品麻豆| 免费在线观看视频国产中文字幕亚洲 | 一区福利在线观看| 人人澡人人妻人| 波野结衣二区三区在线| 在线 av 中文字幕| 久久精品亚洲av国产电影网| 九九爱精品视频在线观看| 18禁裸乳无遮挡动漫免费视频| 欧美变态另类bdsm刘玥| 亚洲精品av麻豆狂野| 久久人妻熟女aⅴ| 99热全是精品| 精品一区二区三区四区五区乱码 | 精品国产乱码久久久久久小说| 极品少妇高潮喷水抽搐|