Patent Grant
11609830
The invention relates to a method for investigating a functional behavior of a component of a technical installation.
It is known that those components of a compnt type which are exhibiting unusual or abnormal behavior can be detected or identified as being non-functional. For example, a current operating parameter value of a component to be investigated is compared with a fixed threshold value. For example, if the operating parameter value of the component to be tested is greater than the fixed threshold, it is detected that the component to be investigated is displaying an abnormal behavior. Determining which threshold value to use is difficult and subject to great uncertainty. Often, the threshold values are specified on the basis of experience values or are arbitrarily chosen. External factors such as ambient conditions (e.g. ambient temperature, humidity, etc.) and/or installation-dependent factors (e.g. running time of the installation, operating temperatures of surrounding components, etc.) are ignored as far as possible.
Thus, for example, the specifications DE 199 61 631 A1, DE 37 34 487 A1, DE 197 32 046 A1 and EP 1 542 108 A1 teach different methods for monitoring a component of a technical installation, in each of which a signal representing the functional behavior of the component to be investigated is compared with a predefined reference signal.
An object of the invention is to specify an improved, more accurate method for investigating a functional behavior of a component of a technical installation.
The object is achieved by a method for investigating a functional behavior of a component of a technical installation, in which according to the invention a signal of the component that represents the functional behavior of the component is compared with a reference signal which describes an average functional behavior of identical components. According to the invention, during the comparison a comparison variable describing the deviation of the signal from the reference signal is determined. This means that the comparison variable preferably describes the deviation of the signal of the component to be investigated from the reference signal. Further, a probability of occurrence of the comparison variable is determined using a definable distribution of a plurality of such comparison variables.
Advantageously, the functional behavior is assessed on the basis of the probability of occurrence of the comparison variable. The probability of occurrence of the comparison variable is preferably a measure of the functional behavior of the component to be investigated, in particular a measure of whether the component to be investigated is behaving normally, in particular whether the component to be investigated is fully functional.
The probability of occurrence of the comparison variable is advantageously dependent on the plurality of such comparison variables, on the basis of which, in particular, the definable distribution is determined. The plurality of such comparison variables can be updated, in particular periodically. In this way, the functional behavior of the component to be investigated can be tested and/or assessed using a dynamic criterion, namely the probability of occurrence.
A further advantage of the method according to the invention is its high sensitivity. This means that even small deviations of the signal of the component to be investigated from the reference signal can be detected.
The signal from the signal component to be investigated can be an unprocessed signal. In addition, the signal from the component to be investigated can be a post-processed signal. In particular, the signal of the component to be investigated can be filtered and/or smoothed. Furthermore, the signal of the component to be investigated can be an average value of multiple signals and/or a convolution of multiple signals.
In an advantageous design of the invention the comparison variable is a maximum cumulative deviation between the signal of the component to be investigated and the reference signal.
The signal from the component to be tested can comprise a plurality of operating parameter values, in particular as a function of time. Furthermore, the reference signal can comprise multiple reference values, in particular as a function of time.
The signal of the components to be tested and/or the reference signal can be smoothed. Furthermore, the signal of the component to be investigated and/or the reference signal can be normalized. In particular, the operating parameter values of the signal, or the operating parameter values for short, and/or the reference values of the reference signal, the reference values for short, can be normalized.
It is advantageous if the operating parameter values are accumulated over time. Further, it is preferable if the reference values are accumulated over time. Preferably, a deviation, also a cumulative deviation, between the cumulative operating parameter values and the cumulative reference values is calculated as a function of time.
Advantageously, a maximum cumulative deviation, also a maximum cumulative distance, between the cumulative operating parameter values and the cumulative reference values is determined as the comparison value. The maximum cumulative deviation corresponds advantageously to the maximum of the above-mentioned deviation between the cumulative operating parameter values and the cumulative reference values.
An operating parameter can be, for example, a temperature, a vibration parameter, a loudness, a voltage, an electric current, a consumed/converted power, a rotational frequency or similar. Further, the operating parameter can be one which can be determined from a plurality of parameters.
Advantageously, the operating parameter values and/or the reference values are time-dependent. In this way, the signal from the component to be investigated and/or the reference signal can be time-dependent.
For example, the operating parameter values and/or the reference values may be or have been determined directly as a function of time. Further, the operating parameter values and/or the reference values can be determined indirectly as a function of time, for example as a function of a running time, which is in turn time-dependent.
In addition, the operating parameter values and/or the reference values may be accumulated directly over time. In addition, the operating parameter values and/or the reference values may be accumulated indirectly over time, for example via the running time.
Preferably the comparison variable is determined using a statistical test, in particular a so-called Kolmogorov-Smirnov test.
In addition, in particular alternatively, the comparison variable can be a maximum Euclidean distance between the signal of the component to be investigated and the reference signal.
Advantageously, in order to determine the maximum Euclidean distance the operating parameter values are compared with the reference values. The maximum Euclidean distance is preferably the maximum deviation between an operating parameter value at any given instant and a reference value at the same instant.
In addition, the comparison variable can be normalized using the time, one or more of the reference values, a constant and/or another value.
The reference signal, in particular the reference values, can be simulated. Furthermore, the reference signal, in particular the reference values, can be estimated using a hypothetical probability density function.
In a preferred configuration of the invention, the reference signal is determined empirically. This means that it is preferable if the reference values are determined empirically.
The signal of the component to be investigated advantageously describes the functional behavior of the component within a specified time interval. In addition, the reference signal preferably describes an average functional behavior of identical components within the same period of time.
In this way, external, in particular time-dependent, influences such as environmental factors can be taken into account in the same way both for the signal of the component to be investigated and for the reference signal. In other words, in this way a normalization of the external influences can be performed.
In principle, the reference signal can also describe an average functional behavior of identical components within a different time interval, for example, within a previous time interval.
Identical components can be understood to mean components which are of the same type as the component to be tested. The reference signal preferably describes an average functional behavior of multiple components which are the same as the component to be investigated.
Further, it is advantageous if the reference signal is an average over multiple signals of multiple identical components of the same technical installation.
The multiple signals of the multiple identical components can each contain multiple signal values. In particular, a respective reference value can be an average of multiple signal values (in particular, of the multiple signals) at a specific time. In particular, a respective reference value can be a mean value, for example an arithmetic mean, a median or a modal value, of the multiple signal values at the given time. This means that in the averaging the mean value of the multiple signal values is calculated, in particular, at a specific point in time. Furthermore, a respective reference value can be determined using a (different) location parameter, e.g. using a quantile of the multiple signal values, at the specific point in time.
In the averaging of the multiple signals of the multiple identical components, a weighting of the signals can be performed.
If the reference signal is an average derived from the multiple signals of the multiple identical components of the same technical installation, then external factors, in particular, system-dependent factors, such as the effects of a running time of the installation, effects of a service implemented by the installation or the like, can be taken into account in the same way both for the signal of the component to be investigated and for the reference signal. In other words, in this way a normalization of the external influences can be performed. This eliminates the need for modeling external factors.
In principle, the reference signal can also be an average over a plurality of signals of a plurality of identical components of at least one other technical installation.
An identical component can be another component of the same component type, in particular of the same technical installation. Preferably, components of the same component type are identical in design.
Also, an identical component can be an identically responding component of a different component type, in particular, of the same technical installation. The phrase ‘identically responding’ can be understood to mean responding in the same way to external factors. Advantageously, the identically responding component responds to external influences in a similar or identical way as the component to be investigated. In particular, the identically responding component can respond in the same manner as the component to be investigated to a running time of the installation, to a service implemented by the installation, to environmental factors (such as an ambient temperature), etc.
Advantageously, a limit value is defined for the probability of occurrence of the comparison variable, in short: for the probability of occurrence. The limit value can specify the probability of occurrence above which a behavior of the component to be investigated is deemed to be normal, or should be detected as normal. The limit value can be, for example, 5% (i.e. 0.05), 1% (i.e. 0.01), 0.5% (i.e. 0.005) or 0.1% (i.e. 0.001).
If the probability of occurrence of the comparison variable is greater than the specified limit value, a normal behavior of the component to be investigated is detected. A normal behavior can be an expected behavior.
This means that in the (specifiable) distribution of the plurality of such comparison variables, in particular using the limit value, a range can be specified for expected behavior or for normal behavior, in which advantageously a normal behavior of the component to be investigated is detected.
If the probability of occurrence of the comparison variable is less than the specified limit value, an abnormal behavior of the component to be investigated is preferably detected. An abnormal behavior can be an unexpected behavior. For example, an abnormally behaving component may be not fully functional and/or deemed to be not fully functional.
The probability of occurrence of the comparison variable is determined using a specifiable distribution of a plurality of such comparison variables.
Preferably, the plurality of such comparison variables already exists. For example, the distribution of the plurality of such comparison variables may be specified. Further, the specifiable distribution of the plurality of such comparison variables may be determined on the basis of the plurality of such comparison variables.
Further, the plurality of such comparison variables can be determined.
In particular, the plurality of such comparison variables is determined, for example, by allowing a respective comparison variable to be determined for multiple components and/or for multiple time intervals in the same way as previously described. On the basis of the determined comparison variables, the distribution (of the plurality of such comparison variables) can be determined.
The probability of occurrence advantageously depends on the specified distribution of the plurality of such comparison variables. In this way, the probability of occurrence is dependent on the plurality of such comparison variables. Advantageously, the more of these comparison variables that are or have been determined, the more stable the method is. Moreover, the plurality of such comparison variables can be (at least partially) updated. In particular, in the updating additional comparison variables can be included in the set of the plurality of such comparison variables. As a result of the updating of the plurality of such comparison variables, the distribution of the plurality of the comparison variables can be updated. In this way, the method can be self-learning.
Advantageously, the specified distribution of the plurality of such comparison variables is preferably determined by, for multiple time intervals, comparing a signal of the respective time interval, representing the functional behavior of at least one identical component, with a respective reference signal which preferably describes an average functional behavior of identical components within the same respective time interval. For example, for multiple time intervals a signal of the respective time interval representing the functional behavior of the component to be investigated can be compared with a respective reference signal, which preferably describes an average functional behavior of identical components within the same respective time interval.
In each comparison, in particular to determine the distribution of the plurality of such comparison variables, a comparison variable describing the deviation of the signal from the reference signal is preferably determined in each case. On the basis of the multiple comparison variables, the distribution of the plurality of such comparison variables can be determined.
Furthermore, it is expedient if the specified distribution of the plurality of such comparison variables is determined by, for multiple components of the same component type, in particular of the same technical installation and/or at least one other technical installation of the same installation type, a signal of the respective component representing the functional behavior of the respective component being compared with a reference signal, which advantageously describes an average functional behavior of identical components, in particular of the respective installation. In each comparison, in particular to determine the distribution of the plurality of such comparison variables, a comparison variable describing the deviation of the signal from the reference signal is preferably determined in each case. On the basis of the multiple comparison variables, the distribution of the plurality of such comparison variables can be determined.
Furthermore, the invention is directed towards a computer program having commands which, when the computer program is executed on a computer, causes said computer to execute the above-mentioned method and/or one of its refinements.
In addition, the invention is directed towards a computer readable storage medium having commands which, when the computer program is executed by a computer, cause said computer to execute the above-mentioned method and/or one of its refinements.
The description of advantageous embodiments of the invention given so far contains numerous features, some of which are combined together in the individual dependent claims. However, these features can also be advantageously considered individually and combined to form other meaningful combinations. In particular, these features can be combined individually and in any suitable combination with the method according to the invention, the computer program according to the invention and the computer-readable storage medium. Thus, method features are also to be regarded as properties relating to the subject matter of the corresponding computer program and/or the corresponding computer-readable storage medium and vice versa.
Even if in the description or in the claims some of the terms are used in the singular or in conjunction with a numeral, the scope of the invention for these terms shall not be restricted to the singular or the respective numeral.
The properties, features and advantages of the present invention and the manner in which these are achieved will become clearer and more comprehensible in conjunction with the following description of the exemplary embodiment, which is explained in more detail in connection with the drawings. The exemplary embodiment serves as the explanation of the invention and does not limit the invention to the combination of features specified therein, which includes functional features. In addition, appropriate features of any exemplary embodiment can also be explicitly considered in isolation, separately from any exemplary embodiment, incorporated into another exemplary embodiment in order to supplement its provisions, and combined with any one of the claims.
Shown are:
On the x-axis 8 of the graph 2 the time t is plotted. On the y-axis 10 of the graph 2 an operating parameter is plotted.
In this example, the component to be investigated is a wheelset bearing of a rail vehicle. The identical components are other wheelset bearings of the same rail vehicle, in particular all the other wheelset bearings of the same rail vehicle. Each wheelset bearing comprises a sensor which detects values of an operating parameter. The operating parameter is sensitive to damage to the wheelset bearing.
In this example, the operating parameter is a temperature of the respective component, in this case the respective wheelset bearing.
The sensors determine (quasi-)continuous values of the operating parameter, operating parameter values for short, in this example with a frequency of f=1 min−1.
Each wheelset bearing can also comprise a plurality of sensors whose signals are then averaged to form a signal of the respective component. Further, the operating parameter values can be smoothed.
For each wheelset bearing Ri with i=1, 2, . . . N each signal comprises operating parameter values xRi(t) as a function of time t.
In this example, the component R1 (i.e. i=1) is the component to be investigated. The signal 4 of the component to be investigated comprises the operating parameters xRi with i=1.
The reference signal 6 is calculated by calculating the expected curve {circumflex over (X)}Ri(t) of all other wheelset bearings Rj with j≠i. To this end a mean value is determined, here an arithmetic mean, of the operating parameters xRj of the wheelset bearings Rj with j≠i:
{circumflex over (x)}R
In principle, the mean value could also be a median, a modal value or a quantile.
In principle, it is possible that statistical and/or systematic fluctuations of at least one signal can be allowed for using a correction factor for the fluctuating signal. For the sake of better clarity, no correction factors have been introduced here.
As an example of the signal 4 of the component to be investigated,
The signal 4 of the component to be investigated is compared with the reference signal 6.
In the comparison, a comparison variable 14 describing the deviation of the signal 4 from the reference signal 6 is determined.
In
To determine the comparison variable 14 the operating parameter values xRi(t) are accumulated over time. The cumulative operating values XRi are shown in the diagram as the solid line 20.
In addition, the reference values {circumflex over (X)}Ri(t) are accumulated over time. The cumulative reference values {circumflex over (X)}Ri are shown in the diagram 12 as the dashed line 22.
A maximum cumulative deviation di, also the maximum cumulative distance di, between the accumulated operating values XRi and the accumulated reference values {circumflex over (X)}Ri is determined as the comparison variable. The maximum cumulative deviation di is determined using a so-called Kolmogorov-Smirnov test statistic. The maximum cumulative deviation di is normalized to the number of measured values N in the respective time interval, wherein here, for example, N=min(Nv,N{circumflex over (V)}). The comparison variable 14, here the maximum cumulative deviation di, is calculated as follows:
The comparison variable 14, here the maximum cumulative deviation di, is indicated in
On the x-axis 26 of the graph 24 the comparison variable 14, here the maximum cumulative deviation di, is plotted. On the y-axis 28 of the graph 24 a frequency is plotted.
The method described under
The distribution 30 of the plurality of such distribution variables is determined by obtaining the frequencies of these determined distribution variables. In this way, an empirical frequency distribution can be determined. The empirical frequency distribution is shown in the diagram 24 in
Further, in this example to determine the distribution 30 of the plurality of such comparison variables, a distribution function of a given distribution type is defined as a parametrized distribution function. The distribution function can be a distribution density function or a cumulative distribution function. In this example, a distribution density function is specified, for example. For example, the specified distribution type can be a logarithmic normal distribution or an exponential function. The parameters of the parameterized distribution function are fitted (using known methods), so that a modified distribution function is determined. The modified distribution function is shown in the diagram 24 in
The modified distribution function is used to calculate the cumulative distribution function. The cumulative distribution function is shown in the diagram 24 in
A probability of occurrence of the comparison variable di with i=1, in short d1, is determined for the component R1 to be investigated using the specified distribution 30 of the plurality of such comparison variables.
The probability of occurrence is then directly a measure of how anomalous the observed temperature profile is on the wheelset bearing Ri to be investigated, here R1.
For the probability of occurrence a limit value G is specified. In particular, the limit value G is specified before the calculation of the probability of occurrence. The limit value G can be, for example, 1% (0.01), 0.5% (0.005) or 0.1% (0.001).
If the probability of occurrence of the comparison variable d1 is greater than the specified limit value G, a normal behavior of the component to be investigated is detected.
If the probability of occurrence of the comparison variable d1 is less than the specified limit value G, an abnormal behavior of the component to be investigated is detected.
If the probability of occurrence of the comparison variable d1 falls below the limit value G, for example, here e.g. 0.005, an alarm is generated which can be forwarded to a maintenance planning and/or maintenance control system, for example. In addition, an alarm can be generated when the probability of occurrence of the comparison variable d1 falls below the limit value G several times within a predefined time interval.
Although the invention has been illustrated and described in greater detail by means of the preferred exemplary embodiments, the invention is not restricted by the examples disclosed and other variations can be derived therefrom by the person skilled in the art without departing from the scope of protection of the invention.
| Number | Date | Country | Kind |
|---|---|---|---|
| 102017215341.5 | Sep 2017 | DE | national |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2018/070319 | 7/26/2018 | WO |
| Publishing Document | Publishing Date | Country | Kind |
|---|---|---|---|
| WO2019/042671 | 3/7/2019 | WO | A |
| Number | Name | Date | Kind |
|---|---|---|---|
| 6108616 | Borchers et al. | Aug 2000 | A |
| 8261132 | Klein et al. | Sep 2012 | B2 |
| 20070115800 | Fonseka | May 2007 | A1 |
| 20070150220 | Bode | Jun 2007 | A1 |
| 20100189051 | Li et al. | Jul 2010 | A1 |
| 20110137432 | Wang et al. | Jun 2011 | A1 |
| 20140006330 | Biem | Jan 2014 | A1 |
| 20140136137 | Tarshish-Shapir | May 2014 | A1 |
| 20150213706 | Bai et al. | Jul 2015 | A1 |
| 20160140009 | Watanabe | May 2016 | A1 |
| 20170288811 | Tosato | Oct 2017 | A1 |
| 20200241979 | Fankhauser | Jul 2020 | A1 |
| Number | Date | Country |
|---|---|---|
| 3734487 | Apr 1989 | DE |
| 19732046 | Jan 1999 | DE |
| 19961631 | Jun 2001 | DE |
| 102005053648 | May 2007 | DE |
| 102006019896 | Oct 2007 | DE |
| 1542108 | Jun 2005 | EP |
| 2012144068 | Aug 2012 | JP |
| 2364911 | Aug 2009 | RU |
| WO2016191711 | Dec 2016 | WO |
| Number | Date | Country | |
|---|---|---|---|
| 20200241979 A1 | Jul 2020 | US |
