Next Article in Journal
Influence of Pre-Radiation and Photo-Bleaching on the Yb-Doped Fiber Laser Radiated with Gamma-ray
Next Article in Special Issue
Switching Control of Closed-Loop Supply Chain Systems with Markov Jump Parameters
Previous Article in Journal
Study on Load–Slip Curve of a PBL Shear Key at a Steel–Concrete Composite Joint
Previous Article in Special Issue
A Cosine-Similarity-Based Deconvolution Method for Analyzing Data-Independent Acquisition Mass Spectrometry Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Three-Dimensional Scanning Path Planning of Casing Parts Based on Industrial Robot

1
Aviation Manufacturing Technology Digital Defense Key Discipline Laboratory, Shenyang Aerospace University, Shenyang 110136, China
2
School of Mechanical and Electrical Engineering, Shenyang Aerospace University, Shenyang 110136, China
3
School of Artificial Intelligence, Shenyang Aerospace University, Shenyang 110136, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2023, 13(10), 6162; https://doi.org/10.3390/app13106162
Submission received: 5 April 2023 / Revised: 24 April 2023 / Accepted: 25 April 2023 / Published: 17 May 2023

Abstract

:

Featured Application

This paper analyzes the automatic scanning path planning method and the scanning quality influence rule for casing parts.

Abstract

According to the requirements for the rapid scanning and measurement of the geometric shape during the process of chemical milling of an aviation engine casing, the scanning path of the casing is planned and studied. This paper introduces the principle and method of the tracking scanner and automatic measuring system and analyzes the scanning area range, approach distance, and wide angle of the field. The casing process is modeled by applying part of the machine, obtaining a series of scanning path point and synthesizing the scanning trajectory. On this basis, the entire scanning process is divided into two alternating actions: scanning measurement and posture adjustment, and the mathematical model of the annular scanning path on the outer surface of the casing part is obtained. Adjusting the scan height was used to solve the repeated scan area problem, and the results show that the adjustment method effectively shortened the scan path’s length and time. The simulation method verifies the planned finite ring-scanning path, which verifies the correctness and feasibility of the mathematical model. Through the automatic scanning reconstruction process test, the reconstruction rate of the ring scanning trajectory reaches 85%, which is 80% higher than the manual detection efficiency.

1. Introduction

With the continuous development of aero-engine manufacturing technology, under the premise of ensuring the high stability of aero-engines [1], the lightweight design and manufacture of aero-engine casings has become one of the most critical and effective technical means to improve the thrust–weight ratio of aero-engines. The casing is the basic skeleton of an aviation engine, comprising not only the central parts but also the main bearing parts of an aviation engine. The engine casing is a critical component of an aviation engine. Under the harsh conditions of high temperature and high pressure for a long time, it bears the high-temperature load and thrust load of the combustion chamber, and its manufacturing precision and performance requirements are incredibly high.
At present, the lightweight problem of aviation engine casing parts is solved mainly by chemical milling after traditional metal cutting [2]. Chemical milling removes the excess metal on the surface of the features, forming reinforcement and structural convex platform characteristics, reducing the weight of the self-structure without affecting the stress state of the casing elements. The chemically milled casing has characteristics of a thin wall, complex reinforced reinforcement structure, high precision, and weak rigid conical shell parts [3]. In the process of chemically milling the casing parts, it is necessary to measure the casing many times during the chemical milling, to evaluate the machining accuracy of the casing parts and the effectiveness of the chemical milling process. The high concentration of acid corrosion liquid remaining on the surface of the pieces will also constantly corrode the parts during the measurement process [4], thus influencing wall thickness. Currently, this enterprise still uses the traditional ultrasonic thickness measurement method to measure the thickness of the critical position of the casing. In addition to the harsh operating environment, there is low detection efficiency and poor repetition accuracy in the results. The accurate measurement for the chemical milling process has become one of the greatest technological limitations in optimizing the chemical milling process and evaluating its lightweight-forming quality.
The surface shape of significant components is complex, and the manufacturing technology is complicated, so the morphology and size measurement has become the bottleneck that restricts the development and technological advancement of large equipment. Under field conditions, measuring the complex surface of significant components must meet the requirements for large scale, high point-cloud density, high precision, and high efficiency. Given the importance of the precision measurement of complex surfaces during the manufacturing of significant components in large equipment, the integrated automatic measurement method applicable to the industrial area has been studied [5]. Realizing the size of complex components requires both measurement accuracy and efficiency. The research background is described in Figure 1.
As early as the 1980s, Germany began studying 3-D measurement technology in surface-structured light. In 1985, Dr. Breckmann introduced phase shift interferometry to three-dimensional morphology measurement, establishing a new three-dimensional morphology measurement technology: the phase-measurement contour technique (PMP) [6,7]. Additionally, Dr. Steinbichler and Professor Reinhold Ritter at Germany’s Braunschweig Technical University introduced the Steinbichler GmbH COMET 5-3D measurement system [8,9] and the GOM Gmbh Atos structural light 3-D measurement system [10]. In China, Tsinghua University, Beijing University of Aeronautics and Astronautics, Shanghai Jiao Tong University, Xi’an Jiaotong University, Sichuan University, Huazhong University of Science and Technology, Tianjin University, and many other universities have also conducted systematic research on the stripe projection technology. At the same time, some commercial measurement systems, such as those produced by the Beijing Tianyuan 3D Technology Co., Ltd. (Beijing, China); Shanghai Shuzao Technology Co., Ltd. (Shanghai, China); Hangzhou Xianlin 3D Technology Co., Ltd. (Hangzhou, China); and Beijing Bowei Hengxin Technology Development Co., Ltd. (Beijing, China), etc., also launched 3-D measurement systems [11] based on stripe projection technology.
In industrial manufacturing, low measurement efficiency will delay the product’s production cycle, and the measurement environment needs to maintain a relatively stable state throughout the product’s production cycle; therefore, efficiency in the measurement process is needed to achieve automatic measurement. For the measurement of large, complex curved-surface components, the stripe projection sensors are required to measure the measured object in multiple poses, to realize the all-around coverage of the measured area of significant components. The high-precision fringe projection sensor needs to project multiple fringe images and keep the measured object stationary during the measurement process [12,13]. In addition, when measuring the global control points with the splicing camera based on close-up photogrammetry, it is also necessary to increase the number of poses to capture the global control points on the measured component to ensure accuracy [14]. The traditional close-view photogrammetry systems require manual operation, and the system is complex. Due to the advantages of high automation, high flexibility, and high efficiency of industrial robots, industrial robots were introduced into the measurement system. The fringe projection sensor and the stitching camera based on close-range photogrammetry are installed at the end of the robot through adapters to achieve automatic measurement within the range of motion [15]. Many researchers have optimized the scanning trajectory according to the parameters of the scanning system itself to obtain better scanning results and scanning efficiency. Generally, the scanning path can be generated by determining the point sequence of the scanning direction with the normal surface vector of the CAD model of the known part [16]. Nguyen [17] proposed a new path planning method that used a conformal map to stretch a 2D surface onto a 3D plane to control the overlap between two adjacent scanning paths. However, this method may cause a sudden change in the direction of the scanner between the two paths. Njaastad et al. [18] combined industrial robots with high-precision 3D cameras and laser distance sensors to form an automated scanning system to measure ship propeller blades. According to the shooting range of the 3D camera and the optimal proximity distance of the laser sensor, the scanning trajectory of the highly inclined ship propeller blade was planned with a resolution of 0.1 mm.
The above research shows that automatic 3D laser scanning is one of the key technologies to improve scanning efficiency and accuracy, and trajectory planning also plays an important role in automated laser-scanning systems. Aiming at the requirement of rapid thickness measurement in the chemical milling of casing parts, in this paper, a general scanning and measuring system composed of a six degrees of freedom industrial robot and tracking 3D scanner was designed. Its core part includes a six degrees of space parallel industrial robot and a three-dimensional tracking scanning system. The specific measurement process is as follows: (1) the three-dimensional point cloud data of the casing were obtained by automatically scanning the casing parts in the working area of the tracker through the six degrees of freedom industrial robot clamping the spherical scanning head; (2) the contoured fitting and post-processing of the point cloud data were carried out to realize the three-dimensional reconstruction of the measured workpiece; (3) the thickness measurement and analysis of the measured position were carried out in the Polyworks measurement software. For complex large-scale structural parts, it is time-consuming to plan the scanning path by manual teaching programming, and the spherical scanning head is easy to collide with the measured part. The accuracy of the scanning path also directly affects the scanning efficiency and the reconstruction rate of the measured parts. Therefore, the planning of the scanning path is particularly important for popularizing automated 3D scanning systems in the measurement of large structural parts [19].
In this paper, the constraints such as scanning area range, proximity distance and scene width angle of the tracking 3D scanner are analyzed. Combined with the geometric structure characteristics of the casing parts, the mathematical model of the annular scanning path was established by using the mathematical discretization modeling method. The change rules of scanning stroke, robot attitude transformation times and scanning time under circular scanning path were analyzed. The Robot Studio software platform was used to simulate the automatic scanning trajectory, which verifies that the scanning trajectory meets the constraints. The influence of the scanning path on the scanning reconstruction rate was analyzed by scanning experimental results.

2. Automated Tracking Type 3D Scanning Measurement System

2.1. System Composition

In the manufacturing process of aeroengine large parts, it is necessary to quickly measure the geometric dimensions of the part to guide the manufacturing process and evaluate the processing quality. This paper designed an automatic tracking 3D scanning-measurement system with a laser 3D scanner and industrial robot as the core, which mainly includes a 3D scanning data acquisition subsystem (3D scanning head, laser positioning tracking instrument), motion control subsystem (industrial robot, safety protection), turntable support subsystem and software system (automatic scanning and detection software, system automatic control software), as shown in Figure 2. The motion control subsystem includes a six degrees of freedom industrial robot (called a scanning robot), which can flexibly scan the 3-D scanning head in the trapezoidal working space of the laser positioning tracker to realize the automatic acquisition of the point cloud of the casing engine. The turntable support subsystem is a six degrees of freedom parallel-motion platform, which can fix the measured workpiece in the measurement space of the laser positioning tracker, to realize the scanning measurement of multiple poses of the large aviation components.
The I/O module software system was developed to connect the scanning robot to the laser scanner through the I/O signal communication. The communication structure is shown in Figure 3. The I/O control signal was written into the robot scanning trajectory program, and the signal received/sent by the I/O module controls the three-dimensional scanning probe, laser positioning tracker and safety protection accessories for automatic scanning operations to complete the acquisition of the point cloud data of the measured workpiece.
The operation process of the automatic tracking-type 3D scanning measurement system to collect the point cloud data and measurement of casing parts: (1) the laser scanner, optical tracker, scanning robot and six degrees of freedom motion platform were calibrated in advance; (2) the measured casing was placed on a six degrees of freedom parallel motion platform, and the platform pose was adjusted to expose the inner and outer surfaces of the casing to the binocular camera of the laser location tracker as much as possible; (3) the laser scanner parameters were configured according to the material and reflectivity of the measured workpiece, and the scanning template was established; (4) the scanning path *.TB program of the scanning robot was established by teaching pendant or off-line programming method; (5) start the automatic control system, invocation of the scanning template and the scanning path *.TB program to scan the measured workpiece area by area to obtain its point cloud data; (6) the point cloud data of the workpiece were post-processed and imported into the measurement software. The casing wall thickness was measured according to the pre-established measurement program, and the measurement report was generated.

2.2. System Scanning Principle and Technical Parameters

The automatic tracking 3D scanning measurement system takes the optical tracker and the three-dimensional scanner as the core components, and its scanning principle is shown in Figure 4. The coordinate system of the three-dimensional scanning probe is O B , the laser generator emits a laser line projected onto the scanned object, and the laser changes with the shape of the measured object. The coordinate system of the binocular camera of the 3D scanner is O B , and the camera can obtain the three-dimensional coordinates of the laser projected onto the scanned object. The point P is a point on the measured surface, and the distance to the surface of the object to be measured is L, which is represented by a blue laser line. A and B are the scan area sizes. During the scanning operation, the shape of the area composed of laser lines changes, as shown in the red laser line in Figure 4. The coordinate system of the optical tracker is O A . Through the reflective marker point-tracking technology, the coordinate position of the reflective target on the laser three-dimensional scanning probe was obtained at a speed of 60 Hz per second. The phase measurement method was used to calculate the three-dimensional coordinate value P i B j of each reflective target in the working space of the optical tracker, and the pose matrix P i B j of the three-dimensional scanner was obtained. If P i C is a point on a workpiece in the O B coordinate system, the three-dimensional coordinates of the laser line projection point in the binocular camera coordinate system O C are P i B , the three-dimensional coordinates in the optical tracker coordinate system O A are P i A :
P i A = T A B i T B C i P i C
where T B C i is the conversion matrix between the coordinate system O B and the binocular camera coordinate system O C of the 3D scanning probe; T A B i is the conversion matrix between the coordinate system O A of the optical tracker and the coordinate system O B of the 3D dimensional scanning probe. It is obtained from the attitude matrix P i B j of tracking target coordinates.
Using the linear three-dimensional information obtained by the three-dimensional scanner and the relative position of the optical tracker in space, when the scanner moves, the three-dimensional information of the position where the laser passes is continuously obtained, thereby forming continuous three-dimensional point cloud data. The free surface data or geometric data of the measured object will be obtained from each image and calculated by the software. The 3D laser scanner model in this paper is Metra SCAN 3D, and its specific parameters are shown in Table 1. The working distance between the scanner and the workpiece surface is 20–40 cm. When the optimal reference distance is 30 cm, the scanning field of the 3D scanner is 27.5 cm × 25 cm, as shown in Figure 5. Due to the use of active shutter technology, it has good adaptability to the surface of workpieces with different degrees of grayness, and can collect enough pictures to achieve efficient scanning.

3. Scan Path Planning

3.1. Casing Characteristics and Scanning Path Analysis

In the automatic tracking 3D scanning measurement system, the 3D scanning probe is used as the effector of the robot. The robot controls the scanner and moves in the measurement area of the optical tracker for scanning measurement. The scanner needs to move along a straight line and adjust its posture during the scanning process. Therefore, the critical factors in determining the scanning efficiency were the total length of the scanning path and the number of scanner attitude adjustments. In this paper, the scanning object is the aeroengine casing. The overall shape is a thin-walled conical barrel, and the surface has a multi-stiffener structure. The radial diameter random-casing shell changes in height, as shown in Figure 6. According to the geometric characteristics and wall-thickness measurement requirements of the part, to achieve high-efficiency and high-quality scanning reconstruction, a circular scanning path-planning method [20,21] was proposed.

3.2. Mathematical Modeling of the Circular Scan Path

To facilitate the 3D scanner to scan the inside and outside of the casing, the casing was fixed on the six degrees of freedom parallel tooling support platform in the direction of the radius. In order to avoid the interference collision between the 3D scanner and the six degrees of freedom support platform, bottom-up layer-by-layer circular scanning was adopted. The whole scanning process was divided into two parts: one part is a circular rotating scan around the inner and outer walls of the magazine, and the other part is a straight translation scan along the outer wall of the engine magazine. In the first part of the circular scanning process, the scanner performs a linear translation motion along the circumcircle tangent of the engine casing. When the scanner moved beyond the acquisition distance, the attitude of the scanner was adjusted. That is, the equidistant offset rotation scanning trajectory was fitted by a multi-segment linear translation motion. Therefore, the total length of the scanner’s linear movement and the number of pose adjustment times were important indicators to evaluate the efficiency of the scanning path-planning strategy.
The scanning starting point of the scanner was set as point A, and the movement mode of the robot within the effective measurement range of the scanner was linear translation. AW is the optimal measurement distance of 30 cm, θ is the angle of the longitudinal direction of the scanning sector, where θ = 50°, A i B i , A i C i is the distance from the laser line to the surface of the workpiece, where A i B i , A i C i 20   cm , 40   cm , H is the height of the casing, R 0 is the radius of the casing bottom and β is the cone of the casing. When the robot reaches the maximum capacity of the scanning sector near point A1, it completes a linear translation. When the scanner moves to A1, the scanning area changes, and the intersection points of the vertex, the boundary and the center O currently are, respectively, represented by red dotted lines in Figure 7a. Scanning along points A, A1, A2, …, An, and completing the first scan, the scanning trajectory is an external polygon, as shown in Figure 7a. The scanner moves points A, E, and F, upward along the outer wall of the casing, and completes the scanning of the A–F ring until the scanning was completed.

The Establishment of the Mathematical Model of Circular Scanning

  • Circular scanning process
It is assumed that the starting point of the scanner is A and OA is the distance from the scanning area to the center of the casing. The translational scanning motion reaches the pose adjustment point A1, OA1 is the distance from the scanning area to the center of the casing and then there is O A i = R i + A W . The deflection angle between OA and OA1 was γ, and the distance between point A and point A1 was L1. According to the geometric position relationship of the three, the length of the scanning trajectory in the measurement range of the scanner is as follows:
L = n i = 1 m tan γ ( R i + A W )
where i is the number of circles scanning; Ri is the scanning radius of circles scanning.
According to the height Hi of the casing and the cone angle β of the casing, the corresponding radius equation of the scanning is as follows:
R i = ( 2 i 1 ) 27.5 2 sin β + R 0
The number of times the scanner adjusts the pose n is as follows:
n = 360 ° γ
In Figure 7, after the scanner reaches the pose adjustment point A1, there are regional boundaries A1B1 and A1C1 with the geometric position generated on the outer surface of the casing. AW is the angular bisector of angular ∠BAC, and the following relationship was obtained:
cos ( α 2 γ ) = A i B i 2 + O A i 2 R i 2 2 A i B i O A i cos ( α 2 + γ ) = A i C i 2 + O A i 2 R i 2 2 A i C i O A i O A i = ( R i + A W ) cos γ
Equation (5) above can obtain the calculation formula of the deflection angle γ:
γ = arccos ( R + A 1 W 1 ) 2 + A C 2 R 2 2 ( R + A 1 W 1 ) A C α 2 γ = α 2 arccos ( R + A 1 W 1 ) 2 + A B 2 R 2 2 ( R + A 1 W 1 ) A B
where cos α 2 = A B 2 + O A 2 R 1 2 2 A B × O A .
2.
Judgment of the scanning boundary
The scanning field of view corresponding to the laser beam of the scanner is shown in Figure 8. When the initial position of the scanner is at point A, the scanning field width angles in both directions are individually α = 45° and θ = 50°. When the AW is 30 cm, the effective field width of the transverse scanning is 25 cm and the effective field width of the vertical scanning is 27.5 cm.
According to the standard geometric position relationship between the scanner and the cartridge in Figure 7, when the deflection angle γ = 0°, the boundary distance to the enclosure when the standard position was obtained according to Equation (6):
A B = A C = ( R + A W ) C O S α 2 ( R + A W ) 2 ( C O S α 2 ) 2 ( R + A W ) 2 + R 2
3.
The boundary judgment of the translation limit
When the scanner moves along the casing arc offset tangent, the laser beams AB and AC on both sides of the scanning area also move. The intersection point of the A1B1 laser beam irradiated on the casing arc is point B1. When A1B1 coincides with the intersection point W of the standard position angle branch to the arc, the length of A1B1min is the smallest, which is the boundary of the reference scanning area:
A 1 B 1 min = A W 2 + ( 25 2 ) 2 = 32.5   cm
According to Equations (7) and (8), when AW is in the effective area of the deflection angle, A1B1 conforms to the boundary constraint condition in the translation process. Then, judge A1C1, let A1C1max= 40 cm. According to Equation (6), the limit boundary judgment equation of scanning translation is as follows:
A 1 C 1 = ( R + A 1 W 1 ) C O S ( α 2 + γ ) ( R + A 1 W 1 ) 2 [ C O S ( α 2 + γ ) ] 2 ( R + A 1 W 1 ) 2 + R 2 A 1 B 1 = ( R + A 1 W 1 ) C O S ( α 2 γ ) ( R + A 1 W 1 ) 2 [ C O S ( α 2 γ ) ] 2 ( R + A 1 W 1 ) 2 + R 2
The variation range of the deflection angle is obtained by simplifying Equation (9):
0 γ arccos ( R + A 1 W 1 ) 2 + A C 2 R 2 2 ( R + A 1 W 1 ) A C α 2
where A C 40   cm ; A 1 W 1 = R + A W C O S γ R .
4.
The rising posture adjustment process
After completing the first ring scan, the scanner returns to the starting point A. The scanner continued to move up along the outer wall of the casing to reach point E, and the distance between the two points is K1. Move up m times in a turn, obtain the final point F, and complete m-cycle scanning:
m = H cos β × 27.5
K 1 = H i H i 1 cos β ; i 2
where 27.5 cm is the effective scanning area of the scanner in the vertical direction.
Assuming that the total distance of the robot moving upward is K, according to Equations (11) and (12) above, the mathematical expression of the upward translation parallel to the outer wall of the casing is as follows:
K = m k 1 = H ( H i H i 1 ) 27.5 ( cos β ) 2
According to Equations (2) and (13), the mathematical expression of the total distance Q of the circular scanning is as follows:
Q = n i = 1 m tan γ ( R i + A W ) + H ( H i H i 1 ) 27.5 ( cos β ) 2
The total number of posture adjustments of the annular scanning method is n:
n = i = 1 m 360 ° γ i

3.3. Repeat the Scan Area Compensation Method

Substitute the measured casing height H = 60 cm into Equation (11) to obtain the number of scans, not the number of full turns. In the scanning process, the number of scans should also be rounded to achieve complete scanning reconstruction. In this paper, each circle’s scanning height was adjusted to decrease the scanning area, and the repeated scanning area was divided into each ring. The corresponding scanning radius Ri will change, and the average coefficient v compensation method was used to compensate for the repeated scanning area. Taking the scanning casing as an example, there are:
v = 3 2.198 3 = 0.267
H i v = H H × v
According to Equations (3), (16) and (17), the distance between the scanner and the casing before and after scanning was obtained by adjusting the scanning height, as shown in Table 2.
According to the data in Table 2, the corresponding position relationship diagram before and after scanning height adjustment is shown in Figure 9. It can be seen from Figure 9 that the adjusted Hiv is lower than the corresponding height H of the original scan, and the corresponding Riv of the scan is lower than the original scan radius R as a whole.
Using the above established circular scanning path-planning method, the relevant data of the scanning path before and after the scanning height adjustment were calculated as shown in Table 3 and Table 4.
If the scanning speed is v = 0.04 m/s and the posture adjustment time is t = 0.5 s/n, the total time to complete the casing scanning before and after the adjustment is as follows:
T before = Q v + N t = 381.4175 + 88 = 469.417   s T after = Q v + N t = 373.375 + 91 = 464.375   s
Through the analysis of the relevant data of the scanning path before and after the above scanning height adjustment it can be seen that the scanning range distribution can be changed by adjusting the scanning height, the complete scanning of the workpiece can be realized and the scanning quality is better than before the adjustment.

3.4. Analysis of the Scan Path Parameters

The casing geometry is shown in Figure 7, where the casing cone angle β is 7.125°. According to the corresponding data after adjusting the scanning height in Table 4, the ring scanning radius R1 of the scanner in the first circle is 43.75 cm. In the process of annular scanning in the limited scanning area, the changes in related parameters such as γ, A1W1, A1B1, A1C1 and L in the scanning path are shown in Table 5. The relationship between the deflection angle and the scanning distance was plotted, as shown in Figure 10.
It can be seen from Table 5 that the deflection angle γ at 0–5° is an effective scanning range in the first circular scanning. When γ = 0°, the scanner is in the standard position, and the distance between A1C1 and A1B1 and the scanner is equal. The laser beam at the boundary of the scanner is an isosceles triangle. It can be seen from Figure 10 that A1W1 and A1C1 increase with the increase in γ. When the deflection angle γ is between 0–4°, A1B1 decreases with the increase in γ. When the deflection angle γ is 4–5°, A1B1 increases with the increase in γ.
During the second and third rounds of scanning, the radius of the casing gradually increases, and the corresponding scanning radius and deflection angle increase to varying degrees. The data of the relationship between the deflection angle of the second and third turns and the scanning distance are shown in Table 6 and Table 7, respectively.
Figure 11 shows the relationship between the deflection angle, scanning radius and scanning distance in the second circular scanning. It can be seen from Figure 11 that when the scanner rises for scanning, the deflection angle γ increases with the increase in the corresponding scanning radius and the corresponding scanning path also increases.

4. Scanning Simulation Analysis

4.1. Description of the Scan Trajectory Motion

In the circular scanning motion, the scanning motion is regarded as the motion of the robot tool coordinate system {T} relative to the workpiece coordinate system {S}, both translational motion and rotational motion. As shown in Figure 12, the workpiece coordinate system {S} is used as the reference coordinate system, and the starting point is set to P0. Due to the limited scanning range of the scanner, it is necessary to adjust the attitude when reaching the coordinate point P1. The scanning trajectory of the first circle passes through the coordinate point Pi, Pi+1, …, Pn in turn. The rising part of the trajectory passes through the coordinate Pj, Pj+1, the scanner will only have a translational motion until the scan was complete. The spatial pose of the scanner’s scanning trajectory coordinates is expressed as:
P 0 S = P 0 T + P S 0 T ; P 1 S = T 1 0 P 0 S ; P 2 S = T 1 0 T 2 1 P 0 S ; P 3 S = T 1 0 T 2 1 T 3 2 P 0 S ; P i S = T 1 0 T 2 1 T 3 2 T i i 1 P 0 S ; P i + 1 S = T 1 0 T 2 1 T 3 2 T i i 1 T i + 1 i P 0 S ; p n s = T 1 0 T 2 1 T 3 2 T i i 1 T i + 1 i T n n 1 P 0 S = T n 0 P 0 s P m s = T m 0 P 0
The coordinate expression of any coordinate point in the scanning trajectory is as follows:
P n m s = ( T s n 0 + T m 0 ) P 0
where T n 0 = r o t ( z , i = 0 n γ i ) t r a n s ( i = 0 n l i ) = cos i = 0 n γ i sin i = 0 n γ i 0 0 sin i = 0 n γ i cos i = 0 n γ i 0 i = 0 n l i 0 0 1 0 0 0 0 1 ; T m 0 = t r a n s Z ( j = 2 m k j ) = 0 0 K j .
The curve circle model is created with (R + AW) as the radius. The coordinate P 0 = [ R i + AW , l i , H i ] T of the starting point of the scanning path is set. Based on Equation (20), the coordinate position of the coordinate point Pi on the scanning trajectory was solved, and the scanning trajectory was fitted.

4.2. Establishment of Scanning Simulation Model

To verify the rationality of the above scanning path planning, the robot studio robot simulation software is used to simulate the virtual scanning trajectory motion under the condition of 3D scanner parameter constraints. The ABB-IRB6700 model robot and the casing part model were imported into the software, and the robot system was created by the ‘layout’ method. Select the frame creation method to establish the tool and workpiece coordinate system, generate the scanner tool and install it on the robot end effector. Adjust the position and angle of the imported model to make it in the robot workspace.
In the workpiece coordinate system, the scanning motion parameters were set, the running speed is v = 0.04 m/s, the acceleration is a = 0.5 m/s and the posture adjustment time is t = 0.5 s. The obtained scanning trajectory coordinate points were used to set the robot’s motion trajectory, generate the robot’s motion path and complete the robot’s motion trajectory offline programming. Run the scanning trajectory program to perform automatic scanning trajectory simulation motion, and the scanning motion process is shown in Figure 13.
In the process of simulation scanning motion, the robot and the 3D scanner did not interfere or collide with the casing. The scanning range conforms to the scanning space constraint condition of the scanner, and the running time is close to the theoretical analysis running time. The off-line scanning trajectory programming based on the geometric characteristics of the measured parts was realized, and the feasibility and correctness of the scanning path planning were also proved.

5. Scanning Measurement Process Experiment

The scanning trajectory of this experiment adopts the scanning trajectory program generated by off-line programming of simulation software, which was imported into the automatic measurement system by downloading and the scanning speed is 0.04 m/s. According to the geometric characteristics and measurement requirements of the casing, the scanning parameter template and the thickness measurement template were established, respectively, to realize the automatic scanning reconstruction and automatic measurement of the casing. The scanning process is shown in Figure 14.
The experimental results show that the whole ring scanning time is about 15 min, including robot start-up and return movement. The casing scanning reconstruction model has a complete quality and no obvious noise. For the measurement of the casing parts, the total time of scanning and measurement using this system was 60 min, including scanning path programming, measurement template preparation, scanning and post-processing. However, manual measurement using an ultrasonic thickness gauge requires two hours for six people to measure at the same time and the scanning efficiency is increased by more than 80%. The least squares algorithm was used to fit the obtained point cloud data and the model was reconstructed after the workpiece was scanned, as shown in Figure 15.
By establishing the thickness measurement template, the thickness measurement was completed for the engine magazine works. The thickness deviation between the scanned model and the design model is shown more intuitively in the form of a chromatogram to evaluate the quality of the casing milling, as shown in Figure 16 and Figure 17. Comparing the measured values of the three-dimensional scanning in the same area with the measured values measured by the ultrasonic thickness gauge, the deviation of the measured values is within 0.1–0.2%. This error is likely to be caused by the deviation of the measurement position. The measurement accuracy of the measurement method can meet the measurement of such workpieces.

6. Conclusions

In this paper, a non-contact automatic three-dimensional scanning measurement process system and process are proposed and established for the rapid and high-precision measurement requirements of aero-engine casing parts. The scanning trajectory is analyzed and verified by mathematical modeling, simulation analysis and process test. The main conclusions are as follows:
(1) According to the relationship between the scanning area constraint of the scanner and the geometric size of the casing, a mathematical model of the annular scanning trajectory is established. The scanning path planning from part model to off-line robot programming was realized. The scanning path-planning method is suitable for the scanning measurement of cylindrical parts with tapering, such as casing parts;
(2) It is proposed to evenly distribute the scanning repetition area by adjusting the scanning height, and introducing the square compensation coefficient to compensate for the repeated scanning area. Through the numerical calculation of the scanning path before and after the adjustment, it is found that the allocation method shortens the scanning path by 1.1%, and the scanning integrity is also improved;
(3) The robot-studio software is used to conduct the motion simulation of the optimized scanning path and verify the scanning process experiment. The results show that the annular trajectory scanning method established in this paper can meet the scanning measurement requirements of casing parts. The efficiency of scanning measurement is more than 80% higher than that of artificial ultrasonic point-by-point measurement.

Author Contributions

Conceptualization, J.L. and M.W.; methodology, L.Q.; formal analysis, G.L.; investigation, J.L.; data curation, J.L.; writing—original draft preparation, J.L.; writing—review and editing, M.W.; supervision, L.Q.; project administration, G.L.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Liaoning Province ‘Xingliao Talent Plan’ project, grant number XLYC2002086; the Discipline Construction Fund, grant number 310121001 and the National Key Research and Development Plan by Liaoning, grant number 2022JH2/101300211.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Srinivas, G.; Kurkal, R.; Shenoy, S.B. Recent developments in turbomachinery component materials and manufacturing challenges for aero engine applications. IOP Conf. Ser. Mater. Sci. Eng. 2018, 314, 012012. [Google Scholar] [CrossRef]
  2. Xu, F.; Yang, M.; Chai, X.; Yan, K.; Ni, X. Lightweight design method of aero-engine combustor case. Aeronaut. Sci. Technol. 2021, 32, 27–34. [Google Scholar] [CrossRef]
  3. Liu, L.; Wang, E.; Zhang, X.; Liang, W.; Li, X.; Xie, H. MEMS-based 3D confocal scanning microendoscope using MEMS scanners for both lateral and axial scan. Sens. Actuators A Phys. 2014, 215, 89–95. [Google Scholar] [CrossRef] [PubMed]
  4. Popov, I.; Onuh, S.; Dotchev, K. Dimensional error analysis in point cloud-based inspection using a non-contact method for data acquisition. Meas. Sci. Technol. 2010, 21, 075303. [Google Scholar] [CrossRef]
  5. Gerde, J.R.; Christens-Barry, W.A. Evaluating a hybrid 3-dimensional metrology system: Merging data from optical and touch probe devices. In Proceedings of the SPIE Optical Engineering+ Applications 2011, San Diego, CA, USA, 21–25 August 2011. [Google Scholar] [CrossRef]
  6. Chu, C.L.; Wu, C.Y. Development of 3D touch trigger probe with real-time observation. In Proceedings of the 6th International Symposium on Precision Engineering Measurements and Instrumentation, Hangzhou, China, 8–10 August 2010. [Google Scholar] [CrossRef]
  7. Elkott, D.; Veldhuis, S. CAD-based sampling for CMM inspection of models with sculptured features. Eng. Comput. 2007, 23, 187–206. [Google Scholar] [CrossRef]
  8. Hot, J.; Dasquec, A.; Topalov, J.; Mazars, V.; Ringot, E. Titanium valorization: From chemical milling baths to air depollution applications. J. Clean. Prod. 2020, 249, 119–344. [Google Scholar] [CrossRef]
  9. Phan, N.D.M.; Quinsat, Y.; Lavernhe, S.; Lartigue, C. Scanner path planning with the control of overlap for part inspection with an industrial robot. Int. J. Adv. Manuf. Technol. 2018, 98, 629–643. [Google Scholar] [CrossRef]
  10. Mahamood, R.M.; Akinlabi, E.T. Chemical Cutting Process. In Advanced Noncontact Cutting and Joining Technologies; Springer: Berlin/Heidelberg, Germany, 2018; pp. 11–25. [Google Scholar] [CrossRef]
  11. Zhang, C.; Zhao, H.; Qiao, J.; Zhou, C.; Zhang, L.; Hu, G.; Geng, H. Three-dimensional measurement based on optimized circular fringe projection technique. Opt. Express 2019, 27, 2465–2477. [Google Scholar] [CrossRef] [PubMed]
  12. ElKott, D.F.; Veldhuis, S.C. Isoperimetric line sampling for the inspection planning of sculptured surfaces. Comput.-Aided Des. 2004, 37, 189–200. [Google Scholar] [CrossRef]
  13. Duan, X.; Wang, C.; Wang, J.; Zhao, H. A new calibration method and optimization of structure parameters under the non-ideal condition for 3D measurement system based on fiber-optic interference fringe projection. Opt. Int. J. Light Electron. 2018, 172, 424–430. [Google Scholar] [CrossRef]
  14. Wang, Y.D. The Applications of Industrial Robots in the Field of Industrial Automation and Control. Int. Conf. Front. Manuf. Sci. Meas. Technol. 2013, 401–403, 1716–1719. [Google Scholar] [CrossRef]
  15. Li, D.; Tian, J.; Yang, X. A new 3D measurement method and its calibration based on the combination of binocular and monocular vision. Proc. SPIE-Int. Soc. Opt. Eng. 2012, 8563, 9. [Google Scholar] [CrossRef]
  16. Boulch, A.; Guerry, J.; Le Saux, B.; Audebert, N. 3D point cloud semantic labeling with 2D deep segmentation networks. Comput. Graph. 2018, 71, 189–198. [Google Scholar] [CrossRef]
  17. Yang, Y.; Li, Z.; Yu, X.; Li, Z.; Gao, H. A trajectory planning method for robot scanning system using mask R-CNN for scanning objects with unknown model. Neurocomputing 2020, 404, 329–339. [Google Scholar] [CrossRef]
  18. Njaastad, E.B.; Munthe-Kaas, N.H.; Egeland, O. Robotic Autoscanning of Highly Skewed Ship Propeller Blades. IFAC-Pap. OnLine 2018, 51, 435–440. [Google Scholar] [CrossRef]
  19. Song, C.K.; Kim, S.W. Reverse engineering: Autonomous digitization of free-formed surfaces on a CNC coordinate measuring machine. Int. J. Mach. Tools Manuf. 1997, 37, 1041–1051. [Google Scholar] [CrossRef]
  20. Fan, K.C.; Cheng, F.; Wang, W.; Chen, Y.; Lin, J.Y. A scanning contact probe for a micro-coordinate measuring machine (CMM). Meas. Sci. Technol. 2010, 21, 054002. [Google Scholar] [CrossRef]
  21. Saadat, M.; Cretin, L. Measurement systems for large aerospace components. Sens. Rev. 2002, 22, 199–206. [Google Scholar] [CrossRef]
Figure 1. Research background for the key components of aero-engines.
Figure 1. Research background for the key components of aero-engines.
Applsci 13 06162 g001
Figure 2. Automatic tracking type 3D scanning measurement system.
Figure 2. Automatic tracking type 3D scanning measurement system.
Applsci 13 06162 g002
Figure 3. System communication structure.
Figure 3. System communication structure.
Applsci 13 06162 g003
Figure 4. Scanner point cloud scanning principle.
Figure 4. Scanner point cloud scanning principle.
Applsci 13 06162 g004
Figure 5. Metra SCAN 3D scanning acquisition area.
Figure 5. Metra SCAN 3D scanning acquisition area.
Applsci 13 06162 g005
Figure 6. Aero engine casing part.
Figure 6. Aero engine casing part.
Applsci 13 06162 g006
Figure 7. Circular scanning trajectory diagram. (a) Circular scanning top view; (b) Circular scanning main view.
Figure 7. Circular scanning trajectory diagram. (a) Circular scanning top view; (b) Circular scanning main view.
Applsci 13 06162 g007
Figure 8. The horizontal and vertical field of view of the scanning area of the scanner.
Figure 8. The horizontal and vertical field of view of the scanning area of the scanner.
Applsci 13 06162 g008
Figure 9. Scanning height adjustment before and after scanning the corresponding position diagram.
Figure 9. Scanning height adjustment before and after scanning the corresponding position diagram.
Applsci 13 06162 g009
Figure 10. Curve diagram of the relationship between deflection angle and scanning distance.
Figure 10. Curve diagram of the relationship between deflection angle and scanning distance.
Applsci 13 06162 g010
Figure 11. Deflection angle and scan corresponding parameters.
Figure 11. Deflection angle and scan corresponding parameters.
Applsci 13 06162 g011
Figure 12. Schematic diagram of the circular scan trajectory.
Figure 12. Schematic diagram of the circular scan trajectory.
Applsci 13 06162 g012
Figure 13. Scanning simulation motion. In the figure, the blue line, the red line and the green line are the x, y, z axes of the trajectory point coordinate system, and the yellow arrow is the trajectory motion direction.
Figure 13. Scanning simulation motion. In the figure, the blue line, the red line and the green line are the x, y, z axes of the trajectory point coordinate system, and the yellow arrow is the trajectory motion direction.
Applsci 13 06162 g013
Figure 14. Casing scanning process.
Figure 14. Casing scanning process.
Applsci 13 06162 g014
Figure 15. After the scanned artifact reconstruction model.
Figure 15. After the scanned artifact reconstruction model.
Applsci 13 06162 g015
Figure 16. Casing surface thickness chromatogram.
Figure 16. Casing surface thickness chromatogram.
Applsci 13 06162 g016
Figure 17. Casing surface thickness deviation chromatogram.
Figure 17. Casing surface thickness deviation chromatogram.
Applsci 13 06162 g017
Table 1. Metra SCAN 3D scanner parameters.
Table 1. Metra SCAN 3D scanner parameters.
Accuracy (mm)Resolution Ratio (mm)Base Distance L (mm)Scan Area A × B (mm)Measuring Rate (r/s)
0.030.05300275 × 250480,000
Table 2. The distance relationship between the scanner and the corresponding position of the casing before and after adjusting the scanning height (unit cm).
Table 2. The distance relationship between the scanner and the corresponding position of the casing before and after adjusting the scanning height (unit cm).
miHiRiHivRiv
i = 1 13.6444.219.9943.75
i = 2 40.9247.6229.9946.25
i = 3 68.2051.0349.9948.75
Table 3. The relevant data of the casing scanning path before adjustment (in cm).
Table 3. The relevant data of the casing scanning path before adjustment (in cm).
Order NumberAngle of DeflectionThe Number of PosturesSurround Path LengthRising Path LengthTotal Path Length
15.8°63472.50472.5
26.2°59487.7427.49487.74
36.7°54510.44727.49510.447
Total1761470.68754.981525.67
Table 4. The relevant data of the casing scanning path after adjustment (in cm).
Table 4. The relevant data of the casing scanning path after adjustment (in cm).
Order NumberAngle of DeflectionThe Number of PosturesSurround Path LengthRising Path LengthTotal Path Length
15.5°664620462
260487.220.15507.35
36.5°5650420.155245.1
Total 1821453.240.31493.5
Table 5. Relationship Data between Deflection Angle and Scan Distance (first circle) (in cm).
Table 5. Relationship Data between Deflection Angle and Scan Distance (first circle) (in cm).
NumberγA1B1A1C1A1W1L
134.7034.70300
234.2335.2530.011.29
333.8335.8930.062.58
433.5336.6730.163.87
533.3637.6430.345.18
633.3238.8530.626.50
733.4540.3931.037.86
833.7642.3731.599.25
Table 6. Deflection angle and scan distance relationship data (second circle) (in cm).
Table 6. Deflection angle and scan distance relationship data (second circle) (in cm).
NumberγA1B1A1C1A1W1L
134.5634.56300
234.1035.0830.011.33
333.7235.7030.062.66
433.4536.4530.164.01
533.2937.3830.355.35
633.2738.5530.646.72
733.4240.0431.068.12
833.7641.9431.649.56
Table 7. Relationship Data of Deflection Angle and Scan Distance (third circle) (in cm).
Table 7. Relationship Data of Deflection Angle and Scan Distance (third circle) (in cm).
NumberγA1B1A1C1A1W1L
134.4334.43300
233.9934.9430.011.37
333.6335.5330.052.75
433.3736.2630.164.13
533.2337.1730.365.531
633.2338.3030.666.94
733.4139.7531.098.38
833.7741.5931.699.87
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, J.; Wang, M.; Qu, L.; Lv, G. Research on Three-Dimensional Scanning Path Planning of Casing Parts Based on Industrial Robot. Appl. Sci. 2023, 13, 6162. https://doi.org/10.3390/app13106162

AMA Style

Li J, Wang M, Qu L, Lv G. Research on Three-Dimensional Scanning Path Planning of Casing Parts Based on Industrial Robot. Applied Sciences. 2023; 13(10):6162. https://doi.org/10.3390/app13106162

Chicago/Turabian Style

Li, Jing, Minghai Wang, Ligang Qu, and Guangming Lv. 2023. "Research on Three-Dimensional Scanning Path Planning of Casing Parts Based on Industrial Robot" Applied Sciences 13, no. 10: 6162. https://doi.org/10.3390/app13106162

APA Style

Li, J., Wang, M., Qu, L., & Lv, G. (2023). Research on Three-Dimensional Scanning Path Planning of Casing Parts Based on Industrial Robot. Applied Sciences, 13(10), 6162. https://doi.org/10.3390/app13106162

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop