CN109116219B - Distributed fault diagnosis method for circuit system - Google Patents

Distributed fault diagnosis method for circuit system Download PDF

Info

Publication number
CN109116219B
CN109116219B CN201811082945.7A CN201811082945A CN109116219B CN 109116219 B CN109116219 B CN 109116219B CN 201811082945 A CN201811082945 A CN 201811082945A CN 109116219 B CN109116219 B CN 109116219B
Authority
CN
China
Prior art keywords
sensor
distributed
subsystem
circuit system
minimum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811082945.7A
Other languages
Chinese (zh)
Other versions
CN109116219A (en
Inventor
郁明
肖晨雨
姜苍华
王海
杨双龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201811082945.7A priority Critical patent/CN109116219B/en
Publication of CN109116219A publication Critical patent/CN109116219A/en
Application granted granted Critical
Publication of CN109116219B publication Critical patent/CN109116219B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/28Testing of electronic circuits, e.g. by signal tracer
    • G01R31/2832Specific tests of electronic circuits not provided for elsewhere
    • G01R31/2836Fault-finding or characterising
    • G01R31/2843In-circuit-testing

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

The invention discloses a distributed fault diagnosis method of a circuit system, which comprises the following steps: modeling the circuit system, and configuring the causal relationship among all parameters to obtain a global diagnosis bonding diagram model of the circuit system; taking each sensor in the circuit system as a basic unit, taking the measured value of the sensor as the local output of the circuit system, extracting the minimum subsystem based on each sensor from the circuit system, and respectively obtaining the local diagnosis bonding diagram model of each minimum subsystem; obtaining a distributed analytic redundancy relation according to a local diagnostic bond map model of a minimum subsystem, and analyzing the distributed analytic redundancy relation to obtain a distributed fault feature matrix; and carrying out fault diagnosis according to the distributed fault feature matrix. The invention effectively improves the fault isolation performance under the conditions of single fault and multiple faults, reduces the complexity of subsequent fault element identification, and improves the diagnosis speed and the diagnosis precision of the faults.

Description

Distributed fault diagnosis method for circuit system
Technical Field
The invention relates to the field of fault diagnosis of a circuit system, in particular to a distributed fault diagnosis method of the circuit system.
Background
In modern industrial production, many complex electrical and mechanical circuit systems are generated, such as complex integrated circuit systems, electromechanical circuit systems, hydraulic circuit systems, and the importance of fault diagnosis and real-time monitoring of the circuit systems is increasingly highlighted in order to improve the safety and reliability of the circuit systems.
With the increase of the complexity of the circuit system, the difficulty of fault diagnosis of the circuit system is increased, and a common fault diagnosis method is to model the circuit system based on a bonding diagram to obtain a global diagnosis bonding diagram model of the circuit system, directly obtain a global fault feature matrix of the circuit system according to the global diagnosis bonding diagram model of the circuit system, and then perform fault diagnosis on the circuit system according to the global fault feature matrix. However, as the complexity of the circuit system increases, the following main problems occur in this method of global fault diagnosis: firstly, under the condition of single fault, because some elements in a circuit system have the same fault characteristics, fault isolation cannot be realized; second, under multiple fault conditions, fault isolation cannot be achieved because the fault characteristics of some components in the circuitry cancel or mask each other. Meanwhile, the complexity of subsequent fault element identification is affected by the performance of fault isolation, and the more elements are included in a possible fault set obtained by fault isolation, the higher the complexity of fault element identification is, and the lower the identification accuracy is.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention provides a distributed fault diagnosis method for a circuit system, which effectively improves the fault isolation performance under the conditions of single fault and multiple faults, reduces the complexity of subsequent fault element identification, and improves the fault diagnosis speed and the fault diagnosis precision.
In order to achieve the purpose, the invention adopts the following technical scheme that:
a distributed fault diagnosis method of a circuit system, comprising the steps of:
s1, modeling the circuit system, and configuring the causal relationship among all parameters to obtain a global diagnosis bonding diagram model of the circuit system;
s2, taking each sensor in the circuit system as a basic unit, taking the measured value of the sensor as the local output of the circuit system, extracting the minimum subsystem based on each sensor from the circuit system, and respectively obtaining the local diagnosis bond map model of each minimum subsystem;
s3, obtaining a distributed analytic redundancy relationship according to the local diagnostic bond map model of the minimum subsystem, and analyzing the distributed analytic redundancy relationship to obtain a distributed fault feature matrix;
and S4, performing fault diagnosis according to the distributed fault feature matrix.
In step S2, the manner of obtaining the local diagnostic bond map model of the smallest subsystem is specifically as follows: taking a certain sensor of the circuit system as an output sensor, taking the measured value of the output sensor as the local output of the circuit system, namely taking the measured value of the output sensor as the output signal of the minimum subsystem based on the output sensor, starting from the output sensor to the end of another sensor, taking the measured value of the sensor which is the end sensor as the input signal of the minimum subsystem, and obtaining a local diagnostic bonding map model of the minimum subsystem based on the output sensor; in a path from the output sensor to the termination sensor, a parameter associated with the measured value of the output sensor is derived from causal relationships between various parameters on a global diagnostic bond map model of the circuitry, and the parameter associated with the measured value of the output sensor is used as a parameter of a local diagnostic bond map model based on a smallest subsystem of the output sensor.
In step S1, the parameters of the global diagnostic bond map model of the circuit system include: potential source UM=VinResistive element { R1,R2,R3}, capacitive element { C1,C2,C3}, voltage sensor { De1,De2}, current sensor Df, switch k1(ii) a Wherein, the potential source VinIs the input voltage of the system; resistance R1And a capacitor C1Series connected, resistance R2And a capacitor C2Series connected, resistance R3And a capacitor C3Are connected in series; will resistance R1And a capacitor C1Formed series circuit, resistor R2And a capacitor C2Formed series circuit, resistor R3And a capacitor C3The formed series circuits are connected in parallel; voltage sensor De1Connected in parallel to a capacitor C2Both ends of (a); voltage sensor De2Connected in parallel to a capacitor C3Both ends of (a); current sensor Df and resistor R1Potential source VinAre connected in series; switch k1And a resistor R3Are connected in series.
Using current sensor Df as basic unit from current sensor Df to voltage sensor De1Terminate, extract based on current sensor DfMinimum subsystem
Figure BDA0001802428360000021
The parameters of the local diagnostic bond map model of the minimal subsystem include: potential source
Figure BDA0001802428360000022
Resistive element { R1,R2}, capacitive element C1A flow sensor Df;
with a voltage sensor De1As a basic unit, a slave voltage sensor De1To start, to the voltage sensor De2And the current sensor Df is terminated, and the voltage sensor De is extracted1Minimum subsystem of
Figure BDA0001802428360000031
The minimum subsystem
Figure BDA0001802428360000032
The parameters of the local diagnostic bond map model of (a) include: potential source
Figure BDA0001802428360000033
Resistive element { R2,R3}, capacitive element { C1,C2}, a voltage sensor De1Switch k1
With a voltage sensor De2As a basic unit, a slave voltage sensor De2To start, to the voltage sensor De1Terminate, extract based on the voltage sensor De2Minimum subsystem of
Figure BDA0001802428360000034
The minimum subsystem
Figure BDA0001802428360000035
The parameters of the local diagnostic bond map model of (a) include: potential source
Figure BDA0001802428360000036
Resistive element R3Capacitive element C3D, a voltage sensor De2Switch k1
In step S3, the method includes the following steps:
s31, obtaining a distributed analytical redundancy relationship according to the local diagnostic bond map models of the three minimum subsystems and the causal relationship among all parameters in the local diagnostic bond map models;
in the above manner, the minimum subsystem is obtained
Figure BDA0001802428360000037
Distributed analytic redundancy relationship DARR1Comprises the following steps:
Figure BDA0001802428360000038
minimum subsystem
Figure BDA0001802428360000039
Distributed analytic redundancy relationship DARR2Comprises the following steps:
Figure BDA00018024283600000310
minimum subsystem
Figure BDA00018024283600000311
Distributed analytic redundancy relationship DARR3Comprises the following steps:
Figure BDA00018024283600000312
wherein, in DARR1In (e)1-1=Vin,e1-8=De1,e1-1And e1-8Respectively representing minimum subsystems
Figure BDA00018024283600000313
Of two inputs, Df*Representing a minimum subsystem
Figure BDA00018024283600000314
The measured value of the electric current sensor Df;
in DARR2In (e)2-1=Vin-R1Df,e2-10=De2,e2-1And e2-10Respectively representing minimum subsystems
Figure BDA00018024283600000315
The value of the two input voltages of (a),
Figure BDA00018024283600000316
representing a minimum subsystem
Figure BDA00018024283600000317
Voltage sensor De1A measured value of (a);
in DARR3In (e)3-1=De1,e3-1Representing a minimum subsystem
Figure BDA00018024283600000318
The value of the input voltage of (a),
Figure BDA00018024283600000319
representing a minimum subsystem
Figure BDA00018024283600000320
Voltage sensor De2A measured value of (a);
s32, according to the distributed analytic redundancy relation DARR of the three minimum subsystems1、DARR2、DARR3Obtaining a distributed fault characteristic matrix, wherein the distributed fault characteristic matrix is as follows:
θ/r dr1 dr2 dr3
R1 1 0 0
R 2 1 1 0
R3 0 k1 k1
C1 1 0 0
C 2 0 1 0
C 3 0 0 k1
wherein dr is1、dr2、dr3Respectively corresponding to the minimum subsystem
Figure BDA0001802428360000041
Residual error of (1), R1、R2、R3Resistive elements, C, being circuitry1、C2、C3Capacitive elements of the circuit system, wherein 1 in the distributed fault signature matrix indicates that the residual error of the column is sensitive to the elements of the row, 0 in the distributed fault signature matrix indicates that the residual error of the column is insensitive to the elements of the row, and k in the distributed fault signature matrix1The switch k represents the sensitivity of the residual of the column to the elements of the row1Determining;
due to the minimal subsystem
Figure BDA0001802428360000042
Distributed analytic redundancy relationship DARR1The element contained in (A) has R1、R2、 C1D is therefore dr1To the element R1、R2、C1Sensitive, so in the distributed fault signature matrix, the drth1R of1Line, R2Line, C1The value of the line is 1 at the dr1R of3Line, C2Line, C3The value of the row is 0;
due to the minimal subsystem
Figure BDA0001802428360000043
Distributed analytic redundancy relationship DARR2The element contained in (A) has R2、R3、 C2D is therefore dr2To the element R2、R3、C2Is sensitive and R is3And also with a switch k1Correlation, therefore dr2To the element R3Is sensitive to by the switch k1Determining, therefore, in the distributed fault signature matrix, the drth2R of2Line, C2The value of the line is 1, dr2R of3The value of a row is k1Dr th2R of1Line, C1Line, C3The value of the row is 0;
due to the minimal subsystem
Figure BDA0001802428360000044
Distributed analytic redundancy relationship DARR3The element contained in (A) has R3、C3D is therefore dr2To the element R3、C3Is sensitive and R is3、C3Are all connected with switch k1Correlation, therefore dr3To the element R3、C3Is sensitive to by the switch k1Determining, therefore, in the distributed fault signature matrix, the drth3R of3Line, C3The value of a row is k1Dr th2R of1Line, R2Line, C1Line, C2The row has a value of 0. .
The invention has the advantages that:
(1) the modeling method of the bonding diagram clearly and quantitatively describes the physical relationship among the elements of the circuit system, and is convenient for quickly and accurately positioning the specific fault element during fault diagnosis.
(2) The invention decomposes a complex circuit system into a plurality of simple and independent minimum subsystems, thereby reducing the complexity of fault isolation.
(3) Under the condition of single fault, the fault isolation of the distributed fault diagnosis method is more accurate.
(4) Under the condition of multiple faults, the distributed fault diagnosis method improves the isolation capability of the multiple faults and simultaneously reduces the complexity of fault identification after fault isolation.
(5) Under the conditions of single fault and multiple faults, the distributed fault diagnosis method effectively improves the fault isolation performance, reduces the complexity of subsequent fault element identification, and improves the diagnosis speed and the diagnosis precision of the faults.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a global diagnostic bond map model of the circuit system of the present invention.
FIG. 3(1) is a diagram of the minimal subsystem of the present invention
Figure BDA0001802428360000051
Local diagnostic bond map model.
FIG. 3(2) is a diagram of the minimal subsystem of the present invention
Figure BDA0001802428360000052
Local diagnostic bond map model.
FIG. 3(3) is a diagram of the minimal subsystem of the present invention
Figure BDA0001802428360000053
Local diagnostic bond map model.
FIG. 4(a) shows an injection resistor R according to the present invention1Global analytical redundancy relationship GARR at fault1The waveform of (2).
FIG. 4(b) shows an injection resistor R according to the present invention1Global analytical redundancy relationship GARR at fault2The waveform of (2).
FIG. 4(c) shows an injection resistor R according to the present invention1Global analytical redundancy relationship GARR at fault3The waveform of (2).
FIG. 4(d) shows an injection resistor R according to the present invention1Distributed analytic redundancy relationship DARR upon failure1The waveform of (2).
FIG. 4(e) shows an injection resistor R according to the present invention1Distributed analytic redundancy relationship DARR upon failure2The waveform of (2).
FIG. 4(f) shows an injection resistor R according to the present invention1Distributed analytic redundancy relationship DARR upon failure3The waveform of (2).
FIG. 5(a) shows an injection resistor R according to the present invention1And a capacitor C3Global analytical redundancy relationship GARR during simultaneous failure1The waveform of (2).
FIG. 5(b) shows an injection resistor R according to the present invention1And a capacitor C3Global analytical redundancy relationship GARR during simultaneous failure2The waveform of (2).
FIG. 5(c) shows an injection resistor R according to the present invention1And a capacitor C3Global analytical redundancy relationship GARR during simultaneous failure3The waveform of (2).
FIG. 5(d) shows an injection resistor R according to the present invention1And a capacitor C3Distributed analytic redundancy relationship DARR during simultaneous failure1The waveform of (2).
FIG. 5(e) shows an injection resistor R according to the present invention1And a capacitor C3Distributed analytic redundancy relationship DARR during simultaneous failure2The waveform of (2).
FIG. 5(f) shows an injection resistor R according to the present invention1And a capacitor C3Distributed analytic redundancy relationship DARR during simultaneous failure3The waveform of (2).
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In this embodiment, the advantages of the distributed fault diagnosis method of the present invention are highlighted by comparing the fault diagnosis effects based on the global fault feature matrix and the distributed fault feature matrix.
As shown in fig. 1, a distributed fault diagnosis method for a circuit system includes the following steps:
and S1, taking a complex circuit system as an example, modeling the circuit system, configuring the causal relationship among all parameters, and obtaining a global diagnosis bonding diagram model of the circuit system.
And S2, taking each sensor in the circuit system as a basic unit, taking the measured value of the sensor as the local output of the circuit system, extracting the minimum subsystem based on each sensor from the circuit system, and respectively obtaining the local diagnosis bonding diagram model of each minimum subsystem.
And S3, obtaining a distributed analytic redundancy relationship according to the local diagnostic bond map model of the minimum subsystem, and analyzing the distributed analytic redundancy relationship to obtain a distributed fault feature matrix.
And S4, performing fault diagnosis according to the distributed fault feature matrix.
As shown in fig. 2, in step S1, the parameters in the global diagnostic bond map model of the circuit system include: potential source UM=VinResistive element { R1,R2,R3}, capacitive element { C1,C2,C3}, voltage sensor { De1,De2}, current sensor Df, switch k1. Wherein, the potential source VinIs the input voltage of the system; resistance R1And a capacitor C1Series connected, resistance R2And a capacitor C2Series connected, resistance R3And a capacitor C3Are connected in series; will resistance R1And a capacitor C1Formed series circuit, resistor R2And a capacitor C2Formed series circuit, resistor R3And a capacitor C3The formed series circuits are connected in parallel; voltage sensor De1Connected in parallel to a capacitor C2Both ends of (a); voltage sensor De2Connected in parallel to a capacitor C3Both ends of (a); current sensor Df and resistor R1Potential source VinAre connected in series; switch k1And a resistor R3Are connected in series.
In the global diagnosis bond-graph model, directional arrows represent keys, 0 and 1 represent two nodes, namely a common potential node and a common current node, wherein the potentials on the keys connected with the common potential node are the same, and the currents on the keys connected with the common current node are the same. In thatDescribing constraints for nodes and parameters in the diagnostic bond graph, eiIndicating the potential, i.e. voltage, f, at the bond numbered iiIndicating a flow, i.e., current, over the key numbered i. In this example, VinAnd R1Df are connected in series, the currents flowing through the three elements are the same, and a common current junction 1 is used in a bonding diagram1Connection, VinAnd C1Parallel connection, using a common potential junction 0 in a bonding diagram2Is connected to R2In series with the front-end circuit, the same current flows as the front-end circuit, and a common current junction 1 is used in the bonding diagram3Connected to the front end, C2、De1Connected in parallel with the front-end circuit and at the same voltage, and the common potential junction 0 is used in the bonding diagram4To the front end, R3、k1The same current flows in series with the front end, and a common current junction 1 is used in the bonding diagram5Connected to a front end, De2、C3Parallel connection with front end, same voltage, common potential junction 0 in bonding diagram6Connected to the front end, thereby forming a bond map model of the circuitry used in the present invention.
In the global diagnostic bond map model of the circuit system, the constraint relationship between each node and key is as follows:
Figure BDA0001802428360000071
in step S2, the manner of obtaining the local diagnostic bond map model of the smallest subsystem is specifically as follows: taking a certain sensor of the circuit system as an output sensor, taking the measured value of the output sensor as the local output of the circuit system, namely taking the measured value of the output sensor as the output signal of the minimum subsystem based on the output sensor, starting from the output sensor to the end of another sensor, taking the measured value of the sensor which is the end sensor as the input signal of the minimum subsystem, and obtaining a local diagnostic bonding map model of the minimum subsystem based on the output sensor; in a path from the output sensor to the termination sensor, a parameter associated with the measured value of the output sensor is derived from causal relationships between various parameters on a global diagnostic bond map model of the circuitry, and the parameter associated with the measured value of the output sensor is used as a parameter of a local diagnostic bond map model based on a smallest subsystem of the output sensor.
As shown in fig. 3(1), the current sensor Df is used as a basic unit from the current sensor Df to the voltage sensor De1Terminate, extract the smallest subsystem based on the current sensor Df
Figure RE-GDA0001861329930000081
The parameters of the local diagnostic bond map model of the minimal subsystem include: potential source
Figure RE-GDA0001861329930000082
Resistive element { R1,R2}, capacitive element C1Flow sensor Df.
Minimum subsystem
Figure RE-GDA0001861329930000083
The local diagnostic bond map model of (2) wherein the constraint relationship between each node and bond is as follows:
Figure RE-GDA0001861329930000084
Figure RE-GDA0001861329930000085
shown in FIG. 3(2), with a voltage sensor De1As a basic unit, a slave voltage sensor De1To start, to the voltage sensor De2And the current sensor Df is terminated, and the voltage sensor De is extracted1Minimum subsystem of
Figure RE-GDA0001861329930000086
The minimum subsystem
Figure RE-GDA0001861329930000087
The parameters of the local diagnostic bond map model of (a) include: potential source
Figure RE-GDA0001861329930000088
Figure RE-GDA0001861329930000089
Resistive element { R2,R3}, capacitive element { C1,C2}, a voltage sensor De1Switch k1
Minimum subsystem
Figure RE-GDA00018613299300000810
The local diagnostic bond map model of (2) wherein the constraint relationship between each node and bond is as follows:
Figure RE-GDA0001861329930000091
Figure RE-GDA0001861329930000092
as shown in the graph of figure 3(3),
with a voltage sensor De2As a basic unit, a slave voltage sensor De2To start, to the voltage sensor De1Terminate, extract based on the voltage sensor De2Minimum subsystem of
Figure BDA0001802428360000093
The minimum subsystem
Figure BDA0001802428360000094
The parameters of the local diagnostic bond map model of (a) include: potential source
Figure BDA0001802428360000095
Resistive element R3Capacitive element C3D, a voltage sensor De2Switch k1
Minimum subsystem
Figure BDA0001802428360000096
The local diagnostic bond map model of (2) wherein the constraint relationship between each node and bond is as follows:
Figure BDA0001802428360000097
Figure BDA0001802428360000098
in step S3, the method includes the following steps:
and S31, obtaining a distributed analytic redundancy relation through a causal path overlay method according to the local diagnostic bond map model of the minimum subsystem and the causal relation among all parameters in the local diagnostic bond map model.
The specific manner of the causal path coverage is as follows: selecting a node in the bonding graph model, listing related analytical redundancy relations through constraint relations between the selected node and the key, replacing unknown variables in the formula with known variables or measurable variables connected with the node through a causal path, removing all the unknown variables, and obtaining an analytical redundancy relation formula after removing all the unknown variables. And after obtaining an analytic redundancy relation, selecting another node, listing the analytic redundancy relation related to the node according to the mode, checking whether the listed analytic redundancy relation is independent from other analytic redundancy relations, if so, keeping the analytic redundancy relation, if not, selecting the next node, and so on until all nodes are selected.
In this embodiment, the minimum subsystem is obtained
Figure BDA0001802428360000101
Distributed resolving of redundant relationshipsDARR1Comprises the following steps:
Figure BDA0001802428360000102
minimum subsystem
Figure BDA0001802428360000103
Distributed analytic redundancy relationship DARR2Comprises the following steps:
Figure BDA0001802428360000104
minimum subsystem
Figure BDA0001802428360000105
Distributed analytic redundancy relationship DARR3Comprises the following steps:
Figure BDA0001802428360000106
wherein, in DARR1In (e)1-1=Vin,e1-8=De1,e1-1And e1-8Respectively representing minimum subsystems
Figure BDA0001802428360000107
Of two inputs, Df*Representing a minimum subsystem
Figure BDA0001802428360000108
The measured value of the electric current sensor Df of (a),
Figure BDA0001802428360000109
represents a differential operation over time;
in DARR2In (e)2-1=Vin-R1Df,e2-10=De2,e2-1And e2-10Respectively representing minimum subsystems
Figure BDA00018024283600001010
The value of the two input voltages of (a),
Figure BDA00018024283600001011
representing a minimum subsystem
Figure BDA00018024283600001012
Voltage sensor De1Measured value of (a), k1Represents the open and close state of the switch in the circuit system, the open state is 0, the close state is 1,
Figure BDA00018024283600001013
represents a differential operation over time;
in DARR3In (e)3-1=De1,e3-1Representing a minimum subsystem
Figure BDA00018024283600001014
The value of the input voltage of (a),
Figure BDA00018024283600001015
representing a minimum subsystem
Figure BDA00018024283600001016
Voltage sensor De2Measured value of (a), k1Represents the open and close state of the switch in the circuit system, the open state is 0, the close state is 1,
Figure BDA00018024283600001017
representing a differential operation over time.
S32, according to the distributed analytic redundancy relation DARR of the three minimum subsystems1、DARR2、DARR3Obtaining a distributed fault characteristic matrix, wherein the distributed fault characteristic matrix is as follows:
Figure BDA00018024283600001018
Figure BDA0001802428360000111
wherein dr is1、dr2、dr3Respectively corresponding to the minimum subsystem
Figure BDA0001802428360000112
Residual error of (1), R1、R2、R3Resistive elements, C, being circuitry1、C2、C3Capacitive elements of the circuit system, wherein 1 in the distributed fault signature matrix indicates that the residual error of the column is sensitive to the elements of the row, 0 in the distributed fault signature matrix indicates that the residual error of the column is insensitive to the elements of the row, and k in the distributed fault signature matrix1The switch k represents the sensitivity of the residual of the column to the elements of the row1And (6) determining.
Due to the minimal subsystem
Figure BDA0001802428360000113
Distributed analytic redundancy relationship DARR1The element contained in (A) has R1、R2、 C1D is therefore dr1To the element R1、R2、C1Sensitive, so in the distributed fault signature matrix, the drth1R of1Line, R2Line, C1The value of the line is 1 at the dr1R of3Line, C2Line, C3The row has a value of 0.
Due to the minimal subsystem
Figure BDA0001802428360000114
Distributed analytic redundancy relationship DARR2The element contained in (A) has R2、R3、 C2D is therefore dr2To the element R2、R3、C2Is sensitive and R is3And also with a switch k1Correlation, therefore dr2To the element R3Is sensitive to by the switch k1Determining, therefore, in the distributed fault signature matrix, the drth2R of2Line, C2The value of the line is 1, dr2R of3The value of a row is k1Dr th2R of1Line, C1Line, C3The row has a value of 0.
Due to the minimal subsystem
Figure BDA0001802428360000115
Distributed analytic redundancy relationship DARR3The element contained in (A) has R3、C3D is therefore dr2To the element R3、C3Is sensitive and R is3、C3Are all connected with switch k1Correlation, therefore dr3To the element R3、C3Is sensitive to by the switch k1Determining, therefore, in the distributed fault signature matrix, the drth3R of3Line, C3The value of a row is k1Dr th2R of1Line, R2Line, C1Line, C2The row has a value of 0.
In step S4, the element parameters and sensor data of the minimum subsystem are substituted into the corresponding distributed analytical redundancy relationship to obtain the residual error of the minimum subsystem, and the residual errors dr of the three minimum subsystems are substituted1、dr2、 dr3A set of residues is constructed.
The residual error is a numerical estimation for analyzing the redundancy relation and is a basis for judging whether the circuit system has a fault, if the circuit system has no fault, all residual values in the circuit system are 0, and if the circuit system has a fault, at least one residual value in the circuit system is 1. The residual error fluctuates in a certain range due to the problems of external noise, disturbance, element parameter tolerance and measurement accuracy of the sensor, and in order to prevent misdiagnosis caused by the fluctuation, the residual errors dr are respectively detected1、dr2、dr3Setting a residual threshold, and if the fluctuation of the residual exceeds the residual threshold, determining that the residual is not in the normal stateThe residual error value is 1, namely at least one element in the minimum subsystem corresponding to the residual error has a fault; and if the fluctuation of the residual is within the residual threshold value, setting the residual value to be 0, namely that all the element parameters in the minimum subsystem corresponding to the residual are normal values. The size of the residual threshold is determined according to the difference of the diagnosis object and the difference of the application environment, and specifically, the residual threshold is determined according to the experiment of the diagnosis object and the application environment.
FIG. 4(d) shows a DARR for distributed resolution of redundancy relationships1The waveform of (1), the fluctuation of which exceeds a residual threshold, the corresponding residual dr1Has a value of 1; FIG. 4(f) shows a DARR for distributed resolution of redundancy relationships3The fluctuation of the waveform does not exceed the residual error threshold, corresponding residual error dr3The value of (d) is 0.
A global fault diagnosis method for a circuit system, comprising the steps of:
and S21, taking a complex circuit system as an example, modeling the circuit system, configuring the causal relationship among all parameters, and obtaining a global diagnosis bonding diagram model of the circuit system.
And S22, obtaining a global analytic redundancy relationship according to the global diagnostic bond diagram model of the circuit system, and analyzing the global analytic redundancy relationship to obtain a global fault feature matrix.
And S23, performing fault diagnosis according to the global fault feature matrix.
Step S21 in the global diagnosis method of a circuit system coincides with step S1 in the distributed diagnosis method of a circuit system.
As shown in fig. 2, in step S21, the parameters in the global diagnostic bond map model of the circuit system include: potential source UM=VinResistive element { R1,R2,R3}, capacitive element { C1,C2,C3}, voltage sensor { De1,De2}, current sensor Df, switch k1. Wherein, the potential source VinIs the input voltage of the system; resistance R1And a capacitor C1Series connected, resistance R2And a capacitor C2Series connected, resistance R3And a capacitor C3Are connected in series; will resistance R1And a capacitor C1Formed series circuit, resistor R2And a capacitor C2Formed series circuit, resistor R3And a capacitor C3The formed series circuits are connected in parallel; voltage sensor De1Connected in parallel to a capacitor C2Both ends of (a); voltage sensor De2Connected in parallel to a capacitor C3Both ends of (a); current sensor Df and resistor R1Potential source VinAre connected in series; switch k1And a resistor R3Are connected in series.
In the global diagnosis bond-graph model, directional arrows represent keys, 0 and 1 represent two nodes, namely a common potential node and a common current node, wherein the potentials on the keys connected with the common potential node are the same, and the currents on the keys connected with the common current node are the same. In describing the constraints for each node and each parameter in the diagnostic bond graph, eiIndicating the potential, i.e. voltage, f, at the bond numbered iiIndicating a flow, i.e., current, over the key numbered i. In this example, VinAnd R1Df are connected in series, the currents flowing through the three elements are the same, and a common current junction 1 is used in a bonding diagram1Connection, VinAnd C1Parallel connection, using a common potential junction 0 in a bonding diagram2Is connected to R2In series with the front-end circuit, the same current flows as the front-end circuit, and a common current junction 1 is used in the bonding diagram3Connected to the front end, C2、De1Connected in parallel with the front-end circuit and at the same voltage, and the common potential junction 0 is used in the bonding diagram4To the front end, R3、k1The same current flows in series with the front end, and a common current junction 1 is used in the bonding diagram5Connected to a front end, De2、C3Parallel connection with front end, same voltage, common potential junction 0 in bonding diagram6Connected to the front end, thereby forming a bond map model of the circuitry used in the present invention.
In the global diagnostic bond map model of the circuit system, the constraint relationship between each node and key is as follows:
Figure BDA0001802428360000131
in step S22, the method includes the following steps:
and S221, obtaining a global analysis redundancy relation through a causal path overlay method according to the global diagnosis bonding diagram model of the circuit system and the causal relation among all parameters in the global diagnosis bonding diagram model.
The specific manner of the causal path coverage is as follows: selecting a node in the bonding graph model, listing related analytical redundancy relations through constraint relations between the selected node and the key, replacing unknown variables in the formula with known variables or measurable variables connected with the node through a causal path, removing all the unknown variables, and obtaining an analytical redundancy relation formula after removing all the unknown variables. And after obtaining an analytic redundancy relation, selecting another node, listing the analytic redundancy relation related to the node according to the mode, checking whether the listed analytic redundancy relation is independent from other analytic redundancy relations, if so, keeping the analytic redundancy relation, if not, selecting the next node, and so on until all nodes are selected.
In this embodiment, the obtained global analytic redundancy relationships of the circuit system are respectively:
Figure BDA0001802428360000141
Figure BDA0001802428360000142
Figure BDA0001802428360000143
wherein, VinIs the input voltage value, R, of the circuit system1,R2,R3Respectively representing the resistance values, C, in the circuit system1,C2,C3Respectively representing the capacitance values, De, in the circuitry1,De2Respectively representing the measured value of a voltage sensor in the circuit system, Df representing the measured value of a current sensor in the circuit system, k1Represents the open and close state of the switch in the circuit system, the open state is 0, the close state is 1,
Figure BDA0001802428360000144
representing a differential operation over time.
S222, according to the global analysis redundancy relation GARR of the circuit system1,GARR2,GARR3Obtaining a global fault characteristic matrix, wherein the global fault characteristic matrix is as follows:
θ/r gr1 gr2 gr3
R1 1 1 0
R 2 1 1 0
R 3 0 0 k1
C1 1 0 0
C 2 0 1 0
C3 0 k1 k1
wherein, gr1、gr2、gr3For three residuals of the circuitry, R1、R2、R3Resistive elements, C, being circuitry1、C2、C3The capacitive elements are circuit systems, 1 in the global fault feature matrix indicates that the residual error of the column is sensitive to the elements of the row, 0 in the global fault feature matrix indicates that the residual error of the column is insensitive to the elements of the row, and k in the global fault feature matrix1The switch k represents the sensitivity of the residual of the column to the elements of the row1And (6) determining.
In step S23, substituting the element parameters and sensor data of the circuit system into the corresponding global analytic redundancy relationship to obtain the residual error of the circuit system, and obtaining the three residual errors gr of the circuit system1、gr2、gr3A set of residues is constructed.
Based on a global fault signature matrix, R, when a single fault occurs in a circuit system1And R2All the same fault signatures in the global fault signature matrix are [ 110 ]]As shown in FIGS. 4(a), (b), and (c), three residuals gr from the circuit system are observed1、gr2、gr3The set of constructed residuals is [ 110 ]]Then the possible faulty element obtained by global fault diagnosis is R1And R2Result in R1And R2Fault isolation cannot be achieved in the circuit system; based on a distributed fault signature matrix, R1And R2Respectively is [100 ]]And [ 110 ]]That is, if the corresponding failure characteristics are different, as shown in fig. 4(d), (e), and (f), the residual dr of the three smallest subsystems is observed1、dr2、dr3The set of constructed residuals is [100 ]]Then the possible failure element obtained by the distributed failure feature matrix is R1. Therefore, the fault isolation based on the distributed fault signature matrix of the invention is more accurate than the fault isolation based on the global fault signature matrix.
When multiple faults occur in the circuit system, e.g. switch k1Closed, i.e. k 11, element R in the circuit system1、 C3The simultaneous failures are based on the global failure signature matrix, as shown in FIGS. 5(a), (b), and (c), from the three residuals gr of the circuitry1、gr2、gr3The set of constructed residuals is [ 111 ]]At this time, it cannot be determined whether other elements in the circuit system are in failure at the same time, so that in the subsequent failure element identification process, all elements in the circuit system need to be identified, thereby increasing the complexity of failure diagnosis; based on the distributed fault feature matrix, shown in the graphs (d), (e) and (f), the residual error dr of the three minimum subsystems1、dr2、dr3The set of constructed residuals is [ 101 ]]The possible failure element obtained by the distributed failure feature matrix is R1、C1、C3. Therefore, the present inventionThe clear distributed fault diagnosis method improves the isolation capability of multiple faults and simultaneously reduces the complexity of fault identification after fault isolation.
The invention is not to be considered as limited to the specific embodiments shown and described, but is to be understood to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention as defined by the appended claims.

Claims (4)

1. A distributed fault diagnosis method for a circuit system, comprising the steps of:
s1, modeling the circuit system, and configuring the causal relationship among all parameters to obtain a global diagnosis bonding diagram model of the circuit system;
s2, taking each sensor in the circuit system as a basic unit, taking the measured value of the sensor as the local output of the circuit system, extracting the minimum subsystem based on each sensor from the circuit system, and respectively obtaining the local diagnosis bond map model of each minimum subsystem;
s3, obtaining a distributed analytic redundancy relationship according to the local diagnostic bond map model of the minimum subsystem, and analyzing the distributed analytic redundancy relationship to obtain a distributed fault feature matrix;
s4, performing fault diagnosis according to the distributed fault feature matrix;
in step S2, the manner of obtaining the local diagnostic bond map model of the smallest subsystem is specifically as follows: taking a certain sensor of the circuit system as an output sensor, taking the measured value of the output sensor as the local output of the circuit system, namely taking the measured value of the output sensor as the output signal of the minimum subsystem based on the output sensor, starting from the output sensor to the end of another sensor, taking the measured value of the sensor which is the end sensor as the input signal of the minimum subsystem, and obtaining a local diagnostic bonding map model of the minimum subsystem based on the output sensor; in a path from the output sensor to the termination sensor, a parameter associated with the measured value of the output sensor is derived from causal relationships between various parameters on a global diagnostic bond map model of the circuitry, and the parameter associated with the measured value of the output sensor is used as a parameter of a local diagnostic bond map model based on a smallest subsystem of the output sensor.
2. The distributed fault diagnosis method for circuit systems according to claim 1, wherein in step S1, the parameters of the global diagnostic bond map model for circuit systems include: potential source UM=VinResistive element { R1,R2,R3}, capacitive element { C1,C2,C3}, voltage sensor { De1,De2}, current sensor Df, switch k1(ii) a Wherein, the potential source VinIs the input voltage of the system; resistance R1And a capacitor C1Series connected, resistance R2And a capacitor C2Series connected, resistance R3And a capacitor C3Are connected in series; will resistance R1And a capacitor C1Formed series circuit, resistor R2And a capacitor C2Formed series circuit, resistor R3And a capacitor C3The formed series circuits are connected in parallel; voltage sensor De1Connected in parallel to a capacitor C2Both ends of (a); voltage sensor De2Connected in parallel to a capacitor C3Both ends of (a); current sensor Df and resistor R1Potential source VinAre connected in series; switch k1And a resistor R3Are connected in series.
3. The distributed fault diagnosis method of circuit system according to claim 2, wherein the current sensor Df is used as a basic unit, starting from the current sensor Df to the voltage sensor De1Terminate, extract the smallest subsystem based on the current sensor Df
Figure FDA0002592273090000021
The parameters of the local diagnostic bond map model of the minimal subsystem include: potential source
Figure FDA0002592273090000022
Resistive element { R1,R2}, capacitive element C1A flow sensor Df;
with a voltage sensor De1As a basic unit, a slave voltage sensor De1To start, to the voltage sensor De2And the current sensor Df is terminated, and the voltage sensor De is extracted1Minimum subsystem of
Figure FDA0002592273090000023
The minimum subsystem
Figure FDA0002592273090000024
The parameters of the local diagnostic bond map model of (a) include: potential source
Figure FDA0002592273090000025
Resistive element { R2,R3}, capacitive element { C1,C2}, a voltage sensor De1Switch k1
With a voltage sensor De2As a basic unit, a slave voltage sensor De2To start, to the voltage sensor De1Terminate, extract based on the voltage sensor De2Minimum subsystem of
Figure FDA0002592273090000026
The minimum subsystem
Figure FDA0002592273090000027
The parameters of the local diagnostic bond map model of (a) include: potential source
Figure FDA0002592273090000028
Resistive element R3Capacitive element C3D, a voltage sensor De2Switch k1
4. The distributed fault diagnosis method for circuit systems according to claim 3, wherein step S3 includes the following specific steps:
s31, obtaining a distributed analytical redundancy relationship according to the local diagnostic bond map models of the three minimum subsystems and the causal relationship among all parameters in the local diagnostic bond map models;
in the above manner, the minimum subsystem is obtained
Figure FDA0002592273090000029
Distributed analytic redundancy relationship DARR1Comprises the following steps:
Figure FDA00025922730900000210
minimum subsystem
Figure FDA00025922730900000211
Distributed analytic redundancy relationship DARR2Comprises the following steps:
Figure FDA00025922730900000212
minimum subsystem
Figure FDA00025922730900000213
Distributed analytic redundancy relationship DARR3Comprises the following steps:
Figure FDA00025922730900000214
wherein, in DARR1In (e)1-1=Vin,e1-8=De1,e1-1And e1-8Respectively representing minimum subsystems
Figure FDA0002592273090000031
Of two inputs, Df*Representing a minimum subsystem
Figure FDA0002592273090000032
The measured value of the electric current sensor Df;
in DARR2In (e)2-1=Vin-R1Df,e2-10=De2,e2-1And e2-10Respectively representing minimum subsystems
Figure FDA0002592273090000033
The value of the two input voltages of (a),
Figure FDA0002592273090000034
representing a minimum subsystem
Figure FDA0002592273090000035
Voltage sensor De1A measured value of (a);
in DARR3In (e)3-1=De1,e3-1Representing a minimum subsystem
Figure FDA0002592273090000036
The value of the input voltage of (a),
Figure FDA0002592273090000037
representing a minimum subsystem
Figure FDA0002592273090000038
Voltage sensor De2A measured value of (a);
s32, according to the distributed analytic redundancy relation DARR of the three minimum subsystems1、DARR2、DARR3Obtaining a distributed fault characteristic matrix, wherein the distributed fault characteristic matrix is as follows:
θ/r dr1 dr2 dr3 R1 1 0 0 R2 1 1 0 R3 0 k1 k1 C1 1 0 0 C2 0 1 0 C3 0 0 k1
wherein dr is1、dr2、dr3Respectively corresponding to the minimum subsystem
Figure FDA0002592273090000039
Residual error of (1), R1、R2、R3Resistive elements, C, being circuitry1、C2、C3Capacitive elements of the circuit system, wherein 1 in the distributed fault signature matrix indicates that the residual error of the column is sensitive to the elements of the row, 0 in the distributed fault signature matrix indicates that the residual error of the column is insensitive to the elements of the row, and k in the distributed fault signature matrix1The switch k represents the sensitivity of the residual of the column to the elements of the row1Determining;
due to the minimal subsystem
Figure FDA00025922730900000310
Distributed analytic redundancy relationship DARR1The element contained in (A) has R1、R2、C1D is therefore dr1To the element R1、R2、C1Sensitive, so in the distributed fault signature matrix, the drth1R of1Line, R2Line, C1The value of the line is 1 at the dr1R of3Line, C2Line, C3The value of the row is 0;
due to the minimal subsystem
Figure FDA00025922730900000311
Distributed analytic redundancy relationship DARR2The element contained in (A) has R2、R3、C2D is therefore dr2To the element R2、R3、C2Is sensitive and R is3And also with a switch k1Correlation, therefore dr2To the element R3Is sensitive to by the switch k1Determining, therefore, in the distributed fault signature matrix, the drth2R of2Line, C2The value of the line is 1, dr2R of3The value of a row is k1Dr th2R of1Line, C1Line, C3The value of the row is 0;
due to the minimal subsystem
Figure FDA0002592273090000041
Distributed analytic redundancy relationship DARR3The element contained in (A) has R3、C3D is therefore dr2To the element R3、C3Is sensitive and R is3、C3Are all connected with switch k1Correlation, therefore dr3To the element R3、C3Is sensitive to by the switch k1Determining, therefore, in the distributed fault signature matrix, the drth3R of3Line, C3The value of a row is k1Dr th2R of1Line, R2Line, C1Line, C2The row has a value of 0.
CN201811082945.7A 2018-09-17 2018-09-17 Distributed fault diagnosis method for circuit system Active CN109116219B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811082945.7A CN109116219B (en) 2018-09-17 2018-09-17 Distributed fault diagnosis method for circuit system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811082945.7A CN109116219B (en) 2018-09-17 2018-09-17 Distributed fault diagnosis method for circuit system

Publications (2)

Publication Number Publication Date
CN109116219A CN109116219A (en) 2019-01-01
CN109116219B true CN109116219B (en) 2021-04-06

Family

ID=64859657

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811082945.7A Active CN109116219B (en) 2018-09-17 2018-09-17 Distributed fault diagnosis method for circuit system

Country Status (1)

Country Link
CN (1) CN109116219B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230259113A1 (en) * 2022-02-11 2023-08-17 Novity, Inc. Subsystem-level model-based diagnosis
CN115083093A (en) * 2022-05-16 2022-09-20 浙江台谊消防股份有限公司 Electrical fire monitoring system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002108949A (en) * 2000-10-04 2002-04-12 Tadano Ltd Bond graph analysis system
WO2012119241A1 (en) * 2011-03-04 2012-09-13 Kylowave Inc. System and method for generating discrete-time model (dtm) of continuous-time model (ctm) for a dynamical system
CN103853889A (en) * 2014-03-12 2014-06-11 南京航空航天大学 Power electronic system fault diagnosis method based on bond graph
CN104376206A (en) * 2014-11-14 2015-02-25 浙江工业大学 Large-scale reaction kettle distributed fault diagnosis method based on sensor network
CN105301470A (en) * 2015-08-12 2016-02-03 南京航空航天大学 IGBT fault diagnosis method for traction drive system of CRH5 type high-speed train
CN108398637A (en) * 2018-01-29 2018-08-14 合肥工业大学 A kind of method for diagnosing faults of non-linear Mechatronic Systems

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IT1395258B1 (en) * 2009-07-23 2012-09-05 Fond Istituto Italiano Di Tecnologia PROCEDURE FOR THE GENERATION OF A SET OF CONFLICTS FOR THE MODEL-BASED DIAGNOSIS OF SYSTEMS, AND CORRESPONDING DIAGNOSIS PROCEDURE

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002108949A (en) * 2000-10-04 2002-04-12 Tadano Ltd Bond graph analysis system
WO2012119241A1 (en) * 2011-03-04 2012-09-13 Kylowave Inc. System and method for generating discrete-time model (dtm) of continuous-time model (ctm) for a dynamical system
CN103853889A (en) * 2014-03-12 2014-06-11 南京航空航天大学 Power electronic system fault diagnosis method based on bond graph
CN104376206A (en) * 2014-11-14 2015-02-25 浙江工业大学 Large-scale reaction kettle distributed fault diagnosis method based on sensor network
CN105301470A (en) * 2015-08-12 2016-02-03 南京航空航天大学 IGBT fault diagnosis method for traction drive system of CRH5 type high-speed train
CN108398637A (en) * 2018-01-29 2018-08-14 合肥工业大学 A kind of method for diagnosing faults of non-linear Mechatronic Systems

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于第一原理的实用故障诊断算法;方敏等;《东南大学学报》;19990520;第29卷(第3期);83-88 *

Also Published As

Publication number Publication date
CN109116219A (en) 2019-01-01

Similar Documents

Publication Publication Date Title
Sunter et al. Test metrics for analog parametric faults
JP4557337B2 (en) Method and system for diagnosing multiple errors and faults based on X list
Huang et al. Fault diagnosis of analog circuits based on machine learning
EP1922554B1 (en) Current measurement circuit and method of diagnosing faults in same
EP3617723B1 (en) Integrated circuit device with integrated fault monitoring system
CN108398637B (en) Fault diagnosis method for nonlinear electromechanical system
CN109116219B (en) Distributed fault diagnosis method for circuit system
CN109032872B (en) Bayesian network-based equipment fault diagnosis method and system
US9773358B2 (en) Device and method for detecting faults in electronic systems
CN112630682A (en) Fault detection method, device and equipment of sensor
CN110619922A (en) Semiconductor fault analysis device and fault analysis method thereof
CN110441669B (en) Gradual fault diagnosis and service life prediction method for uncertain hybrid circuit system
US10683101B2 (en) Method for determining at least one faulty piece of equipment of an aircraft, and corresponding system
Parai et al. Multi-source data fusion technique for parametric fault diagnosis in analog circuits
US10865726B2 (en) Fault diagnosis of electronic control unit (ECU)
CN103941136A (en) Signal sampling circuit resistor-capacitor component short circuit and open circuit fault diagnosis method
CN108845247B (en) Fault diagnosis method for analog circuit module
Chakrabarti et al. Fault diagnosis for mixed-signal electronic systems
Sindia et al. Polynomial coefficient based DC testing of non-linear analog circuits
Yang et al. The effect of time-between-events for sequence interaction testing of a real-time system
Zilch et al. A versatile test set generation tool for structural analog circuit testing
CN107957269B (en) Inertial navigation system fault characteristic judgment and testability prediction method
US8049524B2 (en) Method for detecting component defects of an analog signal processing circuit, especially for a measurement transmitter
CN102915032A (en) Method for detecting state of electronic control system and diagnosing faults
Pekpe et al. Subspace method for sensor fault detection and isolation-application to grinding circuit monitoring

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant