RU2435189C2 - Method of searching for faulty unit in dynamic system - Google Patents
Method of searching for faulty unit in dynamic system Download PDFInfo
- Publication number
- RU2435189C2 RU2435189C2 RU2009123999/08A RU2009123999A RU2435189C2 RU 2435189 C2 RU2435189 C2 RU 2435189C2 RU 2009123999/08 A RU2009123999/08 A RU 2009123999/08A RU 2009123999 A RU2009123999 A RU 2009123999A RU 2435189 C2 RU2435189 C2 RU 2435189C2
- Authority
- RU
- Russia
- Prior art keywords
- integral
- control points
- determined
- dynamic characteristics
- jnom
- Prior art date
Links
Landscapes
- Testing And Monitoring For Control Systems (AREA)
- Complex Calculations (AREA)
Abstract
Description
Изобретение относится к области контроля и диагностирования систем автоматического управления и их элементов.The invention relates to the field of monitoring and diagnosing automatic control systems and their elements.
Известен способ контроля динамического блока в составе системы управления (Патент РФ №2136033, МКИ6 G05B 23/02, 1999), основанный на интегрировании выходного и входного сигналов блока с весом e-αt, где α - вещественная константа.A known method of controlling a dynamic block as part of a control system (RF Patent No. 2136033, MKI 6 G05B 23/02, 1999), based on the integration of the output and input signals of the block with a weight e- αt , where α is a real constant.
Недостатком этого способа является то, что его применение для контроля нескольких блоков системы управления произвольной структуры приводит к необходимости интегрирования входных и выходных сигналов каждого контролируемого блока.The disadvantage of this method is that its use for monitoring several blocks of a control system of an arbitrary structure leads to the need to integrate the input and output signals of each controlled block.
Наиболее близким техническим решением (прототипом) является способ поиска неисправностей в динамических системах (Патент РФ №2199776, МКИ6 G05B 23/02, 2003).The closest technical solution (prototype) is a method for troubleshooting in dynamic systems (RF Patent No. 2199776, MKI 6 G05B 23/02, 2003).
Недостатком этого способа является то, что он предполагает определение функций чувствительности, что является достаточно трудоемкой задачей, связанной со значительными временными и аппаратными затратами.The disadvantage of this method is that it involves the determination of sensitivity functions, which is a rather time-consuming task associated with significant time and hardware costs.
Технической задачей, на решение которой направлено данное изобретение, является уменьшение аппаратных и временных затрат на диагностирование путем замены значений функций чувствительности на величины деформаций интегральных преобразований динамических характеристик системы, полученных в результате пробных отклонений параметров соответствующих блоков, а также обеспечение возможности сравнения результатов различных режимов диагностирования для выбора наилучшего в плане различимости дефектов путем применения нормированных диагностических признаков.The technical problem to which this invention is directed is to reduce the hardware and time costs for diagnosing by replacing the values of the sensitivity functions with the strain values of the integral transformations of the dynamic characteristics of the system obtained as a result of trial deviations of the parameters of the corresponding blocks, as well as providing the possibility of comparing the results of various diagnostic modes to select the best in terms of distinguishability of defects by applying normalized diagnoses ostichnyh signs.
Поставленная задача достигается тем, что регистрируют реакцию заведомо исправной системы fjном(t), j=1, 2, …, k на интервале t∈[0, TK] в k контрольных точках и определяют интегральные оценки выходных сигналов Fjном(α), j=1, …, k, системы, для чего в момент подачи тестового сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование сигналов системы управления в каждой из k контрольных точек с весами e-αt, где , путем подачи на первые входы k блоков перемножения сигналов системы управления, на вторые входы блоков перемножения подают экспоненциальный сигнал e-αt, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Tк, полученные в результате интегрирования оценки выходных сигналов Fjном(α), j=1, …, k регистрируют, определяют интегральные преобразования динамических характеристик модели для каждой из k контрольных точек, полученные в результате пробных отклонений параметров каждого из m блоков, для чего поочередно для каждого блока динамической системы вводят пробное отклонение параметра его передаточной функции и находят интегральные оценки выходных сигналов системы для параметра α и тестового сигнала x(t), полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений Pji(α), j=1, …, k; i=1, …, m регистрируют, определяют деформации интегральных преобразований динамических характеристик модели, полученные в результате пробных отклонений параметров соответствующих блоков ΔPji(α)=Pji(α)-Fjном(α), j=1, …, k; i=1, …, m, определяют нормированные значения деформаций интегральных преобразований динамических характеристик модели, полученные в результате пробных отклонений параметров соответствующих блоков из соотношения , замещают систему с номинальными характеристиками контролируемой, на вход системы подают аналогичный тестовый сигнал x(t), определяют интегральные преобразования динамических характеристик контролируемой системы для k контрольных точек Fj(α), j=1, …, k для параметра α, определяют отклонения интегральных преобразований динамических характеристик контролируемой системы для k контрольных точек от номинальных значении ΔFj(α)=Fj(α)-Fjном(α), j=1, …, k, определяют нормированные значения отклонений интегральных преобразований динамических характеристик контролируемой системы из соотношения ,The problem is achieved by registering the reaction of a known-good system f jnom (t), j = 1, 2, ..., k on the interval t∈ [0, T K ] at k control points and determining the integral estimates of the output signals F jnom (α ), j = 1, ..., k, of the system, for which, at the time the test signal is input to the input of the system with nominal characteristics, the integration of the control system signals at each of k control points with weights e -αt starts at the same time, where , by applying the control system signals to the first inputs k of the multiplication blocks, the exponential signal e -αt is supplied to the second inputs of the multiplication blocks, the output signals of the multiplication blocks are fed to the inputs of the integration blocks k, the integration is completed at time T k , obtained as a result of integration of the estimate jnom output signals F (α), j = 1, ..., k register define integral transform dynamic characteristics model for each of the control points k resulting test parameter deviations ka of m blocks, for which, for each block of the dynamical system, a test deviation of the parameter of its transfer function is introduced and integral estimates of the system output signals for parameter α and test signal x (t) are obtained by integrating the estimate of the output signals for each of the k control points and each of m test deviations P ji (α), j = 1, ..., k; i = 1, ..., m are recorded, the deformations of the integral transformations of the dynamic characteristics of the model are determined, obtained as a result of trial deviations of the parameters of the corresponding blocks ΔP ji (α) = P ji (α) -F jnom (α), j = 1, ..., k ; i = 1, ..., m, determine the normalized strain values of the integral transforms of the dynamic characteristics of the model, obtained as a result of trial deviations of the parameters of the corresponding blocks from the relation , replace the system with the nominal characteristics of the controlled one, apply the same test signal x (t) to the input of the system, determine the integral transformations of the dynamic characteristics of the controlled system for k control points F j (α), j = 1, ..., k for the parameter α, determine the deviations the integral transformations of the dynamic characteristics of the controlled system for k control points from the nominal value ΔF j (α) = F j (α) -F jnom (α), j = 1, ..., k, determine the normalized values of the deviations of the integral transformations of the dynamic characteristics ristik of a controlled system from the relation ,
определяют диагностические признаки из соотношенияdetermine diagnostic signs from the ratio
, i=1, …, m, по минимуму диагностического признака определяют неисправный блок. , i = 1, ..., m, at the minimum of a diagnostic sign, the faulty block is determined.
Сущность предлагаемого способа заключается в следующем.The essence of the proposed method is as follows.
Способ основан на использовании пробных отклонений параметров модели непрерывной динамической системы. Для получения диагностических признаков динамических элементов используются преобразования по Лапласу временных функций:The method is based on the use of trial deviations of the model parameters of a continuous dynamic system. To obtain diagnostic signs of dynamic elements, Laplace transforms of time functions are used:
в области вещественных значений переменной Лапласа p=α, в интервале 0≤α≤∞. Использование преобразования Лапласа позволяет перейти от обработки временных функций к анализу численных значений их функционалов.in the range of real values of the Laplace variable p = α, in the interval 0≤α≤∞. Using the Laplace transform allows you to switch from processing time functions to analyzing the numerical values of their functionals.
Интегральные преобразования находят на временном интервале Tk в k контрольных точках:Integral transformations are found on the time interval T k at k control points:
Прототип (Патент РФ №2199776, МКИ6 G05B 23/02, 2003) предполагает, что для каждого динамического элемента с номером i вычисляется диагностический признак по формуле:The prototype (RF Patent No. 2199776, MKI 6 G05B 23/02, 2003) suggests that for each dynamic element with number i, a diagnostic symptom is calculated by the formula:
где ΔF(α)=(ΔFl(α), ΔF2(α), …, ΔFk(α))T - вектор изображений для вещественного значения переменной Лапласа α отклонений временных характеристик объекта в k контрольных точках;where ΔF (α) = (ΔF l (α), ΔF 2 (α), ..., ΔF k (α)) T is the vector of images for the real value of the Laplace variable α of the deviations of the temporal characteristics of the object at k control points;
- структурная чувствительность для j-й контрольной точки, i-го динамического элемента и значения переменной Лапласа α. - structural sensitivity for the j-th control point, the i-th dynamic element and the value of the Laplace variable α.
По минимуму значения диагностического признака (3) выносят решение о наличии дефекта в динамическом элементе.At a minimum, the values of the diagnostic sign (3) make a decision on the presence of a defect in the dynamic element.
Диагностические признаки одиночных структурных дефектов вида (3) принимают численные значения, которые зависят как от степени отклонения динамических характеристик, так и от свойств модели чувствительности и могут принимать в общем случае любые неотрицательные значения. Для практических целей удобно иметь диагностические признаки, область значений которых была бы ограничена, в частности, интервалом [0, 1], где крайние значения могли бы интерпретироваться как однозначное наличие (отсутствие) дефекта, а промежуточные значения - как вероятность наличия дефекта. Приведение диагностических признаков к единой шкале открывает возможности для сравнительного анализа условий и результатов диагностирования различных объектов диагностирования в терминах количественной различимости дефектов.Diagnostic signs of single structural defects of the form (3) take numerical values, which depend both on the degree of deviation of the dynamic characteristics and on the properties of the sensitivity model and can generally take any non-negative values. For practical purposes, it is convenient to have diagnostic features, the range of which would be limited, in particular, to the interval [0, 1], where extreme values could be interpreted as the unambiguous presence (absence) of a defect, and intermediate values as the probability of a defect. Bringing diagnostic features to a single scale opens up possibilities for a comparative analysis of the conditions and results of diagnosing various diagnostic objects in terms of the quantitative distinguishability of defects.
Возможность получения таких нормированных диагностических признаков дает векторная их интерпретация, рассмотренная ниже. Идея нормировки признака заключается в приведении векторов деформации динамических характеристик к единичной длине.The possibility of obtaining such normalized diagnostic signs provides their vector interpretation, discussed below. The idea of normalizing a sign is to reduce the strain vectors of dynamic characteristics to a unit length.
Произведем нормирование значений коэффициентов структурной чувствительности по всем контрольным точкам:We will normalize the values of the coefficients of structural sensitivity for all control points:
Подставляя эти значения в выражение (3), получим:Substituting these values in expression (3), we obtain:
Коэффициенты структурной чувствительности (4) могут рассматриваться как координаты вектора единичной длины в k-мерном пространстве:The structural sensitivity coefficients (4) can be considered as the coordinates of the vector unit length in k-dimensional space:
, ,
а отклонение изображений реакций объекта - как вектор в пространстве той же размерности с координатами:and the deviation of the images of the reactions of the object - as a vector in space of the same dimension with coordinates:
. .
Тогда формула (5) с учетом введенных векторов запишется в виде:Then formula (5), taking into account the introduced vectors, can be written in the form:
где - скалярное произведение вектора отклонений динамической характеристики объекта на нормированный вектор структурной чувствительности по i-му динамическому элементу.Where - the scalar product of the vector of deviations of the dynamic characteristic of the object by the normalized vector of structural sensitivity with respect to the ith dynamic element.
По определению скалярного произведения:By definition of a scalar product:
, ,
где - означает длину вектора;Where - means the length of the vector;
φ(α) - угол в k-мерном пространстве между этими векторами.φ (α) is the angle in the k-dimensional space between these vectors.
Если направление векторов ΔF(α) и совпадают или противоположны, то cosφ(α)=±1 и , подставляя эти значения в формулу (6), получим:If the direction of the vectors ΔF (α) and are the same or opposite, then cosφ (α) = ± 1 and Substituting these values in the formula (6), we obtain:
. .
Таким образом, в терминах векторной интерпретации поиск одиночного структурного дефекта заключается в подборе такого индекса i, для которого нормированный вектор в k-мерном пространстве в наибольшей степени совпадает или противоположен с направлением соответствующего вектора ΔF(α) деформации динамических характеристик объекта.Thus, in terms of vector interpretation, the search for a single structural defect consists in selecting an index i for which the normalized vector in k-dimensional space to the greatest extent coincides or is opposite to the direction of the corresponding vector ΔF (α) of deformation of the dynamic characteristics of the object.
Сформируем нормированный диагностический признак путем деления каждого слагаемого выражения (5) на величину :We form a normalized diagnostic sign by dividing each term of expression (5) by :
где - элементы нормированного вектора деформации динамической характеристики объекта.Where - elements of the normalized deformation vector of the dynamic characteristics of the object.
Используя векторную интерпретацию выражения (7), запишем его в следующем видеUsing a vector interpretation of expression (7), we write it in the following form
где - угол между векторами единичной длины отклонений динамической характеристики объекта и структурной чувствительности для i-го динамического элемента и параметра интегрального преобразования α.Where - the angle between the unit length vectors of the deviations of the dynamic characteristics of the object and structural sensitivity for the i-th dynamic element and the integral transformation parameter α.
Таким образом, нормированный диагностический признак (7) представляет собой значение квадрата синуса угла, образованного в k-мерном пространстве нормированными векторами структурной чувствительности динамического элемента и деформации интегрального преобразования динамической характеристики объекта диагностирования.Thus, the normalized diagnostic sign (7) represents the value of the square of the sine of the angle formed in the k-dimensional space by the normalized vectors of the structural sensitivity of the dynamic element and the strain of the integral transformation of the dynamic characteristic of the diagnostic object.
Получение модели структурной чувствительности является отдельной, достаточно сложной задачей. Для уменьшения объема требуемых вычислений при получении диагностических признаков заменим в выражении (7) нормированные векторы интегральных преобразований модели структурной чувствительности на нормированные векторы деформаций интегральных преобразований динамических характеристик модели, полученные в результате пробных отклонений параметров соответствующих блоков:Obtaining a model of structural sensitivity is a separate, rather difficult task. To reduce the amount of required calculations when obtaining diagnostic features, we replace the normalized vectors in expression (7) integral transformations of the model of structural sensitivity to normalized vectors deformations of integral transformations of the dynamic characteristics of the model obtained as a result of trial deviations of the parameters of the corresponding blocks:
Выражение (8) позволяет получить значения нормированного диагностического признака, представляющего собой значение квадрата синуса угла, образованного вектором реальной деформации интегральной оценки динамической характеристики, обусловленной наличием дефекта в объекте, и вектором деформации интегральной оценки динамической характеристики, полученной в результате пробного отклонения параметров блока объекта. Пробное отклонение блока, минимизирующее значение диагностического признака (8), указывает на наличие дефекта в этом блоке. Область возможных значений признака лежит в интервале [0, 1].Expression (8) allows one to obtain the values of the normalized diagnostic feature, which is the value of the square of the sine of the angle formed by the vector of real deformation of the integral estimate of the dynamic characteristic due to the presence of a defect in the object, and the deformation vector of the integral estimate of the dynamic characteristic obtained as a result of a test deviation of the parameters of the block of the object. A test deviation of the block, minimizing the value of the diagnostic sign (8), indicates the presence of a defect in this block. The range of possible values of the attribute lies in the interval [0, 1].
Таким образом, предлагаемый способ поиска неисправностей сводится к выполнению следующих операций.Thus, the proposed troubleshooting method is reduced to the following operations.
1. В качестве динамической системы рассматривают систему, состоящую из произвольно соединенных m динамических элементов.1. As a dynamic system, consider a system consisting of arbitrarily connected m dynamic elements.
2. Предварительно определяют время контроля ТК≥ТПП, где ТПП - время переходного процесса системы. Время переходного процесса оценивают для номинальных значений параметров динамической системы.2. Pre-determine the control time T To ≥T PP , where T PP - the transition process of the system. The transient time is estimated for the nominal values of the parameters of the dynamic system.
3. Определяют параметр интегрального преобразования сигналов из соотношения .3. Determine the parameter of the integral signal conversion from the ratio .
4. Фиксируют число контрольных точек k.4. Fix the number of control points k.
5. Предварительно определяют нормированные векторы деформаций интегральных преобразований динамических характеристик модели, полученные в результате пробных отклонений параметров i-го блока каждого из m блоков для номинальных значений параметров передаточных функций блоков и определенного выше параметра α, для чего выполняют пункты 6-10.5. Predefined normalized vectors deformations of integral transformations of the dynamic characteristics of the model obtained as a result of test deviations of the parameters of the i-th block of each of m blocks for the nominal values of the parameters of the transfer functions of the blocks and the parameter α defined above, for which points 6-10 are performed.
6. Подают тестовый сигнал x(t) (единичный ступенчатый, линейно возрастающий, прямоугольный импульсный и т.д.) на вход системы управления с номинальными характеристиками. Принципиальных ограничений на вид входного тестового воздействия предлагаемый способ не предусматривает.6. The test signal x (t) (unit step, linearly increasing, rectangular pulse, etc.) is supplied to the input of the control system with nominal characteristics. The proposed method does not provide fundamental restrictions on the type of input test exposure.
7. Регистрируют реакцию системы fjном(t), j=1, 2, …, k на интервале t∈[0, TK] в k контрольных точках и определяют интегральные оценки выходных сигналов Fjном(α), j=1, …, k системы. Для этого в момент подачи тестового сигнала на вход системы управления с номинальными характеристиками одновременно начинают интегрирование сигналов системы управления в каждой из k контрольных точек с весами e-αt, где , для чего сигналы системы управления подают на первые входы к блоков перемножения, на вторые входы блоков перемножения подают экспоненциальный сигнал e-αt, выходные сигналы к блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Тк, полученные в результате интегрирования оценки выходных сигналов Fjном(α), j=1, …, k регистрируют.7. The reaction of the system f jnom (t), j = 1, 2, ..., k on the interval t∈ [0, T K ] at k control points is recorded and the integral estimates of the output signals F jnom (α), j = 1, ..., k systems. To do this, at the time of supplying a test signal to the input of the control system with nominal characteristics, the integration of the control system signals at each of k control points with weights e -αt starts simultaneously why the control system signals are fed to the first inputs to the multiplication units, the exponential signal e -αt is supplied to the second inputs of the multiplying units, the output signals to the multiplying units are fed to the inputs of the integration blocks k, the integration is completed at time T to , obtained as a result of integration estimates of the output signals F jnom (α), j = 1, ..., k are recorded.
8. Определяют интегральные преобразования динамических характеристик модели для каждой из k контрольных точек, полученные в результате пробных отклонений параметров каждого из m блоков, для чего поочередно для каждого блока динамической системы вводят пробное отклонение параметра передаточной функции и выполняют пункты 6 и 7 для одного и того же тестового сигнала x(t). Полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений Pji(α), j=1, …, k; i=1, …, m регистрируют.8. The integral transformations of the dynamic characteristics of the model for each of the k control points are determined, obtained as a result of test deviations of the parameters of each of the m blocks, for which, for each block of the dynamic system, a test deviation of the transfer function parameter is introduced and steps 6 and 7 are performed for the same same test signal x (t). Estimates of the output signals obtained as a result of integration for each of k control points and each of m test deviations P ji (α), j = 1, ..., k; i = 1, ..., m are recorded.
9. Определяют деформации интегральных преобразований динамических характеристик модели, полученные в результате пробных отклонений параметров соответствующих блоков ΔPji(α)=Pji(α)-Fjном(α), j=1, …, k; i=1, …, m.9. The deformations of the integral transformations of the dynamic characteristics of the model are determined, obtained as a result of trial deviations of the parameters of the corresponding blocks ΔP ji (α) = P ji (α) -F jnom (α), j = 1, ..., k; i = 1, ..., m.
10. Определяют нормированные значения деформаций интегральных преобразований динамических характеристик модели, полученные в результате пробных отклонений параметров соответствующих блоков по формуле .10. Determine the normalized strain values of the integral transforms of the dynamic characteristics of the model, obtained as a result of trial deviations of the parameters of the corresponding blocks by the formula .
11. Замещают систему с номинальными характеристиками контролируемой. На вход системы подают аналогичный тестовый сигнал x(t).11. Replace the system with the rated characteristics controlled. A similar test signal x (t) is supplied to the system input.
12. Определяют интегральные преобразования динамических характеристик контролируемой системы для k контрольных точек Fj(α), j=1, …, k, осуществляя операции, описанные в пунктах 6 и 7 применительно к контролируемой системе.12. Determine the integral transformations of the dynamic characteristics of the controlled system for k control points F j (α), j = 1, ..., k, by performing the operations described in clauses 6 and 7 as applied to the controlled system.
13. Определяют отклонения интегральных преобразований динамических характеристик контролируемой системы для k контрольных точек от номинальных значений ΔFj(α)=Fjном(α), j=1, …, k.13. The deviations of the integral transformations of the dynamic characteristics of the controlled system for k control points from the nominal values ΔF j (α) = F jnom (α), j = 1, ..., k are determined .
14. Вычисляют нормированные значения отклонений интегральных преобразований динамических характеристик контролируемой системы по формуле .14. The normalized values of the deviations of the integral transforms of the dynamic characteristics of the controlled system are calculated by the formula .
15. Вычисляют диагностические признаки наличия неисправного блока по формуле (8).15. Calculate the diagnostic signs of a faulty unit by the formula (8).
16. По минимуму значения диагностического признака определяют дефектный блок.16. At the minimum, the values of the diagnostic sign determine the defective block.
Поскольку диагностические признаки (8) имеют область возможных значений, ограниченную интервалом [0, 1], то разность между ближайшим к минимальному признаку и минимальным признаком (который указывает на дефектный блок) количественно характеризует различимость данного дефекта с учетом расположения блока на структурной схеме, вида и параметров передаточных функций блоков и всех условий диагностирования, при которых получены эти значения диагностических признаков (вид тестового сигнала, величина параметра α, количество и расположение контрольных точек, величина интервала Тк). Наилучшая различимость - когда указанная разность равна единице (в терминах векторной интерпретации нормированные векторы деформаций интегральных преобразований динамических характеристик этих блоков для пробных отклонений ортогональны). Наихудшая различимость - когда указанная разность равна нулю (в терминах векторной интерпретации нормированные векторы деформаций интегральных преобразований динамических характеристик этих блоков для пробных отклонений коллинеарные). Поэтому применение нормированных диагностических признаков позволяет сравнивать результаты диагностирования для выбора оптимальных режимов поиска дефектов.Since the diagnostic signs (8) have a range of possible values limited by the interval [0, 1], the difference between the closest to the minimum sign and the minimum sign (which indicates a defective block) quantitatively characterizes the distinguishability of this defect, taking into account the location of the block on the structural diagram, of the form and parameters of the transfer functions of the blocks and all the diagnostic conditions under which these values of diagnostic signs are obtained (type of test signal, value of parameter α, number and location of role points, the value of the interval T to ). The best distinguishability is when the indicated difference is equal to unity (in terms of vector interpretation, the normalized strain vectors of the integral transforms of the dynamic characteristics of these blocks for trial deviations are orthogonal). The worst distinguishability is when the indicated difference is equal to zero (in terms of vector interpretation, the normalized strain vectors of the integral transforms of the dynamic characteristics of these blocks are collinear for trial deviations). Therefore, the use of normalized diagnostic features allows you to compare the diagnostic results to select the optimal defect search modes.
Рассмотрим реализацию предлагаемого способа поиска одиночного дефекта для системы, структурная схема которой представлена на чертеже.Consider the implementation of the proposed method for finding a single defect for a system whose structural diagram is shown in the drawing.
Передаточные функции блоков:Transfer functions of blocks:
; : , ; : ,
номинальные значения параметров: T1=5 с; К1=1; К2=1; Т2=1 с; К3=1; Т3=5 с. При поиске одиночного дефекта в виде отклонения постоянной времени T1=4 с в первом звене путем подачи ступенчатого тестового входного сигнала единичной амплитуды и интегрального преобразования сигналов для параметра α=0.5 и Тк=10 с получены значения диагностических признаков по формуле (7) при использовании трех контрольных точек, расположенных на выходах блоков, и модели чувствительности: J1=0; J2=0.186; J3=0.018. Минимальное значение признака J1 однозначно указывает на наличие дефекта в первом блоке, а разность между третьим и первым признаками может количественно характеризовать фактическую различимость этого дефекта. Тот же дефект, найденный путем получения пробных отклонений на величину 5% и вычислений по формуле (8), дает следующие значения диагностических признаков: J1=0; J2=0.78; J3=0.074. Анализ значений диагностических признаков показывает, что значения второго и третьего признака, полученные при использовании пробных отклонений, больше, чем при использовании модели чувствительности. Это позволяет сделать вывод, что фактическая различимость дефекта первого блока выше при использовании заявляемого способа. Различимости дефектов второго и третьего блоков при поиске их с использованием пробных отклонений также не хуже, чем при использовании модели чувствительности.nominal values of parameters: T 1 = 5 s; K 1 = 1; K 2 = 1; T 2 = 1 s; K 3 = 1; T 3 = 5 s. When searching for a single defect in the form of a deviation of the time constant T 1 = 4 s in the first link by applying a step test input signal of unit amplitude and integral signal conversion for parameter α = 0.5 and T k = 10 s, the values of diagnostic signs were obtained using formula (7) for using three control points located at the outputs of the blocks, and a sensitivity model: J 1 = 0; J 2 = 0.186; J 3 = 0.018. The minimum value of the sign J 1 unambiguously indicates the presence of a defect in the first block, and the difference between the third and first signs can quantitatively characterize the actual distinguishability of this defect. The same defect found by obtaining test deviations of 5% and calculations by formula (8) gives the following values of diagnostic signs: J 1 = 0; J 2 = 0.78; J 3 = 0.074. An analysis of the values of diagnostic signs shows that the values of the second and third signs obtained using trial deviations are greater than when using the sensitivity model. This allows us to conclude that the actual distinguishability of the defect of the first block is higher when using the proposed method. The distinguishability of defects of the second and third blocks when searching for them using trial deviations is also no worse than when using the sensitivity model.
Поиск неисправного блока согласно предлагаемому способу сводится к выполнению следующих операций.Search for a faulty block according to the proposed method is reduced to performing the following operations.
1. Фиксируют число динамических элементов m=3.1. Fix the number of dynamic elements m = 3.
2. Путем анализа графиков номинальных переходных характеристик определяем время переходного процесса системы. Для данного примера время переходного процесса составляет ТПП=8 с. Фиксируем время контроля Тk≥ТПП. Для данного примера фиксируем Тk=10 с.2. By analyzing the graphs of the nominal transient characteristics, we determine the time of the transition process of the system. For this example, the transient time is T PP = 8 s. We fix the control time T k ≥T PP. For this example, we fix T k = 10 s.
3. Определяем параметр интегрального преобразования сигналов . Для данного примера α=0.5.3. We determine the parameter of integrated signal conversion . For this example, α = 0.5.
4. Фиксируем контрольные точки на выходах блоков: k=3.4. We fix the control points at the outputs of the blocks: k = 3.
5. Предварительно находим элементы векторов ΔP1(α), ΔP2(α), ΔP3(α) деформаций интегральных преобразований динамических характеристик модели, полученные в результате пробных отклонений параметров соответствующих блоков. Величину пробных отклонений выбираем равной 5%:5. Preliminarily we find the elements of the vectors ΔP 1 (α), ΔP 2 (α), ΔP 3 (α) of strains of the integral transformations of the dynamic characteristics of the model, obtained as a result of trial deviations of the parameters of the corresponding blocks. The value of the test deviations is chosen equal to 5%:
6. Находим нормированные векторы деформаций интегральных преобразований динамических характеристик модели, полученные в результате пробных отклонений параметров соответствующих блоков:6. Find normalized vectors deformations of integral transformations of the dynamic characteristics of the model obtained as a result of trial deviations of the parameters of the corresponding blocks:
7. Замещаем систему с номинальными характеристиками контролируемой, в которой введено отклонение параметра Т1 первого блока от номинального на 20%. На вход системы подаем аналогичный тестовый сигнал x(t).7. We replace the system with the nominal characteristics of the controlled one, in which the deviation of the parameter T 1 of the first block from the nominal by 20% is introduced. At the input of the system we apply a similar test signal x (t).
8. Определяем отклонения интегральных преобразований динамических характеристик контролируемой системы для трех контрольных точек от номинальных значений ΔFj(α)=F1(α)-Fjном(α), j=1, 2, 38. We determine the deviations of the integral transformations of the dynamic characteristics of the controlled system for three control points from the nominal values ΔF j (α) = F 1 (α) -F jnom (α), j = 1, 2, 3
ΔF1(α)=0.4003, ΔF2(α)=0.2669, ΔF3(α)=0.07618.ΔF 1 (α) = 0.4003, ΔF 2 (α) = 0.2669, ΔF 3 (α) = 0.07618.
9. Вычисляем нормированные значения отклонений интегральных преобразований динамических характеристик контролируемой системы 9. We calculate the normalized values of the deviations of the integral transforms of the dynamic characteristics of the controlled system
10. Вычисляем диагностические признаки наличия неисправного блока по формуле (8): J1=0, J2=0.7829, j3=0.07399.10. We calculate the diagnostic signs of the presence of a faulty unit by the formula (8): J 1 = 0, J 2 = 0.7829, j 3 = 0.07399.
11. По минимуму значения диагностического признака определяют дефектный блок (в данном случае - №1).11. At a minimum, the values of the diagnostic sign determine the defective block (in this case, No. 1).
Моделирование процессов поиска дефектов во втором и третьем блоках для данного объекта диагностирования при том же параметре α и при единичном ступенчатом входном сигнале дает следующие значения диагностических признаков.Simulation of the defect search processes in the second and third blocks for a given diagnostic object with the same parameter α and with a single step input signal gives the following values of diagnostic signs.
При наличии дефекта в блоке №2 (в виде уменьшения параметра Т2 на 20%): J1=0.7829, J2=0, J3=0.7462.If there is a defect in block No. 2 (in the form of a decrease in the parameter T 2 by 20%): J 1 = 0.7829, J 2 = 0, J 3 = 0.7462.
При наличии дефекта в блоке №3 (в виде уменьшения параметра Т3 на 20%): J1=0.07404, J2=0.7464, J3=0.If there is a defect in block No. 3 (in the form of a decrease in the parameter T 3 by 20%): J 1 = 0.07404, J 2 = 0.7464, J 3 = 0.
Минимальное значение диагностического признака в обоих случаях правильно указывает на дефектный блок.The minimum value of the diagnostic sign in both cases correctly indicates a defective block.
Claims (1)
определяют диагностические признаки из соотношения i=1, …, m, по минимуму диагностического признака определяют неисправный блок. A method for finding a faulty block in a dynamic system, based on the fact that the number m of dynamic elements included in the system is fixed, the monitoring time T K > T PP is determined, the integral signal conversion parameter is determined from the relation use a test signal on the interval t∈ [0, T K ], use the integral estimates obtained for real values of the Laplace variable as dynamic characteristics of the system, fix the number k of control points of the system, record the reaction of the diagnostic object and model, determine the diagnostic sign, at the minimum of a diagnostic sign, a faulty block is determined, characterized in that the reaction of a known-good system F jnom (t), j = 1, 2, ..., k on the interval t∈ [0, T K ] at k control points is recorded and integral estimates of the output signals F jnom (α), j = 1, ..., k, for which, at the time the test signal is input to the input of the system with nominal characteristics, the integration of the control system signals at each of k control points with weights e -αt begins , where by applying the control system signals to the first inputs k of the multiplication blocks, the exponential signal e -t is supplied to the second inputs of the multiplying blocks, the output signals of the multiplying blocks are fed to the inputs of the integration blocks k, the integration is completed at time T K , obtained by integrating the output estimates signals F jnom (α), j = 1, ..., k are recorded, integral transformations of the dynamic characteristics of the model for each of k control points are determined, obtained as a result of test deviations of the parameters of each of m blocks, for which, for each block of the dynamical system, a test deviation of the parameter of its transfer function is introduced and integral estimates of the system output signals for parameter α and test signal x (t) are obtained by integrating the estimate of the output signals for each of the k control points and each of m test deviations P ji (α), j = 1, ..., k; i = 1, ..., m, register, determine the deformations of the integral transformations of the dynamic characteristics of the model, obtained as a result of test deviations of the parameters of the corresponding blocks ΔP ji (α) = P ji (α) -F jnom (α), j = 1, ..., k; i = 1, ..., m, determine the normalized strain values of the integral transforms of the dynamic characteristics of the model, obtained as a result of trial deviations of the parameters of the corresponding blocks from the relation , replace the system with the nominal characteristics of the controlled one, apply the same test signal x (t) to the input of the system, determine the integral transformations of the dynamic characteristics of the controlled system for k control points F j (α), j = 1, ..., k for the parameter α, determine the deviations integral transformations of the dynamic characteristics of the controlled system for k control points from the nominal values ΔF j (α) = F j (α) -F jnom (α), j = 1, ..., k, determine the normalized values of the deviations of the integral transformations of the dynamic characteristics ristik of a controlled system from the relation
determine diagnostic signs from the ratio i = 1, ..., m, at the minimum of a diagnostic sign, the faulty block is determined.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2009123999/08A RU2435189C2 (en) | 2009-06-23 | 2009-06-23 | Method of searching for faulty unit in dynamic system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2009123999/08A RU2435189C2 (en) | 2009-06-23 | 2009-06-23 | Method of searching for faulty unit in dynamic system |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2009123999A RU2009123999A (en) | 2010-12-27 |
RU2435189C2 true RU2435189C2 (en) | 2011-11-27 |
Family
ID=44055421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2009123999/08A RU2435189C2 (en) | 2009-06-23 | 2009-06-23 | Method of searching for faulty unit in dynamic system |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2435189C2 (en) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2473105C1 (en) * | 2011-12-14 | 2013-01-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of detecting faults in units in continuous dynamic system |
RU2506622C1 (en) * | 2012-10-08 | 2014-02-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of finding faulty units in discrete dynamic system |
RU2506623C1 (en) * | 2012-10-08 | 2014-02-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty unit in discrete dynamic system |
RU2513504C1 (en) * | 2013-05-07 | 2014-04-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty unit in continuous dynamic system |
RU2519435C1 (en) * | 2012-12-10 | 2014-06-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of finding faulty units in continuous dynamic system |
RU2528135C1 (en) * | 2013-10-01 | 2014-09-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty unit in continuous dynamic system based on change of position of input signal |
RU2538948C2 (en) * | 2013-02-14 | 2015-01-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty unit in continuous dynamic system |
RU2541857C1 (en) * | 2013-11-06 | 2015-02-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of finding faults in continuous dynamic system based on input of sample deviations |
RU2613402C1 (en) * | 2016-04-12 | 2017-03-16 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Search method of topological defects in continuous dynamic system based on sensitivity functions |
RU2613630C1 (en) * | 2016-03-09 | 2017-03-21 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method for searching faulty unit in continuous dynamic system through introduction of trial deviations |
RU2616512C1 (en) * | 2016-04-20 | 2017-04-17 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method of searching topological defect in continuous dynamic system based on trial deviations introduction |
RU2656923C1 (en) * | 2017-03-10 | 2018-06-07 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method of troubleshooting in continuous dynamic system based on introduction of trial deviations |
RU2676365C1 (en) * | 2018-02-28 | 2018-12-28 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method of searching faulty unit in continuous dynamic system based on introduction of trial deviations |
RU2711000C1 (en) * | 2019-05-27 | 2020-01-14 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method of searching for a faulty unit in a continuous dynamic system based on introduction of trial deviations and analysis of transmission signs |
RU2762532C1 (en) * | 2021-02-11 | 2021-12-21 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method for searching faulty block in continuous dynamical system |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2444774C1 (en) * | 2011-01-13 | 2012-03-10 | Государственное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty module in discrete dynamic system |
RU2451319C1 (en) * | 2011-07-15 | 2012-05-20 | Государственное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty module in dynamic system |
RU2461861C1 (en) * | 2011-10-04 | 2012-09-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty module in continuous dynamic system |
RU2464616C1 (en) * | 2011-11-01 | 2012-10-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faults in dynamic unit in continuous system |
-
2009
- 2009-06-23 RU RU2009123999/08A patent/RU2435189C2/en not_active IP Right Cessation
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2473105C1 (en) * | 2011-12-14 | 2013-01-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of detecting faults in units in continuous dynamic system |
RU2506622C1 (en) * | 2012-10-08 | 2014-02-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of finding faulty units in discrete dynamic system |
RU2506623C1 (en) * | 2012-10-08 | 2014-02-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty unit in discrete dynamic system |
RU2519435C1 (en) * | 2012-12-10 | 2014-06-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of finding faulty units in continuous dynamic system |
RU2538948C2 (en) * | 2013-02-14 | 2015-01-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty unit in continuous dynamic system |
RU2513504C1 (en) * | 2013-05-07 | 2014-04-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty unit in continuous dynamic system |
RU2528135C1 (en) * | 2013-10-01 | 2014-09-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of searching for faulty unit in continuous dynamic system based on change of position of input signal |
RU2541857C1 (en) * | 2013-11-06 | 2015-02-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" | Method of finding faults in continuous dynamic system based on input of sample deviations |
RU2613630C1 (en) * | 2016-03-09 | 2017-03-21 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method for searching faulty unit in continuous dynamic system through introduction of trial deviations |
RU2613402C1 (en) * | 2016-04-12 | 2017-03-16 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Search method of topological defects in continuous dynamic system based on sensitivity functions |
RU2616512C1 (en) * | 2016-04-20 | 2017-04-17 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method of searching topological defect in continuous dynamic system based on trial deviations introduction |
RU2656923C1 (en) * | 2017-03-10 | 2018-06-07 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method of troubleshooting in continuous dynamic system based on introduction of trial deviations |
RU2676365C1 (en) * | 2018-02-28 | 2018-12-28 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method of searching faulty unit in continuous dynamic system based on introduction of trial deviations |
RU2711000C1 (en) * | 2019-05-27 | 2020-01-14 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method of searching for a faulty unit in a continuous dynamic system based on introduction of trial deviations and analysis of transmission signs |
RU2762532C1 (en) * | 2021-02-11 | 2021-12-21 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" | Method for searching faulty block in continuous dynamical system |
Also Published As
Publication number | Publication date |
---|---|
RU2009123999A (en) | 2010-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2435189C2 (en) | Method of searching for faulty unit in dynamic system | |
RU2429518C1 (en) | Method of fault finding in continuous system dynamic unit | |
RU2439648C1 (en) | Method to search for faulty block in dynamic system | |
US11092952B2 (en) | Plant abnormality detection method and system | |
RU2450309C1 (en) | Method of searching for faults in dynamic unit in continuous system | |
Moradi et al. | Intelligent health indicator construction for prognostics of composite structures utilizing a semi-supervised deep neural network and SHM data | |
Boškoski et al. | Bearing fault prognostics based on signal complexity and Gaussian process models | |
RU2444774C1 (en) | Method of searching for faulty module in discrete dynamic system | |
RU2473105C1 (en) | Method of detecting faults in units in continuous dynamic system | |
Wang et al. | Remaining useful life prediction of bearings based on convolution attention mechanism and temporal convolution network | |
RU2613630C1 (en) | Method for searching faulty unit in continuous dynamic system through introduction of trial deviations | |
CN116089870A (en) | Industrial equipment fault prediction method and device based on meta-learning under small sample condition | |
RU2586859C1 (en) | Method of searching for faulty unit in continuous dynamic system based on change of position of input signal | |
Du et al. | RUL prediction based on GAM–CNN for rotating machinery | |
RU2680928C1 (en) | Method for detecting faulty unit in continuous dynamic system based on sensitivity function | |
RU2506623C1 (en) | Method of searching for faulty unit in discrete dynamic system | |
RU2464616C1 (en) | Method of searching for faults in dynamic unit in continuous system | |
RU2451319C1 (en) | Method of searching for faulty module in dynamic system | |
RU2453898C1 (en) | Method of detecting faulty units in dynamic system | |
RU2579543C1 (en) | Method of troubleshooting in discrete dynamical systems based on a change in position input | |
RU2676365C1 (en) | Method of searching faulty unit in continuous dynamic system based on introduction of trial deviations | |
RU2541896C1 (en) | Method of searching for faulty unit in discrete dynamic system based on analysis of signal transmission signs | |
Heindel et al. | A data-driven approach for approximating non-linear dynamic systems using lstm networks | |
CN114970749A (en) | Method, device and readable storage medium for acquiring characteristic parameters of capacitor | |
RU2669035C1 (en) | Method of searching for faulty unit in continuous dynamic system based on change of position of input signal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FA94 | Acknowledgement of application withdrawn (non-payment of fees) |
Effective date: 20110513 |
|
FZ9A | Application not withdrawn (correction of the notice of withdrawal) |
Effective date: 20110610 |
|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20120624 |