Variable Stiffness Identification and Configuration Optimization of Industrial Robots for Machining Tasks

Industrial robots are increasingly being used in machining tasks because of their high flexibility and intelligence. However, the low structural stiffness of a robot significantly affects its positional accuracy and the machining quality of its operation equipment. Studying robot stiffness characteristics and optimization methods is an effective method of improving the stiffness performance of a robot. Accordingly, aiming at the poor accuracy of stiffness modeling caused by approximating the stiffness of each joint as a constant, a variable stiffness identification method is proposed based on space gridding. Subsequently, a task-oriented axial stiffness evaluation index is proposed to quantitatively assess the stiffness performance in the machining direction. In addition, by analyzing the redundant kinematic characteristics of the robot machining system, a configuration optimization method is further developed to maximize the index. For numerous points or trajectory-processing tasks, a configuration smoothing strategy is proposed to rapidly acquire optimized configurations. Finally, experiments on a KR500 robot were conducted to verify the feasibility and validity of the proposed stiffness identification and configuration optimization methods.


Introduction
The application of industrial robots has significantly increased in automatic manufacturing owing to their remarkable advantages such as great operation flexibility, high intelligence, low cost, and low space requirement [1][2][3]. However, the stiffness of industrial robots is only 2% to 5% of that of computer numerical control (CNC) machines because of the inherent characteristics of their series structures [4][5][6]. Such weak stiffnesses of industrial robots significantly affect their positional accuracy and processing quality of products, which severely limits their application in the manufacturing and assembly of high value-added products. Therefore, the stiffness characteristics and strengthening methods must be extensively investigated to improve the operation performance of robot equipment.
Many researchers have investigated robot stiffness characteristics and modeling methods. Dumas and Caro proposed a fast identification method for joint stiffness and established the Cartesian stiffness model of the KUKA KR240-2 robot through experimental results [7]. Zhou et al. identified geometric parameters to revise the robot kinematics model and the Jacobian matrix. Based on this, a more precise joint stiffness matrix and stiffness model could be obtained [8]. Klimchik et al. adopted the method of a virtual axis to couple the connecting rod gravity and external load into the stiffness modeling process, which resulted in a higher accuracy of stiffness modeling [9]. The aforementioned stiffness identification and modeling methods were established based on the assumptions of rigid connecting rods and an elastic rotation axis of the joint. However, the gravity center

Open Access
Chinese Journal of Mechanical Engineering *Correspondence: tw_nj@nuaa.edu.cn of the connecting rod and the torsional stiffness of the motor and reducer vary with the operation configuration. In addition, different sampling configurations result in different identification results because of the inherent nonlinearity of joint stiffness [10]. In other words, a set of joint stiffness identification results cannot satisfy the accuracy requirements of the stiffness model in the entire workspace. Furthermore, the stiffness performance of the robot primarily depends on three factors: (1) the connecting rod structure and material properties, (2) stiffness of the actuator and transmission mechanism, and (3) working configuration. Compared with the first two factors, the optimization of the robot configuration has the advantages of high technical feasibility, good task adaptability, and no need to change the structure and control system of the robot, which indicates a high application value in engineering applications.
Sabourin et al. [11] developed a comprehensive optimization objective function combining kinematic and stiffness performance to optimize the configuration of a robot machining system. Xiong et al. [12] proposed a discrete search algorithm to configure the operation of a milling robot with optimal stiffness performance, and the optimization results significantly increased the trajectory accuracy. Guo et al. [13] used a robot stiffness ellipsoid as the evaluation index to improve robot stiffness performance through configuration optimization, and a higher drilling axial accuracy and countersink depth accuracy were obtained. However, these studies focused on configuration optimization for processing tasks with discrete positions. When many points or trajectory tasks are involved, such as thousands of target positions in drilling tasks of aerospace components, the optimization of individual target points significantly restricts the processing efficiency, which cannot satisfy the requirements of short-cycle manufacturing tasks.
Motivated by these observations, a variable stiffness identification method was investigated for an accurate modeling of industrial robots that overcomes the effects of configurations on stiffness modeling. In addition, a configuration smoothing method is proposed for the rapid acquisition of optimized configurations based on the geometric position features of target machining tasks.
The remainder of this paper is organized as follows. Section 2 presents a planning method for sampling points based on space gridding, and a variable stiffness identification method is proposed. In Section 3, a task-oriented stiffness representation model is established to estimate the stiffness performance in a particular direction, and a configuration optimization strategy is proposed to maximize the performance. Aiming at numerous points or trajectory tasks, Section 4 proposes a robot configuration smoothing method to achieve a rapid acquisition of optimized configurations. In Section 5, verification experiments of the identification results and configuration optimization methods are designed and discussed. Finally, Section 6 concludes the paper.

Identifying Variable Stiffness Based on Space Gridding
The Cartesian stiffness of a robot is directly dependent on the joint stiffness and robot configuration according to the traditional static stiffness model. However, joint stiffness identification and robot configuration are interrelated. The stiffness of each joint varies with the robot configuration because of the complex transmission structure, which primarily consists of a motor and a reducer [14]. However, different configurations result in changes in structural gravity centers and deformations in the kinematic parameters, which affect the identification accuracy of joint stiffness [15]. Thus, variable stiffness identification is an effective method of overcoming the effect of configuration on the stiffness identification results.

Stiffness Identification Principle
Based on the elastic deformation assumption, the relationship between the displacement of the end effector (EE) and the corresponding operating load can be expressed as where F is the generalized force vector acting on the EE, D is the displacement vector, C is the inverse matrix of the stiffness matrix K, i.e., the robot compliance matrix, and J is the robot Jacobian matrix. K θ is a diagonal matrix composed of the joint stiffness: where θ 1 to θ 6 represent the six joints of the industrial robot, and k θ 1 to k θ 6 denote the stiffness of each joint. The deformation caused by the working load is nonlinear with the stiffness of each joint. Thus, when identifying the joint stiffness in experiments, numerous sampling experiments are required. Multiple groups of static loads under different robot configurations should be applied, and the displacements at the end of the robot can be subsequently measured. Simplifying the right-hand side of Eq. (1), the relationship between D and joint compliance can be expressed as where F, A, and C θ can be defined as follows: If the rotation angle and corresponding stiffness of each joint are determined, a robotic stiffness model can be established. If K c1 is the joint stiffness in a particular configuration (C 1 ) and K c2 is the joint stiffness in a closed configuration (C 2 ), then where k 11 to k 16 are the joint stiffnesses in C 1 , k 21 to k 26 are the joint stiffnesses in C 2 , and E is the 2-norm of the difference between K c1 and K c2 . A positive number ξ approaching zero always exists, which satisfies E<ξ, with the angle variation △θ approaching zero. Therefore, the corresponding joint stiffnesses of the two close configurations can be considered highly similar [19,20].
Through the above analysis, the workspace of the robot can be divided, and the joint stiffness corresponding to each subspace can be considered to be approximately fixed values. Owing to the advantage of higher processing efficiency, the reachable space of the robot can be divided into a set of cubic grids (Figure 1(a)).
Eight vertices of a particular grid can be selected as sampling points to envelop a single grid. In addition, the center point of the grid can be selected as a sampling point to ensure that the stiffness parameter in the grid is closer to the actual value. Thus, nine target points (Tags 1 to 9) are used as sample positions to identify the joint stiffness corresponding to a particular grid space ( Figure 1(b)). In the initial sampling posture of a target position, the three axes of the tool coordinate system are consistent with the directions of the three intersecting edges of the grid, such as the sampling posture of Tag 5.  In addition, different sampling postures can be obtained by rotating a certain angle around any axis of the tool coordinate system. Figure 1(c) shows the generation of a new posture by rotating the y-axis of the tool coordinate system. Based on this, different sampling postures of the nine sampling positions can be obtained (Figure 1(b)) [21]. Moreover, the joint stiffness tends to be stable with more than ten groups of sampling configurations in a stiffness identification experiment [22][23][24]. Therefore, in each sampling position of a grid space, more than two configurations should be used to ensure the accuracy of the stiffness identification results.

Variable Stiffness Identification
By combining the static stiffness identification method and the space gridding principle, a variable stiffness identification method is innovatively developed. Sampling positions can be obtained based on space gridding, and the robot stiffness model, which is suitable for an entire workspace, can be built precisely and described as where j is the grid serial number, and K θ (j) is the joint stiffness corresponding to the grid space in which the target position of the current robot is located. Eq. (12) can be described as the variable-stiffness model. Compared with the traditional static stiffness model, the variable stiffness model can more accurately reflect the stiffness characteristics in different robot position intervals; thus, an accurate evaluation of the stiffness performance and precise prediction of load-induced positional error can be achieved. The steps of the variable stiffness identification method are shown in Figure 2.
According to the identification results, the joint stiffness of the grid space in which the target position is located can be obtained. However, when the location of a target point occurs in some special positions (Figure 3), the selecting which sampling points to calculate joint stiffness is a complex problem.
When located in a vertex, side, or plane that belongs to only one grid, the joint stiffness in this grid space can be used. In addition, a target position may also appear in the following special locations: (1) in a common vertex of two, four, or eight grids; (2) in a common line of two or four grids; (3) in the contact plane of two grids. When the grid size is sufficiently small, the aforementioned scenarios can be solved by averaging the stiffness of each joint corresponding to all grids containing the target point.
The proposed stiffness identification method is an innovation over improvement of the static stiffness model. Therefore, although this paper uses KUKA KR500 as a discussion case, the method can be utilized for other types of robots. Based on accurate joint stiffness identification and modeling, the evaluation and optimization methods of robot stiffness performance can be studied for the practical application of the stiffness model in engineering, such as machining stiffness enhancement and high-precision positioning control.

Task-oriented Axial Stiffness Performance Index
In the robot machining process, the working load causes the deformation of the EE, where the deformation in the radial direction of the tool affects the cutting positional accuracy (drilling, boring, etc.) or the trajectory accuracy (milling, grinding, etc.), and the tool axial deformation primarily affects the cutting surface quality (roughness, flatness, etc.). Therefore, the definition of the robot stiffness performance index should be adapted to the load distribution characteristics of the specific task of the robot. Based on the different dimensions of the elements, Eq. (1) can be transformed as follows: where f and m indicate the force and torque matrices, respectively, d and δ represent the translational and rotational displacement matrices, respectively, K fd and K fδ are the force-translational and force-rotational stiffness submatrices, respectively, and K md and K mδ are the torque-translational and torque-translational stiffness submatrices, respectively. The posture error caused by the cutting load is very small, and tracking the working posture of EE using the online detection method is difficult; therefore, the error can be ignored. Subsequently, In addition, the main key factor of machining quality is the linear displacement of the end of the robot caused by the cutting load, and the change in the end posture caused by the torque has a slight effect on the machining accuracy and quality [25,26]. Therefore, Eq. (14) can be simplified similarly as Applying a unit force f to the EE of a robot, we obtain Eq. (16) describes an ellipsoid that changes with the robot operation configuration (Figure 4), whose directions of principal axes are eigenvectors of K T fd K fd . The singular values of K T fd K fd , denoted as λ 1 , λ 2 , and λ 3 , which reflect the lengths of the ellipsoid's semi-axis, represent the square of the stiffness values in three directions, respectively.
During the robot operation, the direction of the machining load at the end of the robot is not consistent with the main axis of the robot stiffness ellipsoid. Therefore, the axial stiffness performance of robots (13) must be studied to realize an accurate evaluation of machining performance. In Figure 4, λ t1 and λ t2 represent the two half axes of the elliptical section of the machining plane in the envelope space of the stiffness ellipsoid, and λ d is the semi-axis, which is the normal vector of the elliptical section. The square roots of the three half-axis lengths represent the stiffness values of the three axes in the robot tool coordinate system [27,28]. The direction of λ d in the stiffness ellipsoid coordinate system can be defined as the unit vector [e x e y e z ], and the mathematical formula of the line of its location is expressed as follows: where t is the proportionality coefficient of the spatial line. Thus, the stiffness ellipsoid can be expressed according to the definition formula of the ellipsoid as where x, y and z represent the coordinates of the intersection point between the straight line shown in Eq. (17) and the stiffness ellipsoid. Combining Eqs. (17) and (18), t can be eliminated, and the coordinates of the intersection point can be obtained as follows: Subsequently, λ d and k x can be calculated as follows: The method can also be extended to calculate k y and k z , which can be expressed as follows: (17) x where [h x h y h z ] and [l x l y l z ] are the unit vectors of λ t1 and λ t2 , respectively. Using these three formulas, the threedimensional stiffness in the robot tool coordinate system can be accurately calculated, which provides a theoretical basis and evaluation standard for precision control and quality improvement of machining tasks.

Stiffness-Oriented Configuration Optimization Method
Redundant degrees of freedom (DOFs) increase the task adaptability of a robot operating system. Consequently, the tasks of the robot operating system can be divided into two subtask levels. The first-level subtask controls the EE to move to the initial position and posture to satisfy the requirements of the operation task. As a second-level subtask, configuration optimization improves machining performance by using redundant DOFs [29]. The 6-DOF industrial robot has a functional redundant DOF in the direction of the tool axis, whether for position processing tasks (drilling, boring, etc.) or trajectory processing tasks (milling, grinding, etc.). In other words, theoretically, an infinite number of poses of the EE can be used to execute machining tasks [30]. Therefore, the change step of the rotation angle of the EE can be set to reduce the amount of calculation in the optimization process and then promote the running speed of the optimization method.
The feasible posture of the EE is described in the Cartesian coordinate system. However, the stiffness performance index of the robot is evaluated in the joint space. Therefore, the inverse kinematics algorithm and the uniqueness principle must be used to solve the joint angles corresponding to the current robot configuration.
The specific steps of robot configuration optimization are as follows: (1) Extract the position coordinates P (p x , p y , p z ) of the target point from the process digital model, and plan the initial posture N(α, β, γ) in the offline programming software. The current configuration C 0 of the robot can be obtained using an inverse kinematics algorithm and the uniqueness principle of the inverse solution.
(2) Use the axial stiffness evaluation index to calculate the target stiffness value in the initial configuration, denoted as k.
(3) Rotate the EE along the tool axis with θ x (Figure 5), where − 180° ≤ θ x ≤ 180°, and define the variable step of θ x as �θ x = 10 • . Thus, θ x can be calculated as follows: (4) Examine the safety and rationality of the new configuration. If the posture of the EE is unreachable or interfered, the configuration is abandoned. Otherwise, calculate the current configuration C ′ and corresponding stiffness performance index k ′ , and update k through the larger value of k and k ′ . (5) Repeat steps 3 and 4 to obtain the rotation angle with the best axial stiffness, and the corresponding optimal configuration can be determined.
The workflow of configuration optimization is shown in Figure 6.

Smooth Processing Strategy
In the aviation manufacturing industry, some parts require numerous drilling tasks. If the operation configuration of each target machining position is optimized individually, the workload of the optimization processing is large and the efficiency is low. Therefore, for processing tasks with numerous discrete positions, according to the geometric characteristics of the product or distribution of these positions, considering the typical operation positions as the optimization objective, the optimization efficiency of the robot configuration in other target positions can be significantly increased using the interpolation smoothing method. Similarly, this method can be extended to milling and other trajectory processing tasks, which are controlled by interpolation points. The above optimization strategies are the extended applications of the configuration optimization method proposed in Section 3. The typical operation points are shown in Figure 7. In addition to the special scenario in which the workpiece surface is a plane, the direction of each axis in the tool coordinate system of each machining position cannot be guaranteed to be in the same direction or plane. Therefore, all target positions must be unified to the same reference coordinate system. In this paper, the base coordinate system is selected as the reference.
The rotation matrix of the posture of the EE can be expressed as an RPY angle: where Rot expresses the rotation matrix, γ i , β i , and α i are the corresponding rotation angles of Base i R in three directions, respectively. Because Base i R is an identity invertible matrix, Base i R x represents a coordinate system obtained by rotating the robot base coordinate system around its x-axis at a certain angle γ ′′ i , whose RPY angle can be expressed as Thus, φ x can be expressed as Therefore, the corresponding posture of the target position can be solved according to the proportional (28) Base i  Figure 9 shows the experimental platform of the study. A KUKA KR500 industrial robot was used as the operation carrier in the identification and machining tasks. An ATI IP60 Omega160 force transducer, which was fixed on the flange, was used to measure the working loads. An API laser tracker was used to establish coordinate systems and measure the positional errors of the robot.
calibration space, and 600 mm, 300 mm, and 150 mm were selected as the side lengths of the cubes to study the effects of different grid sizes. Therefore, the calibration space was divided into 2, 16, and 128 cubic grids, respectively. Based on this, the initial configuration of the samplings could be determined. When rotating the EE at ± 10° along the y-axis of the tool coordinate system, a total of 27 sampling configurations were obtained in a cubic space. The joint stiffness without space gridding was identified as Dividing the calibration space into two symmetrical cubic grids (Figure 10(a)), the joint stiffness in two grids, namely Grids 1 and 2, were identified as By selecting a 300 mm cubic grid as a sampling unit, the calibration space was divided into 16 grids (Figure 10(b)). The entire space was observed as four cuboid spaces: Cuboids 1 to 4, whose longest sides were parallel to the y-axis of the base coordinate system. The joint stiffness of each grid space was identified, and the (38) K θ = diag(1.58 × 10 10 , 6.12 × 10 9 , 5.28 × 10 9 , 4.66 × 10 8 , 2.19 × 10 8 , 3.49 × 10 8 ) N · mm/rad. fluctuation of values was observed according to these four cuboid spaces ( Figure 11). A 150 mm cubic grid was selected as a sampling unit and the calibration space was divided into 128 grid spaces (Figure 10(c)). The entire space was divided into four cuboid spaces, i.e., Cuboids 1 to 4, which could be referred to as first-order cuboid spaces. Each first-order cuboid space was further split into four second-order cuboid spaces, such as Cuboids 1.1 to 1.4, which were formed with eight grid spaces. The joint stiffness of each grid space could be calculated, and the variation regularity of the stiffness value could be investigated based on the second-order cuboid spaces (Figures 12,13,14,15).
Analyzing the variation trend of the joint stiffness in Figures 11 and 12, 13, 14, 15, the stiffness of the second, third, and fifth joints largely maintained the same tendency by dividing the robot calibration space into smaller grid spaces. However, the stiffness of the first, fourth, and sixth joints exhibited different change trends, compared with the tendency with larger grid spaces. Thus, the variation trend with larger grids cannot indicate the local In addition, the joints whose axes were parallel to the y-direction of the base coordinate system, namely the second, third, and fifth joints, exhibited stable change trends of stiffness in the entire calibration space. However, the joints whose axes were perpendicular to the y-direction of the base coordinate system, namely the first, fourth, and sixth joints, had similar stiffness tendencies in the calibration space. We can conclude that the axis direction of a robot joint is related to the stiffness distribution of the joint, which can be used to optimize the machining configuration of industrial robots.
The effectiveness of the variable stiffness identification method was verified through compensation experiments of the load-induced positional error. In the calibration space, a sampling position was randomly selected in each 150 mm grid. Thus, 128 verification points were obtained. A load of 50 kg was fixed to the EE to produce position errors.
The load-induced positional errors were measured by executing control commands before and after loading. According to Eqs. (1) and (15), the load-induced positional errors were calculated as follows: where E is a matrix composed of the first three columns of the first three rows of the objective matrix. By combining the working load and the identified joint stiffness, the load-induced positional errors were predicted, and the control commands were modified through reverse compensation. Subsequently, the compensation effect of load-induced errors was evaluated by executing the modified commands under the loaded state. Defining E x (i), E y (i), and E z (i) as the load-induced positional errors of the ith position on the x-, y-, and z-axes of the base coordinate system, the load-induced absolute errors E(i) were calculated using the following equation: The absolute positional errors before and after compensation with different grid sizes are shown in Figure 16, and Figure 17 shows the error distribution for different grids. The average absolute positional error induced by the working load was 0.2868 mm, and the maximum error In conclusion, the validity of the identification results was verified, and the variable parameter error could better characterize the error model of a robot and obtain better positional accuracy.

Configuration Optimization and Smooth Processing Experiments
The operation configuration optimization strategy and smooth processing method were verified through the simulation layout of a robot operation system in the DELMIA software. Using a workblank fixed on the tooling as the processing object, the machining path was planned in the software ( Figure 18). The path from Tag 1 to Tag 2 was parallel to the y-axis of the base coordinate system, and the path from Tag 3 to Tag 4 also satisfied this scenario. According to the requirements of not using singular and joint-limit configurations, the rotation angle range of the robot EE was selected as − 90° ≤ θ x ≤ 90°.
In the verification experiments, the normal stiffness of the robot motion trajectory in its work plane, defined as k v , and the tool axial stiffness k x , which were directly related to the trajectory accuracy and surface cutting quality, respectively, were considered as the evaluation indexes. The normal stiffness k v was calculated using the following equation: Thus, the fluctuations of k x and k v with different rotation angles (θ x ) in Tags 1 to 4 are shown in Figures 19,20,21,22. Specifically, all the rotation angles corresponding to the configuration with optimal axial stiffness in Tags 1 to 4 were θ x = 0°, that is, the initial robot configuration. The fluctuation of k v in Tags 1 and 4 formed a trough with θ x = 0°, and two peaks were formed with θ x = − 20° and θ x = 10°. The optimal k v was observed at θ x   Tags 2 and 3 at θ x = − 10° and θ x = 20°, respectively, and the optimal k v was observed at θ x =20°.
With the rotation angle farther from θ x corresponding to the configuration with optimal stiffness performance, the stiffness in the target direction gradually decreased. Note that when the EE rotated away from the base coordinate system, the stiffness value decreased gradually. Figures 23 and 24 compare k x and k v at different positions, respectively. The four positions were symmetrically distributed along the y-axis of the robot base coordinate system, and the stiffness change curves at the corresponding positions also exhibited spatial symmetry along the y-axis of the base coordinate system. In addition, the stiffness value exhibited a downward trend along the positive x-axis direction of the base coordinate system. In other words, under the same operational posture of the EE, k x or k v in Tag 1 was lower than the corresponding directional stiffness in Tag 3, and k x or k v in Tag 2 was lower than the corresponding directional stiffness in Tag 4.
To increase the trajectory accuracy in the operation process, the operation configuration with optimal k v was selected first. Based on the optimized configuration of Tags 1 to 4, the smoothing processing of the trajectory interpolation point was conducted to achieve the rapid acquisition of the optimal stiffness configuration of corresponding positions and smoothing results ( Figure 25). The smooth processing strategy simplified the configuration optimization of all target positions into several geometric feature positions, which significantly increased the configuration optimization efficiency. The optimization effect of smooth processing was evaluated by comparing it with the optimal stiffness performance obtained from the configuration optimization of the interpolation points. The interpolation positions of Tags 3 and 4 were selected, and the seven interpolation positions were defined as Tags 5 to 11. The k v obtained using the configuration optimization method and smooth processing are shown in Table 1, and Figure 26 shows a comparison between the smooth processing and configuration optimization.
As Table 1 shows, the stiffness index k v after smooth processing was highly similar to the result of configuration optimization, and the stiffness loss after fairing was less than 0.409% (in Tag 8). The stiffness of the target direction in Tags 5 and 11 after smoothing was slightly higher than that of configuration optimization. This was because the step size of the rotation angle was larger in the process of optimization, which resulted in the optimized configuration obtained by smoothing closer to the configuration with optimal stiffness.
In summary, smooth processing effectively increases the optimization efficiency of a robot operating configuration and ensures that the improvement effect of the axial stiffness through smooth processing is not significantly reduced compared with that of configuration optimization, which fully proves the effectiveness of the fairing method.

Machining Experiment
The effects of configuration optimization and smoothing on milling quality were studied using a robot machining system. In the verification experiment, the milling task was performed in the cylinder head of an automobile engine. The planned trajectory in the machining experiment is shown in Figure 27, where Tags 1 and 2 represent the start and end positions of the trajectory, respectively.
By selecting the tool axial stiffness as the optimization objective, the rotation angles corresponding to the start and end positions of the machining trajectory were obtained using the configuration optimization method; the optimized configurations of Tags 1 and 2 are shown in Table 2, and smooth processing of the interpolation positions was conducted to obtain optimal configurations. Based on this, the trajectory control program of the robot was generated through offline programming software, and the robot could be driven to perform milling tasks. The process parameters of the robot milling experiment are listed in Table 3.
The comparison results of the milling process after configuration optimization are shown in Figure 28. The milling surface quality before configuration optimization was relatively poor, and the blades severely vibrated. The roughness of the surface was R a 2.356, as measured using a roughness meter (Sanfeng SJ-210). After optimizing the operation configuration of the robot to increase the end operation rigidity, the obtained milling surface of the workpiece was relatively smooth, the surface roughness was R a 0.597, and the accuracy and quality of milling were significantly improved.   Tag3  Tag5  Tag6 Tag7  Tag8  Tag9 Tag10 Tag11 Tag4 Configuration smooth processing Configuration optimization one by one