• <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/).

    91字幕亚洲| 又黄又粗又硬又大视频| 91麻豆av在线| 真人一进一出gif抽搐免费| 国产av精品麻豆| 一区福利在线观看| 亚洲中文字幕一区二区三区有码在线看 | 男人舔女人下体高潮全视频| 曰老女人黄片| 黄色毛片三级朝国网站| 亚洲国产中文字幕在线视频| 无人区码免费观看不卡| 这个男人来自地球电影免费观看| 国产三级黄色录像| 精品一品国产午夜福利视频| 亚洲伊人色综图| 久久精品91蜜桃| 老司机深夜福利视频在线观看| 国产成人av激情在线播放| 精品久久蜜臀av无| 美女 人体艺术 gogo| 老鸭窝网址在线观看| 亚洲人成77777在线视频| 午夜福利一区二区在线看| 91成人精品电影| 少妇粗大呻吟视频| 国产99白浆流出| 一边摸一边抽搐一进一出视频| 热re99久久国产66热| 久久国产精品人妻蜜桃| 国产精品久久电影中文字幕| 美女 人体艺术 gogo| 青草久久国产| 国产精品久久久久久人妻精品电影| 久久久久久久午夜电影| 欧美黄色片欧美黄色片| 一a级毛片在线观看| 18禁观看日本| 国产一级毛片七仙女欲春2 | 日日夜夜操网爽| 午夜成年电影在线免费观看| 精品国产一区二区久久| 男人操女人黄网站| 亚洲五月天丁香| 一夜夜www| 99在线人妻在线中文字幕| 精品国产美女av久久久久小说| 欧美日本亚洲视频在线播放| 成人av一区二区三区在线看| 一边摸一边做爽爽视频免费| 国产精品乱码一区二三区的特点 | 日本 av在线| 久久国产乱子伦精品免费另类| 悠悠久久av| 变态另类成人亚洲欧美熟女 | 女性生殖器流出的白浆| 婷婷六月久久综合丁香| 国产在线观看jvid| 无遮挡黄片免费观看| 成人国产一区最新在线观看| 妹子高潮喷水视频| 麻豆成人av在线观看| 岛国在线观看网站| 色婷婷久久久亚洲欧美| 国产伦一二天堂av在线观看| 国产99久久九九免费精品| 成人18禁高潮啪啪吃奶动态图| 午夜久久久在线观看| 99久久精品国产亚洲精品| 别揉我奶头~嗯~啊~动态视频| 亚洲五月婷婷丁香| 一级片免费观看大全| 男女之事视频高清在线观看| 亚洲狠狠婷婷综合久久图片| 国产三级黄色录像| 欧美激情极品国产一区二区三区| 亚洲成av片中文字幕在线观看| 色综合亚洲欧美另类图片| 一进一出好大好爽视频| 一进一出好大好爽视频| 亚洲精品久久成人aⅴ小说| 亚洲片人在线观看| 亚洲第一av免费看| 巨乳人妻的诱惑在线观看| av中文乱码字幕在线| 一区二区三区国产精品乱码| 国产激情欧美一区二区| 成人18禁在线播放| 中出人妻视频一区二区| 精品一区二区三区四区五区乱码| 精品国产一区二区久久| 日韩一卡2卡3卡4卡2021年| 69精品国产乱码久久久| 欧美人与性动交α欧美精品济南到| 后天国语完整版免费观看| 一本久久中文字幕| 国产精品久久视频播放| 亚洲精华国产精华精| 国产精品久久久久久人妻精品电影| 满18在线观看网站| 高清黄色对白视频在线免费看| 午夜亚洲福利在线播放| 69av精品久久久久久| 91av网站免费观看| 久久中文字幕人妻熟女| 黄色视频,在线免费观看| 大陆偷拍与自拍| 精品免费久久久久久久清纯| а√天堂www在线а√下载| 制服丝袜大香蕉在线| 大陆偷拍与自拍| 俄罗斯特黄特色一大片| 大陆偷拍与自拍| 两个人免费观看高清视频| 村上凉子中文字幕在线| 91成人精品电影| 国产精品久久久av美女十八| 精品国产一区二区三区四区第35| 国产真人三级小视频在线观看| 久久国产乱子伦精品免费另类| 男人舔女人下体高潮全视频| 91麻豆精品激情在线观看国产| 亚洲第一青青草原| 精品国产美女av久久久久小说| 亚洲av美国av| 国产精品二区激情视频| 国产精品日韩av在线免费观看 | 丰满人妻熟妇乱又伦精品不卡| 搡老熟女国产l中国老女人| 久久 成人 亚洲| 亚洲 欧美一区二区三区| 久久久久久久久中文| 欧美乱色亚洲激情| 欧美国产精品va在线观看不卡| 不卡av一区二区三区| 九色亚洲精品在线播放| 中文亚洲av片在线观看爽| 亚洲精品中文字幕在线视频| 丝袜在线中文字幕| 国产男靠女视频免费网站| 国产亚洲精品第一综合不卡| 巨乳人妻的诱惑在线观看| 精品一品国产午夜福利视频| 免费av毛片视频| 国产精品二区激情视频| 久久香蕉国产精品| 欧美日韩乱码在线| 国产高清有码在线观看视频 | 久久久国产成人精品二区| 亚洲第一欧美日韩一区二区三区| 日本在线视频免费播放| 久久久久久久久免费视频了| 国产激情欧美一区二区| 国产亚洲精品av在线| 岛国在线观看网站| 黄网站色视频无遮挡免费观看| 黄色 视频免费看| 欧美+亚洲+日韩+国产| 村上凉子中文字幕在线| 日韩 欧美 亚洲 中文字幕| 少妇被粗大的猛进出69影院| 亚洲成人国产一区在线观看| 亚洲成av片中文字幕在线观看| www.999成人在线观看| 亚洲精品粉嫩美女一区| 国产精品自产拍在线观看55亚洲| 夜夜爽天天搞| 久久性视频一级片| 日本vs欧美在线观看视频| 首页视频小说图片口味搜索| 色av中文字幕| 亚洲精品国产色婷婷电影| 日日爽夜夜爽网站| 老司机深夜福利视频在线观看| 久热爱精品视频在线9| 欧美日韩瑟瑟在线播放| 国产精品99久久99久久久不卡| 亚洲熟女毛片儿| 久久久久久人人人人人| 午夜免费成人在线视频| 男女做爰动态图高潮gif福利片 | 亚洲中文字幕一区二区三区有码在线看 | 亚洲一区高清亚洲精品| 午夜久久久在线观看| 人妻丰满熟妇av一区二区三区| 成人精品一区二区免费| 色婷婷久久久亚洲欧美| 久久久久久免费高清国产稀缺| 国产av在哪里看| 一a级毛片在线观看| 亚洲成人国产一区在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 欧美激情 高清一区二区三区| 久久久久久国产a免费观看| 国内毛片毛片毛片毛片毛片| 丝袜人妻中文字幕| 好看av亚洲va欧美ⅴa在| 免费无遮挡裸体视频| 亚洲精品国产区一区二| 欧美av亚洲av综合av国产av| 色综合亚洲欧美另类图片| 黄色视频不卡| 免费人成视频x8x8入口观看| 成年人黄色毛片网站| 性少妇av在线| 成人18禁高潮啪啪吃奶动态图| 精品国产国语对白av| 女同久久另类99精品国产91| 人人妻,人人澡人人爽秒播| av在线播放免费不卡| 国产成年人精品一区二区| 大陆偷拍与自拍| 亚洲精品久久国产高清桃花| 无人区码免费观看不卡| 亚洲精品国产一区二区精华液| 怎么达到女性高潮| 熟妇人妻久久中文字幕3abv| 国产视频一区二区在线看| 好男人电影高清在线观看| av在线播放免费不卡| 黄色视频不卡| 免费搜索国产男女视频| 日本一区二区免费在线视频| 人妻久久中文字幕网| 啦啦啦韩国在线观看视频| 亚洲国产精品合色在线| 不卡一级毛片| 日韩精品免费视频一区二区三区| 亚洲av日韩精品久久久久久密| 国产精品 国内视频| 99久久精品国产亚洲精品| 怎么达到女性高潮| 欧美精品啪啪一区二区三区| 亚洲五月色婷婷综合| 日本欧美视频一区| 国产一区二区在线av高清观看| 日韩av在线大香蕉| 欧美日韩一级在线毛片| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文日韩欧美视频| 黄网站色视频无遮挡免费观看| av天堂在线播放| 欧美黑人欧美精品刺激| 精品久久久久久久毛片微露脸| 99riav亚洲国产免费| 亚洲成国产人片在线观看| 国产成人av激情在线播放| 亚洲专区字幕在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品久久成人aⅴ小说| 国产av精品麻豆| 男人操女人黄网站| 欧美国产日韩亚洲一区| 色综合婷婷激情| 一区二区三区激情视频| 啦啦啦免费观看视频1| 久久香蕉国产精品| 给我免费播放毛片高清在线观看| 国产亚洲欧美精品永久| 亚洲精品av麻豆狂野| 又紧又爽又黄一区二区| 搡老岳熟女国产| 国产一区二区激情短视频| 午夜福利在线观看吧| 女人精品久久久久毛片| 国产99久久九九免费精品| 97人妻天天添夜夜摸| 久久国产乱子伦精品免费另类| 变态另类丝袜制服| 中亚洲国语对白在线视频| 久久这里只有精品19| 人妻丰满熟妇av一区二区三区| 欧美精品亚洲一区二区| 午夜福利18| 欧美激情极品国产一区二区三区| 久久精品91无色码中文字幕| 婷婷丁香在线五月| 村上凉子中文字幕在线| 欧美日本中文国产一区发布| 日本 欧美在线| 在线av久久热| 精品国产国语对白av| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美一区视频在线观看| 精品国产乱码久久久久久男人| 18禁黄网站禁片午夜丰满| 国产高清videossex| 老司机深夜福利视频在线观看| 午夜精品国产一区二区电影| 性欧美人与动物交配| 午夜免费鲁丝| 午夜两性在线视频| av在线播放免费不卡| 国产精品精品国产色婷婷| 中文字幕人妻熟女乱码| 多毛熟女@视频| 天天一区二区日本电影三级 | 中国美女看黄片| 黄色片一级片一级黄色片| 免费观看精品视频网站| 亚洲中文字幕日韩| 免费av毛片视频| 一个人观看的视频www高清免费观看 | 久久国产精品影院| av在线播放免费不卡| 久久人妻福利社区极品人妻图片| av福利片在线| 脱女人内裤的视频| 亚洲电影在线观看av| 亚洲视频免费观看视频| 国产高清有码在线观看视频 | 丝袜美足系列| ponron亚洲| 变态另类丝袜制服| 99香蕉大伊视频| 精品久久久久久久久久免费视频| av免费在线观看网站| 精品国产一区二区三区四区第35| 性色av乱码一区二区三区2| 免费在线观看黄色视频的| 啦啦啦观看免费观看视频高清 | av网站免费在线观看视频| 成人三级做爰电影| 欧美 亚洲 国产 日韩一| 久久久久亚洲av毛片大全| 日韩三级视频一区二区三区| 如日韩欧美国产精品一区二区三区| 午夜成年电影在线免费观看| 亚洲午夜精品一区,二区,三区| 在线观看一区二区三区| 欧美在线一区亚洲| 色综合欧美亚洲国产小说| 亚洲自偷自拍图片 自拍| 人人妻,人人澡人人爽秒播| 久久久久久久午夜电影| 国产极品粉嫩免费观看在线| 人妻丰满熟妇av一区二区三区| 妹子高潮喷水视频| 欧美黑人欧美精品刺激| 男人的好看免费观看在线视频 | 午夜a级毛片| 日韩av在线大香蕉| 亚洲avbb在线观看| 曰老女人黄片| avwww免费| 电影成人av| 国产精品永久免费网站| 国产精品免费视频内射| 男女做爰动态图高潮gif福利片 | 在线av久久热| 好看av亚洲va欧美ⅴa在| 人人妻人人澡欧美一区二区 | 脱女人内裤的视频| 亚洲中文字幕一区二区三区有码在线看 | 9热在线视频观看99| 国产av一区二区精品久久| 麻豆av在线久日| 正在播放国产对白刺激| 精品不卡国产一区二区三区| 美女高潮喷水抽搐中文字幕| 给我免费播放毛片高清在线观看| 黄色a级毛片大全视频| 少妇粗大呻吟视频| 女性生殖器流出的白浆| 精品国产一区二区久久| 9191精品国产免费久久| 午夜影院日韩av| 亚洲精品中文字幕一二三四区| 久久精品国产99精品国产亚洲性色 | 高潮久久久久久久久久久不卡| 欧美日韩黄片免| 国产精品一区二区三区四区久久 | 久久精品aⅴ一区二区三区四区| 高清黄色对白视频在线免费看| 大码成人一级视频| 欧美性长视频在线观看| 99re在线观看精品视频| 午夜福利在线观看吧| 欧美国产日韩亚洲一区| 亚洲欧美激情在线| 人人妻人人爽人人添夜夜欢视频| 给我免费播放毛片高清在线观看| 国产成人av激情在线播放| 黄色 视频免费看| 国产一区二区三区在线臀色熟女| 久久久国产成人精品二区| 美女国产高潮福利片在线看| 久久精品国产亚洲av香蕉五月| 亚洲人成电影观看| 50天的宝宝边吃奶边哭怎么回事| 欧美一区二区精品小视频在线| 成人亚洲精品一区在线观看| 久久久久国产一级毛片高清牌| 9热在线视频观看99| 国产高清激情床上av| 久久天躁狠狠躁夜夜2o2o| av天堂久久9| 老汉色∧v一级毛片| 琪琪午夜伦伦电影理论片6080| 极品人妻少妇av视频| 欧美精品啪啪一区二区三区| 久久婷婷人人爽人人干人人爱 | 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产美女av久久久久小说| a在线观看视频网站| 亚洲 欧美一区二区三区| www.自偷自拍.com| 国内毛片毛片毛片毛片毛片| 精品午夜福利视频在线观看一区| 成年版毛片免费区| 18禁国产床啪视频网站| 国产亚洲av嫩草精品影院| 一级毛片精品| 成人av一区二区三区在线看| 91成年电影在线观看| 91麻豆av在线| 欧美成人午夜精品| 国产一区二区激情短视频| 欧美中文日本在线观看视频| 手机成人av网站| 亚洲视频免费观看视频| bbb黄色大片| 国产免费男女视频| 日韩欧美一区二区三区在线观看| a在线观看视频网站| 91在线观看av| 日韩 欧美 亚洲 中文字幕| 国产亚洲av嫩草精品影院| 丝袜美足系列| av在线播放免费不卡| 十八禁人妻一区二区| 一级毛片高清免费大全| 村上凉子中文字幕在线| 亚洲五月天丁香| 中文字幕久久专区| 涩涩av久久男人的天堂| 国产成人精品久久二区二区91| 国产av一区在线观看免费| 午夜视频精品福利| 欧美激情高清一区二区三区| 午夜久久久在线观看| 身体一侧抽搐| 日日干狠狠操夜夜爽| 久久国产精品人妻蜜桃| 日韩免费av在线播放| 香蕉丝袜av| xxx96com| 人人妻,人人澡人人爽秒播| 国产精品 国内视频| 窝窝影院91人妻| 巨乳人妻的诱惑在线观看| 亚洲色图 男人天堂 中文字幕| 丝袜美足系列| 国产免费av片在线观看野外av| 国产麻豆成人av免费视频| 不卡av一区二区三区| 亚洲成人国产一区在线观看| 亚洲成a人片在线一区二区| 亚洲电影在线观看av| 嫩草影院精品99| 高清在线国产一区| 亚洲 国产 在线| 老鸭窝网址在线观看| 真人做人爱边吃奶动态| 97人妻精品一区二区三区麻豆 | 男人舔女人下体高潮全视频| 91国产中文字幕| 看黄色毛片网站| 精品国产美女av久久久久小说| 国产极品粉嫩免费观看在线| 老熟妇乱子伦视频在线观看| 日韩三级视频一区二区三区| 十分钟在线观看高清视频www| 一区二区三区激情视频| 亚洲九九香蕉| 亚洲国产精品久久男人天堂| 一进一出抽搐gif免费好疼| 国产伦人伦偷精品视频| 久久影院123| 无遮挡黄片免费观看| 日韩大码丰满熟妇| 亚洲一区二区三区不卡视频| 亚洲激情在线av| av片东京热男人的天堂| 亚洲三区欧美一区| 精品免费久久久久久久清纯| 久久国产亚洲av麻豆专区| 91麻豆精品激情在线观看国产| 精品无人区乱码1区二区| 日韩中文字幕欧美一区二区| 天天添夜夜摸| 亚洲在线自拍视频| 欧美在线黄色| 精品熟女少妇八av免费久了| 搞女人的毛片| 国产精品九九99| 精品久久久久久成人av| 自线自在国产av| 电影成人av| 韩国av一区二区三区四区| 国产91精品成人一区二区三区| 人人妻人人澡欧美一区二区 | 级片在线观看| 97人妻精品一区二区三区麻豆 | 国产aⅴ精品一区二区三区波| 午夜福利高清视频| 一级,二级,三级黄色视频| 亚洲国产日韩欧美精品在线观看 | 久99久视频精品免费| 国产精品一区二区精品视频观看| 亚洲狠狠婷婷综合久久图片| 精品第一国产精品| 精品一品国产午夜福利视频| 性欧美人与动物交配| 18禁观看日本| 久久午夜亚洲精品久久| 亚洲av五月六月丁香网| 免费看美女性在线毛片视频| 波多野结衣巨乳人妻| 深夜精品福利| √禁漫天堂资源中文www| 精品人妻1区二区| 女警被强在线播放| 亚洲电影在线观看av| 免费人成视频x8x8入口观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲av日韩精品久久久久久密| 动漫黄色视频在线观看| 久久久久久久久久久久大奶| 午夜福利影视在线免费观看| 琪琪午夜伦伦电影理论片6080| 午夜久久久久精精品| 亚洲片人在线观看| 亚洲欧美精品综合久久99| av天堂在线播放| 亚洲人成电影免费在线| 黄网站色视频无遮挡免费观看| 一个人观看的视频www高清免费观看 | 99热只有精品国产| 国产激情久久老熟女| 国产又爽黄色视频| 亚洲一区高清亚洲精品| www.自偷自拍.com| 国产欧美日韩综合在线一区二区| 亚洲国产高清在线一区二区三 | 国产精品香港三级国产av潘金莲| e午夜精品久久久久久久| 老汉色∧v一级毛片| 给我免费播放毛片高清在线观看| 国产视频一区二区在线看| 国产精品久久视频播放| 欧美不卡视频在线免费观看 | 丝袜美足系列| 女性生殖器流出的白浆| 中文字幕av电影在线播放| 日本精品一区二区三区蜜桃| 色尼玛亚洲综合影院| 日本 av在线| 熟女少妇亚洲综合色aaa.| 日韩成人在线观看一区二区三区| 亚洲欧美激情综合另类| 一进一出抽搐gif免费好疼| 久久午夜亚洲精品久久| 精品国产超薄肉色丝袜足j| 日韩大码丰满熟妇| 麻豆成人av在线观看| 国产精品 国内视频| 亚洲va日本ⅴa欧美va伊人久久| 久久中文字幕一级| 美国免费a级毛片| 99国产综合亚洲精品| 亚洲国产日韩欧美精品在线观看 | 国产精品一区二区免费欧美| 动漫黄色视频在线观看| 男女之事视频高清在线观看| 精品久久久精品久久久| 久久久精品欧美日韩精品| 啦啦啦韩国在线观看视频| 欧美精品亚洲一区二区| 国产成人精品无人区| 少妇熟女aⅴ在线视频| 精品免费久久久久久久清纯| 看黄色毛片网站| 亚洲第一青青草原| 国产区一区二久久| 99国产极品粉嫩在线观看| 午夜老司机福利片| 国产一级毛片七仙女欲春2 | 国产又色又爽无遮挡免费看| 99久久国产精品久久久| 午夜视频精品福利| 免费搜索国产男女视频| 国产精品98久久久久久宅男小说| 日本 欧美在线| 国产精品免费视频内射| 视频在线观看一区二区三区| 久久国产精品影院| 麻豆av在线久日| 在线十欧美十亚洲十日本专区| 女性生殖器流出的白浆| 国产av一区在线观看免费| 午夜免费观看网址| 日韩欧美国产在线观看| 村上凉子中文字幕在线| 两人在一起打扑克的视频| 久久香蕉精品热| 亚洲精华国产精华精| 777久久人妻少妇嫩草av网站| 叶爱在线成人免费视频播放| 美国免费a级毛片| 日韩欧美免费精品| 热re99久久国产66热| 视频区欧美日本亚洲| 午夜免费鲁丝| 亚洲精品美女久久av网站|