• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A new interpretable fault diagnosis method based on belief rule base and probability table

    2023-04-22 02:04:52ZhichoMINGZhijieZHOUYouCAOShuiwenTANGYunCHENXioxiHANWeiHE
    CHINESE JOURNAL OF AERONAUTICS 2023年3期

    Zhicho MING, Zhijie ZHOU,*, You CAO, Shuiwen TANG, Yun CHEN,Xioxi HAN, Wei HE

    a High-Tech Institute of Xi’an, Xi’an 710025, China

    b School of Computer Science and Information Engineering, Harbin Normal University, Harbin 150025, China

    KEYWORDS Aerospace relay;Belief rule base;Expert knowledge;Fault diagnosis;Interpretability constraints

    Abstract It is vital to establish an interpretable fault diagnosis model for critical equipment.Belief Rule Base (BRB) is an interpretable expert system gradually applied in fault diagnosis.However,the expert knowledge cannot be utilized to establish the initial BRB accurately if there are multiple referential grades in different fault features.In addition, the interpretability of BRB-based fault diagnosis is destroyed in the optimization process, which reflects in two aspects: deviation from the initial expert judgment and over-optimization of parameters.To solve these problems, a new interpretable fault diagnosis model based on BRB and probability table, called the BRB-P, is proposed in this paper.Compared with the traditional BRB,the BRB-P constructed by the probability table is more accurate.Then, the interpretability constraints, i.e., the credibility of expert knowledge, the penalty factor and the rule-activation factor, are inserted into the projection covariance matrix adaption evolution strategy to maintain the interpretability of BRB-P.A case study of the aerospace relay is conducted to verify the effectiveness of the proposed method.

    1.Introduction

    Critical equipment,such as the electro-hydraulic servo and the aerospace relay, plays a significant role in the operation of complex systems.Once the fault occurs, if unable to be ruled out timely, it may cause considerable economic losses.Fault diagnosis is an effective measure to diagnose abnormal working states.

    According to different information used in fault diagnosis,three types of methods are usually adopted, including the model-based method,1–2the data-based method3and the expertise-based method.4In the model-based methods, the fault diagnosis model is established by analyzing the fault mechanism of equipment.5Good transparency,precise mechanism and strong interpretability constitute its main advantages.Nonetheless, this method depends on establishing an accurate fault mechanism model.If there are a large number of fault samples, a data-based method such as the neural network can be employed to describe the causal relationship between fault features and fault modes.6–9The model established in this method works in black-box perspective and is entirely independent from system mechanism.10When there is a lack of fault samples, and the system mechanism model is challenging to obtain, the expert knowledge and historical experience can be utilized to construct the fault model, which is the expertise-based method.11However,the diagnostic accuracy is relevant to the expert knowledge due to the subjective deviation of experts.

    Based on the above analysis, it is necessary to establish an excellent fault diagnosis model for critical equipment that generally has the following characteristics:

    (1) Considering the complex structure, the accurate mechanism model is hard to be established.

    (2) In most cases, the equipment works in the normal state with limited fault samples,and the sample characteristic tends to be‘‘low density”.12

    (3) Due to the particularity, there is a demand to establish an understandable fault diagnosis model for critical equipment, which increases the trust between humans and the model.13

    Based on the above characteristics,this paper aims to establish an interpretable fault diagnosis model for critical equipment.Interpretability in this paper means that the operation of a model can be understood by humans.14–15

    Yang et al.16proposed the Belief Rule Base (BRB) inference methodology using the evidential reasoning approach,which is a typical expert system.The rule of BRB has the understandable knowledge form of‘‘IF-THEN”.17BRB is developed based on the evidence theory.18The uncertainty can be measured by belief function rather than probability in evidence theory.There is no need for a large number of historical samples to determine parameters such as the prior distribution.The quantitative information and qualitative knowledge can be effectively combined to construct the model.The evidence theory has been applied in information fusion,expert system and decision analysis based on the above advantages.There is a general knowledge representation scheme based on the belief structure for various types of inputs, such as quantitative data and qualitative information.Thus, BRB has the ability of complex nonlinear modeling based on semantic information.19When the fault samples are limited, or the mechanism model is difficult to be obtained, the BRB can be established based on the expert knowledge or historical experience.The parameters of BRB can be modified during the optimization process.The reasoning engine of BRB is the Evidential Reasoning (ER) approach that is used to fuse the activated rules.The ER approach is a generalization form of the traditional Bayesian inference rule.20BRB has been widely employed in various fields,such as disease diagnosis,engineering manufacturing, industrial management,21–25etc.

    However, there are still two problems in the BRB-based fault diagnosis.

    (1) In the modeling process,many fault features and modes exist in the system mentioned above.The ability of experts to deal with complex systems is limited.The expert knowledge cannot be accurately embedded in BRB when there are a large number of rules.

    (2) In the optimization process, the interpretability of BRB may be destroyed,and the trust degree between humans and the model is reduced.On the one hand, the initial BRB is given and effectively explained by experts.Inevitably, the optimized BRB may deviate from the initial expert judgment so that the optimized rules cannot be understood by experts.On the other hand, the parameters of BRB are over-optimized.26In other words, all parameters of a fault diagnosis model are optimized previously whether they work or not during the training process.

    To accurately establish the model of complex systems,Chang et al.17proposed a disjunctive belief rule base modeling method, in which BRB model was constructed under the disjunctive assumption rather than the conjunctive assumption.A new approximate BRB model,named ABRB,was proposed by Cao et al.27The input of a rule is converted into a single input form by introducing the independency factor in ABRB.Reducing the number of rules in these methods is convenient to embed expert knowledge.In addition, Zhang et al.28proposed a method of automatically generating initial parameters for BRB by means of cloud model and partial rules given by experts.The establishment of BRB has been conducted with in-depth research, but the interpretability of the model has not been mentioned in the above literature.With the attention paid to the interpretability of algorithms, the research on the interpretability of BRB has made some progress.Yang et al.29put forward a viewpoint that the parameters of BRB should be fine-tuned in optimization.Zhou et al.19stated the intrinsic interpretability of BRB and proposed a new health-state assessment model based on BRB with interpretability.30Cao et al.26summarized the interpretability of BRB and proposed a modified optimization algorithm with the interpretability constraints.The interpretability of BRB was defined, and the evaluation indicator of interpretability was proposed by You et al.31

    In this paper,a new fault diagnosis model,called BRB-P,is constructed based on BRB with a probability table and three interpretability constraints.It considers both the establishment of BRB and the maintenance of model interpretability.The expert knowledge is expressed as the probability table to establish the initial fault diagnosis model.In this way, the number of initialized parameters is decreased to be convenient for applying the expert knowledge.Meanwhile, the completeness of the rule base can be preserved.Then, the parameters of the initial BRB-P are optimized by the limited fault samples to reduce the error brought by the subjectivity of experts.In this paper,the Projection Covariance Matrix Adaption Evolution Strategy (P-CMA-ES) optimization algorithm is utilized to optimize the parameters of BRB-P,and interpretability constraints are added to maintain the interpretability of the fault diagnosis model.Specifically, there are three interpretability constraints, introducing the credibility of expert knowledge,the penalty factor and rule-activation factor to the optimization algorithm, respectively.The credibility of expert knowledge and the penalty factor can avoid deviating from the initial judgment of experts in optimization.The problem of parameters over-optimization can be solved by introducing the rule-activation factor.

    The remainder of this paper is organized as follows.In Section 2, the relevant problems of BRB for fault diagnosis are discussed, and the framework of BRB-P is briefly introduced.In Section 3,the establishment and reasoning method of BRBP is proposed.The optimization of BRB-P with interpretability constraints is conducted in Section 4.A case study of the relay fault diagnosis is performed to verify the effectiveness of the proposed method in Section 5.This paper is concluded in Section 6.

    2.Problem formulation and framework of BRB-P model

    In this section, the problem formulation of fault diagnosis based on BRB is presented.Then, the basic framework of BRB-P is constructed.

    2.1.Problem formulation of fault diagnosis based on BRB

    Suppose there are M fault features and N fault modes in the system.The BRB is utilized to describe the relationship between them.A rule base is composed of a series of rules,which are the comprehensible form of‘‘IF-THEN”.The structure of the k th rule is expressed as follows:

    Problem 1.All the referential values of fault features should be traversed when establishing a BRB-based fault diagnosis model.Thus,the number of rules in BRB can be calculated by

    where cmis the number of referential grades of the m th fault feature.As can be seen in Eq.(2), K increases exponentially with M.The initialized parameters given by experts include the belief degrees, fault feature weights and rule weights.The total number of initialized parameters can be calculated by

    Taking the fault diagnosis of the relay as an example, 4 fault features and 6 fault modes are the input and output variables of the fault diagnosis model.3 referential grades are assigned for each fault feature.The k th rule is shown in Eq.(4).There are 81 (34=81) rules in BRB and 571 initialized parameters in the traditional BRB, which should be given by experts to maintain the interpretability of BRB.However,the ability to process a great deal of information for experts is limited.32–33Zhang et al.28pointed out that experts were hard to apply knowledge and historical experience to construct the BRB-based model when there were a large number of parameters.In other words, BRB cannot be accurately established by experts.It is necessary to propose a new mechanism to transform the expert knowledge into the initial BRB.

    Problem 2.After the initial BRB is established by experts, the optimization process is employed to reduce the model error brought by the subjectivity of experts.34The target of the current optimization is to improve the accuracy of the model while maintaining interpretability is ignored.In this process,the interpretability of BRB is destroyed in two aspects: the deviation from the initial judgment of experts and the overoptimization of parameters.The initial judgment of experts is a key element of the interpretability of BRB.In other words,the model result is consistent with the experts’cognition.In this paper,we assume that the experts are authoritative.It is noted that the target of optimization should be to fine-tune the parameters of the BRB model based on the initial judgment.However,the optimized rules may be different from the initial rules to a large extent,which means that the deviation from the expert knowledge occurs.35–36The final result shows that experts are unable to understand the optimized BRB.For instance, Liu et al.37established the BRB model for fault diagnosis of an airborne missile,whose partial rules are shown in Table 1.Refrigeration noise voltage (x1) and refrigerating flow (x2) are selected as fault features.Taking the 8th rule as an example, the belief degrees given by experts are (0.8600,0.0000, 0.1400).Fault D1occurs when the input is exactly the reference values of fault features in the 8th rule.However, the corresponding result of the optimized BRB indicates that the fault D3occurs in the same case.It is noted that the optimized rules are inconsistent with the experts’cognition.

    In addition, the over-optimization problem of parameters means that the parameters of unactivated rules are optimized,which weakens the interpretability of BRB.30The existing normal or fault states are included in the fault diagnosis model.However, the limited samples may contain the partial states of a system.In other words, some rules unactivated in the training process never play roles in the fusion process.It is irrational that the parameters of a fault diagnosis model are totally optimized previously whether they work or not.

    2.2.Framework of fault diagnosis model based on BRB-P

    To solve the problems mentioned above,a new fault diagnosis model based on BRB-P is proposed, whose framework is depicted in Fig.1.The following part highlights the functions of each module in Fig.1.

    In the modeling process,the initial BRB-P is established by the probability table rather than being given directly by expertsto solve Problem 1.The factors of a probability table are produced by the expert knowledge,which decreases the number of initialized parameters and preserves the completeness of BRB.The reasoning process of BRB-P includes the transformation of input information, calculation of rule activation weights and rule fusion by ER.38In this process, the mapping relationship between fault features (x) and fault modes (D) is presented as follows:

    Table 1 Partial rules of BRB in fault diagnosis model of an airborne missile.

    Fig.1 Interpretable fault diagnosis model based on BRB-P.

    It is necessary to utilize the existing fault samples to improve the model’s performance since the parameters may be subjectively given by experts.The P-CMA-ES optimization algorithm is introduced to enhance the accuracy of the initial BRB-P.In allusion to Problem 2, interpretability constraints are inserted into P-CMA-ES to maintain the interpretability of BRB.12–13Both the deviation from the initial judgment of experts and the over-optimization problem of parameters are avoided.

    The detailed fault diagnosis process is presented in Sections 3 and 4.

    3.Establishment and reasoning of BRB-P model

    In this section, BRB-P is exploited to describe the mapping relationship between fault features and fault modes.The probability table is introduced to establish the initial BRB-P model,and the ER approach is used to fuse the activated rules to generate diagnostic results.

    3.1.Establishment of BRB-P model

    In this subsection, the transition process of expert knowledge to the initial BRB-P is conducted by a probability table.The factors of the probability table are transformed into belief degrees of the initial BRB-P.

    (1) Generation of probability table

    The expert knowledge is aggregated into a probability table.The factors of the probability table are the values of probability.The implication of each factor is the probability of fault feature near a specific referential value when a fault mode occurs.For instance, the probability of Pull-in Voltage(PV) near referential grade L is 0.1 when fault contact wear(D1)occurs in the fault diagnosis of a relay.The specific meaning can be expressed as

    Table 2 Probability table.

    The number of the initialized parameters in BRB-P can be calculated by

    A simple example is cited to illustrate the effect of reducing the number of initialized parameters given by experts.It is assumed that there are 6 fault modes and 3 fault features.Each fault feature has 3 referential grades.There are 192(6×33+3+33=192) parameters set by experts in the traditional BRB.In comparison, experts only specify 84(6×(3+3+3)+3+33=84) parameters in BRB-P.With the number of fault features increasing,the number of the initialized parameters is pictured in Fig.2.

    (2) Establishment of initial BRB-P

    The target of transforming a probability table to rules of BRB is to calculate the belief degrees of each rule.This process can be cognitively expressed by Eq.(10), where ^βn,krepresents the belief degree of the n th fault mode when the input information is the referential values of fault features in the k th rule.

    The ER approach has the advantage of effectively fusing qualitative knowledge and quantitative data.Consequently,factors of the probability table and the weights of fault features are fused into belief degrees of rules in BRB-P by ER.To briefly present the ER fusion process, the first rule is taken as an example to reflect the process from the probability given by experts to belief degrees of rules.

    Fig.2 Comparison of the number of initialized parameters in BRB and BRB-P.

    where μ is an intermediate variable to calculate belief degrees,and ^βj,1is the belief degree of Djin the first rule.According to the above analysis, the belief degrees of the first rule are obtained.Other rules can be generated in the same way.The form of the initial BRB-P is shown in Table 4.

    3.2.Reasoning process of BRB-P

    In this subsection,the reasoning process of BRB-P is presented based on the initial BRB-P,which includes the transformation of input information, calculation of activation weights, and fusion of activated rules.

    Step 1.Transformation of input information

    In this paper, the utility-based method is adopted to transform the input information.The quantitative input is trans-

    Table 3 Evidence collection.

    Table 4 Rule base of initial BRB-P.

    where βjrepresents the belief degree of the j th fault mode.The principle of ER utilized in the process of evidence fusion is essentially the same as that used in the activated rules fusion.

    4.Optimization of BRB-P with interpretability constraints

    In this section, three interpretability constraints are firstly introduced to maintain the interpretability of BRB-P model in the optimization process.26,30Then, the initial BRB-P fault diagnosis model is optimized by P-CMA-ES with three interpretability constraints.The complete process of fault diagnosis based on BRB-P is illustrated at the end of this section.

    4.1.Interpretability constraints

    The definition and the specific application of three interpretability constraints are illustrated as follows:

    (1) The credibility of expert knowledge is introduced to describe the modeling ability of the initial BRB-P.

    The information in the probability table is not entirely reliable due to the subjectivity of experts.The optimization process is exploited to improve the accuracy of the fault diagnosis model.The credibility of expert knowledge26is utilized to calculate the step size of optimization.The step size of optimization is relatively large when the accuracy of the initial BRB-P model is lower.The search range of the optimization is determined by the step size.The higher the credibility of expert knowledge is, the more accurate the initial BRB-P is,and the smaller the optimization step size is.To represent this relationship, the transformation of Euclidean distance is used to define the credibility of expert knowledge in this paper,which is expressed as

    (2) The penalty factor is proposed to enable the optimized BRB-P as close as the initial BRB-P.

    The parameters of BRB-P should be fine-tuned in optimization to ensure the effective utilization of the initial judgment of experts.16As such,in this paper,the objective function of optimization is improved with the penalty factor.The penalty factor is given by

    (3) The rule-activation factor is utilized to mark the activated rules.

    The unactivated rules in the training process should preserve the initial judgment given by experts due to the limited fault samples.In this paper, the rule-activation factor Wkis used to express the activation states of the k th rule.If the k th rule has been activated at least once by the training samples,Wkis equal to 1.Otherwise, Wkis equal to 0, and the corresponding rule cannot participate in the optimization process.It should be pointed out that the k th rule is activated when ωk>0.001.26

    4.2.Parameter optimization with interpretability constraints

    P-CMA-ES algorithm has been well applied to weaken the subjectivity of experts to enhance the accuracy of BRB.However,the purpose of optimization is to improve the accuracy of the model, and the optimized rules are considerably different from the initial rules.Moreover, the parameters may be over-optimized without constraints.Hence, in this paper, the interpretability constraints are introduced to P-CMA-ES to improve the accuracy of the initial BRB-P under the requirement of interpretability.The steps are shown below.

    where f(Ψ) is the traditional objective function, such as the accuracy rate, mean square error, etc.C(β*,β~*) represents the penalty factor.The symbol of objective function depends on the specific application.

    Step 1.Initialization operation

    The initial parameters of the optimization algorithm are presented in Table 5.The initial values of parameters are given by experts.

    Step 2.Sampling operation

    Ψ0is taken as the center point to generate the initial populations with normal distribution:

    Table 5 Initial parameters of P-CMA-ES.

    where c1and c2are the learning rate.p is the evolution path,whose updating rule is

    Step 6.Termination criterion

    Go back to Step 2 until the maximum evolution generation is achieved.

    Step 7.Output optimization result

    If Wk=1, the k th rule optimized is preserved.Otherwise,the parameters of the k th rule are reserved.

    4.3.Process of fault diagnosis based on BRB-P

    In this subsection,the implementation process of the proposed interpretable fault diagnosis method based on BRB-P is presented in Fig.3.The specific steps are shown as follows:

    Step 1.Determination of fault features and fault modes

    The fault features and fault modes are the input and output variables in the process of fault diagnosis respectively, which can be selected by the mechanism analysis or the feature extraction method.The referential grades of fault features can be determined by expert experience or specific methods,such as Receiver Operating Characteristic (ROC) curve.

    Step 2.Establishment and reasoning of BRB-P fault diagnosis model

    The probability table of the system is given by experts based on expert knowledge and historical experience, which has the unified and meaningful form to express the expert knowledge.The factors of the probability table are collected as evidence to generate the belief degrees of rules through the process of normalization and fusion.Then, the input information is transformed into belief distribution of referential grades of fault features by Eqs.(16)–(18).As profiled in Eqs.(19)–(21), the activation weights of rules are calculated.Activated rules are fused to generate the diagnostic results by Eqs.(22)–(23).

    Step 3.Optimization of parameters under interpretability requirement

    Interpretability constraints are generated based on the practical system.The credibility of expert knowledge is counted based on Eqs.(24)–(25) as the step size of optimization to control the search range; the penalty factor shown as Eq.(26) is employed to construct the new objective function to enable the optimized fault diagnosis model as close as the initial one; the rule-activation factors are calculated based on the training samples to pick out the rules that need to be optimized.The interpretability constraints mentioned above are inserted into the P-CMA-ES optimization algorithm shown as Eqs.(29)–(36).The parameters of the initial BRB-P model are optimized to weaken the expert subjectivity under the interpretability requirement.

    Step 4.Output the diagnostic result

    So far,the fault diagnosis model based on BRB-P has been established.When the observation data of fault features are input into BRB-P, the diagnostic result is available.

    5.Case study

    Relay as a control component plays a significant role in signal transmission,load switching,etc.The performance of a relay is considerably essential for the normal operation of the whole system.Taking the relay (JZX-22F) as an example, the fault diagnosis model based on BRB-P is constructed to monitor whether the relay is in a fault state.

    5.1.Preparatory work

    Some normal relays are divided into six groups.One group is not processed.The other five groups are injected Contact Wear(CW), Adhesion of Normally Open Contact (ANOC), Adhesion of Normally Closed Contact (ANCC), Loose Spring(LS) and Coil Fault (CF).In this subsection, the STS2104A relay test system shown in Fig.4 is selected to test the above 6 groups of relays.

    First,the upper computer,tester,adapter and relay are connected to make up the whole system.Then, the equipment parameters in the upper computer are written according to the electrical parameters of the JZX-22F relay.Finally,the test key is pressed repeatedly to obtain the experimental samples which are stored on the upper computer.1800 experimental samples have been obtained,in which 600 groups are test samples and the rest are training samples.The training samples are shown in Fig.5.

    In BRB-P model,four variables:Pull-in Voltage(PV),Pullin Time(PT),Resistance of No Contact(RNC)and Resistance of Contact(RC)are selected as fault features.Three referential grades are assigned for PV, PT, RNC and RC.They are Low(L),Medium(M)and High(H).The referential values of fault features are profiled in Table 6, which are given by experts according to the manual and experimental samples.266 working states of a relay are shown in Table 7.

    5.2.Establishment of initial BRB-P

    In this subsection, the initial BRB-P is automatically generated by the probability table given by experts.The first task is to set up the probability table according to the expert knowledge.On the premise of known fault modes and fault features, the probability of the fault feature near a specific referential value when a fault mode occurs is given by experts.The probability table of the relay is shown in Table 8.

    The initial BRB-P is established after the process of evidence acquisition, normalization and fusion.The partial rules of the initial BRB-P are shown in Table 9, and the complete rule base is presented in Appendix A.

    Fig.3 Process of fault diagnosis based on BRB-P.

    Fig.4 Relay test system.

    5.3.Optimization process

    The iterative optimization process is shown in Fig.6.

    The partial optimized rules are shown in Table 10, and the complete rule base of BRB-P is given in Appendix B.

    Fig.5 Training samples of relays.

    Table 6 Referential value of fault feature.

    Table 7 Working states of a relay.

    Table 8 Probability table of relay fault diagnosis.

    Table 9 Partial rules of initial BRB-P.

    Fig.6 Iterative process of optimization.

    5.4.Comparative study

    In this subsection, the traditional BRB, the BRB-P without interpretability constraints, the back propagation neural network and support vector machine are conducted for comparison.

    (1) Traditional BRB

    According to the expert knowledge and historical experience, the belief degrees, the fault feature weights and the rule weights can be set by experts.In contrast with the proposed method, the initial BRB is established through the traditional method.In other words,81 rules are provided by experts based on the mechanism in the process of relay fault diagnosis.Please refer to Appendix C for the specific BRB, named BRB0.

    (2) BRB-P without interpretability constraints

    The initial BRB-P is optimized by the P-CMA-ES without interpretability constraints to verify the effectiveness of interpretability constraints by comparison.The specific parameters of this optimization include: the search step size is 0.5; the maximum evolution generation is 200; the objective function is f(Ψ).In this optimization, the parameters of all rules are optimized based on the initial BRB-P.Please refer to Appendix D for the result of optimization, named BRB1.

    (3) Back Propagation Neural Network (BPNN) method

    Table 10 Partial rules of BRB-P.

    Table 11 Parameters of BPNN.

    Fig.7 Structure of BPNN.

    Table 12 Accuracy comparison of different models.

    BPNN is used to establish the relationship of fault features and fault modes, which is a typical black-box algorithm.The structure of BPNN is shown in Fig.7, where fi(i=1,2,3) is the transfer function.The structure of BPNN includes the input layer, hidden layer and output layer, which are connected through weight parameters.Similarly, the values of fault features are regarded as the input of the neural network model.However, the output is one of the working states of the relay represented by a digital label.In the training process,the weight parameters are constantly updated to obtain better diagnostic results.The parameters of BPNN are shown in Table 11.

    (4) Support Vector Machine (SVM) method

    SVM is a typical method used in discriminative model,which has the advantages of less dependence on samples and preventing overfitting.The main target of SVM is to search the classification hyperplane to classify samples.This experiment is carried out with the help of the LibSVM toolbox.Radial Basis Function (RBF) is adopted as kernel function.

    5.5.Result analysis

    In this subsection,the comparative study of the initial BRB-P,BRB0,BRB-P,BRB1,BPNN and SVM is made to analyze the accuracy and interpretability of models.

    (1) Analysis of model accuracy and the number of initialized parameters

    The accuracy of the models mentioned above is concluded in Table 12, which is divided into two parts.

    In Part A, the accuracy of the initial BRB-P is desirable.Less subjectivity of experts is contained in the initial BRB-P owing to fewer parameters given by experts.The accuracy of the initial BRB-P is higher than that of BRB0.The results of the relay fault diagnosis in Part A are shown in Fig.8.By comparing BRB0 with BRB-P, the number of initialized parameters reduces from 571 to 157, calculated by Eqs.(8)–(9).It can be seen that the number of parameters in BRB-P is reduced, and the accuracy of BRB-P is effectively enhanced by introducing the probability table.

    The diagnostic results of Part B are shown in Fig.9.Although BPNN and SVM can achieve high accuracy, they are data-driven models whose reasoning processes cannot be understood intuitively.It is obvious that the modeling ability of SVM is better than that of BPNN when there are limited fault samples.Initializing the parameters by experts releases the burden of optimization,and the accuracy of the optimized models based on BRB-P is higher in these models.Compared with BRB1, the accuracy of BRB-P is decreased due to the constrained optimization process.However,more importantly,the interpretability of the model is retained under the premise of the same maximum evolution generation.The specific analysis of interpretability is presented below.

    Fig.8 Comparison of model accuracy in Part A.

    Fig.9 Comparison of model accuracy in Part B.

    Table 13 Euclidean distance of belief degrees for different models.

    (2) Analysis of model interpretability

    BRB is essentially an interpretable model.However, the excellent property may be destroyed in the optimization process, which mainly reflects in two aspects: the deviation from the initial judgment of experts and over-optimization of parameters.The following will focus on these two aspects.

    To avoid deviating from experts’initial judgment,the optimized parameters are required to be close to the initial BRB-P.As can be seen from Tables 9 and 10, the parameters of rules are fine-tuned based on the initial rules.In other words,although the parameters are optimized to improve the accuracy, most of the rules are still in line with the cognition of experts.More specifically, the Euclidean distance shown as Eq.(37)is introduced to describe this closeness.From Table 13,compared with BRB1, the Euclidean distance between the initial BRB-P and BRB-P is shorter, which indicates that the result of optimization with interpretability constraints is closer to the initial judgment of experts.

    Fig.10 Belief degrees of rules.

    To further reveal the level of closeness,the belief degrees of models and the difference between them are drawn in Figs.10 and 11 respectively.It can be seen from Fig.10 that the belief degrees of BRB-P are closer to the belief degrees of the initial BRB-P than those of the BRB1.Fig.11 also presents that the differences between the initial belief degrees and the optimized belief degrees of BRB-P are closer to 0.

    To avoid the over-optimization of parameters, the ruleactivation factor is proposed.The value of the ruleactivation factor is shown in Fig.12.The sum of Wkis 76,which means that the remaining 5 rules are not involved in the fusion process.

    The unactivated rules are usually rare conditions or bad inputs.To ensure the completeness of the rule base, the rules in all cases are set by experts to deal with new situations or error signals.Here, 5 unactivated rules are presented in Table 14, in which working states of relays are recognized as normal states by experts.The rule-activation factor is introduced into the optimization process to judge whether a rule is activated, so as to determine whether to optimize the rule.In the fault diagnosis model based on BRB-P, the parameters of the 5 rules mentioned above are not optimized and reserve the initial judgment given by experts to avoid the overoptimization of parameters.However, these rules are updated in BRB1 without basis.Obviously, the optimized results for these 5 rules are irrational in BRB1.For example,the diagnostic results of BRB-P and the initial BRB-P are that the relay is in normal state when the input is consistent with the antecedent of the 27th rule.However, the diagnostic result of BRB1 is the adhesion of normally closed contact in the same case.Consequently, the over-optimization of parameters is avoided, and the interpretability of fault diagnosis model is preserved in this proposed method.

    Fig.11 Difference of belief degrees.

    Fig.12 Rule-activation factor.

    Table 14 Unactivated rules.

    Table A1 Weights of fault features in initial BRB-P.

    Table A2 Rule base of initial BRB-P.

    Table B1 Probability table of BRB-P.

    Table B2 Weights of fault feature in BRB-P.

    6.Conclusions

    Fault diagnosis plays an irreplaceable role in the operation and maintenance of the system.An interpretable fault diagnosis approach can build good trust between the model and human beings.As an interpretable expert system, BRB has been widely applied in the field of fault diagnosis.

    In this paper, a new interpretable fault diagnosis model based on BRB-P is proposed, where the probability table is exploited to transform the expert knowledge into rules of BRB-P.Three interpretability constraints are inserted into the P-CMA-ES optimization algorithm to improve the accuracy of the model under the requirement of interpretability.Moreover, the fault diagnosis model of the JZX-22F relay based on BRB-P is established.The accuracy and interpretability of BRB-P model can be maintained by comparing with other fault diagnosis models.

    There are mainly two contributions in this paper.On the one hand, the probability table is proposed to express the expert knowledge and establish the initial BRB-P, which reduces the number of initialized parameters under the premise of higher accuracy.The problem of embedding expert knowledge to BRB when many fault features have multiple referential grades is solved.On the other hand, three interpretability constraints are introduced to P-CMA-ES.The credibility of expert knowledge can describe the modeling ability of the initial BRB-P.The penalty factor can enable the optimized BRBP as close as the initial BRB-P.The problems of the deviation from the experts’initial judgment can be solved by these twointerpretability constraints.The rule-activation factor is utilized to mark the activated rules during the training process,which avoids over-optimization of parameters.The interpretability of BRB-P is maintained in optimization.

    Table D1 Weights of fault features in BRB1.

    In this paper, a new interpretable fault diagnosis model based on expert knowledge is proposed.Many data-based models are used in the industry when mass fault samples are available.How to establish an interpretable fault diagnosis model in terms of mass fault samples and how to conduct the performance analysis of BRB-P should be further discussed.The real-time data can be disturbed by the actual environment and the reliability of sensors in practical applications.One of the future works is how to process data to improve the accuracy of fault diagnosis.

    Table C2 Rule base of BRB0.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (No.61833016), the Shaanxi Outstanding Youth Science Foundation, China (No.2020JC-34), the Shaanxi Science and Technology Innovation Team, China(No.2022TD-24),and the Natural Science Foundation of Heilongjiang Province of China (No.LH2021F038).

    Appendix A.

    See Table A1 and Table A2.

    Appendix B.

    See Table B1, Table B2 and Table B3.

    Appendix C.

    See Table C1 and Table C2.

    Appendix D.

    See Table D1 and Table D2.

    Table D2 Rule base of BRB1.

    亚洲国产毛片av蜜桃av| 日韩大片免费观看网站| 十八禁网站网址无遮挡| 极品少妇高潮喷水抽搐| 色播在线永久视频| 免费不卡的大黄色大毛片视频在线观看| 日韩av不卡免费在线播放| 黑丝袜美女国产一区| 亚洲男人天堂网一区| 如何舔出高潮| av在线播放精品| 日韩不卡一区二区三区视频在线| 日本免费在线观看一区| 高清在线视频一区二区三区| 欧美精品亚洲一区二区| 黄色怎么调成土黄色| 99久久中文字幕三级久久日本| 成人手机av| 中文字幕人妻丝袜制服| 不卡视频在线观看欧美| 亚洲第一av免费看| 男人舔女人的私密视频| 国产亚洲一区二区精品| 丝袜脚勾引网站| 国语对白做爰xxxⅹ性视频网站| 多毛熟女@视频| 精品久久蜜臀av无| 亚洲天堂av无毛| 视频在线观看一区二区三区| 日本免费在线观看一区| 欧美日本中文国产一区发布| 久久毛片免费看一区二区三区| 久久精品国产鲁丝片午夜精品| 黄色 视频免费看| 水蜜桃什么品种好| 国产精品久久久久成人av| 看免费av毛片| 青春草亚洲视频在线观看| 99久久中文字幕三级久久日本| 国产精品偷伦视频观看了| 高清黄色对白视频在线免费看| 九草在线视频观看| 尾随美女入室| 一区福利在线观看| 丝袜喷水一区| 日韩中文字幕欧美一区二区 | 叶爱在线成人免费视频播放| 精品国产一区二区三区久久久樱花| 久久女婷五月综合色啪小说| 国产精品不卡视频一区二区| 亚洲三级黄色毛片| 午夜福利视频在线观看免费| 中文字幕亚洲精品专区| 最近的中文字幕免费完整| www日本在线高清视频| 男男h啪啪无遮挡| xxxhd国产人妻xxx| 少妇猛男粗大的猛烈进出视频| 国产成人aa在线观看| 久久久久久久大尺度免费视频| 亚洲国产看品久久| 欧美黄色片欧美黄色片| 在线天堂中文资源库| 精品国产乱码久久久久久男人| 久久人妻熟女aⅴ| 亚洲av中文av极速乱| 精品一区二区免费观看| 久久久国产一区二区| 亚洲情色 制服丝袜| 又黄又粗又硬又大视频| 国产精品欧美亚洲77777| 多毛熟女@视频| 亚洲,一卡二卡三卡| 日本vs欧美在线观看视频| 99re6热这里在线精品视频| av视频免费观看在线观看| 久久精品国产亚洲av天美| 热re99久久精品国产66热6| 国产片特级美女逼逼视频| 嫩草影院入口| 美女xxoo啪啪120秒动态图| 免费不卡的大黄色大毛片视频在线观看| 国产一区二区三区av在线| 男女啪啪激烈高潮av片| 免费观看在线日韩| 久久鲁丝午夜福利片| 精品一区二区三卡| 我要看黄色一级片免费的| 黄片播放在线免费| 久久99热这里只频精品6学生| 亚洲av欧美aⅴ国产| 91午夜精品亚洲一区二区三区| 午夜免费男女啪啪视频观看| 精品亚洲乱码少妇综合久久| 妹子高潮喷水视频| 亚洲五月色婷婷综合| 丰满少妇做爰视频| 久久97久久精品| 日韩av免费高清视频| 久久这里只有精品19| 国产熟女午夜一区二区三区| 一级片免费观看大全| 免费观看在线日韩| 久久国产精品大桥未久av| 深夜精品福利| 免费看不卡的av| 巨乳人妻的诱惑在线观看| 十八禁高潮呻吟视频| 老汉色∧v一级毛片| 九色亚洲精品在线播放| 亚洲欧美清纯卡通| 少妇猛男粗大的猛烈进出视频| 久久综合国产亚洲精品| 午夜激情久久久久久久| 99热国产这里只有精品6| 大香蕉久久成人网| 边亲边吃奶的免费视频| 久久久久久久久久人人人人人人| 不卡视频在线观看欧美| 亚洲视频免费观看视频| 精品一区二区三卡| 97在线人人人人妻| www.自偷自拍.com| 国产一级毛片在线| 亚洲精品在线美女| 日韩中文字幕欧美一区二区 | 欧美激情 高清一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲国产毛片av蜜桃av| 亚洲少妇的诱惑av| 免费黄网站久久成人精品| 飞空精品影院首页| 麻豆乱淫一区二区| 9热在线视频观看99| 免费在线观看视频国产中文字幕亚洲 | 久久国产精品大桥未久av| 老汉色av国产亚洲站长工具| 满18在线观看网站| 在线天堂最新版资源| 大陆偷拍与自拍| 免费大片黄手机在线观看| 亚洲精品一二三| 日韩欧美一区视频在线观看| 国产精品久久久av美女十八| 老女人水多毛片| 新久久久久国产一级毛片| 午夜精品国产一区二区电影| 另类精品久久| 久久久久久人人人人人| 亚洲av综合色区一区| 亚洲av综合色区一区| 国产一级毛片在线| 亚洲av综合色区一区| 伦精品一区二区三区| 男人爽女人下面视频在线观看| 黄色毛片三级朝国网站| 日韩成人av中文字幕在线观看| 欧美日韩视频高清一区二区三区二| 亚洲av电影在线观看一区二区三区| 欧美少妇被猛烈插入视频| 亚洲精品国产一区二区精华液| 久久久国产一区二区| 亚洲国产欧美在线一区| 熟女少妇亚洲综合色aaa.| 26uuu在线亚洲综合色| 精品国产一区二区三区久久久樱花| 精品少妇黑人巨大在线播放| 在线观看免费视频网站a站| 日韩视频在线欧美| 日韩精品免费视频一区二区三区| 在线观看一区二区三区激情| 国产野战对白在线观看| 男女边吃奶边做爰视频| av免费在线看不卡| 伦理电影免费视频| 26uuu在线亚洲综合色| 99香蕉大伊视频| 26uuu在线亚洲综合色| 亚洲欧洲国产日韩| 99热国产这里只有精品6| 少妇人妻 视频| 国产精品秋霞免费鲁丝片| av福利片在线| 日韩制服骚丝袜av| 18禁裸乳无遮挡动漫免费视频| 日韩一卡2卡3卡4卡2021年| 国产精品秋霞免费鲁丝片| 亚洲综合色惰| 国产男人的电影天堂91| 午夜久久久在线观看| 91久久精品国产一区二区三区| 精品一区在线观看国产| 国产精品一区二区在线不卡| 我的亚洲天堂| 日本爱情动作片www.在线观看| 中文字幕av电影在线播放| 毛片一级片免费看久久久久| 久久精品国产综合久久久| 亚洲欧美清纯卡通| 亚洲精品日本国产第一区| 亚洲五月色婷婷综合| 欧美成人精品欧美一级黄| 18禁观看日本| 男女午夜视频在线观看| 日韩视频在线欧美| 另类亚洲欧美激情| 欧美+日韩+精品| 99国产精品免费福利视频| 蜜桃在线观看..| 国产淫语在线视频| 1024香蕉在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲五月色婷婷综合| 久久人人爽人人片av| 婷婷色av中文字幕| 黑人巨大精品欧美一区二区蜜桃| 美女中出高潮动态图| 亚洲视频免费观看视频| 九草在线视频观看| 中文字幕另类日韩欧美亚洲嫩草| 精品福利永久在线观看| 久久精品亚洲av国产电影网| 久久久久精品久久久久真实原创| 色吧在线观看| 国产视频首页在线观看| 国精品久久久久久国模美| 人妻人人澡人人爽人人| 午夜激情av网站| 美女视频免费永久观看网站| videos熟女内射| 午夜免费观看性视频| 十八禁网站网址无遮挡| 国产熟女欧美一区二区| 天天躁夜夜躁狠狠久久av| 男女国产视频网站| 超碰成人久久| 精品亚洲成a人片在线观看| 不卡视频在线观看欧美| 亚洲精品自拍成人| 少妇人妻久久综合中文| 欧美日韩视频高清一区二区三区二| 亚洲精品第二区| av在线app专区| 日韩一区二区视频免费看| 国产在线视频一区二区| 国产福利在线免费观看视频| 男人操女人黄网站| 有码 亚洲区| 午夜日韩欧美国产| 婷婷色av中文字幕| 久久精品aⅴ一区二区三区四区 | 精品人妻熟女毛片av久久网站| 在线看a的网站| 欧美人与性动交α欧美精品济南到 | kizo精华| 香蕉丝袜av| 亚洲欧美成人综合另类久久久| 中文字幕色久视频| 9色porny在线观看| 大片免费播放器 马上看| 欧美精品亚洲一区二区| 亚洲国产精品国产精品| 赤兔流量卡办理| 久久鲁丝午夜福利片| 免费观看av网站的网址| 亚洲色图 男人天堂 中文字幕| 在线观看三级黄色| √禁漫天堂资源中文www| 国产亚洲欧美精品永久| av在线播放精品| 国产黄色视频一区二区在线观看| 亚洲一区中文字幕在线| 亚洲av免费高清在线观看| 99九九在线精品视频| 我要看黄色一级片免费的| 精品国产超薄肉色丝袜足j| 自线自在国产av| 美女中出高潮动态图| 精品酒店卫生间| 国产黄色免费在线视频| 爱豆传媒免费全集在线观看| 熟女少妇亚洲综合色aaa.| 欧美日韩视频高清一区二区三区二| 久久热在线av| 亚洲婷婷狠狠爱综合网| 亚洲精品成人av观看孕妇| 亚洲av在线观看美女高潮| 秋霞伦理黄片| 黄色一级大片看看| 亚洲精品aⅴ在线观看| 国产国语露脸激情在线看| 国产亚洲最大av| 我要看黄色一级片免费的| 久久久久久人人人人人| 秋霞伦理黄片| 国产一区二区激情短视频 | 精品人妻熟女毛片av久久网站| 最近最新中文字幕大全免费视频 | av不卡在线播放| 男人舔女人的私密视频| 不卡视频在线观看欧美| 欧美亚洲 丝袜 人妻 在线| 免费久久久久久久精品成人欧美视频| 欧美日韩国产mv在线观看视频| 三上悠亚av全集在线观看| 日韩在线高清观看一区二区三区| 久久99一区二区三区| 国产欧美日韩综合在线一区二区| 久久影院123| a 毛片基地| 天天影视国产精品| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩视频精品一区| 欧美中文综合在线视频| 热99国产精品久久久久久7| 国产欧美日韩一区二区三区在线| 两个人免费观看高清视频| 欧美xxⅹ黑人| 2018国产大陆天天弄谢| 9热在线视频观看99| 欧美黄色片欧美黄色片| 久久精品国产亚洲av高清一级| 视频区图区小说| av天堂久久9| 国产成人aa在线观看| 国产精品久久久av美女十八| 国产精品av久久久久免费| av福利片在线| 久久ye,这里只有精品| 免费在线观看视频国产中文字幕亚洲 | 欧美亚洲 丝袜 人妻 在线| 日韩三级伦理在线观看| 免费观看性生交大片5| 少妇人妻精品综合一区二区| 男人舔女人的私密视频| 99国产精品免费福利视频| 国产精品国产三级专区第一集| 国产精品人妻久久久影院| 午夜91福利影院| 精品99又大又爽又粗少妇毛片| 尾随美女入室| 亚洲av中文av极速乱| 亚洲精品国产一区二区精华液| 精品国产乱码久久久久久男人| 一级毛片电影观看| 99久久综合免费| 老熟女久久久| 丝袜美腿诱惑在线| 中国三级夫妇交换| 又黄又粗又硬又大视频| 欧美成人午夜精品| 欧美97在线视频| 国产成人精品久久久久久| 欧美日韩亚洲国产一区二区在线观看 | 涩涩av久久男人的天堂| 国产免费现黄频在线看| 如日韩欧美国产精品一区二区三区| 亚洲av电影在线观看一区二区三区| 97人妻天天添夜夜摸| 国产av一区二区精品久久| 中文天堂在线官网| 欧美国产精品一级二级三级| 国产片特级美女逼逼视频| 91成人精品电影| 最近最新中文字幕大全免费视频 | 亚洲精品国产av蜜桃| 免费在线观看完整版高清| 超色免费av| 日本av手机在线免费观看| 国产一级毛片在线| 久久鲁丝午夜福利片| 国产 一区精品| 综合色丁香网| 中文字幕亚洲精品专区| 看十八女毛片水多多多| 亚洲国产av影院在线观看| 在线 av 中文字幕| 亚洲精品一二三| 777米奇影视久久| 晚上一个人看的免费电影| 国产男女超爽视频在线观看| 亚洲av福利一区| 亚洲欧美日韩另类电影网站| 国产在线免费精品| 2018国产大陆天天弄谢| 人妻系列 视频| a级片在线免费高清观看视频| 欧美bdsm另类| 啦啦啦视频在线资源免费观看| 国产成人欧美| 国产探花极品一区二区| 国产熟女午夜一区二区三区| 国产精品久久久久久av不卡| 亚洲国产看品久久| 久久免费观看电影| 搡老乐熟女国产| 一边摸一边做爽爽视频免费| 日韩在线高清观看一区二区三区| 只有这里有精品99| 天美传媒精品一区二区| 国产精品久久久久久久久免| 国产一区二区激情短视频 | 老鸭窝网址在线观看| 国产精品三级大全| 99热全是精品| 亚洲,欧美,日韩| 精品国产乱码久久久久久小说| 黄色毛片三级朝国网站| 自线自在国产av| 国产av精品麻豆| 亚洲精品乱久久久久久| 国产精品嫩草影院av在线观看| 男女无遮挡免费网站观看| 老鸭窝网址在线观看| 成人18禁高潮啪啪吃奶动态图| 久久精品久久久久久久性| 美女高潮到喷水免费观看| 日本vs欧美在线观看视频| 欧美人与性动交α欧美软件| 日韩视频在线欧美| 麻豆av在线久日| 黑人欧美特级aaaaaa片| 国产成人欧美| 波野结衣二区三区在线| 少妇的丰满在线观看| 男人爽女人下面视频在线观看| 午夜福利一区二区在线看| 免费观看a级毛片全部| 亚洲 欧美一区二区三区| 亚洲精品国产av蜜桃| 成人亚洲欧美一区二区av| 可以免费在线观看a视频的电影网站 | 成人亚洲精品一区在线观看| 在线观看国产h片| 国产又爽黄色视频| 美女主播在线视频| 99久久精品国产国产毛片| 国产淫语在线视频| 久久精品亚洲av国产电影网| 亚洲 欧美一区二区三区| 高清在线视频一区二区三区| 午夜福利在线观看免费完整高清在| 日本猛色少妇xxxxx猛交久久| 亚洲三区欧美一区| 国产白丝娇喘喷水9色精品| 一级毛片 在线播放| 国产成人aa在线观看| 少妇熟女欧美另类| 欧美黄色片欧美黄色片| 亚洲欧美成人综合另类久久久| 婷婷色综合www| 少妇的逼水好多| 欧美日韩一级在线毛片| 99九九在线精品视频| 久久精品熟女亚洲av麻豆精品| 街头女战士在线观看网站| 在线观看免费高清a一片| 人人澡人人妻人| 纯流量卡能插随身wifi吗| 一级毛片黄色毛片免费观看视频| 久久久国产一区二区| 久久久精品区二区三区| 成人亚洲欧美一区二区av| 亚洲精品国产色婷婷电影| h视频一区二区三区| 又大又黄又爽视频免费| 亚洲欧美中文字幕日韩二区| 捣出白浆h1v1| 多毛熟女@视频| 亚洲美女搞黄在线观看| 欧美精品人与动牲交sv欧美| 日日摸夜夜添夜夜爱| 精品人妻在线不人妻| 国产成人91sexporn| 亚洲一区中文字幕在线| 亚洲在久久综合| 啦啦啦在线观看免费高清www| 久久久久久人妻| 咕卡用的链子| 蜜桃在线观看..| 免费不卡的大黄色大毛片视频在线观看| 黄网站色视频无遮挡免费观看| 大陆偷拍与自拍| 人人澡人人妻人| 免费少妇av软件| 欧美日韩精品成人综合77777| 国产老妇伦熟女老妇高清| 日韩熟女老妇一区二区性免费视频| 男女国产视频网站| 丝袜喷水一区| 少妇被粗大猛烈的视频| 亚洲精品成人av观看孕妇| 99国产综合亚洲精品| 亚洲情色 制服丝袜| 日本欧美视频一区| 精品国产乱码久久久久久小说| 亚洲精品日本国产第一区| 国产成人精品福利久久| 国产精品欧美亚洲77777| 在线亚洲精品国产二区图片欧美| 中文天堂在线官网| 啦啦啦中文免费视频观看日本| 人妻 亚洲 视频| 最新中文字幕久久久久| 黄网站色视频无遮挡免费观看| 免费观看av网站的网址| 亚洲成人手机| 亚洲成av片中文字幕在线观看 | 18禁国产床啪视频网站| 久久久亚洲精品成人影院| 啦啦啦在线观看免费高清www| 午夜老司机福利剧场| 乱人伦中国视频| 久久久国产精品麻豆| 搡女人真爽免费视频火全软件| 国产人伦9x9x在线观看 | videos熟女内射| 亚洲国产精品一区二区三区在线| 夫妻午夜视频| 三级国产精品片| 国产一区二区三区综合在线观看| 日韩精品有码人妻一区| 99九九在线精品视频| 日韩精品免费视频一区二区三区| 满18在线观看网站| 日韩大片免费观看网站| 丰满少妇做爰视频| 美女午夜性视频免费| 街头女战士在线观看网站| 免费观看性生交大片5| 伊人久久国产一区二区| 亚洲av日韩在线播放| 久久热在线av| h视频一区二区三区| 一本色道久久久久久精品综合| 熟妇人妻不卡中文字幕| 午夜av观看不卡| 免费看av在线观看网站| 精品福利永久在线观看| 高清欧美精品videossex| 人体艺术视频欧美日本| 久久精品国产综合久久久| 亚洲精品美女久久久久99蜜臀 | 伊人久久大香线蕉亚洲五| 日产精品乱码卡一卡2卡三| 99久国产av精品国产电影| 久久女婷五月综合色啪小说| 综合色丁香网| 国产又爽黄色视频| 十分钟在线观看高清视频www| 国产精品久久久久久久久免| 免费高清在线观看视频在线观看| 欧美精品亚洲一区二区| 欧美日韩一区二区视频在线观看视频在线| 久久青草综合色| 亚洲男人天堂网一区| 伦理电影免费视频| 两个人免费观看高清视频| 高清av免费在线| 国产xxxxx性猛交| 各种免费的搞黄视频| 国产日韩欧美亚洲二区| 中文字幕人妻熟女乱码| 韩国av在线不卡| 国产精品偷伦视频观看了| 激情视频va一区二区三区| 久久久亚洲精品成人影院| 咕卡用的链子| 亚洲欧美精品自产自拍| 黄色视频在线播放观看不卡| 国产成人精品无人区| 日韩 亚洲 欧美在线| 老司机影院毛片| 免费观看性生交大片5| 熟女av电影| 国产人伦9x9x在线观看 | 国产精品一国产av| 一边摸一边做爽爽视频免费| 中文天堂在线官网| 曰老女人黄片| 国产成人精品无人区| 99久国产av精品国产电影| 精品国产一区二区三区久久久樱花| 热re99久久国产66热| 亚洲中文av在线| 亚洲,一卡二卡三卡| 欧美日本中文国产一区发布| 在线观看美女被高潮喷水网站| 免费看不卡的av| 男女国产视频网站| 精品人妻在线不人妻| 免费久久久久久久精品成人欧美视频| 欧美av亚洲av综合av国产av | 超色免费av| 99热全是精品| 亚洲av中文av极速乱| 最近2019中文字幕mv第一页| 国产女主播在线喷水免费视频网站| 热99国产精品久久久久久7| 日本免费在线观看一区| av电影中文网址| 国产av码专区亚洲av| 日韩免费高清中文字幕av| 国产av精品麻豆| 最近手机中文字幕大全| 在线天堂最新版资源| 日韩免费高清中文字幕av| 国产精品av久久久久免费| 亚洲国产欧美网| 人体艺术视频欧美日本| 人妻人人澡人人爽人人| 熟女电影av网| 2021少妇久久久久久久久久久| 欧美激情高清一区二区三区 | 啦啦啦在线免费观看视频4|