A zone-layered trimming method for ceramic core of aero-engine blade based on an advanced reconfigurable laser processing system

Ceramic structural parts are one of the most widely utilized structural parts in the industry. However, they usually contain defects following the pressing process, such as burrs. Therefore, additional trimming is usually required, despite the deformation challenges and difficulty in positioning. This paper proposes an ultrafast laser processing system for trimming complex ceramic structural parts. Opto-electromechanical cooperative control software is developed to control the laser processing system. The trimming problem of the ceramic cores used in aero engines is studied. The regional registration method is introduced based on the iterative closest point algorithm to register the path extracted from the computer-aided design model with the deformed ceramic core. A zonal and layering processing method for three-dimensional contours on complex surfaces is proposed to generate the working data of high-speed scanning galvanometer and the computer numerical control machine tool, respectively. The results show that the laser system and the method proposed in this paper are suitable for trimming complex non-datum parts such as ceramic cores. Compared with the results of manual trimming, the method proposed in this paper has higher accuracy, efficiency, and yield. The method mentioned above has been used in practical application with satisfactory results.


Introduction
Ceramic cores, which are employed to form cavities in the turbine blade, are one of the main components for casting hollow high-performance aircraft engine turbine blades. The surface quality of the ceramic core directly affects the forming process as well as the casting quality of the monocrystalline silicon blade [1]. Ceramic cores are manufactured via pressing and sintering processes. Thus, they usually contain manufacture-related defects such as burrs and blocked holes, as presented in Fig. 1. These defects have to be trimmed to meet the requirements of blade casting.
However, the ceramic material is brittle, the surface is rough, and the contact pressure between the tool and the workpiece is relatively high. Thus, it is often difficult to position and process the ceramic core on the traditional machining systems. Currently, ceramic core trimming relies on manual hand-held diamond abrasive tools or small electric drills, resulting in low efficiency, poor product consistency, and low yield. Laser processing, one of the non-contact processing methods, can be employed to process hard and brittle materials such as glass or single-crystal silicon [2][3][4]. In recent years, laser processing technology has been rapidly expanding to many industrial application fields, such as laser drilling [5][6][7], laser-assisted milling [8][9][10], laser polishing [11][12][13][14], and laser surface texturing [15][16][17][18]. Although lasers can be used to trim ceramic cores, there are still some key issues that need to be resolved, which are summarized as follows: 1) The surface is rough, and the core is deformed compared to the CAD model, making it impossible to effectively locate by using a fixture and trimming the burrs according to the path generated by the CAD.
2) The trimming profile is a set of complex 3D spatial curves, so an advanced laser processing system with fiveaxis processing functions and unique processing methods is needed to achieve trimming.
The trimming problem of the ceramic core can be summarized as the laser trimming problem for complex non-datum components. Researches on this problem are rare at present, and similar research mainly focuses on the repair of damaged areas of aero-engine blades [19][20][21][22][23]. Yilmaz et al. [19] investigated the repair methodology of thin-curved blades. In this methodology, non-contact optical measurement equipment is used to obtain 3D data on the surface of the damaged blades. Furthermore, a three-axis machining strategy is employed to remove the excess weld material. Both smoothness and error requirements of the repair path are investigated in depth. Bremer [21] described the importance of automatic blade repair and provided a complete set of blade automatic repair processes, including inspection, weld preparation, welding and reprofiling. In another study [20], the authors detailed reverse engineering technology application in blade repair. Liu et al. [22] investigated in detail the application of additive manufacturing technology for repairing military aircraft components. Gao et al. [23] described 3D non-contact measurement technology in blade repair and equipment integration. Chen et al. [24] proposed a repair method based on non-damaged blades in service at the same time instead of the original design model so that the reconstructed target repairing surface will more accord with the deformation state of damaged blades. However, the research above mainly focuses on the repair of metal components. Compared with ceramic cores, the repair of metal components such as blades focuses on the reconstruction of the damaged part of the blade, and rarely considers the deformation of blades, and there is no positioning problem. At the same time, most of the utilized processing systems are traditional cutting machine tools.
For the problem that the trimming path extracted by the computer-aided design (CAD) model cannot be used due to the core's deformation, the laboratory researchers have proposed a corresponding solution. Based on the iterative closest point (ICP) registration algorithm and the subregional operation of the CAD model and the scanned point cloud model, Hou et al. [25] realized the registration of CAD extracted path and deformed ceramic core. There are many similar kinds of research [26][27][28][29][30][31] on estimating the spatial pose of objects or matching edges based on the point cloud registration algorithm. Fan et al. [26] pointed out that edges are more common than feature points or feature planes, so edges can be used to register point clouds with CAD models. Vock et al. [27] proposed a sampling strategy with hypothesis validation to drastically reduce runtime. They demonstrated that the point cloud-based pose estimation was helpful for fast matching. Xie et al. [30] proposed an iterative varianceminimization matching method to address measured points that have measuring defects. Furthermore, a distance based on point to tangent distance and first-order point-to-point distance is presented to construct the objective function to speed up convergence. The results show that the method has excellent stability in dealing with measured points with measuring defects. Li et al. [31] proposed a new partial shape matching method via constructing an objective function considering the requirements associated with different grinding allowances. In addition to the matching method, Xiao et al. [32] proposed a tool path generation method directly with point clouds without surface fitting. Machining areas could be recognized from the entire blade model by splitting the point cloud and analyzing geometric parameters of points. The cutter location point is generated by extending on the normal vector direction of the corresponding point.
In addition to pose measurement and path generation problems, advanced laser processing systems are also essential for the laser trimming of complex components. Nowadays, there are many studies on the function expansion of the laser processing system. Such investigations aim to achieve higher equipment functions utilized in special material processing, composite processing, and complex component processing. Many researchers [33][34][35] have conducted investigations in this  field. Erkorkmaz et al. [33] introduced an advanced fiveaxis laser drilling system to achieve flying drilling of groups of holes on curved surfaces. Jiang et al. [34] introduced a laser lithography system for surface texturing. In addition, the authors achieved texture marking on large-scale surfaces and proposed a new laser layered lithography method for large-scale surface texturing. Cuccolini et al. [35] introduced a five-axis laser milling system composed of a two-dimensional (2D) high-speed galvanometer and a five-axis computer numerical control (CNC) machine tool. This milling system can achieve pattern-like laser milling on complex surfaces. Aiming at the trimming problem of the ceramic core, an advanced ultrafast laser processing system is developed in this paper. A laser processing head with a high-speed scanning galvanometer inside is installed on a conventional five-axis CNC machine tool. In order to scan the clamped deformed ceramic core to obtain the point cloud model and achieve registration of CAD extracted path and deformed ceramic core, a stereo vision measuring device is integrated within the system. An opto-electromechanical cooperative control (OEMCC) software is developed to realize the collaborative work of the CNC machine tool, laser, galvanometer and stereo vision measuring device. A zonal and layering processing algorithm for contour clusters on freeform surfaces is proposed to deal with the complex spatial path. The algorithm generates the processing data of both the highspeed scanning galvanometer and the CNC machine tool.
The content of this paper is arranged as follows. In Section 2, the zonal and layering processing algorithm is introduced in detail. In Section 3, the hardware and software of the laser system are introduced, and the trimming process of the ceramic core is given. In Section 4, the experiments and the results of the ceramic core trimming are presented. Finally, the conclusions are made in Section 5.
2 Laser zone-layered trimming method for surface contour cluster 2.1 Principle of the laser layering process Due to the 3D processing capability of the laser focal depth, the 3D contours can be layer-processed by a 2D laser galvanometer [34]. According to the energy distribution characteristics in the proximity of the laser beam focus, the principle of 3D contour layering processing is provided. Following the focus of the laser beam via field lens, the energy distribution near the focus can be expressed as Eq. (1) [36]: where r is the laser spot radius, is the energy density at the center of the spot, which can be expressed by the laser power and the spot radius at the focal, is the radius of the laser beam at the waist and is the radius at various positions along the optical axis, and z represents the distance to the center of the focus along the optical axis. Parameters and can be expressed as Eqs.
(2)-(6) [36]: where parameter P represents the laser power, h is the repetition frequency, E is the single pulse energy, is the laser wavelength, is the focal length of the field lens, M 2 is the mode parameter that characterizes the beam quality and the value is 1.1, D is the diameter of the laser beam, and is the Rayleigh length. In this research, a picosecond laser with a power of 53.59 W, a repetition frequency of 250 kHz, and a wavelength of 532 nm is used for ceramic core trimming, an f-theta lens with an effective focal length of 259.4 mm is used for laser focus. The spot radius at the focal plane can be obtained from Eq. (6) after the experimentally measured beam diameter . The measurement results of the beam quality at the entrance of the field lens are shown in Fig. 2. The ellipticity of the spot is 89.5%. The minor ellipse diameter is 4.468 mm, and the major ellipse diameter is 4.994 mm. The beam diameter D takes the average of the minor and the major ellipse diameters, the value is 4.7313 mm. Then the relevant laser parameters can be calculated by Eqs. (1)- (6). The single pulse energy E = 214.36 µJ, the Rayleigh length mm, the spot radius 20.4 µm, and the energy density .
The focus energy distribution is depicted in Fig. 3, where the energy density at the focal (z = 0) is 32.7248 J/cm 2 . The distribution of laser energy on the focal plane at various positions near the focus is shown in Fig. 3(a). The variation of the maximum energy density near the focus is presented in Fig. 3(b). It can be observed that the attenuation of the maximum energy density increases with the increase in distance from the focus.
According to Fig. 3, the variety of laser energy density can be ignored within a specific range near the focus. The laser energy density within this range can meet the processing requirements, so the range is defined as laser focus processing length (LPL) [34]. The size of LPL is I 0 z related to the energy density that meets the processing requirements. For example, if an energy density higher than 97% of the maximum energy density can meet the processing requirements, the value can reach 0.4 mm. Therefore, even if a scanning galvanometer is often employed as a 2D galvanometer in practical applications, the laser has 3D processing ability, which is the basis of layered processing.

Theoretical basis for laser zone trimming of ceramic core contour clusters
It is challenging to ensure that the incident direction of the laser is consistent with the ideal direction during the actual processing when processing the contours of the curved surface. Generally, there is an angle between the direction of the optical axis and the ideal processing direction of the contour. As shown in Fig. 4, this angle is defined as the incident angle of the laser, and the tilt of the optical axis will affect the shape and size of the spot near the laser focus, which in turn affects the line width of laser engraving. Therefore, to ensure that the size of the laser engraving is within the allowable range of the design requirements, the range of the laser incident angle needs to be limited.
Generally, the energy distribution of the laser beam can be described by Gaussian distribution. In the case of unbiased incidence, the laser spot radius r can be expressed as Eq. (7) [37]: I th θ where the parameter represents the ablation threshold of the material. When the incident angle of the laser is , the changes of the laser spot size is as shown in Fig. 5. The laser spot on the curved surface is approximately an ellipse, and the spot size can be further expressed as Eq. (8) [38]: r min r max where and respectively represent the radius size along the ellipse's minor axis and major axis. When the material and laser parameters are fixed, the change rate of the laser spot size with the angle of incidence is shown in Fig. 5.
It can be observed in Fig. 5 that when the laser incident angle is less than a certain angle, the change rate of the laser spot size is close to 1, that is, the laser spot size is approximately unchanged. For example, when the value  I 0 / I th is 2, and the incident angle is 5°, the change rates of the minor and major semi-axes of the spot are 99.72% and 100.11%, respectively. With the laser incidence angle increasing, the change rate will change drastically, so in actual processing, the incident angle of the laser should be limited to a specific range. In order to explain the influence of the incidence angle on the line width, the laser grooving experiment under different incidence angles is designed. The experimental device is shown in Fig. 6. The laser processing head is configured on a five-axis motion platform, and the silicon wafer is clamped by the plane fixture. The rotation of the five-axis motion platform is used to realize the laser grooving with different incident angles.
It is necessary to find the position of the laser focus before grooving. For this purpose, an experiment of finding the focus is carried out. First, the five-axis motion platform moves along the Y-axis and scans a groove on the silicon wafer with a galvanometer at a certain distance. The change of the groove width near the focus is shown in Fig. 7. When the coordinate of the Y-axis is -155. 16, the width of the groove is the smallest, which is 47 µm. Take this position as the focus position under the current clamping condition.
Engrave the first groove with the focus position, then rotate the turntable, and use the laser with different incident angles to engrave the groove. The groove width is shown in Fig. 8. It can be seen from the measurement data that when the incident angle is less than 5°, the groove width of the main ablation zone changes little, which is consistent with the results given in Fig. 5. Therefore, a specific incident angle can be allowed in actual processing, which is also the basis of the zonal laser processing of the surface contours.   Typically, the position coordinates and the normal vectors of the contours can be extracted by computer-aided manufacturing (CAM) software and the corresponding CAD models. This information enables the preprocessing of the ceramic core trimming path. The geometric information of the surface contours is shown in Fig. 9, including the contour position data and the contour feature vectors. These contours can be divided into three types: The first type is a cluster of circular contours on the curved surface, while the second and the third type are two parts of the same strip contour. The reason for dividing a complete strip contour into two smaller parts with corresponding feature vectors is that the direction of the feature vectors should be more conducive to the layering operation.
In Fig. 10, the extracted contour clusters of the ceramic core designated for trimming and the corresponding feature vectors are presented. There are a total of 60 feature vectors corresponding to 60 contours. The corresponding data is provided in detail in Electronic Supplementary Materials. The maximum angle between the feature vectors can be calculated for each type of contour. The maximum included angle of the feature vectors in area 1 is 13.9697°. The maximum included angle of the feature vectors in area 2 is 15.82°. And the maximum included angle of the feature vectors in area 3 is 13.03°. Section 2.2 shows that when the laser incident angle is relatively small, the line width of laser engraving is basically the same. For a higher trimming quality, the maximum incident angle is limited to 5°, and the contours in areas 1-3 are partitioned. The maximum incident angle between the feature vectors of each cluster after the partition is less than 5°.   In order to improve the processing efficiency of surface contours with similar features, a partition method based on the k-means clustering algorithm is presented. This method groups the contours with similar machining characteristics. First, the position data and feature vectors of the trimming contours are extracted via CAM software. The feature vector is defined as the ideal direction of the laser axis when the contours are trimmed, which is manually specified according to the contour characteristics. Then, the k-means clustering algorithm is used to cluster these contours. k The k-means clustering algorithm is an iterative clustering algorithm that uses distance as a similarity measure. The algorithm sets the number of clusters in advance, and the clustering center is obtained based on the sample average within the cluster. The algorithm takes the error square sum function as the objective function, and the objective function minimization criterion as the basis for evaluating whether the clustering converges. In this paper, the contour feature vectors are used as the clustering data, and the objective function (J) can be written as Eq. (9): where is the number of feature vectors, is the preset number of categories, and is the six-dimensional data representing the feature vector, is the coordinate of the cluster center. For the ceramic core used in this paper, is 60. Since the data representing the geometric positions and the feature vectors have different scales, position data is normalized to obtain a more accurate clustering result. The essence of normalization is to assign data weights in different dimensions and to improve the clustering effects by adjusting those weights. Thus, the clustering effects under various weights can be further studied. In this paper, the clustering result via normalization operations is achieved. The formula for normalization is shown in Eq. (10): x max y max z max where , , and represent the maximum values of the feature vector data in the X-, Y-, and Z-axis  x iŷiẑî directions, respectively, and , , and represent normalized data. The normalized feature vector data is expressed as .

k k
For the k-means clustering, the desired result is that contours with close distances and similar feature vectors are grouped into one category, and the contours in the same category meet the requirement that the maximum included angle between feature vectors is less than 5°. Since areas 1-3 have their characteristics, and the maximum angles are all greater than 5°, each area must be divided into at least two regions, that is, the value must be greater than 6. Through verification, when is equal to 10, the ideal clustering result can be obtained. The normalized clustering result is shown in Fig. 11(a), and the non-normalized clustering effect is shown in Fig. 11(b).
According to Fig. 11, when the clustering data is not normalized, the clustering results tend to be based on the position data. It is due to the weight of the feature vector data being relatively insignificant compared to the weight of the position data in the objective function. As a result, the partition effect of the contour cluster is incorrect. However, the contour cluster can be correctly classified following the normalization process. In Table 1, the clustering centers and the maximum included angles between the feature vectors in each category are shown. All included angles are less than 5°.

3D contour layering algorithm and processing path generation
The clustering results of the ceramic core contours designated for trimming are provided in Section 2.4. In this paper, the 5 + 2 axis processing method is employed to trim the ceramic core. In other words, positioning is done via the five-axis machine tool, while simultaneously scanning with 2D high-speed galvanometer. Thus, it is necessary to further process the clustered 3D contours to meet the processing requirements of the 2D galvanometer. The specific steps are: (1) transformation of the contour data corresponding to the same cluster from the workpiece coordinate system (WCS) to the local coordinate system (LCS), (2) the layering operation of the contour data, and (3) the projection transformation from 3D data to 2D data. After clustering, the vector direction of the clustering center corresponds to the positioning direction of the five-axis machine tool during processing. In other words, it corresponds to the direction of the laser beam, i.e., the Z-axis direction in LCS, as demonstrated in Fig. 12. The contour data is transformed from the WCS to where parameter represents the data of the ith cluster in WCS, represents the data in LCS after the coordinate transformation of , and represents the transformation matrix from WCS to LCS . Here, i represents the index of the contour cluster.
The contour point data in the LCS can be determined by Eq. (11). , , and represent the unit vector along the coordinate axis in the ith LCS. Then , , and are defined as follows: is the feature vector of the cluster after contour clustering; is determined as the difference between the unit vector in the Z-axis direction in WCS and the projection in the Z-axis direction in LCS, as shown in Fig. 12. The parameter can be determined by the cross product of and , as shown in Eq. (12): (12) The LCS calculation results of some contour clusters are shown in Fig. 13, such as contour cluster 2, contour cluster 8, and contour cluster 10. The results of the LCS defined in the above manner are satisfactory.
The value of the unit vectors , , and can be obtained via Eq. (12). Then, the transformation matrix can be obtained via Eq. (13). The contour data in the WCS can be converted to the LCS by applying Eq. (11).
where represents the coordinate of the coordinate origin of LCS in WCS . After the contour data is transformed to LCS, layering and projection operations are required, as shown in Fig. 14. With being the direction of the Z-axis for LCS and simultaneously the direction of the contour layering, the layering operation can be executed based on the Z-axis coordinates of the data. z max z min First, the maximum and minimum coordinates and of the points in the current layering are obtained. Then, the contour data is divided into several layers along the direction with interval. When finished, the contours of each layer are still represented as 3D data. Since the galvanometer is 2D, it is necessary to project the 3D data in each layer along the direction. The 2D projected contour data required for the galvanometer can The positioning data of the five-axis machine tool corresponding to each layer can be obtained in the following way. Assuming that the clustering center of the ith cluster is expressed as , , then the position data of the five-axis machine tool corresponding to each layer can be calculated by Eq. (14): floor (x) x N layer where the function represents the downward rounding of the variable , is the layer number of the current cluster, represents the positioning coordinates of the mth layer of the ith contour cluster, and is the layering distance, the value is 0.3 mm.  3 Process and control of ceramic core trimming Galvanometer scanning and machine tool positioning data can be obtained by layering and projecting the 3D contours. The five-axis ultra-fast laser processing system employed within this paper is shown in Fig. 16. The system is integrated with a high-precision five-axis motion platform, a 2D scanning galvanometer, an ultrafast picosecond laser and a 3D structured light vision measurement module. The five-axis motion platform, the scanning galvanometer and the 3D visual scanner are high-precision equipment and instruments. The positioning accuracy and repeat positioning accuracy of the linear axes of the five-axis motion platform are higher than 5 µm, and the rotating axes are higher than 3″. The repeatability accuracy of the scanning galvanometer is higher than 1 µrad. The measurement accuracy of the 3D scanner is 15-8 µm.
In Fig. 17, the hardware integration scheme of the system is shown. The industrial computer, visual measurement computer, five-axis numerical control system and scanning galvanometer are connected through an Ethernet switch, and communication and data transmission between them is realized based on socket communication technology. The laser is connected with industrial control computer and numerical control system through RS232 serial communication technology. The components above are mutually combined through the OEMCC software. The software is developed via C# language, using WinForm interface library, OpenTK graphics display library, and point cloud processing library. During ceramic core trimming, the point cloud model of the ceramic core is first obtained via a 3D vision measurement device. Then, the standard template trajectory exported by CAM is registered with the actual point cloud model through the registration algorithm to obtain the registered 3D trajectory. According to the contour feature vectors, the contours are clustered, layered and projected to machine tool positioning data as well as the galvanometer scanning data. The ceramic core trimming process is completed by the OEMCC software. The specifics of the process are shown in Fig. 18.
The direction of the optical axis for machine tool positioning is the direction of the feature vector at the center of the ith cluster. In actual processing, when the machine tool locates the corresponding position of a specific layer of data, the high-speed galvanometer is used to scan the contours of the layer. Then the machine tool is located to the position of the next layer, and the Xiaodong WANG et al. Ceramic core trimming based on a zone-layered method galvanometer completes the scanning of the contours in the current layer again. After the contours of all layers in the current region have been scanned, the contour processing of the following area is performed.

Experimental validation
In this paper, online trimming of the ceramic core is achieved. Following the clamping of the ceramic core, the point cloud model is obtained via a stereo vision measuring device. The trimming trajectory following the position correction is obtained via the ICP registration algorithm, which registers the template trajectory and the point cloud model. The registration process between CAD extracted trajectory and the edges to be trimmed of the measured point cloud is shown in Fig. 19. Firstly, the template trajectory is extracted from the CAD model. Then, taking the points of the CAD model as the source point cloud and the measured point cloud of the ceramic core as the target point cloud, the coarse registration is carried out based on the ICP registration algorithm. The template trajectory is transformed with the obtained transformation matrix to realize the rough registration between the template trajectory and the edges of the measured point cloud. Then, taking the template trajectory points as the reference points, the point cloud data is clipped within a specific range of the reference points to obtain the sub-point clouds near the reference points in the source point cloud and the target point cloud, respectively. These sub-point clouds represent the feature contour on the ceramic core. The ceramic core is deformed compared with the standard CAD model, the coarse registration cannot achieve the complete alignment between the template trajectory and the edges to be trimmed, so it is necessary to further fine registration based on sub-point clouds. The registration in this step is based on the scale ICP registration algorithm, that is, the point cloud registration algorithm with scale scaling. The accurate registration of template trajectory and trimming edges can be realized through the fine registration. Finally, the fine registration operation is carried out for all sub-point clouds to obtain the modified trajectory that can be used for practical ceramic core trimming.
The detailed process is shown in Ref.
[25]. The visual ceramic core measurement is shown in Fig. 20(a). The visual measurement system is installed on the five-axis motion platform to measure the ceramic core from various angles, thus generating the point cloud model as shown in Fig. 20(b). Template trajectory results before and after the registration are depicted in Fig. 21. It can be observed that the post-registration modified trajectory is well aligned with the ceramic core edges. After obtaining the registered trajectory, the trajectory is partitioned and layered according to the method proposed in this paper. Both scanning data of the high-speed scanning galvanometer and the movement path data of the five-axis motion platform are obtained.
In the experiment, the picosecond laser is used. The  output laser wavelength is 532 nm, the repetition frequency is 1 Hz-2 MHz, the maximum output power is 60 W, the maximum single pulse energy is 30 µJ, the pulse width is 10 ps, and the beam quality is equal to M 2 = 1.1. After a triple beam expander expands the laser beam, it is transmitted through a series of reflectors, a galvanometer and a field lens to focus on the edge of the ceramic core designated for trimming. The optical system is shown in Fig. 22. Laser trajectory during ceramic core trimming is shown in Fig. 23. The laser scanning trajectory after layered processing of the contour cluster 2 is shown in Fig. 23(a). The laser scanning trajectory of the contour cluster 8 is shown in Fig. 23(c), and the laser scanning trajectory of    the contour cluster 10 is shown in Fig. 23(b). The results before and after ceramic core trimming are shown in Fig.  24. There are defects such as burrs and hole blockages in the characteristic structure of the ceramic core before trimming, as shown in Figs To guarantee the quality of trimming for manual trimming is challenging, and there are often phenomena such as uneven burr trimming and burr residue. Although the manual trimming method has been used in the ceramic core of aero-engine blades, in contrast, the effect of using the laser processing equipment and the methods proposed in this paper has been significantly improved, not only has a higher trimming accuracy and efficiency but also dramatically improves the yield.
The experimental results show that the method proposed in this paper can effectively realize the trimming of ceramic cores. The method makes full use of the characteristics of the galvanometer and the machine tool, that is, it avoids the vibration caused by the frequent start and stop in pure machine tool processing and makes full use of the high speed and high dynamic response characteristics of the galvanometer. This paper proposes a zonal and layering laser processing method of 3D surface contour clusters, and the experimental validation is carried out via ceramic core trimming. The overall benefits of the proposed method can be summarized as follows: (1) Based on the feature vector and the position data of the contours, the k-means clustering algorithm is utilized to partition the contour cluster. By normalizing the position data and the feature vector data, the problem of contour misclassification caused by using only position data is solved.
(2) To construct the LCS, the clustering center is employed as the origin point, and the vector direction of the clustering center is taken as the Z-axis direction. To Fig. 24 Trimming results of the ceramic core: (a) ceramic core before trimming, (b) burr and hole blockage of long groove, (c) burr and hole blockage of hole, (d) ceramic core after laser trimming, (e) long groove after laser trimming, (f) hole after laser trimming, (g) ceramic core after manual trimming, (h) long groove after manual trimming, and (i) hole after manual trimming.
obtain the scanning data of the galvanometer, the contour cluster data in the WCS is transformed into the LCS by the coordinate transformation algorithm. The definition method of the LCS greatly facilitates the subsequent layering and projection operations.
(3) In order to achieve the online ceramic core trimming, a laser processing system with the reconfiguration of the multi-axis motion platform and the galvanometer is implemented. A laser processing software OEMCC is employed. This software can realize functions such as vision measurement, motion control, trajectory registration and generation, and laser scanning. The software is suitable for laser processing systems with similar configurations.
The processing technology and equipment proposed within this paper are suitable for online laser trimming of workpieces with curved contour clusters, especially for non-metallic, difficult-to-process materials such as ceramic cores. In the future, further optimization of process parameters will be implemented to achieve better trimming quality.  Normalization of Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution, and reproduction in any medium or format as long as appropriate credit is given to the original author(s) and source, a link to the Creative Commons license is provided, and the changes made are indicated.

Nomenclature
The images or other third-party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
Visit http://creativecommons.org/licenses/by/4.0/ to view a copy of this license.