• <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涩爱 | 一进一出抽搐gif免费好疼| 免费无遮挡裸体视频| 又爽又黄无遮挡网站| 国产成人福利小说| 午夜影院日韩av| 欧美高清成人免费视频www| 色av中文字幕| 国产精品久久电影中文字幕| 嫩草影院入口| 超碰av人人做人人爽久久| 久久久国产成人免费| 久久草成人影院| 亚洲内射少妇av| 少妇丰满av| 国产伦精品一区二区三区视频9| 国产精品,欧美在线| 免费大片18禁| 99久久精品国产国产毛片| 91麻豆av在线| 亚洲图色成人| 大又大粗又爽又黄少妇毛片口| 日韩一区二区视频免费看| 欧美3d第一页| 国产aⅴ精品一区二区三区波| 啦啦啦观看免费观看视频高清| 久久久久久久久中文| 国产69精品久久久久777片| 亚洲av电影不卡..在线观看| 日韩欧美在线二视频| 直男gayav资源| 高清在线国产一区| 黄色丝袜av网址大全| 久久精品久久久久久噜噜老黄 | 久久精品国产清高在天天线| 深夜精品福利| 少妇被粗大猛烈的视频| 国产不卡一卡二| 日韩欧美三级三区| 99九九线精品视频在线观看视频| 国模一区二区三区四区视频| 99热6这里只有精品| 校园春色视频在线观看| 联通29元200g的流量卡| 国内精品一区二区在线观看| 日韩欧美 国产精品| av.在线天堂| 国产免费一级a男人的天堂| 日韩欧美 国产精品| 久久欧美精品欧美久久欧美| 欧美bdsm另类| 他把我摸到了高潮在线观看| 日韩精品有码人妻一区| 成人毛片a级毛片在线播放| 亚洲精品久久国产高清桃花| 欧美一级a爱片免费观看看| 免费人成在线观看视频色| 日日夜夜操网爽| 成人精品一区二区免费| 少妇高潮的动态图| 女的被弄到高潮叫床怎么办 | 日本五十路高清| 国产精品免费一区二区三区在线| 亚洲最大成人av| 成人无遮挡网站| 两人在一起打扑克的视频| 九九爱精品视频在线观看| 熟妇人妻久久中文字幕3abv| 亚洲中文字幕一区二区三区有码在线看| 一a级毛片在线观看| 亚洲自偷自拍三级| АⅤ资源中文在线天堂| 亚洲五月天丁香| 国产中年淑女户外野战色| 午夜精品一区二区三区免费看| 亚洲精华国产精华液的使用体验 | 亚洲精品色激情综合| 国产国拍精品亚洲av在线观看| 丰满人妻一区二区三区视频av| 欧美日韩中文字幕国产精品一区二区三区| 日日干狠狠操夜夜爽| 无人区码免费观看不卡| 三级毛片av免费| 黄色丝袜av网址大全| 成人综合一区亚洲| 免费无遮挡裸体视频| 日日啪夜夜撸| 黄色女人牲交| 国国产精品蜜臀av免费| 欧美一区二区亚洲| 久9热在线精品视频| 麻豆av噜噜一区二区三区| 窝窝影院91人妻| 九九在线视频观看精品| av中文乱码字幕在线| 在线观看一区二区三区| 亚洲熟妇中文字幕五十中出| 欧美黑人巨大hd| 18禁裸乳无遮挡免费网站照片| 国产视频内射| 国产精品一区二区免费欧美| 国内精品一区二区在线观看| 伦理电影大哥的女人| 欧美三级亚洲精品| 国产女主播在线喷水免费视频网站 | 精品午夜福利在线看| 天美传媒精品一区二区| 春色校园在线视频观看| 久久久久性生活片| 国产aⅴ精品一区二区三区波| 男女那种视频在线观看| 日韩,欧美,国产一区二区三区 | 少妇被粗大猛烈的视频| 搡老妇女老女人老熟妇| 99久久精品国产国产毛片| 欧美色欧美亚洲另类二区| 九九爱精品视频在线观看| 99热精品在线国产| 白带黄色成豆腐渣| 免费看光身美女| 成年版毛片免费区| .国产精品久久| 国语自产精品视频在线第100页| 精品一区二区三区视频在线观看免费| 男人舔奶头视频| av在线蜜桃| 女人被狂操c到高潮| 我要搜黄色片| 啦啦啦啦在线视频资源| 赤兔流量卡办理| 夜夜看夜夜爽夜夜摸| 欧美日韩精品成人综合77777| 成年女人永久免费观看视频| av女优亚洲男人天堂| 日日夜夜操网爽| 1000部很黄的大片| 亚洲一级一片aⅴ在线观看| 亚洲国产精品成人综合色| xxxwww97欧美| 亚洲三级黄色毛片| 国产精品久久久久久久电影| 最好的美女福利视频网| 成人三级黄色视频| 日本五十路高清| 草草在线视频免费看| 色哟哟·www| 免费搜索国产男女视频| 亚洲国产色片| 女的被弄到高潮叫床怎么办 | 国产午夜精品久久久久久一区二区三区 | 欧美一级a爱片免费观看看| 不卡视频在线观看欧美| 在线天堂最新版资源| 国产三级中文精品| 国产精品一区二区三区四区久久| 在线观看午夜福利视频| 欧美高清成人免费视频www| 免费一级毛片在线播放高清视频| 亚洲成人中文字幕在线播放| 一本一本综合久久| 亚洲性久久影院| av在线观看视频网站免费| 偷拍熟女少妇极品色| 99久国产av精品| 成人特级黄色片久久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 黄色一级大片看看| 白带黄色成豆腐渣| 3wmmmm亚洲av在线观看| 丰满乱子伦码专区| 97热精品久久久久久| 国内久久婷婷六月综合欲色啪| 看黄色毛片网站| 高清毛片免费观看视频网站| 两个人视频免费观看高清| 亚洲成人久久爱视频| 中文字幕久久专区| 国产高清不卡午夜福利| 在线观看av片永久免费下载| 悠悠久久av| 一本一本综合久久| 一进一出抽搐动态| 久久午夜福利片| 免费观看在线日韩| 日韩欧美 国产精品| 日韩中字成人| 国产精品爽爽va在线观看网站| 毛片一级片免费看久久久久 | 又爽又黄a免费视频| 赤兔流量卡办理| 国产熟女欧美一区二区| 99国产极品粉嫩在线观看| 精品乱码久久久久久99久播| 免费观看人在逋| 国产亚洲欧美98| 亚洲三级黄色毛片| 精品一区二区三区视频在线观看免费| 国产一区二区在线av高清观看| 少妇人妻精品综合一区二区 | 亚洲aⅴ乱码一区二区在线播放| 直男gayav资源| 老熟妇仑乱视频hdxx| 国产一区二区三区在线臀色熟女| 女的被弄到高潮叫床怎么办 | 给我免费播放毛片高清在线观看| 亚洲av成人精品一区久久| 国内少妇人妻偷人精品xxx网站| 日韩欧美三级三区| 简卡轻食公司| 亚洲成av人片在线播放无| 日本免费一区二区三区高清不卡| 国产亚洲91精品色在线| 成熟少妇高潮喷水视频| 午夜福利在线在线| 国产亚洲av嫩草精品影院| 看黄色毛片网站| 国产精品久久久久久精品电影| 热99re8久久精品国产| h日本视频在线播放| 男女边吃奶边做爰视频| 男女那种视频在线观看| 色噜噜av男人的天堂激情| 日韩 亚洲 欧美在线| 免费av不卡在线播放| 九九热线精品视视频播放| 日韩欧美在线二视频| 亚洲av中文av极速乱 | 我的老师免费观看完整版| 亚洲七黄色美女视频| 啦啦啦韩国在线观看视频| 精品国产三级普通话版| 国产单亲对白刺激| 午夜福利在线观看免费完整高清在 | 亚洲精品亚洲一区二区| 久久久久国内视频| 99久久九九国产精品国产免费| 九色成人免费人妻av| 欧洲精品卡2卡3卡4卡5卡区| 免费看a级黄色片| 亚洲精品久久国产高清桃花| 免费无遮挡裸体视频| a级毛片a级免费在线| 老师上课跳d突然被开到最大视频| 亚洲乱码一区二区免费版| 俺也久久电影网| 尾随美女入室| 九色成人免费人妻av| 亚洲人成网站在线播| 成人av一区二区三区在线看| 精品久久久久久久末码| 少妇熟女aⅴ在线视频| 少妇的逼好多水| 一卡2卡三卡四卡精品乱码亚洲| 悠悠久久av| 日本 欧美在线| 日韩在线高清观看一区二区三区 | 亚洲avbb在线观看| 少妇人妻一区二区三区视频| xxxwww97欧美| 97碰自拍视频| 国产精品一区www在线观看 | 久99久视频精品免费| av福利片在线观看| 69av精品久久久久久| 午夜日韩欧美国产| 精品久久国产蜜桃| 免费观看精品视频网站| 亚州av有码| 国产精品久久电影中文字幕| www日本黄色视频网| 别揉我奶头~嗯~啊~动态视频| 小说图片视频综合网站| 毛片一级片免费看久久久久 | 高清日韩中文字幕在线| 中文字幕免费在线视频6| 一个人看的www免费观看视频| 丰满人妻一区二区三区视频av| 亚洲一区高清亚洲精品| 12—13女人毛片做爰片一| 欧美xxxx黑人xx丫x性爽| 高清在线国产一区| 日韩精品青青久久久久久| 午夜福利视频1000在线观看| 亚洲天堂国产精品一区在线| 此物有八面人人有两片| 日韩 亚洲 欧美在线| 午夜免费成人在线视频| videossex国产| 搞女人的毛片| a级毛片a级免费在线| 国产v大片淫在线免费观看| 男女啪啪激烈高潮av片| 1024手机看黄色片| 久久精品91蜜桃| 男女视频在线观看网站免费| netflix在线观看网站| 天堂√8在线中文| 精品人妻一区二区三区麻豆 | 亚洲一级一片aⅴ在线观看| 欧美极品一区二区三区四区| 亚洲精品在线观看二区| 午夜激情欧美在线| 日本黄大片高清| 中国美女看黄片| 黄色配什么色好看| 亚洲av第一区精品v没综合| 人人妻人人看人人澡| 亚洲一级一片aⅴ在线观看| 欧美高清性xxxxhd video| 最好的美女福利视频网| 日本黄色片子视频| 国产精品福利在线免费观看| 国产成人aa在线观看| 精品一区二区三区视频在线观看免费| 男人的好看免费观看在线视频| 婷婷亚洲欧美| 免费人成视频x8x8入口观看| 亚洲成av人片在线播放无| 欧美黑人欧美精品刺激| 最近最新中文字幕大全电影3| 国产伦精品一区二区三区视频9| 久久精品国产亚洲网站| 婷婷精品国产亚洲av在线| 春色校园在线视频观看| 日韩中文字幕欧美一区二区| 久久午夜福利片| 久久久久国产精品人妻aⅴ院| 日韩欧美免费精品| 国产伦在线观看视频一区| 综合色av麻豆| 日日夜夜操网爽| 欧美黑人欧美精品刺激| 午夜精品一区二区三区免费看| 91麻豆精品激情在线观看国产| 国产精品久久电影中文字幕| 久久亚洲精品不卡| 亚洲av不卡在线观看| 欧美又色又爽又黄视频| 99热这里只有是精品在线观看| 日韩一本色道免费dvd| 欧美色视频一区免费| 亚洲中文字幕日韩| 18禁裸乳无遮挡免费网站照片| 婷婷精品国产亚洲av| 亚洲欧美日韩卡通动漫| 亚洲色图av天堂| 99国产精品一区二区蜜桃av| 免费看光身美女| 国产黄色小视频在线观看| 91狼人影院| 久久久精品欧美日韩精品| 精品久久久久久久人妻蜜臀av| 亚洲内射少妇av| 淫妇啪啪啪对白视频| 亚洲综合色惰| 亚洲在线自拍视频| 国产亚洲91精品色在线| 国产欧美日韩精品一区二区| 欧美一区二区国产精品久久精品| 窝窝影院91人妻| 哪里可以看免费的av片| 免费人成在线观看视频色| 久久人人爽人人爽人人片va| 日韩一区二区视频免费看| 久久精品影院6| 免费观看精品视频网站| 国产精品三级大全| 小蜜桃在线观看免费完整版高清| 舔av片在线| 尾随美女入室| 91在线精品国自产拍蜜月| 内射极品少妇av片p| 免费大片18禁| 看免费成人av毛片| 嫩草影院入口| av.在线天堂| 国产亚洲av嫩草精品影院| 国产午夜福利久久久久久| 免费观看的影片在线观看| 999久久久精品免费观看国产| 禁无遮挡网站| 91久久精品电影网| 国内精品久久久久久久电影| 美女大奶头视频| 国产v大片淫在线免费观看| 精品人妻视频免费看| 国产淫片久久久久久久久| 国产一区二区三区视频了| 在现免费观看毛片| 日韩一本色道免费dvd| bbb黄色大片| 大型黄色视频在线免费观看| 国产日本99.免费观看| 午夜福利18| 午夜激情福利司机影院| 亚洲真实伦在线观看| 国产69精品久久久久777片| 亚洲自拍偷在线| 国产v大片淫在线免费观看| 3wmmmm亚洲av在线观看| 在线播放国产精品三级| 国产精品人妻久久久影院| 男女那种视频在线观看| 麻豆国产97在线/欧美| 欧美最黄视频在线播放免费| 九色成人免费人妻av| 麻豆一二三区av精品| 麻豆国产97在线/欧美| 赤兔流量卡办理| 亚洲性夜色夜夜综合| 国产亚洲精品av在线| 观看美女的网站| 国产精品乱码一区二三区的特点| 夜夜爽天天搞| 国产极品精品免费视频能看的| 男人和女人高潮做爰伦理| 日本-黄色视频高清免费观看| 精品欧美国产一区二区三| 1000部很黄的大片| 国产视频内射| 少妇丰满av| 国产在视频线在精品| 在线免费观看的www视频| 啦啦啦韩国在线观看视频| 久久久久久大精品| 伊人久久精品亚洲午夜| 欧美成人a在线观看| 中文亚洲av片在线观看爽| 午夜精品在线福利| 波多野结衣高清无吗| 色视频www国产| 3wmmmm亚洲av在线观看| 又爽又黄a免费视频| 国产精品人妻久久久影院| 看片在线看免费视频| 床上黄色一级片| 亚洲成av人片在线播放无| 一夜夜www| 一个人看的www免费观看视频| 久久午夜福利片| 亚洲成人免费电影在线观看| 51国产日韩欧美| 亚洲乱码一区二区免费版| 欧美成人免费av一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 尤物成人国产欧美一区二区三区| 亚洲av成人av| 91狼人影院| 日韩av在线大香蕉| 麻豆成人av在线观看| 男女视频在线观看网站免费| 麻豆久久精品国产亚洲av| 无人区码免费观看不卡| 美女大奶头视频| 高清在线国产一区| 国产黄a三级三级三级人| 小蜜桃在线观看免费完整版高清| 舔av片在线| 日韩精品中文字幕看吧| 国产亚洲av嫩草精品影院| 亚洲精品亚洲一区二区| 中出人妻视频一区二区| 久久这里只有精品中国| 午夜久久久久精精品| 国产精品一及| 校园春色视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 国产精品久久久久久久电影| 国产真实乱freesex| 99久久久亚洲精品蜜臀av| 色精品久久人妻99蜜桃| 看免费成人av毛片| 日本 欧美在线| 麻豆成人av在线观看| 男女视频在线观看网站免费| 亚洲人成网站在线播放欧美日韩| 久久九九热精品免费| 日韩中文字幕欧美一区二区| 国产成人aa在线观看| 不卡一级毛片| 少妇人妻精品综合一区二区 | 岛国在线免费视频观看| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲av嫩草精品影院| 日韩人妻高清精品专区| 久久草成人影院| 日本与韩国留学比较| 免费看日本二区| 美女大奶头视频| 亚洲图色成人| 欧美黑人欧美精品刺激| 人人妻人人看人人澡| 校园人妻丝袜中文字幕| .国产精品久久| 国内精品宾馆在线| 男女啪啪激烈高潮av片| 如何舔出高潮| 搡女人真爽免费视频火全软件 | 国产精品伦人一区二区| 婷婷色综合大香蕉| 老师上课跳d突然被开到最大视频| 变态另类丝袜制服| 波多野结衣巨乳人妻| 在线观看美女被高潮喷水网站| 欧美一级a爱片免费观看看| 欧美日韩中文字幕国产精品一区二区三区| 欧美激情国产日韩精品一区| 99久久无色码亚洲精品果冻| 日本黄色视频三级网站网址| 国产成人av教育| 少妇被粗大猛烈的视频| 日本一本二区三区精品| 有码 亚洲区| 亚洲va日本ⅴa欧美va伊人久久| 久久久久九九精品影院| av中文乱码字幕在线| 国产在视频线在精品| 亚洲最大成人av| 在线观看66精品国产| 三级毛片av免费| 亚洲专区中文字幕在线| 老司机福利观看| 男女做爰动态图高潮gif福利片| 成人亚洲精品av一区二区| 久久中文看片网| 精品国内亚洲2022精品成人| 欧美性感艳星| 少妇高潮的动态图| 久久久久国内视频| 在线观看免费视频日本深夜| 桃红色精品国产亚洲av| 午夜福利在线在线| 精品久久久久久成人av| 亚洲成人中文字幕在线播放| 老女人水多毛片| 久久久久久久午夜电影| 精品国内亚洲2022精品成人| 超碰av人人做人人爽久久| 亚洲专区国产一区二区| 很黄的视频免费| 一进一出好大好爽视频| 在线观看舔阴道视频| 高清毛片免费观看视频网站| 亚洲美女黄片视频| 日韩欧美国产在线观看| 国产精品精品国产色婷婷| 麻豆成人av在线观看| 99热这里只有是精品在线观看| 97人妻精品一区二区三区麻豆| 欧美极品一区二区三区四区| 丰满人妻一区二区三区视频av| 色5月婷婷丁香| 国内精品宾馆在线| av国产免费在线观看| 中文字幕久久专区| 欧美在线一区亚洲| 国产乱人视频| 中文字幕人妻熟人妻熟丝袜美| 国产一区二区激情短视频| 免费看光身美女| 精品久久久久久久末码| 国产三级中文精品| 国产 一区 欧美 日韩| 听说在线观看完整版免费高清| 动漫黄色视频在线观看| 免费看日本二区| 哪里可以看免费的av片| 网址你懂的国产日韩在线| 在线免费观看不下载黄p国产 | a级一级毛片免费在线观看| 亚洲精品影视一区二区三区av| 神马国产精品三级电影在线观看| 老女人水多毛片| 校园人妻丝袜中文字幕| 国产高清视频在线播放一区| 久久精品国产鲁丝片午夜精品 | 又紧又爽又黄一区二区| 最新在线观看一区二区三区| 国产精品一区二区三区四区久久| 又黄又爽又免费观看的视频| 99精品久久久久人妻精品| 男女做爰动态图高潮gif福利片| 国产精品不卡视频一区二区| 禁无遮挡网站| 男女下面进入的视频免费午夜| 日日摸夜夜添夜夜添小说| 国模一区二区三区四区视频| 国产 一区精品| 亚洲在线自拍视频| 男女视频在线观看网站免费| 自拍偷自拍亚洲精品老妇| 国产成人a区在线观看| 色播亚洲综合网| av天堂在线播放| 精品久久久噜噜| 色播亚洲综合网| av中文乱码字幕在线| 免费观看精品视频网站| 91久久精品国产一区二区成人| 午夜免费成人在线视频| 成人综合一区亚洲| 日本 欧美在线| 国产真实伦视频高清在线观看 | 欧美成人性av电影在线观看| 在线观看av片永久免费下载| 国产 一区 欧美 日韩| 一区二区三区四区激情视频 | 日韩欧美 国产精品| 久久精品国产鲁丝片午夜精品 | 欧美日韩黄片免|