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

    An electrochemical machining method for aeroengine blades based on four-directional synchronous feeding of cathode tools

    2023-10-25 12:13:08JiaLIULibingHUIDongqianJIAYanLIUDiZHU
    CHINESE JOURNAL OF AERONAUTICS 2023年9期

    Jia LIU, Libing HUI, Dongqian JIA, Yan LIU, Di ZHU

    College of Mechanical and Electrical Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

    KEYWORDS

    Abstract Blades are critical components of modern aero-engine.Among the many characteristic structures of the blade,the leading/trailing edges are key structures that have the greatest influence on the aerodynamic effect and power conversion of the blade.Electrochemical machining(ECM)is regarded as one of the most important techniques in blade manufacturing due to its process characteristics of high material removal rates, virtually no tool wear, and no areas of thermal or mechanical damage to the workpiece rim zones.Herein, an ECM method based on the fourdirectional synchronous feeding of four cathode tools is proposed to improve the machining accuracy of the leading and trailing edges of the blade.During blade ECM using this method,four cathode tools feed toward the basin/back surfaces and leading/trailing edges respectively.The dynamic processing simulation and flow field simulation results of the ECM process show that the proposed method eliminates the sharp changes in the electric field and electrolyte flow field at the leading and trailing edges seen in traditional machining methods.Thus, the electric field and flow field stability of the leading and trailing edges at the final stage of machining is greatly improved.Experimental comparison of the conventional and proposed ECM methods showed that four-directional synchronous feeding results in improved profile accuracy over repeated machining processes and good surface quality.

    1.Introduction

    Blades are core components of aero-engines and take the pole of energy transformation.The leading and trailing edges are the most important functional areas of a blade.During the operation of an aero-engine, the airflow divides and merges at the leading and trailing edges, respectively.A small profile error in the leading or trailing edge will significantly affect the aerodynamic performance of the aero-engine.1–3Therefore,among all areas of the blade, the leading and trailing edges have the highest requirements for profile accuracy.However,the leading and trailing edges are complex, weak, and rigid structures;they are also ultra-thin and twisty with small curvature.4In addition, aero-engine blades are usually made of materials that are difficult to cut(e.g.,nickel-based superalloys and titanium alloys).5–7In the process of milling,the vibration and deformation of the workpiece can seriously affect the machining accuracy.Occasionally, subsurface cracks caused by cutting heat also affect the surface integrity of the blade.8Thus, it is very difficult to control the machining accuracy of the leading and trailing edges of blades.

    Electrochemical machining (ECM) is a non-traditional machining technique9in which materials are removed via a controlled electrochemical anodic reaction.ECM has the advantages of no tool wear and good surface quality regardless of the material hardness and toughness.10,11Based on these characteristics, ECM has become one of the main processing methods in aero-engine blade manufacturing.12Experts and scholars have conducted detailed research on blade ECM.Demirtas et al.studied the effects of two ECM parameters,feed velocity and the electrical conductivity of the electrolyte,on freeform surface machining.13Brusilovski analyzed the influence of pulse voltage, pulse width, electrolyte pressure,and other parameters on machining accuracy in pulse ECM.14Paczkowski et al.established a mathematical model of ECM and evaluated the distribution of the machining gap.15Fujisawa et al.established a three-dimensional simulation model of an ECM blade and predicted the distribution of machining clearance.16Li and Niu proposed a model and a finite element numerical approach for cathode design based on the potential distribution of the inter-electrode gap.17Ernst et al.developed a rapid method for design method based on simulation and optimization.18Liu et al.proposed a multidirectional auxiliary electrolyte flow mode for the ECM of blisk.19Klocke et al.established a multi-physical coupled material removal model20and a surface micro-morphology prediction model for ECM21and numerically simulated the ECM process of the leading and trailing edges of the blade.22Zhu et al.studied the design strategy of cathode tool trajectory for the ECM of blisk23and proposed a tool design method to improve the machining accuracy of the leading and trailing edges of the blade.24Wang and Zhu proposed an ECM method based on tangential feeding in which the leading/trailing edges are electrochemically processed by the cathode tools,which are fed along the tangential direction of the mean camber line of the blade.25.

    In the conventional ECM of blades, the processing mode usually involves two cathode tools feeding toward the basin and back surfaces of the blade.The basin/back surfaces and leading/trailing edges of the blade are thus processed at the same time.In this processing mode, the ECM accuracy of the basin and back surfaces can be well controlled.However,the leading and trailing edges are at both ends of the opening gap between the two cathode tools.The continuous feeding of the two cathode tools causes drastic changes in the electric field and electrolyte flow field at the leading and trailing edges,which seriously affects the accuracy of the ECM.In addition,the current density at the leading/trailing edges increases rapidly in the final stage of the ECM process, which is detrimental to maintaining process stability and consistency.For aero-engine blades, the accuracy requirement for the leading and trailing edges is often higher than that for the basin and back surfaces.To improve the ECM accuracy of the leading and trailing edges, a new ECM method for blades based on the four-directional synchronous feeding of four cathode tools is proposed in this study.In the proposed method, four cathode tools are arranged in a cross shape to form a closed curved surface.A specially designed slip surface is located between each pair of adjacent cathode tools.When the four cathode tools are fed according to the set feed velocity ratio,the sliding surface between any two adjacent cathodes always remains in contact,and the surface surrounded by the four tools is always closed.During blade ECM,four cathode tools feed toward the basin/back surfaces and leading/trailing edges.All the profile surfaces of the blade are processed simultaneously with the contraction of the closed surface surrounded by the tools.The proposed method eliminates the opening electric and flow field in the leading/trailing edges, which improves the stability of the electric field and electrolyte flow field,thereby improving the controllability of machining accuracy.In order to verify the feasibility of the proposed method, the simulation and experiment were carried out.

    2.Principle of four-directional synchronous feeding method

    The schematic of the traditional ECM method for blades is shown in Fig.1.In this method, the workpiece is connected to the positive pole of power supply, and the cathode tools are connected to the negative poles.A high-velocity, highpressure electrolyte passes through the machining gaps between the cathode tools and the basin/back surfaces.During processing, the two cathode tools feed toward the basin and back surfaces, respectively.The material of the workpiece is gradually dissolved into ions via an electrochemical anodic dissolution reaction.After a period of time,the dissolution velocity of the workpiece material in the basin and back surfaces is roughly balanced with the feed velocity of the cathode tools.In this processing mode,the ECM accuracy of the basin and back surfaces can be tightly controlled.The leading and trailing edges, however, are located at the opposite ends of the gap between the two cathode tools.The leading and trailing edges do not have cathode tools feeding along their direction.Due to the large machining gap,the materials at the leading and trailing edges are always in a free dissolution state.It is hard to precisely control the shape of the leading and trailing edges without the cathode tool profile constraint.In addition, with the continuous face-to-face feeding of the two cathode tools,the shape of the machining gaps at the leading and trailing edges changes dramatically, especially in the final stage of machining.This will cause significant changes in the electric field and the electrolyte flow field in this region, resulting in random fluctuations of current density in the machining gaps.In ECM,the corrosion removal rate of the material is proportional to the current density on the workpiece surface.Thus,continuous and drastic changes in the electric field and electrolyte flow field will seriously affect repeatability.As a result,in the traditional ECM process, it is difficult to control the ECM accuracy of the leading and trailing edges.

    Fig.1 Schematic diagram of traditional blade ECM process.

    In order to improve the ECM accuracy of the leading and trailing edges, a new ECM method for blades based on the four-directional synchronous feeding of four cathode tools is proposed.A schematic of the proposed ECM method is shown in Fig.2.The workpiece is connected with the positive pole of the power supply.Four cathode tools are arranged in a cross shape and connected to the negative pole of the power supply.The four tools form a closed surface,and a sliding surface separates every-two adjacent tools.When the four tools are fed at a certain velocity ratio, sliding motion only occurs between adjacent tools, and the surface surrounded by the tools remains closed at all times.During blade ECM, the four tools are fed toward the basin/back surfaces and the leading/trailing edges, respectively.All the surfaces of the blade are processed at the same time.In this ECM method,two specially designed cathode tools are fed to the leading and trailing edges respectively.With the continuous feeding of the cathode tools, the electric field distribution in the machining gap of the leading and trailing edges will enter the equilibrium state.In this ECM state, the controllability of the leading and trailing edge shapes will be greatly improved.Furthermore, since the surface surrounded by the tools remain closed during the ECM process, the drastic changes in the electric field and electrolyte flow field of the leading and trailing edges are caused by the gap between the cathode tools, which exist in the traditional ECM method,are eliminated.The proposed method improves the stability of the electric field and electrolyte flow field,which helps to improve the controllability of the machining accuracy.

    Fig.2 Schematic diagram of proposed blade ECM process.

    To ensure that the sliding surfaces of two adjacent tools remain in contact at all times during the synchronous movement of four tools, the relationship between the dip angle of the sliding surface and the feed velocity ratio of two adjacent tools must be determined.The method used to determine the dip angle of the sliding surface is shown in Fig.3.Take the sliding surface between the basin-surface tool and the leading-edge tool as an example.The feed directions of the two tools are V1and V2,the angles between the sliding surface and the two tools are α and β, the angle between α and β is γ,and the feed velocity ratio V1/V2is K.The included angles α and β can be calculated using the Eqs.(1)and(2),respectively.By using this method to design four sliding surfaces,the structures of the four synchronous feeding tools can be determined.

    3.Simulation of the forming processes in blade ECM

    3.1.Simulation model and boundary conditions

    The simulation geometric models of the blade machining gaps for the traditional and proposed ECM methods are shown in Fig.4.The workpiece blank is a rectangle with a length of 48 mm, a width of 8 mm and a height of 36 mm.The cathode tool profiles of the two methods are designed according to the typical stator blade model of an aero-engine.In the proposed four-directional feeding method,the cathode tools are perpendicular to each other and are fed at the same velocity, and the four dip angles of the sliding surfaces are all 45°.In both ECM modes, the electrolyte flows from the tip surface to the root surface of the blade,and the direction of electrolyte flow is perpendicular to the simulation section.Therefore, the effects of electrolytic products on electrolyte conductivity were ignored in the simulation, and the conductivity distribution in the machining gap was assumed to be uniform and constant.

    Fig.3 Method for determining dip angle of a slip surface.

    Fig.4 Simulation geometric models of blade machining gap.

    In the simulation models shown in Fig.4, the boundaries Γc1and Γp1are the workpiece surfaces,and the boundary conditions are φa= U.The boundaries Γc2to Γc3and Γp2to Γp9are the cathode tool surfaces,and the boundary conditions are φc= 0.Other boundaries such as Γc4and Γp5are the insulation surfaces of the fixture, and the boundary conditions are?φ/?n = 0.In the ECM process, the dissolution of the workpiece material is dominated by the electric field, which is described by the Laplace equation:

    In ECM, the material dissolution velocity νnobeys Faraday’s law, as shown in Eq.(4):

    where η is the electrolytic processing efficiency,ω is the electrochemical equivalent volume of the workpiece material,and i is the current density on the surface of the workpiece.

    During the simulation of the traditional and proposed ECM methods, all tools were fed to the workpiece at the constant velocity of ν.The entire ECM processing period was divided into small time intervals (Δt).For each time interval Δt,the electric field distribution in the machining gap was first calculated using Eq.(3).Next,according to the distribution of current density on the anode surface, the displacement of the anode boundary Δl was calculated using Eqs.(4) and (5).At the same time, the tool boundary also translates the distance Δd along the feed direction (Fig.5), as shown in Eq.(6).

    Fig.5 Simulation calculation method of ECM process.

    The new workpiece and tool boundaries are then taken to the next iteration, and the above calculation is repeated.Repeat the above simulation steps until the profile of the workpiece enters the tolerance range of the theoretical contour,and the final profile of the workpiece is obtained.The simulation parameters and conditions are shown in Table 1.The current efficiency curve of the workpiece material (SS304 stainless steel) used in the simulations was obtained from the study by Wang et al.26.

    3.2.Simulation results and analysis

    The total time of the simulation process is 1200 s.Since the leading and trailing edges molding occurs at the last stage of the machining process,the most drastic changes in the leading and trailing edges profiles occur during this time, so the simulation results within the last 120 s are selected for processing and analysis.

    Fig.6 shows the simulation results (current density distribution and anode formation process)for the leading and trailing edges of the traditional ECM method at 1080 s,1100 s,???,1200 s during the machining process.With the continuous feeding of the two cathode tools,the electric field distributions and profiles at the leading and trailing edges of the blade change drastically during traditional ECM.In order to quantify the simulation results, two measurement boundaries were intercepted at the leading edge (points A to C in Fig.6) and trailing edge (points D to F).Eleven reference points were evenly selected within these two boundaries, where points B and E represent the vertices of the leading and trailing edges respectively, as shown in Fig.7.

    Fig.8 shows the current density distribution of the leading and trailing edges at the moments (t = 1160 s, 1180 s and 1200 s) during ECM.In the final stage of ECM, the electric field is obviously concentrated near the vertices, and the current density near points B and E increases rapidly.The rapid increase in local current density will accelerate the dissolutionof local workpiece material.This will significantly reduce the controllability of machining accuracy.

    Table 1 Simulation conditions of forming processes.

    Fig.6 Simulation results of current density distribution for traditional ECM method.

    Fig.9 shows the simulation results of the anode formation process based on the proposed ECM method.With the synchronous feeding of the four cathode tools, the feed velocity of the cathode boundaries is approximately equal to the corrosion velocity of the anode boundary.Thus,the machining process occurs in an approximately equilibrium state, and the electric field distribution in the entire machining gap is relatively stable.As for the traditional ECM simulation,two measurement boundaries were selected in the areas of the leading and trailing edges (points U to W and points X to Z, respectively, where points V and Y are vertices).During the last 40 seconds of ECM, the current density at each point in the two research boundaries basically remains stable, as shown in Fig.10.The uniform variation in the anode profile will significantly improve the controllability of ECM accuracy at the leading and trailing edges.

    Fig.7 Location of sampling points.

    Fig.8 Current density distribution in final stage of traditional ECM method.

    Fig.9 Simulation results of current density distribution for proposed ECM method.

    4.Simulation of the electrolyte flow field in blade ECM

    4.1.Simulation model and boundary conditions

    To analyze the effect of different processing modes on the stability of the electrolyte flow field, the simulation of the flow field in traditional and proposed methods during the last 40 seconds was carried out.In the simulation,the electrolyte flow fields at two moments(t=1160 s and t=1200 s)were chosen as the object of analysis.Since the velocity of electrolyte is much higher than the velocity of tool feed and workpiece material dissolution,the motion of tool and workpiece boundary is neglected.Based on the simulation results of the forming process, the geometric models of the electrolyte flow channel for these two moments are designed, as shown in Fig.11.

    Fig.10 Current density distribution in final stage of the proposed ECM method.

    In ECM,both the traditional and proposed methods adopt the flow mode of electrolyte from the blade tip surface to the root surface of the blade.In the simulation model, the inlet and outlet boundaries of the electrolyte are set as pressure type.The pressures of electrolyte inlet and outlet are set to 0.8 MPa and 0.1 MPa,respectively.Other surfaces in the electrolyte flow channel are set to be walls.The RNG k-ε model is used to accurately describe the turbulent flow of electrolyte.Its turbulent kinetic energy k and dissipation rate ε are shown in Eq.(7) and Eq.(8).27.

    Fig.11 Geometric models of electrolyte flow field.

    where ρ is the density of the electrolyte; t is the time; u is the flow velocity; μ is the effective viscosity coefficient; Gkis the turbulent kinetic energy generation term caused by the average velocity gradient; and μtis the turbulent viscosity.Additionally, the constant= C1ε= 1.42, C2ε= 1.68,Cμ= 0.0845, αk= 1, αε= 1.

    4.2.Simulation results and analysis

    Fig.12 shows the simulation results of the electrolyte velocity distributions at the leading and trailing edges of the traditional and proposed methods at the two moments.In the traditional ECM method, the electrolyte velocity at the leading and trailing edges of the blade decreases significantly as the gap between the two tools decreases.By contrast, the apparent shape change of the electrolyte flow channel near the leading and trailing edges is eliminated in the proposed ECM method due to the synchronous feeding of the four cathode tools.

    In order to quantify the simulation results, eleven measuremnt lines were selected on the cross-section of the electrolyte flow channel at the leading and trailing edges.1100 reference points were set uniformly on these lines, as shown in Fig.13.Using the flow velocity distribution cloud maps of the leading and trailing edges at two moments (t = 1160 s and t = 1200 s) of both methods, the flow velocity data of 100 reference points on each measurement line were extracted and the mean values were taken for analysis, and the results were shown in Fig.14.In the traditional ECM method, the average flow velocity is about 16 m/s.Additionally, the flow velocity fluctuates very sharply at the leading and trailing edges.The flow velocity at the leading edge ranged from 20.6 m/s to 15.14 m/s at 1160 s, however, it is only 16.37 m/s to 15.81 m/s at 1200 s.On the contrary,in the proposed ECM method, the average flow velocity is about 19 m/s and there is little fluctuation in the flow velocity at the leading and trailing edges.The flow velocity at the leading edge varies from 19.46 m/s to 18.36 m/s at 1160 s and from 19.18 m/s to 17.92 m/s at the last moment.The flow velocity fluctuation at the trailing edge is similar to those at the leading edge.Furthermore, in the traditional ECM method, the variance of the electrolyte flow velocity for the leading edge is 23.5 m2/s2and 21.7 m2/s2, respectively, at 1160 s and 1200 s.However, in the proposed method,the variance of the flow velocity is 6.5 m2/s2and 6.8 m2/s2,respectively.The stable electrolyte flow field can reduce the random error in the ECM and significantly improve the consistency of repeated processing.

    Fig.12 Simulation results of electrolyte velocity distribution.

    Fig.13 Location of measurement line.

    Fig.14 Electrolyte velocity analysis results.

    Table 2 Experimental conditions.

    Fig.15 ECM system and corresponding process devices.

    Fig.16 Specimens machined using traditional and proposed ECM methods.

    Fig.17 Blade middle profiles of traditional and proposed ECM methods.

    Fig.18 Deviation distribution of each profile with respect to average profile.

    5.Experiment

    To verify the feasibility of the proposed method, comparison experiments of the two ECM methods were carried out.The experiments for both methods were simplified as two halfblade tests on the leading-edge and trailing-edge parts.Fig.15 shows the ECM system and corresponding devices for the traditional and proposed methods.In the newly proposed ECM method, three cathode tools are fed at a constant velocity toward the basin, back, and leading-edge (or trailingedge) surfaces of the blade, respectively.A high-speed, highpressure electrolyte flows in from the blade tip and out from the blade root.During traditional ECM, two cathode tools are fed toward the basin and back surfaces of the blade,respectively, and the gap above the leading edge (or trailing edge) is filled with insulating blocks.The workpiece is a rectangle with a length of 28 mm, width of 8 mm, and height of 36 mm.The workpieces and tools are made of SS304 stainless steel.The fixture is made of glass fiber-reinforced plastic.The experimental parameters are consistent with the simulation parameters, as shown in Table 1.The experimental electrolyte and other experimental parameters are given in Table 2.

    Fig.16 shows the leading- and trailing-edge parts of the specimens machined using the traditional and proposed ECM methods.Each experiment was repeated three times.On each specimen, three section lines located at the blade tip,middle of the blade,and blade root were selected for measurement and analysis.In ECM, the profile repeatability between multiple machined specimens is an important index for evaluating the accuracy of the machining method.Therefore, the analysis of the experimental results primarily focuses on the repeatability of the profiles of the machined specimens.

    The blade middle profiles near the leading and trailing edges of the specimens are shown in Fig.17.The three profiles obtained for the same ECM method and position were averaged to obtain the average profiles at the leading and trailing edges (shown as dotted lines in Fig.16).The distributions of deviation from the average profiles are shown for both ECM methods in Fig.18.For the traditional ECM method,the deviation from the average profile ranges from 0.036 to - 0.025 mm (from the positive direction to the negative direction) for the leading edge and from 0.061 to-0.040 mm for the trailing edge.In comparison, the deviations obtained using the proposed ECM method range from 0.011 to - 0.013 mm at the leading edge and from 0.012 to - 0.017 mm at the trailing edge.Thus, the proposed method significantly improved the machining accuracy at the leading and trailing edges.These results demonstrate that the deviation in the profile at the middle of the blade from the average profile is obviously smaller for the proposed method than for the traditional method.

    Fig.19 Repeated deviation distributions for traditional and proposed ECM methods.

    To evaluate the repeated ECM accuracy of the three section lines on the blade,the difference between the positive deviation and negative deviation of the measuring point from the average profile was defined as the evaluation index.Fig.19 shows the repeated deviation distributions for the blade tip, blade middle, and blade root of the specimens processed using the traditional and proposed ECM methods.For the traditional ECM method, the maximum deviations on the leading and trailing edges were 0.064 and 0.086 mm, respectively.For the proposed ECM method,the maximum deviations on the leading and trailing edges were 0.036 and 0.038 mm, respectively.Thus,the four-directional synchronous feeding method significantly improved the machining accuracy at the leading and trailing edges compared to the traditional ECM method.

    Fig.20 shows the measured surface roughness values obtained for the two ECM methods.The surface roughness of the specimens processed using both methods were better than Ra0.4 μm.Overall, the experimental results demonstrate the feasibility and effectiveness of the proposed ECM method.

    Fig.20 Surface 3D profile and roughness.

    6.Conclusions

    A method for the ECM of blades based on the four-directional synchronous feeding of four cathode tools is presented in this paper.In this method, four cathode tools are fed toward the basin/back surfaces and leading/trailing edges, respectively,and all profile surfaces of the blade are processed at the same time.The new method was evaluated through simulation and experiment, and the conclusions are summarized as follows:

    1) In the conventional two-directional feeding method for the ECM of blades, the materials at the leading and trailing edges are always in a free dissolution state,and the electric field changes dramatically during the ECM process.The simulation results indicate that the distributions of current density and electrolyte velocity at the leading/trailing edges change sharply in the final stage of the ECM process, which is detrimental to the control of machining accuracy.

    2) The proposed ECM method eliminates the sharp changes in the electric field and electrolyte flow field at the leading and trailing edges.In the final stage of blade ECM,the distributions of electric field and flow field are approximately in a steady state,which improves the controllability of machining accuracy.

    3) The experimental comparison of the traditional and proposed ECM methods demonstrated that the profiles of the leading and trailing edges have higher accuracy when machined using four-directional synchronous feeding compared to traditional ECM.The surfaces obtained using the proposed ECM method also have good surface quality.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (Grant 52075253) and the National Science and Technology Major Project (Grant 2017-VII-0004-0097).

    久久久久久久亚洲中文字幕 | 又粗又爽又猛毛片免费看| 天天躁日日操中文字幕| 成人永久免费在线观看视频| 亚洲精品粉嫩美女一区| 国产色婷婷99| 国产私拍福利视频在线观看| www.熟女人妻精品国产| 久久精品国产亚洲av香蕉五月| 麻豆av噜噜一区二区三区| 欧美日本亚洲视频在线播放| 老司机午夜十八禁免费视频| 波多野结衣高清无吗| 亚洲人成电影免费在线| 日韩免费av在线播放| 人妻久久中文字幕网| 一区二区三区激情视频| 男人的好看免费观看在线视频| 国产av一区在线观看免费| 大型黄色视频在线免费观看| 熟女人妻精品中文字幕| 91午夜精品亚洲一区二区三区 | 国产精品精品国产色婷婷| 特级一级黄色大片| 嫩草影院新地址| 欧美不卡视频在线免费观看| 欧美在线一区亚洲| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精华国产精华精| 亚洲精品在线美女| 男人和女人高潮做爰伦理| 国内少妇人妻偷人精品xxx网站| 老司机午夜十八禁免费视频| 亚洲欧美日韩无卡精品| 老司机福利观看| 久久精品国产亚洲av香蕉五月| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 琪琪午夜伦伦电影理论片6080| 在线国产一区二区在线| 在线播放无遮挡| 好看av亚洲va欧美ⅴa在| 成人鲁丝片一二三区免费| 中国美女看黄片| 欧美一区二区国产精品久久精品| eeuss影院久久| 欧美最黄视频在线播放免费| 一级毛片久久久久久久久女| 不卡一级毛片| 亚洲人成伊人成综合网2020| 国产一区二区在线观看日韩| 亚洲av成人精品一区久久| 搡老熟女国产l中国老女人| 小蜜桃在线观看免费完整版高清| eeuss影院久久| 99热这里只有是精品50| 性色av乱码一区二区三区2| 久久性视频一级片| 欧美色视频一区免费| 天美传媒精品一区二区| 97超视频在线观看视频| 亚洲欧美日韩无卡精品| 色哟哟哟哟哟哟| 小说图片视频综合网站| 动漫黄色视频在线观看| 婷婷色综合大香蕉| 久久精品影院6| 免费看a级黄色片| 十八禁国产超污无遮挡网站| 日本免费a在线| 欧美三级亚洲精品| 首页视频小说图片口味搜索| 国产免费av片在线观看野外av| 日日夜夜操网爽| 69人妻影院| 一进一出抽搐gif免费好疼| 性欧美人与动物交配| 丰满人妻一区二区三区视频av| 国产熟女xx| 欧美成狂野欧美在线观看| 黄色女人牲交| 亚洲午夜理论影院| 搡女人真爽免费视频火全软件 | 成熟少妇高潮喷水视频| 99在线人妻在线中文字幕| 伊人久久精品亚洲午夜| 久久九九热精品免费| 欧美日本亚洲视频在线播放| 男女那种视频在线观看| 免费观看人在逋| 国产高潮美女av| 午夜福利成人在线免费观看| 国产精品三级大全| 国产毛片a区久久久久| 午夜精品久久久久久毛片777| 男女下面进入的视频免费午夜| av专区在线播放| 日本 欧美在线| av中文乱码字幕在线| 观看免费一级毛片| 51午夜福利影视在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产精品亚洲av一区麻豆| 五月玫瑰六月丁香| 黄色配什么色好看| 成人av在线播放网站| 国产主播在线观看一区二区| 欧美日本视频| 给我免费播放毛片高清在线观看| 亚洲国产精品sss在线观看| 日韩大尺度精品在线看网址| 亚洲国产精品成人综合色| 久久精品人妻少妇| 中文在线观看免费www的网站| 制服丝袜大香蕉在线| 成人毛片a级毛片在线播放| 精品人妻一区二区三区麻豆 | 国产黄色小视频在线观看| 12—13女人毛片做爰片一| 亚洲国产精品久久男人天堂| 亚洲性夜色夜夜综合| 日韩欧美 国产精品| 久久久久精品国产欧美久久久| 免费大片18禁| 国产精品综合久久久久久久免费| 一本精品99久久精品77| 最近中文字幕高清免费大全6 | 成人鲁丝片一二三区免费| 国产亚洲欧美在线一区二区| 永久网站在线| 国产大屁股一区二区在线视频| 亚洲专区中文字幕在线| 亚洲成a人片在线一区二区| 国产av不卡久久| 日韩中文字幕欧美一区二区| 午夜福利成人在线免费观看| 免费大片18禁| 丰满乱子伦码专区| 色精品久久人妻99蜜桃| 欧美日韩乱码在线| 琪琪午夜伦伦电影理论片6080| 国产高清视频在线观看网站| 高清日韩中文字幕在线| 免费高清视频大片| 日韩中文字幕欧美一区二区| 国产综合懂色| 给我免费播放毛片高清在线观看| 欧美高清成人免费视频www| 日本免费一区二区三区高清不卡| 日本精品一区二区三区蜜桃| 亚洲精品影视一区二区三区av| 99久国产av精品| 欧美+日韩+精品| 一级a爱片免费观看的视频| 我要看日韩黄色一级片| 桃色一区二区三区在线观看| 在线播放无遮挡| 日本五十路高清| 亚洲av成人av| 一区二区三区免费毛片| 亚洲五月天丁香| 国产亚洲欧美98| netflix在线观看网站| 人妻夜夜爽99麻豆av| 国产精品,欧美在线| 少妇裸体淫交视频免费看高清| 国产亚洲欧美98| 老司机深夜福利视频在线观看| 在线观看一区二区三区| 欧美性感艳星| 一级作爱视频免费观看| 中文在线观看免费www的网站| 波野结衣二区三区在线| 99久久精品热视频| 午夜福利在线在线| 国产三级黄色录像| 嫩草影院入口| 精品午夜福利视频在线观看一区| 免费看光身美女| 午夜a级毛片| 男插女下体视频免费在线播放| 日韩大尺度精品在线看网址| 亚洲 国产 在线| av在线天堂中文字幕| 国产v大片淫在线免费观看| 黄色配什么色好看| 男人舔女人下体高潮全视频| 亚洲国产精品成人综合色| 国产综合懂色| 日韩免费av在线播放| 最近视频中文字幕2019在线8| 能在线免费观看的黄片| 国产精品人妻久久久久久| 欧美日本亚洲视频在线播放| 91在线精品国自产拍蜜月| 欧美中文日本在线观看视频| 欧美不卡视频在线免费观看| 免费电影在线观看免费观看| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区免费欧美| 身体一侧抽搐| 九九久久精品国产亚洲av麻豆| 日韩欧美精品免费久久 | 中文字幕av成人在线电影| 最近视频中文字幕2019在线8| 最近最新免费中文字幕在线| 亚洲精品456在线播放app | 网址你懂的国产日韩在线| 国产成年人精品一区二区| 一进一出抽搐动态| 国产亚洲精品久久久久久毛片| 97超视频在线观看视频| 欧美成人免费av一区二区三区| 亚洲电影在线观看av| 免费大片18禁| 真人做人爱边吃奶动态| 在线免费观看的www视频| 麻豆av噜噜一区二区三区| 波野结衣二区三区在线| www.色视频.com| 免费看美女性在线毛片视频| 久久精品国产99精品国产亚洲性色| 女同久久另类99精品国产91| 91在线精品国自产拍蜜月| 看免费av毛片| 久久精品国产清高在天天线| 一级av片app| 精品一区二区三区人妻视频| 啪啪无遮挡十八禁网站| 久久久久久九九精品二区国产| 欧美zozozo另类| 成人一区二区视频在线观看| 色吧在线观看| 成人亚洲精品av一区二区| 在线观看av片永久免费下载| 天天躁日日操中文字幕| 国产午夜精品论理片| 久久久久久久久久黄片| 久久伊人香网站| 亚洲av成人不卡在线观看播放网| 久9热在线精品视频| 日韩精品中文字幕看吧| 色综合站精品国产| 中文在线观看免费www的网站| 国产午夜精品久久久久久一区二区三区 | 亚洲人成网站高清观看| 男女做爰动态图高潮gif福利片| 老司机午夜十八禁免费视频| 久久热精品热| 好男人在线观看高清免费视频| 伦理电影大哥的女人| 午夜福利免费观看在线| 怎么达到女性高潮| 99热只有精品国产| 级片在线观看| 日本 av在线| 中文字幕精品亚洲无线码一区| av天堂中文字幕网| 亚洲国产精品成人综合色| 国产精品人妻久久久久久| 久久久久久久久久黄片| 午夜久久久久精精品| 亚洲成人精品中文字幕电影| 国内精品久久久久久久电影| 日韩免费av在线播放| 高清毛片免费观看视频网站| 白带黄色成豆腐渣| 淫秽高清视频在线观看| 天堂av国产一区二区熟女人妻| 香蕉av资源在线| 在线十欧美十亚洲十日本专区| 成年女人永久免费观看视频| 99热这里只有是精品50| 日韩欧美一区二区三区在线观看| 99热这里只有精品一区| 美女免费视频网站| 久久精品国产99精品国产亚洲性色| 一进一出抽搐gif免费好疼| 97热精品久久久久久| 99国产精品一区二区三区| 宅男免费午夜| 直男gayav资源| 精品人妻视频免费看| 亚洲人与动物交配视频| 久久热精品热| 欧美日韩福利视频一区二区| 日本黄大片高清| 免费搜索国产男女视频| 少妇高潮的动态图| 成人午夜高清在线视频| 午夜免费男女啪啪视频观看 | 亚洲美女视频黄频| 亚洲人成网站在线播| 国产三级在线视频| 99热这里只有是精品在线观看 | 最新在线观看一区二区三区| 亚洲真实伦在线观看| 一本久久中文字幕| 国产精品野战在线观看| 午夜福利18| 精品熟女少妇八av免费久了| 身体一侧抽搐| 欧美色视频一区免费| 欧美日韩综合久久久久久 | www.www免费av| 少妇高潮的动态图| 日韩av在线大香蕉| 老熟妇乱子伦视频在线观看| 性色av乱码一区二区三区2| 18+在线观看网站| 亚洲自偷自拍三级| 又黄又爽又刺激的免费视频.| 草草在线视频免费看| 又爽又黄无遮挡网站| 国产三级中文精品| 国产中年淑女户外野战色| 一进一出抽搐gif免费好疼| 免费观看的影片在线观看| 国语自产精品视频在线第100页| 嫩草影院新地址| 中文资源天堂在线| 免费看美女性在线毛片视频| 欧美日韩国产亚洲二区| 18禁在线播放成人免费| 国产精品亚洲美女久久久| 欧美日韩中文字幕国产精品一区二区三区| 天堂av国产一区二区熟女人妻| www.999成人在线观看| 国产又黄又爽又无遮挡在线| 欧美午夜高清在线| 国产又黄又爽又无遮挡在线| 美女高潮喷水抽搐中文字幕| 亚洲av成人精品一区久久| 91麻豆精品激情在线观看国产| 免费黄网站久久成人精品 | 亚洲在线自拍视频| 日韩欧美三级三区| 国产在线精品亚洲第一网站| 97超视频在线观看视频| 欧美一级a爱片免费观看看| 国产精品久久久久久精品电影| 婷婷丁香在线五月| 亚洲片人在线观看| 色综合欧美亚洲国产小说| 久久国产乱子免费精品| 精品日产1卡2卡| 久久人人精品亚洲av| 97超级碰碰碰精品色视频在线观看| 成人毛片a级毛片在线播放| 色精品久久人妻99蜜桃| 少妇的逼水好多| 免费av毛片视频| 国产日本99.免费观看| 国产爱豆传媒在线观看| 非洲黑人性xxxx精品又粗又长| 大型黄色视频在线免费观看| 亚洲欧美激情综合另类| 久久99热这里只有精品18| 欧美日韩亚洲国产一区二区在线观看| av欧美777| 成人亚洲精品av一区二区| 精品久久久久久久久亚洲 | 欧美激情久久久久久爽电影| 欧美在线一区亚洲| 又紧又爽又黄一区二区| 性色avwww在线观看| 在线观看免费视频日本深夜| 欧美乱妇无乱码| 午夜福利在线观看免费完整高清在 | 制服丝袜大香蕉在线| 国产成+人综合+亚洲专区| 美女xxoo啪啪120秒动态图 | 亚洲最大成人手机在线| 成人av一区二区三区在线看| 中文字幕精品亚洲无线码一区| 男人舔奶头视频| bbb黄色大片| 91字幕亚洲| 别揉我奶头~嗯~啊~动态视频| 成人av一区二区三区在线看| 国产中年淑女户外野战色| 人人妻,人人澡人人爽秒播| aaaaa片日本免费| 免费看a级黄色片| 欧美黄色淫秽网站| av国产免费在线观看| 成年女人毛片免费观看观看9| 国产真实伦视频高清在线观看 | 成年版毛片免费区| 日韩欧美精品免费久久 | 97超级碰碰碰精品色视频在线观看| 人人妻人人澡欧美一区二区| 色综合欧美亚洲国产小说| 少妇丰满av| 日韩欧美在线二视频| 精品国产三级普通话版| 久久国产精品人妻蜜桃| 长腿黑丝高跟| 亚洲成av人片在线播放无| 深爱激情五月婷婷| 黄色一级大片看看| 久久午夜福利片| 一区二区三区激情视频| 欧美日韩福利视频一区二区| 自拍偷自拍亚洲精品老妇| 美女高潮的动态| 丁香欧美五月| 黄色丝袜av网址大全| 国产91精品成人一区二区三区| 久久久久久久亚洲中文字幕 | 少妇熟女aⅴ在线视频| 91久久精品国产一区二区成人| 亚洲av二区三区四区| 一边摸一边抽搐一进一小说| 国产极品精品免费视频能看的| 亚洲国产色片| 日日摸夜夜添夜夜添小说| 亚洲av日韩精品久久久久久密| 国产精品人妻久久久久久| 我要搜黄色片| 波野结衣二区三区在线| 精品欧美国产一区二区三| 两个人的视频大全免费| 日本撒尿小便嘘嘘汇集6| 国产三级在线视频| 精品久久久久久成人av| 国产91精品成人一区二区三区| 夜夜躁狠狠躁天天躁| 欧美日韩中文字幕国产精品一区二区三区| 国产精品99久久久久久久久| 久久亚洲真实| 欧美日韩黄片免| 成人性生交大片免费视频hd| 亚洲av电影在线进入| 亚洲欧美精品综合久久99| 久久精品国产亚洲av香蕉五月| 国产 一区 欧美 日韩| 国产成年人精品一区二区| av专区在线播放| 久久久精品欧美日韩精品| 伊人久久精品亚洲午夜| 午夜老司机福利剧场| 精品国产亚洲在线| 99热6这里只有精品| 丁香欧美五月| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 色视频www国产| 少妇被粗大猛烈的视频| 国产高清三级在线| 综合色av麻豆| 人妻制服诱惑在线中文字幕| 欧美xxxx黑人xx丫x性爽| 久久精品91蜜桃| 国产美女午夜福利| 国产真实伦视频高清在线观看 | 亚洲人成网站在线播放欧美日韩| 日韩有码中文字幕| 91av网一区二区| 国产av一区在线观看免费| www.999成人在线观看| 日韩欧美在线乱码| 国产精品久久久久久人妻精品电影| 日韩中文字幕欧美一区二区| 日韩 亚洲 欧美在线| 深夜a级毛片| 99久久成人亚洲精品观看| 97超视频在线观看视频| 18+在线观看网站| 岛国在线免费视频观看| 亚洲第一区二区三区不卡| 在线免费观看不下载黄p国产 | 欧美激情在线99| 亚洲 欧美 日韩 在线 免费| 中亚洲国语对白在线视频| 日韩精品中文字幕看吧| 国产高清有码在线观看视频| 一个人看视频在线观看www免费| 一区二区三区免费毛片| 在线观看av片永久免费下载| 精品午夜福利视频在线观看一区| 97超视频在线观看视频| 国产精品不卡视频一区二区 | 国产精品亚洲一级av第二区| 简卡轻食公司| 国产一区二区在线观看日韩| 欧美一级a爱片免费观看看| 国产精品精品国产色婷婷| h日本视频在线播放| 麻豆国产av国片精品| 最好的美女福利视频网| 精品久久久久久久久久久久久| 国产三级黄色录像| 黄色一级大片看看| 青草久久国产| 午夜免费成人在线视频| 高清在线国产一区| 九九热线精品视视频播放| 少妇高潮的动态图| 99国产精品一区二区蜜桃av| 757午夜福利合集在线观看| 性色av乱码一区二区三区2| 亚洲最大成人中文| 日日干狠狠操夜夜爽| 网址你懂的国产日韩在线| www.www免费av| 亚洲第一区二区三区不卡| 国产欧美日韩一区二区精品| 亚洲中文字幕一区二区三区有码在线看| 桃红色精品国产亚洲av| 日本黄色视频三级网站网址| 国产精华一区二区三区| 精品不卡国产一区二区三区| 中国美女看黄片| 在线观看66精品国产| 亚洲精华国产精华精| 欧美日韩综合久久久久久 | 午夜福利免费观看在线| 99久久99久久久精品蜜桃| 国产精品精品国产色婷婷| 亚洲一区二区三区不卡视频| 免费大片18禁| 亚洲va日本ⅴa欧美va伊人久久| av视频在线观看入口| 又黄又爽又免费观看的视频| 久久6这里有精品| 人妻制服诱惑在线中文字幕| 免费搜索国产男女视频| 一级黄色大片毛片| 高清毛片免费观看视频网站| 亚洲精品在线美女| 国产欧美日韩一区二区精品| 伦理电影大哥的女人| 亚洲精品久久国产高清桃花| 又紧又爽又黄一区二区| 国产亚洲av嫩草精品影院| 精品一区二区三区人妻视频| 熟女人妻精品中文字幕| 一进一出抽搐gif免费好疼| 麻豆成人午夜福利视频| 国产三级中文精品| 亚洲人与动物交配视频| 波多野结衣高清无吗| 美女 人体艺术 gogo| 最好的美女福利视频网| 性插视频无遮挡在线免费观看| 在线观看66精品国产| 国产亚洲av嫩草精品影院| 欧美三级亚洲精品| 国产淫片久久久久久久久 | 嫩草影院精品99| 欧美日韩国产亚洲二区| 搡老妇女老女人老熟妇| 757午夜福利合集在线观看| 免费搜索国产男女视频| 亚洲人成网站在线播放欧美日韩| 毛片女人毛片| ponron亚洲| 99久久精品国产亚洲精品| 亚洲欧美日韩无卡精品| 蜜桃久久精品国产亚洲av| 国产91精品成人一区二区三区| 日韩大尺度精品在线看网址| 噜噜噜噜噜久久久久久91| 亚洲一区二区三区不卡视频| 免费搜索国产男女视频| 亚洲成人免费电影在线观看| 中文字幕精品亚洲无线码一区| 久久精品国产亚洲av香蕉五月| 9191精品国产免费久久| 午夜免费激情av| 搡老岳熟女国产| 丁香六月欧美| 亚洲精品久久国产高清桃花| 日韩欧美精品v在线| 99热精品在线国产| 久久亚洲真实| 人人妻人人澡欧美一区二区| 欧美一级a爱片免费观看看| 91久久精品国产一区二区成人| 亚洲 国产 在线| 国产综合懂色| 日日干狠狠操夜夜爽| 日韩免费av在线播放| 黄色丝袜av网址大全| 日韩欧美免费精品| 日韩欧美 国产精品| 丰满人妻熟妇乱又伦精品不卡| 美女xxoo啪啪120秒动态图 | 中出人妻视频一区二区| 成年女人永久免费观看视频| 在线观看免费视频日本深夜| 免费看日本二区| 精品久久久久久久久av| 在线a可以看的网站| 亚洲一区高清亚洲精品| 精品久久久久久久久av| 日本成人三级电影网站| 国产单亲对白刺激| 国产亚洲欧美在线一区二区| 亚洲成人免费电影在线观看| 欧美成人免费av一区二区三区| 亚洲无线在线观看| 动漫黄色视频在线观看| 中文字幕免费在线视频6| 丰满乱子伦码专区| 精品国产亚洲在线| 免费黄网站久久成人精品 | 国产视频一区二区在线看| 中出人妻视频一区二区| 亚洲无线在线观看| 永久网站在线| 精品久久久久久成人av| 18禁在线播放成人免费| 一进一出抽搐gif免费好疼| 午夜久久久久精精品| 久久6这里有精品|