An effective similarity determination model for case-based reasoning in support of low-carbon product design

Abstract

In the initial stage of low-carbon product design, design information is always uncertain and incomplete, as well as the coupling properties between design attributes, thus it requires retrospective coordination for design conflicts resulting from the inclusion of low-carbon requirements. Reusing the prior design knowledge can promote design efficiency, however, the acquisition of similar cases knowledge not only needs to consider the similarity of design problems, but also the adaptability of candidate cases. This study presents an effective similarity determination model to support low-carbon product design, and targets of the proposed model are (1) to reasonably determine design ranges of attribute values for product cases retrieval by representing the uncertain design attributes with fuzzy set theory; (2) to construct an efficient indexing structure to generate the index set of similar cases based on the improved discretized highest similarity method by proposing two effective strategies; and, (3) to establish similarity estimation models for different types of attributes, and it calculates the information content of each attribute to evaluate the adaptability of cases based on the Information Axiom. The applicability of the proposed model is demonstrated through a case study of similar cases retrieval for the vacuum pump low-carbon design.

Introduction

In the 21^st century, human society confronts with significant challenges from climate change. The IPCC’s fourth report of environmental assessment pointed out that climate change and global warming are mainly caused by the large amount of greenhouse gases generated by human activities, where the industrial production and manufacturing account for a large proportion. ¹ The International Academy for Production Engineering (CIRP) established the low-carbon manufacturing working group (CWG) in Pairs, France, in January 2008, it aims to guide enterprises to take measures to reduce carbon emissions throughout product life cycle, including energy saving and emission reduction at manufacturing stage, green supply chain, services at use stage, and components recycling at the end of life stage. ² The preliminary design stage of the product influences about 80% economic cost, environmental and social impacts in its whole life cycle. ³ Low-carbon product design takes account of the environmental indicators during product planning and conceptual design stages, and puts forward the strategies and ideas for the generation of new design schemes with the purpose of carbon emission reduction, and it has attracted a heated research in both academia and industry.^4–6

The research of low-carbon product design mainly focuses on the carbon footprint estimation,^7–9 low-carbon design theory and methodology,^10,11 and low-carbon optimization.^12–14 Carbon footprint, the sum of greenhouse gases emissions at each life cycle stage, is a wildly accepted quantitative metric to evaluate the environmental impact of the product. Song and Lee ¹⁵ estimated the life cycle carbon footprint of parts based on the design information obtained from the bill of material (BOM), and established the greenhouse gases-BOM (g-BOM) to identify parts with high carbon emissions, which indicates the direction for subsequent improvement and optimization of the parts. Devanathan et al. ¹⁶ combined the quality function deployment (QFD) and life cycle assessment (LCA) method to construct a function-environment impact matrix, which enables designers to evaluate the environmental impact of the product function at design stage. However, in the low-carbon product design, since the coupling properties of design quality characteristics, it will inevitably break the original design equilibrium system when low-carbon requirements are incorporated, which results in the complex design conflicting problems. For instance, modifying the component to reduce its carbon emission, but the economic cost is increased and the mechanical property is deteriorated. Therefore, it is necessary to employ the conflicting problem-solving method to coordinate the design contradictions and generate design strategies to make a trade-off between carbon footprint, cost and function performance. The common used contradictory problems solving methods include the theory of inventive problems solving (TRIZ, a Russian acronym),^17,18 Extenics,^19–21 constraint satisfaction problems (CSP), ²² the negotiation method.^23,24

From the abovementioned studies in low-carbon product design, it indicates that whether it is the carbon footprint estimation, the generation of low-carbon design schemes through solving the contradictory problems, and the low-carbon optimization method, they are all dependent on the preceding design knowledge. That is, it is a reasoning innovation on the basis of the product case knowledge. Thus, integrating the case-based reasoning (CBR) technique and design conflicting problem-solving methods can promote the design efficiency, ²⁵ and to accurately grasp the similar cases knowledge corresponding to the low-carbon requirements is the prerequisite for low-carbon design. In our research work, it presents an effective similarity determination model to obtain similar product cases in consideration of the similarity and adaptability for knowledge reuse and adaptation for the low-carbon product design.

The CBR method reuses the prior case knowledge to solve the similar design problems, it stores the product cases with the solution to problems previously solved in the case base, rather than extracting the complex rule knowledge, such as the rule-based method.^26,27 Whereas CBR method requires less knowledge extraction in the knowledge acquisition task, it needs to collect a significant number of product cases to wildly cover the design problem and solution spaces in order to handle the new emergent design problems, which contributes to the retrieval utility problem of the case base. Retrieval utility problem refers to the opposite property between the accuracy and the efficiency, ²⁸ that is, when the number of product cases reaches the saturation point, the retrieval accuracy is basically stable, and if it continually adds new cases, the retrieval efficiency decreases. On the other hand, it is unreasonable to select the candidate cases by the single metric, similarity. Because the most similar candidate case may be not the most appropriate for the adaptation purpose. Consequently, it requires additional knowledge to construct the similarity determination model for the candidate cases in consideration of the similarity and adaptability.^29,30

The most commonly used cases retrieval method is the nearest neighbor algorithm and its various improved variants,^31,32 nearest neighbor method is wildly adopted for the flat case base as its simplicity and practicability. However, it needs to calculate the distance between each case in the case base (memory) and the query case, and then makes the final decision, thus, it will consume extensive searching time when facing the large size case base. Compared with the flat case base, the structure tree indexing method is more suitable for the cases retrieval in a hierarchical structure case base; it hierarchically divides the case base into multiple subsets based on the attribute nodes, thus it only needs to search the subset, in which the attribute values of the cases falling within the value range of the node attributes. Hjaltason and Samet ³³ summarized various structure tree indexing methods, and Stephane et al. ³⁰ pointed out the complexity in case base maintenance caused by the structural expansion under dynamic cases addition, and revealed the limitation in processing the overlapping cases. In addition, some useful machine learning methods are applied to the case retrieval, such as the neural network classification method, ³⁴ genetic algorithm, ³⁵ the clustering method.^36,37

Aiming at addressing the retrieval utility problem of the case base, the following three measures are generally adopted: ²⁸ (1) construct the indexing structure for the case base, in this approach, the number of cases is not limited, and the constructed indexing structure can improve the retrieval efficiency without reducing the accuracy; however, it is not easy to construct an appropriate indexing structure according to the attribute properties of the cases; (2) formulate the principle of cases addition, that is, control the number of cases in the case base by assessing whether the new case satisfies the addition principle; (3) formulate the principle of cases deletion, namely, control the number of cases by evaluating whether the existing cases should be removed from the case base based on the deletion principle. The difficulty of the latter two measures lies in the formulation of evaluation criteria for the principles. In terms of the adaptability of the similar cases, Smyth and Keane ²⁹ pointed out that it requires additional knowledge to select the candidate case in addition to the similarity, and they proposed the concept of adaptation-guided retrieval scheme. Leake et al. ³⁸ further studied the relationship between the similarity and adaptability. Rousu and Aarts ³⁹ incorporated the adaptation cost into the evaluation indicators of the similarity and adaptability. Roldan Reyes et al. ⁴⁰ reported that the CSP and CBR depend on the general knowledge and contextual knowledge respectively, and they integrated these two methods to address adaptation problems in the process of engineering design. Qi et al. ⁴¹ incorporated the support vector machine (SVM) and adaptability-related knowledge to obtain adaptable cases, the knowledge of adaptability contains the adaptability characteristics of the old cases and the guideline to the construction of the adaptability-involving SVM adaptation model with respect to retrieved cases.

In addition, in the earlier low-carbon design stage for products, facing the fuzzy low-carbon requirements, the incomplete design information, and the potential complex design conflicting problems, even if the designers are experienced domain experts, they cannot grasp the deep knowledge to determine the precise values for the design attributes. But they roughly know the possible design ranges, such as the lower and upper boundaries of the attribute values to guarantee that the decision making for the design parameters is still satisfying the design targets. Giachetti et al. ⁴² revealed this uncertainty in the expression of the design information at each design stage, from the linguistic variables, fuzzy numbers, to the final real numbers. It also exactly demonstrates that the essence of product design is a complex backtracking process of clarifying the “ill-defined” target problems. ⁴³ The rest of this article is organized as follows: in the next section, low-carbon requirements are detailed analyzed, and the fuzzy set theory is used to represent the uncertain design attributes; the discretized highest similarity indexing method is improved with two proposed strategies to construct the effective indexing structure for the low-carbon product case base; the comprehensive estimation model by integrating the similarity and adaptability is presented to evaluate the candidate cases; a case study about the similar cases retrieval for the vacuum pump low-carbon design is provided to demonstrate the effectiveness of our proposed method; discussion is given, and conclusion is finally drawn along with the recommendations for the further research.

Requirement analysis and representation for low-carbon product design

For low-carbon product design, it is essential to obtain similar cases from the case base according to the low-carbon requirements, and reuse the prior knowledge and conduct the adaptation to generate a new design scheme. Thus, it needs to extract the retrieval attributes from the requirement information. In addition, as the design information is incomplete and unascertained at earlier design stage, in this paper, different kinds of descriptions of the requirement information are converted into the standardized fuzzy numbers.

Low-carbon requirement analysis

Low-carbon requirement comprises the traditional customer preference, that is the appearance attributes of the product, the general quality performance, etc., which is named the voice of customers (VOC), and the low-carbon demand, such as less material usage, low energy consumption, less carbon emissions. Mausi et al. ⁴⁴ summarized the customer requirements for the environment, defined as voice of environment (VOE), and the quality characteristics (QCs) for the product environmental conscious design (Eco-design), as listed in Table 1. Low-carbon design is a branch of the Eco-design, while the difference in VOC with different products, the VOE and corresponding QCs involved in the design requirements for different products are basically similar. Thus, this study incorporates those low-carbon related VOE and QCs into the design requirement analysis.

Table 1.

The customer requirements and quality characteristics for the environment. ⁴⁴

Requirements for the environment	Quality characteristics
Less material usage	Mass
Less energy consumption	Volume
Use of environmental friendly materials	Number of parts
Use of environmental friendly technology	Number of material types
Easy to transport and retain	Mechanical properties
Easy to process and assemble	Physical lifetime
High durability	Amount of energy consumption
Easy to reuse	Rate of recycled materials
Easy to disassemble	Toxicity of materials
Easy to sort	Degree of air pollutant
Harmless to the living environment	Degree of water pollutant
Safe to incinerate	Degree of soil pollutant
Safe to landfill	Biodegradability
Safe to emission	Noise, electromagnetic wave

In addition, QFD method is a useful design tool to collect vaguely described needs and deploy them to the actual design work. The relationship of low-carbon requirements and quality characteristics can be established based on the house of quality (HoQ), as shown in Figure 1. According to the deployment for the low-carbon requirements, designers can identify the important quality characteristics with high relative weight, and extract them as the retrieval attributes to search the similar cases knowledge for low-carbon design.

Figure 1.

Mapping relationship of low-carbon requirements and QCs based on the HoQ.

Representation of the uncertain requirement information with fuzzy set theory

In the life cycle-oriented low-carbon design process, the requirement information not only comes from the product users, but also includes the feedback information from the workers who are in charge of the product manufacturing, assembly, sales, maintenance, and recycling. However, it should be noted that some of this requirement information can be quantitatively expressed, for instance, specify the numerical attribute value with a crisp number, or designate it in a required range; and some requirement information is quantitative immeasurable, it should be expressed with a linguistic term. For example, it adopts the linguistic term to describe the low-carbon and environmental impact (LCEI) of the product as “LCEI is poor”, “LCEI is fair”, “LCEI is good”.

Thus, the different types of descriptions of the requirement information, including the precise numerical value, a range of the numerical value, and the linguistic term, make it difficult to construct the evaluation criteria of the similarity and adaptability for the cases. Fortunately, the above requirement information can be quantitative measured based on the fuzzy set theory.^45,46 The fuzzy set in a domain U is defined by a membership function μ, which has the values in range [0, 1]. μ(x) indicates the degree to which x is the possible value in a fuzzy subset. In this paper, the uncertain attribute values of the requirement information are transformed into the standardized triangular fuzzy number, named STFN (Triangle), or the standardized trapezoidal fuzzy number, named STFN (Trapezoid), as illustrated in Figure 2.

Figure 2.

Membership functions for the standardized fuzzy number conversion: (a) membership function of STFN (Triangle) and (b) membership function of STFN (Trapezoid).

Here, it takes the STFN (Trapezoid) as an example to describe its transformation effect. Suppose the attribute value A, the transformation result of A based on the fuzzy set theory is represented as STFN (A) = (a ^l , a ^m , a ⁿ , a ^u ), a ^l ≤a ^m ≤a ⁿ ≤a ^u ≤u, as shown in Figure 2(b), where the domain U of the fuzzy set is U = [0, u], and the membership function is as equation (1).

$$\mu (x)=\{\begin{array}{cc}(x-a^l)/(a^m-a^l)\hfill & a^l\le x\le a^m\hfill \\ \multicolumn{1}{c}{1}& a^m\le x\le a^n\hfill \\ \multicolumn{1}{c}{(a^u-x)/(a^u-a^n)}& a^n\le x\le a^u\hfill \\ \multicolumn{1}{c}{0}& \mathrm{for}\mathrm{otherwise}\hfill \end{array}$$

(1)

When a ^l = a ^m = a ⁿ = a ^u , A is a numerical value and it is transformed into a STFN(Trapezoid); when a ^l = a ^m , a ⁿ = a ^u , A is a range value; and when a ^m = a ⁿ , the STFN (Trapezoid) is equal to the STFN (Triangle).

Set U = [0, 1], the linguistic term can be quantitatively expressed by STFN (Triangle) or STFN (Trapezoid) as described in Figure 2(a) and Table 2, and the fuzzy numerical value can also be converted into the STFN (Triangle) or STFN (Trapezoid) as expressed in Table 3.

Table 2.

Descriptions of the linguistic term (take the LCEI as the example).

The linguistic term	Interpretation	STFN (Triangle)	STFN (Trapezoid)
Very poor (VP)	LCEI of the candidate case is very poor	0.0, 0.0, 0.25	0.0, 0.0, 0.0, 0.25
Poor (P)	LCEI of the candidate case is poor	0.0, 0.25, 0.50	0.0, 0.25, 0.25, 0.50
Fair (F)	LCEI of the candidate case is fair	0.25, 0.50, 0.75	0.25, 0.50, 0.50, 0.75
Good (G)	LCEI of the candidate case is good	0.50, 0.75, 1.0	0.50, 0.75, 0.75, 1.0
Very good (VG)	LCEI of the candidate case is very good	0.75, 1.0, 1.0	0.75, 1.0, 1.0, 1.0

Table 3.

Descriptions of the fuzzy numerical values.

Fuzzy numerical values	Interpretation	STFN (Triangle)	STFN (Trapezoid)
About 0.7	the attribute value is about 0.7	0.6, 0.7, 0.8	0.6, 0.7, 0.7, 0.8
Range (0.4, 0.8)	the attribute value is in range (0.4, 0.8)	0.4, 0.6, 0.8	0.4, 0.6, 0.6, 0.8
At least 0.8	the attribute value is at least 0.8	0.8, 1.0, 1.0	0.8, 1.0, 1.0, 1.0
At most 0.4	the attribute value is at most 0.4	0.0, 0.0, 0.4	0.0, 0.0, 0.0, 0.4

The improved D-HS indexing method for the index set of similar cases

Indexing method is used to gather a group of similar cases for similarity calculation, and thus it can avoid the extensive searching for the whole memory. In this paper, an efficient indexing method named Discretized Highest Similarity (D-HS) ⁴⁷ is introduced, however, this general approach presents its weakness in application, two effective strategies are proposed to construct the improved indexing structure of the case base for low-carbon design.

Principle of the D-HS indexing method

With the principle of D-HS method, it distributes the intervals of each attribute value of the product cases, and designates the spatial coordinates for the cases in the memory. Cases with the same coordinates as that of the query case constitute the similar group. The simple indexing structure of D-HS method facilitates the maintenance for the case base, and it can effectively address the utility problem.

Let C represent the low-carbon case base, it is composed of N product cases (c _i ), and each case consists of M attributes (c_j) in a vector, as expressed in equations (2) and (3).

$$\mathrm{C}=\{{\mathrm{c}}_i,i=1,2,\dots ,N\}$$

(2)

$$\mathrm{c}=[c_j],j=1,2,\dots ,M$$

(3)

where, it defines the expression c_j(c _i ) as the value of attribute c_j of the c _i .

Suppose the low-carbon query case q, that is the target case, which consists of M attributes q_i, and thus the query case can be described as below.

$$\mathrm{q}=\left[q_i\right],i=1,2,\dots ,M$$

(4)

where, it defines the expression q_i(q) as the value of the attribute q_i of the query case q.

Let S represent the index set of similar cases, generated by the indexing rule R based on the D-HS principle as expressed in equation (5).

$$\begin{array}{c}S=R(\mathrm{q})=\underset{i=1}{\overset{M}{\cap }}{\mathrm{C}}^i,\hfill \\ {\mathrm{C}}^i=\{{\mathrm{c}}_j\in \mathrm{C}:c_i({\mathrm{c}}_j)\in [q_i(\mathrm{q})-\frac{\mathrm{\Delta }{\delta }_i}{2},q_i(\mathrm{q})+\frac{\mathrm{\Delta }{\delta }_i}{2}]\}\hfill \end{array}$$

(5)

Where, △δ_i denotes the interval width of the i-th attribute value; C ⁱ is the subset of cases, in which the attribute value of each case falls within the same interval as that of the query case for the i-th attribute.

Based on the principle, the indexing rule R can be described in detail as equation (6).

$$\{\begin{array}{c}\multicolumn{1}{c}{S=R(\mathrm{q})=\underset{i=1}{\overset{M}{\cap }}{\mathrm{C}}_{l_i}^i}\\ \multicolumn{1}{c}{{\mathrm{C}}_{l_i}^i=\{{\mathrm{c}}_j\in \mathrm{C}:c_i({\mathrm{c}}_j),q_i(\mathrm{q})\in [(l_i-1)\mathrm{\Delta }{\delta }_i,l_i\mathrm{\Delta }{\delta }_i]\}}\\ \multicolumn{1}{c}{\mathrm{\Delta }{\delta }_i=\frac{1}{L_i},l_i=1,...,L_i,j=1,2,...,N}\end{array}$$

(6)

Where, L_i denotes the number of intervals of the i-th attribute, l_i is the index number of L_i.

In equation (6), each attribute is discretized by designating the intervals. The number of intervals of the nominal attributes equals to the number of discrete attribute values, while the number of intervals of the numerical attributes is predefined by the designer. Take the two dimensional indexing space as an example, M = 2, c = [c₁, c₂] = [c_x, c_y], define L_x = L_y = 4, then the discrete intervals of c_x are described as: [0, 0.25], (0.25, 0.5], (0.5, 0.75], (0.75, 1]; the corresponding subsets of cases are $\{{\mathrm{C}}_1^x,{\mathrm{C}}_2^x,{\mathrm{C}}_3^x,{\mathrm{C}}_4^x\}$ . The discrete intervals of c_y are described as: [0, 0.25], (0.25, 0.5], (0.5, 0.75], (0.75, 1]; and the corresponding subsets of cases are $\{{\mathrm{C}}_1^y,{\mathrm{C}}_2^y,{\mathrm{C}}_3^y,{\mathrm{C}}_4^y\}$ . The two dimensional indexing structure is given as in Figure 3. Suppose the query case falls in the interval with index 3 of the attribute c_x, and the interval with the index 3 of the attribute c_y, then the cases falling within the intersection space constitute the index set S, that is $S={\mathrm{C}}_3^x\cap {\mathrm{C}}_3^y$ = {c₄, c₇, c₁₂,c₁₅}.

Figure 3.

The indexing structure of D-HS method adapted from the literature 47: (a) the distribution of cases in two dimensional space and (b) the scaled up view of the index set S.

D-HS method makes the similarity evaluation limited to the cases in the index set S, thus it greatly promotes the retrieval efficiency. Besides, it is easy to keep the maintenance for the case base when new cases are added or removed. However, the indexing structure of D-HS in Figure 3 has the limitations in retrieval accuracy, which mainly come from two aspects: (1) while the interval width of each attribute is equally designated, the attribute values of cases cannot be evenly distributed. (2) when the location of the query case is close to the boundaries of intervals, the similar cases are easily neglected.

Strategies to improve the D-HS method

Strategy that keep approximately equal number of cases for attribute intervals division

From the equation (6), it can be seen that the widths of intervals of one certain attribute are equally divided, that is, each interval width of the i-th attribute is 1/L_i. However, it cannot make sure the even distribution of the i-th attribute values of cases in each interval, which significantly impacts the retrieval result. For instance, the most cases fall within one certain interval, while few cases are sparsely distributed in other intervals. In this context, the index set S may gathers most of the cases, and thus it cannot achieve the effect of reducing the similarity calculation for searching the whole memory by the indexing approach. On the contrary, S may contain no one case.

Fayyad and Irani ⁴⁸ proposed a multi-interval discretized method for continuous numerical attributes based on the information entropy, and constructed the MDLP interval division criteria. Galushka and Patterson ⁴⁷ applied the above method to construct a new indexing structure, name D-HS^E, and proposed a selection mechanism between the D-HS^M and the D-HS^E according to the distribution of the cases, the former has better retrieval efficiency, and the latter has higher accuracy. Based on the principle of D-HS method, this study proposes the strategy that keep the number of cases approximately equal in each attribute interval to make cases evenly distributed in a new indexing structure, the procedure is as below.

For the i-th attribute, cases are sorted by the attribute value from smallest to largest, and it defines the number of cases in each interval of the i-th attribute as Int(N/L_i), where, Int() is the rounding function. Thus, the lower boundary of each interval is the minimum attribute value of cases falling within the interval, and the upper boundary of the interval is the maximum attribute value. The indexing rule is improved as equation (7).

$$\{\begin{array}{c}S=R(\mathrm{q})=\underset{i=1}{\overset{M}{\cap }}{\mathrm{C}}_{l_i}^i\hfill \\ {\mathrm{C}}_{l_i}^i=\{\mathrm{c}\in \mathrm{C}:c_i({\mathrm{c}}_j),q_i(\mathrm{q})\in [[c_i({\mathrm{c}}_{((l_i-1)\mathrm{\Delta }{\delta }_i)}),c_i({\mathrm{c}}_{(l_i\mathrm{\Delta }{\delta }_i)})]\cup [c_i({\mathrm{c}}_{((L_i-1)\mathrm{\Delta }{\delta }_i)}),1]]\}\hfill \\ \mathrm{\Delta }{\delta }_i=\mathrm{Int}(N/L_i)\hfill \\ c_i({\mathrm{c}}_0)=0,l_i=1,...,(L_i-1),j=1,2,...,N\hfill \end{array}$$

(7)

Where, N is the number of cases. L_i is the predefined number of intervals of the i-th attribute; $c_i({\mathrm{c}}_{(l_i\mathrm{\Delta }{\delta }_i)})$ denotes the attribute value of the (l_i△δ_i)-th case in the sequence.

However, the number of cases in each attribute interval is usually not equal as the overlapping attribute values of cases, and thus, it should adjust the number of cases in each attribute interval, and even redefine the number of attribute intervals. For this problem, the criterion to identify whether the attribute values of the cases are evenly distributed in one certain attribute intervals is constructed, as expressed in equation (8).

$$A_i^{\mathrm{Well}}=\{\begin{array}{cc}1,\hfill & \frac{{\sigma }_i}{{\sigma }_{max(i)}}<\sigma *\hfill \\ \multicolumn{1}{c}{0,}& \frac{{\sigma }_i}{{\sigma }_{max(i)}}\ge \sigma *\hfill \end{array}$$

(8)

In equation (8), when $A_i^{\mathrm{Well}}$ = 1, it indicates the values of the i-th attribute of cases are basically evenly distributed; and when $A_i^{\mathrm{Well}}$ = 0, it reveals that the interval division for the i-th attribute is inappropriate, and thus it needs to adjust the number of intervals L_i, and the number of cases in each interval. σ_i denotes the standard deviation of counts of values of the i-th attribute falling within the predefined L_i intervals, it is calculated as equation (9); σ_max(i) denotes the standard deviation in the scenario where all the values of the i-th attribute falls into a single interval, based on equation (9), ${\sigma }_{max(i)}=N/\sqrt{L_i}$ ; σ* is a threshold to identify the acceptable balance between the well and poor distribution of the attribute values.

$${\sigma }_i=\sqrt{\frac{L_i\sum _{l_i=1}^{L_i}{n}_{l_i}^2-{\left(\sum _{l_i=1}^{L_i}{n}_{l_i}\right)}^2}{L_i(L_i-1)}}$$

(9)

Where, n_li denotes the number of cases falling within the interval with index l_i of the i-th attribute.

Strategy of similarity matching for attributes

When the query case q does not fall within the center of the intervals, as shown in Figure 3, the location of q is near to the boundaries of the intersection. In this scenario, cases c₁, c₃, c₅, and c₁₉ are not included in the index set S, however they are more closer to the query case q than some cases in S, which results in the retrieval inaccuracy.

Brodu ⁴⁹ proposed a spherical indexing method for neighboring cases, such as the two dimensional attributes in Figure 3(b), it takes the location of query case q as the center to draw a circle, then cases in the circle constitute the index set S. However, this query sphere algorithm is not applicable to the hyper polyhedron in case that the lower and upper boundary values of one or more attribute intervals are asymmetric with respect to the recommended points. In this paper, a simple and effective strategy of similarity matching for attributes is proposed to obtain the similar edge candidate cases.

Suppose there are 10 cases in the case base with four dimensional attributes, and each attribute is designated four intervals. Set the attribute values of candidate cases and the query case: c₁(0.30, 0.55, 0.40, 0.15), c₂(0.30, 0.85, 0.20, 0.65), c₃(0.15, 0.35, 0.30, 0.95), c₄(0.20, 0.60, 0.35, 0.90), c₅(0.80, 0.60, 0.65, 0.70), c₆(0.40, 0.45, 0.65, 0.80), c₇(0.80, 0.30, 0.35, 0.85), c₈(0.20, 0.90, 0.70, 0.65), c₉(0.80, 0.70, 0.20, 0.20), c₁₀(0.80, 0.30, 0.40, 0.80), q = (0.20, 0.30, 0.30,0.80), which can be arranged in Table 4.

Table 4.

Attribute values matching between the candidate cases and query case.

IntervalsAttributes	0.0–0.25	0.25–0.50	0.50–0.75	0.75–1.0
A 1	q, c3, c4, c8	c1, c2, c6		c5, c7, c9, c10
A 2		q, c3, c6, c7, c10	c1, c4, c5, c9	c2, c8
A 3	c2, c9	q, c1, c3, c4, c7, c10	c5, c6, c8
A 4	c1, c9		c2, c5, c8	q, c3, c4, c6, c7, c10

In Table 4, each attribute value of case c₃ is matched with that of the query case q; cases c₄, c₇, c₁₀ all have three attribute values falling within the same intervals as that of the query case q, respectively; for case c₆, two attribute values are matched, that is, the values of attributes A₂ and A₄; cases c₁ and c₈ only have single attribute value matched with that of the query case q. It concludes that cases having more matched attribute values may close to the query case q, such as the cases c₃, c₄, c₇, and c₁₀. Based on this principle, the strategy of similarity matching for attributes can be quantitatively modeled as equation (10).

$$\mathrm{SMA}({\mathrm{c}}_i,\mathrm{q})=\sum _{j=1}^M\mathrm{Match}(c_j({\mathrm{c}}_i),q_j(\mathrm{q}))$$

(10)

Where, SMA(c_i, q) indicates the number of matched attributes between case c _i and the query case q; Match() is a matching function, based on the strategy that keep approximately equal number of cases for attribute intervals division as expressed in equation (7), then the matching function is constructed as below.

If $c_j({\mathrm{c}}_i),q_j(\mathrm{q})\in [[c_j({\mathrm{c}}_{((l_j-1)\mathrm{\Delta }{\delta }_j)}),c_j({\mathrm{c}}_{(l_j\mathrm{\Delta }{\delta }_j)})]\cup [c_j({\mathrm{c}}_{((L_j-1)\mathrm{\Delta }{\delta }_j)}),1]]$ ,

then Match (c_j(c _i ), q_j(q)) = 1;

If $\{\begin{array}{c}[c_j({\mathrm{c}}_i)\in [c_j({\mathrm{c}}_{((l_{j1}-1)\mathrm{\Delta }{\delta }_j)}),c_j({\mathrm{c}}_{(l_{j1}\mathrm{\Delta }{\delta }_j)})]]\cap [c_j({\mathrm{c}}_i)\notin [c_j({\mathrm{c}}_{((L_j-1)\mathrm{\Delta }{\delta }_j)}),1]]\hfill \\ \multicolumn{1}{c}{[q_j(\mathrm{q})\in [c_j({\mathrm{c}}_{((l_{j2}-1)\mathrm{\Delta }{\delta }_j)}),c_j({\mathrm{c}}_{(l_{j2}\mathrm{\Delta }{\delta }_j)})]]\cap [q_j(\mathrm{q})\notin [c_j({\mathrm{c}}_{((L_j-1)\mathrm{\Delta }{\delta }_j)}),1]]}\end{array},$ and l_j1≠l_j2,

then Match (c_j(c _i ), q_j(q)) = 0.

Based on the strategy of similarity matching for attributes, a new index set of similar cases S_b is constructed as in equation (11), which is a supplement to the index set S obtained from equation (7), and it can effectively avoids the neglect of similar cases when the query case q locates near the boundaries of the intervals by setting the threshold Const of the number of matched attributes. In final, the index set of similar cases S₀ in equation (12) is obtained.

$$S_b=\{{\mathrm{c}}_i:\mathrm{SMA}({\mathrm{c}}_i,\mathrm{q})\ge \mathrm{Const},\mathrm{Const}\le M\}$$

(11)

$$S_0=\{{\mathrm{c}}_i:{\mathrm{c}}_i\in S\cup S_b\}$$

(12)

Construction of the indexing structure for the low-carbon product case base

Based on the improved D-HS method, the indexing structure of the product case base for low-carbon design is established as illustrated in Figure 4.

• Step 1: sequence the cases for M dimensional attributes, respectively; designate the number of intervals of each attribute; and set up the attributes coordinate system for the product case base, labeled as ((0, L_i), (0, L_j), (0, L_k),…).

• Step 2: designate the coordinates for each product case according to the values of M dimensional attributes, labeled as (l_i(c _i ), l_j(c _i ), l_k(c _i ),…).

• Step 3: Input the requirement information, namely the values of M dimensional attributes of the query case q, and designate the attribute coordinates for the query case, labeled as (l_i(q), l_j(q), l_k(q),…).

• Step 4: generate the index set S, in which the coordinates of cases are same as that of the query case based on the indexing rule S = R(q).

• Step 5: define the threshold Const and gather the cases which satisfy the condition SMA(c _i , q)≥Const, these cases constitute the index set S_b, and it obtains the final index set S₀ = S∪S_b.

• Step 6: evaluate the similarity and adaptability of each product case in S₀, and output the acceptable cases for knowledge reuse and adaptation. The evaluation models for the similarity and adaptability are described in the following sections.

• Step 7: if the output case is successfully adapted to satisfy the design requirements, then it is considered as a new product case added to the case base with updated attribute coordinates.

• Step 8: the total number of cases will increase when new adapted cases are continuously added into the case base, thus, if the current number of cases N₁ satisfies the condition N₁≥(1+10%) N, then the coordinate system of the case base should be updated.

Figure 4.

The procedure of the construction of the indexing structure for the case base.

Similarity and adaptability evaluation models for cases retrieval

Similarity evaluation

In this section, it establishes the similarity evaluation models for different kinds of attributes in consideration of the fuzzy requirement information for the cases in the index set S₀.

Similarity evaluation model for the precise numerical attribute

When the requirement information for one certain numerical attribute is ascertained, namely, the numerical attribute value of the query case q is a crisp number. For the precise numerical attribute, the classical distance measurement model is used to calculate the local similarity as expressed in equation (13).

$$\mathrm{si}{m}_1(c_j,q_j)=1-\frac{|c_j({\mathrm{c}}_i)-q_j(\mathrm{q})|}{D_j}$$

(13)

Where, D_j = maximum value –minimum value of the j-th attribute.

Similarity evaluation model for the fuzzy numerical attribute

When the requirement information for the numerical attribute is unascertained, that is, the numerical attribute value is a fuzzy number, and usually it is a range number. In this paper, it sets up the feasible interval and desired interval for each fuzzy numerical attribute to make the decision making of the parameters more flexible for the designer; and the dependent function presented in Extenics^19,50 is adopted to establish the similarity evaluation model for fuzzy numerical attributes, as expressed in equation (14).

$$\mathrm{si}{m}_2(c_k,q_k)=\{\begin{array}{cc}1,\hfill & K(c_k({\mathrm{c}}_i))\ge 1\hfill \\ \multicolumn{1}{c}{e^{K(c_k({\mathrm{c}}_i))-1},}& K(c_k({\mathrm{c}}_i))<1\hfill \end{array}$$

(14)

Where, K(c_k(c _i )) is the dependent function for the value of the k-th attribute of case c _i .

Suppose the feasible interval X = [a, b], and the desired interval X₀ = [c, d], satisfying X₀⊆X; given the ideal point x₀ of interval X₀, then the dependent function of an arbitrary point x is calculated as in equation (15), here, x can be taken as the attribute value of case c _i , namely x = c_k(c _i ).

$$K(x)=\frac{\rho (x,x_0,X)}{\rho (x,x_0,X)-\rho (x,x_0,X_0)}$$

(15)

Where, ρ(x, x₀, X) is named extension distance, which describes the distance of the point x to interval X; and ρ(x, x₀, X₀) is the distance of the point x to interval X₀, accordingly. However, the calculation of the extension distance depends on the location of the ideal point x₀. Here, it takes the feasible interval X as the example, when x₀ locates at the left side of the middle point of interval X, then it should calculate left side extension distance ρ_l(x, x₀, X), otherwise, it should calculate the right side extension distance ρ_r(x, x₀, X), as expressed in equation (16).

$$\begin{array}{c}{\rho }_l(x,x_0,X)=\{\begin{array}{cc}a-x,\hfill & x\le a\hfill \\ \multicolumn{1}{c}{\frac{b-x_0}{a-x_0}(x-a),}& x\in (a,x_0)\hfill \\ \multicolumn{1}{c}{x-b,}& x\ge x_0\hfill \end{array}\hfill \\ {\rho }_r(x,x_0,X)=\{\begin{array}{cc}a-x,\hfill & x\le x_0\hfill \\ \multicolumn{1}{c}{\frac{a-x_0}{b-x_0}(b-x),}& x\in (x_0,b)\hfill \\ \multicolumn{1}{c}{x-b,}& x\ge b\hfill \end{array}\hfill \end{array}$$

(16)

The dependent function of each attribute value of cases has the properties as below:

when x∈X₀, then K(x)≥1, it reveals that the attribute value completely meets the requirement;

when x∈X-X₀, then 0≤K(x)<1, it indicates that the attribute value is acceptable;

and when x∉X, then K(x)<0, it reveals that the attribute value cannot satisfy the requirement.

Similarity evaluation model for the nominal attribute

The l-th nominal attribute value of case c _i and query case q are firstly transformed into the standardized triangular fuzzy number STFN(c_l(c _i )) and STFN(q_l(q)), respectively; and then the local similarity of the l-th nominal attribute is estimated based on the ratio of area of the common range A_CR to the area of design range A_STFN(DR), as expressed below.

$$\mathrm{si}{m}_3(c_l,q_l)=\frac{A_{\mathrm{CR}}}{A_{\mathrm{STFN}(\mathrm{DR})}}$$

(17)

Where, for the l-th nominal attribute, the area of STFN(DR) is equal to the area of STFN(q_l(q)); A_CR is the overlapping section between the area of STFN(c_l(c _i )) and the area of STFN(DR).

Therefore, the global similarity of case c _i and query case q is calculated as in equation (18).

$$\begin{array}{c}\multicolumn{1}{c}{\mathrm{SIM}({\mathrm{c}}_i,\mathrm{q})=\sum _{j=1}^{m_1}{\omega }_j\mathrm{si}{m}_1(c_j,q_j)}\\ +\sum _{k=1}^{m_2}{\omega }_k\mathrm{si}{m}_2(c_k,q_k)+\sum _{l=1}^{m_3}{\omega }_l\mathrm{si}{m}_3(c_l,q_l)\hfill \end{array}$$

(18)

Where, m₁ is the number of precise numerical attributes, m₂ is the number of fuzzy numerical attributes, and m₃ is the number of nominal attributes; ω_j, ω_k, ω_l are the weight of the corresponding attributes, ∑ω_j+∑ω_k+∑ω_l = 1. In addition, each local similarity should be normalized in order to avoid the distorted result because of the different amplitudes caused by the different similarity evaluation models.

Adaptability evaluation

In this paper, it construct a novel adaptability evaluation model based on the fuzzy axiomatic design (FAD).^51,52 Axiomatic design (AD) ⁵³ is wildly used for the product design and development as its two foundational axioms, the Independence Axiom and the Information Axiom. In Independence Axiom, it requires that each functional requirement (FR) should be independent in order to construct the stable product system with an uncoupled design environment. The information in Information Axiom is referred to the probability of satisfying FRs, which provides a theoretical support for the robust design, where the design system can accommodate the large variations in design attributes values and yet satisfy the FRs. Based on the Information Axiom, in this section, it construct the evaluation model for the adaptability of candidate cases in S₀ by estimating the information content of each attribute, the smaller the information content of the case, the better the adaptability.

The concept of information content is defined to quantitatively estimate the information of the design system, suppose the information content I_k for the given FR _k , then I_k is expressed as below.

$$I_k={\log }_2(\frac{1}{P_k})$$

(19)

Where, P_k is the probability of achieving the FR _k ; and when there are more than one functional requirements, then the information content I of the design system is the sum of all probabilities. If one or more probabilities are zero, then I approaches to the infinity; in this scenario, the design system will not work although other design attributes are completely satisfied. Conversely, if all the probabilities are equal to one, then I approaches to zero, and this is an ideal design system.

In a design project, the probability of success for the design task is given by the tolerance specified by the designer, that is, the design range (DR), and the tolerance that the product system is capable of satisfying, that is, the system range (SR). Suppose the functional requirement FR _k follows the uniform distribution, and the system probability density function (pdf) is clear, as shown in Figure 5, it illustrates that the intersection between the design range and system range is the region where the acceptable solution exists. Therefore, P_k can be expressed as equation (20).

$$P_k=\frac{A_{\mathrm{CR}}}{A_{\mathrm{SR}}}=\frac{L_{\mathrm{CR}}}{L_{\mathrm{SR}}}$$

(20)

Where, A_CR, A_SR are the area of common range and system range respectively; and L_CR, L_SR are the interval width of common range and system range; accordingly, the information content for the FR _k is written as below.

$$I_k={\log }_2(\frac{L_{\mathrm{SR}}}{L_{\mathrm{CR}}})$$

(21)

When FR _k follows the random distribution, as shown in Figure 6, the probability of achieving the FR _k in the common range (CR) is expressed as equation (22).

$$P_k=\underset{\mathrm{C}{\mathrm{R}}^l}{\overset{\mathrm{C}{\mathrm{R}}^u}{\int }}\mathrm{Ps}(\mathrm{F}{\mathrm{R}}_k)d\mathrm{F}{\mathrm{R}}_k$$

(22)

Figure 5.

Identification for DR, SR, CR and A_CR with the uniform distribution.

Figure 6.

Identification for DR, SR, CR, and A_CR with the random distribution.

Where, CR ^l , CR ^u are the lower boundary and upper boundary of CR, respectively; Ps denotes the system probability density function (pdf). In this context, the A_SR is equal to one, thus, the information content I_k for FR _k is calculated as in equation (23).

$$I_k={\log }_2(\frac{1}{A_{\mathrm{CR}}})$$

(23)

However, in the earlier design stage, the information of design system is always incomplete and unascertained, thus it is hard to obtain the precise probability density function of the design system. In this paper, we adopt the membership function to transform the DR and SR to the STFN(DR) and STFN(SR) with crisp numbers as illustrated in Figure 7, and thus, the information content I_k is expressed as in equation (24).

$$I_k={\log }_2(\frac{A_{\mathrm{STFN}(\mathrm{SR})}}{A_{\mathrm{CR}}})$$

(24)

Where, A_STFN(SR) is the area of STFN(SR).

Figure 7.

The standardized triangular fuzzy number transformation for the DR, SR.

According to the evaluation model, it is crucial to obtain the DR, SR, and A_CR of each attribute for the information content calculation. Here, we also specify the rules for Benefit type and Cost type attributes to construct the design range.

For Benefit type attributes, that is, the attribute value is like “at least a”, or “more than a”, the design range can be constructed as in the Figure 8(a). where, a ^l = DR ^l , a ^m = a ^u = max(U(x)), μ(a ^l ) = 0, μ(a ^u ) = 1; DR ^l denotes the lower boundary of the design range, U(x) denotes the whole design domain.

Figure 8.

The STFN (triangle) transformation for DR of benefit and Cost type attributes: (a) STFN(DR)of the Benefit type attribute and (b) STFN (DR) of the Cost type attribute.

For Cost type attributes, that is, the attribute value is like “at most a”, or “less than a”, then the design range can be constructed as in the Figure 8(b). where, a ^l = a ^m = min(U(x)), a ^u = DR ^u , μ(a ^l ) = 1, μ(a ^u ) = 0; DR ^u denotes the upper boundary of the design range.

In addition, when the system range of the design attribute is a crisp value, that is, x_j = c_j(c _i ), in this scenario, the system range is constructed as in Figure 9, here, the Benefit type attribute is taken as the example; and the information content is calculated as in equation (25). ⁵²

$$\begin{array}{c}\multicolumn{1}{c}{I_j={\log }_2(\frac{1}{P(c_j(c_i))}),P(c_j({\mathrm{c}}_i))=u(x_j)}\\ =\{\begin{array}{cc}\frac{x_j-a^l}{a^u-a^l},\hfill & \mathrm{Benefit}\mathrm{type}\mathrm{attributes}\hfill \\ \multicolumn{1}{c}{\frac{a^u-x_j}{a^u-a^l},}& \mathrm{Cost}\mathrm{type}\mathrm{attributes}\hfill \end{array}\hfill \end{array}$$

(25)

Figure 9.

Construction of A_CR for the crisp SR (the Benefit type attribute).

A case study

Low-carbon requirement analysis for the vacuum pump

The CNC cutting machine is widely used to automatically cut the leather and fabric materials, however, unlike the rigid material processing, the materials should be fixed to the working table by the adsorption force with the pressure difference. In the design task, it needs to configure a suitable vacuum pump for the working requirement of the cutting machine in consideration of low-carbon impact for the environment throughout the product life cycle as the demand of the government policy.

For the task, we concludes twelve design requirements for the vacuum pump after detailed discussion with the staff from the cutting machine industry, which include traditional performance needs and low-carbon demand, and related quality characteristics, as illustrated in Figure 10. Based on the HoQ for the requirements analysis, designers can clearly revels the important QCs and pay more attentions to these attributes during the redesign process, and this study takes the QCs with high relative weight as the retrieval attributes to obtain the similar cases knowledge.

Figure 10.

The mapping analysis for low-carbon requirements and QCs based on the HoQ.

There are 151 product cases in the case base, to obtain the similar cases for knowledge reuse and adaptation, it extracts nine QCs as the retrieval attributes, the total sum of their relative weight are more than 80%. The nine QCs are labeled as A₁-A₉: rate of suction and exhaust (A₁), ultimate pressure (A₂), rated power (A₃), noise (A₄), cost (A₅), low-carbon and environmental impact (A₆), failure rate (A₇), maintenance convenience (A₈), and material recycling and reuse (A₉); A₁-A₅ are numerical attributes, A₆-A₉ are nominal attributes, and each attribute value is provided with the feasible interval, desired interval, and the recommended value according to the fuzzy low-carbon requirements, as shown in Table 5.

Table 5.

Attribute values of the vacuum pump with fuzzy low-carbon requirements.

Retrieval attributes	Feasible interval X	Desired interval X0	Recommended value x
A 1/(L/s)	[50, 250]	[120, 180]	150
A 2/Pa	[0, 9000]	[3000, 7000]	5000
A 3/kW	[7.5, 37]	[11, 18.5]	15
A 4/dB	[68, 82]	[72, 78]	75
A 5/Yuan	[2000, 20000]	[4000, 12000]	8000
A 6	(0.75, 1.0, 1.0)	(0.75, 1.0, 1.0)	(0.75, 0.75, 0.75)
A 7	(0.0, 0.0, 0.30)	(0.0, 0.0, 0.30)	(0.30, 0.30, 0.30)
A 8	(0.75, 1.0, 1.0)	(0.75, 1.0, 1.0)	(0.75, 0.75, 0.75)
A 9	(0.75, 1.0, 1.0)	(0.75, 1.0, 1.0)	(0.75, 0.75, 0.75)

For example, the determination of the attribute value of the rate of suction and exhaust A₁, during the cutting process, the air in the working chamber needs to be exhausted within 5–10 s, and thus the value of A₁ can be derived from the equation (26). ⁵⁴

$$t=2.3K\frac{V}{S}\lg (\frac{P_0}{P})$$

(26)

Where, t is the time to exhaust the air from the working chamber, P₀ is the initial pressure of the working chamber, and P is the final pressure, V is the working chamber volume, S is the value of A₁, K is the coefficient. Based on the calculation result, the rate of the suction and exhaust varies in range [50, 250]L/s. However, for the practical situation, factors such as pumping efficiency, air leakage, cost, and the demand of energy saving need to be considered, thus, in addition to the feasible interval X = [50, 250]L/s, the desired interval is also set as X₀ = [120, 180]L/s, and the recommended value is provided x = 150L/s.

The ultimate pressure (A₂) is determined by the required pressure difference to produce the adsorption force to guarantee that the leather and fabric materials are stably fixed during the cutting process. Suggestion knowledge from the cutting machine customers is that when the pressure difference is more than 10,000 Pa, there is no relative motion between the cutting materials. Thus the working pressure of the vacuum pump is less than 90,000 Pa (set 1atm = 100 kPa); and according to the design principle that the ultimate pressure is one order of magnitude smaller than the working pressure, thus the required ultimate pressure should be less than 9000 Pa. It also should take account of the air leakage factor and the harmful friction in the cutting process caused by the excessive pressure difference; thus, in addition to the feasible interval X = [0, 9000]Pa, the desired interval is set as X₀ = [3000, 7000]Pa, and the recommended value is provided x = 5000 Pa.

For the nominal attributes, the vacuum pump is required with a good low-carbon and environmental performance, low failure rate, good property in maintenance, and high rate of material recycling and reuse. These linguistic terms for the requirements are transformed into the standardized triangular fuzzy numbers, and according to the actual needs from the working situation of the cutting machine, the feasible and desired intervals for attributes can also be designated along with the recommended values.

Generation of the index set of similar cases

The indexing structure of the product case base is constructed based on the strategy that keep approximately equal number of cases for attribute intervals division. Values of each attribute are normalized and sequenced from smallest to largest, respectively. For attributes A₁-A₆, the number of intervals are predefined as 5, and for attributes A₇-A₉, the number of the intervals are predefined as 3, as illustrated in Figures 11 and 12.

Figure 11.

The distribution of cases of attributes A₁-A₅ (η_i = σ_i/σ_(max(i), σ^*= 0.20).

Figure 12.

The distribution of cases of attributes A₆-A₉ (η_i = σ_i/σ_(max(i), σ^*= 0.20).

Take A₁ as an example, attribute A₁ is split into five intervals, in the interval with index 1, there are 33 cases, the attribute value varies in a range (0, 0.019); in the interval with index 2, there are 30 cases, the attribute value varies in a range (0.019, 0.063); in the interval with index 3, there are 32 cases, the attribute value varies in a range (0.063, 0.125); in the interval with index 4, there are 30 cases, the attribute value varies in a range (0.125, 0.25); and in the interval with index 5, there are 26 cases, the attribute value varies in a range (0.25, 1). The distribution of cases for each attribute is estimated based on the criterion in equation (8). It indicates the ratio of σ_i to σ_max(i) is smaller than the threshold σ^* of each attribute, especially for A₇-A₉, although many attribute values are overlapping, the splitting result still make sure these attributes belong to the “good” attributes with basically uniform distribution.

In order to make a clear comparison with the original method for intervals division in the D-HS, that is each interval width is equal to 1/L_i, L_i is the interval number of the i-th attribute. The result of the cases distribution is illustrated in Figure 13, where the interval number of each numerical attribute A₁-A₅ is predefined as 5, and the interval number of each nominal attribute A₆-A₉ is predefined as 10. It indicates that most of cases are falling within the interval with index 1 for A₁-A₅, and the uniformity of cases distribution of attributes A₆-A₉ is also not ideal, which will badly impact the accuracy of the similar cases retrieval.

Figure 13.

The distribution of cases by the original attribute intervals division method: (a) distribution of cases of A₁-A₅ and (b) distribution of cases of A₆-A₉.

The distribution of cases in Figures 11 and 12 demonstrates that the predefined interval number of each attribute is acceptable, and accordingly, the coordinates system of the case base can be established as shown in Table 6; and then the coordinates of each case should be designated according to the attribute values. Take the product case c₇ as the example, the coordinates are given as described in Table 7.

Table 6.

The coordinates system of the case base.

Ai	Interval 1		Interval 2		Interval 3		Interval 4		Interval 5
	l 1	Range	l 2	Range	l 3	Range	l 4	Range	l 5	Range
A 1	1	0, 0.019	2	0.019, 0.063	3	0.063, 0.125	4	0.125, 0.25	5	0.25, 1
A 2	1	0, 3E-6	2	3E-6, 0.001	3	0.001, 0.130	4	0.130, 0.20	5	0.20, 1
A 3	1	0, 0.019	2	0.019, 0.05	3	0.05, 0.094	4	0.094, 0.188	5	0.188, 1
A 4	1	0, 0.175	2	0.175, 0.25	3	0.25, 0.40	4	0.40, 0.50	5	0.50, 1
A 5	1	0, 0.022	2	0.022, 0.049	3	0.049, 0.089	4	0.089, 0.179	5	0.179, 1
A 6	1	0, 0.50	2	0.50, 0.65	3	0.65, 0.70	4	0.70, 0.75	5	0.75, 1
A 7	1	0, 0.25	2	0.25, 0.30	3	0.30, 1	–	–	–	–
A 8	1	0, 0.65	2	0.65, 0.70	3	0.70, 1	–	–	–	–
A 9	1	0, 0.65	2	0.65, 0.70	3	0.70, 1	–	–	–	–

Table 7.

The coordinates l_i(c₇) of the product case c₇.

Attributes	A 1	A 2	A 3	A 4	A 5	A 6	A 7	A 8	A 9
Values	150	6666	15	74	3200	0.75	0.30	0.75	0.75
Normalization	0.188	0.333	0.188	0.350	0.027	0.75	0.30	0.75	0.75
li(c7)	4	5	4	3	2	4	2	3	3

Query case is the target case in the retrieval algorithm, the coordinates of query case are also designated as shown in Table 8 based on the recommended value of each attribute provided in Table 5.

Table 8.

The coordinates of l_i(q) of the query case q.

Attributes	A 1	A 2	A 3	A 4	A 5	A 6	A 7	A 8	A 9
Values	150	5000	15	75	8000	0.75	0.30	0.75	0.75
Normalization	0.188	0.25	0.188	0.370	0.075	0.75	0.30	0.75	0.75
li(q)	4	5	4	3	3	4	2	3	3

The cases labeled with the same coordinates as the query case are extracted to constitute the index set S, and the number of attributes which are matched with the query case q for each case c _i , SMA(c _i , q) is estimated, as shown in Table 9. From Table 9, it indicates that there is no one case can completely match with the query case in each attribute, that is, the index set S is an empty set. In this situation, the threshold Const is set, Const = 4, and to obtain the index set S_b = {c _i : SMA(c _i , q)≥4}; thus the index set S₀ = {c _i : c _i ∈SUS_b} = S_b; there are 25 cases in S₀, the attributes of which are listed in Table 10.

Table 9.

The result of SMA(c _i , q) between the product case and query case.

SMA(c i , q)	Number of cases	Cases corresponding to the SMA(c i , q)
9	0	–
8	1	c7;
7	5	c6, c8, c39, c43, c44;
6	1	c38;
5	4	c9, c17, c23, c134;
4	14	c5, c18, c22, c24, c25, c47, c87, c114, c123, c130, c133, c144, c145, c150;
3	29	c1, c2, c3, c4, c10, c11, c26, c27, c28, c36, c37, c45, c46, c88, c94, c97, c98, c102, c103, c106, c108, c109, c110, c115, c126, c127, c138, c139, c149;
2	37	c12, c13, c14, c15, c16, c29, c30, c31, c32, c34, c41, c42, c48, c51, c52, c55, c59, c61, c63, c66, c67, c75, c85, c91, c104, c105, c107, c111, c112, c113, c124, c129, c135, c136, c141, c143, c151;
1	39	c19, c21, c33, c35, c40, c53, c54, c56, c57, c60, c62, c64, c65, c68, c69, c70, c71, c74, c76, c77, c78, c84, c86, c93, c95, c96, c100, c101, c117, c118, c121, c122, c125, c128, c131, c132, c140, c141, c142;
0	21	c20, c49, c50, c58, c72, c73, c79, c80, c81, c82, c83, c89, c90, c92, c99, c116, c119, c120, c137, c146, c147, c148.

Table 10.

Cases in the index set S₀.

Cases/c i	A 1/(L/s)	A 2/Pa	A 3/kW	A 4/dB	A 5/Yuan	A 6	A 7	A 8	A 9
c7	150	6666	15	74	3200	0.75*	0.3*	0.75*	0.75*
c6	100	6666	11	72	2500	0.75*	0.3*	0.75*	0.75*
c8	200	6666	18.5	75	4900	0.75*	0.3*	0.75*	0.75*
c39	138	3300	15	74	6800	0.75*	0.3*	0.65*	0.75*
c43	111	3300	11	73	5500	0.75*	0.3*	0.65*	0.75*
c44	138	3300	15	74	7200	0.75*	0.3*	0.65*	0.75*
c38	111	3300	11	73	4050	0.75*	0.3*	0.65*	0.75*
c9	250	6666	30	78	6000	0.7*	0.3*	0.75*	0.75*
c17	100	3333	15	74	2700	0.65*	0.3*	0.75*	0.75*
c23	95	3300	11	76	3800	0.75*	0.25*	0.75*	0.75*
c134	70	2600	7.5	75	7500	0.5*	0.3*	0.75*	0.75*
c5	50	6666	5.5	67	1950	0.75*	0.25*	0.75*	0.75*
c18	200	3333	22	78	5400	0.65*	0.3*	0.75*	0.75*
c22	65	3300	7.5	72	3600	0.75*	0.25*	0.75*	0.75*
c24	113	3300	15	80	4500	0.75*	0.25*	0.7*	0.75*
c25	130	3300	15	80	4850	0.75*	0.25*	0.7*	0.75*
c47	191	8106	22	82	9400	0.65*	0.3*	0.7*	0.65*
c87	150	1	15	84	18000	0.5*	0.3*	0.65*	0.75*
c114	120	10	15	78	85000	0.85*	0.3*	0.7*	0.75*
c123	150	0.05	2.2	80	6200	0.7*	0.3*	0.75*	0.7*
c130	166	5489.9	15	82	7150	0.7*	0.35*	0.7*	0.7*
c133	50	2600	5.5	72	5200	0.5*	0.3*	0.75*	0.75*
c144	150	2600	11	80	14000	0.7*	0.3*	0.7*	0.75*
c145	200	2600	15	80	18300	0.7*	0.3*	0.7*	0.75*
c150	50	300	7.5	76	8200	0.7*	0.3*	0.7*	0.75*

Note: 0.25* = (0.1, 0.25, 0.4); 0.3* = (0.2, 0.3, 0.4); 0.35* = (0.2, 0.35, 0.5); 0.5* = (0.4, 0.5, 0.6); 0.65* = (0.5, 0.65, 0.8); 0.7* = (0.6, 0.7, 0.8); 0.75* = (0.65, 0.75, 0.85); 0.85* = (0.75, 0.85, 0.95).

Estimation of the similarity and adaptability

When the index set S₀ is generated, cases in the S₀ still need to be estimated from two perspectives, the similarity and adaptability according to the low-carbon requirements. For the estimation of similarity, A₁-A₅ belong to fuzzy numerical attributes, thus it adopts the dependent function to calculate the local similarity of each attribute; and for the nominal attributes A₆-A₉, the estimation model based on the membership function as constructed in equation (17) is used to calculate the local similarity of each attribute. The result of the similarity estimation for each case is listed in Table 11.

Table 11.

The result of similarity estimation for cases in S₀.

c i	sim2(Ai, qi), i  = 1, 2, 3, 4, 5;					sim3(Ai, qi), i  = 6, 7, 8, 9				SIM(c i , q)
	A 1	A 2	A 3	A 4	A 5	A 6	A 7	A 8	A 9
c25	1.000	1.000	1.000	0.607	1.000	0.229	0.592	0.057	0.229	0.659
c23	0.700	1.000	1.000	1.000	0.949	0.229	0.592	0.229	0.229	0.654
c39	1.000	1.000	1.000	1.000	1.000	0.229	0.167	0.050	0.229	0.653
c44	1.000	1.000	1.000	1.000	1.000	0.229	0.167	0.050	0.229	0.653
c24	0.912	1.000	1.000	0.607	1.000	0.229	0.592	0.057	0.229	0.646
c7	1.000	1.000	1.000	1.000	0.780	0.229	0.167	0.229	0.229	0.644
c43	0.882	1.000	1.000	1.000	1.000	0.229	0.167	0.050	0.229	0.635
c38	0.882	1.000	1.000	1.000	1.000	0.229	0.167	0.050	0.229	0.635
c8	0.753	1.000	1.000	1.000	1.000	0.229	0.167	0.229	0.229	0.631
c114	1.000	0.369	1.000	1.000	0.0001	0.914	0.167	0.057	0.229	0.589
c6	0.751	1.000	1.000	1.000	0.550	0.229	0.167	0.229	0.229	0.583
c18	0.753	1.000	0.827	1.000	1.000	0.050	0.167	0.229	0.229	0.579
c17	0.751	1.000	1.000	1.000	0.619	0.050	0.167	0.229	0.229	0.562
c144	1.000	0.875	1.000	0.607	0.779	0.057	0.167	0.057	0.229	0.559
c130	1.000	1.000	1.000	0.368	1.000	0.057	0.148	0.057	0.057	0.558
c22	0.452	1.000	0.368	1.000	0.896	0.229	0.592	0.229	0.229	0.526
c145	0.753	0.875	1.000	0.607	0.455	0.057	0.167	0.057	0.229	0.488
c9	0.368	1.000	0.537	1.000	1.000	0.057	0.167	0.229	0.229	0.483
c47	0.844	0.641	0.827	0.368	1.000	0.050	0.167	0.057	0.050	0.472
c134	0.490	0.875	0.368	1.000	1.000	0.000	0.167	0.229	0.229	0.456
c87	1.000	0.368	1.000	0.223	0.472	0.000	0.167	0.050	0.229	0.432
c133	0.368	0.875	0.208	1.000	1.000	0.000	0.167	0.229	0.229	0.416
c123	1.000	0.368	0.082	0.607	1.000	0.057	0.167	0.229	0.057	0.400
c150	0.368	0.407	0.368	1.000	1..000	0.057	0.167	0.057	0.229	0.384
c5	0.368	1.000	0.208	0.287	0.36	0.229	0.592	0.229	0.229	0.378

Make a comparison between Tables 9 and 11, it indicates that the cases with high number of matched attributes cannot guarantee high similarity. For instance, case c₇ has the highest number of matched attributes, but the similarity of c₇ is not the best. This situation is same to the cases c₆, c₈. Conversely, cases c₂₅, c₂₃, and c₂₄ have relative high similarity although the number of matched attributes of which are not very high. Therefore, the strategy of similarity matching for attributes can effectively prevent similar cases with high similarity from being missed.

Estimation of the adaptability for cases is based on the information content of each attribute, which is determined by the design range and system range. The design range of each attribute is obtained from the requirement analysis, as explained in Table 5, and the range values of numerical attributes need to be transformed into the standardized triangular fuzzy number. The system range of each attribute reveals the maximum range for adaptation of the series product cases under the existing technical conditions. For example, the ultimate pressure of the SK series vacuum pump is generally at the range 6000–8000 Pa, but it can reach 3000 Pa or lower by installing a fore-stage nozzle. Here, it takes case c₂₅ as the example to describe the estimation of the adaptability for each attribute.

In Table 12, the normalized attribute values of design range and system range of the case c₂₅ are provided. For attributes A₁(rate of suction and exhaust), A₂(ultimate pressure), A₃(rated power), they follow the fuzzy triangular distribution; A₂ belongs to the Cost type attribute, and the system range is a crisp value, thus the system area of A₂ is constructed based on the principle in Figure 9, here with the Cost type attribute. In Figure 14, the calculation models of information content I_Ai with respect to A₁, A₂, A₃ of case c₂₅ are illustrated, respectively.

• I_{A 1} = log₂(1/ (A_CR /A_STFN(SR) )) = log₂(1/0.28) = 1.836;

• I_{A 2} = log₂(1/ (A_CR /A_STFN(SR) )) = log₂(1/0.747) = 0.421;

• I_{A 3} = log₂(1/ (A_CR /A_STFN(SR) )) = log₂(1/0.679) = 0.558.

Table 12.

Design range of the requirements and system range of c₂₅.

Ai	Design range	System range
A 1	0.15, 0.188, 0.225	0.081, 0.163, 0.20
A 2	0.15, 0.25, 0.35	0.165, 0.165, 0.165
A 3	0.138, 0.188, 0.231	0.094, 0.188, 0.231
A 4	0.30, 0.45	0.375, 0.75
A 5	0.035, 0.116	0.023, 0.067
A 6	0.75, 1.0, 1.0	0.65, 0.75, 0.85
A 7	0.0, 0.0, 0.30	0.10, 0.25, 0.40
A 8	0.75, 1.0, 1.0	0.60, 0.70, 0.80
A 9	0.75, 1.0, 1.0	0.65, 0.75, 0.85

Figure 14.

Information content estimation for attributes A₁-A₃: (a) information content I_A1, (b) information content I_A2, and (c) Information content I_A3.

For attributes A₄(noise), A₅(cost), they follow the uniform distribution, the calculation models of the information content for A₄ and A₅ are illustrated in Figure 15.

• I_{A 4} = log₂(1/ (L_CR/L_SR )) = log₂(1/0.20) = 2.322;

• I_{A 5} = log₂(1/ (L_CR/L_SR )) = log₂(1/0.727) = 0.460.

Figure 15.

Information content estimation for attributes A₄ and A₅: (a) information content I_A4, and (b) information content I_A5.

For nominal attributes A₆-A₉, the linguistic term of each attribute value is transformed into the standardized triangular fuzzy number, and the calculation models of the information content for A₆-A₉ are illustrated in Figure 16.

• I_{A 6} = log₂(1/ (A_CR/_ASTFN(SR) )) = log₂(1/0.143) = 2.806;

• I_{A 7} = log₂(1/ (A_CR/A_STFN(SR) )) = log₂(1/0.296) = 1.756;

• I_{A 8} = log₂(1/ (A_CR /A_STFN(SR) )) = log₂(1/0.036) = 4.796;

• I_{A 9} = log₂(1/ (A_CR/A_STFN(SR) )) = log₂(1/0.143) = 2.806.

Figure 16.

Information content estimation for attributes A₆-A₉: (a) information content I_A6, (b) information content I_A7, (c) information content I_A8, and (d) information content I_A9.

Then the information content of case c₂₅ is obtained:

• $\begin{array}{c}{I}_{c25}=({\omega }_{A1},{{\omega }_A}_2,{{\omega }_A}_3,{{\omega }_A}_4,{{\omega }_A}_5,{{\omega }_A}_6,{{\omega }_A}_7,{{\omega }_A}_8,{{\omega }_A}_9)\hfill \\ •{({I_A}_1,{I_A}_2,{I_A}_3,{I_A}_4,{I_A}_5,{I_A}_6,{I_A}_7,{I_A}_8,{I_A}_9)}^T\hfill \\ =(0.149,0.106,0.136,0.080,0.105,0.157,0.087,0.082,0.098)\hfill \\ •{(1.836,0.421,0.558,2.322,0.460,2.806,1.756,4.796,2.806)}^T\hfill \\ =1.890.\hfill \end{array}$

The information content of the remaining cases in S₀ also can be calculated with the same method, the result is shown in Table 13. In addition, in order to consider the comprehensive influence of similarity and adaptability of the product cases, the comparison analysis histogram of the similarity and information content of each case in the index set S₀ is generated, as illustrated in Figure 17.

Table 13.

The result of information content estimation for cases in S₀.

C i	IA 1	IA 2	IA 3	IA 4	IA 5	IA 6	IA 7	IA 8	IA 9	I ci
c25	1.836	0.421	0.558	2.322	0.460	2.806	1.756	4.796	2.806	1.890
c23	2.988	0.421	1.910	1.114	0.460	2.806	1.756	2.806	2.806	1.986
c39	1.966	0.421	0.558	0.737	0.615	2.806	3.000	5.573	2.806	1.970
c44	1.966	0.421	0.558	0.737	0.474	2.806	3.000	5.573	2.806	1.956
c24	2.988	0.421	0.558	1.806	0.460	2.806	1.756	4.796	2.806	2.020
c7	1.221	2.556	0.558	1.586	0.480	2.806	3.000	2.806	2.806	1.913
c43	3.059	0.421	1.910	1.224	0.571	2.806	3.000	5.573	2.806	2.352
c38	3.059	0.421	1.910	1.000	0.648	2.806	3.000	5.573	2.806	2.342
c8	3..023	2.556	0.548	1.224	0.460	2.806	3.000	2.806	2.806	2.149
c114	2.857	1.977	0.910	0.737	α	0.806	3.000	4.796	2.806	α
c6	3.265	2.556	1.910	1.224	0.927	2.806	3.000	2.806	2.806	2.419
c18	3.023	0.421	1.514	1.737	0.460	4.158	3.000	2.806	2.806	2.307
c17	3.265	0.421	0.910	1.586	0.808	4.158	3.000	2.806	2.806	2.286
c144	1.224	0.816	1.910	1.486	1..552	4.796	3.000	4.796	2.806	2.493
c130	1.214	0.396	0.910	3.699	1.351	4.796	3.756	4.796	4.796	2.728
c22	3.427	0.421	2.434	1.415	0.619	2.806	1.756	2.806	2.806	2.163
c145	3.023	0.816	0.558	1.486	2.523	4.796	3.000	4.796	2.806	2.679
c9	3.857	2.556	3.878	1.224	0.514	4.796	3.000	2.806	2.806	3.044
c47	2.699	α	3.718	2.322	0.876	5.573	3.000	4.796	5.573	α
c134	3.251	0.816	2.434	1.221	0.533	α	3.000	2.806	2.806	α
c87	0.941	1.977	0.558	3.171	3.699	α	3.000	5.573	2.806	α
c133	3.836	0.816	3.012	1.224	0.514	α	3.000	2.806	2.806	α
c123	0.941	1.977	α	2.322	0.456	4.796	3.000	2.806	4.796	α
c150	3.836	1.862	1.565	1..000	0.626	4.796	3.000	4.796	2.806	2.810
c5	3.836	2.556	3.012	1.114	1.806	2.806	1.756	2.806	2.806	2.629

Note: α denotes that the corresponding information content approaches infinity.

Figure 17.

Comparison analysis of the similarity and information content for cases in S₀.

In Figure 17, it indicates that cases c₂₅, c₂₃, c₃₉, and c₄₄ have high similarity and small information content; and for the case c₇, although the similarity is not high, the information content is small; therefore, these cases are considered as the similar cases knowledge for reuse. For cases c₄₃, c₃₈, c₈, although the similarity is not low, the information content is relative large; however, the information content of cases c₁₁₄, c₄₇, c₁₃₄, c₈₇, c₁₃₃, c₁₂₃ are all approaching infinity, it indicates that there are one or more requirement attributes cannot be satisfied by adaptation under the current technical conditions.

Meanwhile, the local attribute similarity and information content can guide the designer to make targeted improvement for the candidate cases, and it also could indicate the direction of technology updates to some extent. For instance, the attributes A₄ (noise) and A₈ (maintenance convenience) of c₂₅ cannot be completely satisfied, and the related information content of A₄ and A₈ are relative high. The designer maybe adopt new scheme to enhance the local similarity rather than make adaption on the original structure, such as install the noise reduction device, and replace the mechanical sealing component with a simple and effective packing sealing component to reduce the failure rate and facilitate the maintenance. Case c₁₁₄ is a screw vacuum pump, which has good environmental performance and high work efficiency, but its cost is far beyond the scope of the cost constraint; and thus reducing the cost with new technologies is a key issue. The information content of attribute A₆ (low-carbon and environmental impact) of cases c₈₇, c₁₃₃, c₁₃₄ are all approaching infinity, it reveals that the traditional technology can no longer meet the customer’s low-carbon and environmental needs, and it is urgent to implement the green manufacturing technologies to update the product cases to meet the market demand and the government environmental regulations.

Discussion

In this paper, an effective similarity determination model which integrates the representation of fuzzy requirement information, the efficient indexing approach, and the evaluation method for the similarity and adaptability into a holistic similar cases retrieval methodology is proposed to support low-carbon product design. A case study concerning the similar cases retrieval for the vacuum pump low-carbon design is conducted to verify the applicability of the proposed methodology; and the theoretical and practical implications of our research work are discussed below.

In the earlier design stage, designers can basically specify the accurate values for the very familiar requirement attributes. However, when taking account of the low-carbon requirements, it requires the reconfiguration for corresponding quality characteristics, which possibly contributes to complex design conflicting problems as the coupling property between the quality characteristics, namely, one design factor about the low-carbon requirements is improved, while the related factor is deteriorated. Consequently, even if there are experienced domain experts for the design task, sometimes they cannot grasp deep knowledge and precisely address the complex low-carbon design problems with incomplete design information, but they can roughly know the possible range of design attribute values, and the tolerable differences between the query case attribute values and that of the candidate cases to ensure the designated parameter ranges are still satisfying the design purpose. Therefore, for the fuzzy requirement attributes, this study adopts the fuzzy set theory to represent the incomplete and uncertain design information, and it defines the feasible interval and desired interval to make flexible for the designers in decision-making to the design parameters.

For the improved D-HS method, it adopts the recommended attribute values to assign the indexing coordinates for the query case rather than uses the design range of each attribute to construct the index set mainly due to the following considerations: if it uses design ranges of attributes to search the similar cases, then it must identify all cases in the memory whether they are falling within the index set every time with a new requirement, and thus it loses the meaning of the intervals division. Conversely, with the recommended values, the coordinates of the query case can be designated, and cases having the same coordinates with that of the query case are extracted to build the index set; in addition, when the coordinates of the current query case are same as that of the previous one, similar cases in the index set saved in the database can be immediately reused.

In our study, we propose the strategy that the number of cases is approximately equal for each interval of attributes to make sure the cases can be evenly distributed. In case study, it predefine the interval number of the i-th attribute L_i, however, the determination of L_i is actually an iterative process to satisfy the criteria in equation (8). Besides, for two neighboring cluster cases, if the average values of the i-th attribute of the two clusters are quite different, then it is necessary to construct the attribute interval for two clusters separately, and re-designate the L_i, even though only the sum of the number of cases of two clusters can satisfy the condition of uniform distribution of the cases. The strategy of similarity matching for attributes proposed in the paper provides a supplementary index set S_b by defining the required number of matched attributes, namely the threshold Const, which effectively avoids the omission of the similar cases. For example, in the case study, cases c₂₅,c₂₄,c₂₃ are included in the initial index set S₀, the number of the matched attributes of these three cases are not high, they can be seen as the edge cases as illustrated in Figure 3(b); although they do not fall within the same intervals with the query case, they have the high similarity, and luckily, they are not neglected.

The effective estimation models of the similarity and adaptability of the cases are respectively established in our study. For the similarity calculation of precise numerical attributes, it adopts the traditional distance measurement to calculate the local similarity; for the fuzzy numerical attributes, it uses the dependent function to construct the local similarity evaluation model, as it can better distinguish the degree to which the attribute values of the candidate case meet the design requirements. That is, when the value of dependent function of the attribute is equal or more than one, it indicates the attribute value is completely meet the demand; when the value of dependent function is between zero and one, then the attribute value is acceptable; while when the value of dependent function is less than zero, the attribute value cannot be satisfied. For the nominal attributes, the traditional method is to use the binary estimation model, 0 and 1. For instance, the linguistic term of the attribute low-carbon and environmental impact (A₆) of the query case is required “good”, when the linguistic term of A₆ of the candidate case is “good”, then the similarity is equal to 1, otherwise is 0; however, when the linguistic term of A₆ of the candidate case is “very good”, “fair”, “poor”, or “very poor”, it cannot clearly distinguish the similarity of these values. In our study, the linguistic terms of the nominal attribute are transformed into the standardized fuzzy triangular numbers, and the similarity estimation model is constructed based on the ratio of area of the common range to the area of the design range.

The estimation model of the adaptability of the candidate case is proposed based on the information content presented in the Information Axiom. For one certain attribute of the case, the smaller the information content, the better the adaptability performance, because the system range of this attribute has a high degree of overlap with the design range, which indicates it has a high probability to meet the design requirements through modification; and the information content of the case is the sum of the information content of each attribute. To estimate the information content of each attribute, it is important to identify its system range that is based on the existing manufacturing technologies, design specification, related laws and regulations, etc. Thus it is hard to delineate a clear system range of the attribute with the incomplete design information. In this study, for the numerical attribute, the upper and lower boundary values that the attribute can reach under the current technical conditions of the series product cases are defined as the system range. For example, the rate of suction and exhaust of the SK series vacuum pump can be promoted by enhancing the rated power or increasing the number of impellers, and thus to obtain the maximum and minimum value of the rate of suction and exhaust to construct the system range. For the nominal attribute, the system range is defined by evaluating the corresponding sub-factors. For instance, to define the system range of the attribute A₆, it needs to evaluate its sub-factors, such as the life cycle carbon footprint of the product, assemble and disassemble for the components reuse, the material recycling and disposal methods; the linguistic term of the attribute is final obtained by scoring the each sub-factor, and then the system range of the attribute is specified by comparing series product cases.

Conclusion and further research

In low-carbon product design, it takes account of the low-carbon requirements at preliminary design stage in consideration of the reduction of the carbon emissions throughout the product life cycle. However, as the coupling properties between the quality characteristics, it inevitably causes the design conflicting problems when incorporating the low-carbon indicators, which needs to go through multiple backtracking design processes to coordinate the contradictions. Thus, it is a good strategy for exploiting the prior design knowledge and experience to solve the similar design problems. In our study, we propose an effective similarity determination method to obtain acceptable cases knowledge for reuse, similarity here is not only referred to the similarity between the attribute values of the query case and that of the candidate case, but the similarity of the system range of the candidate case and the design range of the query case.

The incomplete and uncertain requirements information is represented based on the fuzzy set theory. The efficient D-HS indexing approach is improved with two valid strategies, the strategy that keep approximately equal number of the cases for intervals division, and the strategy of similarity matching for attributes; the improved D-HS method can effectively guarantee the basically even distribution of cases in each attribute interval, and avoid the omission of the similar cases when the location of the query case is close to the boundaries of the design ranges. In terms of the similarity estimation, it adopts the classical distance measurement and dependent function to calculate the local similarity for the precise numerical attribute and the fuzzy numerical attribute; and for the fuzzy nominal attribute, the similarity is based on the ratio of area of the common range to the area of the design range. For the adaptability estimation, it adopts the information content presented in the Information Axiom to construct the evaluation model for the adaptability of each attribute, the smaller the total information content, the better the adaptability of the candidate case. Thus, the cases having high similarity and small information content can be considered as the acceptable cases for knowledge reuse and adaptation.

In our proposed methodology, the information content of each attribute depends on the corresponding design range and system range, in the further research, it is necessary to use the sensitivity analysis and Monte Carlo analysis methods to more reasonably determine the boundaries of the design range and system range of each attribute, as well as the distribution of its fuzzy membership function. In addition, it is essential to adopt the contradictory problems solving methods, such as the TRIZ, Extenics, to address the design conflicts during the adaptation for the similar cases, which facilitate the generation of the low-carbon innovative design schemes.