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

    Cutting tool temperature prediction method using analytical model for end milling

    2016-11-23 01:57:32WuBohiCuiDiHeXiodongZhngDinghuTngKi
    CHINESE JOURNAL OF AERONAUTICS 2016年6期

    Wu Bohi,Cui Di,He Xiodong,Zhng Dinghu,*,Tng Ki

    aLaboratory of Contemporary Design and Integrated Manufacturing Technology,Ministry of Education,Northwestern Polytechnical University,Xi’an 710072,China

    bHong Kong University of Science and Technology,Hong Kong,China

    Cutting tool temperature prediction method using analytical model for end milling

    Wu Baohaia,Cui Dia,He Xiaodonga,Zhang Dinghuaa,*,Tang Kaib

    aLaboratory of Contemporary Design and Integrated Manufacturing Technology,Ministry of Education,Northwestern Polytechnical University,Xi’an 710072,China

    bHong Kong University of Science and Technology,Hong Kong,China

    Dramatic tool temperature variation in end milling can cause excessive tool wear and shorten its life,especially in machining of difficult-to-machine materials.In this study,a new analytical model-based method for the prediction of cutting tool temperature in end milling is presented.The cutting cycle is divided into temperature increase and decrease phases.For the temperature increase phase,a temperature prediction model considering real friction state between the chip and tool is proposed,and the heat flux and tool-chip contact length are then obtained through finite element simulation.In the temperature decrease phase,a temperature decrease model based on the one-dimension plate heat convection is proposed.A single wire thermocouple is employed to measure the tool temperature in the conducted milling experiments.Both of the theoretical and experimental results are obtained with cutting conditions of the cutting speed ranging from 60 m/min to 100 m/min,feed per tooth from 0.12 mm/z to 0.20 mm/z,and the radial and axial depth of cut respectively being 4 mm and 0.5 mm.The comparison results show high agreement between the physical cutting experiments and the proposed cutting tool temperature prediction method.

    1.Introduction

    Cutting heat is a fundamental physical phenomenon in machining and causes high temperature in local cutting zone,which results in many serious problems,including excessive toolwear, shortened tool life,and low machining accuracy.1–3End milling,which plays a critical role in manufacturing,has been extensively employed to finishing-cut parts of planar or curved shapes,as well as in rough machining.During an end milling process,enormous heat is produced because of the severe deformation of metal in the cutting zone and intensive friction in the chip-tool interface,especially for difficult-to-cut materials.The resulting excessive temperature inevitably worsens the tool wear and shortens its life,while atthesametimedetrimentsthesurfaceintegrityand machining quality.Understandably,it is necessary to study the temperature of cutting tool in order to improve the cutting condition in end milling.

    There are many theoretical or experimental research works on interrupted cutting temperature,especially on end milling temperature.Radulescu and Kapoor4proposed an analytical model to predict the tool temperature field in the metal cutting process,which can be applied to both continuous and interrupted three-dimensional cutting.A time-dependent heat flux model was introduced to precisely represent the heating and cooling cycle in interrupted cutting.Stephen and Ali5analyzed tool temperature in interrupted cutting,both theoretically and through physical experiments.In their work,a theoretical model of semi-infinite rectangular corner heated by heat flux varying on time with different spacious distributions is employed to investigate the tool temperature distribution;the theoretical results are compared with the measurement resultsfrom infrared and tool-chip thermocouple.Lin6researched the tool-workpiece interface temperature problem in end milling by an inverse heat conduction approach,where the machine surface temperature measured by IR pyrometer is regarded as a boundary condition and an inverse finite element method is employed to estimate the tool-work interface temperature.Ueda et al.7used a two-color pyrometer to measure the temperature of tool flank face in high speed milling and investigated the effect of cutting parameters on temperature at a carbide tool flank face.The results show that the cutting speed is the most important factor in causing the temperature rise.Lazoglu and Altintas8proposed a temperature prediction model of the tool and chip in continuous machining and for time varying milling processes,based on the finite difference method.Firstly,a heat transfer model between the chip and tool rake face is employed to study the steady cutting operation,especially the orthogonal cutting.Then this model is extended to analyze a time-varying milling process considering the chip thickness varying with time.Sato et al.9,10invented an infrared radiation pyrometer with two optical fibers connected by a fiber coupler to measure the chip-tool interface temperature in end milling with brittle CBN tool.Their method is proved very practical for measuring the chip-tool temperature during chip formation.Later they used the measurement equipment to study the cyclic temperature variation beneath the rake face of tool in end milling,and compared the measurement results with the theoretical results obtained by the same theoretical model with Stephen and Ali.The comparison results show good agreement and it is validated that the temperature variation in up milling is inverse with that of down milling.Coz et al.11proposed a temperature measuring system for rotating tools that is made of a thermocouple integrated into the milling or drilling tool near the cutting edge and a wireless transmission unit and data conditioning system incorporated into the tool-holder.Jen et al.12obtained a numerical solution of nonlinear heat conduction by a volume control method to study the time-varying cutting temperature.Yang and Zhu13analyzed the cutting temperature in milling process of titanium alloy Ti6Al4V by a finite element model of helix double-edge cutting based on a new material constitutive model.Their analysis results suggest that the temperature at the rake face is higher than that on the flank and the high temperature is closer to the cutting edge.Jen et al.14proposed a temperature prediction model applied under transient conditions,which improved Stephenson’model by a fixed-point iteration process in quasi-steady energy partitioning.Feng et al.15analyzed the workpiece temperature by a heat conduction model based on the fundamental characteristics of the milling process.Specifically,the method of weight particle swarm optimization was used to evaluate the heat flux.Their experimental results on AISI1045 show that the global maximum and minimum heat are 2.856× 106W/(m2·°C)and 2.823×106W/(m2·°C)respectively;further,the interfacial heat flux was apparently divided into three non-linear stages.Cui et al.16analyzed transient average temperature in face milling.Chen et al.17investigated heat flux and temperature distribution on tool-work interface based on a three-dimension transient model of inverse heat conduction in high speed milling process.

    Notwithstanding many valuable results aforementioned,the actual friction state of the tool-chip interface and the temperature dropping phase are not considered in these studies,resulting in inaccuracy of the analytical models of the cutting temperature.To address this issue,in this paper a new analytical model-based method for predicting the cutting tool temperature in end milling is presented,with both the aforementioned missing factors considered.

    The rest of this paper is organized as follows.The analytical cutting temperature model is presented in Section 2,including the models for both the temperature rising and temperature dropping phase.Section 3 discusses how to determine the heat flux and the chip-tool contact length using the finite element method,which is necessary in order to use the analytical model established in Section 2.The results of physical cutting experiments on Inconel718 milling are given in Section 4,and the comparison analysis between the theoretical and experimental results is discussed in Section 5.Finally,the conclusions are summarized in Section 6.

    2.Methodology

    Compared with turning,end milling is a discontinuous machining operation,with two totally opposite phases.In the temperature rising phase,the insert performs a cutting action and is heated by the heat source from the secondary and the tertiary deformation zone,causing its temperature gradually reach the top temperature.In the secondary deformation zone,the heat comes from the work done in the chip deformation and the slide friction between the insert and chip.Besides,in the tertiary deformation zone,the heat is produced from overcoming the friction between the insert flank face and the newly machined surface.Nevertheless,after accomplishing the cutting action,the insert becomes completely exposed in the ambient air,which causes its temperature to drop until the next cutting action begins,thus comprising the temperature dropping phase.These two phases make up one complete cycle of an end milling operation.

    2.1.Model of temperature rising phase

    During the temperature rising phase,the temperature of insert varies with time and it can be described as a non-steady heat conduction process.Generally,in end milling,the insert on the cutter is roughly of rectangular shape.Moreover,the insert is heated at a corner by the heat source produced by contact with workpiece.The heat source,which causes the temperature increment of insert,is produced by friction between chip and rake face.Because of a roughly rectangular friction contact area between chip and rake face,its shape can approximately be regarded as a rectangular.Before the insert performs cutting action,the conductor consists of insert and cutter has steady temperature.With the cutting action perform,it begins to be heated by heat source over insert surface.According to heat transfer,the conductor like this can be assumed as a semi-infinite body in analysis temperature change.Therefore,as an approximation,the insert is modeled as a semi-infinite rectangular corner heated by a time-varying heat flux on the tool-chip interface,as shown in Fig.1.

    Fig.1 End milling diagram and approximation heating model.

    In this paper,with the heat radiation and temperature effect on the thermal properties of the tool material neglected,the governing equation and boundary condition for the temperature T of the insert are set as where α is the thermal diffusivity,t is the time,λ the thermal conductivity,Q the heat flux,and Lxand Lyare the dimensions of the heat source.Except the interface between the tool and chip,the rest surface of tool is regarded as insulated and the initial temperature is set to be 0°C.

    The green function can be used to solve the above heat conduction equation.5The green function G for the temperature of a semi-infinite corner excited by an instantaneous heat source at time τ at the surface point(x=xp,y=yp,z=0)is G(x,y,z,xp,yp,0,D)

    where D is a parameter with units of length,and its physical meaning is characteristic dimension for the penetration of temperature field at time t resulted from instantaneous heat source at time τ.According to Osakis’theory,18the solution of the temperature in a semi-infinite corner is

    Fig.2 Stress distribution on rake face.

    The magnitude of heat flux is strongly affected by the friction along the tool-chip interface.However,there are two different friction states in this interface in the process of end milling,respectively called stick zone and slip zone,as shown in Fig.2.

    At the beginning of chip formation,the workpiece is squeezed from the tool and in turn forms a chip.As a result,the shear stress ε in the area that the chip closely contacts with the tool tip is constant kchip.Moreover,there is very high normal stress σ in this area.The chip is able to separate from the workpiece on the condition that the shear stress is increasingly promoted.At this moment,the shear stress matches the normal stress by the friction coefficient μ.In this study,the heat flux Q is considered to come from the uniform heat source and the liner heat source based on the following empirical equations:

    where β is length ratio of stick zone over slip zone,q1(x)reflects the type of heat source distribution,and q(τ)is the source heat flux varying with time.

    Eqs.(5)and(6)then will be utilized to solve the temperature equation.First,the inner double integral of the green function GR needs to be calculated:

    An error function erf(x)is introduced in the inner double integral calculation:

    Then GU which represents the integral of variable x is calculated as

    Similarly,the integral of variable y can be obtained.So the inner double integral of the green function GR is

    At last,the temperature solution of the insert on the cutter is

    2.2.Model of temperature dropping phase

    This phase refers to the time when the insert is entirely exposed in ambient air,causing its temperature to drop.The process of temperature dropping can be modeled as one of heat convection.Therefore,by approximating the insert as a thin plate with certain thickness(see Fig.3),based on the heat conduction theory,the governing equation and boundary condition for solving its temperature T can be obtained by analysis of one-dimensional thermal conduction of plate.Specifically,the following PDEs are in the following order:

    where h is the convection transfer rate and Tffluid temperature.Excess temperature θ =T(L,τ)-Tfis introduced,L is the thickness of insert,the method of separation of variables then be employed to solve the problem,i.e.

    Fig.3 Insert heat convection diagram.

    3.Heat flux and tool-chip contact length evaluation

    Apparently,the heat flux Q and the chip-tool contact length Lxmust be obtained in order to get the tool temperature.But to measure them by experiment is too difficult to realize;as an effective alternative,the finite element simulation is fixed,while the workpiece approaches at a particularly constant speed.

    3.1.Finite element model

    Deformation software is employed in the simulation of the milling Inconel718 process.In the simulation,the tool is regarded as a rigid body,while the workpiece is regarded as isotropic elastic–plastic with isotropic strain-hardening.Because of high strain rate and high temperature deformation occurring in machining Inconel718,the Jason-Cook constitutive equation,which is used by many researchers,is used to describe the high strain and temperature deformation behavior,as given in Eq.(17):

    where Td,Tr,and Tmare the deformation temperature,room temperature and melt point temperature,respectively. ε,ˉε and ˙ε are shear stress,shear strain and shear strain rate,respectively.The constants A and B stand for the material yield stress and strain hardening.The mechanics characteristics of the material are represented by the stress rate sensitive coefficient C,the strain hardening exponent n and the thermal softening coefficient m.These parameters Inconel718 are listed in Table 1.19

    The mesh function in the software is based on the Lagrangian-Euler techniques and adaptive re-meshing.Specifically,the mesh could be remeshed in the vicinity to form a normal mesh again after the initial mesh has become distorted.Furthermore,the mesh density in the cutting zone should be higher than that in the rest regions of the tool and workpiece.In the present study,the numbers of quadrangle elements in the insert and workpiece are 28197 and 49079 respectively.

    Table 1 Jason-cook constitutive of Inconel718 parameters.

    The simulation is performed at ambient temperature provided that the initial temperature of the tool and workpiece is 20°C.The convection coefficient between the workpiece and air is set as 20 W/(m2·°C).According to the Coulomb law,the friction factor is taken to be 0.4.During the simulation,the insert rake angle 5°and clearance angle 5°are fixed,while the workpiece approaches at a particularly constant speed.

    3.2.Chip-tool contact length Lxevaluation

    The whole simulation is divided into three groups based on the cutting speed,i.e.,the cutting speed V is 60 m/min,80 m/min,and 100 m/min.At each of the three,the simulation is performed by altering feed per tooth fzin the range from 0.12 mm/z to 0.20 mm/z.The radial depth of cut aeand axial depth of cut apare set to be 4 mm and 0.5 mm respectively.The fitting formula about the relationship between Lxand fzat any particular cutting speed is obtained by the leastsquare curve fitting method.The specific results from the simulation and the fitting formulas are listed in Tables 2 and 3.

    3.3.Heat flux evaluation

    The same finite element method is used in the evaluation of heat flux Q.However,compared with the Lxevaluation,the coefficient of heat flux distribution γ must be taken into account in the Q evaluation,since not all heat flux generatedin the chip-tool interface can reach the tool.20Equation below provides a way to calculate it:

    Table 2 Lxand Q obtained from simulation.

    Table 3 Fitting formula of chip-tool contact length.

    Table 4 Fitting formula of heat flux.

    where λ,ρ and c respectively represent thermal conductivity,density,and specific heat,with the subscript t indicating whether the particular coefficient is for the tool(with t)or the workpiece(without t).The specific results of the heat flux Q and the fitting formulas are listed in Table 4.

    4.Experiments

    To verify the established temperature model,a physical milling experiment is carried out.It is performed at a four-axis machine center,as shown in Fig.4.The workpiece material is Inconel718,which is extensively used in aeroengines.A double-tooth end milling cutter with a carbide insert is selected in this experiment and its specific parameters are shown in Table 5.

    Fig.4 Temperature measuring experiment.

    Table 5 Tool parameters.

    Fig.5 Single wire thermocouple and tool.

    Table 6 Critical cutting parameters.

    Table 7 Theoretical and experimental results.

    Single wire thermocouple is used in the milling experiment to measure the temperature.First,the workpiece is divided into two parts across the line of milling direction.Next,two pieces of megahit,as insulted conductor,are introduced between the two parts.Besides,a nickel–chromium wire is employed as thermocouple and located in between them.Finally,the three components mentioned above are combined as a whole sample by binder,as shown in Fig.5.In the process of milling experiment,the insulted conductor is destroyed by the milling action,which makes the nickel–chromium wire exposed and contact with the workpiece.When the tool is approaching,the temperature signal is produced.

    Fig.6 Flowchart for temperature rise calculation.

    Those cutting parameters,such as cutting speed,axial depth of cut,radial depth of cut,and feed per tooth,do affect the temperature.However,the cutting speed,and feed per tooth are more crucial and are considered to be dominating ones in insert wear.So the focus of our milling experiment is placed on these two parameters and their ranges are chosen in accordance with the optimal values from the SANDVIK recommendation,as listed in Table 6.The rest of parameters are simply set to be some constants,e.g.,the axial depth of cut is set to be 0.5 mm,and the radial depth of cut is set to be 4 mm.

    5.Results and discussion

    Table 7 shows the theoretical results calculated by our tool prediction model as the process shown in Fig.6 and the experimental results measured by single wire thermocouple under different cutting conditions with a fixed radial depth of cut and axial depth of cut.

    The trends of both theoretical and experimental results in temperature change are verified to be in good agreement.However,the experimental results tend to be lower than the theoretical calculated ones.The reason behind this deviation may be conjectured from the following aspects.First,the coated layer which can prevent heat from entering the tool has not been taken into account in the theoretical model.Secondly,the cutting speed in our experiment belongs to high speed milling for Inconel718.So,compared to conventional milling,there may not be enough heat conducting into the tool during the short contact time.Thirdly,the analysis assumption that the tool is regarded as a rectangular corner ignores the toolrakeangle,which may howeveraffectthetool temperature.

    The temperature change has an inverse trend against the cutting speed.This is partly because that the increase of cutting speed may reduce the heating time of the cutter during the machining process.The temperature is also found to increase with the per tooth feed;this agrees with the common knowledge that,with the increase of the per tooth feed,there will be more material to be removed in a cutting cycle and hence more energy is consumed and converted into heat.

    6.Conclusions

    (1)A new analytical model for predicting the cutting tool temperature in end milling is presented in this paper with consideration of real friction state on the tool-chip interface and the temperature dropping phase.

    (2)With different cutting parameters,both computer simulation using the proposed analytical model and physical cutting experiment are carried out,and the two results show good agreement.

    (3)The experimental results suggest that the tool temperature increases with the feed per tooth but decreases with the cutting speed.The proposed theoretical model for cutting tool temperature prediction can be further enhanced and then used to optimize the cutting condition to prevent excessive tool wear in end milling.

    Acknowledgments

    This work was supported by the National Basic Research Program of China(No.2013CB035802)and National Natural Science Foundation of China(No.51475382).

    1.Abukhshim NA,Mativenga PT,Sheikh MA.Heat generation and temperature prediction in metal cutting:a review and implication for high speed machining.Int J Mach Tools Manuf 2006;46(7–8):782–800.

    2.Young HT.Cutting temperature responses to flank wear.Wear 1996;201(1–2):117–20.

    3.Pradip M,Jayaramachandran R,Ganesan S.Finite element analysis of temperature rise in metal cutting process.Appl Therm Eng 2005;25(14–5):2152–68.

    4.Radulescu R,Kapoor SG.An analytical model for prediction of tool temperature fields during continuous and interrupted cutting.J Eng Ind 1994;116(2):135–43.

    5.Stephenson DA,Ali A.Tool temperatures in interrupted metal cutting.J Eng Ind 1992;114(2):127–36.

    6.Lin J.Inverse estimation of the tool-work interface temperature in end milling.Int J Mach Tools Manuf 1995;35(5):751–60.

    7.Ueda T,Hosokawa A,Oda K,Yamada K.Temperature on flank face of cutting tool in high speed milling.CIRP Ann–Manuf Technol 2001;50(1):37–40.

    8.Lazoglu I,Altintas Y.Prediction of tool and chip temperature in continuous and interrupted machining.Int J Mach Tools Manuf 2002;42(9):1011–22.

    9.Sato M,Ueda T,Tanaka H.An experimental technique for the measurement of temperature CBN tool face in end milling.Int J Mach Tools Manuf 2007;47(14):2071–6.

    10.Sato M,Tamura N,Tanaka H.Temperature variation in the cutting tool in end milling.J Manuf Sci Eng 2011;133(2):1–7.

    11.Coz GL,Marinescu M,Devillez A,Dudzinski D,Velnom L.Measuring temperature of rotating cutting tools:application to MQL drilling and dry milling of aerospace alloys.Appl Therm Eng 2012;36(4):434–41.

    12.Jen TC,Eapen S,Gutierrez G.Nonlinear numerical analysis in transient cutting tool temperatures.J Manuf Sci Eng 2003;125(1):48–56.

    13.Yang Y,Zhu WW.Study on cutting temperature during milling of titanium alloy based on FEM and experiment.Int J Adv Manuf Technol 2014;73(9):1511–21.

    14.Jen TC,Aloysius U,Anagonye AU.An improved transient model of tool temperatures in metal cutting.J Manuf Sci Eng 2001;123(1):30–7.

    15.Feng Y,Zheng L,Wang ML,Wang BS,Hou J,Yuan TJ.Research on cutting temperature of work-piece in milling process based on WPSO.Int J Adv Manuf Technol 2015;79(1):427–35.

    16.Cui XB,Zhao J,Pei ZQ.Analysis of transient average tool temperatures in face milling.Int Commun Heat Mass Transfer 2012;39(6):786–91.

    17.Chen M,Sun FH,Wang HL,Yuan RW,Qu ZH,Zhang SQ.Experimental research on the dynamic characteristics of the cutting temperature in the process of high-speed milling.J Mater Process Technol 2003;138(1–3):468–71.

    18.Ozisik MN.Heat conduction.3rd ed.New York:Wiley;1980.p.300–44.

    19.Uhlmann E,Graf von der Schulenburg M,Zettier R.Finite element modeling and cutting simulation of Inconel718.CIRP Ann– Manuf Technol 2007;56(1):61–4.

    20.Sekhon GS,Chenot JL.Numerical simulation of continuous chip formation during non-steady orthogonal cutting.Eng Comput 1993;10(1):31–48.

    Wu Baohai received the B.S.,M.S.and Ph.D.degrees from Xi’an Jiaotong University in 1997,2000,and 2005,respectively,and is now an associate professor there.His main research interests are multi-axis machining and intelligent machining.

    Cui Di received the B.S.and M.S.degrees from Xi’an Technological University,Northwestern Polytechnical University in 2011 and 2015,respectively,and is now a Ph.D.candidate there.His main research interests are multi-axis machining and intelligent machining.

    16 November 2015;revised 4 January 2016;accepted 11 January 2016

    Available online 21 October 2016

    Analytical model;

    Cutting tool;

    End milling;

    Temperature prediction;Tool temperature

    ?2016 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.This is anopenaccessarticleundertheCCBY-NC-NDlicense(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    *Corresponding author.

    E-mail address:dhzhang@nwpu.edu.cn(Z.Dinghua).

    Peer review under responsibility of Editorial Committee of CJA.

    Production and hosting by Elsevier

    http://dx.doi.org/10.1016/j.cja.2016.03.011

    1000-9361?2016 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.

    This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    天堂8中文在线网| 欧美三级亚洲精品| 中文资源天堂在线| 日韩成人伦理影院| 国产黄片美女视频| 久久国产精品大桥未久av | 日本av免费视频播放| 香蕉精品网在线| 国产精品嫩草影院av在线观看| 久久韩国三级中文字幕| 插阴视频在线观看视频| 亚洲四区av| 欧美xxxx性猛交bbbb| 少妇人妻精品综合一区二区| 国产精品一二三区在线看| 国产精品久久久久久精品古装| 九九爱精品视频在线观看| av在线老鸭窝| 日韩精品免费视频一区二区三区 | 美女国产视频在线观看| 国产免费一区二区三区四区乱码| 涩涩av久久男人的天堂| 欧美xxⅹ黑人| 中文精品一卡2卡3卡4更新| 美女内射精品一级片tv| 亚洲精品一区蜜桃| 一级毛片黄色毛片免费观看视频| 亚洲av.av天堂| 99久久综合免费| 纯流量卡能插随身wifi吗| 午夜福利,免费看| 亚洲精品久久久久久婷婷小说| 一边亲一边摸免费视频| 久久人人爽人人片av| 两个人免费观看高清视频 | 国产视频内射| 内射极品少妇av片p| 免费高清在线观看视频在线观看| 成人影院久久| 成人二区视频| 性高湖久久久久久久久免费观看| 亚洲一级一片aⅴ在线观看| 成人无遮挡网站| 免费黄网站久久成人精品| 看免费成人av毛片| 极品人妻少妇av视频| 嫩草影院入口| 亚洲国产精品国产精品| 亚洲欧美精品自产自拍| 日韩中字成人| 狂野欧美激情性bbbbbb| 午夜av观看不卡| 亚洲真实伦在线观看| 久久精品久久久久久噜噜老黄| 国产亚洲欧美精品永久| 国产成人午夜福利电影在线观看| 一个人看视频在线观看www免费| 纵有疾风起免费观看全集完整版| 精品熟女少妇av免费看| 街头女战士在线观看网站| 丰满人妻一区二区三区视频av| 国产男人的电影天堂91| 一本色道久久久久久精品综合| 国产免费福利视频在线观看| 永久免费av网站大全| 视频区图区小说| 欧美国产精品一级二级三级 | 在线观看一区二区三区激情| 国产精品一区二区三区四区免费观看| 亚洲电影在线观看av| 久久久久久久久久久久大奶| 亚洲一区二区三区欧美精品| 少妇人妻 视频| 中文精品一卡2卡3卡4更新| 久久青草综合色| 2021少妇久久久久久久久久久| 国产视频内射| 欧美精品高潮呻吟av久久| 欧美老熟妇乱子伦牲交| 久久女婷五月综合色啪小说| 欧美精品一区二区大全| 搡女人真爽免费视频火全软件| 成人漫画全彩无遮挡| 国产亚洲一区二区精品| 国产白丝娇喘喷水9色精品| 人人妻人人澡人人看| 国产成人a∨麻豆精品| 亚洲美女搞黄在线观看| 午夜福利网站1000一区二区三区| 成人漫画全彩无遮挡| 香蕉精品网在线| 亚洲经典国产精华液单| 九色成人免费人妻av| 国产成人aa在线观看| 午夜免费男女啪啪视频观看| 视频区图区小说| 亚洲va在线va天堂va国产| 欧美精品高潮呻吟av久久| 国产女主播在线喷水免费视频网站| 我的女老师完整版在线观看| 成人影院久久| 国产精品女同一区二区软件| 欧美xxⅹ黑人| 日产精品乱码卡一卡2卡三| 一级,二级,三级黄色视频| 在线观看一区二区三区激情| 欧美一级a爱片免费观看看| 国产av一区二区精品久久| 久久久久久久久久久免费av| 国产女主播在线喷水免费视频网站| 又大又黄又爽视频免费| 国产精品免费大片| 国产精品久久久久久av不卡| 秋霞在线观看毛片| 丁香六月天网| 欧美人与善性xxx| 国产成人精品久久久久久| 日韩不卡一区二区三区视频在线| 青春草亚洲视频在线观看| 欧美精品亚洲一区二区| 精品亚洲成国产av| 内射极品少妇av片p| 国产欧美亚洲国产| 亚洲精品中文字幕在线视频 | 女的被弄到高潮叫床怎么办| 赤兔流量卡办理| 免费观看av网站的网址| 一区在线观看完整版| 久久久精品94久久精品| 能在线免费看毛片的网站| 国产一级毛片在线| 国产在线男女| 永久网站在线| 久久久久精品久久久久真实原创| 国产精品一区二区在线不卡| 狂野欧美激情性bbbbbb| 成人亚洲精品一区在线观看| 国产亚洲午夜精品一区二区久久| 久久久久久久大尺度免费视频| .国产精品久久| 自拍欧美九色日韩亚洲蝌蚪91 | 国产一区二区在线观看日韩| 国产片特级美女逼逼视频| 亚洲精华国产精华液的使用体验| 亚洲国产精品专区欧美| 99热6这里只有精品| 国产在线视频一区二区| av又黄又爽大尺度在线免费看| 水蜜桃什么品种好| 夫妻性生交免费视频一级片| 777米奇影视久久| 97精品久久久久久久久久精品| 人人妻人人爽人人添夜夜欢视频 | 热re99久久精品国产66热6| 精品酒店卫生间| 成人国产av品久久久| 亚洲精品乱码久久久v下载方式| 精品熟女少妇av免费看| 午夜免费鲁丝| 美女中出高潮动态图| 免费高清在线观看视频在线观看| 最近手机中文字幕大全| 亚洲,欧美,日韩| 国产一区二区在线观看av| 国产淫片久久久久久久久| 日韩一区二区视频免费看| 哪个播放器可以免费观看大片| av在线观看视频网站免费| 日韩成人av中文字幕在线观看| 欧美 亚洲 国产 日韩一| 国产一区有黄有色的免费视频| 国产av码专区亚洲av| 婷婷色综合大香蕉| 美女cb高潮喷水在线观看| 久久毛片免费看一区二区三区| 精品久久久久久久久亚洲| 婷婷色综合www| 两个人免费观看高清视频 | 又黄又爽又刺激的免费视频.| 日韩av在线免费看完整版不卡| 高清不卡的av网站| 精品久久久噜噜| 国产亚洲5aaaaa淫片| 嫩草影院入口| 18+在线观看网站| 亚洲美女视频黄频| 天堂俺去俺来也www色官网| 欧美最新免费一区二区三区| 人人澡人人妻人| 亚洲不卡免费看| 高清黄色对白视频在线免费看 | 99精国产麻豆久久婷婷| 人人澡人人妻人| 在线观看人妻少妇| 少妇被粗大猛烈的视频| 亚洲精品第二区| 色94色欧美一区二区| 男女边吃奶边做爰视频| 91成人精品电影| 国产亚洲5aaaaa淫片| 亚洲va在线va天堂va国产| 中文资源天堂在线| 2018国产大陆天天弄谢| 国产91av在线免费观看| av在线观看视频网站免费| 最后的刺客免费高清国语| 少妇的逼水好多| 一区二区av电影网| 欧美3d第一页| 久久免费观看电影| 美女cb高潮喷水在线观看| 午夜福利,免费看| 各种免费的搞黄视频| 国产69精品久久久久777片| 99精国产麻豆久久婷婷| 国精品久久久久久国模美| 国产成人精品一,二区| 国产真实伦视频高清在线观看| 亚洲av免费高清在线观看| 免费看光身美女| 黑人高潮一二区| 日本黄色片子视频| 亚洲电影在线观看av| 校园人妻丝袜中文字幕| 少妇精品久久久久久久| 久久久久精品久久久久真实原创| 国产亚洲一区二区精品| 不卡视频在线观看欧美| 亚洲欧美一区二区三区黑人 | 国产色婷婷99| 国产有黄有色有爽视频| 日韩中文字幕视频在线看片| 我要看黄色一级片免费的| 丝袜喷水一区| 天堂中文最新版在线下载| 欧美精品国产亚洲| 国产 精品1| 亚洲精品一二三| 在线观看免费日韩欧美大片 | 我要看黄色一级片免费的| 欧美97在线视频| 日韩av不卡免费在线播放| 亚洲情色 制服丝袜| 晚上一个人看的免费电影| 又黄又爽又刺激的免费视频.| 国产片特级美女逼逼视频| 精品亚洲成a人片在线观看| 国产在线视频一区二区| 日韩中字成人| 国产免费一区二区三区四区乱码| 新久久久久国产一级毛片| 高清午夜精品一区二区三区| 久久久久久久久大av| 亚洲av免费高清在线观看| 秋霞伦理黄片| 日日啪夜夜撸| 亚洲精品久久午夜乱码| 亚洲美女搞黄在线观看| 一级毛片黄色毛片免费观看视频| 少妇的逼好多水| 麻豆精品久久久久久蜜桃| 日本黄色日本黄色录像| 视频区图区小说| 99久久精品国产国产毛片| 国产成人freesex在线| 亚洲内射少妇av| 亚洲国产精品国产精品| 少妇人妻 视频| 啦啦啦啦在线视频资源| 亚洲真实伦在线观看| 日韩亚洲欧美综合| 九九在线视频观看精品| 大片免费播放器 马上看| 亚洲性久久影院| 九色成人免费人妻av| 新久久久久国产一级毛片| 国产熟女午夜一区二区三区 | 精品亚洲乱码少妇综合久久| 国产一区二区在线观看日韩| 国产精品国产三级国产专区5o| 亚洲国产精品一区二区三区在线| 国产精品久久久久久av不卡| 精品酒店卫生间| 欧美性感艳星| 国产精品不卡视频一区二区| 黑人猛操日本美女一级片| 少妇精品久久久久久久| 日韩大片免费观看网站| 国产爽快片一区二区三区| 亚洲欧美精品专区久久| 国产美女午夜福利| 一级毛片黄色毛片免费观看视频| 午夜91福利影院| 亚洲欧美日韩东京热| 久久ye,这里只有精品| 大又大粗又爽又黄少妇毛片口| 国产成人91sexporn| 纯流量卡能插随身wifi吗| av在线观看视频网站免费| 国产 精品1| 永久网站在线| 日本欧美国产在线视频| 午夜福利在线观看免费完整高清在| 大片免费播放器 马上看| 免费看光身美女| 成人18禁高潮啪啪吃奶动态图 | 久久99热这里只频精品6学生| h视频一区二区三区| 亚洲av男天堂| 涩涩av久久男人的天堂| av免费在线看不卡| 精品亚洲乱码少妇综合久久| 少妇人妻久久综合中文| av黄色大香蕉| 日日撸夜夜添| av线在线观看网站| 国产一区二区在线观看av| 18禁动态无遮挡网站| 在线观看免费日韩欧美大片 | 久久av网站| 国产亚洲午夜精品一区二区久久| 日韩一区二区视频免费看| 国产精品人妻久久久影院| 亚洲av在线观看美女高潮| 日本午夜av视频| 韩国高清视频一区二区三区| 免费人妻精品一区二区三区视频| 中文乱码字字幕精品一区二区三区| 91精品一卡2卡3卡4卡| 国产亚洲91精品色在线| 九九爱精品视频在线观看| 亚洲av日韩在线播放| 黄色毛片三级朝国网站 | 日本欧美国产在线视频| 久久人妻熟女aⅴ| 中文字幕av电影在线播放| 99国产精品免费福利视频| 亚洲精品国产色婷婷电影| 亚洲丝袜综合中文字幕| 国产亚洲91精品色在线| 国产欧美日韩精品一区二区| 欧美另类一区| 国产欧美日韩精品一区二区| 97在线视频观看| 中文字幕久久专区| 国产精品免费大片| 午夜免费男女啪啪视频观看| 亚洲国产精品一区三区| 亚洲国产精品999| 伦理电影免费视频| 久久亚洲国产成人精品v| 噜噜噜噜噜久久久久久91| 午夜免费观看性视频| 九色成人免费人妻av| 男人添女人高潮全过程视频| 一级a做视频免费观看| 亚洲三级黄色毛片| 欧美另类一区| 免费看不卡的av| 亚洲精品国产成人久久av| 欧美变态另类bdsm刘玥| 亚洲精品国产色婷婷电影| 欧美变态另类bdsm刘玥| 一级二级三级毛片免费看| 欧美老熟妇乱子伦牲交| 天美传媒精品一区二区| 伦理电影大哥的女人| 啦啦啦中文免费视频观看日本| 一个人看视频在线观看www免费| 在线观看美女被高潮喷水网站| 高清视频免费观看一区二区| 亚洲色图综合在线观看| 九草在线视频观看| 好男人视频免费观看在线| 亚洲av二区三区四区| 六月丁香七月| 全区人妻精品视频| 国产淫片久久久久久久久| 男的添女的下面高潮视频| 在线 av 中文字幕| 两个人的视频大全免费| 99热6这里只有精品| 国产亚洲5aaaaa淫片| 丰满饥渴人妻一区二区三| 亚洲美女黄色视频免费看| 日本黄色片子视频| 26uuu在线亚洲综合色| 在线亚洲精品国产二区图片欧美 | 国产成人a∨麻豆精品| h日本视频在线播放| 国产成人免费无遮挡视频| a级片在线免费高清观看视频| www.色视频.com| 啦啦啦视频在线资源免费观看| 欧美日韩国产mv在线观看视频| 在线观看免费视频网站a站| 欧美bdsm另类| 日韩亚洲欧美综合| 伊人久久精品亚洲午夜| 国产亚洲91精品色在线| 亚洲精品中文字幕在线视频 | 成人国产麻豆网| 韩国高清视频一区二区三区| 人妻夜夜爽99麻豆av| 亚洲婷婷狠狠爱综合网| 欧美人与善性xxx| a级毛片免费高清观看在线播放| 日日啪夜夜爽| 深夜a级毛片| 国产黄色视频一区二区在线观看| 国产精品一区二区在线观看99| 伦理电影免费视频| 成人综合一区亚洲| 色婷婷av一区二区三区视频| 美女视频免费永久观看网站| 国产永久视频网站| 插阴视频在线观看视频| 夜夜骑夜夜射夜夜干| 观看免费一级毛片| 3wmmmm亚洲av在线观看| 国产精品熟女久久久久浪| 国产毛片在线视频| 国产亚洲欧美精品永久| 涩涩av久久男人的天堂| 精品久久久久久久久亚洲| 99热全是精品| 中文天堂在线官网| 午夜91福利影院| 国产精品一区二区三区四区免费观看| 搡老乐熟女国产| 免费人成在线观看视频色| 99热网站在线观看| 嘟嘟电影网在线观看| 亚洲精品久久久久久婷婷小说| 国产淫片久久久久久久久| 你懂的网址亚洲精品在线观看| 精品少妇内射三级| 国产在线视频一区二区| 99九九线精品视频在线观看视频| 日韩电影二区| 极品人妻少妇av视频| 精品人妻熟女av久视频| 纯流量卡能插随身wifi吗| 天天操日日干夜夜撸| 亚洲精品日本国产第一区| 国产有黄有色有爽视频| 最近最新中文字幕免费大全7| 少妇的逼好多水| av在线老鸭窝| 久久久久久久久大av| 中文字幕制服av| 亚洲美女黄色视频免费看| 99热6这里只有精品| 大话2 男鬼变身卡| 国产淫片久久久久久久久| 22中文网久久字幕| 免费在线观看成人毛片| 天天操日日干夜夜撸| 大香蕉97超碰在线| 黄色配什么色好看| 亚洲国产精品一区三区| 80岁老熟妇乱子伦牲交| 亚州av有码| 综合色丁香网| videos熟女内射| 最黄视频免费看| 国产精品伦人一区二区| freevideosex欧美| 亚洲国产色片| 亚洲经典国产精华液单| 黑人猛操日本美女一级片| 亚洲图色成人| 狠狠精品人妻久久久久久综合| 蜜臀久久99精品久久宅男| 99热全是精品| 日日啪夜夜撸| a 毛片基地| 日韩av不卡免费在线播放| 日韩 亚洲 欧美在线| 午夜日本视频在线| 国产视频内射| 欧美日韩视频精品一区| 国产精品久久久久久久久免| 一区在线观看完整版| 亚洲精品国产av蜜桃| 国产伦精品一区二区三区视频9| 精品亚洲成a人片在线观看| 亚洲精品色激情综合| 一个人免费看片子| 777米奇影视久久| 久久婷婷青草| a级毛片免费高清观看在线播放| h视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 欧美 日韩 精品 国产| 国产精品福利在线免费观看| 午夜激情福利司机影院| 国产成人精品久久久久久| 久久6这里有精品| 80岁老熟妇乱子伦牲交| 亚洲国产精品专区欧美| 香蕉精品网在线| 深夜a级毛片| 最近的中文字幕免费完整| 国产精品一区二区性色av| 春色校园在线视频观看| 欧美日韩综合久久久久久| 亚洲欧美日韩卡通动漫| 国产在线免费精品| 少妇猛男粗大的猛烈进出视频| 五月开心婷婷网| 99久久综合免费| 久久国产精品大桥未久av | 久久毛片免费看一区二区三区| 男女国产视频网站| 亚洲精华国产精华液的使用体验| 久久久久人妻精品一区果冻| 国产伦理片在线播放av一区| 观看美女的网站| 亚洲国产日韩一区二区| 亚洲,一卡二卡三卡| 性高湖久久久久久久久免费观看| 在线观看免费高清a一片| 夫妻性生交免费视频一级片| 精品久久久精品久久久| 国产熟女欧美一区二区| 2021少妇久久久久久久久久久| 久久99一区二区三区| 黄色视频在线播放观看不卡| 美女大奶头黄色视频| 狂野欧美激情性bbbbbb| 国产精品秋霞免费鲁丝片| 国产在线男女| 99九九在线精品视频 | 久久久久久人妻| 国产精品三级大全| 国产黄片美女视频| 亚洲精品成人av观看孕妇| www.色视频.com| 在线播放无遮挡| 亚洲精品国产av成人精品| 国产精品福利在线免费观看| 亚洲av成人精品一区久久| 校园人妻丝袜中文字幕| 两个人的视频大全免费| 99热这里只有精品一区| 一级毛片电影观看| 成人亚洲欧美一区二区av| 青春草国产在线视频| 18禁裸乳无遮挡动漫免费视频| 亚洲,欧美,日韩| 热99国产精品久久久久久7| 日本午夜av视频| 精品一品国产午夜福利视频| 91精品国产国语对白视频| 男女边摸边吃奶| 亚洲精品一区蜜桃| 一级二级三级毛片免费看| 欧美三级亚洲精品| 晚上一个人看的免费电影| 一级a做视频免费观看| 人人妻人人澡人人爽人人夜夜| 99热6这里只有精品| 成人国产av品久久久| 亚洲av综合色区一区| 久热这里只有精品99| 黑人猛操日本美女一级片| 大又大粗又爽又黄少妇毛片口| √禁漫天堂资源中文www| av专区在线播放| 成人影院久久| 伊人久久精品亚洲午夜| 大片免费播放器 马上看| 日本-黄色视频高清免费观看| 曰老女人黄片| 男的添女的下面高潮视频| 在线观看www视频免费| 亚洲欧洲国产日韩| 丝袜喷水一区| 黑人猛操日本美女一级片| 久久久久久久久久成人| 国产精品久久久久久av不卡| 亚洲电影在线观看av| 最近的中文字幕免费完整| 哪个播放器可以免费观看大片| 久久久久久久久久久丰满| 午夜激情久久久久久久| 日日摸夜夜添夜夜添av毛片| 一区二区三区四区激情视频| 99精国产麻豆久久婷婷| 欧美精品亚洲一区二区| 国产一区二区三区av在线| 亚洲精品一区蜜桃| 国产精品免费大片| 欧美激情国产日韩精品一区| 好男人视频免费观看在线| 两个人的视频大全免费| 国产日韩欧美视频二区| 精品久久久久久久久亚洲| 国产淫片久久久久久久久| 两个人免费观看高清视频 | 久久久久久人妻| av免费在线看不卡| 精品一品国产午夜福利视频| 成人漫画全彩无遮挡| 亚洲经典国产精华液单| 五月天丁香电影| 国产午夜精品久久久久久一区二区三区| 欧美成人午夜免费资源| 又粗又硬又长又爽又黄的视频| 亚洲无线观看免费| 久久青草综合色| 欧美三级亚洲精品| 国产视频内射| 老司机影院毛片|