CZ303481B6 - Zpusob stanovení míry opotrebení prevodovky a zarízení k provádení tohoto zpusobu - Google Patents
Zpusob stanovení míry opotrebení prevodovky a zarízení k provádení tohoto zpusobu Download PDFInfo
- Publication number
- CZ303481B6 CZ303481B6 CZ20110665A CZ2011665A CZ303481B6 CZ 303481 B6 CZ303481 B6 CZ 303481B6 CZ 20110665 A CZ20110665 A CZ 20110665A CZ 2011665 A CZ2011665 A CZ 2011665A CZ 303481 B6 CZ303481 B6 CZ 303481B6
- Authority
- CZ
- Czechia
- Prior art keywords
- load
- gearbox
- test
- particles
- gear
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000012360 testing method Methods 0.000 claims abstract description 71
- 239000002245 particle Substances 0.000 claims abstract description 58
- 238000001228 spectrum Methods 0.000 claims abstract description 17
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 13
- 229910052742 iron Inorganic materials 0.000 claims abstract description 6
- 230000002441 reversible effect Effects 0.000 claims abstract description 3
- 230000005540 biological transmission Effects 0.000 claims description 10
- 238000013517 stratification Methods 0.000 claims description 9
- 238000000611 regression analysis Methods 0.000 claims description 5
- 238000005299 abrasion Methods 0.000 abstract description 16
- 238000004458 analytical method Methods 0.000 abstract description 2
- 230000001186 cumulative effect Effects 0.000 abstract 2
- 239000003921 oil Substances 0.000 description 36
- 238000005070 sampling Methods 0.000 description 16
- 238000004519 manufacturing process Methods 0.000 description 12
- 238000002474 experimental method Methods 0.000 description 11
- 238000000926 separation method Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 239000000314 lubricant Substances 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000011109 contamination Methods 0.000 description 3
- 239000013527 degreasing agent Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000001479 atomic absorption spectroscopy Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 101100459912 Caenorhabditis elegans ncs-1 gene Proteins 0.000 description 1
- 241001442234 Cosa Species 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000005237 degreasing agent Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000003517 fume Substances 0.000 description 1
- 239000012208 gear oil Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000007970 homogeneous dispersion Substances 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000009349 indirect transmission Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 239000006249 magnetic particle Substances 0.000 description 1
- 238000007885 magnetic separation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000011022 opal Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Míra opotrebení prevodovky (5) podle obsahu oterových cástic v olejové náplni se stanoví tak, že se nejprve v reálném provozu vozidla zjistí zátežové spektrum pro každý prevodový stupen, tj. pocet otácek vykonaných tímto stupnem v po sobe následujících intervalech hodnoty krouticího momentu prenášeného tímto stupnem vztaženo na jednotku dráhy vozidla, címž se získá pomerná hodnota práce prenesené stupnem v jednotlivých intervalech - momentovýchhladinách. Pro jednotlivé stupne se sestaví náhradní zátežové bloky tak, že se vybere jedna kladná a jedna záporná momentová hladina s nejvetší prevedenou prací a do ní se pripocte práce všech ostatních momentových hladin se stejným znaménkem. Sestava takto vytvorených dílcích zátežových bloku pro jednotlivé stupne vytvorí souhrnný zátežový blok pro prevodovku, urcí se délka souhrnného zátežového bloku odpovídající zvolené dráze vozidla nebo prenesené práci, na zkušebním stendu s reverzibilním pohonem a brzdným zarízením se testovaná prevodovka provozuje tak, že se pri konstantních otáckách v hodnote 0,5 až 0,8 maximálních otácek podrobuje zatížení po sobe následujícími zátežovými bloky, pricemž mezi jednotlivými zátežovými bloky, vždy ihned po jejich ukoncení, se z prevodovky odeberou vzorky olejové náplne, které se podrobí analýze ke zjištení koncentrace cástic železa v jednotlivých fázích zkoušky. Príslušné zarízení je tvoreno stejnosmerným hnacím motorem (1) s tyristorovým regulátorem (2) otácek, redukcní prevodovkou (6), zkušebním stojanem (4) s hnacím vretenem, na který se upíná zkoušená prevodovka (5) a dále vírivým dynamometrem (7) propojeným s regulátorem (8) momentu.
Description
Způsob stanovení míry opotřebení převodovky a zařízení k provádění tohoto způsobu
Oblast techniky
Vynález se týká způsobu stanovení míry opotřebení převodovky podle obsahu otěrových částic v olejové náplni a zařízení k provádění tohoto způsobu.
Dosavadní stav techniky
Provozní zkoušky převodovek testují skutečný agregát v běžném provozu. Je to spojeno se značnými provozními náklady při poměmě nízké zkušební produktivitě. Další nevýhodou provozních zkoušek je při jejich omezeném počtu nízká statistická významnost výsledků. Chyba může být vnesena nejen rozptylem provozních podmínek, ale zejména rozptylem vstupních technologických parametrů. Při sledování opotřebení může například rozdíl v počátečních podmínkách způsobený nepřesností výroby, tepelného zpracování i montáže překrýt očekávaný vliv různých maziv.
Oproti tomu laboratorní modelové zkoušky využívají typizované převodovky (např. FZG). Jejich hlavní výhodou je vysoká produktivita a velmi dobrá reprodukovatelnost výsledků. Systém zkoušky většinou umožňuje sledovat závislost výsledku na změně jediného vstupního parametru. Pozitivně se zde projevuje možnost zkrácení zkoušky preferováním té části zátěžového spektra, která sledované účinky ovlivňuje nejvíce. Otázkou je, jak dalece odpovídají zkušební podmínky standardizovaných zkoušek provozním podmínkám skutečným. Vynález si klade za úkol navrhnout způsob laboratorní zkoušky opotřebení převodovky, který by co nejvíce odpovídal reálným provozním podmínkám a přinášel výsledky s maximální výpovědní schopností.
Podstata vynálezu
Uvedený úkol řeší způsob stanovení míry opotřebení převodovky podle obsahu otěrových částic v olejové náplni, jehož podstata spočívá v tom, že se nejprve v reálném provozu vozidla zjistí zátěžové spektrum pro každý převodový stupeň, tj. počet otáček vykonaných tímto stupněm v po sobě následujících intervalech hodnoty kroutícího momentu přenášeného tímto stupněm vztaženo na jednotku dráhy vozidla, čímž se získá poměrná hodnota práce přenesené stupněm v jednotlivých intervalech - momentových hladinách, načež se sestaví pro jednotlivé stupně náhradí zátěžové bloky tak, že se vybere jedna kladná a jedna záporná momentová hladina s největší převedenou prací a do ní se připočte práce všech ostatních momentových hladin se stejným znaménkem, přičemž sestava takto vytvořených dílčích zátěžových bloků pro jednotlivé stupně vytvoří souhrnný zátěžový blok pro převodovku, určí se délka souhrnného zátěžového bloku odpovídající zvolené dráze vozidla nebo přenesené práci, na zkušebním stendu s reverzibilním pohonem a brzdným zařízením se testovaná převodovka provozuje tak, že se při konstantních otáčkách v hodnotě 0,5 až 0,8 maximálních otáček podrobuje zatížení po sobě následujícími zátěžovými bloky, přičemž mezi jednotlivými zátěžovými bloky, vždy ihned po jejich ukončení, se z převodovky odeberou vzorky olejové náplně, které se podrobí analýze ke zjištění koncentrace částic železa v jednotlivých fázích zkoušky.
Ke zpřesnění výsledků se koncentrace částic železa ve vzorcích koriguje započtením poměrné části částic odebraných v předchozích vzorcích a částic usazených na stěnách převodovky a rovněž se zohlední rozvrstvení částic v olejové náplni.
Ke zjištění středního nárůstu koncentrace částic železa během zkoušky je výhodné podrobit hodnoty koncentrace regresní analýze.
- 1 CZ 303481 B6
Zařízení k provádění způsobu podle nároku 1 je tvořeno stejnosměrným elektromotorem s tyristorovou regulací otáček a převodovou skříní a dále vířivým dynamometrem propojeným s regulátorem momentu.
Přehled obrázků na výkresech
Vynález bude dále objasněn pomocí výkresu, na němž obr. 1 představuje schéma sestavení zátěžového bloku pro 3. stupeň příkladné převodovky, obr. 2 je vývojový diagram iteračního procesu ío ke korekci zohledňující rozvrstvení částic v olejové náplni, obr. 3 představuje porovnání vývoje koncentrace částic železa zjištěné ve vzorcích s vypočtenou průměrnou koncentrací v náplni, na obr. 4 je průběh koncentrací částic podrobený regresní analýze a na obr. 5 je schéma zkušebního zařízení.
Příklady provedení vynálezu
Způsob stanovení míry opotřebení podle vynálezu byl aplikován na pětistupňovou převodovku osobního automobilu. Na základě zkušeností byla pro test převodovky určena délka celé zkoušky odpovídající ujetí 25 000 km. Jako podklad pro sestavení náhradního zatížení bylo použito jednotkové otáčkové spektrum vytížení hnacího ústrojí při jízdě automobilu. Toto spektrum bylo získáno při jízdních zkouškách na zkušební dráze a představuje počty otáček v jednotlivých momentových intervalech/hladinách přepočtené na 1 km zkušební dráhy. V tabulce 1 je uvedeno toto spektrum pro každý z 5 převodových stupňů rozdělené do 17 momentových hladin podle velikosti vstupního momentu převodovky.
| Hl. | M í Nm 1 | 1 . p- st. | 2.p.st. | 3. p.et. | 4. p.st. | 5.p.xt. |
| 17 | teo. o | O. OO | 0. oc | 0. 13 | . 0. 20 | 0. 57 |
| 15 | 144.0 | 0. Ol | 0. 43 | 0. 87 | 0. S2 | 0. 11 |
| 15 | 13®. O | 0.13 | O. 34 | 0. 07 | 0. 74 | O. 23 |
| 14 | 112.0 | 1.02 | 0. 8Θ | 5. 22 | 2. OS | 2. 25 |
| 13 | 9«. O | 3. 01 | 35. oa | 213. 42 | 85. 05 | 152. 41 |
| 12 | 30.0 | 12. 24 | 110.25 | 23S. 85 | 24. 40 | ISO. Sl |
| 11 | 54.0 | ie. 74 | 48. 07 | 13- 31 | 5. 15 | 40. 31 |
| 10 | 43. 0 | 4. 31 | 33. 14 | 5. 02 | 70. 51 | 1 4. 24 |
| 9 | 33. O | 2. 33 | 32.25 | 13. SS | 53. 14 | 11.44 |
| 9 | ie. o | 3. 55 | 37. 55 | 236.11 | 72. 12 | 22. 00 |
| 7 | 0.0 | 3. 34 | 52. 47 | 30. 07 | 95. 94 | 5, 75 |
| S | -1 e. o | 12. 03 | 65. 44 | 117. 84 | 92. 07 | 4. 52 |
| 3 | -32. 0 | 1 . £7 | 1S. 72 | 8. 15 | 32. 94 | O. 13 |
| 4 | -43. 0 | O. 03 | 2.S4 | 1. 17 | 27. 10 | 0. OO |
| 3 | -54. 0 | 0. 00 | 0. 25 | 2. 32 | 3. 24 | 0. 02 |
| a | -30. 0 | O. OO | 0. OO | O. 15 | 0. 05 | O. 03 |
| 1 | -05, 0 | 0.00 | o, oo | O. 12 | 0. 00 | 0. 00 |
| Sux» | otáček | Ol. 35 | 300. 47 | 933. 83 | 505. 40 | 433. 53 |
| Dráha íkml | 0. OOO | 0.114 | 0. 32S | O. 288 | O. 254 | |
| Čas | C s! | 1.234 | 7.534 | 13. 014 | 11.100 | G. 409 |
Tab. 1
-2CZ 303481 B6
Následně byl ze zátěžových spekter jednotlivých převodových stupňů sestaven pro každý stupeň zátěžový blok. Na obr. 1 je graficky znázorněno zátěžové spektrum pro 3. převodový stupeň podle hodnot z tabulky 1. Na ose x jsou vstupní kroutící momenty na převodovce resp. sled momen5 tových hladin. Šířka sloupce histogramu je dána zvolenou hodnotou hladiny, tedy intervaly vstupního kroutícího momentu a výška poetem otáček na 1 km dráhy. Plocha sloupce tedy představuje práci přenesenou na dané momentové hladině, tj. v daném rozmezí momentů.
Rychlostní stupeň není namáhán jen při akceleraci a ustáleném pohybu vozidla - kladná část io spektra, nýbrž i při brzdění motorem - záporná část spektra. Při sestavování zátěžového bloku je třeba zohlednit obě části spektra, takže zátěžový blok bude mít dvě části. Zátěžový blok je znázorněn na obr. 1 šrafovanými plochami. Dospělo se k němu tak, že se vybrala jedna kladná a jedna záporná momentová hladina s největší převedenou prací a do ní se připočetla práce všech ostatních momentových hladin se stejným znaménkem. Zátěžový blok znázorněný šrafovanou oblastí tak představuje náhradu práce celého spektra přepočtenou na hladinu resp. hladiny přenášející největší procento celkové práce v kladném i záporném smyslu.
Zátěžový blok navržený tímto postupem je tedy pro každý převodový stupeň dvouhladinový, jedna hladina je s kladnou a jedna se zápornou hodnotou kroutícího momentu. Protože práce převe20 děná prvním převodovým stupněm představuje pouze 2,5 % z celkové převedené práce, je tento převod ze zátěžového bloku vyřazen (H-l, H+l). Ze stejných důvodů je vypuštěna záporná hladina pátého převodového stupně (H-5). Řazením zátěžových bloků jednotlivých rychlostních stupňů za sebou byl vytvořen souhrnný zátěžový blok pro celou převodovku.
Velký význam pro urychlení zkoušky má volba zkušebních otáček. S ohledem na růst vnitřních dynamických sil byly jako mezní přijaty otáčky rovnající se 3/4 maximálních otáček motoru. V konkrétním případě testované převodovky to představuje frekvenci otáčení 3900 n/min.
Bylo také nutno vyřešit problém délky (celkové práce) jednoho zátěžového bloku. Snaha o získá30 ní co největšího množství dat kjedné zkoušce vede ke stanovení poměrně krátkých zátěžových bloků a vyšší hustotě odebírání vzorků oleje. To ovšem vyvolává značné nároky na obsluhu, protože zkušební doby na jednotlivých hladinách jsou poměrně krátké. Snaha o kompromis v tomto směru vedla k vytvoření zátěžového bloku odpovídajícího 2500 km, což podle naměřeného otáčkového spektra odpovídá převedení práce ve výši asi 1900 MJ, Jeho rozpis pro zkoušenou pře35 vodovku je v tabulce 2,
| Hladina bi oku | Hladi na ot. sp. | t Nm 1 | prev. st. | Čas l h 1 | Práce t MJ I |
| H + 2 | 12 | SO | 2 | 2. 5 | 294 |
| .4+4 | 13 | 95 | 4 | 2. O | 2se |
| H+S | 13 | 93 | 5 | 3. 5 | soe |
| H-4 | Ξ | -32 | 4 | l.S | 71 |
| K~2 | S | -16 | 2 | 1.0 | 24 |
| H-3 | 6 | -16 | 3 | 1.5 | 35 |
| H*3 | 13 | 96 | 3 | s. o | 706 |
Tab. 2
Popsaný postup při sestavení zátěžového bloku vychází z těchto poznatků a úvah:
-3 CZ 303481 B6
Převodovka je na zkušebním stavu zatěžována při konstantních otáčkách cyklicky se opakujícími zátěžovými bloky, jejichž účinek je adekvátní v provozu zjištěnému momentovému spektru. Hodnota zkušebních otáček má pochopitelně velký význam pro urychlení zkoušky. S ohledem na růst vnitřních dynamických sil však nemohou být použity větší zkušební otáčky, než jsou maximální provozní. Urychlení zkoušky lze dosáhnout při sestavování náhradního zatížení tím, že se přihlédne k fyzikální podstatě vzniku opotřebení ozubených kol. Metoda výpočtu opotřebení ozubených kol s pomocí hodnoty lineární intenzity opotřebení A je rozpracována v [1]. Mírou relativního prokluzování profilů zubů jsou tzv. měrné skluzy 3, které se rovnají poměru skluzové rychlosti v bodě doteku k tečným složkám obvodových rychlostí v tomto bodě. Tloušťka opotřebené vrstvy u soukolí ze stejného materiálu za dobu t potom bude:
ho — . & 'Ppr O)» kde w'n, je měrná výpočetní síla na bok zubu, E modul pružnosti, v je hodnota Poissonovy konstanty ppr je vstupní poloměr křivosti kontaktních ploch, n je frekvence otáčení a i je součinitel vícepárového záběru kola. Měrnou výpočetní sílu na bok zubu lze vyjádřit vztahem:
b^áw,.cosav.<x>sfi' (2).
kde Mt je hnací moment, Kha součinitel podílu zatížení jednotlivých zubů, Kup součinitel nerovnoměrnosti zatížení zubu, K!!v součinitel vnitřních dynamických sil, bw pracovní šířka ozubení, dw! valivý průměr hnacího kola, valivý záběrový úhel v čelní rovině a úhel sklonu ozubení. Při konstrukci náhradního zatížení, tzv. zátěžového bloku, je cílem sestavit při konstantních vstupních otáčkách takový sled zátěžových hladin, který způsobí v co nej kratší době stejné opotřebení převodovky jako reálný zátěžový cyklus. Jako podklad je vhodné použít naměřené jednotkové otáčkové spektrum zatížení převodovky. Například u hnacího ústrojí osobního automobilu může být toto spektrum získáno při jízdních zkouškách na vhodně zvolené zkušební dráze [2]. Při náhradě reálného zatížení zkušebním blokem lze předchozí vztahy formálně zjednodušit na tvar:
ho
(3).
kde Kj představuje obecně závislost na kinematice záběru. Součinitel A lze vyjádřit jako:
HB
K2> (4), kde HB]q tvrdost boku zubu podle Brinella, K2 vyjadřuje závislost na parametrech charakterizujících kvalitu mazání a σΗ je kontaktní tlak v ozubení:
(5), když K3 vyjadřuje závislost kontaktního tlaku na geometrických poměrech v kontaktu. Je-li potom K4 funkce nezávislá na Mt, n a t lze výše uvedené shrnout do vztahu:
= Ki-K2-KíMiM = Κ.,Μ,Μ = k4.a. HB (6).
-4CZ 303481 B6
Opotřebení převodu je tedy úměrné práci A převodem převáděné. Náhradní zátěžový blok byl sestaven tak, aby pro každý převodový stupeň (u vícestupňových převodovek) obsahoval pouze dvě zátěžové hladiny při konstantních otáčkách (jednu v zátěžovém a druhou v brzdném smyslu zatížení). Vybrána je vždy momentová hladina s největší převedenou prací, do které se připočte práce všech ostatních momentových hladin. To velmi urychlí zkoušku, protože odpadnou dlouhé fáze s malým hnacím momentem. Postupné střídání těchto bloků pro všechny převodové stupně tak bude působit stejné opotřebení všech ozubených kol jako reálný provoz převodovky.
Znečištění oleje otěrovými částicemi při zkoušce způsobují především ozubení kola a tvoří je v naprosté většině částice železa. Protože jsou zátěžové bloky krátké, musí být metoda citlivá, ale z provozních důvodů nesmí být příliš nákladná. Ideální je proto analyzátor pracující na principu atomové absorpční spektrometrie (AAS), který je běžnou součástí tribotechnických laboratoří. Výhodou je, že hořák analyzátoru spaluje i směs s poměrně velkými otěrovými částicemi (<10pm).
Vždy po absolvování jednoho zátěžového bloku je proveden vzorkovací odběr. Množství odebíraného oleje je důležitou volbou. Za odebíraný objem je nutno neustále doplňovat čerstvý olej, takže je tímto procesem narušeno modelování skutečných provozních podmínek. Je tedy snaha toto množství minimalizovat. Na druhé straně však metody analýzy exploatovaného oleje vyžadují konkrétní minimální objem.
Velmi významně závisí výsledek zkoušky na zjištění vzorků oleje po ukončení zkoušky. Celá náplň je po odebrání posledního vzorku vypuštěna z převodovky, poté promíchána a odebrán vzorek. Vypouštění má být prováděno se snahou, aby v převodovce zbylo co nejméně reziduálního oleje. Ještě za provozní teploty je skříň naplněna odmaŠťovadlem a za pomalého protáčení proplačhována. Po vypuštění je skříň demontována a vymyta znovu v lázni odmašťovadla. Po odpaření odmašťovadla v digestoři jsou tyto úsady naředěný známým množstvím čistého oleje a rovněž analyzovány na otěr.
Nejjednodušším modelem experimentu by bylo po odejití jistého počtu zátěžových cyklů podrobit analýze vypuštěnou olejovou náplň i vymyté ůsady. V tomto případě však nelze ve výsledcích podchytit různou intenzitu tvorby otěrových částic v jednotlivých etapách zkoušky (např. fázi záběhu). V případě havárie převodovky v závěru experimentu jsou prakticky znehodnoceny výsledky časově i ekonomicky velmi nákladného procesu. Z těchto důvodů bylo přistoupeno ke vzorkování olejové náplně v průběhu experimentu. Pravidelné vzorkování olejové náplně také poskytuje vyšší počet naměřených hodnot a zvyšuje tak statistickou významnost vyhodnocení zkoušky, kterou často nelze opakovat.
S ohledem na skutečnost, že analýza odebraných vzorků popisuje pouze stav olejové náplně, nikoliv usazeniny a odběry předchozích vzorkuje nutno v metodice zohlednit usazování a odebírání otěrových částic. Základními příčinami odlučování otěru z oleje jsou usazovací procesy na pohyblivých a pevných dílech převodovky, mající za následek přilnutí částic zejména na nepohyblivých dílech. Magnetická separace se při experimentu neprojevuje, protože magnetický sepa45 rátor částic musí být před zkouškou vždy odstraněn.
Pro stanovení závislosti mezi generovaným množstvím otěrových částic a jeho dispergovaného podílu je potřebné definovat následující zjednodušující předpoklady:
· otěrové částice mezi sebou navzájem nereagují, tj. nevytvářejí se chemické sloučeniny ani aglomeráty částic, • neuvažují se ztráty v náplni (únik netěsnostmi, opal atd.) exploatovaného oleje,
-5CZ 303481 B6 • výše uvedené procesy odlučování částic jsou proporcionální počtu částic v mazivu a tato proporce je neměnná v průběhu celé zkoušky, • ve výše uvedeném smyslu je předpokládáno chování částic všech sledovaných velikostí «5 μηη).
Je-li za časový úsek Δτ vyprodukováno množství částic o stejné hmotnosti mx a z původní hmotnosti m: částic v oleji se odloučil ar-tý díl, je nová hmotnost částic dispergovaných v oleji za úsek Δγ.
m} = m,.(l -a) + mx (7).
Bude-li produkce otěrových částic při konstantním momentovém zatížení neměnná, bude v čase /.Jr množství volného otěru v oleji:
(8).
Za předpokladu, že zkouška nebude narušena vzorkováním, je možno tedy ze znalosti dispergovaného množství otěru v oleji po zkoušce počátečního znečištění m2 a množství vymytého ve formě úsad tnu vyjádřit hmotnost vyprodukovaného otěru:
(9).
Množství vyprodukovaných částic, nacházejících se neusazených v oleji bude:
Nyní je možné stanovit intenzitu produkce částic mx za časový úsek Δτ.
(10).
(11) a dosud neznámý součinitel odlučování a\
^.(7-(/-α)ϋ (12) a ~ -35
Tento implicitní vztah lze řešit numerickou iterací a je jádrem analytického programu, který simuluje rozvoj opotřebení při zkouškách, které se mohou výrazně lišit v časové ose. Pokud nastane situace, že není známa hmotnost usazených (vymytých) částic opotřebení, lze na základě znalosti součinitele a z předchozích experimentů stanovit:
.m (13)
Zavedeme-li z praktických důvodů pojmy průměrné koncentrace otěrových částic N = m/Vf, a měrné intenzity zvyšování koncentrace částic X = kde VL je celkový objem oleje v pře45 vodovce, budou platit výše uvedené vztahy i po naměřené koncentrace (s tím rozdílem, že hodno-6CZ 303481 B6 ta X bude mít jiný rozměr než mx). Na je v tomto případě fiktivní koncentrace představující hmotnost zvážených úsad přepočtenou na objem olejové lázně.
Při vzorkování znečištění olejové náplně je odebráno reprezentativní množství Vv oleje a je 5 zjišťována koncentrace nečistot N. Olejový systém převodovky je po odběru ihned doplněn stejným množstvím čistého maziva. Celková koncentrace znečištění v převodovce N tedy poklesne na hodnotu Ν':
/^ = ^/-^ = ^.(/-^), (14),
IO kde iv je poměr objemu vzorku a_celkového objemu. Bylo-li v průběhu zkoušky provedeno q odběrů o různých koncentracích NL, je nutno toto odebrané množství otěru zohlednit jak v hodnotě vyprodukované koncentrace otěru:
_ _ _ _ q-i_ _ _ _
Ν'=Ν,+Ν„- N; + /,..£ Nl =N, + Nu- N, + /K.5„ (15), tak jejího dispergovaného podílu:
Nl = N, + (7-a)ř (16), kde Z/, reprezentuje čas odběru vzorku měřený od počátku zkoušky. Další postup stanovení X a součinitele a je stejný.
Během zkoušky není intenzita produkce otěrových částic stále konstantní. Výrazně se mění v závislosti na vstupním točivém momentu při postupném střídání jednotlivých momentových hladin zátěžového bloku. Svůj vliv na proměn nost produkce otěru uplatní také různé provozní podmínky v průběhu zkoušky (záběh, havarijní stav). Má-li zůstat měrná intenzita zvyšování koncentrace částic X po celý experiment nezávislá na těchto vlivech, je nutno reálnou intenzitu generování otěru v každém úseku zkoušky Xr vyjádřit ve tvaru:
Χ,=Μ^ζ2ί, (17), kde Mi je hodnota hnacího momentu a mění se plánovitě několikrát mezi následujícími odběry, zatímco k. je parametr charakterizující opotřebení mezi dvěma následujícími odběry a v tomto úseku bude pro zjednodušení považován za konstantní. Po zavedení indexování m, označující chronologické poradí zátěžové hladiny zvyšující svou hodnotu při změně této hladiny, a index L, označující chronologické pořadí odběru. Index m dosahuje hodnot l ažp, a L se pohybuje v intervalu 7 až q. Označíme-li odpovídající hodnoty Mtm a kzi. lze za předpokladu, že tp = t() = /, odvodit:
NcS^l-a)' s4.n.
(18) (19)
-7CZ 303481 B6 pro známé Nea 7Vř/. Známe-li pouze N(l a parametr a, potom:
a.Nd
Ss.(J-dť (20), kde hodnoty X a S4 zohledňují proměnlivost produkce otěru:
M^P
L-J ^4 ~ ~ m)
L=-i (21) (22).
Je nutno uvažovat s přesně danou časovou souvislostí mezi indexováním m a L. Po ukončení ío experimentu nejsou známy hodnoty součinitelů k7L. V prvním krokuje tedy nutno při výpočtu X a parametru a použít odhadnutých hodnot (£z/ = kZ2 — = kZg = 1) a poté přistoupit k vyčíslení součinitelů kZL na základě naměřených koncentrací odebíraných vzorků.
Jsou-li v časech daných tL atL,} vzorkováním naměřeny koncentrace NL a NL h lze při znalosti 15 a, Xa časového průběhu změn M/m (v časech daných tm začínají momentové hladiny o Mfm) stanovit hodnotu kzi, mezi těmito odběry:
a_ N^-NL.il-iMl-a)'··^ (23) h (1-a)'-----20 Nyní je možno znovu přikročit k opětnému stanovení parametrů a a X. Iterační cyklus je nutno opakovat do splnění cílových podmínek. Je vhodné, vzhledem k závislosti X na a, dodatečně upravit po každém výpočtu součinitelů jejich velikost na hodnoty k'Zi, = c. kZL, kde konstanta c je volena tak, aby při iteračních krocích docházelo pouze ke změnám parametru a (X zůstává neměnné). Má-li zůstat měrná produkce otěrových částic Y po každé iteraci stejná, musí být koe25 fícienty Az/j upraveny tak, aby i součet násobků kZL. \ί!η,.Λτ pro celou zkoušku (představující celkovou převedenou práci) byl stále stejný. Konstantu c proto vypočteme ze vztahu:
m=p
ΣμιΛι.., -d m=í (24)
L-l
-d
Za provozu převodovky dochází k výraznému rozvrstvení otěrových částic v olejové náplni. Jde o dynamické rozvrstvení způsobené prouděním maziva ve skříni. Protože vzorek oleje je odebírán vzorkovacím kohoutem vždy ze stejného prostoru skříně, je potřebné znát závislost mezi koncentrací částic N v tomto prostoru a průměrnou koncentrací Ϋ v celé náplni. Po myšlenkovém rozdělení olejového prostoru na k vrstev byl proto do matematického modelu zaveden souči-8CZ 303481 B6 nitel rozvrstvení b, který vyjadřuje, jaký podíl částic opotřebení z celkového množství ve vrstvě přestoupí za časový úsek Jrdo vrstvy následující.
V každé vrstvě se tedy za tento okamžik zvyšuje koncentrace o množství otěru vytvořeného třecí5 mi kontakty a o množství infiltrované ze sousední vrstvy. Současně však nastává úbytek vlivem odlučování a přechodem částic do vrstvy sousedící na druhé straně. Koncentraci y—té vrstvy v čase daném /-tým násobkem Δτ lze vyjádřit vztahem (pro
N,j = -a-b) + + X (25).
Pro stanovení koncentrace pro j=l a j=k použijeme logických okrajových podmínek. První vrstva je charakterizována tím, že do ní nepřechází částice z žádné jiné vrstvy; bude tedy:
N,,i = N,_,,.(l-a-b) + X (26)
Poslední vrstva se vyznačuje tím, že z této oblasti již nemají částice kam přecházet a proto:
S.^N^b + X (27).
Problém se řeší zejména z toho důvodu, že koncentrace zjištěná ve vzorku odebraném v čase zdaleka nemusí odpovídat průměrné koncentraci v tomto okamžiku Af. Za předpokladu, že na počátku zkoušky je disperze částic dostatečně homogenní (tedy N: = N2), bude se rozdíl mezi nimi postupně zvětšovat. Většinou jsou známé obě hodnoty koncentrace pri posledním odběru po ukončení zkoušky, kdy se průměrná koncentrace A^ vypočte z koncentrace odebraného vzorku
N(, a koncentrace homogenizované vypuštěné lázně N'q:
N,=i,.Nq+(l-iv)^\ (28).
Vzhledem k tomu, že do výše uvedených vztahů pro zjištění intenzity produkce Částic X a odlu30 čovacího parametru a je nutno dosazovat hodnoty průměrných koncentrací, musí se tyto hodnoty stanovit. Byla zvolena metoda postupného modelování rozvrstvování za stanovených podmínek s krokem Δτ. V první fázi se vychází z hodnot a, kZi. a X stanovených ze zjednodušeného předpokladu, že koncentrace je vždy homogenní. Ze získaných výpočtových hodnot průměrných a vzorkovaných koncentrací ~ŇLv a M, a známé skutečné vzorkované koncentrace Nfj se úměrou stanoví vstupní hodnoty průměrných koncentrací Nf. do dalšího iteračního kroku:
y<”'> = ,v — L ν'·νϊ>’ (29), kde <v>, <w>, <u> jsou příznaky kroku iterace. Pozornost je nutno věnovat také stanovení 40 hodnoty součinitele rozvrstvení b. Pro první iterační krok je volen odhadem a pro jeho korekci v průběhu iteračního cyklování se použije vztah:
£<w+/> = b(w\^
N„.Nqv
Na.N (30).
qv
Cílem je tedy dosáhnout stejného rozložení výpočtové a naměřené koncentrace k okamžiku ukončení zkoušky. Při modelování vzorkování olejové náplně a s ním spojeném poklesu koncentrace
-9CZ 303481 B6 předpokládáme, že otěrové částice jsou odebírány ze všech vrstev, úměrně jejich zastoupení. Kritériem pro ukončení celého iteraěního procesuje stabilizace součinitele rozvrstvení b.
Za situace, kdy z nej různějších důvodů není k dispozici koncentrace vypuštěné náplně Níf (např. ? havárie, průběžné vyhodnocení experimentu) a v předcházejících zkouškách byla ověřena hodnota parametru b pro stejné podmínky, lze výpočet realizovat náhradním způsobem. Na počátku iterace je za tuto hodnotu dosazena koncentrace posledního vzorku Nif, v dalších krocích je tato hodnota korigována pomocí vztahu:
A v,.
(31)
Iterace je ukončena požadovanou stabilizací této hodnoty.
Aplikací výše uvedených vztahů vznikl analytický program, který na základě výsledků zjištěných 15 po ukončení experimentu (koncentrace vypuštěné náplně, koncentrace naředěných úsad atd.) a údajů získaných vzorkováním umožňuje odhadnout skutečnou intenzitu produkce otěrových částic Xr v libovolném okamžiku zkoušky. Vývojový diagram programuje na obr. 2.
Na obr. 3 je ukázka modelování vývoje koncentrace Fe částic v olejové náplni. Naměřené vzor20 kované hodnoty jsou vyznačeny kroužky, a čarami jsou znázorněny průběhy vypočtených průměrných koncentrací - plná čára 1 a koncentrace v blízkosti vzorkovacího kohoutu - Čárkovaná čára 2.
Jsou-li známy hodnoty intenzity produkce otěru Xr ve všech úsecích zkoušky a chronologie zatě25 žování agregátu, lze s pomocí měrné intenzity tvorby otěru, určené vztahem Xp = Xr/ sestrojit závislost měrného nárůstu koncentrace otěru v olejové lázni (při jednotkovém hnacím momentu) na čase zkoušky - viz obr. 4.
K tomu byla provedena regresní analýza této závislosti. S přihlédnutím na přijatý fyzikální model lineární závislosti opotřebovaného objemu na čase zkoušky, vyhoví nejlépe přímková regresní funkce. V této fázi je zachována interaktivita programu, protože obsluha se musí na základě grafického znázornění závislosti rozhodnout, zda podrobí analýze celý rozsah experimentu, či bude sledovat a vyěíslovat pouze fázi záběhu nebo jiné časové údobí. Je také možno cenzurovat chybně odebrané či chybně analyzované vzorky, způsobující často klamnou informací o zápor35 něm opotřebení (viz obr. 2). Regresní analýza se^provádí pomocí standardních vzorců. Získaná hodnota regresní měmé intenzity tvorby otěru X představuje odhad střední hodnoty produkce otěrových Částic v převodovce při hnacím momentu 1 Nm za časový úsek Δτ. Je-li zkušební metodika založena na cyklickém opakování zátěžového bloku, lze potom hodnoty intenzity tvorby částic za časový úsek přepočítat na intenzitu tvorby částic za jednotku převedené práce.
Postup při zajišťování míry opotřebení převodovky probíhal na zařízení podle obr. 5. Tvořil je stejnosměrný hnací motor 1 s cizím buzením: Pmax = 38,5 kW, nmax = 1500 ot/min tyristorovým regulátorem 2 otáček, nepřímý převod 3 klínovými řemeny o převodovém poměru i - 0,4, zkušební stojan 4 s hnacím vřetenem, na který se upíná zkoušená převodovka 5, přičemž výstup zkoušené převodovky 5 je spojen kloubovou hřídelí s redukční převodovkou 6, která zvyšuje vstupní otáčky do dynamometru ]_ s ohledem na jeho charakteristiku a disponuje převodovými stupni: it = 0,37, i2 ~ 0,61,13 = 1, a konečný vířivý dynamometr 7 o parametrech: Pmax = 125 kW, nmax = 8000 ot/min propojený s regulátorem 8 momentu.
Vířivý dynamometr 7 a redukční převodovka 6 jsou chlazeny vodou, zkoušená převodovka 5 vzduchem pomocí axiálního ventilátoru 9. Teplota redukční i zkoušené převodovky 6, 5 je udržována ve zvoleném rozmezí. Na obou převodovkách 5, 6 jsou za tímto účelem umístěna na několi- 10CZ 303481 B6 ka místech teplotní čidla, spolu se snímači hladiny chvění pro bezdemontážní diagnostiku technického stavu agregátů.
Způsob podle vynálezu na základě výsledků zjištěných po ukončení experimentu a údajů získa5 ných vzorkováním modeluje odlučovací a rozvrstvovací procesy olejové lázně převodovky.
Umožňuje zpětně určit nejen průměrnou koncentraci otěrových částic v libovolném okamžiku zkoušky, ale i odhadnout intenzitu produkce otěrových částic v daném úseku mezi vzorkovacími odběry, vyjádřenou v mg na 1 GJ převedené práce. Tato hodnota je potom objektivním měřítkem otěruvzdomosti testovaného převodového oleje nebo převodové skříně ve vybrané fázi testu.
Claims (4)
1. Způsob stanovení míry opotřebení převodovky (5) podle obsahu otěrových částic v olejové náplni, vyznačující se tím, že se nejprve v reálném provozu vozidla zjistí zátěžové spektrum pro každý převodový stupeň, tj. počet otáček vykonaných tímto stupněm v po sobě
20 následujících intervalech hodnoty kroutícího momentu přenášeného tímto stupněm vztaženo na jednotku dráhy vozidla, čímž se získá poměrná hodnota práce přenesené stupněm v jednotlivých intervalech - momentových hladinách, načež se sestaví pro jednotlivé stupně náhradní zátěžové bloky tak, že se vybere jedna kladná a jedna záporná momentová hladina s největší převedenou prací a do ní se připočte práce všech ostatních momentových hladin se stejným znaménkem, při25 čemž sestava takto vytvořených dílčích zátěžových bloků pro jednotlivé stupně vytvoří souhrnný zátěžový blok pro převodovku, určí se délka souhrnného zátěžového bloku odpovídající zvolené dráze vozidla nebo přenesené práci, na zkušebním stendu s reverzibilním pohonem a brzdným zařízením se testovaná převodovka provozuje tak, že se při konstantních otáčkách v hodnotě 0,5 až 0,8 maximálních otáček podrobuje zatížení po sobě následujícími zátěžovými bloky, přičemž
30 mezi jednotlivými zátěžovými bloky, vždy ihned po jejich ukončení, se z převodovky odeberou vzorky olejové náplně, které se podrobí analýze ke zjištění koncentrace částic železa v jednotlivých fázích zkoušky.
2. Způsob podle nároku 1, vyznačující se tím, že se zjištěná koncentrace částic
35 koriguje započtením poměrné části částic odebraných v předchozích vzorcích a částic usazených na stěnách převodovky (5) a rovněž se zohlední rozvrstvení částic v olejové náplni.
3. Způsob podle nároku 2, vyznačující se tím, že se ke zjištění středního nárůstu koncentrace částic během zkoušky podrobí hodnoty koncentrace regresní analýze.
4. Zařízení k provádění způsobu podle nároku 1, vyznačující se tím, že je tvořeno stejnosměrným hnacím motorem (1) s tyristorovým regulátorem (2) otáček, redukční převodovkou (6), zkušebním stojanem (4) s hnacím vřetenem, na který se upíná zkoušená převodovka (5) a dále vířivým dynamometrem (7) propojeným s regulátorem (8) momentu.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CZ20110665A CZ303481B6 (cs) | 2011-10-18 | 2011-10-18 | Zpusob stanovení míry opotrebení prevodovky a zarízení k provádení tohoto zpusobu |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CZ20110665A CZ303481B6 (cs) | 2011-10-18 | 2011-10-18 | Zpusob stanovení míry opotrebení prevodovky a zarízení k provádení tohoto zpusobu |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CZ2011665A3 CZ2011665A3 (cs) | 2012-10-10 |
| CZ303481B6 true CZ303481B6 (cs) | 2012-10-10 |
Family
ID=46964886
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CZ20110665A CZ303481B6 (cs) | 2011-10-18 | 2011-10-18 | Zpusob stanovení míry opotrebení prevodovky a zarízení k provádení tohoto zpusobu |
Country Status (1)
| Country | Link |
|---|---|
| CZ (1) | CZ303481B6 (cs) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10588268B2 (en) | 2016-06-28 | 2020-03-17 | Deere & Company | Baler with wear reducing controller |
| CN113062901A (zh) * | 2021-06-03 | 2021-07-02 | 盛瑞传动股份有限公司 | 自动变速器测试台液压控制系统及测试方法 |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109307594B (zh) * | 2018-09-28 | 2024-04-26 | 浙江海天机械有限公司 | 一种摩托齿轮防磨损检测装置 |
| CN111198001B (zh) * | 2018-11-16 | 2022-09-30 | 中国北方车辆研究所 | 一种液力减速器应力应变的测试方法及系统 |
| CN115952693B (zh) * | 2023-03-13 | 2023-08-04 | 中国重汽集团济南动力有限公司 | 基于大数据的变速箱载荷谱转化方法、装置、设备、介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2116896U (zh) * | 1992-02-29 | 1992-09-23 | 济南交通高等专科学校 | 汽车变速器磨合试验台 |
| JP2002286585A (ja) * | 2001-03-28 | 2002-10-03 | Ishikawajima Harima Heavy Ind Co Ltd | 歯車の摩耗試験装置 |
| CN201653713U (zh) * | 2010-03-28 | 2010-11-24 | 荆州恒盛汽车系统有限公司 | 一种动力转向器磨损试验的液压装置 |
-
2011
- 2011-10-18 CZ CZ20110665A patent/CZ303481B6/cs not_active IP Right Cessation
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2116896U (zh) * | 1992-02-29 | 1992-09-23 | 济南交通高等专科学校 | 汽车变速器磨合试验台 |
| JP2002286585A (ja) * | 2001-03-28 | 2002-10-03 | Ishikawajima Harima Heavy Ind Co Ltd | 歯車の摩耗試験装置 |
| CN201653713U (zh) * | 2010-03-28 | 2010-11-24 | 荆州恒盛汽车系统有限公司 | 一种动力转向器磨损试验的液压装置 |
Non-Patent Citations (1)
| Title |
|---|
| (Sledovani mno×stvi oterovych kovu v mazacich olejich, prednaÜka na seminari Tribotechnika v provozu a ·dr×be, Ing.M.Holec a kol. [http://www.crc.cz/data/publications/sledovani_oterovych_kovu_v_mazacich_olejich.pdf]), 2.11.2010 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10588268B2 (en) | 2016-06-28 | 2020-03-17 | Deere & Company | Baler with wear reducing controller |
| CN113062901A (zh) * | 2021-06-03 | 2021-07-02 | 盛瑞传动股份有限公司 | 自动变速器测试台液压控制系统及测试方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CZ2011665A3 (cs) | 2012-10-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CZ2011665A3 (cs) | Zpusob stanovení míry opotrebení prevodovky a zarízení k provádení tohoto zpusobu | |
| MX2014014463A (es) | Medios y metodos que aplican sflt-1/plgf o relacion de endoglina/plgf para descartar el inicio de preeclampsia en un determinado periodo de tiempo. | |
| WolAk et al. | The kinetics of changes in kinematic viscosity of engine oils under similar operating conditions | |
| Sepehrnia et al. | Rheological behavior of MWCNT-SnO2/SAE50 hybrid nanolubricant: experimental evaluation and viscosity prediction using optimized machine learning model | |
| Sharma et al. | Experimental study of tribological behavior of casted aluminium-bronze | |
| CN113916843A (zh) | 基于灰色模型的机械设备润滑油性能预测方法 | |
| CH708013A2 (de) | Vorrichtung zur Überwachung eines Schmieröls in einem Vorratsbehälter einer Turbomaschine. | |
| CN102519721B (zh) | 变速器齿轮接触斑点实验装置 | |
| Nakhaei et al. | Prediction of XRF analyzers error for elements on-line assaying using Kalman Filter | |
| Wilson et al. | Development of a modular in-situ oil analysis prognostic system | |
| JP7363692B2 (ja) | 劣化推定装置、劣化推定方法、及び劣化推定プログラム | |
| Sitzmann et al. | The integration of predictive and preventive maintenance strategies based on oil and vibration analyses techniques for industrial gearboxes | |
| Hein et al. | Test method for time-scaled fatigue tests of gear transmission systems | |
| Seo | A study on improving the prediction accuracy of cold forging die life based on quantitative evaluation of phosphate film damage | |
| Mohammed et al. | Control system design of a Pin-on-Disc wear machine using PLC | |
| Xing et al. | Fatigue Wear Analysis and Life Prediction of Surgical Instrument Driving Wires | |
| JP4836773B2 (ja) | 潤滑状態評価装置、潤滑状態評価方法、プログラム及び記録媒体 | |
| Furch et al. | The Use of Portable Tribodiagnostic Equipment to Determine the Technical Condition of an Internal Combustion Engine within the Life Cycle | |
| Baron et al. | Research and correlation of diagnostic methods for assessment of the state of oil filling in cycloid gearbox | |
| Karajagikar et al. | Investigate the condition of single speed worm gearbox using ferrography and vibration analysis,' | |
| Iqbal et al. | Brake Pad Life Monitoring System Using Machine Learning | |
| Stodola | The results of ferrography tests and their evaluation | |
| Ori et al. | Vibration Diagnostic Methods of Automatic Transmission Service Requirement Prediction. Infrastructures 2023, 8, 31 | |
| Lee et al. | Accounting for driver variability in vehicle powertrain durability testing | |
| WO2026105369A1 (ja) | 寿命推定システム、寿命推定方法、及び寿命推定プログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | Patent lapsed due to non-payment of fee |
Effective date: 20181018 |