CZ2018517A3 - Diagnostický systém strojů - Google Patents

Diagnostický systém strojů Download PDF

Info

Publication number
CZ2018517A3
CZ2018517A3 CZ2018-517A CZ2018517A CZ2018517A3 CZ 2018517 A3 CZ2018517 A3 CZ 2018517A3 CZ 2018517 A CZ2018517 A CZ 2018517A CZ 2018517 A3 CZ2018517 A3 CZ 2018517A3
Authority
CZ
Czechia
Prior art keywords
data
machine
model
sensors
diagnostic system
Prior art date
Application number
CZ2018-517A
Other languages
English (en)
Inventor
Jan OtoupalĂ­k
Josef Burian
Original Assignee
4Dot Mechatronic Systems S.R.O.
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 4Dot Mechatronic Systems S.R.O. filed Critical 4Dot Mechatronic Systems S.R.O.
Priority to CZ2018-517A priority Critical patent/CZ2018517A3/cs
Priority to SK50075-2018A priority patent/SK500752018A3/sk
Priority to US17/280,412 priority patent/US11543329B2/en
Priority to PCT/CZ2019/050039 priority patent/WO2020064033A1/en
Publication of CZ2018517A3 publication Critical patent/CZ2018517A3/cs

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • G01M99/005Testing of complete machines, e.g. washing-machines or mobile phones
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0221Preprocessing measurements, e.g. data collection rate adjustment; Standardization of measurements; Time series or signal analysis, e.g. frequency analysis or wavelets; Trustworthiness of measurements; Indexes therefor; Measurements using easily measured parameters to estimate parameters difficult to measure; Virtual sensor creation; De-noising; Sensor fusion; Unconventional preprocessing inherently present in specific fault detection methods like PCA-based methods
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0243Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
    • G05B23/0254Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model based on a quantitative model, e.g. mathematical relationships between inputs and outputs; functions: observer, Kalman filter, residual calculation, Neural Networks
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2223/00Indexing scheme associated with group G05B23/00
    • G05B2223/06Remote monitoring

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • Medical Informatics (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Physics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

Předmětem vynálezu je způsob provádění technické diagnostiky strojů pomocí diagnostického systému strojů využívajícího alespoň dvou senzorů pro umístění na strojích vybraných ze skupiny snímače vibrací, snímače přetvoření, snímače polohy a snímače vzdálenosti, kdy naměřená data jsou vyhodnocována vyhodnocovacím procesem zahrnujícím krok spárování naměřených dat a krok porovnání zpracovaných dat s modelovými stavy.

Description

Diagnostický systém strojů
Oblast techniky
Vynález se týká diagnostického systému strojů sloužícího ke snížení rizika poškození stroje, zvýšení využití strojů, snížení rizika vzniku odchylek na výrobcích, zvýšení produktivity a monitoringu procesu.
Dosavadní stav techniky
Během provozu strojů dochází k působení různých sil, které mohou způsobovat vznik značných deformací částí stroje. Během provozu může také docházet ke vzniku vysokých vibrací způsobených různými procesy během provozu. Z tohoto důvodu je vhodné využít diagnostických systémů, které by umožňovaly měření fyzikálních veličin působících v různých částech strojů.
V současném stavu techniky jsou známy systémy využívající různé druhy senzorů jako ochranných prvků. Tyto senzory měří pouze překročení určité mezní hodnoty a neumožňují tak diagnostiku strojů.
Dále jsou známy diagnostické systémy využívající různé druhy senzorů, přičemž jednotlivé druhy senzorů jsou využívány pro měření určitých jevů a z nich vyplývajících poruchových stavů. Toto řešení je však nevhodné, jelikož může docházet k chybnému rozpoznání poruchových stavů z důvodu rozpoznávání poruchových stavů pouze na základě dat naměřených jedním typem senzoru.
Bylo by tedy vhodné poskytnout diagnostický systém, který by umožňoval dostatečně přesné vyhodnocování různých poruchových stavů částí strojů, případně nástroje strojů nebo výrobku. Zároveň by řešení mělo umožňovat identifikaci různých fází provozu stroje. Toto řešení by zároveň mělo umožňovat identifikaci nových, dříve neznámých, poruchových stavů a jejich ukládání do množiny modelových stavů.
Podstata vynálezu
Výše uvedené nedostatky odstraňuje způsob provádění technické diagnostiky strojů pomocí diagnostického systému strojů zahrnujícího informační rozhraní a řídicí část, přičemž řídicí část obsahuje paměť, jehož podstata spočívá vtom, že diagnostický systém stroje dále zahrnuje alespoň dva senzory pro umístění na stroji vybrané ze skupiny snímače vibrací, snímače přetvoření, snímače polohy a snímače vzdálenosti, v paměti řídicí části je dále uložen modelový bezporuchový stav obsahující průběhy odpovídající průběhům a veličinám měřeným pomocí alespoň dvou vybraných senzorů při bezporuchovém provozu, a dále je v paměti řídicí části uložen alespoň první modelový poruchový stav a druhý modelový poruchový stav, přičemž jak první modelový poruchový stav, tak druhý modelový poruchový stav obsahují průběhy odpovídající průběhům a veličinám měřeným pomocí alespoň dvou vybraných senzorů při poruše během provozu, naměřená data ze senzorů jsou po naměření odeslána do řídicí části, následně jsou naměřená data z alespoň dvou vybraných senzorů v řídicí části zpracována vyhodnocovacím procesem obsahujícím krok spárování naměřených dat tak, aby si data z jednotlivých senzorů časově odpovídala, přičemž tato zpracovaná data jsou následně v řídicí části porovnána s modelovými stavy a na základě tohoto porovnání vydá řídicí část na informační rozhraní informaci, že se stroj nachází v bezporuchovém stavu nebo v alespoň jednom z modelových poruchových stavů. Výše popsaný způsob provádění technické diagnostiky strojů prováděný pomocí diagnostického systému strojů s alespoň dvěma senzory, jejichž naměřená data jsou dále zpracována vyhodnocovacím procesem zahrnujícím krok spárování naměřených dat, umožňuje
- 1 CZ 2018 - 517 A3 dosažení cíle dostatečně přesného vyhodnocování různých poruchových stavů, jelikož pomocí takto zpracovaných dat je možné určit poruchové stavy, které by nebylo možné určit s využitím dat naměřených pouze jedním senzorem, nebo s využitím dat naměřených více senzory, avšak bez kroku spárování takto naměřených dat, neboť při vyhodnocování se přihlíží k datům naměřeným všemi senzory.
Pokud zpracovaná data neodpovídají modelovému bezporuchovému stavu a pokud zároveň neodpovídají žádnému modelovému poruchovému stavu, vydá řídicí část na informační rozhraní informaci, že se stroj nachází v neznámém stavu. Vyhodnocení, že se jedná o neznámý stav, napomáhá dosažení cíle identifikace neznámých stavů, jelikož prostřednictvím informačního rozhraní může být tento stav uživatelem identifikován a uložen jako nový modelový stav.
Vyhodnocovací proces dále obsahuje krok úpravy dat pomocí přenosové funkce. Krok úpravy dat pomocí přenosové funkce umožňuje odstranění vlivu umístění jednotlivých senzorů v případě, že je umístění senzorů odlišné od toho umístění senzorů, pro které byly stanoveny průběhy modelových stavů.
Diagnostický systém dále zahrnuje alespoň jeden snímač teploty a vyhodnocovací proces dále obsahuje krok provedení teplotní kompenzace dat naměřených senzory vybranými ze skupiny snímače vibrací, snímače přetvoření a snímače vzdálenosti a snímače polohy, přičemž teplotní kompenzace je provedena na základě dat naměřených snímačem teploty. Teplotní kompenzací dat naměřených jinými senzory, než je snímač teploty, je možné odstranit vliv změny teploty okolí senzorů, a tím zvýšit přesnost vyhodnocování různých poruchových stavů částí stroje, nástroje nebo tvářeného výrobku.
Vyhodnocovací proces dále obsahuje krok filtrace zpracovaných dat, kdy data z alespoň jednoho senzoru jsou použita pro nastavení parametrů filtru pro filtraci dat z alespoň jednoho senzoru. Využití dat z jednoho senzoru pro nastavení parametrů filtru pro filtraci dat z alespoň jednoho senzoru umožňuje přesnější identifikaci relevantní části dat.
Diagnostický systém strojů zahrnuje informační rozhraní a řídicí část, přičemž řídicí část obsahuje paměť, a jehož podstata spočívá v tom, že diagnostický systém strojů dále zahrnuje alespoň dva senzory pro umístění na stroji vybrané ze skupiny snímače vibrací, snímače přetvoření, snímače polohy a snímače vzdálenosti, v paměti řídicí části je dále uložen modelový bezporuchový stav obsahující průběhy odpovídající průběhům a veličinám měřeným pomocí alespoň dvou vybraných senzorů při bezporuchovém provozu, a dále je v paměti řídicí části uložen alespoň první modelový poruchový stav a druhý modelový poruchový stav, přičemž jak první modelový poruchový stav, tak druhý modelový poruchový stav obsahují průběhy odpovídající průběhům a veličinám měřeným pomocí alespoň dvou vybraných senzorů při poruše během provozu. Využití diagnostického systému strojů umožňuje dosažení cíle dostatečně přesného vyhodnocování různých poruchových stavů, jelikož pomocí dvou senzorů je možné určit poruchové stavy, které by nebylo možné určit s využitím dat naměřených pouze jedním senzorem.
Strojem může být výrobní stroj.
Strojem může být dopravní prostředek.
Modelovým bezporuchovým stavem v případě, že strojem je výrobní stroj, je modelový bezporuchový stav alespoň jednoho prvku ze skupiny části stroje, nástroje nebo výrobku.
Modelovým poruchovým stavem v případě, že strojem je výrobní stroj, je kterýkoliv modelový poruchový stav ze skupiny modelového poruchového stavu části stroje, modelového poruchového stavu nástroje nebo modelového poruchového stavu výrobku.
-2 CZ 2018 - 517 A3
Modelovým bezporuchovým stavem v případě, že strojem je dopravní prostředek, je modelový bezporuchový stav části stroje.
Modelovým poruchovým stavem v případě, že strojem je dopravní prostředek, je modelový poruchový stav části stroje. Řídicí část zahrnuje monitorovací jednotku a vzdálený server, přičemž senzory jsou datově spojeny s monitorovací jednotkou a monitorovací jednotka je datově spojena se vzdáleným serverem.
V jednom případě je alespoň jedním ze senzorů snímač vibrací a alespoň jedním dalším senzorem je snímač polohy.
V jednom případě je alespoň jedním ze senzorů snímač vibrací a alespoň jedním dalším senzorem je snímač přetvoření.
V jednom případě jsou alespoň dvěma senzory snímače vibrací.
V jednom případě jsou senzory tenko vrstvě senzory.
Snímačem polohy je lineární enkodér nebo rotační enkodér.
Alespoň jedním ze senzorů je snímač teploty, přičemž snímač teploty je umístěn v těsné blízkosti alespoň jednoho senzoru jiného typu. Využitím snímače teploty je možné provádět teplotní kompenzaci dat naměřených senzory jiného typu.
Objasnění výkresů
Podstata vynálezu je dále objasněna na příkladech jeho provedení, které jsou popsány s využitím připojených výkresů, kde:
Obr. 1 znázorňuje diagnostický systém na stroji.
Obr. 2 znázorňuje vetknutý nosník se dvěma vruby.
Obr. 3 znázorňuje časový okamžik tváření 1.
Obr. 4 znázorňuje časový okamžik tváření 5.
Obr. 5 znázorňuje časový okamžik tváření 7.
Obr. 6 znázorňuje časový okamžik tváření 10.
Obr. 7 znázorňuje časový okamžik tváření 15.
Obr. 8 znázorňuje časový okamžik tváření 20.
Obr. 9 znázorňuje časový okamžik tváření 25.
Obr. 10 znázorňuje časový okamžik tváření 30.
Obr. 11 znázorňuje časový okamžik tváření 35.
Obr. 12 znázorňuje časový okamžik tváření 40.
Obr. 13 znázorňuje časový okamžik tváření 43.
-3 CZ 2018 - 517 A3
Obr. 14 znázorňuje časový okamžik tváření 45.
Obr. 15 znázorňuje časový okamžik tváření 50.
Obr. 16 znázorňuje průběh přetvoření spodní zápustky při procesu tváření.
Obr. 17 znázorňuje průběh vibrací spodní zápustky při procesu tváření.
Obr. 18 znázorňuje průběh kovací síly při procesu tváření.
Obr. 19 znázorňuje průběh polohy beranu při procesu tváření.
Obr. 20 znázorňuje průběh zrychlení horního upínače ve směru pohybu beranu při procesu tváření.
Obr. 21 znázorňuje průběh zrychlení spodní zápustky kolmo ke směru pohybu beranu pň procesu tváření.
Obr. 22 znázorňuje měřený průběh přetvoření a měřený průběh přetvoření po teplotní kompenzaci.
Obr. 23 znázorňuje změnu citlivosti piezoelektrického akcelerometru v závislosti na teplotě.
Obr. 24 znázorňuje horizontální klikový lis.
Obr. 25 znázorňuje modelový bezporuchový stav.
Obr. 26 znázorňuje modelový poruchový stav vyššího tvárného odporu přetvářeného polotovaru.
Obr. 27 znázorňuje modelový poruchový stav trhliny v upínači.
Obr. 28 znázorňuje modelový poruchový stav trhliny v rámu.
Obr. 29 znázorňuje porovnání zpracovaných dat prvního konkrétního příkladného provedení a modelového poruchového stavu vyššího tvárného odporu přetvářeného polotovaru.
Obr. 30 znázorňuje průběhy vibrací při porovnání zpracovaných dat prvního konkrétního příkladného provedení a modelového bezporuchového stavu.
Obr. 31 znázorňuje průběhy přetvoření při porovnání zpracovaných dat prvního konkrétního příkladného provedení a modelového bezporuchového stavu.
Obr. 32 zpracovaná data druhého konkrétního příkladného provedení.
Obr. 33 znázorňuje průběhy vibrací při porovnání zpracovaných dat druhého konkrétního příkladného provedení a modelového bezporuchového stavu.
Obr. 34 znázorňuje průběhy přetvoření při porovnání zpracovaných dat druhého konkrétního příkladného provedení a modelového bezporuchového stavu.
Obr. 35 znázorňuje porovnání zpracovaných dat druhého konkrétního příkladného provedení a modelového bezporuchového stavu.
-4 CZ 2018 - 517 A3
Obr. 36 znázorňuje porovnání zpracovaných dat druhého konkrétního příkladného provedení a modelového poruchového stavu trhliny v rámu.
Obr. 37 znázorňuje zpracovaná data třetího konkrétního příkladného provedení.
Obr. 38 znázorňuje průběhy vibrací při porovnání zpracovaných dat třetího konkrétního příkladného provedení a modelového bezporuchového stavu.
Obr. 39 znázorňuje průběhy přetvoření při porovnání zpracovaných dat třetího konkrétního příkladného provedení a modelového bezporuchového stavu.
Obr. 40 znázorňuje zpracovaná data prvního konkrétního příkladného provedení.
Obr. 41 znázorňuje porovnání zpracovaných dat prvního konkrétního příkladného provedení s modelovým bezporuchovým stavem.
Obr. 42 znázorňuje podmnožinu F120 modelového poruchového stavu trhliny v upínači.
Obr. 43 znázorňuje podmnožinu G130 zpracovaných dat třetího konkrétního příkladného provedení.
Obr. 44
Obr. 45 znázorňuje porovnání podmnožin F120 a G130.
znázorňuje absolutní hodnotu rozdílu funkcí podmnožin F120 a G130.
Obr. 46 příkladného provedení a modelového bezporuchového stavu.
znázorňuje průběhy vzdálenosti při porovnání zpracovaných dat čtvrtého konkrétního
Obr. 47 znázorňuje průběhy vibrací při porovnání zpracovaných dat čtvrtého konkrétního příkladného provedení a modelového poruchového stavu trhliny v upínači.
Obr. 48 znázorňuje průběhy vzdálenosti při porovnání zpracovaných dat čtvrtého konkrétního příkladného provedení a modelového poruchového stavu trhliny v upínači.
Obr. 49
Obr. 50
Obr. 51
Obr. 52
Obr. 53
Obr. 54
Obr. 55
Obr. 56
Obr. 57
Obr. 58 znázorňuje modelový poruchový stav vůle vedení.
znázorňuje modelový poruchový stav vůle zápustky.
znázorňuje modelový poruchový stav vůle upínače.
znázorňuje modelový poruchový stav vnitřního kroužku ložiska ve frekvenční oblasti.
znázorňuje modelový poruchový stav vnitřního kroužku ložiska v časové oblasti.
znázorňuje naměřená data prvním akcelerometrem v časové oblasti.
znázorňuje naměřená data třetím akcelerometrem v časové oblasti.
znázorňuje naměřená data prvním akcelerometrem ve frekvenční oblasti.
znázorňuje zpracovaná data ve frekvenční oblasti.
znázorňuje zpracovaná data v časové oblasti.
-5 CZ 2018 - 517 A3
Obr. 59 znázorňuje porovnání zpracovaných dat a modelového poruchového stavu vnitřního kroužku ložiska ve frekvenční oblasti.
Obr. 60 znázorňuje porovnání zpracovaných dat a modelového poruchového stavu vnitřního kroužku ložiska v časové oblasti.
Obr. 61 znázorňuje modelový bezporuchový stav.
Obr. 62 znázorňuje zpracovaná data pátého konkrétního příkladného provedení.
Obr. 63 znázorňuje porovnání zpracovaných dat pátého konkrétního příkladného provedení s modelovým bezporuchovým stavem.
Obr. 64 znázorňuje modelový bezporuchový stav ve frekvenční oblasti.
Obr. 65 znázorňuje zpracovaná data pátého konkrétního příkladného provedení ve frekvenční oblasti.
Obr. 66 znázorňuje porovnání zpracovaných dat pátého konkrétního příkladného provedení s modelovým bezporuchovým stavem ve frekvenční oblasti.
Obr. 67 znázorňuje válcovačku se dvěma nástroji.
Obr. 68 znázorňuje válcovačku se dvěma nástroji v řezu.
Obr. 69 znázorňuje modelový bezporuchový stav válcovačky.
Obr. 70 znázorňuje modelový bezporuchový stav válcovačky.
Obr. 71 znázorňuje modelový poruchový stav poškození zubu nástroje.
Obr. 72 znázorňuje modelový poruchový stav trhliny v polotovaru.
Obr. 73 znázorňuje zpracování nasimulovaných dat při poruchovém stavu poškození zubu nástroje.
Obr. 74 znázorňuje zpracování nasimulovaných dat při poruchovém stavu trhliny v polotovaru.
Obr. 75 znázorňuje zpracovaná data sedmého konkrétního příkladného provedení.
Obr. 76 znázorňuje zpracovaná data sedmého konkrétního příkladného provedení.
Obr. 77 znázorňuje porovnání zpracovaných dat sedmého konkrétního příkladného provedení s modelovým bezporuchovým stavem válcovačky.
Obr. 78 znázorňuje porovnání zpracovaných dat sedmého konkrétního příkladného provedení s modelovým poruchovým stavem poškození zubu nástroje.
Obr. 79 znázorňuje porovnání zpracovaných dat sedmého konkrétního příkladného provedení s modelovým poruchovým stavem trhliny v polotovaru.
Obr. 80 znázorňuje porovnání zpracovaných dat sedmého konkrétního příkladného provedení s modelovým bezporuchovým stavem válcovačky.
Obr. 81 znázorňuje zpracovaná data osmého konkrétního příkladného provedení.
-6 CZ 2018 - 517 A3
Obr. 82 znázorňuje zpracovaná data osmého konkrétního příkladného provedení.
Obr. 83 znázorňuje porovnání zpracovaných dat osmého konkrétního příkladného provedení s modelovým bezporuchovým stavem.
Obr. 84 znázorňuje porovnání zpracovaných dat osmého konkrétního příkladného provedení s modelovým bezporuchovým stavem.
Obr. 85 znázorňuje porovnání zpracovaných dat osmého konkrétního příkladného provedení s modelovým poruchovým stavem poškození zubu nástroje.
Obr. 86 znázorňuje porovnání zpracovaných dat osmého konkrétního příkladného provedení s modelovým poruchovým stavem trhliny v polotovaru.
Obr. 87 znázorňuje zpracovaná data devátého konkrétního příkladného provedení.
Obr. 88 znázorňuje zpracovaná data devátého konkrétního příkladného provedení.
Obr. 89 znázorňuje porovnání zpracovaných dat devátého konkrétního příkladného provedení s modelovým bezporuchovým stavem válcovačky.
Obr. 90 znázorňuje porovnání zpracovaných dat devátého konkrétního příkladného provedení s modelovým bezporuchovým stavem válcovačky.
Obr. 91 znázorňuje porovnání zpracovaných dat devátého konkrétního příkladného provedení s modelovým poruchovým stavem poškození zubu nástroje.
Obr. 92 znázorňuje porovnání zpracovaných dat devátého konkrétního příkladného provedení s modelovým poruchovým stavem trhliny v polotovaru.
Obr. 93 je vidět schéma diagnostikované brusky.
Obr. 94 zobrazuje modelový bezporuchový stav, modelový poruchový stav opotřebení části nástroje a modelový poruchový stav geometrie obrobku.
Obr. 95 zobrazuje spárovaná data z piezoelektrického akcelerometru a prvního rotačního enkodéru a druhého rotačního enkodéru.
Obr. 96 zobrazuje průběh zpracovaných dat $42 9 (x
Obr. 97 zobrazuje průběh zpracovaných dat $43 o (x
Obr. 98 zobrazuje trendy zpracovaných dat $42o(X), $421 (x)· $430(x) a $431 (x
Obr. 99 zobrazuje porovnání průběhů zpracovaných dat a modelových stavů.
Obr. 100 je schéma diagnostikovaného karuselového soustruhu.
Obr. 101 zobrazuje modelový bezporuchový stav, modelový poruchový stav lineárního vedení ve směru osy x a modelový poruchový stav obrobku.
Obr. 102 zobrazuje spárovaná data z piezoelektrických akcelerometrů a snímačů polohy.
Obr. 103 zobrazuje trendy zpracovaných dat $474(x). $475(^)- $476(x) a $477(x
-7 CZ 2018 - 517 A3
Obr. 104 zobrazuje porovnání průběhů zpracovaných dat a modelových stavů.
Obr. 105 zobrazuje dvoustopý dopravní prostředek.
Obr. 106 zobrazuje modelový bezporuchový stav dopravního prostředku.
Obr. 107 zobrazuje modelový poruchový stav ložiska první nápravy.
Obr. 108 zobrazuje modelový poruchový stav ložiska třetí nápravy.
Obr. 109 zobrazuje modelový poruchový stav brzdy třetí nápravy.
Obr. 110 zobrazuje modelový poruchový stav ozubeného převodu.
Obr. 111 zobrazuje zpracovaná data jedenáctého konkrétního příkladného provedení.
Obr. 112 zobrazuje porovnání zpracovaných dat jedenáctého příkladného provedení a bezporuchového stavu dopravního prostředku.
Obr. 113 zobrazuje porovnání zpracovaných dat jedenáctého příkladného provedení a poruchového stavu ložiska třetí nápravy.
Obr. 114 je schéma letadla s umístěním senzorů.
Obr. 115 zobrazuje modelový bezporuchový stav a modelový poruchový stav trhliny v patě křídla.
Obr. 116 jsou zobrazena spárovaná data.
Obr. 117 jsou zobrazena filtrovaná zpracovaná data
Příklad uskutečnění vynálezu
Příkladem vynálezu je diagnostický systém strojů. Stroji jsou buď výrobní stroje, nebo dopravní prostředky. Výrobním strojem je kterýkoliv výrobní stroj ze skupiny tvářecích strojů, obráběcích strojů nebo jakýchkoliv jiných výrobních strojů. Dopravním prostředkem je kterýkoliv dopravní prostředek ze skupiny pozemních dopravních prostředků, leteckých dopravních prostředků, vodních dopravních prostředků, železničních dopravních prostředků.
Diagnostický systém strojů zahrnuje alespoň dva senzory umístěné na stroji. Senzory jsou vybrány ze skupiny snímače vibrací, snímače přetvoření, snímače vzdálenosti a snímače polohy. Snímačem vibrací je akcelerometr. V tomto příkladném provedení snímače vibrací je použitým snímačem vibrací piezoelektrický akcelerometr s integrovanou elektronikou. Ten má oproti ostatním druhům akcelerometrů vyšší odolnost v průmyslovém prostředí a větší rozsah měření. Piezoelektrický akcelerometr může měřit pouze změny ve zrychlení, tudíž nelze měřit stacionární gravitační pole země (zrychlení „g“), a slouží tedy ke snímání absolutního zrychlení dynamických procesů. V alternativním příkladném provedení snímače vibrací je snímačem vibrací jakýkoliv druh akcelerometrů jiný než piezoelektrický akcelerometr. Snímačem přetvoření je tenzometr. Tenzometrem je kterýkoliv snímač přetvoření ze skupiny foliového, polovodičového, piezorezistivního, MEMS, optického tenzometru, či tenkostěnného tenzometru. V příkladném provedení je použit polovodičový tenzometr, který je nej výhodnějším typem pro použití v průmyslovém prostředí díky vysoké citlivosti a odolnosti proti vlivům prostředí. Snímače přetvoření jsou zde použity pro měření relativních dynamický a statických procesů.
CZ 2018 - 517 A3
Snímačem polohy je kterýkoliv snímač ze skupiny lineárního enkodéru a rotačního enkodéru. Snímačem vzdálenosti je kterýkoliv snímač ze skupiny kapacitního snímače vzdálenosti, laserového snímače polohy, konfokálního snímače vzdálenosti, optického dálkoměru, induktivního snímače, magneto-indukčního snímače nebo lankového snímače. Snímače vzdálenosti jsou použity pro měření relativních dynamických a statických procesů. Senzory jsou umístěny na stroji. Diagnostický systém strojů dále obsahuje řídicí část. V jednom z příkladných provedeních řídicí části zahrnuje řídicí část monitorovací jednotku a výpočetní zařízení. Senzory jsou datově spojeny s monitorovací jednotkou. Datové spojení senzorů s monitorovací jednotkou je drátové nebo bezdrátové. Spojení senzorů s monitorovací jednotkou je dále přímé nebo přes sběrnici nebo přes jinou výpočetní jednotku. Výpočetní jednotkou je myšlen například řídicí systém stroje. Senzory jsou dále s monitorovací jednotkou spojeny analogově nebo digitálně. V případě analogového spojení je monitorovací jednotka vybavena analogově-digitálním převodníkem, kterým je převeden signál ze senzoru tak, aby byl dále digitálně zpracovatelný. V případě digitálního spojení je senzor vybaven analogově-digitálním převodníkem a do monitorovací jednotky diagnostického systému se posílá digitální signál. Monitorovací jednotka je umístěna v blízkosti monitorovaného stroje. Monitorovací jednotka je datově spojena s výpočetním zařízením. Datové spojení monitorovací jednotky s výpočetním zařízením je drátové nebo bezdrátové. Výpočetním zařízením je vzdálený server. V alternativním provedení je výpočetním zařízením řídicí systém stroje nebo monitorovací jednotka nebo jakékoliv jiné relevantní výpočetní zařízení. Diagnostický systém strojů dále obsahuje informační rozhraní datově spojené s řídicí částí. Informačním rozhraním je jakékoliv zařízení schopné předat informaci o tom, že se stroj nachází v určitém stavu. V prvním příkladném provedení informačního rozhraní je informačním rozhraním osobní počítač. V alternativním příkladném provedení je informačním rozhraním jakékoliv elektronické zařízení s displejem nebo autonomní spolupracující systém mající vliv na provoz stroje. Datové spojení informačního rozhraní s řídicí částí je drátové nebo bezdrátové. Informační rozhraní je v jednom z příkladných provedeních součástí výpočetního zařízení.
V jednom z příkladných provedeních je monitorovací jednotka dále datově spojena s periferními zařízeními. Periferní zařízení je jakékoliv zařízení ze skupiny řídicí jednotky stroje, systému plánování výroby, měniče motoru stroje, externí databáze, dalších senzorů a firemních informačních systémů, systémů sledování kvality, různých dalších strojů majících vliv na diagnostikovaný stroj, skladových systémů, bezpečnostních systémů, měřidel nebo jakéhokoliv jiného relevantního periferního zařízení.
V jednom z příkladných provedeních dále zahrnuje diagnostický systém alespoň jeden další senzor ze skupiny snímače teploty, tenkostěnného snímače teploty, snímače síly, snímače tlaku, tenkostěnného snímače tlaku, snímače polohy, snímače rychlosti, gyroskopu, napěťového měřiče, měřiče momentu nebo jakéhokoliv jiného senzoru, který zpřesní vyhodnocování stavu stroje.
Příklad diagnostického systému strojů je viditelný na obrázku 1, kde strojem je stroj 1, dvěma senzory jsou první senzor 2 a druhý senzor 3, monitorovací jednotkou je monitorovací jednotka 4, výpočetním zařízením je výpočetní zařízení 5, řídicí částí je řídicí část 6, periferním zařízením je periferní zařízení 7 a informačním rozhraním je informační rozhraní 8.
Ve všech příkladných provedeních řídicí část dále obsahuje paměť. V paměti řídicí části jsou uloženy modelové stavy. Modelovými stavy jsou modelové poruchové stavy nebo modelový bezporuchový stav. Modelovými poruchovými stavy jsou kterékoliv modelové poruchové stavy ze skupiny modelových poruchových stavů části stroje, modelových poruchových stavů nástroje stroje nebo modelových poruchových stavů výrobku. Modelovým bezporuchovým stavem je modelový bezporuchový stav alespoň jednoho prvku ze skupiny části stroje, nástroje stroje nebo výrobku. V jednom z příkladných provedení modelového bezporuchového stavu v případě, že je diagnostikována část strojů a výrobek, je modelovým bezporuchovým stavem modelový bezporuchový stav části stroje a výrobku. V paměti řídicí části je uložen modelový bezporuchový stav a alespoň první modelový poruchový stav a alespoň druhý modelový poruchový stav. Žádný
-9 CZ 2018 - 517 A3 z modelových poruchových stavů není stejný jako jiný modelový poruchový stav. Všechny modelové stavy obsahují veličiny a průběhy, které se při provozu mění a odpovídají průběhům a veličinám měřeným při bezporuchovém provozu nebo při poruše během provozu. Fyzikálními parametry ovlivňujícími modelové stavy jsou parametry týkající se vlastností materiálu jednotlivých částí stroje, jejich vzájemného spojení, průběhů fyzikálních veličin, jako je průběh síly, geometrie stroje, průběh přetvoření, průběh napětí, dále dynamika a kinematika součástí, parametry přenosu tepla, akustické a elektromagnetické vlastnosti. V prvním příkladném provedení odvození průběhů modelových stavů jsou průběhy modelových stavů odvozeny na základě teoretického fýzikálního popisu chování stroje během provozu stroje. Takto odvozené průběhy modelových stavů umožňují použití na nových strojích bez nutnosti předchozích měření.
V alternativním provedení odvození průběhů modelových stavů jsou průběhy modelových stavů odvozeny na základě teoretického fýzikálního popisu chování stroje během provozu stroje upraveného na základě předcházejících měření reálného chování stroje během provozu stroje. Takto odvozené průběhy modelových stavů umožňují úpravy teoretického fyzikálního popisu chování stroje na základě měření reálného chování stroje během provozu stroje a jedná se tedy o zpřesněný fyzikální popis. Příkladné provedení zpřesnění fýzikálního popisu na základě měření během reálného provozu je možné na vidět obrázku 2. V tomto příkladném provedení je popisováno chování vetknutého nosníku 15 se dvěma vruby. Tento vetknutý nosník je součástí rámu stroje, kterým je v tomto příkladném provedení tvářecí stroj a fyzikální popis chování popisuje přetvoření 13 nosníku s vrubem během provozu stroje. Na základě měření vibrací prvním snímačem 11 vibrací na nosníku a druhým snímačem 12 vibrací na nosníku dojde k úpravě fýzikálního popisu chování vetknutého nosníku, a to tak, že jeden z vrubů vetknutého nosníku je virtuálně zvětšen. Virtuálním zvětšením vrubu vetknutého nosníku a fyzikálním popisem chování je popisované přetvoření během provozu stroje větší a odpovídá tak reálnému chování. Model byl zpřesněn a původní vrub 9 byl změněn na zvětšený vrub JO, který je větší než vrub 9. Na základě zpřesněného modelu je přetvoření 14 nosníku se zvětšeným vrubem od síly F větší a je blíže reálnému stroji. V dalším alternativním provedení odvození průběhů modelových stavů jsou průběhy modelových stavů odvozeny na základě matematického modelu vytvořeného strojově na základě naučených znalostí využívajících předcházejících měření reálného chování stroje během provozu. Strojové vytváření matematického modelu na základě naučených znalostí je v tomto příkladném provedení realizováno prostřednictvím neuronové sítě. Takto odvozené průběhy modelových stavů umožňují použití na strojích bez výrobní dokumentace nebo strojích, které jsou v provozu již delší dobu a různé části tak mají vůle, které nejsou zaneseny ve výrobní dokumentaci. V dalším alternativním provedení odvození průběhů modelových stavů jsou průběhy modelových stavů změřeny na základě předcházejících měření reálného chování stroje během provozu stroje. Takto odvozené průběhy modelových stavů umožňují použití bez nutnosti vytvářet jakékoliv matematické modely. Při odvozování různých modelových stavů lze využít více způsobů odvozování. Výše uvedená příkladná provedení odvozování průběhů modelových stavů je možné kombinovat a získávat tak různé modelové stavy různými způsoby. Poruchou části stroje je porucha kterékoliv části stroje, takovouto poruchou je například porucha motoru, hydraulického válce, brzdy, spojky, převodů, ložisek, vedení nebo rámu stroje. Poruchou nástroje stroje je například prasklý nástroj, opotřebení nástroje, vůle v nástroji, poškození povrchu nástroje nebo poškození geometrie nástroje. Poruchou výrobku je porucha vstupního polotovaru nebo finálního výrobku stroje, takovouto poruchou je například špatná geometrie, struktura materiálu, chemické vlastnosti materiálu nebo struktura povrchu.
Ve všech příkladných provedeních jsou data ze senzorů po naměření odeslána do řídicí části.
V řídicí části jsou naměřená data zpracována vyhodnocovacím procesem obsahujícím krok spárování naměřených dat. V kroku spárování naměřených dat jsou naměřená data alespoň ze dvou senzorů spárována tak, aby tato data měla stejnou časovou značku. Spárování dat tak, aby měla data stejnou časovou značku, znamená, že dochází ke vzájemnému posunu naměřených průběhů tak, že na časové ose naměřených průběhů začínají data vztahující se k totožnému podnětu naměřenému senzory ve stejný okamžik. V alternativním provedení kroku spárování naměřených dat jsou naměřená data alespoň ze dvou senzorů spárována tak, že jsou řazena dle určitého opakujícího se vzorce za sebe.
- 10CZ 2018 - 517 A3
Po vykonání kroku spárovaní naměřených dat se z naměřených dat stávají zpracovaná data, přičemž zpracovaná data mohou dále být zpracována jakýmikoliv dalšími kroky.
Takto zpracovaná data jsou v dalším kroku v řídicí části porovnána s modelovými stavy a na základě tohoto porovnání řídicí část vydá na informační rozhraní informaci, že se stroj nachází v bezporuchovém stavu nebo v alespoň jednom z modelových poruchových stavů. Řídicí část vydá na informační rozhraní informaci, že se stroj nachází v určitém stavu, pokud jsou zpracovaná data podobná průběhům modelových stavů. Určování podobnosti je provedeno buď prostřednictvím informačního rozhraní, nebo strojově na základě vložených znalostí, nebo na základě naučených znalostí. Strojové určování podobnosti na základě vložených znalostí je v jednom z příkladných provedení určování podobnosti realizováno prostřednictvím stavového automatu nebo vícehodnotové logiky. U stavového automatu je nutné mít jasně definovaná kritéria, je tedy vhodný pro určování podobnosti jednoznačně definovaných průběhů. Vícehodnotová logika je naopak vhodná pro práci s neurčitými výstupy. Při určování podobnosti lze využít více metod zpracování. Řídicí část vydá informaci o tom, v jakém stavu se na nachází stroj, do informačního rozhraní, které ji zobrazí.
V jednom z příkladných provedeních určování podobnosti se v případě strojového určování podobnosti na základě vložených znalostí nejprve definuje operace určení podobnosti COM, která určuje, zda si jsou průběhy podobné. Operace určení podobnosti COM je v definována vzorcem:
COM x e (A, B) [F @G*] <E, kde F je skupina modelových stavů, G* jsou data, u kterých zjišťujeme podobnost na intervalu <A,B> s podmínkami E. Operace @ zde zastupuje operaci porovnání fůnkci, jako je například rozdíl a podíl, případně rozdíl derivací, integrací a další. Operace COM zastupuje statistické vyhodnocení rozdílů například maximální hodnoty, průměrné hodnoty a dalších.
V jednom z příkladných provedení určování podobnosti se v případě strojového určování podobnosti na základě naučených znalostí vytvoří rozhodovací algoritmus na základě zpracovaných dat a označení příslušných stavů těchto dat. Určování podobnosti dle tohoto příkladného provedení zahrnuje krok přípravy dat pro učení, kdy jsou označeny příslušné stavy zpracovaných dat. Následně je proveden krok vytvoření rozhodovacího algoritmu. Rozhodovací algoritmus je v tomto příkladném provedení reprezentován neuronovou sítí a vytvoření rozhodovacího algoritmu tak zahrnuje nastavení vstupních a výstupních parametrů neuronové sítě. Následně je proveden krok učení rozhodovacího algoritmu, kdy jsou do rozhodovacího algoritmu zasílána zpracovaná data spolu s označením příslušných stavů zpracovaných dat do té doby, než je odstraněna chyba v určování stavu zpracovaných dat. Výsledkem tohoto kroku je vytvoření modelových stavů a rozhodovacího algoritmu. V dalším kroku je rozhodovací algoritmus využíván pro rozpoznání příslušných stavů stroje. V tomto příkladném provedení určování podobnosti pomocí strojového určování podobnosti na základě naučených znalostí je spolu s výše uvedenými kroky s výhodou možné provést krok vylepšení rozhodovacího algoritmu, kdy jsou do stávajícího rozhodovacího algoritmu zasílána další zpracovaná data spolu s označením příslušných stavů zpracovaných dat, přičemž tedy dochází k opětovnému učení rozhodovacího algoritmu. V tomto příkladném provedení určování podobnosti pomocí strojového určování podobnosti na základě naučených znalostí je spolu s výše uvedenými kroky s výhodou možné provést krok spárování naměřených dat, přičemž tento krok předchází kroku přípravy dat pro učení. Strojové určování podobnosti na základě naučených znalostí dle tohoto příkladného provedení je vhodné pro určování podobnosti složitých procesů s nepřesnými modely a s velkým množstvím naměřených dat k učení.
- 11 CZ 2018 - 517 A3
V dalším z příkladných provedeních určování podobnosti se v případě strojového určování podobnosti na základě naučených znalostí vytvoří rozhodovací algoritmus na základě zpracovaných dat a modelových stavů, kde si fyzikální význam zpracovaných dat a modelových stavů navzájem odpovídají. Určování podobnosti dle tohoto příkladného provedení zahrnuje krok přípravy dat pro učení, kdy ke zpracovaným datům jsou přiřazovány odpovídající modelové stavy. Následně je proveden krok vytvoření rozhodovacího algoritmu. Rozhodovací algoritmus je v tomto příkladném provedení reprezentován neuronovou sítí a vytvoření rozhodovacího algoritmu tak zahrnuje nastavení vstupních a výstupních parametrů neuronové sítě. Následně je proveden krok učení rozhodovacího algoritmu, kdy jsou do rozhodovacího algoritmu zasílána zpracovaná data spolu s modelovými stavy do té doby, než je odstraněna chyba v určování stavu zpracovaných dat. Výsledkem tohoto kroku je vytvoření rozhodovacího algoritmu. V dalším kroku je rozhodovací algoritmus využíván pro rozpoznání příslušných stavů stroje. V tomto příkladném provedení určování podobnosti pomocí strojového určování podobnosti na základě naučených znalostí je spolu výše uvedenými kroky s výhodou možné provést krok vylepšení rozhodovacího algoritmu, kdy jsou do stávajícího rozhodovacího algoritmu zasílána další zpracovaná data spolu s odpovídajícími modelovými stavy, přičemž tedy dochází k opětovnému učení rozhodovacího algoritmu. V tomto příkladném provedení určování podobnosti pomocí strojového určování podobnosti na základě naučených znalostí je spolu s výše uvedenými kroky s výhodou možné provést krok spárování naměřených dat, přičemž tento krok předchází kroku přípravy dat pro učení. Strojové určování podobnosti na základě naučených znalostí dle tohoto příkladného provedení je vhodné pro určování podobnosti bez nutnosti využití předchozích naměřených dat a je ho tak možné použít pro nové stroje.
V jednom z příkladných provedení vyhodnocovacího procesu je v kroku porovnávání zpracovaných dat a modelových stavů v případě, že zpracovaná data neodpovídají modelovému bezporuchovému stavu stroje a zároveň ani neodpovídají žádnému modelovému poruchovému stavu, řídicí částí stanoveno, že se jedná o neznámý stav. Řídicí část vydá informaci o tom, že se stroj nachází v neznámém stavu, do informačního rozhraní, které ji zobrazí a zároveň uživateli umožní přiřadit zobrazený neznámý stav bezporuchovému nebo určitému poruchovému stavu stroje. Po přiřazení zobrazeného neznámého stavu bezporuchovému nebo určitému poruchovému stavu části stroje, nástroje stroje nebo výrobku, se tato informace uloží do paměti řídicí části jako nový modelový bezporuchový nebo poruchový stav.
V jednom z příkladných provedení dále vyhodnocovací proces zahrnuje krok úpravy dat pomocí přenosové funkce. Průběhy modelových stavů jsou stanoveny pro určité rozmístění senzorů na stroji. Přenosová funkce popisuje přenos signálu z oblasti vzniku do oblasti jeho snímání.
V jednom z příkladných provedení využití přenosové funkce, upraví přenosová funkce data tak, že odstraní vliv umístění jednotlivých senzorů v případě, že je umístění senzorů odlišné od toho umístění senzorů, pro které byly stanoveny průběhy modelových stavů. V dalším příkladném provedení využití přenosové funkce umožňuje přenosová funkce nahrazení měření na více místech pomocí více senzorů téhož druhu, a to měřením na jednom místě pomocí jednoho senzoru prostřednictvím toho, že data naměřená pomocí jednoho senzoru se upraví aplikováním několika různých přenosových funkcí, přičemž každá z těchto přenosových funkcí upraví data tak, že data budou srovnatelná s daty jinak naměřenými pomocí vícera senzorů na více místech. Umístění senzorů má vliv na amplitudu a na fázový posuv jednotlivých frekvencí dat, přičemž přenosová funkce je popsána v komplexním tvaru rovnicí:
y(t) _ yoe-' <v)t+lP) _ yo_ ίφ u(ť) u0 ’ kde — je poměr amplitud a φ je fázový posuv.
u0
Přenosová funkce se získává analyticky ze znalosti mechanického chování stroje a jeho analytických rovnic. V alternativních provedeních získání přenosové funkce se přenosová funkce
- 12CZ 2018 - 517 A3 získá za pomoci parametrického fyzikálního modelu nebo modelu vytvořeného umělou inteligencí nebo experimentálně pomocí impulsní odezvy.
V jednom z příkladných provedeních kroku úpravy dat pomocí přenosové funkce jsou upravovanými daty naměřená data. V dalším z příkladných provedeních kroku úpravy dat pomocí přenosové funkce jsou upravovanými daty zpracovaná data.
V jednom z příkladných provedení obsahuje vyhodnocovací proces proces normalizace zpracovaných dat. Pro určení stavu je potřeba, aby si analyzovaná data a modelový průběh navzájem odpovídaly. To je realizováno procesem výběru intervalu zpracovaných dat a jeho normalizováním na požadovaný definiční obor odpovídající definičnímu oboru hodnot modelového průběhu. Tento proces je nezbytný pro proces strojového porovnání dat. V jednom z příkladných provedení může existovat stav který je popsán fúnkcemi Λ(χ) až /n(x) a měřená data „G“ popsána fúnkcemi g^x) až gn(x), kde n udává počet signálu popisujících stav. Při příkladném porovnání funkcí /i(x) a g^x) s různými definičními obory, se normalizuje funkce gr (x) na funkci g^(x) se stejným definičním oborem, jako má funkce (x) dle vztahu:
í X “i í J A i
9i [b. _ a- * ^di ~ ci) + ci)~ Yi g* (x) =------ --------------------------* (βι - at) + aif x e (ab bj); i = 1
Si - Yi
V jednom z příkladných provedení dále vyhodnocovací proces zahrnuje krok filtrace zpracovaných dat. V kroku filtrace zpracovaných dat jsou data alespoň jednoho senzoru použita pro nastavení parametrů filtru pro filtraci dat z alespoň jednoho senzoru.
V prvním příkladném provedení filtrace zpracovaných dat je filtrace provedena tak, že parametrem filtru je časový interval stanovený na základě dat z jednoho senzoru. Takto nastavený filtr se následně aplikuje na data z druhého senzoru a po filtraci tak z těchto dat zůstane pouze určitý časový interval. Ve druhém příkladném provedení filtrace zpracovaných dat je filtrace provedena tak, že parametrem filtruje časový průběh dat z jednoho senzoru, který pak určuje zesílení nebo útlum filtrovaných dat. Takto nastavený filtr se následně aplikuje na data z druhého senzoru a po filtraci jsou v datech potlačeny informace nenesoucí informaci o stavu. Příkladem informací nenesoucích informace o stavu, jsou rázy, posloupnost rázů, signál z jiných dílů, signál z procesu, tření součástí, určité frekvence zachycené jedním senzorem. Ve třetím příkladném provedení filtrace zpracovaných dat je filtrace provedena tak, že parametrem filtruje samotný průběh dat z jednoho senzoru. Takto nastavený filtr se následně aplikuje na data z druhého senzoru a dojde tak k odečtení dat z prvního senzoru od dat z druhého senzoru. Ve čtvrtém příkladném provedení filtrace zpracovaných dat je filtrace provedena tak, že parametrem filtru je samotný průběh dat z jednoho senzoru, kterým je snímač polohy a kinematický model diagnostikovaného stroje. Takto nastavený filtr se následně aplikuje na data z druhého senzoru, kde se tvoří aritmetický průměr ze všech měření v bodech nastavených filtrem. Tím se docílí eliminace šumu a dojde ke zvýraznění nosné informace. Měření je nutno provádět tak dlouho, dokud neobdržíme alespoň minimální počet hodnot měření v každém měřeném bodě. Ten se může měnit v závislosti na poměru šumu a nosné informace v datech.
V dalším příkladném provedení diagnostického systému strojů a vyhodnocovacího procesuje na stroji dále umístěn snímač teploty. Data naměřená snímačem teploty jsou odesílána do řídicí částí. Snímač teploty je umístěn v těsné blízkosti alespoň jednoho senzoru jiného typu. Při provozu stroje dochází k zahřívání okolí senzorů od výrobku, uvolňováním energie při provozu stroje, třením ložisek při provozu nebo přenosem tepla z dalších částí stroje. Zahřívání od zmiňovaných částí může způsobovat změny teploty v okolí senzorů, čímž dochází ke změnám citlivosti těchto senzorů. Pokud je snímač přetvoření vyroben z materiálu s odlišnou teplotní roztažností, než jakou má v daném místě materiál stroje, je tedy nutné provádět teplotní kompenzaci naměřených dat. Na obrázku 22 je tedy možné vidět rozdíl mezi stejnými naměřenými daty, přičemž průběh 26 představuje naměřená data bez provedené teplotní
- 13 CZ 2018 - 517 A3 kompenzace a průběh 25 představuje stejná naměřená data s provedenou teplotní kompenzací. V tomto příkladném provedení vyhodnocovací proces tedy zahrnuje krok teplotní kompenzace. Řídicí část poté v kroku teplotní kompenzace, na základě dat ze snímače teploty, provede teplotní kompenzaci dat naměřených senzory vybranými ze skupiny snímače vibrací, snímače přetvoření, snímače vzdálenosti a snímače polohy.
Příkladné provedení kroku teplotní kompenzace dat naměřených snímačem přetvoření se provede následujícím způsobem. Data naměřená snímačem přetvoření jsou data týkající se silového namáhání. Silové namáhání je pomocí snímače přetvoření změřeno měřením změny napětí vůči budícímu napětí, přičemž ke změně napětí dojde vlivem změny odporu v závislosti na deformaci. Změna napětí vůči budícímu napětí se převede na změnu délky vyjádřenou poměrným prodloužením ε pomocí vzorce:
kde AU/U reprezentuje změnu napětí vůči budícímu napětí a K reprezentuje převodní konstantu. Teplotní kompenzace dat se poté provádí změnou převodní konstanty K v závislosti na teplotě.
Na obrázku 23 je možné vidět průběh změny citlivosti snímače vibrací, konkrétně piezoelektrického akcelerometru. Příkladné provedení kroku teplotní kompenzace dat naměřených snímačem vibrací se provede tedy následujícím způsobem. Daty naměřenými snímačem vibrací jsou data týkající se dynamických procesů. Tyto procesy jsou měřeny změnou náboje v piezoelektriku. Tato změna náboje se poté převede na napětí, a to je následně přepočítáno na zrychlení. Teplotní kompenzace dat se poté provádí změnou náboje Bqa:
V dalším příkladném provedení vyhodnocovacího procesu obsahuje vyhodnocovací proces krok úpravy dat. V kroku úpravy dat se provede oddělení relevantní informace od zbytku dat. Relevantní informace je taková informace, která popisuje fyzikální procesy vztahující se k danému stavu. Úpravou dat je myšleno provedení alespoň jedné operace ze skupiny operací odstranění šumu, převodu dat do jiného prostoru, integrace, derivace, modulace signálu, demodulace signálu, převzorkování, tvorby trendů nebo zhodnocení statistického. Příkladem operace odstranění šumu je filtrace filtrem typu horní, dolní nebo pásmové propusti, obálková metoda a další. Převod do jiného prostoru umožňuje výhodnější zpracování a zobrazení některých poruchových stavů. Příkladem převodu do jiného prostoru je Fourierova transformace, kdy signál z časové oblasti je převeden na oblast frekvenční. Trendy jsou vytvářeny na základě výsledku určitých operací prováděných v jednotlivých časových okamžicích nad úsekem signálu. Příkladem takových operací je vyhodnocení efektivní hodnoty měřených veličin, velikosti amplitudy měřených veličin, počtu špičkových hodnot nebo délky vlny.
V jednom z příkladných provedení kroku úpravy dat jsou upravovanými daty naměřená data.
V dalším z příkladných provedení kroku úpravy dat jsou upravovanými daty zpracovaná data.
Vyhodnocovací proces vždy obsahuje minimálně kroky spárování naměřených dat a krok porovnání zpracovaných dat s modelovými stavy.
V prvním příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok teplotní kompenzace, krok úpravy dat pomocí přenosové funkce, krok spárování naměřených dat, krok filtrace zpracovaných dat, krok úpravy dat, krok normalizace zpracovaných dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
- 14CZ 2018 - 517 A3
Ve druhém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok teplotní kompenzace, krok úpravy dat pomocí přenosové funkce, krok spárování naměřených dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
Ve třetím příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok teplotní kompenzace, krok spárování naměřených dat, krok úpravy dat pomocí přenosové funkce, krok úpravy dat, krok filtrace zpracovaných dat, krok úpravy dat, krok normalizace zpracovaných dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
Ve čtvrtém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok teplotní kompenzace, krok spárování naměřených dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
V pátém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok úpravy dat pomocí přenosové funkce, krok spárování naměřených dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
V šestém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok spárování naměřených dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
V sedmém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok teplotní kompenzace, krok spárování naměřených dat, krok úpravy dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
V osmém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok úpravy dat pomocí přenosové funkce, krok spárování naměřených dat, krok úpravy dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
V devátém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok spárování naměřených dat, krok úpravy dat, krok filtrace zpracovaných dat, krok úpravy dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
V desátém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok úpravy naměřených dat krok teplotní kompenzace, krok úpravy dat pomocí přenosové funkce, krok spárování naměřených dat, krok filtrace zpracovaných dat, krok úpravy dat, krok normalizace zpracovaných dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
V jedenáctém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok teplotní kompenzace, krok úpravy dat pomocí přenosové funkce, krok spárování naměřených dat, krok úpravy dat, krok filtrace zpracovaných dat, krok normalizace zpracovaných dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
Ve dvanáctém příkladném provedení sekvence kroků vyhodnocovacího procesu je sekvence kroků vyhodnocovacího procesu následující: krok teplotní kompenzace, krok spárování naměřených dat, krok úpravy dat pomocí přenosové funkce, krok filtrace zpracovaných dat, krok úpravy dat, krok normalizace zpracovaných dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy.
- 15 CZ 2018 - 517 A3
V prvním konkrétním příkladném provedení je strojem tvářecí stroj, konkrétně klikový kovací lis zahrnující pohon, tvářecí část a rám 37. Tvářecí část zahrnuje upínač 39 sloužící k upínání nástroje zahrnujícího spodní a horní zápustku. Použitými senzory dle tohoto příkladného provedení jsou první piezoelektrický akcelerometr 27, polovodičový tenzometr 31, první snímač 43 teploty a snímač 50 síly. Všechny senzory kromě snímače 50 síly jsou umístěny na tvářecí části klikového kovacího lisu 54, konkrétně ve spodní zápustce 41. Řídicí část v tomto příkladném provedení zahrnuje monitorovací jednotku. Monitorovací jednotka je přímo datově spojena s prvním piezoelektrickým akcelerometrem 27. polovodičovým tenzometrem 31 a snímačem 43 teploty. Dále je datově spojena přes řídicí systém stroje se snímačem 50 síly umístěným na rámu 37. Snímač 50 síly umístěný na rámu 37 umožňuje získávání dat o velikosti tvářecí síly. Řídicí část diagnostického systému se skládá z monitorovací jednotky a výpočetního zařízení, kterým je vzdálený server. Dále jsou v tomto příkladném provedení v paměti řídicí části uloženy čtyři modelové stavy: modelový bezporuchový stav, který zahrnuje množinu funkcí {/ioi(%)< -/iosU)}, jehož průběh je viditelný na obrázku 25, 30 a 31, modelový poruchový stav vyššího tvárného odporu přetvářeného polotovaru, který zahrnuje množinu funkcí {/in(%)< > fns(.x)}, jehož příkladný průběh je viditelný na obrázku 26, modelový poruchový stav trhliny v rámu, který zahrnuje množinu funkcí {/i2i(x), > frzs(x)} Jehož příkladný průběh je viditelný na obrázku 28 a modelový poruchový stav trhliny na upínači, který zahrnuje množinu funkcí , fí3s(.x)}, jehož příkladný průběh je viditelný na obrázku 27. Modelový bezporuchový stav a modelový poruchový stav vyššího tvárného odporu přetvářeného polotovaru je vytvořen na základě předcházejících měření reálného chování stroje během provozu. Modelový poruchový stav trhliny v rámu a modelový poruchový stav trhliny na upínači je vytvořen na základě teoretického fyzikálního popisu chování stroje během provozu. Průběhy modelových stavů jsou vytvořeny pomocí stejného postupu zpracovaní, jak jsou zpracována aktuálně naměřená data v popisu níže. Senzory naměřená data jsou naměřena při procesu tváření.
V upínači je při procesu tváření vysoká teplotní dilatace a je tedy nejprve nutné provést krok teplotní kompenzace. Krok teplotní kompenzace polovodičového tenzometru 31 spočívá v měření teploty snímačem 43 teploty v místě měření polovodičovým tenzometrem 31 a následném využití těchto dat k přepočtu součinitele deformační citlivosti polovodičového tenzometru 31. Přepočtený součinitel deformační citlivosti polovodičového tenzometru 31 dle okamžité teploty v místě měření se následně využívá k přepočtu hodnot naměřených polovodičovým tenzometrem 31 na hodnoty reálného přetvoření. Závislost odporu polovodičového tenzometru 31 na deformaci a teplotě není lineární a je vyjádřena vztahem:
R£,t = Ro.t + ^0.25 [fi(e + («mat ~ %í)(t ~ 25)) + ύ2(ε + (amat - asi)(t - 25))2].
Odpor dále vstupuje do vztahu (platí pro 25 °C):
Λ Z?
Kg25 = Q + 2ύ2ε = Finální rovnice vyjadřující závislost součinitele deformační citlivosti na teplotě a deformaci je následující:
^,ř = ^,25(l+y^(t-25)
Konstanty Cý, C2,B jsou udávány výrobcem a jejich přesný výpočet se provádí empiricky. Pomocí výše uvedených rovnic je provedena teplotní kompenzace dat naměřených polovodičovým tenzometrem 31. Druhým krokem je spárování naměřených dat. Třetím krokem je krok úpravy dat, jejichž časové průběhy íJnyCx), áfiieCx) jsou na obrázcích 30, 31. Data z polovodičového tenzometru 31 jsou filtrována filtrem typu dolní propusti s mezní frekvencí 100 Hz, výsledná data tedy nesou informaci o přetvoření upínače v průběhu procesu tváření. Data z piezoelektrického akcelerometrů 27 jsou použita pro dva nezávislé druhy úpravy. První druh
- 16CZ 2018 - 517 A3 úpravy dat z piezoelektrického akcelerometru 27 je filtrace filtrem typu dolní propusti s mezní frekvencí 100 Hz, výsledná data tedy nesou informace o absolutním zrychlení upínače 39. Druhý druh úpravy dat z piezoelektrického akcelerometru 27 je filtrace filtrem typu horní propusti s mezní frekvencí 20 kHz, výsledná data tedy nesou informaci o rázech, kontaktu tvářecí části a výrobku a struktuře výrobku. Čtvrtým krokem je krok filtrace zpracovaných dat, kdy data naměřená piezoelektrickým akcelerometrem 27 a zároveň polovodičovým tenzometrem 31 jsou použita jako parametr filtru, kterým jsou vyfiltrována data naměřená polovodičovým tenzometrem 31. V průběhu vibrací je možné vidět v časovém okamžiku 20 impuls, který je způsoben vymezením vůle mezi rámem stroje a tvářecí částí. V průběhu přetvoření je možné vidět v časovém okamžiku 35 maximální hodnotu přetvoření. Jako parametr filtruje tedy zvolen časový interval mezi časovými okamžiky 20 a 35. Výsledkem filtrace je tedy plné utlumení průběhu přetvoření ve všech časových okamžicích mimo interval mezi časovými okamžiky 20 až 35. Pro další zpracování je použit i původní průběh přetvoření mezi okamžiky 0 až 50. Pátým krokem je krok úpravy dat. V pátém kroku jsou zpracovaná data dále statisticky vyhodnocena, přičemž data z polovodičového tenzometru 31 jsou vyhodnocena rozdílem minimální a maximální amplitudy, data z piezoelektrického akcelerometru 27 upravená filtrem typu dolní propusti s mezní frekvencí 100 Hz jsou vyhodnocena rozdílem minimální a maximální amplitudy a data z piezoelektrického akcelerometru 27 upravená filtrem typu horní propusti s mezní frekvencí 20 kHz jsou vyhodnocena pomocí výpočtu efektivní hodnoty signálu (RMS). Z takto upravených dat je následně vytvořena množina trendů {gm(x), , 5n6(+))· které jsou viditelné na obrázku 40. Šestým krokem je krok vyhodnocení podobnosti zpracovaných dat {5111(+),---5118(+)} s množinou modelových stavů {/111(+)---/118(+))- Vyhodnocení podobnosti je v tomto příkladném provedení prováděno pomocí informačního rozhraní, jímž je v tomto příkladném provedení rozhraní tvářecího stroje 54 zahrnující obrazovku a uživatelský vstup. Jak je vidět na obrázcích 30, 31 a 41 zpracovaná data {5m(+), -·,5ιΐ8(+)) neodpovídají modelovému bezporuchovému stavu určeného množinou {/101(+), --,/108(+))- jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 29 zpracovaná data {5m(x), - ,5118(+)) a modelový poruchový stav vyššího tvárného odporu přetvářeného polotovaru, určeného množinou {/111(+), ---/118(+)), si jsou podobné. Vyšší tvárný odpor přetvářeného polotovaru je způsoben špatným nastavením parametru procesu tváření, konkrétně teploty tváření.
Druhé konkrétní příkladné provedení se kromě naměřených dat a šestého kroku shoduje s prvním konkrétním příkladným provedením. Šestým krokem je krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Vyhodnocení podobnosti je v tomto příkladném provedení prováděno pomocí informačního rozhraní, jímž je v tomto příkladném provedení osobní počítač. Jak je vidět na obrázku 33, 34, 35, zpracovaná data {5121(+),----5128(+)} zobrazená na obrázku 32 neodpovídají modelovému bezporuchovému stavu určeného množinou {/101(+), - ,/108(+)/jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 36, zpracovaná data {5131(+), ---,5130(+)} a modelový poruchový stav trhliny v rámu {/131(+), ---,/130(+)), si odpovídají, jelikož průběhy si jsou podobné. Trhlina v rámu 37 je v tomto modelovém poruchovém stavu v místě rámu 37 přibližně pod tvářecí částí.
Třetí konkrétní příkladné provedení se kromě naměřených dat a pátého kroku shoduje s prvním konkrétním příkladným provedením. Pátým krokem je krok normalizace zpracovaných dat. Vstupy normalizace jsou průběhy zpracovaných dat G130 a množina funkcí F120. Množina funkcí U.20 = {/121(+),/124(+),/125(+)} zobrazená na obrázku 42 je podmnožinou množiny {/121(+), --,/120(+)} , přičemž tyto funkce reprezentují průběhy modelového poruchového stavu trhliny na upínači 39. Průběhy zpracovaných dat G130 = {5131(+), 5134(+), 5135(+)} na obrázku 43 jsou podmnožinou množiny {5131(+),, 5i3s(+)), zobrazené na obrázku 37, 38 a 39. Funkce 5131(+) reprezentuje zpracovaná data naměřená polovodičovým tenzometrem 31 mezi časovými okamžiky 20 až 35 a upravená filtrem typu dolní propusti s mezní frekvencí 100 Hz. Funkce 5134(+) reprezentuje zpracovaná data naměřená polovodičovým tenzometrem 31 mezi časovými okamžiky 0 až 50 a upravená filtrem typu dolní propusti s mezní frekvencí 100 Hz. Funkce 5i3s(+) reprezentuje zpracovaná data naměřená pomocí piezoelektrického akcelerometru 27 mezi časovými okamžiky 0 až 50 a upravená filtrem typu dolní propusti s mezní frekvencí 100
- 17CZ 2018 - 517 A3
Hz. Průběhy zpracovaných dat Guo jsou transformovány na průběhy dat Gno*, které mají stejný definiční obor a obor hodnot, jako mají funkce z množiny funkcí F120. Funkce z množiny funkcí F120 mají definiční obor D(f)t = {a^bi) a obor hodnot H(f)t = {α^, βι), kde iE{l,n). Příkladem pro transformaci ¢/134(%) na ί)(34 (x) je rovnice:
/7134 — •^134(h *U134 c134) + c134) /134
V°134 a134__________________________ ^134_ri34
()5134 - «134) + «134; x θ (a, b).
Šestým krokem je krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Vyhodnocení podobnosti je v tomto příkladném provedení prováděno strojově na základě vložených znalostí, konkrétně pomocí funkce COM. Vstupy funkce COM jsou v tomto příkladném provedení průběhy zpracovaných dat Guo* a množina funkcí F120. Funkce z množiny funkcí F120 reprezentují průběhy modelových poruchových stavů části tvářecího stroje. Funkce průběhů zpracovaných dat Guo* reprezentují zpracovaná data se stejným definičním oborem a oborem hodnot, jako mají funkce z množiny funkcí F 120, jak je zobrazeno na obrázku 44. Funkce COM je pro toto příkladné provedení definována rovnicemi:
COM rc. r „_ %e (AB)[F12O@G13° ] <E x E^bf^121^ @ £i3i(%)] < eí AND x e(a^,b2)^12^ @ #1340)] < e2 AND AND x E {a3,b3}^12^X) @ £1350)] < e3
Na základě vložených znalostí si musí být podobné všechny funkce z množiny funkcí F120 a funkce průběhů zpracovaných dat Guo*, a to tak, že funkční obor absolutní hodnoty rozdílů funkcí, zobrazený na obrázku 45, omezený kvantilem 95 % všech nejnižších hodnot na definičním oboru rozdílů funkcí je menší než 5 % definičního oboru. Na základě tohoto lze upravit rovnice funkce COM:
max0.95 z . * , .. , ,rrA max0.95 z . * , .. , ,rrA x E (a^b^^121^ ~ 9í3í^ < eiAND xe (a2,b2)^124^ ~ 9í34^ < 62 AND AND %T(u°3'9b3)l/125 W £ΐ35(χ)! < ^3 = 1 AND 1 AND 1
Porovnání funkcí je viditelné na obrázku 36, zpracovaná data {^131(^), > Vi36(x)} a modelový poruchový stav trhliny vupínači 39 {/i3i(%), -..,/130(%)} si odpovídají, jelikož průběhy si jsou podobné. Na základě výsledku vyhodnocení je do informačního rozhraní zaslána informace o tom, že poruchou je trhlina na upínači 39.
Čtvrté konkrétní příkladné provedení se shoduje se třetím konkrétním příkladným provedením s tím rozdílem, že polovodičový tenzometr 31 je nahrazen snímačem 51 vzdálenosti, v tomto konkrétním příkladném provedení indukčním, který je umístěn ve spodním upínači 39 a měří vzdálenost mezi spodním upínačem 39 a rámem 37. Modelový bezporuchový stav jev tomto příkladném provedení rozšířen o průběh /iog(%), modelový poruchový stav vyššího tvárného odporu přetvářeného polotovaru byl rozšířen o průběh /j19(x), modelový poruchový stav trhliny na upínači 39 byl rozšířen o průběh /i29(%), modelový poruchový stav trhliny v rámu 37 byl rozšířen o průběh /139(%). Doplněné průběhy popisují změnu ve vzdálenosti mezi upínačem 39 a rámem 37 během procesu tváření. Modelový bezporuchový stav je tvořen množinou dat /107(%). /109(%). Modelový poruchový stav vyššího tvárného odporu přetvářeného polotovaru je tvořen množinou dat /47(^), /119(%). Modelový poruchový stav trhliny na upínači 39 je tvořen množinou dat /i2?(%), /129(%). Modelový poruchový stav trhliny v rámu 37 je tvořen množinou
- 18 CZ 2018 - 517 A3 dat /13 7 (x)· /13 9 (x)· Zpracovanými daty využívanými v kroku vyhodnocení podobnosti zpracovaných dat s modelovými stavy jsou v tomto příkladném provedení zpracovaná data $137 (x) a $139(X)· Vyhodnocení podobnosti je v tomto příkladném provedení prováděno pomocí informačního rozhraní, jímž je v tomto příkladném provedení osobní počítač. Jak je vidět na obrázku 38, 46, zpracovaná data neodpovídají modelovému bezporuchovému stavu, jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 47 a 48, zpracovaná data a modelový poruchový stav trhliny na upínači 39 si odpovídají, jelikož průběhy si jsou podobné.
V pátém konkrétním příkladném provedení je strojem tvářecí stroj, konkrétně klikový kovací lis 54 zahrnující pohon, tvářecí část, vedení 32 a rám. Pohon zahrnuje beran 36 a převodové ústrojí. Tvářecí část zahrnuje upínač sloužící k upínání nástroje zahrnujícího spodní zápustku 41a horní zápustku 40. Použitými senzory dle tohoto příkladného provedení je pět akcelerometrů, v tomto konkrétním příkladném provedení piezoelektrických, a dva snímače teploty. První akcelerometr 27 je umístěn ve spodní zápustce 41 a snímá vibrace ve směru osy y, která je shodná s osou pohybu beranu 36. Druhý akcelerometr 28 je umístěn v horním upínači 42 nástroje a snímá vibrace ve směru osy y. Třetí akcelerometr 29 je umístěn na uložení 38 valivého ložiska, přičemž valivým ložiskem je valivé ložisko převodového ústrojí, a snímá vibrace ve směru osy y. Čtvrtý akcelerometr 30 je umístěn v horním upínači 42 a snímá vibrace kolmo ke směru osy pohybu beranu 36, tedy ve směru osy x. Pátý akcelerometr 53 je umístěn ve spodní zápustce 41 a snímá vibrace ve směru osy x. První snímač 43 teploty je umístěn ve spodní zápustce 41 a druhý snímač 44 teploty je umístěn na uložení 38 valivého ložiska. Řídicí část zahrnuje monitorovací jednotku a výpočetní zařízení, kterým jev tomto příkladném provedení řídicí systém stroje. Monitorovací jednotka je přímo datově spojena s piezoelektrickými akcelerometry a snímači teploty. Dále je v tomto příkladném provedení v paměti řídicí části uloženo pět modelových stavů: modelový bezporuchový stav, jehož průběh je viditelný na obrázku 61, 20 a 21, modelový poruchový stav vůle vedení 32, jehož průběh je viditelný na obrázku 49, modelový poruchový stav vůle upínače 39, jehož průběh je viditelný na obrázku 51, modelový poruchový stav vnitřního kroužku ložiska, jehož průběh je viditelný na obrázku 52, 53, a modelový poruchový stav vůle spodní zápustky 41, jehož průběh je viditelný na obrázku 50. Modelový bezporuchový stav je vytvořen strojově na základě naučených znalostí, přičemž data ktomu použitá jsou získána měřením při bezporuchovém provozu. Strojové vytvoření modelového bezporuchového stavu na základě naučených znalostí je realizováno prostřednictvím neuronové sítě. Prostřednictvím neuronové sítě jsou v tomto příkladném provedení tedy realizovány krok spárovaní naměřených dat, krok úpravy zpracovaných dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Modelový poruchový stav vůle vedení 32, modelový poruchový stav vůle upínače 39 a modelový poruchový stav vůle spodní zápustky 41 jsou vytvořeny stejným postupem popsaným níže. Modelový poruchový stav vůle vedení 32, modelový poruchový stav vůle upínače 39 a modelový poruchový stav vůle spodní zápustky 39 jsou odvozeny základě teoretického fýzikálního popisu chování tvářecího stroje 54 během provozu upraveného na základě předcházejících měření reálného chování tvářecího stroje 54 během provozu. Odvození probíhá vnesením vůle beranu 36 nebo vůle upínače 39 nebo vůle spodní zápustky 41 do zpřesněného fýzikálního popisu chování tvářecího stroje 54 a následnou simulací jednotlivých poruchových stavů metodou konečných prvků. Nasimulované průběhy poruchových stavů jsou následně dále zpracovány obdobně, jako jsou zpracovávána naměřená data, zpracování tedy zahrnuje krok spárování naměřených dat, krok úpravy zpracovaných dat, krok filtrace zpracovaných dat. V kroku úpravy zpracovaných dat jsou průběhy vibrací použity pro dva nezávislé druhy úprav. První druh úpravy průběhů vibrací je filtrace filtrem typu dolní propusti s mezní frekvencí 100 Hz. Druhý druh úpravy průběhů vibrací je filtrace filtrem typu horní propust s mezní frekvencí 20 kHz. Z průběhů vibrací upravených prvním druhem úpravy jsou vytvořeny trendy z rozdílů minimálních a maximálních hodnot amplitudy. Z průběhů vibrací upravených druhým druhem úpravy jsou vytvořeny trendy výpočtem efektivních hodnot signálů (RMS). Takto upravená data jsou následně filtrována krokem filtrace dat. Průběh vibrací odpovídající průběhu měřenému pátým akcelerometrem 54 je použit ke dvěma nezávislým filtracím. První filtrací jsou odfiltrována veškerá data mimo časové okamžiky 10 až 35. Druhou filtrací jsou odfiltrována veškerá data mimo časové okamžiky 7 až 43. Dále z průběhu vibrací odpovídajícímu průběhu měřeným čtvrtým akcelerometrem 30 jsou
- 19CZ 2018 - 517 A3 odfiltrována veškerá data mimo časové okamžiky 1 až 7. Takto nasimulované a zpracované průběhy jsou následně uloženy jako modelový poruchový stav vůle vedení 32. modelový poruchový stav vůle upínače 39 a modelový poruchový stav vůle spodní zápustky 41. Modelový poruchový stav vnitřního kroužku ložiska je odvozen na základě teoretického fýzikálního popisu chování tvářecího stroje 54 během provozu, přičemž tento fýzikální popis je v tomto příkladném provedení reprezentován analytickým vztahem pro chybové frekvence ložisek. Chybové frekvence zohledňují geometrii ložiska a rychlost otáčení jednotlivých kroužků. K výpočtům chybových frekvencí slouží níže uvedené rovnice:
f _ rie 'e 60’ /íp = lí I1 + ΐ cos “] _ í kde ni [RPM] označuje otáčky vnitřního kroužku, ne [RPM] označuje otáčky vnějšího kroužku, fip [Hz] označuje poruchovou frekvenci vnitřního kroužku, P [mm] označuje rozteč tělísek, Dw [mm] označuje průměr valivého tělíska, z [-] označuje počet valivých tělísek v řadě, a [°] označuje stykové úhly valivého tělíska, fl [Hz] označuje frekvence otáčení vnitřního kroužku, fe [Hz] označuje frekvence otáčení vnějšího kroužku. Za účelem zjištění aktuálních otáček vnitřního kroužku nebo vnějšího kroužku ložiska je monitorovací jednotka dále spojena s měničem motoru. Naměřená data pomocí akcelerometrů 27, 29 a snímačů 43, 44 teploty jsou vyhodnocena vyhodnocovacím procesem zahrnujícím krok teplotní kompenzace, krok spárování naměřených dat, krok úpravy dat pomocí přenosové fůnkce, krok filtrace zpracovaných dat, krok úpravy dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy je v tomto příkladném provedení realizován strojově na základě naučených znalostí prostřednictvím neuronové sítě. Rozhodovací algoritmus neuronové sítě je vytvořen na základě zpracovaných dat a modelových stavů, kde si fyzikální význam zpracovaných dat a modelových stavů navzájem odpovídají. Jak je vidět na obrázku 63 a 66, v tomto příkladném provedení jsou naměřená data {^oiCO, >9207^x)} zobrazená na obrázku 62 a 65 podobná modelovému bezporuchovému stavu {/201 (%)< ---,/207(^)), který je zobrazen na obrázku 61 a 64. Neuronová síť tak na informační rozhraní vydá informaci o tom, že se tvářecí stroj 54 nachází v bezporuchovém stavu.
Šesté konkrétní příkladné provedení se kromě průběhů naměřených dat a vyhodnocovacího procesu shoduje s pátým konkrétním příkladným provedení. Data naměřená prvním akcelerometrem 27 jsou ve tvaru m201(/cT) = {m201(lT),.,,m201(nT)}, viditelná na obrázku 54, a po Fourierově transformaci jsou tatáž data m2oi(k/F) = {m201(lF),..,m201(mF)}. viditelná na obrázku 56, data naměřená třetím akcelerometrem 29 jsou ve tvaru m203(/cT) = {m203(lT),. .,m203(nT)}, viditelná na obrázku 55 a ve tvaru ^203(^/^) = {m203(lF),..,m203(mF)}, data naměřená prvním snímačem 43 teploty jsou ve tvaru m218(kT) = {m218(lT),.., m218(nT)}. data naměřená druhým snímačem 44 teploty jsou ve tvaru m219(/cT) = {m219(lT),..,m219(nT)}, přičemž n je počet naměřených hodnot, T je perioda měření, m je počet diskrétních frekvencí a F je perioda mezi jednotlivými frekvencemi. Naměřená data pomocí akcelerometrů 29, 17 a snímačů 44, 43 teploty jsou vyhodnocena vyhodnocovacím procesem zahrnujícím postupně krok teplotní kompenzace, krok spárování naměřených dat, krok úpravy dat pomocí přenosové fůnkce, krok filtrace zpracovaných dat, krok úpravy zpracovaných dat a krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Pro úpravu dat pomocí přenosové funkce je nutné nejprve přenosovou fůnkci získat. Přenosová fůnkce TF13(z) se v tomto příkladném provedení získává experimentálně konkrétně pomocí impulsní odezvy, přičemž impulsem jev tomto příkladném provedení ráz způsobený nárazem beranu 36 skrz polotovar do rámu 37. Alternativně může být ráz vytvořen uměle. Pro určení přenosové funkce jsou užita zpracovaná data m201(/cT) = {m201(lT),.,,m201(nT)}, kde n je
-20CZ 2018 - 517 A3 počet naměřených hodnot a T je perioda měření, přičemž tato perioda měření odpovídá periodě měření pomocí prvního akcelerometru 27. Pro určení přenosové funkce jsou dále užita data m203(kT) = {m203(lT), ..,m203(nT)}, kde n je počet naměřených hodnot a T je perioda měření, přičemž tato perioda měření odpovídá periodě měření pomocí třetího akcelerometru 29. Měření dat m201(/íT) a m203(/cT) začíná před impulsem a končí s koncem jeho šíření materiálem. Měření probíhá v bezporuchovém stavu tvářecího stroje 54. Diskrétní naměřené signály jsou převáděny Z-transformací do komplexního oboru. Nejprve je použito přímé Ztransformace dle předpisu:
n
Mt(z) = Z{mi(kT)] = mi(kT)z~k fc=0 = m.i(0) + m.i(T)z 1 + m.i(2T)z 2 + ···
Přímou Z-transformací je získána spojitá funkce komplexní proměnné ML(z) z funkce mL(kT). a tím je získána funkce M201(z) a M203(z). Pro odvození přenosové funkce TF13(z) popisující přenos signálu tvářecím strojem 54 mezi prvním akcelerometrem 27, a třetím akcelerometrem 29 je využito naměřených diskrétních dat z obou akcelerometrů. Měřená data z prvního akcelerometru 27 M201(z) představují impuls. Měřená data z druhého akcelerometru 29 M203 (z), který je umístěn na uložení 38 valivého ložiska, představují tentýž impuls ovlivněný průchodem tvářecím strojem 54. Přenosová funkce TF13(z) se získá ze znalosti obou diskrétních signálů a lze vyjádřit jako
EF13(z) ^20 3 (Z) kde M201(z) je vstupní signál, M203(z) výstupní signál a TF13(z) je přenosová funkce v oblasti Z-transformace. Prvním krokem vyhodnocovacího procesu je krok teplotní kompenzace dat naměřených prvním akcelerometrem 27 m201(/cT) a třetím akcelerometrem 29 ηι203(/τΓ) pomocí dat naměřených prvním snímačem 43 teploty m218(/cT) a druhým snímačem 44 teploty m219(kT). Změna teplotní citlivosti je popsána vztahem CA(K), kde Kje teplota. Pomocí teplotní kompenzace jsou získána data a ητ2θ3(/τΓ). Teplotní kompenzaci lze popsat rovnicemi:
m^iCkT) = m20i(/cT)CA(m218(/cr)), m203 (kT) = m203(/cT)CA(m219(/cr)).
Druhým krokem vyhodnocovacího procesu je krok spárování naměřených dat m°201(kT') a m203(/cr), kde ke každé hodnotě m^^kT) je přiřazena odpovídající hodnota ητ2θ3(/τΓ). Zpracovaná data m^^kT) jsou pomocí přenosové funkce transformována na data M23(z). Transformace je provedena násobením Μ^ζζ) a TF13(z). transformace lze tedy vyjádřit rovnicí:
Μ2 τ 073(ζ) = (M^z^TF^z).
Transformovaná data M2q^3(z) v komplexním oboru lze do časového oboru převést pomocí zpětné Z-transformace dle rovnice:
mJ^3(kT) = Z 1{M^3(z)J = M^3(z)zk rdz.
J c
-21 CZ 2018 - 517 A3
Křivka C zahrnuje všechny póly výrazu Mjoi13 (z)zfe_1. Třetím krokem je krok filtrace zpracovaných dat, kde jsou od sebe odečtena data m^oiCkT) a mJo^CkT), výsledkem jsou zpracovaná data mfů-Ř/cT) na obrázku 58. Krok filtrace zpracovaných dat lze tedy vyjádřit rovnicí:
^203(^) = ^o O3(kT) - m^3(kT).
Čtvrtým krokem je krok úpravy zpracovaných dat, kdy ve zpracovaných datech mfosíkT) jsou pomocí filtrace potlačeny jiné frekvence, než je poruchová frekvence vnitřního kroužku ftp. Pátým krokem je krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. V rámci pátého kroku je provedena Fourierova transformace a data 03(^/^) na obrázku 57 jsou porovnávána ve frekvenční oblasti, přičemž porovnání je strojové na základě vložených znalostí. Chybová frekvence vnitřního kroužku ložiska se projeví až při jeho poruše. Jak je vidět na obrázku 59 a 60, v tomto příkladném provedení jsou zpracovaná data a modelový poruchový stav vnitřního kroužku ložiska si podobná a neuronová síť tak na informační rozhraní vydá informaci o tom, že se tvářecí stroj 54 nachází v modelovém poruchovém stavu vnitřního kroužku ložiska.
V sedmém konkrétním příkladném provedení je strojem tvářecí stroj, konkrétně válcovačka se dvěma nástroji viditelná na obrázku 67 a zahrnující tvářecí část. Použitými senzory dle tohoto příkladného provedení jsou dva piezoelektrické akcelerometry 18. 19 a dva rotační enkodéry 24. 47. První piezoelektrický akcelerometr 18 je umístěn na prvním valivém ložisku 20 prvního nástroje 21 a snímá tak vibrace prvního nástroje 21, druhý piezoelektrický akcelerometr 19 je umístěn na druhém valivém ložisku 22 druhého nástroje 23 a snímá tak vibrace druhého nástroje 23, první rotační enkodér 24 je umístěn na prvním nástroji 21 a snímá polohu prvního nástroje 21, druhý rotační enkodér 47 je umístěn na druhém nástroji 23 a snímá polohu druhého nástroje 23. Řídicí část zahrnuje monitorovací jednotku a výpočetní zařízení. Výpočetním zařízením je v tomto příkladném provedení vzdálený server. Monitorovací jednotka je přímo datově spojena s piezoelektrickými akcelerometry 18, 19 a přes řídicí systém stroje je datově spojena s rotačními enkodéry 24. 47. Přímé datové spojení monitorovací jednotky s piezoelektrickými akcelerometry 18. 19 je v tomto příkladném provedení realizováno drátově. Pohyb nástrojů 21. 23 se při procesu tváření v tomto příkladném provedení skládá z pohybu rotačního a z pohybu translačního. Rotační pohyb nástrojů 21, 23 je rotačním pohybem kolem vlastní osy. Translační pohyb nástrojů 21. 23 je translačním pohybem směrem vůči sobě. Kombinací rotačního a translačního pohybu nástrojů 21, 23 dochází ke vtlačování jednotlivých zubů 48 nástroje do polotovaru 16, viditelném na obrázku 68. Při vtlačování každého zubu 48 dochází ke vzniku vibrací v bodě 49 kontaktu nástroje a tvářeného výrobku. Dále jsou v tomto příkladném provedení v paměti řídicí části uloženy tři modelové stavy: modelový bezporuchový stav zahrnující {/333(^),/334(^),/335(^),/330(^)), jehož průběh je viditelný na obrázku 69, 70, modelový poruchový stav poškození zubu nástroje zahrnující {/331(^),/332(^)), jehož průběh je viditelný na obrázku 71, modelový poruchový stav trhliny v polotovaru zahrnující {/337(^),/333(^)), jehož průběh je viditelný na obrázku 72. Modelový poruchový stav poškození zubu nástroje je odvozen pomocí teoretického fyzikálního popisu chování tvářecího stroje během provozu, přičemž nasimulovaný průběh m301(kT), kde k je krok měření a T je perioda měření a který je viditelný na obrázku 73, odpovídá průběhům naměřeným prvním piezoelektrickým akcelerometrem 18 a nasimulovaný průběh m333(/cT), kde k je krok měření a T je perioda měření a který je viditelný na obrázku 73, odpovídá průběhům naměřeným prvním rotačním enkodérem 24. Takto nasimulované průběhy jsou dále zpracovány kroky, které jsou totožné s kroky vyhodnocovacího procesu. Nasimulované průběhy jsou nejprve upraveny krokem úpravy dat. V tomto kroku je nasimulovaný průběh m331(/cT) vyfiltrován filtrem typu horní propusti s mezní frekvencí 20 kHz, přičemž výsledkem této filtrace je funkce m331(/cT). Druhým krokem je krok spárovaní naměřených dat, kde je ke každé hodnotě nasimulovaného průběhu m333(/cT) přiřazena odpovídající hodnota funkce m°331(kT). Třetím krokem je krok filtrace zpracovaných dat, přičemž parametry filtru určuje zpracovaný nasimulovaný průběh m333(/cT) a filtrovanými daty je zpracovaná funkce m°331(kT). Filtrace je provedena tak, že do časové proměnné kT
-22CZ 2018 - 517 A3 zpracované funkce m331(kT') je dosazena časová proměnná kT ze zpracovaného nasimulovaného průběhu m333(/cT), čímž je získána funkce tnj31(ypTe). kde j je pořadí otáčky nástroje, p je krok polohy a Te je perioda měření polohy. Poté je pro každé p určen aritmetický průměr absolutní hodnoty vdaném bodě a je získána funkce ητ|31(ρΓ6), kde p je krok polohy a Te je perioda měření polohy, což lze vyjádřit rovnicí:
^33ΐ(Ρύ) = jZ7ilw33i(/'pTe)|.
Získaná funkce m331(pTe) má potlačený šum a je v ní zvýrazněna informace o stavu. Čtvrtým krokem je krok úpravy dat. V rámci tohoto krokuje provedena demodulace funkce m331(pTe) za pomoci filtrace pásmovou propustí obsahující frekvenci odpovídající počtu zubů z na nástroji 21. výsledkem tohoto krokuje zpracovaná funkce m331(pTe). Frekvencí odpovídající počtu zubu zt na nástroji 21 je myšlena v časové oblasti převrácená hodnota počtu zubů nástroje 21, které se dotknou výrobku 16 za jednu sekundu. Zpracovaná funkce m331(pTe) je v rámci pátého kroku úpravy dat dále statisticky vyhodnocena tak, že je určeno minimum funkce φ. poté se z daného 2π bodu minima φ funkce rozdělí na z stejných intervalů o velikost —. Na každém intervalu zL se zt určí lokální maximum. Výsledkem je množina lokálních maxim. Tato množina je dále statisticky vyhodnocena a je určeno minimum a směrodatná odchylka. Z těchto hodnot jsou poté sestaveny trendy, reprezentované funkcemi minim z maxim /331(%) a směrodatných odchylek maxim /332(%). Tyto trendy tvoří modelový poruchový stav poškození zubu nástroje. Bezporuchový modelový stav je vytvořen obdobným postupem a je tedy reprezentován množinou dvou funkcí, přičemž první funkcí je funkce nej menších maxim /333(%) a druhou funkcí je funkce směrodatných odchylek maxim /334(%). Modelový poruchový stav trhliny v polotovaru je odvozen pomocí teoretického fyzikálního popisu chování tvářecího stroje během provozu, přičemž nasimulovaný průběh m336(/cT), kde k je krok měření a T je perioda měření a který je viditelný na obrázku 74, odpovídá průběhům naměřeným prvním piezoelektrickým akcelerometrem 18, nasimulovaný průběh m332(/cT). kde k je krok měření a T je perioda měření a který je viditelný na obrázku 74, odpovídá průběhům naměřeným druhým piezoelektrickým akcelerometrem 19. nasimulovaný průběh m338(kT), kde k je krok měření a T je perioda měření a který je viditelný na obrázku 74, odpovídá průběhům naměřeným prvním rotačním enkodérem 24 a nasimulovaný průběh m334(/cT). kde k je krok měření a T je perioda měření a který je viditelný na obrázku 74, odpovídá průběhům naměřeným druhým rotačním enkodérem 47. Takto nasimulované průběhy jsou dále zpracovány kroky, které jsou totožné s kroky vyhodnocovacího procesu. Nasimulované průběhy jsou v prvním kroku upraveny krokem úpravy dat pomocí přenosové funkce. Na nasimulovaný průběh m334 (/cT) je aplikována přenosová funkce TF334337, která se získá na základě teoretického fýzikálního popisu chování tvářecího stroje během provozu. Tento krok je možné vyjádřit rovnicí:
mTF334,337 (/cr) = TF334 337
Druhým krokem je krok úpravy dat, kdy nasimulovaný průběh m332(/cT) je upraven pomocí filtru typu horní propusti s mezní frekvencí 20 kHz a výsledkem je funkce m332(kT). Třetím krokem je krok spárovaní naměřených dat, kde ke každé hodnotě funkce m|32(/cT) je přiřazena odpovídající hodnota funkce m33|34,337(/cT). Čtvrtým krokem je krok filtrace zpracovaných dat, přičemž parametry filtru určuje zpracovaná funkce m334 ’ 7 (kT) a filtrovanými daty je zpracovaná funkce m332(kT). Filtrace je provedena tak, že do časové proměnné kT zpracované funkce m°332(kT) je dosazena časová proměnná kT ze zpracované funkce m33|34337, čímž je získána funkce m332(jpTg). kde j je pořadí otáčky nástroje 47. p je krok polohy a Te je perioda měření polohy. Poté je pro každé p určen aritmetický průměr absolutní hodnoty v daném bodě a je získána funkce m332(pTe). kde p je krok polohy a Te je perioda měření polohy, což lze vyjádřit rovnicí:
-23 CZ 2018 - 517 A3 i
W-hÁPTe) =-Y |w332(/pTe)|
Získaná funkce m332(pTe) má potlačený šum a je v ní zvýrazněna informace o stavu. Pátým krokem je krok úpravy dat. V rámci tohoto krokuje provedena demodulace funkce m232(pTe) za pomoci filtrace pásmovou propustí obsahující frekvenci odpovídající počtu zubů zs výrobku, výsledkem tohoto krokuje fúnkce m332(pTe)· Frekvencí odpovídající počtu zubu zs výrobku je myšlena v časové oblasti převrácená hodnota počtu zubů nástroje 23, které se dotknou výrobku za jednu sekundu. Dále je první až pátý krok opakován pro nasimulované průběhy m336(/cT) a m338(/cT). přičemž výsledkem pátého kroku je fúnkce ιη336(ρΤβ)· Šestým krokem je krok filtrace zpracovaných dat, přičemž parametry filtru určuje zpracovaná funkce m3j2(pTe) a filtrovanými daty je zpracovaná fúnkce m^CpTg). Přičemž filtrace je provedena tak, že zpracované funkce m332(pTe) a m336(pTe) jsou od sebe odečteny a výsledkem je zpracovaná fúnkce m332(pTe). což lze vyjádřit rovnicí:
m332 (pTe) = m332(pTe) - m^36(pTe)
Zpracovaná fúnkce m332(pTg) je v rámci kroku úpravy dat dále statisticky vyhodnocena tak, že je určeno minimum fúnkce a maximum a je sledována úhlová vzdálenost mezi maximem a minimem a rozdíl maxima a minima zpracované fúnkce m332(pTg) a z těchto hodnot jsou následně vytvořeny trendy, kde úhlová vzdálenost mezi minimem a maximem je reprezentována fůnkcí /337(%) a velikost rozdílu maxima a minima je reprezentována fůnkcí /333(%). Tyto trendy tvoří modelový poruchový stav trhliny v polotovaru. Naměřená data jsou pro určení poškození zubu nástroje v tomto příkladném provedení zpracována stejnými kroky, jakými jsou zpracovány nasimulované průběhy m331(/cT) a m333(/cT) pro odvození modelového poruchového stavu poškození zubu nástroje. Zpracovaným průběhům poté odpovídají funkce $333(%) a $334(^) na obrázku 75. Naměřená data jsou pro určení trhliny v polotovaru v tomto příkladném provedení zpracována stejnými kroky, jakými jsou zpracovány nasimulované průběhy m332(/cT). m334(/cT). m336(/cT) «m338(fcf) pro odvození modelového poruchového stavu trhliny v polotovaru. Zpracovaným průběhům poté odpovídají funkce $335(%) a ¢/330(%) na obrázku 76. Krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy je v tomto příkladném provedení realizován prostřednictvím informačního rozhraní, konkrétně prostřednictvím dotykového displeje řídicího systému stroje. Jak je vidět na obrázku 77 a 80, zpracovaná data g333(x) a $334(^), ^335(^) a g336(x), modelový bezporuchový stav /333(%) a /334(^), /335(A) a /33ó(%) si odpovídají, jelikož průběhy si jsou podobné. Jak je vidět na obrázku 78, zpracovaná data g333(x). $334(%) neodpovídají modelovému poruchovému stavu poškození zubu nástroje /s3i(%), /332(A), jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 79, zpracovaná data g333(x). g33e(x) neodpovídají modelovému poruchovému stavu trhliny v polotovaru f337(x). f33s(x)· jelikož průběhy si nejsou podobné. Na základě výsledků porovnání je na informační rozhraní vydána informace, že se válcovačka nachází v bezporuchovém stavu.
Osmé konkrétní příkladné provedení se kromě průběhů naměřených dat a výsledku kroku vyhodnocení podobnosti zpracovaných dat s modelovými stavy shoduje se sedmým konkrétním příkladným provedení. Naměřená data jsou pro určení poškození zubu nástroje v tomto příkladném provedení zpracována stejnými kroky, jakými jsou zpracovány nasimulované průběhy m331(/cT) a m333(/cT) pro odvození modelového poruchového stavu poškození zubu nástroje v sedmém konkrétním příkladném provedení. Zpracovaným průběhům poté odpovídají fúnkce $331(%) a $332(^) na obrázku 81. Naměřená data jsou pro určení trhliny v polotovaru v tomto příkladném provedení zpracována stejnými kroky, jakými jsou zpracovány nasimulované průběhy m332(/cT), m334(/cT), m336(/cT) a m338(kT) pro odvození modelového poruchového stavu trhliny v polotovaru. Zpracovaným průběhům poté odpovídají funkce $337(%) a $338(%) na
-24CZ 2018 - 517 A3 obrázku 82. Krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy je v tomto příkladném provedení realizován prostřednictvím informačního rozhraní, konkrétně prostřednictvím dotykového displeje řídicího systému stroje. Jak je vidět na obrázku 83, 84, zpracovaná data $33i(x) a $332(^)- $337(^) a $338(%), modelový bezporuchový stav f333(x) a /334(A). /335(//) a/33ó(x) si neodpovídají, jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 85, zpracovaná data $33i(%), $332(^) odpovídají modelovému poruchovému stavu poškození zubu nástroje f33í(x)· f333(.x), jelikož průběhy si jsou podobné. Jak je vidět na obrázku 86, zpracovaná data $337(//), g333(%) neodpovídají modelovému poruchovému stavu trhliny v polotovaru f337(x)· f33s(.x), jelikož průběhy si nejsou podobné. Na základě výsledků porovnání je na informační rozhraní vydána informace, že se válcovačka nachází v poruchovém stavu poškození zubu nástroje, kde je porucha zubu 17.
Deváté konkrétní příkladné provedení se kromě průběhů naměřených dat a výsledku kroku vyhodnocení podobnosti zpracovaných dat s modelovými stavy shoduje se sedmým konkrétním příkladným provedení. Naměřená data jsou pro určení poškození zubu nástroje v tomto příkladném provedení zpracována stejnými kroky, jakými jsou zpracovány nasimulované průběhy m331(líT~) a m333(kT~) pro odvození modelového poruchového stavu poškození zubu nástroje v sedmém konkrétním příkladném provedení. Zpracovaným průběhům poté odpovídají funkce $33g(%) a $340(x) na obrázku 87. Naměřená data jsou pro určení trhliny v polotovaru v tomto příkladném provedení zpracována stejnými kroky, jakými jsou zpracovány nasimulované průběhy m332(kT), m334.(kT), m336(kT) a m338(kT) pro odvození modelového poruchového stavu trhliny v polotovaru. Zpracovaným průběhům poté odpovídají funkce $341(%) a $342(^) na obrázku 88. Krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy je v tomto příkladném provedení realizován prostřednictvím informačního rozhraní, konkrétně prostřednictvím dotykového displeje řídicího systému stroje. Jak je vidět na obrázku 89, 90, zpracovaná data $33g(%) a $340(//), $341(x) a $342(//) a modelový bezporuchový stav /333(A) a /334(//) , /335(//) a /330(/0 si neodpovídají, jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 91, zpracovaná data $33g(A), $340(x) neodpovídají modelovému poruchovému stavu poškození zubu nástroje f33í(x)· f333(.x), jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 92, zpracovaná data $341 (A). $342(^) odpovídají modelovému poruchovému stavu trhliny v polotovaru f337(x)· f33s(.x), jelikož průběhy si jsou podobné. Na základě výsledků porovnání je na informační rozhraní vydána informace, že se válcovačka nachází v poruchovém stavu trhliny v polotovaru.
V desátém konkrétním příkladném provedení je strojem výrobní stroj, konkrétně bruska 401, viditelná na obrázku 93. Použitými senzory jsou piezoelektrický akcelerometr 408 a dva rotační enkodéry 409. 410. Piezoelektrický akcelerometr 408 je umístěn na ložiskovém domku ložiska 406 brusného vřetene 407 brusky 401. První rotační enkodér 409 je umístěn na brusném vřetení 407 brusky 401 a snímá polohu brusného vřetene 407, respektive polohu brusného nástroje 405. Druhý rotační enkodér 410 je umístěn na unášecím vřetení 403 brusky 401 a snímá polohu unášecího vřetene 403. respektive polohu obrobku 404. Řídicí část diagnostického systému se skládá z monitorovací jednotky a výpočetního zařízení, kterým je vzdálený server. Monitorovací jednotka je přímo datově spojena s piezoelektrickým akcelerometrem 408 a oběma rotačními enkodéry 409. 410. V tomto příkladném provedení jsou v paměti řídicí části uloženy tři modelové stavy: modelový bezporuchový stav, viditelný na obrázku 94 a popsaný křivkami {fni(x), fn2(x)}, modelový poruchový stav opotřebení části brusného nástroje 405, viditelný na obrázku 5 a popsaný křivkami {f42i(x),Í422(x)} a modelový poruchový stav geometrie obrobku 404. viditelný na obrázku 94 a popsaný křivkami {f43i(x),fi32(x)}. Průběhy modelového bezporuchového stavu, modelového poruchového stavu opotřebení části brusného nástroje 405 a modelového poruchového stavu geometrie obrobku 404 jsou nasimulovány na základě teoretického fyzikálního popisu chování brusky 401 během provozu upraveného na základě předcházejících měření reálného chování brusky během provozu a jsou zpracovány stejnými kroky, kterými jsou zpracována naměřená data dle tohoto příkladného provedení. Modelový poruchový stav geometrie obrobku 404 je v tomto příkladném provedení způsoben poruchou vzdáleného převodového ústrojí hydraulického agregátu brusky 401. U modelového poruchového
-25 CZ 2018 - 517 A3 stavu opotřebení brusného nástroje 405 se porucha projevuje změnou velikosti 418 amplitudy modelového průběhu křivky £417 v konkrétním úhlovém natočení 419 hřídele. U modelového poruchového stavu geometrie obrobku 404 se porucha projevuje na frekvenci 426 poruchy geometrie obrobku 404, která je vizualizovaná na modelovém průběhu křivky f423. přičemž frekvence 426 poruchy geometrie obrobku 404 je určena největší amplitudou ve frekvenčním spektru křivky f423 určené křivkou f425 a velikostí 424 amplitudy průběhu modelové křivky /423. Naměřená data z piezoelektrického akcelerometru 408 jsou označena m1(kT) = {mj/lT),.. naměřená data z prvního rotačního enkodéru 409 jsou označena m2(kT) = {m2(IT),..,m2(nT)} a naměřená data z druhého rotačního enkodéru 410 jsou označena m3(kT) = {m3(lT),. ,,m3(nT)}, přičemž n je počet naměřených hodnot a T je perioda měření. Tato naměřená data jsou vyhodnocena vyhodnocovacím procesem, kde prvním krokem je krok spárování naměřených dat. V kroku spárování naměřených dat je ke každé hodnotě m1(kT) přiřazena odpovídající hodnota m2(kT) a m3(/cT). Takto spárovaná data již nesou informaci o hodnotě vibrací v každé poloze hřídele 407, 403. Spárovaná data jsou znázorněna na obrázku 95. Druhým krokem je krok filtrace zpracovaných dat. V tomto kroku se data z piezoelektrického akcelerometru 408 sčítají v závislosti na poloze brusného vřetene 407 a unášecího vřetene 403 a kinematice diagnostikované soustavy. Tím vzniknou data, která odpovídají přesně jedné periodě diagnostikované části stroje. V tomto kroku je tedy provedena filtrace dat m1 (kT) z piezoelektrického akcelerometru 408 filtrem daným daty m2 (kT) z prvního rotačního enkodéru 409 umístěného na brusném vřetení 407. Data m1(kT) se tedy sčítají přesně v každém měřeném bodě T409 = tn + 2π brusného vřetene 407, kde tn je daným bodem měření. Zároveň je provedena filtrace dat m^kT) z piezoelektrického akcelerometru 408 filtrem daným daty m3(/cT) druhého rotačního enkodéru 410 umístěného na unášecím vřetení 403. Data m^kT) se tedy zároveň sčítají přesně v každém měřeném bodě T410 = ίη + 2π unášecího vřetene 403, kde tn je daným bodem měření. Pro každý měřený bod polohy p je dále určen aritmetický průměr zpracovaných dat dle vzorců:
i mí2(pT2) = g429(x) = -Y m1 (jpT2) mi3(pT3) = g430(x) = yZllii(/pT3)|, kde j je pořadí otáčky dané hřídele. U takto upravených dat dojde při sečtení dat z dostatečného množství otáček k potlačení šumu, který je u měřených dat na obráběcích strojích vysoký, a tím dojde ke zvýraznění nosné informace. Třetím krokem je krok úpravy dat, kdy jsou zpracovaná data filtrována filtrem typu dolní propusti s mezní frekvencí 100 Hz. Průběhy zpracovaných dat $429 a g430 jsou zobrazeny na obrázcích 96 a 97. Dále je v tomto kroku provedena úprava dat g429 , ze kterých se tvoří množiny trendů {¢/421,//422) a úprava dat g430 , ze kterých se tvoří množiny trendů {//431,//432)- Zpracovaná data g4so jsou převedena rychlou Fourierovou transformací do frekvenční oblasti a je sledováno, zda je porucha na frekvenci poruchy modelového poruchového stavu brusného nástroje 405. Hodnota frekvence poruchy tvoří křivku g432. Velikost amplitudy poruchy je analyzována metodou Min-Max, která spočívá v odečtení nejnižší hodnoty zpracovaných dat od té nejvyšší. Trend těchto amplitud tvoří křivku g434. U zpracovaných dat g429 se sleduje, zdaje porucha ve stejném úhlovém natočení brusného vřetene 407. Trend hodnot těchto úhlových natočení tvoří křivku g422. Velikost amplitudy poruchy se analyzuje metodou Min-Max. Trend těchto amplitud tvoří křivku g424, jež je vidět na obrázku 98. Čtvrtým krokem je krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Vyhodnocení podobnosti je v tomto příkladném provedení prováděno pomocí informačního rozhraní, jímž je v tomto příkladném provedení osobní počítač. Jak je vidět na obrázku 99, zpracované skupiny dat {//421,//431) neodpovídají modelovému bezporuchovému stavu popsanému křivkami [f4íí, f4í2] , jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 99 zpracované množiny dat {//421,//431) neodpovídají modelovému poruchovému stavu opotřebení brusného nástroje 405 popsanému křivkami {/411,/412). jelikož průběhy si nejsou podobné. Jak
-26CZ 2018 - 517 A3 je vidět na obrázku 99 zpracované množiny dat {5431,5432} odpovídají modelovému poruchovému stavu geometrie obrobku 404 popsanému křivkami {/411,/412), jelikož průběhy si jsou podobné.
V jedenáctém konkrétním příkladném provedení je strojem výrobní stroj, konkrétně karuselový soustruh 445 zobrazený na obrázku 100. Použitými senzory jsou dva piezoelektrické akcelerometry 446, 447 s integrovanými snímači teploty. První piezoelektrický akcelerometr 446 je umístěn na vřetenu 448 lože. Druhý piezoelektrický akcelerometr 447 je umístěn na lineárním vedení 449 příčníkové podpory. Řídicí část diagnostického systému obráběcích strojů se skládá z monitorovací jednotky a výpočetního zařízení, kterým je vzdálený server. Monitorovací jednotka je přímo datově spojena s piezoelektrickým akcelerometry 446, 447 a dále je datově spojena s řídicím systémem karuselového soustruhu 445. Z řídicího systému karuselového soustruhu 445 jsou získávány informace o poloze příčníkové podpory na lineárním vedení 449 ve směru osy x, o úhlové poloze obrobku 452 vůči nástroji 453 a o poloze nástroje 453 vůči obrobku 452 v ose z. V tomto příkladném provedení jsou v paměti řídicí části uloženy tři modelové stavy: modelový bezporuchový stav, viditelný na obrázku 101 a popsaný množinou křivek Fx(x) = {/407(+)-/468(+)), modelový poruchový stav lineárního vedení 449 v ose x, viditelný na obrázku 101 a popsaný množinou křivek F2(x) = {/459(+),/400(+)) a modelový poruchový stav obrobku 452. viditelný na obrázku 101 a popsaný křivkami F3(x) = {/461(+),/402(+))- Modelový poruchový stav lineárního vedení 449 v ose x je odvozen na základě teoretického fýzikálního popisu chování karuselového soustruhu 445 během provozu, kde se porucha projevuje změnou velikosti amplitudy 464 modelového průběhu křivky fw, v konkrétní poruchové poloze 465 lineárního vedení 449. Modelový poruchový stav obrobku 452 je odvozen na základě teoretického fyzikálního chování karuselového soustruhu 445 během provozu, konkrétně metodou konečných prvků. Naměřená data z prvního piezoelektrického akcelerometru 446 jsou označena m4(kT) = ..,m1(nT)j, naměřená data z druhého piezoelektrického akcelerometru 447 jsou označena m2(kT) = (m2(lT),..,m2(nT)j, naměřená data z prvního snímače teploty integrovaného na prvním piezoelektrickém akcelerometru 446 jsou označena m3(kT) = (m3(lT),..,m3(nT)j, naměřená data ze druhého snímače teploty integrovaného na druhém piezoelektrickém akcelerometru 447 jsou označena m4(kT) = (m4(lT),..,m4(nT)j, data o úhlovém natočení obrobku 452 vůči nástroji 453 jsou označena m5(kT) = {m5(lT), ..,m5(nT)}, data o poloze nástroje 453 vůči obrobku 452 vose z jsou označena m6(kT) = {m6(lT),..,m6(nT)} a data o poloze nástroje 452 vůči lineárnímu vedení 449 ve směru osy x jsou označena m7(/cT) = {m7(lT),..,m7(nT)}, přičemž n je počet naměřených hodnot a T je perioda měření. Prvním krokem je krok teplotní kompenzace, přičemž teplotně kompenzovanými daty jsou naměřená data m1(kT) a m2(/cT) a jsou kompenzována pomocí naměřených dat m3(kT) a m4(kT). Změna teplotní citlivosti je popsána vztahem Cý(K), kde K je teplota. Pomocí teplotní kompenzaci jsou získána data m°4(kT) a m°2 (kT).
m4°(kT) = m^kT) * Cý(m3(/cT)) m2°(kT) = m2(/cT)Cý(m4(/cT)).
Druhým krokem je krok úpravy dat pomocí přenosové funkce. Vzhledem ke vzdálenému umístění pouze prvního piezoelektrického akcelerometru 446 od zdroje vibrací vzniklých během procesu obrábění, je přenosová funkce aplikována pouze na naměřená data prvního piezoelektrického akcelerometru 446. Přenosová funkce je získána experimentálně pomocí impulsní odezvy, kde impulsem je impuls od procesu obrábění. V alternativním příkladném provedení je možné impuls vytvořit i uměle a sledovat jeho šíření materiálem. Pro získání přenosové funkce jsou využita naměřená data m4 (kT) = [m4 (1T),.., m4 (nT)}, kde n je počet naměřených hodnot a T je perioda měření. Snímání dat začíná před impulsem a končí s koncem jeho šíření materiálem. Diskrétní naměřená data jsou převáděna Z-transformací do komplexního oboru. Nejprve je použito přímé Z-transformace dle jejího předpisu:
-27CZ 2018 - 517 A3 n
Mi(zj = Z{mi(kT)] = mi(kT)z~k = mi(0) + m^Tjz^1 + mi(2T)z~2 fc=0
Přímou Z-transformací je získána spojitá funkce komplexní proměnné M^zj z funkce m^kT). Dosazením m*(kTj a mp(kT) do vzorce pak získáme funkce Mf(z) a Mp(z). Data Mp(zj představují impuls. Měřená data ze snímače Mý (z), který je umístěn na uložení ložiska vřetene 448 lože představují tentýž impuls, který je ovlivněn charakterem systému, tedy přenosovou funkcí. Přenosová funkce TFip (z) se získá znalostí obou signálů a lze vyjádřit jako:
Mp(z)
TF1D(z) = — lp 7 M*(z)
Aplikace přenosové funkce je provedena násobením Μ^(ζ) a TFlp(z). Všechna data m^(kT) je tedy nutné nejprve převést Z-transformací na M^(z) a aplikaci pak lze vyjádřit rovnicí:
<Flp(z) = (M^TFlp(z)
Z komplexního oboru se přejde do časového oboru pomocí zpětné Z-transformace dle předpisu:
m^lp(kT) = Z 1{M^l Flp(z)] = Mmlp(z)zk rdz J c
Křivka C zahrnuje všechny póly výrazu M^Flp(z)zfe 1. Třetím krokem je krok spárování naměřených dat. V kroku spárování naměřených dat je ke každé hodnotě mFFlp (kT) je přiřazena odpovídající hodnota m5(kT) a m6(kT) a ke každé hodnotě m2(kT) je přiřazena odpovídající hodnota m7(kT). Takto spárovaná data již nesou informaci o hodnotě vibrací v každém měřeném bodě. Spárovaná data jsou znázorněna na obrázku 102. Čtvrtým krokem je krok filtrace zpracovaných dat, kde se zpracovaná data mFFlp(kT) průměrují ve stejných bodech prostoru určených zpracovanými daty polohy m5(kTj a m6(kT). Zpracovaná data m2(k.Tj jsou průměrována ve stejných bodech určených daty polohy m7(kT). U takto filtrovaných dat dojde při sečtení dat z dostatečného množství otáček k potlačení šumu, který je u měřených dat na obráběcích strojích vysoký. Pátým krokem je krok úpravy dat, kdy zpracovaná data jsou dále filtrována filtrem typu horní propusti s mezní frekvencí 20 kHz. Z takto zpracovaných dat se tvoří množina trendů hodnot G postupem shodujícím se s postupem tvorby modelových stavů. Množina G obsahuje množiny dat ($474(x), $476(^)1(^474^)-0475WM0476W-0477 W) zobrazené na obrázku 103. Šestým krokem je krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Vyhodnocení podobnosti jev tomto příkladném provedení prováděno strojově na základě vložených znalostí, přičemž vstupy je množina modelových stavů F složená z modelových stavů {Ž459W-/460WK/461W-/462 W1ÍČ467 W-AósW) a množina upravených dat G složená z upravených dat {$474 W-0476W1Í0474W-0475 W1Í0476W-0477 W}· Krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy v tomto příkladném provedení zahrnuje krok normalizace zpracovaných dat, kdy množina G je transformována na množinu G*, která má stejný definiční obor a funkční obor, jako mají funkce F. Funkce F mají definiční obor Dfj = {ct-i, bi) a funkční obor Hfi = {αι,βι}, kde i = 1,2,3. Příkladem transformace je transformace g474(t) na $474(1:) dle rovnice:
^474(h *Ur74 c474) + c474) /474 $474* = ----- o-------------------(/?474 - «474) + «474; t e (a,b).
°474 /474
Porovnání funkcí je realizováno prostřednictvím operace COM, která určuje podobnost funkcí:
-28 CZ 2018 - 517 A3
COM
X e (A, B) [F @ G*] <E, kde E je množina podmínek {e1; e3, e3}. Porovnání funkcí je viditelné na obrázku 104, přičemž v tomto příkladném provedení je funkcí COM vyhodnoceno, že se karuselový soustruh 445 nachází ve stavu poruchového stavu geometrie obrobku 452 a vydá tuto informaci na informační rozhraní.
Ve dvanáctém konkrétním příkladném provedení je strojem dvoustopý pozemní dopravní prostředek 508. konkrétně dvoustopý dopravní prostředek se samonosnou karoserií a třemi nápravami 504. 505 a 506. Každá z náprav 504. 505 a 506 obsahuje ložisko a mechanickou brzdu, viditelné na obrázku 105. Dvoustopý dopravní prostředek 508 je poháněn pohonnou jednotkou 507 s převodovým ústrojím 509. Použitými senzory dle tohoto příkladného provedení jsou tři piezoelektrické akcelerometry 501. 502. 503. První piezoelektrický akcelerometr 501 je umístěn na první nápravě 504, přičemž první náprava 504 obsahuje první ložisko, druhý piezoelektrický akcelerometr 502 je umístěn na druhé nápravě 505, přičemž druhá náprava 505 obsahuje druhé ložisko, a třetí piezoelektrický akcelerometr 503 je umístěn na třetí nápravě 506. přičemž třetí náprava 506 obsahuje třetí ložisko. Řídicí část zahrnuje monitorovací jednotku 510 a výpočetní zařízení, kterým jev tomto příkladném provedení vzdálený server. Monitorovací jednotka 510 je přímo datově spojena s piezoelektrickými akcelerometry 501. 502. 503 a dále je datově spojena s řídicí jednotkou dvoustopého dopravního prostředku 508. Dále je v tomto příkladném provedení v paměti řídicí části uloženo pět modelových stavů: modelový bezporuchový stav reprezentovaný funkcemi f501 (x). /502(+) a /503(+), jehož průběh je viditelný na obrázku 106, modelový poruchový stav prvního ložiska reprezentovaný fůnkcemi /511(+), /512(+), /513(+) a kde zpoždění signálu Δχ je menší než zpoždění signálu Δ2, jehož průběh je viditelný na obrázku 107, modelový poruchový stav druhého ložiska reprezentovaný fůnkcemi /521(+), /522(+), /523(+) a kde zpoždění signálu Δ3 je menší než zpoždění signálu Δ4, jehož průběh je viditelný na obrázku 108, modelový poruchový stav brzdy třetí nápravy 506 reprezentovaný funkcemi /531(x), /532(+), /533(+) a kde zpoždění signálu Δ5 je menší než zpoždění signálu Δ6, jehož průběh je viditelný na obrázku 109, a modelový poruchový stav převodového ústrojí 509 reprezentovaný funkcemi /541(x), /542(+), /543(+) a kde zpoždění signálu Δ7 je menší než zpoždění signálu Δ8,jehož průběh je viditelný na obrázku 110. Všechny modelové stavy byly vytvořeny na základě teoretického fyzikálního popisu chování stroje 508 během provozu stroje upraveného na základě předcházejících měření reálného chování dvoustopého dopravního prostředku 508 během provozu dvoustopého dopravního prostředku 508. Naměřená data jsou zpracována vyhodnocovacím procesem obsahujícím krok spárování naměřených dat a krok úpravy dat. Zpracovaná data jsou reprezentována fůnkcemi #101(+)#102(+), #103(+) viditelnými na obrázku 111, přičemž zpracování odpovídá zpracování průběhů modelových stavů. Dalším krokem vyhodnocovacího procesu je krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy jev tomto příkladném provedení realizován pomocí stroje s vloženými znalostmi. Jak je vidět na obrázku 112, zpracovaná data neodpovídají modelovém bezporuchovému stavu, jelikož průběhy si nejsou podobné. Jak je vidět na obrázku 113, zpracovaná data odpovídají modelovému poruchovému stavu třetího ložiska, jelikož průběhy si jsou podobné. Na základě výsledků porovnání je na informační rozhraní vydána informace, že se dvoustopý dopravní prostředek nachází v poruchovém stavu druhého ložiska.
Ve třináctém příkladném provedení je strojem letecký dopravní prostředek, konkrétně letadlo
601. viditelné na obrázku 114. Použitými senzory dle tohoto příkladného provedení jsou tenkovrstvý snímač 603 přetvoření, piezoelektrický akcelerometr 604 a snímač 605 teploty. Tenkovrstvý snímač 603 přetvoření je umístěn na nosném profilu 606 křídla 602 v patě křídla
602, piezoelektrický akcelerometr 604 je umístěn na nosném profilu 606 křídla 602 ve středu křídla 602 a snímač 605 teploty je umístěn v těsné blízkosti tenkovrstvého snímače 603 přetvoření. Řídicí část zahrnuje monitorovací jednotku a výpočetní zařízení, kterým je v tomto
-29CZ 2018 - 517 A3 příkladném provedení řídicí systém letadla 601. Monitorovací jednotka je přímo datově spojena se senzory. Dále jsou v tomto příkladném provedení v paměti řídicí části uloženy tři modelové stavy: modelový bezporuchový stav, modelový poruchový stav trhliny v patě křídla 602 a modelový poruchový stav trhliny ve středu křídla 602. Modelový bezporuchový stav, modelový poruchový stav trhliny v patě křídla 602 a modelový poruchový stav trhliny středu křídla 602 jsou odvozeny na základě simulace teoretického fyzikálního popisu chování stroje 601 během provozu, konkrétně metodou konečných prvků, kde jako proměnné vstupují předpětí křídla 602 a velikost změn síly vyvozovaných na křídlo 602. velikost trhliny v patě křídla 602 a velikost trhliny ve středu křídla 602. Výstupem této simulace je tuhost křídla 602 reprezentovaná odezvou systému měřenou piezoelektrickým akcelerometrem 604 a popsaná fůnkci fn = f(F,AF, t1; t2). Funkce fn je dále rozdělena na množiny jednotlivých křivek modelových poruchový stavů {/óoiCEnAF^tn^tzi), ...,/6nn(Fn,AFn,tln,t2n)}. Na obrázku 115 je znázorněn modelový bezporuchový stav pro hodnotu síly předpětí F02 a změnu síly působící na křídlo 602 AF2. Na obrázku 115 je dále znázorněn modelový poruchový stav trhliny v patě křídla 602 pro hodnotu síly předpětí F01, změnu síly působící na křídlo AFr a velikost trhliny t. Naměřená data z tenkovrstvého snímače 603 přetvoření jsou označena m^kT) = , m^nT}}, naměřená data z piezoelektrického akcelerometrů 604 jsou označena m2(/cT) = {m2(lT),..,m2(nT)} a naměřená data ze snímače 605 teploty jsou označena m3(kT) = {m3 (1T), ..,m3 (nT)}, kde n je počet naměřených hodnot a T je perioda měření. Naměřená data jsou vyhodnocena vyhodnocovacím procesem, kde prvním krokem je krok teplotní kompenzace, druhým krokem je krok spárování naměřených dat a třetím krokem je krok úpravy dat. Krokem spárování naměřených dat jsou zpracovaná data viditelná na obrázku 116. Ve třetím kroku jsou ve zpracovaných datech z tenkovrstvého snímače 603 přetvoření vyhledány úseky dat g6nn délky dx. Úseky dat g^nn zachycují změnu přetvoření o velikosti alespoň AFmin původní hodnoty a hodnoty předcházející této změně, která má setrvalý stav, který je charakteristický jeho maximálním směrodatnou odchylkou od průměru AFmm/10 alespoň po dobu dx/2 a který je charakteristický jeho maximální směrodatnou odchylkou od průměru AFmin/10 po dobu dx/2, jak je vidět na obrázku 117. K úsekům dat g^nn jsou následně přiřazena zpracovaná data z piezoelektrického akcelerometrů 604, jak je vidět na obrázku 117. Čtvrtým krokem je krok vyhodnocení podobnosti zpracovaných dat s modelovými stavy. Vyhodnocení podobnosti je realizováno strojově na základě s vložených znalostí. Během vyhodnocení podobnosti je porovnána množina zpracovaných dat {g601(F1,AF1), , genn(Fn, AFn~)} s množinou dat popisující modelový bezporuchový stav, modelový poruchový stav trhliny v patě křídla a modelový poruchový stav trhliny středu křídla {/eoiÚT, AFn t11; t21),, f6nn(Fn, AFn, tln, t2n)}. Na základě výsledku vyhodnocení podobnosti vydá řídicí část na informační rozhraní informaci o tom, zda se stroj nachází v bezporuchovém stavu, v poruchovém stavu trhliny v patě křídla 602, v poruchovém stavu trhliny ve středu křídla 602 nebo v neznámém stavu.

Claims (19)

  1. PATENTOVÉ NÁROKY
    1. Způsob provádění technické diagnostiky strojů pomocí diagnostického systému strojů zahrnujícího informační rozhraní a řídicí část, přičemž řídicí část obsahuje paměť, vyznačující se tím, že diagnostický systém strojů dále zahrnuje alespoň dva senzory pro umístění na stroji vybrané ze skupiny snímače vibrací, snímače přetvoření, snímače polohy a snímače vzdálenosti, v paměti řídicí části je dále uložen modelový bezporuchový stav obsahující průběhy odpovídající průběhům a veličinám měřeným pomocí alespoň dvou vybraných senzorů při bezporuchovém provozu, a dále jev paměti řídicí části uložen alespoň první modelový poruchový stav a druhý modelový poruchový stav, přičemž jak první modelový poruchový stav, tak druhý modelový poruchový stav obsahují průběhy odpovídající průběhům a veličinám měřeným pomocí alespoň dvou vybraných senzorů při poruše během provozu, naměřená data ze senzorů jsou po naměření odeslána do řídicí části, následně jsou naměřená data z alespoň dvou vybraných senzorů v řídicí části zpracována vyhodnocovacím procesem obsahujícím krok spárování naměřených dat tak, aby si data z jednotlivých senzorů časově odpovídala, přičemž tato zpracovaná data jsou následně
    -30CZ 2018 - 517 A3 v řídicí části porovnána s modelovými stavy a na základě tohoto porovnání vydá řídicí část na informační rozhraní informaci, že se stroj nachází v bezporuchovém stavu nebo v alespoň jednom z modelových poruchových stavů.
  2. 2. Způsob provádění technické diagnostiky strojů pomocí diagnostického systému strojů podle nároku 1 vyznačující se tím, že pokud zpracovaná data neodpovídají modelovému bezporuchovému stavu a pokud zároveň neodpovídají žádnému modelovému poruchovému stavu, vydá řídicí část na informační rozhraní informaci, že se stroj nachází v neznámém stavu
  3. 3. Způsob provádění technické diagnostiky strojů pomocí diagnostického systému strojů podle kteréhokoliv z předchozích nároků vyznačující se tím, že vyhodnocovací proces dále obsahuje krok úpravy dat pomocí přenosové funkce.
  4. 4. Způsob provádění technické diagnostiky strojů pomocí diagnostického systému strojů podle kteréhokoliv z předchozích nároků vyznačující se tím, že diagnostický systém dále zahrnuje alespoň jeden snímač teploty a vyhodnocovací proces dále obsahuje krok provedení teplotní kompenzace dat naměřených senzory vybranými ze skupiny snímače vibrací, snímače přetvoření, snímače vzdálenosti a snímače polohy, přičemž teplotní kompenzace je provedena na základě dat naměřených snímačem teploty.
  5. 5. Způsob provádění technické diagnostiky strojů pomocí diagnostického systému strojů podle kteréhokoliv z předchozích nároků vyznačující se tím, že vyhodnocovací proces dále obsahuje krok filtrace zpracovaných dat, kdy data z alespoň jednoho senzoru jsou použita pro nastavení parametrů filtru pro filtraci dat z alespoň jednoho senzoru.
  6. 6. Diagnostický systém strojů zahrnuje informační rozhraní a řídicí část, přičemž řídicí část obsahuje paměť, vyznačující se tím, že diagnostický systém strojů dále zahrnuje alespoň dva senzory pro umístění na stroji vybrané ze skupiny snímače vibrací, snímače přetvoření, snímače polohy a snímače vzdálenosti, v paměti řídicí části je dále uložen modelový bezporuchový stav obsahující průběhy odpovídající průběhům a veličinám měřeným pomocí alespoň dvou vybraných senzorů při bezporuchovém provozu, a dále je v paměti řídicí části uložen alespoň první modelový poruchový stav a druhý modelový poruchový stav, přičemž jak první modelový poruchový stav, tak druhý modelový poruchový stav obsahují průběhy odpovídající průběhům a veličinám měřeným pomocí alespoň dvou vybraných senzorů při poruše během provozu.
  7. 7. Diagnostický systém strojů podle nároku 6 vyznačující se tím, že strojem je výrobní stroj.
  8. 8. Diagnostický systém strojů podle nároku 6 vyznačující se tím, že strojem je dopravní prostředek.
  9. 9. Diagnostický systém strojů podle nároku 6 nebo 7 vyznačující se tím, že modelovým bezporuchovým stavem je modelový bezporuchový stav alespoň jednoho prvku ze skupiny části stroje, nástroje nebo výrobku.
  10. 10. Diagnostický systém strojů podle kteréhokoliv z nároků 6, 7 a 9 vyznačující se tím, že modelovým poruchovým stavem je kterýkoliv modelový poruchový stav ze skupiny modelového poruchového stavu části stroje, modelového poruchového stavu nástroje nebo modelového poruchového stavu výrobku.
  11. 11. Diagnostický systém strojů podle nároku 6 nebo 8 vyznačující se tím, že modelovým bezporuchovým stavem je modelový bezporuchový stav části stroje.
  12. 12. Diagnostický systém strojů podle kteréhokoliv z nároků 6, 8 a 11 vyznačující se tím, že modelovým poruchovým stavem je modelový poruchový stav části stroje.
    -31 CZ 2018 - 517 A3
  13. 13. Diagnostický systém strojů podle kteréhokoliv z nároků 6 až 12 vyznačující se tím, že řídicí část zahrnuje monitorovací jednotku a vzdálený server, přičemž senzory jsou datově spojeny s monitorovací jednotkou a monitorovací jednotka je datově spojena se vzdáleným serverem.
    5
  14. 14. Diagnostický systém stroje podle kteréhokoliv z nároků 6 až 13 vyznačující se tím, že alespoň jedním ze senzorů je snímač vibrací a alespoň jedním dalším senzorem je snímač polohy.
  15. 15. Diagnostický systém stroje podle kteréhokoliv z nároků 6 až 13 vyznačující se tím, že alespoň jedním ze senzorů je snímač vibrací a alespoň jedním dalším senzorem je snímač ίο přetvoření.
  16. 16. Diagnostický systém stroje podle kteréhokoliv z nároků 6 až 13 vyznačující se tím, že alespoň dvěma senzory jsou snímače vibrací.
    15
  17. 17. Diagnostický systém stroje podle kteréhokoliv z nároků 6 až 16 vyznačující se tím, že senzory jsou tenko vrstvě senzory.
  18. 18. Diagnostický systém strojů podle kteréhokoliv z nároků 6 až 17 vyznačující se tím, že snímačem polohy je lineární enkodér nebo rotační enkodér.
  19. 19. Diagnostický systém strojů podle kteréhokoliv z nároků 6 až 18 vyznačující se tím, že alespoň jedním ze senzorů je snímač teploty, přičemž snímač teploty je umístěn v těsné blízkosti alespoň jednoho senzoru jiného typu.
CZ2018-517A 2018-09-30 2018-09-30 Diagnostický systém strojů CZ2018517A3 (cs)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CZ2018-517A CZ2018517A3 (cs) 2018-09-30 2018-09-30 Diagnostický systém strojů
SK50075-2018A SK500752018A3 (sk) 2018-09-30 2018-12-21 Diagnostický systém strojov
US17/280,412 US11543329B2 (en) 2018-09-30 2019-09-26 Diagnostic system of machines
PCT/CZ2019/050039 WO2020064033A1 (en) 2018-09-30 2019-09-26 Diagnostic system of machines

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CZ2018-517A CZ2018517A3 (cs) 2018-09-30 2018-09-30 Diagnostický systém strojů

Publications (1)

Publication Number Publication Date
CZ2018517A3 true CZ2018517A3 (cs) 2020-04-08

Family

ID=68653293

Family Applications (1)

Application Number Title Priority Date Filing Date
CZ2018-517A CZ2018517A3 (cs) 2018-09-30 2018-09-30 Diagnostický systém strojů

Country Status (4)

Country Link
US (1) US11543329B2 (cs)
CZ (1) CZ2018517A3 (cs)
SK (1) SK500752018A3 (cs)
WO (1) WO2020064033A1 (cs)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117784710B (zh) * 2024-02-23 2024-06-04 江西秦锐智能科技有限公司 数控机床远程状态监控***及其方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6539343B2 (en) * 2000-02-03 2003-03-25 Xerox Corporation Methods for condition monitoring and system-level diagnosis of electro-mechanical systems with multiple actuating components operating in multiple regimes
US20100284546A1 (en) 2005-08-18 2010-11-11 Debrunner Victor Active noise control algorithm that requires no secondary path identification based on the SPR property
CN102308194B (zh) * 2008-12-22 2014-10-22 S.P.M.仪器公司 分析***
ITCO20120008A1 (it) 2012-03-01 2013-09-02 Nuovo Pignone Srl Metodo e sistema per monitorare la condizione di un gruppo di impianti
WO2015187220A2 (en) 2014-05-07 2015-12-10 Sikorsky Aircraft Corporation Rotor system structural fault estimation

Also Published As

Publication number Publication date
US20220034756A1 (en) 2022-02-03
US11543329B2 (en) 2023-01-03
SK500752018A3 (sk) 2020-04-02
WO2020064033A1 (en) 2020-04-02

Similar Documents

Publication Publication Date Title
EP3629118B1 (en) Diagnostic system of machine tools
Chen et al. CPS modeling of CNC machine tool work processes using an instruction-domain based approach
CN104808585B (zh) 一种机床健康状态快速检查方法
Fines et al. Machine tool positioning error compensation using artificial neural networks
Diez et al. Feasibility study of in-process compensation of deformations in flexible milling
CN107111297A (zh) 用于由至少一台计算机数控机器加工的工件的部件分析的计算机实现方法
El Badaoui et al. Modeling and detection of localized tooth defects in geared systems
CN108803499A (zh) 控制装置以及机器学习装置
Bauer et al. Analysis of industrial robot structure and milling process interaction for path manipulation
Indri et al. Friction modeling and identification for industrial manipulators
Liu et al. In-process identification of machine tool dynamics
D'Emilia et al. Mechatronics applications of measurements for smart manufacturing in an industry 4.0 scenario
Sparham et al. ANFIS modeling to predict the friction forces in CNC guideways and servomotor currents in the feed drive system to be employed in lubrication control system
Peng et al. Prediction of milling force based on spindle current signal by neural networks
Zhao et al. Robotic grinding process monitoring by vibration signal based on LSTM method
CZ2018517A3 (cs) Diagnostický systém strojů
CN104699905B (zh) 基于频域响应的调速***齿轮传动机构辨识建模方法
EP3856507B1 (en) Diagnostic system of forming machines
Benker et al. Experimental derivation of a condition monitoring test cycle for machine tool feed drives
Castaño et al. Application of hybrid incremental modeling for predicting surface roughness in micromachining processes
Aslan Integration of virtual and on-line machining process control and monitoring using cnc drive measurements
Giovannitti et al. A virtual sensor for backlash in robotic manipulators
Sicard et al. Design Considerations for Building an IoT Enabled Digital Twin Machine Tool Sub-System
Lu et al. In-process surface quality monitoring of the slender workpiece machining with digital twin approach
Čereška Vibrodiagnostics and monitoring of the mechanical-dynamic elements of mechatronic systems