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

    国产精品综合久久久久久久免费 | 国内精品久久久久精免费| 欧美激情 高清一区二区三区| 亚洲全国av大片| 91麻豆精品激情在线观看国产| 首页视频小说图片口味搜索| 黄色视频,在线免费观看| a在线观看视频网站| 成人国产一区最新在线观看| 丁香六月欧美| 国产国语露脸激情在线看| 啦啦啦 在线观看视频| 久久欧美精品欧美久久欧美| 看片在线看免费视频| 亚洲人成伊人成综合网2020| 制服丝袜大香蕉在线| 亚洲精品久久成人aⅴ小说| 亚洲欧美日韩另类电影网站| 久久人妻av系列| 久久午夜亚洲精品久久| 午夜免费成人在线视频| 老司机深夜福利视频在线观看| 欧美最黄视频在线播放免费| 我的亚洲天堂| 久久草成人影院| 99久久久亚洲精品蜜臀av| 国产免费男女视频| 亚洲成国产人片在线观看| 国产精品久久久人人做人人爽| 女人被躁到高潮嗷嗷叫费观| 午夜福利成人在线免费观看| 在线观看午夜福利视频| 国产精品 国内视频| 亚洲人成77777在线视频| 亚洲国产毛片av蜜桃av| 亚洲专区国产一区二区| 亚洲欧洲精品一区二区精品久久久| 欧美乱妇无乱码| 国产精品爽爽va在线观看网站 | 操出白浆在线播放| 这个男人来自地球电影免费观看| 97碰自拍视频| 亚洲三区欧美一区| 九色亚洲精品在线播放| 不卡av一区二区三区| 国产精品乱码一区二三区的特点 | 久久久精品欧美日韩精品| 日本精品一区二区三区蜜桃| 欧美日韩亚洲国产一区二区在线观看| 视频在线观看一区二区三区| 18禁美女被吸乳视频| 中文字幕av电影在线播放| 成人免费观看视频高清| 久久国产精品男人的天堂亚洲| 成人特级黄色片久久久久久久| 99香蕉大伊视频| 9191精品国产免费久久| 一级毛片精品| 国产成人精品久久二区二区免费| 男男h啪啪无遮挡| 欧美最黄视频在线播放免费| 国产aⅴ精品一区二区三区波| av视频在线观看入口| 日本五十路高清| 免费看a级黄色片| 午夜精品国产一区二区电影| 亚洲全国av大片| 91在线观看av| 免费av毛片视频| 中文亚洲av片在线观看爽| 我的亚洲天堂| 亚洲成国产人片在线观看| 丝袜人妻中文字幕| or卡值多少钱| 精品一品国产午夜福利视频| 在线播放国产精品三级| 级片在线观看| 久久久久久久久免费视频了| 欧美日韩福利视频一区二区| 啪啪无遮挡十八禁网站| 大陆偷拍与自拍| 此物有八面人人有两片| 老司机深夜福利视频在线观看| 国内精品久久久久精免费| 最新美女视频免费是黄的| 久久久水蜜桃国产精品网| 久久青草综合色| 亚洲一区高清亚洲精品| 成人手机av| 99国产综合亚洲精品| 亚洲片人在线观看| 久久人人精品亚洲av| 久热爱精品视频在线9| 国产精品国产高清国产av| 久久久国产成人免费| 久久香蕉国产精品| 天天添夜夜摸| 国产精品九九99| 俄罗斯特黄特色一大片| av视频在线观看入口| 母亲3免费完整高清在线观看| 欧美一区二区精品小视频在线| 国产成人啪精品午夜网站| 欧美成人免费av一区二区三区| 亚洲免费av在线视频| 免费久久久久久久精品成人欧美视频| 亚洲五月天丁香| 日本免费a在线| 欧美成人性av电影在线观看| 国产免费av片在线观看野外av| 日本欧美视频一区| 欧美午夜高清在线| 午夜福利,免费看| 99国产精品免费福利视频| 色在线成人网| 大型黄色视频在线免费观看| 亚洲人成网站在线播放欧美日韩| 国产伦一二天堂av在线观看| 黑人欧美特级aaaaaa片| 欧美成人免费av一区二区三区| 精品熟女少妇八av免费久了| av网站免费在线观看视频| 99国产综合亚洲精品| 久久久水蜜桃国产精品网| 久久精品国产清高在天天线| 亚洲少妇的诱惑av| 黄色 视频免费看| 精品少妇一区二区三区视频日本电影| 一级毛片高清免费大全| 国产欧美日韩精品亚洲av| 女同久久另类99精品国产91| 国产成人精品久久二区二区免费| 久久久久久久久免费视频了| 亚洲av电影不卡..在线观看| av片东京热男人的天堂| 中文字幕人成人乱码亚洲影| 国产精品久久久av美女十八| 国产欧美日韩一区二区三| 久久久水蜜桃国产精品网| 中文字幕另类日韩欧美亚洲嫩草| 麻豆av在线久日| 日本一区二区免费在线视频| 熟女少妇亚洲综合色aaa.| 丝袜人妻中文字幕| 日本精品一区二区三区蜜桃| 久久久久久久午夜电影| 精品国产乱码久久久久久男人| 露出奶头的视频| 天堂动漫精品| 免费av毛片视频| 国产av一区二区精品久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人系列免费观看| 亚洲第一青青草原| 免费高清在线观看日韩| 88av欧美| 精品乱码久久久久久99久播| 黄色女人牲交| 一区二区日韩欧美中文字幕| av天堂在线播放| 天天躁夜夜躁狠狠躁躁| 国产亚洲精品久久久久久毛片| 亚洲色图av天堂| 欧美性长视频在线观看| 一本综合久久免费| 女人爽到高潮嗷嗷叫在线视频| 亚洲中文日韩欧美视频| 91精品国产国语对白视频| 亚洲精品国产色婷婷电影| 伦理电影免费视频| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久电影中文字幕| 老司机福利观看| x7x7x7水蜜桃| 免费不卡黄色视频| 女人精品久久久久毛片| 亚洲成a人片在线一区二区| www.精华液| 热99re8久久精品国产| 夜夜躁狠狠躁天天躁| 欧美日韩瑟瑟在线播放| 午夜精品国产一区二区电影| 亚洲精品中文字幕在线视频| 大型黄色视频在线免费观看| av在线播放免费不卡| 亚洲av成人av| 国内毛片毛片毛片毛片毛片| 久久香蕉国产精品| 99国产精品一区二区蜜桃av| 这个男人来自地球电影免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲在线自拍视频| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩高清在线视频| 国产亚洲精品一区二区www| 女性被躁到高潮视频| 精品国产亚洲在线| 国产午夜福利久久久久久| 亚洲国产精品999在线| 在线免费观看的www视频| 免费看美女性在线毛片视频| 亚洲人成伊人成综合网2020| 亚洲精品国产色婷婷电影| 久久国产亚洲av麻豆专区| 一级毛片高清免费大全| 国产精品98久久久久久宅男小说| 国产欧美日韩一区二区三区在线| 久久精品国产亚洲av香蕉五月| 欧美激情高清一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 伦理电影免费视频| 精品久久久久久久久久免费视频| 色老头精品视频在线观看| 日本免费一区二区三区高清不卡 | 久久中文字幕人妻熟女| www国产在线视频色| 窝窝影院91人妻| 777久久人妻少妇嫩草av网站| 神马国产精品三级电影在线观看 | 久久性视频一级片| 久久这里只有精品19| 欧洲精品卡2卡3卡4卡5卡区| 欧美另类亚洲清纯唯美| 亚洲在线自拍视频| 日日爽夜夜爽网站| 亚洲国产精品合色在线| 动漫黄色视频在线观看| 国产精品久久久av美女十八| 在线观看免费视频网站a站| 琪琪午夜伦伦电影理论片6080| 亚洲七黄色美女视频| 国产精品久久视频播放| 久久午夜亚洲精品久久| 日韩一卡2卡3卡4卡2021年| 人妻久久中文字幕网| 高清在线国产一区| 亚洲自偷自拍图片 自拍| 色哟哟哟哟哟哟| 国产精品久久久人人做人人爽| 美女大奶头视频| 91精品三级在线观看| av网站免费在线观看视频| 一区二区三区精品91| 大码成人一级视频| 精品久久久久久久人妻蜜臀av | 午夜久久久久精精品| 男男h啪啪无遮挡| 国产一区二区三区视频了| 精品人妻1区二区| 97人妻精品一区二区三区麻豆 | 国产成人精品久久二区二区免费| 黑人巨大精品欧美一区二区mp4| 91老司机精品| 免费无遮挡裸体视频| 男人舔女人的私密视频| 久久精品国产清高在天天线| 啦啦啦 在线观看视频| 免费人成视频x8x8入口观看| 国产成人av教育| 精品日产1卡2卡| 免费观看精品视频网站| 久久香蕉精品热| 亚洲片人在线观看| 欧美成人性av电影在线观看| tocl精华| 亚洲欧美精品综合久久99| 色综合婷婷激情| 99久久久亚洲精品蜜臀av| 男女下面进入的视频免费午夜 | 麻豆一二三区av精品| 女性被躁到高潮视频| 麻豆成人av在线观看| 欧美黑人精品巨大| 免费高清在线观看日韩| 色播亚洲综合网| 国产亚洲欧美98| 色综合婷婷激情| 中文亚洲av片在线观看爽| 精品人妻在线不人妻| 深夜精品福利| 大型av网站在线播放| 麻豆国产av国片精品| 可以免费在线观看a视频的电影网站| 亚洲欧美一区二区三区黑人| 男人操女人黄网站| 99久久综合精品五月天人人| 搡老熟女国产l中国老女人| av中文乱码字幕在线| 精品久久久久久久毛片微露脸| 国产亚洲av嫩草精品影院| 亚洲男人的天堂狠狠| 又黄又爽又免费观看的视频| 婷婷六月久久综合丁香| 久久天堂一区二区三区四区| 亚洲美女黄片视频| 欧美午夜高清在线| 欧美在线一区亚洲| 国产成人av激情在线播放| 日韩免费av在线播放| 久久香蕉国产精品| 两人在一起打扑克的视频| 91在线观看av| 黑人巨大精品欧美一区二区蜜桃| 18禁美女被吸乳视频| 午夜成年电影在线免费观看| av电影中文网址| 亚洲欧美日韩无卡精品| 天天添夜夜摸| 欧美日韩一级在线毛片| 亚洲国产欧美一区二区综合| 操美女的视频在线观看| 极品人妻少妇av视频| 在线观看免费午夜福利视频| 热99re8久久精品国产| 久久青草综合色| 999久久久精品免费观看国产| 久久精品国产清高在天天线| 美女午夜性视频免费| 淫秽高清视频在线观看| tocl精华| 伦理电影免费视频| 亚洲va日本ⅴa欧美va伊人久久| 成人av一区二区三区在线看| 色综合亚洲欧美另类图片| 青草久久国产| 伦理电影免费视频| 黑人巨大精品欧美一区二区mp4| 在线观看一区二区三区| 日本免费a在线| 国产精品 欧美亚洲| 亚洲成人免费电影在线观看| 国产一卡二卡三卡精品| 久久精品91无色码中文字幕| 脱女人内裤的视频| 在线观看66精品国产| 999久久久精品免费观看国产| 亚洲成人国产一区在线观看| 欧美成人性av电影在线观看| 一区在线观看完整版| 欧美国产日韩亚洲一区| 欧美日韩瑟瑟在线播放| 老司机靠b影院| 欧美日本中文国产一区发布| tocl精华| 色综合婷婷激情| 最好的美女福利视频网| 最近最新中文字幕大全免费视频| 欧美国产日韩亚洲一区| 成人特级黄色片久久久久久久| 亚洲精品久久成人aⅴ小说| 天堂√8在线中文| 精品国产乱码久久久久久男人| 女人精品久久久久毛片| 国产熟女xx| 国产亚洲精品第一综合不卡| 亚洲性夜色夜夜综合| 中文字幕色久视频| 黑人巨大精品欧美一区二区mp4| 欧美激情 高清一区二区三区| 精品欧美国产一区二区三| 亚洲中文日韩欧美视频| 久久国产精品男人的天堂亚洲| 免费人成视频x8x8入口观看| 欧美久久黑人一区二区| 成人国产一区最新在线观看| 琪琪午夜伦伦电影理论片6080| 欧美丝袜亚洲另类 | 波多野结衣高清无吗| 制服丝袜大香蕉在线| 欧洲精品卡2卡3卡4卡5卡区| 女同久久另类99精品国产91| 老汉色av国产亚洲站长工具| 精品一区二区三区四区五区乱码| 成人永久免费在线观看视频| 免费少妇av软件| 中国美女看黄片| 久久久精品国产亚洲av高清涩受| 久久青草综合色| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久久久亚洲av鲁大| 天堂动漫精品| 老鸭窝网址在线观看| 国产国语露脸激情在线看| 久久久久精品国产欧美久久久| 熟妇人妻久久中文字幕3abv| 精品电影一区二区在线| 宅男免费午夜| 亚洲精品久久国产高清桃花| 一级毛片精品| 国产亚洲精品第一综合不卡| 亚洲av熟女| 亚洲男人的天堂狠狠| 夜夜看夜夜爽夜夜摸| 女性生殖器流出的白浆| 无人区码免费观看不卡| 色婷婷久久久亚洲欧美| 国产私拍福利视频在线观看| 一级作爱视频免费观看| 国产亚洲精品综合一区在线观看 | 可以在线观看的亚洲视频| 搞女人的毛片| 久久中文看片网| 色精品久久人妻99蜜桃| 日本vs欧美在线观看视频| 怎么达到女性高潮| 久久婷婷成人综合色麻豆| 一a级毛片在线观看| 亚洲精品中文字幕一二三四区| 日本vs欧美在线观看视频| 99国产精品99久久久久| 日本黄色视频三级网站网址| av中文乱码字幕在线| 一a级毛片在线观看| 亚洲无线在线观看| 18禁裸乳无遮挡免费网站照片 | 国产成+人综合+亚洲专区| 国产精品1区2区在线观看.| 99久久久亚洲精品蜜臀av| 国产精品久久久av美女十八| 中文字幕人妻丝袜一区二区| 妹子高潮喷水视频| 男人舔女人下体高潮全视频| www.www免费av| 久久青草综合色| 人人妻,人人澡人人爽秒播| 免费观看人在逋| 91大片在线观看| 久久精品91蜜桃| 国产精品精品国产色婷婷| 国产精品久久电影中文字幕| 老熟妇仑乱视频hdxx| 午夜福利免费观看在线| 欧美在线一区亚洲| 欧美中文综合在线视频| 精品一区二区三区视频在线观看免费| 19禁男女啪啪无遮挡网站| 大码成人一级视频| 久久精品91蜜桃| 成人免费观看视频高清| 亚洲欧美精品综合久久99| 亚洲av五月六月丁香网| 九色国产91popny在线| 国产精品久久久久久亚洲av鲁大| 色播亚洲综合网| 日日爽夜夜爽网站| 在线视频色国产色| 青草久久国产| 午夜福利影视在线免费观看| 91麻豆精品激情在线观看国产| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美乱色亚洲激情| 黄色女人牲交| av福利片在线| 亚洲五月色婷婷综合| 亚洲av电影在线进入| 亚洲欧美激情综合另类| 国产熟女xx| 国产成人系列免费观看| 精品免费久久久久久久清纯| 久久精品国产99精品国产亚洲性色 | 国产又爽黄色视频| 亚洲av第一区精品v没综合| 中出人妻视频一区二区| 欧美精品亚洲一区二区| 午夜福利一区二区在线看| 国产精品美女特级片免费视频播放器 | 丝袜美腿诱惑在线| 黑人欧美特级aaaaaa片| 亚洲国产精品999在线| 大香蕉久久成人网| 少妇熟女aⅴ在线视频| 美女高潮到喷水免费观看| 欧美日韩亚洲国产一区二区在线观看| 国产高清视频在线播放一区| 欧美在线一区亚洲| 一区二区三区国产精品乱码| 在线永久观看黄色视频| 精品欧美一区二区三区在线| 国产亚洲av高清不卡| 美女高潮到喷水免费观看| 操美女的视频在线观看| 亚洲欧美日韩另类电影网站| 久久中文看片网| 首页视频小说图片口味搜索| 亚洲精品av麻豆狂野| 人人妻,人人澡人人爽秒播| 两个人视频免费观看高清| 两人在一起打扑克的视频| 欧美精品亚洲一区二区| 国产精品一区二区三区四区久久 | 国产精品久久久av美女十八| 久久午夜亚洲精品久久| 日韩国内少妇激情av| 亚洲,欧美精品.| 首页视频小说图片口味搜索| 国产99白浆流出| 免费女性裸体啪啪无遮挡网站| 国产高清videossex| 18禁观看日本| 日韩欧美免费精品| 人人澡人人妻人| 巨乳人妻的诱惑在线观看| 国产区一区二久久| 国产欧美日韩一区二区三区在线| 国产色视频综合| 侵犯人妻中文字幕一二三四区| 国产成人精品久久二区二区91| 免费在线观看视频国产中文字幕亚洲| 国产精品久久久久久精品电影 | 色综合站精品国产| 亚洲国产中文字幕在线视频| 日本免费一区二区三区高清不卡 | 日日夜夜操网爽| 精品人妻在线不人妻| 中文字幕色久视频| 黑人巨大精品欧美一区二区蜜桃| 国产亚洲av嫩草精品影院| 午夜激情av网站| 亚洲欧美一区二区三区黑人| 亚洲av片天天在线观看| 一二三四社区在线视频社区8| 亚洲第一电影网av| 搡老岳熟女国产| 女人精品久久久久毛片| 999久久久国产精品视频| 制服诱惑二区| 亚洲午夜精品一区,二区,三区| 日韩视频一区二区在线观看| 久久人妻av系列| avwww免费| aaaaa片日本免费| 亚洲中文字幕一区二区三区有码在线看 | 免费观看人在逋| 伦理电影免费视频| 97人妻精品一区二区三区麻豆 | 亚洲av成人av| 欧美国产日韩亚洲一区| 久久狼人影院| 亚洲av电影不卡..在线观看| 女人高潮潮喷娇喘18禁视频| 老司机午夜十八禁免费视频| 一进一出抽搐gif免费好疼| 深夜精品福利| 久久人妻福利社区极品人妻图片| 国产亚洲av嫩草精品影院| 久久久久国产一级毛片高清牌| 女人高潮潮喷娇喘18禁视频| 99久久国产精品久久久| 纯流量卡能插随身wifi吗| 欧美乱妇无乱码| 亚洲 欧美 日韩 在线 免费| a在线观看视频网站| 亚洲专区国产一区二区| 亚洲精华国产精华精| 欧美激情久久久久久爽电影 | www.熟女人妻精品国产| 在线观看免费视频网站a站| 午夜福利成人在线免费观看| 国产成人系列免费观看| 欧美黄色淫秽网站| av中文乱码字幕在线| 成人国语在线视频| 男女午夜视频在线观看| 国内精品久久久久久久电影| 久久人人97超碰香蕉20202| 操出白浆在线播放| 亚洲三区欧美一区| 如日韩欧美国产精品一区二区三区| 精品电影一区二区在线| 很黄的视频免费| 国产亚洲av高清不卡| 十八禁网站免费在线| 99久久精品国产亚洲精品| 亚洲中文字幕日韩| 电影成人av| 精品午夜福利视频在线观看一区| 免费在线观看影片大全网站| 欧美成狂野欧美在线观看| 夜夜看夜夜爽夜夜摸| 99riav亚洲国产免费| 一级黄色大片毛片| 久久精品91无色码中文字幕| 老司机在亚洲福利影院| av视频在线观看入口| 又黄又爽又免费观看的视频| av视频在线观看入口| 欧美激情久久久久久爽电影 | 午夜免费鲁丝| 免费高清视频大片| 国产激情欧美一区二区| 国产三级黄色录像| 热re99久久国产66热| 90打野战视频偷拍视频| 色播在线永久视频| 亚洲性夜色夜夜综合| 国产亚洲精品一区二区www| 亚洲狠狠婷婷综合久久图片| 国产精品野战在线观看| 国产精品九九99| 国产1区2区3区精品| 深夜精品福利| 啦啦啦韩国在线观看视频| 亚洲精品粉嫩美女一区| 国产99久久九九免费精品| 日韩三级视频一区二区三区| 啦啦啦韩国在线观看视频| 琪琪午夜伦伦电影理论片6080| 色播在线永久视频| 欧美激情高清一区二区三区| 欧美国产精品va在线观看不卡| 亚洲国产精品成人综合色| 侵犯人妻中文字幕一二三四区| 桃红色精品国产亚洲av| 12—13女人毛片做爰片一| 99久久国产精品久久久|