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

    Development of an In-Situ Laser Machining System Using a Three-Dimensional Galvanometer Scanner

    2020-06-01 02:06:20XiaoLiBinLiuXuesongMeiWenjunWangXiaodongWangXunLi
    Engineering 2020年1期

    Xiao Li, Bin Liu,*, Xuesong Mei, Wenjun Wang, Xiaodong Wang, Xun Li

    a State Key Laboratory for Manufacturing Systems Engineering, Xi’an Jiaotong University, Xi’an 710049, China

    b Shaanxi Key Laboratory of Intelligent Robots, Xi’an Jiaotong University, Xi’an 710049, China

    c Xi’an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences, Xi’an 710119, China

    Keywords:

    In-situ laser machining

    Three-dimensional galvanometer scanner Line structured light

    Three-dimensional measurement

    ABSTRACT In this study,a three-dimensional(3D)in-situ laser machining system integrating laser measurement and machining was built using a 3D galvanometer scanner equipped with a side-axis industrial camera.A line structured light measurement model based on a galvanometer scanner was proposed to obtain the 3D information of the workpiece.A height calibration method was proposed to further ensure measurement accuracy, so as to achieve accurate laser focusing. In-situ machining software was developed to realize time-saving and labor-saving 3D laser processing. The feasibility and practicability of this in-situ laser machining system were verified using specific cases.In comparison with the conventional line structured light measurement method,the proposed methods do not require light plane calibration,and do not need additional motion axes for 3D reconstruction; thus they provide technical and cost advantages. The insitu laser machining system realizes a simple operation process by integrating measurement and machining, which greatly reduces labor and time costs.

    1. Introduction

    Over the past 20 years, laser machining technology has found its way into many industrial applications, such as welding [1-3],cutting [4,5], drilling [6,7], texturing [8-10], and micro-structures manufacturing [11-13]. With the development of laser machining technology and equipment, the application of this technology has gradually become popularized, from mass production to small batch customized production. It is essential to develop a laser machining technique that is more suitable for diversified products,and the advent of the galvanometer scanner greatly assists in this endeavor.A galvanometer scanner can quickly and accurately control the position of the laser spot. Compared with traditional processing methods, laser processing with a galvanometer scanner has many advantages, including high dynamic performance and processing speed, no wear on mechanical tools, contact-free processing, and exceptional flexibility without retooling [14-16].

    The major challenge of three-dimensional(3D)laser machining is to continually keep the laser focus spot on a specific position in 3D space[17,18].Noh et al.[19]proved that the microgrooves produced by an on-focus laser possess higher quality in terms of depth and uniformity compared with those produced by an off-focus laser. Cao et al. [17,20] used machine vision to compensate for a defocused laser and to continually keep the laser focused on the surface. Wang et al. [21] divided a large freeform surface into sub-areas, sub-layers, and sub-blocks to ensure that the laser was always focused. Based on the 3D surface of the workpiece,dynamic adjustment of the laser focus spot in the spatial position is the key to realizing non-planar laser machining. The 3D galvanometer scanner improves the laser beam by allowing it to change its focus position rapidly in 3D space; the 3D machining capability of the galvanometer scanner benefits from its dynamic focusing unit[22,23].Xiao et al.[24]used a 3D galvanometer scanner to directly mark patterns on freeform curved surfaces. Diaci et al. [25] used a 3D galvanometer scanner to mark complete patterns on a curved surface and inclined surface,which could not be done using a general 2D galvanometer scanner.

    Most current 3D galvanometer scanner machining systems are unable to perform 3D measurement; therefore, it is necessary to obtain the 3D information of the workpiece by means of 3D software modeling or special scanning equipment before laser machining, and the workpiece must be positioned by fixtures.However, when workpieces are produced in a small batch and are diversified, a great deal of time will be spent on preparation work,resulting in increased hardware and time costs.At the same time, multiple processes may introduce more errors. In-situ processing technology generally refers to direct in-situ measurement and processing on an object;this technology has been investigated in the fields of industry, bioengineering, and more [26-29]. When this technology is used in 3D laser machining, its main advantage lies in the integration of measurement and machining, which can simplify the traditional workflow—including individual modeling,clamping, adjustment, and alignment—and avoid the workload and errors caused by redundant processes.The galvanometer scanner can project line structured light for 3D measurement [30-32].Diaci et al. [25] attempted to use a laser system for measurement and machining in order to complete rapid and flexible laser marking and engraving;however,they did not provide a specific implementation method for measurement or a detailed evaluation of the results.

    Traditional 3D measurement using line structured light is mainly based on the direct projection of laser lines onto the surface of an object[33].Zhou et al.[34]used line structured light to measure rails; they found that the errors in length, width, and thickness were about 0.3 mm. Liu et al. [35] measured the diameter of an axis (the diameter was about 30 mm) by means of line structured light. The diameter was obtained by fitting the measured points, and the error was found to be less than 20 μm. Li et al.[36] used line structured light to measure a cube with a height of 100 mm,and observed an error of about 1 mm.In fact,the absolute measurement accuracy is determined by factors such as the equipment, measurement distance, and algorithm. For example, a small camera field of view and a short measurement distance can result in relatively high measurement accuracy.

    The current research on line structured light measurement method mainly includes measurement modeling,the optical center extraction method, and the calibration method. Measurement modeling is the basis of this technology, and is mainly based on geometric derivation. The measurement modeling using a galvanometer scanner is different from the typical method due to the use of the galvanometer scanner and different types of focusing lens. Optical center extraction is an important factor in measurement accuracy.The optical center extraction method with subpixel precision can achieve relatively high measurement accuracy, and there are mainly three types of typical subpixel methods including the Gaussian fitting method [37], gray centroid method [5], and Hessian matrix method [38]. Calibration mainly involves camera parameters calibration and light plane equation calibration. The calibration of the light plane equation is more difficult, and the problem of how to obtain high-precision control points to fit the light plane is a key issue in calibration.

    In this study, a 3D in-situ laser machining system integrating laser measurement and machining was built using a 3D galvanometer scanner and an industrial camera. The main contribution of this paper is the realization of the integration of in-situ measurement and machining, including: ①the proposal of a line structured light measurement model based on a galvanometer scanner; ②the proposal of a height calibration method to ensure measurement accuracy; and ③ the development of in-situ machining software to realize time-saving and labor-saving 3D laser processing. The feasibility and practicability of this in-situ laser machining system were verified through specific cases. In comparison with the conventional line structured light measurement method, the proposed methods do not require light plane calibration,and do not require additional motion axes for 3D reconstruction. The proposed in-situ laser machining system realizes a simple operation process,which reduces labor and time costs.

    2. Methods

    2.1. 3D measurement method

    Fig. 1(a) shows the 3D measurement model based on the line structured light projected by the galvanometer scanner. There are four coordinate systems in this model: the world coordinate system Ow-xwywzw,the galvanometer coordinate system Og-xgygzg,the camera coordinate system Oc-xcyczc, and the image coordinate system O-uv. To simplify the model, the galvanometer coordinate system is coincident with the world coordinate system. In this research, a non-telecentric f-θ field lens was used. In Fig. 1(a),point a is the center of the galvanometer scanner’s y-axis mirror,whose world coordinates are a(0, 0, H), where H is the height of point a in the world coordinates. Points b and c are the two endpoints of the laser line projected by the galvanometer scanner.Their world coordinates can be precisely controlled by the galvanometer scanner, and are b(x1, y1, 0) and c(x2, y2, 0). To obtain the 3D coordinates of a point on the object, at least three constraints should be known.According to the typical pinhole camera model [39], Eq. (1) can be derived.

    Fig. 1. Principle and equipment diagram of the 3D measurement subsystem. (a) Schematic diagram of the measurement method based on line structured light using the galvanometer scanner; (b) equipment diagram of the measurement and machining system.

    where fxand fyare the focal distance expressed in units of horizontal and vertical pixels; and cxand cyare the principal point coordinates. Suppose that a space point d, (xc, yc, zc) has the coordinates of d in the camera coordinate system, and (u, v) are the coordinates in the pixel coordinate system.In order to uniquely determine the coordinate value of the point, at least one more constraint is required—that is,the light plane equation.The conventional line structured light measurement method requires calibration of the light plane to determine the light plane equation, but the light plane generated by the galvanometer scanner can be precisely controlled. The laser is focused by the field lens, which includes the non-telecentric f-θ field lens and the telecentric f-θ field lens. When using a non-telecentric f-θ field lens, the light plane in the world coordinate system passes through the three points a, b, and c (Fig. 1(a)); when using a telecentric f-θ field lens,the light plane is always perpendicular to the xy plane in the world coordinate system,and the world coordinates of the point a in Fig.1 are considered to be a(x1, 0, H). Therefore, the three known noncollinear points (a, b, and c) in the world coordinate system can be converted to the camera coordinate system according to the camera’s external parameters, and the light plane equation in the camera coordinate system can be determined as follows:

    where A,B,C,and D are the coefficients of the light plane equation.This is synchronized with Eq. (1) to obtain Eq. (3):

    Thus, coordinates on the pixel coordinate system (u,v) can be converted to the camera coordinate system (xc, yc, zc) according to Eq. (3), and then to the world coordinate system (xw, yw, zw)according to Eq. (4), where R is the rotation matrix and T is the translation matrix in the camera’s external parameters.

    2.2. Height calibration method

    On the basis of the methods described above, this study introduces the use of a height calibration method to ensure the measurement accuracy of the height direction. The principle of this method is to find the relationship among the actual height, the laser line offset distance, and the pixel distance; this relationship is shown in Eq. (5).

    where Z is the actual height, x is the laser line offset distance, p is the pixel distance, and f represents the function of x and p. In this study, the laser line is parallel to the y-axis direction in the world coordinate system,and its offset direction is along the x-axis direction, as shown in Fig. 1(a). The three-axis motion platform can be controlled to move a certain distance in the z-axis direction and simultaneously take multiple laser line images for height calibration.

    The center coordinate of the laser line in the pixel coordinate system is the optical center. The optical center with a height of 0 and a laser line offset distance of 0 is defined as the reference optical center,and the average pixel distance between any other optical center and the reference optical center is defined as P. In the plane with a height of 0,the average pixel distance between other optical centers and the reference optical center is defined as P1.When the laser line offset distance is the same and the height is different,the average pixel distance between optical centers at different heights and the optical center at a height of 0 is defined as P2. Therefore, Eq. (6) is obtained in the pixel coordinate system.

    Fig. 2(a) shows the relationship between P1and x (-30 to 30 mm, calibrated every 5 mm). Fig. 2(b) shows the relationship between P2and Z (-15 to 15 mm, calibrated every 5 mm). It can be seen from the figures that P1has a good linear relationship with x, and P2has a good linear relationship with Z.

    Fig.3 plots the relationship among Z,P2,and x,which is almost fitted to a plane.Therefore,when calculating the actual height Z,P1is first calculated by linear interpolation according to the value of x in Fig. 2(a); next, according to the calculated value of P, P2can be obtained using Eq.(6).Finally,Z is obtained by interpolation using P2and x, as shown in Fig. 3.

    Fig.2. (a)Relationship between the average pixel distance P1 and the laser line offset distance x in the height of 0;(b)relationship between the average pixel distance P2 and the height Z in the laser line offset distance of 0.

    Fig. 3. 3D map of the relationship among the actual height Z, pixel distance P2, and laser line offset distance x.

    2.3. Optical center extraction method

    In this method, the laser line images are preprocessed by median filtering, and the gray centroid method is then used to extract the optical center. First, the point gmaxwith the highest light intensity in a column of pixels is found, and then a gray threshold K =k×gmax(0 <k <1) is determined. Here, k is a scale factor, and the pixel coordinates of all points whose gray value is larger than K are set to ui(i=0,1,2···); and the corresponding gray values of these points are set to gi(i=0,1,2···). Finally, the optical center U is obtained by Eq. (7):

    During the extraction of the optical center by the gray centroid method, the thresholds of each column of pixels are different, so even if the gray distribution is not uniform, it will not bring great errors to the determination of the optical center position, thus improving the extraction accuracy.

    2.4. In-situ machining method

    After the extraction of the optical center, the 3D point cloud surface of the object can be obtained according to the methods described in Sections 2.1 and 2.2. The density of the point cloud is determined by the distance between the projected laser lines and the resolution of the industrial camera.The point cloud surface can reflect the actual 3D shape and pose of the object.

    After obtaining the surface point cloud of the object,machining patterns or paths are projected to the corresponding position on the point cloud to generate the 3D galvanometer scanner processing instructions. In this study, the machining patterns were projected onto the 3D point cloud in the form of pixels. Each pixel in the machining pattern could find its adjacent points in the point cloud, and the average height value of these points was taken as the height value of the pixel. In order to realize in-situ machining,in-situ machining software was developed,which achieves the precise projection of patterns at any position of the workpiece by zooming, translating, and rotating the imported patterns.

    3. Experiments and evaluation

    3.1. Equipment condition

    The 3D galvanometer scanner used in this study was a SCANLAB intelliSCANse14 equipped with a 255 mm focal length nontelecentric f-θ field lens. The industrial camera was a Vieworks VH-5MG 5-million-pixel black-and-white camera, equipped with a U-TRON FV1520 lens, and the camera was fixed on the side of the 3D galvanometer scanner. The field of view of the industrial camera was about 200 mm×200 mm when the working distance was 350 mm.The laser source was a 1064 ns fiber laser(RFL-P30Q)with a 680 nm indicator light. The 3D galvanometer scanner was mounted on a three-axis motion platform, as shown in Fig. 1(b).

    3.2. 3D measurement accuracy evaluation

    The measurement accuracy directly decides the laser focusing accuracy during in-situ processing. A standard-height step block was selected (Fig. 4(a)) with a length of 80 mm, a width of 20 mm, and four standard heights of 3, 5, 10, and 15 mm for the measurement accuracy test. The scanned laser line spacing was 0.5 mm, the scanning speed of the galvanometer scanner was 5000 mm·s-1, and the exposure time of the industrial camera was 50 ms. 3D point clouds were generated using Method 1 (the method using the camera parameters and the light plane equation)and Method 2 (the method that includes height calibration).Fig.4(b)shows the point cloud obtained by Method 2,which accurately reconstructs the height information of the standard-height step block. Fig. 5(a) indicates the average value of the measured height. Fig. 5(b) shows that the error of Method 1 varies with change in the y value. The error at both ends is large, and reaches a maximum at a height of 15 mm,where it exceeds 1 mm.In contrast, the error of Method 2 is stable in the height direction. It is always within the range of ±0.2 mm, and does not change significantly with change in the y value.

    The y value of the first and last points of the laser lines scanned onto the workpiece were extracted, and the length of the workpiece obtained by the scanning was calculated to be(79.900 ± 0.154) mm, which is consistent with the actual workpiece length of 80 mm.

    3.3. Implementation and examples

    3.3.1. Case 1: In-situ machining on a standard-height step block

    In order to verify the function of the 3D in-situ laser machining system developed in this study, 3D measurement, pattern projection, and in-situ machining were carried out on a standard-height step block (Fig. 4(a)). This process included the following steps:

    Step 1: Pre-calibration work. The Matlab camera calibration toolbox was used to perform camera calibration in order to obtain the camera’s internal and external parameters.When obtaining the external parameters, the x-axis and y-axis direction lines were projected to the chessboard to coincide the galvanometer coordinate system with the world coordinate system, as shown in Fig. 6(a). Height calibration was performed using the method described in Section 2.2.

    Fig. 4. (a) The standard-height step block with standard heights of 3, 5, 10, and 15 mm; (b) 3D point cloud of the standard-height step block obtained by Method 2.

    Fig. 5. (a) Average height values obtained by Method 1 and Method 2; (b) height error of Method 1 and Method 2.

    Step 2:3D in-situ measurement.The scanning speed of the galvanometer scanner was 5000 mm·s-1, the exposure time of the camera was 50 ms, and the laser line was scanned every 0.5 mm to traverse the object. One of the laser line images is shown in Fig. 6(b). The optical center extracted by the gray center method is illustrated in Fig. 6(c). The 3D point cloud of the workpiece surface was generated according to the methods described in Sections 2.1 and 2.2.

    Step 3: Machining pattern projection. The generated 3D point cloud and the imported machining patterns were displayed by OpenGL.After a certain amount of zooming,translating,and rotating, the pattern was projected and displayed on the point cloud.

    Step 4: 3D in-situ machining. In order to protect the standardheight step block, thermal paper was attached to it, and the laser was used to mark the patterns on the thermal paper. After in-situ machining, these patterns were consistent with the design patterns, as shown in Fig. 6(d).

    In this study,the in-situ machining software(Fig.7)was developed on a C# WinForm platform. The main functions of the software include scanning laser lines, importing the point cloud,importing machining patterns, translating and rotating the patterns, generating projection patterns, and controlling the galvanometer scanner and the three-axis motion platform.

    3.3.2. Case 2: Local scanning and in-situ machining

    Fig. 6. Implementation process and results of Case 1. (a) The galvanometer coordinate system is coincided with the world coordinate system during camera calibration;(b)one of the laser line images during 3D measurement;(c)the optical center extracted from the laser line image, where the red line is the extracted optical center; (d) in-situ machining results on the thermal paper.

    When the workpiece is large, it is only necessary to scan the local part that needs to be processed in order to further reduce the unnecessary workload. In this case, a larger non-standard workpiece was selected, as shown in Fig. 8(a), the machining results were shown in Fig. 8(b). It is only necessary to ensure that the local part to be processed is within the processing range of the galvanometer scanner and the visual field of the camera.Figs.8(c)and (d) show the measurement results of a local corner of the workpiece.Two quick response(QR)code patterns were projected on it, and the machining was directly processed after the projection. Designed patterns can be accurately marked at the corresponding position of the workpiece, which is very convenient.From the machining results, it can be seen that the processed patterns are complete and have no obvious flaws, whether on the slope or on the plane;the position and size of the in-situ machined QR code patterns are consistent with the design in the software.

    Fig.7. The in-situ 3D machining software interface displays the 3D point cloud(brown dots)of the standard-height step block and of the machining patterns before and after projection.

    Fig.8. Implementation process and results of Case 2.(a)A larger size non-standard workpiece on the machining platform;(b)in-situ machining results of two QR code patterns at the corner of the workpiece; (c) projection results on the slope;(d) projection results on the plane.

    3.3.3. Case 3: Machining on a curved surface with different pixel distance

    In order to verify the machining capacity on a curved surface,an aluminum alloy arc surface workpiece(40 mm×40 mm×11.5 mm)was selected as the processing object, and the pattern was projected after equal arc length transformation.The machining line spacing was artificially adjusted by the pixel distance in the in-situ machining software. The measurement and projection results are presented in Fig. 9(a). A machining line spacing of 40 and 120 μm was used, respectively. The machining results are shown in Figs. 9(b) and (c). When the machining line spacing was 40 μm, the machined pattern was filled (Fig. 9(b)); when it was changed to 120 μm,there was a gap between each machining line(Fig. 9(c)). Optical microscopy observation in the range of the red dotted box (Fig. 9) showed that the line width produced by the laser was(46.08±1.22)μm,and the actual machining line spacing was (40.04 ± 0.72) and (120.11 ± 0.02)μm, respectively, which is consistent with the design values; this indicates that the system has high processing accuracy.

    3.3.4.Case 4:Tool in-situ measurement and micro-textures processing

    Fig. 9. Implementation process and results of Case 3. (a) Measurement and projection results of the aluminum alloy arc surface workpiece. In-situ machining results with(b) 40 μm machining line spacing and (c) 120 μm machining line spacing. Optical microscopy observation results with (d) 40 μm machining line spacing and (e) 120 μm machining line spacing.

    Fig.10. Implementation process and results of Case 4.(a)Triangular turning tool with a particular rake angle;(b)reconstruction results of the tool;(c)top view of the tool profile (the dimensions in the figure are measured by DigitalMicrograph software); (d) in-situ machining software screenshot (bottom right corner is the top view of the projection of the micro-textures); (e) in-situ machining results of the micro-textures. L1, L2, L3: the side length without considering the tool nose radius.

    Table 1 The in-situ measured value and the Vernier caliper measured value of the tool.

    4. Discussion

    Unlike the conventional line structured light measurement method, in which the laser line is directly generated by the line laser,the laser line in this study was generated by a galvanometer scanner. Because the position of the laser line in the common method cannot be accurately controlled,it is necessary to calibrate the light plane to obtain the light plane equation; a great deal of research has been conducted on this topic. For example, Liu et al.[40] used a ball target and nonlinear optimization to perform the light plane calibration, while Kiddee et al. [41] projected crossline structured light(CLSL)on the chessboard to calculate the light plane equation. These calibration methods are complicated and require accurate calibration targets. In this study, the position of the laser line can be precisely controlled by the galvanometer scanner, thus eliminating the need for light plane calibration, which greatly simplifies the process and reduces the workload.3D reconstruction by means of the typical line structured light measurement method requires additional linear or rotary motion axes to achieve laser line scanning of the workpiece. However, the galvanometer scanner can realize self-scanning, thereby eliminating the need to add additional motion axes, and consequently saving hardware cost. Therefore, line structured light measurement method using a galvanometer scanner has technical and cost advantages.

    Compared with a 2D galvanometer scanner,a 3D galvanometer scanner must acquire the 3D information of the workpiece in order to accurately focus the laser. However, most laser machining systems equipped with a 3D galvanometer scanner lack the 3D measurement capability, making it necessary to obtain the 3D information of the workpiece through other methods or equipment in advance;this is inconvenient and will increase the time cost and workload. The proposed methods add 3D measurement capability to the 3D laser machining system with just one industrial camera.The typical 3D laser machining process mainly includes 3D modeling, importing the models, and positioning the workpieces by the fixture. In comparison, the proposed methods do not require fixture positioning. In-situ measurement can obtain the 3D information of the workpiece and match the machining coordinate system, so that the laser can be focused on a certain position of the workpiece.The above method realizes a simple operation process for laser 3D machining, which greatly reduces time and labor costs.

    When the workpiece size is large, the farther away a point on the workpiece surface is from the origin of the world coordinate system, the more obvious the measurement error will be. This error is due to inaccurate calibration of the camera parameters. If the measured height is inaccurate,the laser focus will be defocused during the 3D laser processing, which will directly affect the processing quality. In order to ensure measurement accuracy in the height direction,the height calibration method was proposed.This method ensures that the height error is always within a certain range, which is generally smaller than the laser Rayleigh length in the laser machining system, thereby ensuring that the laser is always focused.The error of the height calibration method mainly comes from the line structured light not being thin enough and the resolution of the camera not being high; however, the error is stable and within an acceptable range.

    5. Conclusion

    In this study,an a 3D in-situ laser machining system integrating laser measurement and machining was built by using a 3D galvanometer scanner equipped with a side-axis industrial camera.The line structured light was projected onto the workpiece surface by the galvanometer scanner, which increased the 3D measurement capability for the 3D laser machining system. A line structured light measurement model based on the galvanometer scanner was proposed. A method using camera parameters and the light plane equation was used to obtain the 3D surface information. A height calibration method was proposed to further ensure the measurement accuracy. The measurement error of these methods is stable and acceptable in the height direction.The feasibility and practicability of this in-situ laser machining system were verified using specific cases.

    The advantages of the proposed methods are as follows:①Compared with the conventional line structured light measurement method, the proposed methods do not require light plane calibration and additional motion axes to achieve 3D reconstruction,and thus provide technical and cost advantages;②the in-situ laser machining system realizes a simple operation process by the integration of 3D in-situ measurement and laser machining,which greatly reduces labor and time costs.

    For tasks that require high machining accuracy,the in-situ laser machining system can realize higher precision laser processing on a macro-micro scale by configuring a more precise coaxial vision device (such as SCANLAB’s camera adapter), which has great research value and application potential.

    Acknowledgements

    This work was supported by the National Key Research and Development Program of China (2017YFB1104602 and 2016YFB1102502) and the Program for Changjiang Scholars and Innovative Research Team in University (IRT_15R54).

    Compliance with ethics guidelines

    Xiao Li, Bin Liu, Xuesong Mei, Wenjun Wang, Xiaodong Wang,and Xun Li declare that they have no conflict of interest or financial conflicts to disclose.

    白带黄色成豆腐渣| 亚洲精品粉嫩美女一区| 日本免费a在线| www.www免费av| 啪啪无遮挡十八禁网站| 亚洲无线观看免费| 久久精品国产综合久久久| 波多野结衣巨乳人妻| 中文字幕久久专区| 美女被艹到高潮喷水动态| 天天躁日日操中文字幕| 欧美性猛交黑人性爽| 宅男免费午夜| 90打野战视频偷拍视频| 淫秽高清视频在线观看| 琪琪午夜伦伦电影理论片6080| 国产精品女同一区二区软件 | 午夜a级毛片| 国产伦精品一区二区三区视频9 | 免费在线观看影片大全网站| 欧美+亚洲+日韩+国产| 亚洲自拍偷在线| 国产精品香港三级国产av潘金莲| 热99在线观看视频| 久久精品国产清高在天天线| 五月玫瑰六月丁香| 舔av片在线| 免费大片18禁| 男女视频在线观看网站免费| 国产日本99.免费观看| 亚洲精品乱码久久久v下载方式 | 夜夜看夜夜爽夜夜摸| 色在线成人网| 日韩欧美精品v在线| 精品国产三级普通话版| 国产一区二区激情短视频| 亚洲精品美女久久av网站| 淫秽高清视频在线观看| 亚洲欧美日韩卡通动漫| 麻豆国产av国片精品| 亚洲成a人片在线一区二区| 欧洲精品卡2卡3卡4卡5卡区| xxx96com| 亚洲va日本ⅴa欧美va伊人久久| 麻豆久久精品国产亚洲av| 欧美午夜高清在线| 露出奶头的视频| 亚洲国产精品成人综合色| 亚洲五月婷婷丁香| 特大巨黑吊av在线直播| 亚洲中文日韩欧美视频| 一本综合久久免费| 黄色成人免费大全| 999久久久国产精品视频| 18禁国产床啪视频网站| 欧美激情在线99| 精品无人区乱码1区二区| 亚洲av电影不卡..在线观看| 可以在线观看毛片的网站| 久久精品91蜜桃| 日本 av在线| 久久中文字幕一级| 日韩有码中文字幕| 亚洲欧美日韩高清在线视频| 九九在线视频观看精品| 日韩精品青青久久久久久| 亚洲片人在线观看| 日日干狠狠操夜夜爽| 桃色一区二区三区在线观看| 啦啦啦韩国在线观看视频| 动漫黄色视频在线观看| 国产精品免费一区二区三区在线| cao死你这个sao货| 精品日产1卡2卡| 99久久99久久久精品蜜桃| 国产乱人视频| 免费电影在线观看免费观看| 亚洲精品中文字幕一二三四区| 久久久久性生活片| 国产成人av激情在线播放| av天堂中文字幕网| 亚洲国产精品999在线| 长腿黑丝高跟| 中国美女看黄片| 欧美成人免费av一区二区三区| 亚洲最大成人中文| 久久中文看片网| 两性夫妻黄色片| 高潮久久久久久久久久久不卡| 老熟妇仑乱视频hdxx| 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 不卡一级毛片| 国产精品女同一区二区软件 | 97超视频在线观看视频| 国内精品久久久久久久电影| 白带黄色成豆腐渣| 变态另类丝袜制服| 免费看日本二区| 久久久久国产精品人妻aⅴ院| 亚洲成人中文字幕在线播放| 母亲3免费完整高清在线观看| svipshipincom国产片| 女同久久另类99精品国产91| 国产三级黄色录像| 精品日产1卡2卡| 在线a可以看的网站| 国产成年人精品一区二区| 午夜精品久久久久久毛片777| 亚洲中文日韩欧美视频| 热99在线观看视频| 香蕉国产在线看| 国产高清videossex| 精品国产乱子伦一区二区三区| 12—13女人毛片做爰片一| 成人特级av手机在线观看| av片东京热男人的天堂| 黄色日韩在线| 九色成人免费人妻av| 美女 人体艺术 gogo| 欧美乱码精品一区二区三区| 免费av毛片视频| 两性夫妻黄色片| 亚洲中文日韩欧美视频| 亚洲成人中文字幕在线播放| 欧美成狂野欧美在线观看| 综合色av麻豆| 在线观看免费午夜福利视频| 老司机深夜福利视频在线观看| 国产伦精品一区二区三区四那| 日韩精品青青久久久久久| 一个人免费在线观看电影 | 宅男免费午夜| e午夜精品久久久久久久| 国内精品一区二区在线观看| 一二三四社区在线视频社区8| 在线免费观看不下载黄p国产 | 日韩高清综合在线| 精品国产乱子伦一区二区三区| av在线蜜桃| 国产av不卡久久| 国产日本99.免费观看| 色老头精品视频在线观看| 午夜免费观看网址| 成年免费大片在线观看| 全区人妻精品视频| 国产一区二区在线av高清观看| aaaaa片日本免费| 国产精品乱码一区二三区的特点| 丰满人妻一区二区三区视频av | 99久国产av精品| av福利片在线观看| 嫩草影院精品99| 一卡2卡三卡四卡精品乱码亚洲| 久久国产精品影院| 国产精品精品国产色婷婷| 人妻久久中文字幕网| 午夜两性在线视频| 午夜亚洲福利在线播放| 99热精品在线国产| 一区福利在线观看| 久99久视频精品免费| 老司机深夜福利视频在线观看| 色综合亚洲欧美另类图片| 成人性生交大片免费视频hd| 亚洲色图av天堂| www.999成人在线观看| x7x7x7水蜜桃| 亚洲电影在线观看av| a级毛片a级免费在线| 校园春色视频在线观看| 国产爱豆传媒在线观看| 亚洲国产精品成人综合色| 欧美一级a爱片免费观看看| 人人妻,人人澡人人爽秒播| 国产黄色小视频在线观看| 精品久久久久久久久久久久久| 亚洲一区二区三区不卡视频| 亚洲国产看品久久| 亚洲成av人片免费观看| 亚洲男人的天堂狠狠| 日本黄色视频三级网站网址| 日本免费a在线| 欧美大码av| 九色成人免费人妻av| 成人鲁丝片一二三区免费| 小说图片视频综合网站| 我要搜黄色片| 观看免费一级毛片| 美女午夜性视频免费| 两性夫妻黄色片| 久久精品亚洲精品国产色婷小说| 波多野结衣高清作品| 欧美3d第一页| 国产高清视频在线播放一区| 欧美成人性av电影在线观看| 日日夜夜操网爽| 亚洲欧美日韩无卡精品| 精品一区二区三区视频在线 | 在线观看一区二区三区| 精品熟女少妇八av免费久了| 欧美日韩综合久久久久久 | 久久国产精品人妻蜜桃| 啦啦啦免费观看视频1| 真人一进一出gif抽搐免费| 搡老熟女国产l中国老女人| 国产精品一及| 老司机福利观看| 精品国产超薄肉色丝袜足j| 亚洲九九香蕉| 99热精品在线国产| 18禁国产床啪视频网站| 熟女少妇亚洲综合色aaa.| 久久久久久九九精品二区国产| av视频在线观看入口| 可以在线观看毛片的网站| 久久亚洲精品不卡| 女人被狂操c到高潮| 精品国产超薄肉色丝袜足j| 我要搜黄色片| 在线国产一区二区在线| 美女高潮喷水抽搐中文字幕| 校园春色视频在线观看| 午夜影院日韩av| 国产亚洲av嫩草精品影院| 熟妇人妻久久中文字幕3abv| 视频区欧美日本亚洲| 国产精品 欧美亚洲| 日韩精品青青久久久久久| 久久久久国内视频| 又黄又爽又免费观看的视频| 国产精品精品国产色婷婷| 国产高清三级在线| 日韩成人在线观看一区二区三区| 中文字幕精品亚洲无线码一区| 嫁个100分男人电影在线观看| 熟女少妇亚洲综合色aaa.| 国产乱人视频| 久久午夜亚洲精品久久| 黑人巨大精品欧美一区二区mp4| 噜噜噜噜噜久久久久久91| 久久午夜综合久久蜜桃| 国产高清激情床上av| 欧美色视频一区免费| 91老司机精品| 国内揄拍国产精品人妻在线| 国产欧美日韩一区二区精品| 在线观看舔阴道视频| 久久久成人免费电影| 久久久久性生活片| 午夜视频精品福利| 日本一二三区视频观看| 免费av不卡在线播放| 国产精品精品国产色婷婷| 怎么达到女性高潮| 日韩欧美三级三区| 日本在线视频免费播放| 日韩 欧美 亚洲 中文字幕| 午夜激情欧美在线| 国产高潮美女av| 成人无遮挡网站| 熟妇人妻久久中文字幕3abv| 一个人看视频在线观看www免费 | 一本一本综合久久| 午夜两性在线视频| 成人无遮挡网站| 日本熟妇午夜| 午夜亚洲福利在线播放| 中文资源天堂在线| 十八禁网站免费在线| 国产成人啪精品午夜网站| 国产av在哪里看| 男人和女人高潮做爰伦理| 亚洲色图av天堂| 俄罗斯特黄特色一大片| 中文字幕av在线有码专区| 亚洲av美国av| 天堂影院成人在线观看| 热99在线观看视频| 美女扒开内裤让男人捅视频| 少妇人妻一区二区三区视频| 免费在线观看亚洲国产| 男人舔奶头视频| 欧美乱妇无乱码| 国产精品1区2区在线观看.| 制服丝袜大香蕉在线| 波多野结衣巨乳人妻| 亚洲av美国av| 观看免费一级毛片| 免费电影在线观看免费观看| 亚洲欧洲精品一区二区精品久久久| 少妇人妻一区二区三区视频| 国产97色在线日韩免费| 最好的美女福利视频网| 久久久久性生活片| 午夜福利在线在线| x7x7x7水蜜桃| 欧美又色又爽又黄视频| x7x7x7水蜜桃| 国产精品影院久久| 午夜激情欧美在线| 97碰自拍视频| 一本一本综合久久| 亚洲无线观看免费| 国产精品久久久久久精品电影| 亚洲狠狠婷婷综合久久图片| 在线免费观看的www视频| 精品一区二区三区视频在线 | 青草久久国产| 美女被艹到高潮喷水动态| 久久精品人妻少妇| 我要搜黄色片| www.熟女人妻精品国产| 欧美一区二区精品小视频在线| 综合色av麻豆| www.精华液| 1000部很黄的大片| 久久中文字幕一级| 99久久成人亚洲精品观看| 999久久久国产精品视频| 国产高清有码在线观看视频| 久久人妻av系列| 国产精品久久久久久亚洲av鲁大| 日韩av在线大香蕉| 啦啦啦观看免费观看视频高清| 国产av在哪里看| 18美女黄网站色大片免费观看| 成人精品一区二区免费| 国产又色又爽无遮挡免费看| tocl精华| 国产免费男女视频| 日韩欧美三级三区| 性欧美人与动物交配| 成年女人看的毛片在线观看| 国产精品女同一区二区软件 | 每晚都被弄得嗷嗷叫到高潮| 一夜夜www| av黄色大香蕉| 欧美激情久久久久久爽电影| 国内精品美女久久久久久| av在线天堂中文字幕| 色尼玛亚洲综合影院| 精品久久久久久,| 亚洲九九香蕉| 欧美一级a爱片免费观看看| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲av嫩草精品影院| 韩国av一区二区三区四区| 亚洲乱码一区二区免费版| 色精品久久人妻99蜜桃| 麻豆一二三区av精品| 国产一区在线观看成人免费| 国模一区二区三区四区视频 | 99国产极品粉嫩在线观看| 亚洲五月婷婷丁香| 亚洲无线观看免费| 18禁国产床啪视频网站| 欧美丝袜亚洲另类 | 亚洲精品在线观看二区| 久久久久国内视频| 欧美日本视频| 日韩欧美三级三区| 特大巨黑吊av在线直播| 成人av一区二区三区在线看| 欧美zozozo另类| 99国产精品99久久久久| 午夜免费观看网址| 欧美3d第一页| 悠悠久久av| 国产精品久久视频播放| 91在线观看av| 精品乱码久久久久久99久播| 欧美色欧美亚洲另类二区| 999久久久精品免费观看国产| 久久中文字幕一级| 中文资源天堂在线| 国产一区二区三区在线臀色熟女| 9191精品国产免费久久| 免费观看人在逋| 老司机深夜福利视频在线观看| 男女之事视频高清在线观看| 色老头精品视频在线观看| 欧美在线一区亚洲| 三级男女做爰猛烈吃奶摸视频| 成人无遮挡网站| 国产精品久久久久久久电影 | 午夜精品一区二区三区免费看| 十八禁网站免费在线| 看免费av毛片| 露出奶头的视频| 搡老熟女国产l中国老女人| 我要搜黄色片| 亚洲av五月六月丁香网| 狂野欧美白嫩少妇大欣赏| 日本免费一区二区三区高清不卡| 在线看三级毛片| 久久久国产欧美日韩av| 人妻久久中文字幕网| 热99在线观看视频| 国产精品一区二区免费欧美| 亚洲熟妇熟女久久| 在线国产一区二区在线| 一夜夜www| 国产精品国产高清国产av| 国产精品久久久久久人妻精品电影| 久久久久国产精品人妻aⅴ院| 国产精品久久久久久人妻精品电影| 亚洲欧美精品综合一区二区三区| 美女cb高潮喷水在线观看 | 美女高潮的动态| 一级a爱片免费观看的视频| 欧美另类亚洲清纯唯美| 精品乱码久久久久久99久播| 成人欧美大片| 香蕉丝袜av| 日韩国内少妇激情av| 免费人成视频x8x8入口观看| 小说图片视频综合网站| 麻豆成人午夜福利视频| av中文乱码字幕在线| 首页视频小说图片口味搜索| 亚洲片人在线观看| 淫妇啪啪啪对白视频| 国产三级黄色录像| 国产免费男女视频| 成人无遮挡网站| 身体一侧抽搐| xxxwww97欧美| 国产高清激情床上av| 国产aⅴ精品一区二区三区波| 欧美黄色片欧美黄色片| 少妇的逼水好多| 午夜亚洲福利在线播放| 一级黄色大片毛片| 好男人在线观看高清免费视频| 伊人久久大香线蕉亚洲五| 欧美3d第一页| 色综合亚洲欧美另类图片| 欧美一级a爱片免费观看看| 欧美日韩中文字幕国产精品一区二区三区| 搞女人的毛片| 一边摸一边抽搐一进一小说| 99久久久亚洲精品蜜臀av| 国内精品久久久久久久电影| 国产不卡一卡二| 国产 一区 欧美 日韩| 久久久久久久久久黄片| 一级毛片高清免费大全| 日本a在线网址| 变态另类成人亚洲欧美熟女| 欧美日韩一级在线毛片| 久久久国产欧美日韩av| 国产高清videossex| 天堂√8在线中文| 麻豆一二三区av精品| 观看美女的网站| 国产黄片美女视频| 国产淫片久久久久久久久 | 国产亚洲精品一区二区www| cao死你这个sao货| 国内精品美女久久久久久| 精品无人区乱码1区二区| 精华霜和精华液先用哪个| 中文字幕最新亚洲高清| 亚洲无线在线观看| 欧美不卡视频在线免费观看| tocl精华| 亚洲自拍偷在线| 国产精品永久免费网站| 国产午夜精品久久久久久| 国产伦一二天堂av在线观看| 亚洲五月天丁香| 12—13女人毛片做爰片一| 日本黄色片子视频| 久久精品国产亚洲av香蕉五月| 看黄色毛片网站| 国产一区在线观看成人免费| 亚洲欧美精品综合久久99| 精品久久久久久,| 欧美丝袜亚洲另类 | 久久久久久久久久黄片| 国产人伦9x9x在线观看| 黄色女人牲交| 欧洲精品卡2卡3卡4卡5卡区| 亚洲片人在线观看| 国内久久婷婷六月综合欲色啪| 最新中文字幕久久久久 | 超碰成人久久| 在线观看美女被高潮喷水网站 | 欧美性猛交黑人性爽| 中文字幕高清在线视频| 久久精品影院6| 美女被艹到高潮喷水动态| 亚洲精品粉嫩美女一区| 真人做人爱边吃奶动态| 淫秽高清视频在线观看| 91在线观看av| 国产精品精品国产色婷婷| 一本综合久久免费| 一个人免费在线观看电影 | 日本a在线网址| 国内精品一区二区在线观看| 91久久精品国产一区二区成人 | 特级一级黄色大片| 亚洲aⅴ乱码一区二区在线播放| 成人高潮视频无遮挡免费网站| 久久亚洲真实| 99视频精品全部免费 在线 | 亚洲国产精品合色在线| 亚洲九九香蕉| 免费在线观看视频国产中文字幕亚洲| 成年女人永久免费观看视频| 久久久久久久久中文| 香蕉丝袜av| 曰老女人黄片| 男女之事视频高清在线观看| 少妇熟女aⅴ在线视频| 99热这里只有是精品50| 亚洲熟妇中文字幕五十中出| 少妇的逼水好多| 中文字幕熟女人妻在线| 熟妇人妻久久中文字幕3abv| 欧美日韩乱码在线| 一本精品99久久精品77| 久久精品国产综合久久久| 久久精品综合一区二区三区| 亚洲欧洲精品一区二区精品久久久| 国产欧美日韩一区二区三| 亚洲第一电影网av| 岛国视频午夜一区免费看| 级片在线观看| 久久婷婷人人爽人人干人人爱| 男女那种视频在线观看| 国产欧美日韩精品亚洲av| 啦啦啦观看免费观看视频高清| 女生性感内裤真人,穿戴方法视频| 色尼玛亚洲综合影院| 亚洲精品美女久久久久99蜜臀| 亚洲欧美日韩高清在线视频| 在线永久观看黄色视频| 久久人妻av系列| 国产午夜精品论理片| 亚洲国产欧美人成| 国产免费男女视频| 国产精品日韩av在线免费观看| 国产亚洲精品av在线| 两个人视频免费观看高清| 99国产精品一区二区三区| 19禁男女啪啪无遮挡网站| 亚洲在线自拍视频| 哪里可以看免费的av片| 最近最新中文字幕大全电影3| www.999成人在线观看| 91在线观看av| 亚洲,欧美精品.| 国产伦精品一区二区三区四那| avwww免费| 三级毛片av免费| 日韩欧美国产一区二区入口| 国产欧美日韩一区二区精品| 999精品在线视频| 国产97色在线日韩免费| 日韩 欧美 亚洲 中文字幕| 亚洲av成人精品一区久久| 69av精品久久久久久| h日本视频在线播放| 这个男人来自地球电影免费观看| 日本黄色视频三级网站网址| 日韩有码中文字幕| 日韩成人在线观看一区二区三区| 人妻久久中文字幕网| 国产高清激情床上av| 黑人欧美特级aaaaaa片| 国产成人福利小说| 久久久久免费精品人妻一区二区| 天天躁日日操中文字幕| 久久精品91蜜桃| 青草久久国产| 午夜福利在线观看免费完整高清在 | 日本免费一区二区三区高清不卡| 国产亚洲欧美98| 色噜噜av男人的天堂激情| 制服人妻中文乱码| 麻豆国产97在线/欧美| 国产成人精品久久二区二区免费| 亚洲九九香蕉| 国产精品99久久99久久久不卡| 在线观看午夜福利视频| 九九在线视频观看精品| 美女cb高潮喷水在线观看 | 这个男人来自地球电影免费观看| 日本黄大片高清| 最新在线观看一区二区三区| 九九久久精品国产亚洲av麻豆 | 老熟妇仑乱视频hdxx| 亚洲精品456在线播放app | 夜夜爽天天搞| 老汉色∧v一级毛片| 成人鲁丝片一二三区免费| 国产又色又爽无遮挡免费看| 国产在线精品亚洲第一网站| 长腿黑丝高跟| 亚洲自偷自拍图片 自拍| 一边摸一边抽搐一进一小说| 亚洲 国产 在线| 老司机午夜福利在线观看视频| 少妇丰满av| 两个人看的免费小视频| 变态另类丝袜制服| 黑人操中国人逼视频| 国产精品 国内视频| 黄色 视频免费看| 少妇的丰满在线观看| 欧美绝顶高潮抽搐喷水| 免费看美女性在线毛片视频|