Integrated Design of Process-Tolerance for Remanufacturing Based on Failure Feature

: The uncertainty failure of the used part leads the complexity selection of remanufacturing processes. The different remanufacturing process combinations among used parts of used products also make the formulation of tolerance schemes more difﬁcult. It is hard to guarantee the optimality of process-tolerance schemes by traditional serial production modes in general, in which tolerance design is followed by process formulation. In order to generate the optimal remanufacturing scheme of process and tolerance for used products, an optimization method to integrate designs of process-tolerance (IDP-T) based on fault features was presented. In this work, the failure description set of used parts was constructed by combining the attribute characteristics and failure characteristics. Case-based reasoning (CBR) was ﬁrst utilized to generate the feasible remanufacturing process plans of used parts. Then, based on the feasible process plans, the factors of cost, quality loss, closed-loop accuracy and machining ability of remanufacturing were comprehensively considered to construct the optimization model of IDP-T. The Beetle Antennae Search algorithm (BAS) was used for the optimal alternative selection. Finally, a used gearbox was taken as an example to illustrate the validity and practicality of the proposed method. The results showed that the proposed method was effective in the optimization of IDP-T for remanufacturing.


Introduction
With the development of the scale of global industrialization and the shortening of the technological iteration cycle, the elimination quantity and speed of mechanical products increases. This causes countries to pay more attention to the circular economy. Remanufacturing as a technology to restore the performance, and finally reach the "as new" state of used products, can fully release the residual value of waste products [1]. It has become an important part of the circular economy [2,3]. Compared to new product manufacturing, remanufacturing can reduce 60% of energy, 70% of materials, 50% of cost and 80% of air pollution while achieving the same product quality [4]. It is considered to be a technology with both environmental and commercial value [5,6]. For the great advantage to be gained from remanufacturing, to a large degree this depends on the remanufacturing process and tolerance scheme.
The process and tolerance scheme formulated for used parts was a critical step in remanufacturing of the remanufactured products, which established a key connection between used parts and remanufactured products. As a kind of mechanical processing, an optimal process and tolerance scheme of remanufacturing created more benefits in optimization model based on feasible process plans was constructed. In Section 3, the effectiveness and feasibility of the proposed method was verified through the case study of a used gearbox. Finally, Section 4 summarized the main conclusions and future research points.

The Framework of IDP-T for Remanufacturing
IDP-T aims to generate the process and tolerance schemes simultaneously in the design stage. The optimization is considered as the pathways to obtain the optimum scheme. The process is shown in Figure 1.

Feasible process plans
Multi-objective function construction CBR   This optimization framework consists of three steps. 1. Feasible process plan generation. In order to obtain the feasible remanufacturing process plans of the failure surface, FTA (fault tree analysis) and FCE (fuzzy comprehensive evaluation) are used to extract and quantify the failure features of the failure surface, and combine the attributes of the failure surface to construct the failure surface feature description set. Then, CBR (Case-based reasoning) is utilized to match the case base to generate the feasible process plans of the failure surface. 2. Multi-objective function construction. IDP-T is a multi-objective optimization problem. In this step, the minimization of the remanufacturing cost and quality loss of used products are regarded as optimization objectives, and the processing capacity and assembly precision requirements as constraints to construct the optimization model. 3. Optimal scheme generation utilizing BAS. BAS is used to solve the problem of multi-objective optimization. At last, the optimum value is obtained, and the optimal plan of IDP-T is selected

Feasible Process plans generation
With the continuous accumulation of remanufacturing production experience, enterprises have accumulated amounts of similar cases in historical production activities, and there are usually multiple feasible process plans for the same failure surface. Therefore, how to reasonably describe the failure surface features and efficiently retrieve similar cases from the case base is the key step of IDP-T. This optimization framework consists of three steps.

1.
Feasible process plan generation. In order to obtain the feasible remanufacturing process plans of the failure surface, FTA (fault tree analysis) and FCE (fuzzy comprehensive evaluation) are used to extract and quantify the failure features of the failure surface, and combine the attributes of the failure surface to construct the failure surface feature description set. Then, CBR (Case-based reasoning) is utilized to match the case base to generate the feasible process plans of the failure surface.

2.
Multi-objective function construction. IDP-T is a multi-objective optimization problem. In this step, the minimization of the remanufacturing cost and quality loss of used products are regarded as optimization objectives, and the processing capacity and assembly precision requirements as constraints to construct the optimization model.

3.
Optimal scheme generation utilizing BAS. BAS is used to solve the problem of multiobjective optimization. At last, the optimum value is obtained, and the optimal plan of IDP-T is selected.

Feasible Process Plans Generation
With the continuous accumulation of remanufacturing production experience, enterprises have accumulated amounts of similar cases in historical production activities, and there are usually multiple feasible process plans for the same failure surface. Therefore, how to reasonably describe the failure surface features and efficiently retrieve similar cases from the case base is the key step of IDP-T.
The characteristics set of the failure surface i was constructed as a kind of data structure F i = { f i , a i }, which is comprised of failure characteristics set { f i } and attribute feature set {a i } of failure surface i. In this data structure f i = {t i , q i } represents failure type t i and failure quantification q i of surface i; a i = {m i , s i , d i } represents material m i , shape s i and dimension d i of surface i. Once the characteristics set of the failure surface is obtained, the feasible process plans can be generation by CBR. The specific construction process is shown in Figure 2. The different working conditions of products cause the uncertainty of parts failure types, which affects the processes and parameters formulation. In addition, various fault types also exist on different surfaces of a used part. Therefore, failure types extracted from surfaces should be carried out first. The FTA is used in this procedure.
FTA is a deductive failure analysis method from top to bottom, which uses Boolean Logic to combine low-order events to analyze failure scenarios in the system [20]. In this paper, FTA was utilized to identify the failure type of surface i, as shown in Figure 2.  In Figure 2, " "in (a) is the input events, representing the main failure performance of the used parts of this type of used product; " " and" "in (b) is the input events and condition events, representing the specific failure type of used parts and the cause of the failure type, respectively; " " in (c) is the output events, representing the specific failure parts and failure types of the used product; " "and" "are AND gate events and OR gate events, respectively.
When the failure type of surface i was obtained, the damage volume of the failure could be confirmed by the CAD model. Traditionally, the damage degree of failure was often evaluated by "high", "medium" and "low" attributes, which was a fuzzy evaluation method based on a designers' subjective interval. In order to overcome the shortcomings of this evaluation mode, the FCE was used to quantify the failure characteristics unifying the evaluation interval according to the specific situation of remanufactured

Failure Characteristics Set Construction
The different working conditions of products cause the uncertainty of parts failure types, which affects the processes and parameters formulation. In addition, various fault types also exist on different surfaces of a used part. Therefore, failure types extracted from surfaces should be carried out first. The FTA is used in this procedure.
FTA is a deductive failure analysis method from top to bottom, which uses Boolean Logic to combine low-order events to analyze failure scenarios in the system [20]. In this paper, FTA was utilized to identify the failure type of surface i, as shown in Figure 2.
In Figure 2 shape si and dimension di of surface i. Once the characteristics set of the failure surface is obtained, the feasible process plans can be generation by CBR. The specific construction process is shown in Figure 2.
The different working conditions of products cause the uncertainty of parts failure types, which affects the processes and parameters formulation. In addition, various fault types also exist on different surfaces of a used part. Therefore, failure types extracted from surfaces should be carried out first. The FTA is used in this procedure.
FTA is a deductive failure analysis method from top to bottom, which uses Boolean Logic to combine low-order events to analyze failure scenarios in the system [20]. In this paper, FTA was utilized to identify the failure type of surface i, as shown in Figure 2. In Figure 2, " "in (a) is the input events, representing the main failure performance of the used parts of this type of used product; " " and" "in (b) is the input events and condition events, representing the specific failure type of used parts and the cause of the failure type, respectively; " " in (c) is the output events, representing the specific failure parts and failure types of the used product; " "and" "are AND gate events and OR gate events, respectively.
When the failure type of surface i was obtained, the damage volume of the failure could be confirmed by the CAD model. Traditionally, the damage degree of failure was often evaluated by "high", "medium" and "low" attributes, which was a fuzzy evaluation method based on a designers' subjective interval. In order to overcome the shortcomings of this evaluation mode, the FCE was used to quantify the failure characteristics unifying the evaluation interval according to the specific situation of remanufactured and dimension di of surface i. Once the characteristics set of t obtained, the feasible process plans can be generation by CBR. The s process is shown in Figure 2.

Failure Characteristics Set Construction
The different working conditions of products cause the uncerta types, which affects the processes and parameters formulation. In ad types also exist on different surfaces of a used part. Therefore, fail from surfaces should be carried out first. The FTA is used in this proc FTA is a deductive failure analysis method from top to bottom, Logic to combine low-order events to analyze failure scenarios in the paper, FTA was utilized to identify the failure type of surface i, as sho In Figure 2, " "in (a) is the input events, representing the m mance of the used parts of this type of used product; " " and" events and condition events, representing the specific failure type of cause of the failure type, respectively; " " in (c) is the output even specific failure parts and failure types of the used product; " "and events and OR gate events, respectively.
When the failure type of surface i was obtained, the damage vo could be confirmed by the CAD model. Traditionally, the damage d often evaluated by "high", "medium" and "low" attributes, which w tion method based on a designers' subjective interval. In order to o comings of this evaluation mode, the FCE was used to quantify the fa unifying the evaluation interval according to the specific situation shape si and dimension di of surface i. Once the characteristics set of the failure surface is obtained, the feasible process plans can be generation by CBR. The specific construction process is shown in Figure 2.
The different working conditions of products cause the uncertainty of parts failure types, which affects the processes and parameters formulation. In addition, various fault types also exist on different surfaces of a used part. Therefore, failure types extracted from surfaces should be carried out first. The FTA is used in this procedure.
FTA is a deductive failure analysis method from top to bottom, which uses Boolean Logic to combine low-order events to analyze failure scenarios in the system [20]. In this paper, FTA was utilized to identify the failure type of surface i, as shown in Figure 2. In Figure 2, " "in (a) is the input events, representing the main failure performance of the used parts of this type of used product; " " and" "in (b) is the input events and condition events, representing the specific failure type of used parts and the cause of the failure type, respectively; " " in (c) is the output events, representing the specific failure parts and failure types of the used product; " "and" "are AND gate events and OR gate events, respectively.
When the failure type of surface i was obtained, the damage volume of the failure could be confirmed by the CAD model. Traditionally, the damage degree of failure was often evaluated by "high", "medium" and "low" attributes, which was a fuzzy evaluation method based on a designers' subjective interval. In order to overcome the shortcomings of this evaluation mode, the FCE was used to quantify the failure characteristics unifying the evaluation interval according to the specific situation of remanufactured " in (b) is the input events and condition events, representing the specific failure type of used parts and the cause of the failure type, respectively; " and dimension di of surface i. Once the characteristics set of the failure surfa obtained, the feasible process plans can be generation by CBR. The specific construc process is shown in Figure 2.

Failure Characteristics Set Construction
The different working conditions of products cause the uncertainty of parts fai types, which affects the processes and parameters formulation. In addition, various f types also exist on different surfaces of a used part. Therefore, failure types extra from surfaces should be carried out first. The FTA is used in this procedure.
FTA is a deductive failure analysis method from top to bottom, which uses Boo Logic to combine low-order events to analyze failure scenarios in the system [20]. In paper, FTA was utilized to identify the failure type of surface i, as shown in Figure 2. In Figure 2, " "in (a) is the input events, representing the main failure per mance of the used parts of this type of used product; " " and" "in (b) is the in events and condition events, representing the specific failure type of used parts and cause of the failure type, respectively; " " in (c) is the output events, representing specific failure parts and failure types of the used product; " "and" "are AND events and OR gate events, respectively.
When the failure type of surface i was obtained, the damage volume of the fai could be confirmed by the CAD model. Traditionally, the damage degree of failure often evaluated by "high", "medium" and "low" attributes, which was a fuzzy eva tion method based on a designers' subjective interval. In order to overcome the sh comings of this evaluation mode, the FCE was used to quantify the failure characteri unifying the evaluation interval according to the specific situation of remanufactu shape si and dimension di of surface i. Once the characteristics set of the failure surface is obtained, the feasible process plans can be generation by CBR. The specific construction process is shown in Figure 2.
The different working conditions of products cause the uncertainty of parts failure types, which affects the processes and parameters formulation. In addition, various fault types also exist on different surfaces of a used part. Therefore, failure types extracted from surfaces should be carried out first. The FTA is used in this procedure.
FTA is a deductive failure analysis method from top to bottom, which uses Boolean Logic to combine low-order events to analyze failure scenarios in the system [20]. In this paper, FTA was utilized to identify the failure type of surface i, as shown in Figure 2. In Figure 2, " "in (a) is the input events, representing the main failure performance of the used parts of this type of used product; " " and" "in (b) is the input events and condition events, representing the specific failure type of used parts and the cause of the failure type, respectively; " " in (c) is the output events, representing the specific failure parts and failure types of the used product; " "and" "are AND gate events and OR gate events, respectively.
When the failure type of surface i was obtained, the damage volume of the failure could be confirmed by the CAD model. Traditionally, the damage degree of failure was often evaluated by "high", "medium" and "low" attributes, which was a fuzzy evaluation method based on a designers' subjective interval. In order to overcome the shortcomings of this evaluation mode, the FCE was used to quantify the failure characteristics unifying the evaluation interval according to the specific situation of remanufactured shape si and dimension di of surface i. Once the characteristics set of the failure surface is obtained, the feasible process plans can be generation by CBR. The specific construction process is shown in Figure 2.
The different working conditions of products cause the uncertainty of parts failure types, which affects the processes and parameters formulation. In addition, various fault types also exist on different surfaces of a used part. Therefore, failure types extracted from surfaces should be carried out first. The FTA is used in this procedure.
FTA is a deductive failure analysis method from top to bottom, which uses Boolean Logic to combine low-order events to analyze failure scenarios in the system [20]. In this paper, FTA was utilized to identify the failure type of surface i, as shown in Figure 2. In Figure 2, " "in (a) is the input events, representing the main failure performance of the used parts of this type of used product; " " and" "in (b) is the input events and condition events, representing the specific failure type of used parts and the cause of the failure type, respectively; " " in (c) is the output events, representing the specific failure parts and failure types of the used product; " "and" "are AND gate events and OR gate events, respectively.
When the failure type of surface i was obtained, the damage volume of the failure could be confirmed by the CAD model. Traditionally, the damage degree of failure was often evaluated by "high", "medium" and "low" attributes, which was a fuzzy evaluation method based on a designers' subjective interval. In order to overcome the shortcomings of this evaluation mode, the FCE was used to quantify the failure characteristics unifying the evaluation interval according to the specific situation of remanufactured " are AND gate events and OR gate events, respectively.
Failure Quantification (q i ) When the failure type of surface i was obtained, the damage volume of the failure could be confirmed by the CAD model. Traditionally, the damage degree of failure was often evaluated by "high", "medium" and "low" attributes, which was a fuzzy evaluation method based on a designers' subjective interval. In order to overcome the shortcomings of this evaluation mode, the FCE was used to quantify the failure characteristics unifying the evaluation interval according to the specific situation of remanufactured products. A series of quantitative equations were then used to map the damage volume to the score interval of [0-10]. Once the type and damage volume of the failure was obtained, the failure could be quantified by the corresponding quantitative equation.
A failure quantification standard for a used gearbox was taken as an example, as shown in Table 1.
The used parts were affected by attributive character such as material, shape and dimension of the failure surface when making the remanufacturing process plan. Therefore, in order to establish the feasible process plans more accurately, it was necessary to describe the attributes of the failure surface after its failure type, and quantification was obtained. The attribute feature set of failure surface i was established as Finally, the characteristics set of failure surface i could be described as Feasible Process Plans Generation Based on CBR CBR realizes retrieval by calculating the similarity between new cases and existing cases. The similarity can be divided into local similarity and global similarity, where the local similarity is the similarity between two attributes, and the global similarity is the similarity between two cases. The relationship between local similarity and global similarity is as Equation (1).
where N is the new case and E is the existing case; n is the number of influence factors; c N i and c E i are the ith factors of new case and existing case respectively; sim(c N i , c E i ) is the local similarity of ith factors between new case and existing case; ω i is the weight of ith factors, which obtained by experts through AHP, the specific process can be referred to [21].
For local similarity, different data types require different calculation methods. Therefore, data of matching set G i are divided into three types: numerical data, enumerated data, and character data, and their corresponding similarity calculation models are established.
Numerical data consist of a series of continuous values in an interval (such as failure quantification). The similarity calculation model of numerical data is given by Equation (2). where max(c i ) and min(c i ) are the maximum and the minimum values that can be obtained in all cases of c i , respectively.
Enumerated data consist of a series of discrete values in an interval (such as surface dimension). The size of the failure surface (very small, small, medium, large, very large) can be represented by discrete data (0.1, 0.3, 0.5, 0.7, 0.9, respectively). The similarity calculation model of enumerated data is given by Equation (3).
where H i is the maximum discrete values that can be obtained in all cases of c i .

3.
Character data similarity calculation model.
Character data calculate the local similarity through the complete matching of keywords (such as surface shape, material, etc.). When the keywords are matched, the local similarity is 1, otherwise the local similarity is 0. The similarity calculation model of enumerated data is given by Equation (4).
After the global similarity is obtained through Equation (1), the similarity of case base instances is screened and matched by a preset similarity threshold, and the remanufacturing process corresponding to the case with the global similarity that meets the threshold is output. Then, the expert determines whether the feasible plans need to be revised, and the new plans are imported to upgrade and expand the case base. If there is no case that meets the similarity threshold, the process plans will be formulated by the expert and then imported into the case base.
The framework of feasible process plans generation of failure surface based on CBR is shown in Figure 3.

Multi-Objective Function Construction
IDP-T for remanufacturing needs to comprehensively consider the information of design, manufacturing, and assembly stage. Therefore, after obtaining the feasible process plans of each feature surface, the information of manufacturing stage and assembly stage should be quantified to consider in-scheme formulation, and a multi-objective optimization model needs be used to select an optimum scheme.
Remanufacturing cost and quality loss in IDP-T are considered to be the variables in multi-objective functions. The combination of different process-tolerance schemes of each failure surface will lead to different remanufacturing costs, which determine the profit and investment willingness of the company; the variety of tolerance leads to different

Multi-Objective Function Construction
IDP-T for remanufacturing needs to comprehensively consider the information of design, manufacturing, and assembly stage. Therefore, after obtaining the feasible process plans of each feature surface, the information of manufacturing stage and assembly stage should be quantified to consider in-scheme formulation, and a multi-objective optimization model needs be used to select an optimum scheme.
Remanufacturing cost and quality loss in IDP-T are considered to be the variables in multi-objective functions. The combination of different process-tolerance schemes of each failure surface will lead to different remanufacturing costs, which determine the profit and investment willingness of the company; the variety of tolerance leads to different quality losses, which will affect the subsequent processing of remanufacturing. As a result, the optimization model of IDP-T for remanufacturing is established, in which minimizing cost and quality loss is the objective, and is constrained by the requirements of assembly precision and the processing ability of the company. Before introducing the optimization model, some assumptions are provided as follows: 1.
The dimensional chain is in plane.

2.
Tolerances of remove process are symmetrically distributed.

3.
The failure surface needs to be processed as a whole. Partially repaired failure surface needs to coincide with the base plane after remanufacturing, so it can be regarded as a new part and the tolerances conform to the original dimensions.

Objective Functions Cost Function
The cost of IDP-T for remanufacturing includes the cost of machine tools, cutting tools and additional materials; remanufacturing process of failure surface can be divided into removal processes and additive processes. When there are n processes in the remanufacturing of failure surface i, the total cost is shown in Equation (5).
where C ij − and C ij + are the cost of removal process and additive process of surface i using process j, respectively; a and b are the numbers of removal process and additive process of surface i, respectively. The cost of removal process C ij − will increase with the improvement of the machining accuracy. The negative square model is used to construct the tolerance-cost function of the removal process, as in Equation (6).
where a 0ij − and a 1ij − are the fixed cost coefficient and tolerance cost coefficient of surface i using removal process j, respectively; T ij− is the tolerance of surface i using removal process j. The values of a 0ij − and a 1ij − can be obtained by fitting the enterprise costtolerance statistical samples by the least square method. For additive processes, different equipment and parameters will lead to greater difference in precision. According to literature [22,23], for the same equipment, when the main parameters remain stable, the tolerance is relatively fixed in unit thickness and meets the normal distribution. Therefore, when economic parameters are adopted, the additive process cost can be regarded as a fixed value related to the thickness of the additive, as shown in Equation (7).
where a 0ij + and f h ij + are fixed cost, and cost functions are related to the thickness of the additive when dimensions are restored under economic processing parameters, respectively. When the used product, having N failure surface needs remanufacturing tolerance design, the total cost is shown in Equation (8).

Quality Loss Function
Quality loss is an economic loss due to a less than expected performance of the product. According to Taguchi theory, size deviation of the parts will cause quality loss and ultimately affect the subsequent processing. Exponential equations are used to construct the quality loss function of remanufacturing processes, shown as Equations (9) and (10).
where K ij is the quality loss coefficient of surface i using process j; A 0ij and A ij are the original design dimension and remanufacturing actual machining dimension of surface i using process j, respectively; B ij and ∆ ij are the lost cost by out of tolerance and allowance tolerance of surface i using process j, respectively. In process level, the physical meaning of quality loss is the economic loss of subsequent process and service caused by the tolerance of the surface i using process j.
When the tolerance of removal process T ij − is symmetrically distributed, the quality loss of surface i using process j is shown as Equation (11).
For the additive process, quality loss varies greatly with different equipment; the tolerance is regarded as a fixed value related to the thickness of the additive. Therefore, the quality loss of additive processes is shown in Equations (12) and (13).
where σ ij + is the one-way maximum dimensional deviation coefficient per unit thickness with economic machining parameters of surface i using additive process j. Therefore, the total quality loss of the used product is shown as Equation (14).

Constraints Functions
IDP-T should be conformed to the actual machining capacity of the remanufacturing companies and assembly precision requirements of the product. Based on this, the constraint function should be established.

Tolerance Constraint of Closed-Loop
The tolerance of a closed-loop reflects the precision of the manufacturing stage and the assembly stage of the dimension chain. For remanufacturing the dimension chain, the tolerance of the closed-loop should meet the requirements of the closed-loop tolerance of the original dimension chain. However, due to the uncertainty of the part reused mode of remanufacturing the dimension chain, the original tolerance recovery of a single part cannot meet the requirements. It is necessary to include this uncertainty in the integrated tolerance design of remanufacturing.
The reused mode of used parts can be divided into new parts, reused parts, and remanufactured parts. The tolerance of the new part meets the requirements of new product suppliers; the reused parts had absolute dimensions, and the tolerance is 0; the tolerance of remanufactured parts needs to be optimized. Therefore, when considering different part reused modes of remanufacturing dimension chains, the tolerance constraint of a closed-loop is as in Equation (15).
where T i and T new are the tolerance of the remanufacturing surface i and new parts, respectively; ξ i and ξ new are the transfer coefficient of surface i and new parts, respectively; T 0 is the closed-loop tolerance of the original dimension chain.

Constraint of Machining Capability
There are differences in machining capability among different remanufacturing companies, and the actual machining capability of companies needs to satisfy Equation (16).
where T minij − and T maxij − are the minimum and the maximum machining precisions of the removal process of the company in economic processing parameters.

Multi-Objective Function Transformation
The multi-objective function is utilized to select the best scheme. The IDP-T optimization problem hinges on two objective functions. The objective optimization of IDP-T is to minimize cost ( f 1 ) and minimize quality loss ( f 2 ), shown in Equation (17). Since the disparate units in objectives, a weighted mean method is proposed to simplify this process, shown in Equation (18). The constraint quantization is summarized in Equation (19).

Optimization Implementation Process of BAS
IDP-T needs to comprehensively consider cost and quality loss in the design stage, while being constrained by closed-loop accuracy and machining capacity, which is a typical multi-objective and multi-constraint optimization problem. The general method to deal with this kind of problem is to transform it into a single-objective problem and find the optimal value of the model through multiple solutions; the convergence rate is slow and difficult to guarantee the optimality of the solution. Considering the shortcomings of traditional methods, we use the Beetle Antennae Search (BAS) to optimize the model.
Beetle Antennae Search (BAS) is an intelligent algorithm proposed in 2017 [24], which simulates the foraging behavior of beetles in nature. Beetles find food by following odor information received through left and right antennae. By adjusting its direction according to odor intensity difference between antennas, the beetle finally finds food after several iterations.
In order to increase the solving efficiency of BAS, some simplifications should be applied as follows: 1.
The beetle is simplified as a particle, and the left and right antennae are located on both sides of the particle.

2.
The ratio of the movement step size L step to the distance between the antennas d0 is constant, that is L step = c × d0.

3.
The head orientation of the beetle is random after one step.
After completing the above simplification, the optimization process framework of BAS for IDP-T includes four steps:

1.
Beetle position initialized. Select a group of solution in the optimization model of IDP-T solution space as the initial position of beetle x.

2.
Set left and right antennae and distance between two antennae. The fitness value F min is obtained by the left and right antennae of beetles, represented by xl and xr respectively, and the difference value between the two antennae is represented by d0. Since x, xl, xr are located in the same straight line, d0 = norm(xl − xr).
According to the simplified condition 3, the beetle faces in any direction, which means the direction vector that xl points to xr is arbitrary; thus, use the random vector dir = (x 1 , . . . , x i , . . . , x n ) to represent this direction vector, and normalize by dir = dir/norm(dir). Therefore, xl − xr = d0 × dir, when using the particle in simplifications condition 1, xl and xr can be expressed as xl = x + d0 × dir/2 and xr = x − d0 × dir/2, respectively. 3.
The orientation and step size of a beetle. The corresponding fitness values of left and right antennae are represented by f le f t = f (xl) and f right = f (xr) respectively. Then, the beetle moves with a step size of L step . The forward direction of the beetle is shown as Equation (20).
4. Output the optimal scheme and fitness. After comparing the fitness of the left and right tentacles for one step, the step size is updated in L step = L step × eta, where eta is the step size update factor [25]. When the precision requirement p or the specified number of steps n is reached, the beetle stops moving. When thee position of the beetle is the optimal solution of the model, then output position x with best scheme and corresponding fitness value F(x).
The IDP-T optimization framework based on BAS is shown in Figure 4. Then, the beetle moves with a step size of step L . The forward direction of the beetle is shown as Equation (20).
4. Output the optimal scheme and fitness. After comparing the fitness of the left and right tentacles for one step, the step size is updated in step is the step size update factor [25]. When the precision requirement p or the specified number of steps n is reached, the beetle stops moving. When thee position of the beetle is the optimal solution of the model, then output position x with best scheme and corresponding fitness value ( ) F x .
The IDP-T optimization framework based on BAS is shown in Figure 4.

Case Study
As a widely used mechanical product, a gearbox has great remanufacturing value and potential. Taking the remanufacturing of a used gearbox as an example to verify the model, the structure and composition number of the used gearbox is shown in Figure 5, and the dimension chain is shown in Figure 6.

Case Study
As a widely used mechanical product, a gearbox has great remanufacturing v and potential. Taking the remanufacturing of a used gearbox as an example to verify model, the structure and composition number of the used gearbox is shown in Figu and the dimension chain is shown in Figure 6.  According to Figure 6, the closed-loop is formed by the surface 9 A , 10 A , and precision range is 0.2 mm-0.4 mm. After the remanufacturing reuse decision of each the reuse mode of the right bearing bush is "replaced by new", and the size and toler of the new component is +0.050 0.050 6 − mm, the left bearing bush is reused, the left box, r box and gear shaft need remanufacturing, and the failure surfaces which need integr design are 1 3 4 5 8 9 , , , , ,

Remanufacturing Process Plan Formulation of Failure Surface
According to Section 2.2.1, FTA is used to identify the failure type of the fa surface, as shown in Figure 7, in which a represents the failure performance of used box, b the failure type and cause of the used parts, and c the extracted failure type o failure surface.

Case Study
As a widely used mechanical product, a gearbox has great remanufacturing value and potential. Taking the remanufacturing of a used gearbox as an example to verify the model, the structure and composition number of the used gearbox is shown in Figure 5, and the dimension chain is shown in Figure 6.  According to Figure 6, the closed-loop is formed by the surface 9 A , 10 A , and the precision range is 0.2 mm-0.4 mm. After the remanufacturing reuse decision of each part, the reuse mode of the right bearing bush is "replaced by new", and the size and tolerance of the new component is +0.050 0.050 6 − mm, the left bearing bush is reused, the left box, right box and gear shaft need remanufacturing, and the failure surfaces which need integrated design are 1 3 4 5 8 9 , , , , ,

Remanufacturing Process Plan Formulation of Failure Surface
According to Section 2.2.1, FTA is used to identify the failure type of the failure surface, as shown in Figure 7, in which a represents the failure performance of used gear box, b the failure type and cause of the used parts, and c the extracted failure type of the failure surface. According to Figure 6, the closed-loop is formed by the surface A 9 , A 10 , and the precision range is 0.2-0.4 mm. After the remanufacturing reuse decision of each part, the reuse mode of the right bearing bush is "replaced by new", and the size and tolerance of the new component is 6 +0.050 −0.050 mm, the left bearing bush is reused, the left box, right box and gear shaft need remanufacturing, and the failure surfaces which need integrated design are A 1 , A 3 , A 4 , A 5 , A 8 , A 9 .

Remanufacturing Process Plan Formulation of Failure Surface
According to Section 2.2.1, FTA is used to identify the failure type of the failure surface, as shown in Figure 7, in which a represents the failure performance of used gear box, b the failure type and cause of the used parts, and c the extracted failure type of the failure surface.
The volume damage of the failure surface is obtained with a CAD model. Then, the failure quantification is obtained according to the corresponding criteria in Table 1. Once the failure type and quantification of failure surface is obtained, the failure characteristics set, and attribute feature set of each failure surface, is shown in Table 2.
Taking the surface A 1 as an example to illustrate the formulation process of feasible process plans, through AHP the weight coefficients of each indicator are obtained as {material, shape, size, failure type, failure quantification} with ω = (0.12, 0.25, 0.15, 0.25, 0.23). The similarity threshold is set to 80%. After matching, three cases are obtained, as shown in Table 3.  The volume damage of the failure surface is obtained with a CAD model. Then, the failure quantification is obtained according to the corresponding criteria in Table 1. Once the failure type and quantification of failure surface is obtained, the failure characteristics set, and attribute feature set of each failure surface, is shown in Table 2. Taking the surface 1 A as an example to illustrate the formulation process of feasible process plans, through AHP the weight coefficients of each indicator are obtained as {material, shape, size, failure type, failure quantification} with ω = (0.12, 0.25, 0.15, 0.25, 0.23). The similarity threshold is set to 80%. After matching, three cases are obtained, as shown in Table 3.   After expert revision, A 1 is difficult to be grinded. The feasible process plans of A 1 is shown as Equation (21). Feasible processes of other failure surfaces of this used gearbox are shown as Table 4.

Multi-Objective Function Construction
According to the companies experiment, the cost and quality loss parameters of each feasible process plan for failure surfaces is shown in Table 5 (some of the cost parameters is dummy due to the sensitivities), the machining capability and kind of each process is shown in Table 6; the tolerance transfer coefficient is ξ i = ξ new = 1, and the weight of the cost and quality loss is ω 1 = ω 2 = 0.5. Therefore, the objective function of IDP-T is established according to Equations (8) and (14), and the constraint function is established according to Equations (15) and (16).

Optimization by BAS
According to Table 4, there are 16 (2 × 1 × 1 × 2 × 2 × 2) remanufacturing process combinations of this used gearbox that can be obtained. The process is coded as Pabc, in which a represents the serial number of failure surface, b represents the serial number of feasible process plan (shown in Table 4), c represents the serial number of process. The parameters of the BAS algorithm are set, as shown in Table 7. Finally, the optimal IDP-T remanufacturing scheme for the used gearbox is shown in Table 8.

Comparison with Different Scheme
The original remanufacturing scheme of process and tolerance of this used gearbox is "highest matching process + original precision recovery". With this scheme, the cost and quality loss is shown in Table 9. When optimal, the tolerance of the original scheme in the proposed method is achieved without changing process scheme. The cost and quality loss is shown in Table 10.   By comparison with Table 8, the scheme optimized in this paper has significantly reduced the cost and quality loss. The results indicate that the IDP-T optimization model is more effective in optimizing the process. Tolerance for remanufacturing is integrated.

Summary and Conclusions
An IDP-T method based on fault features is presented for remanufacturing, including failure description set establishment, feasible process plan generation, and processtolerance scheme optimization. FTA and FCE are used to extract and quantify the failure types of the failure surfaces, and combines them with the attributes of the failure surfaces to construct the failure description set. Then, CBR is applied to generate feasible process plans for each of the failure surfaces; finally, IDP-T optimization model is established, and BAS is employed to select the optimal scheme. A used gear box is taken as an example, which proves the validity of the proposed methods.
In general, it may be concluded that the multi-objective optimization model with IDP-T, based on fault features for remanufacturing products, is more feasible and more effective than the traditional serial design in minimizing cost and quality loss. Without failure feature description and its integration with IDP-T, as well as optimization of this process, an optimization plan of process and tolerance for remanufacturing products may not sufficiently minimize cost and quality loss. To be specific, the same types of used product should be restored by different process-tolerance schemes due to the different failure features and damage degrees etc., by which the cost and quality loss of remanufacturing for used product will be decreased greatly. Moreover, other factors can influence the optimization plan of process-tolerance too, such as remaining life and environmental factors. Future research of IDP-T should take these factors into consideration.