MISKOLCI EGYETEM
GÉPÉSZMÉRNÖKI ÉS INFORMATIKAI KAR
PNEUMATIKUS ELZÁRÓ-SZERKEZETTEL ÜRÍTETT ÖMLESZTETT ANYAG TÁROLÓ SILÓK ADAGOLÁSI MODELLJE Ph.D. értekezés KÉSZÍTETTE:
Telek Péter okleveles gépészmérnök
SÁLYI ISTVÁN GÉPÉSZETI TUDOMÁNYOK DOKTORI ISKOLA GÉPEK ÉS SZERKEZETEK TERVEZÉSE TÉMATERÜLET ANYAGMOZGATÓ GÉPEK TERVEZÉSE TÉMACSOPORT DOKTORI ISKOLA VEZETŐ: Dr. Tisza Miklós egyetemi tanár, DSc TÉMACSOPORT VEZETŐ: Dr. Lévai Imre Professzor Emeritus TÉMAVEZETŐ: Dr. Szaladnya Sándor Professor Emeritus
Miskolc, 2012.
TARTALOMJEGYZÉK A TÉMAVEZETŐ AJÁNLÁSA ……...………………………………………………………4 BEVEZETÉS ……………………………...…………………………………………………..5 AZ ÉRTEKEZÉSBEN ALKALMAZOTT JELÖLÉSEK ………...…………………………..7 1. ÖMLESZTETT ANYAGOK ADAGOLÁSA ……………...……………………………..12 1.1. Közvetlen adagolás ...………..…….…………………………..……………………13 1.2. Közvetett adagolás ...………………..………………………………………….…..14 1.3. Kombinált adagolás ...…..……………..………………………………………….…..15 1.4. Az adagolt mennyiség ……..………..……………………………………………….16 1.4.1. Az adagolt mennyiség közelítése ……..……………………………………..18 1.4.2. Az adagolás komplex összefüggés rendszere ...………………………..………18 2. SILÓK ELZÁRÓ-SZERKEZETEI ...…………...…….…………………………............20 2.1. Tolózárak ……………..…………………….………………………………………..20 2.2. Szegmens-zárak ……….…...………..………..……………………………………...23 2.3. Elzáró-szerkezetek méretezése ...…………………..…………………………………25 3. ÖMLESZTETT ANYAGOK ADAGOLÁSI JELLEMZŐI ……......................................27 3.1. Statikus állapot-jellemzők ……..…………………………………………………......28 3.2. Állapot-jellemzők kiömlés közben ………...…………………………………………32 3.3. Kiömlő anyagmennyiség ...…………..……...……………………………….……….34 4. PNEUMATIKUS MUNKAHENGEREK MŰKÖDÉSI JELLEMZŐI ……..………........39 4.1. Mechanikai modell ..........…………………………………………………….............40 4.2. Hőtani modell ...……………...……………………………………………………….41 4.3. Áramlástani modell ...…………………………..…………………………….............43 4.4. Pneumatikus munkahengerek működését leíró egyenletrendszer .....………………...45 5. PNEUMATIKUS MUNKAHENGEREK ÁLTALÁNOS MŰKÖDÉSI MODELLJE .......48 5.1. A működés fázisai ...………………………………………………………………….48 5.2. Működési modell ...……………………………………...……………………………49 5.3. A differenciálegyenlet-rendszer megoldása ...………………………………………..51 5.4. A modell eredményeinek ellenőrzése ....……………………………………...............54 5.4.1. A vizsgált munkahenger ...……………………………………………………..54 5.4.2. A mérés elrendezési vázlata és elemei ...………………………………………56 5.4.3. A mérés eredményei ...…………………………………………………………57 5.4.4. A mért és számított eredmények összehasonlítása ...…….………………...…..58 6. PNEUMATIKUS MŰKÖDTETÉSŰ ELZÁRÓ-SZERKEZETEK ERŐJÁTÉKA ............62 6.1. Tolózárak erőjátéka ...…………………………………………………………...........62 6.2. Kéttagú szegmens-zárak erőjátéka ...……...…………………………….....................64
-2-
7. ÖMLESZTETT ANYAGOK VERTIKÁLIS KIÖMLÉSE ...…..........................................67 7.1. Adagolt ömlesztett anyagok ………………………………...………………..............67 7.2. Adagolási környezet ……………………………………………...…………………..69 7.3. A kiömlő mennyiség számítására alkalmas összefüggés kiválasztása …..…………...71 8. AZ ADAGOLÁSI FOLYAMAT ÁLTALÁNOS MODELLJE ...…..………….................75 8.1. A modell elemei ...……………………………………………..…………..................75 8.2. A modell eredményeinek ellenőrzése ...………………………...…………………….78 8.2.1. Mérési környezet …………………………………………..…………………..78 8.2.2. A mérés eredményei ...……………………………………...………………….79 8.2.3. A mért és számított eredmények összehasonlítása ...………...………………...81 .
9. AZ ADAGOLÁSI MODELL ALKALMAZÁSA ...……………………..………………..83 9.1. Az adagolási folyamat paramétereinek meghatározása ….…………………………...83 9.2. Elzáró-szerkezetek kiválasztása ...……………………………….……………...........84 10. ÚJ TUDOMÁNYOS EREDMÉNYEK ...………………………….…………………….87 ÖSSZEFOGLALÁS ...…………………………………………………..……………………89 SUMMARY ...………………………………………………………………………………..90 KÖSZÖNETNYILVÁNÍTÁS ...……………………………………………………………..91 FELHASZNÁLT IRODALMAK JEGYZÉKE ...………………………...………….………92 AZ ÉRTEKEZÉSHEZ KAPCSOLÓDÓ PUBLIKÁCIÓK ...………………………………..94 MELLÉKLETEK ...…………………………………………………………………………..96 M1. Munkahenger paraméterek mérési eredményei ...……………………..…..………....97 M2. Mért és számított munkahenger paraméterek összehasonlítása ...……...….....……..100 M3. Ömlesztett anyagok kiömlésével kapcsolatos mérések ...……………..…………....106 M4. Az adagolási folyamat mérési és számítási eredményei ...………….…..…………..109 M5. A mérésekhez használt berendezések ...…………………………………..………...112
-3-
A TÉMAVEZETŐ AJÁNLÁSA Telek Péter Pneumatikus elzáró-szerkezettel ürített ömlesztett anyag tároló silók adagolási modellje című Ph.D. értekezéséhez
Első ránézésre az ömlesztett anyagok adagolása egyszerű feladat. Ez igaz is akkor, ha az adagolás alárendelt szerepet játszik és az anyagmennyiség ellenőrzésére külön berendezéseket alkalmaznak. Ha viszont egyszerű és pontos adagolást szeretnénk megvalósítani, akkor figyelembe kell venni a szerkezet, az anyag és a működtető elemek komplex rendszerének minden elemét, ami bizonyos esetekben rendkívül bonyolult feladatot eredményez. Különösen igaz ez a pneumatikus munkahengerrel működtetett adagolási folyamatok esetén, ahol a feladat komplexitását növeli a munkahengert működtető közeg összenyomhatósága is. Ezt a szakirodalomban eddig nem tárgyalt témát dolgozta fel doktori értekezésében Telek Péter, amellyel nem könnyű feladatot vállalt magára. Kutatásainak célja egy olyan általános modell kidolgozása volt, amellyel az adagolási folyamat leírható. Értekezésében bemutatja a modell megalkotásához szükséges tématerületek szakirodalmi hátterét, és az adagolási folyamatot leszűkíti bizonyos szerkezeti megoldásokra és ömlesztett anyagokra. Az adagolási modell fontos eleme a pneumatikus munkahenger működése, melynél Telek Péter kidolgozott egy modellt az adagoláshoz szükséges paraméterek (pl. dugattyú elmozdulás, rúderő, stb.) meghatározására. Az adagolási modell vizsgálatára számítógépes programot készített és méréseket végzett, amellyel igazolta az eredmények helytállóságát. A gyakorlat szempontjából az értekezésben foglalt eredményeknek fontos szerepe lehet azokban az ipari folyamatokban, ahol a termelékenység, az egyszerűség és a pontos adagolás magas követelményeket támasztanak az adagolóeszközökkel szemben. Telek Péter doktori tevékenysége során a doktori témájához kapcsolódóan 21 darab publikációt jelentetett meg, melyeket jó szakmai ismeretekkel, önálló, értékes gondolatokkal alkotott. Emellett eredményeit számos hazai és nemzetközi konferencián bemutatta. Szakmai tevékenysége, hozzáállása kiemelkedő, az értekezés fejezetei megfelelő színvonalúak. Közös munkáinkban tevékenyen, önálló elképzelésekkel vett részt. Az értekezésben megjelenő új eredmények Telek Péter önálló eredményei. A tézisekben megfogalmazott új eredmények a doktori (Ph.D.) fokozat elnyeréséhez szükséges kívánalmakat messzemenően kielégítik.
…………………………………. Dr. Szaladnya Sándor Professor Emeritus témavezető Miskolc, 2012. június 14.
-4-
BEVEZETÉS Napjainkban, a technika ugrásszerű fejlődésének következtében, az ipari folyamatok technikai színvonala nagyon fejlett, hatékonyságuk növeléséhez a meglévő berendezések, eszközök kapacitásainak maximális kihasználására van szükség. Ehhez sok esetben felül kell vizsgálni a berendezések működését, és olyan módszereket kell alkalmazni azok méretezésére, amelyeket eddig bonyolultságuk miatt nem használtak. Jó példa erre a számítógépes szimulációk, véges elemes módszerek, és más korszerű megoldások terjedése. Az ömlesztett anyagok adagolási folyamatára is elmondható ugyanez, ahol a pneumatikus elzáró-szerkezettel történő silóürítés méretezése hagyományosan tapasztalati, közelítő összefüggések alapján történik. Az így elérhető adagolási pontosság azokban a technológiai folyamatokban, ahol a minőséget alapvetően meghatározza az egyes anyagok aránya (pl. üveggyártás), nem megfelelő. A pontosság növeléséhez új, a folyamatot jobban közelítő méretezési eljárásra van szükség, ami csak az adagolási folyamatra felírt általános modell alapján készíthető el. Doktori kutatásaim célja: az adagolási folyamat általános leírására alkalmas modell kidolgozása pneumatikus elzáró-szerkezettel ürített ömlesztett-anyag tároló silók esetén. A modell alkalmazásával lehetővé válik: a kiömlő anyagmennyiség megfelelő pontosságú meghatározása az adagolási folyamat paramétereinek függvényében, az adagolási folyamat időszükségletének csökkentése, a követelményekhez jobban illeszkedő elzáró-szerkezetek tervezése, stb. A gyakorlatban az ömlesztett anyagok adagolása egy külön adagoló-szerkezet közbeiktatásával történik, ezért a kapcsolódó szakirodalom alig foglalkozik az elzárószerkezet ilyen irányú méretezésével. Pneumatikusan működtetett elzáró-szerkezetekre vonatkozóan pedig egyáltalán nem ismert ilyen szakirodalom. Kutatási célom eléréséhez a következő feladatokat kell megoldanom: 1. Az adagolási folyamat definiálása és a modellben alkalmazható, pneumatikusan működtethető elzáró-szerkezetek kiválasztása 2. Az adagolási folyamatot befolyásoló ömlesztett anyagjellemzők vizsgálata, a szükséges paraméterek meghatározása 3. Pneumatikus munkahengerek működését leíró általános modell kidolgozása 4. Az adagolási folyamatot leíró általános modell kidolgozása 5. A modellek vizsgálatára alkalmas számítógépes programok elkészítése 6. Mérési folyamatok összeállítása a modellek eredményeinek ellenőrzésére Doktori értekezésem felépítését ahhoz igazítottam, hogy a kutatási téma három különböző szakterületet kapcsol össze, ezért négy különálló fejezet (1-4.) tartalmazza az általános folyamatleírást és az egyes területekhez kapcsolódó szakirodalmak áttekintését és négy fejezetben (5-8.) kerülnek bemutatásra a kutatás során elvégzett feladatok. Az 1. fejezet tartalmazza az adagolási folyamat általános leírását, elemeit, jellemzését és az adagolt anyagmennyiség meghatározásának módját. Mivel az adagolás megvalósítására különböző megoldásokat lehet alkalmazni, ezért részletesen bemutatom az egyes változatok jellemzőit. A fejezet végén összefoglalom a folyamatot befolyásoló paramétereket és a közöttük fennálló kapcsolatokat.
-5-
Az elzáró-szerkezet szerkezeti megoldása jelentős hatást gyakorol az adagolási folyamat idejére és intenzitására. A 2. fejezetben ismertetem a gyakorlatban előforduló siló-elzáró szerkezetek csoportosítását, főbb típusait és az egyes változatok jellemzőit. Az adagolási folyamat jellemzői alapján kiválasztom az általam vizsgált esetben alkalmazható szerkezeti megoldásokat. Az ömlesztett anyagok viselkedése kétféle módon befolyásolja az adagolási folyamatot: o az anyag nyomása az elzáró-szerkezeten keresztül erőt fejt ki a mozgató pneumatikus munkahenger dugattyújára, o az anyag áramlási tulajdonságai befolyásolják a kiömlő mennyiséget. A fenti két tényező meghatározásához szükség van az ömlesztett anyagok tulajdonságait, viselkedését leíró összefüggések felírására, amely a 3. fejezetben található. Mivel az ömlesztett anyagok viselkedése meglehetősen bonyolult, ezért külön bemutatom a statikus állapotra és a kiömlés közben érvényes összefüggéseket. Az elzáró-szerkezet működtetésére egy pneumatikus munkahengert alkalmazok, ezért a 4. fejezetben bemutatom a munkahenger működési paramétereinek meghatározására alkalmas összefüggéseket. A működés leírásánál a levegő tulajdonságaiból, illetve a levegőre érvényes hőtani, illetve áramlási összefüggésekből, illetve a mozgásegyenletből kell kiindulni. A fejezet végén definiálom azt a differenciálegyenlet-rendszert, amivel a működés általánosan leírható. A pneumatikus munkahenger működését leíró differenciálegyenlet-rendszer analitikusan nem oldható meg, ezért az 5. fejezetben bemutatok egy új, általam kidolgozott általános modellt, amelynek alkalmazásával az egyenlet-rendszer numerikusan megoldható. A modell alapján elkészített számítógépes program eredményeit mérésekkel ellenőrzöm. A munkahenger dugattyújára ható erő nagyságát és a kiömlési keresztmetszet méretét az elzáró-szerkezet geometriája határozza meg. A 6. fejezetben mutatom be a működés során fellépő erőket és a kiömlési keresztmetszet változását leíró matematikai összefüggéseket, a kiválasztott elzáró-szerkezet megoldások esetén. Az egyes ömlesztett anyagok, különböző belső tulajdonságaik miatt, teljesen eltérő folyási viselkedést mutathatnak, ezért az adagolási modell nem vonatkozhat minden anyagra. A 7. fejezetben kiválasztok egy anyagtípust, amire a modell alkalmazható és a 3. fejezetben leírt kiömlési összefüggések közül, mérések alapján kiválasztom az adott anyagtípusra használható változatokat. A 8. fejezetben ismertetem a teljes adagolási folyamat általam kidolgozott általános modelljét. A modell alapján elkészített számítógépes program eredményeit mérésekkel ellenőrzöm. A 9. fejezetben bemutatom az általános modell alkalmazhatóságát az adagolt anyagmennyiség és más paraméterek meghatározására, illetve az adagolás vezérlésére. A fejezet végén példát mutatok be az elzáró-szerkezetek adott feladathoz történő kiválasztására, összehasonlító vizsgálatára.
-6-
AZ ÉRTEKEZÉSBEN ALKALMAZOTT JELÖLÉSEK Latin betűvel jelölt változók: A+ AAcf+ AcfAe Af Ak Akh Amax An As As+ AsAsf+ AsfAt Az b B cd cp C Cb dc dc+, dc-, dcf+, dcf-, dg dh dhf dp Dekv Dg Dh Dk dAM dhm dp/dr dQ dU dW
- a dugattyú felülete a plusz kamrában [m2] - a dugattyú felülete a mínusz kamrában [m2] - az átömlési keresztmetszet fékezéskor a plusz kamra esetén [m2] - az átömlési keresztmetszet fékezéskor a mínusz kamra esetén [m2] - az elzáró-szerkezet anyaggal érintkező felülete [m2] - a beömlési keresztmetszet [m2] - a kiömlési keresztmetszet [m2] - a kiömlési keresztmetszet horizontális vetülete [m2] - a kifolyónyílás maximális keresztmetszete [m2] - a kifolyónyílás nyitott keresztmetszete [m2] - a silószár keresztmetszete [m2] - a fojtószelep átömlési keresztmetszete nyitott irányban (plusz kamra) [m2] - a fojtószelep átömlési keresztmetszete nyitott irányban (mínusz kamra) [m2] - a fojtószelep átömlési keresztmetszete fojtott irányban (plusz kamra) [m2] - a fojtószelep átömlési keresztmetszete fojtott irányban (mínusz kamra) [m2] - a kifolyógarat keresztmetszete a kúp alakú térfogatelemen [m2] - a kifolyónyílás zárt keresztmetszetének mérete [m2] - az elzáró-szerkezet működési sebességétől függő tényező - a silószár hosszúsága (téglalap alak) [m] - a külső rugalmassági tényező [N/m] - izobár fajhő [J/(kg∙K)] - az anyag kohéziójának mértéke [Pa] - kifolyási tényező - a csapátmérő [m] - a plusz kamra levegőcsatlakozásának átmérője [m] - a mínusz kamra levegőcsatlakozásának átmérője [m] - kiáramlási átmérő fékezéskor a plusz kamránál [m] - kiáramlási átmérő fékezéskor a mínusz kamránál [m] - a görgőcsapok átmérője [m] - a dugattyúrúd átmérője [m] - a fékeződugattyú átmérője [m] - az átlagos szemcseátmérő [m] - az ekvivalens kifolyónyílás átmérő [m] - a görgők átmérője [m] - a dugattyú átmérője [m] - a kifolyónyílás átmérője (kör alakú) [m] - a kúp alakú térfogatelem palástjának felülete [m2] - a kiáramló térfogatelem magassága [m] - a kifolyógaratban, a szemcsék között kialakuló nyomásváltozás - a rendszerbe bevitt hőmennyiség [J] - a rendszer belső energia változása [J] - a rendszer által végzett munka [J]
-7-
dzs dzt ez f1-f6 ff ffakt Fa Fd Fdug Fe Fh Fm Fmax FR Fs Fszeg Fv g Ge Ges h h0 ht Hs kd kp kt k Kh lg L Lc Lk Lm Ln m md mg ms mt mz ma
- a hengeres térfogatelem magassága [m] - a kúp alakú térfogatelem magassága [m] - a tolózár elmozdulása [m] - függvénykapcsolatok - a folyási faktor - a folyási faktor aktuális értéke - az ömlesztett anyag nyomásából a zárólapra ható erő [N] - a dugattyúra ható összes külső terhelőerő [N] - a munkahenger mozgó részeinek súlyából származó erő [N] - az anyagból a zárólapra ható eredő terhelés [N] - az anyag horizontális nyomásából az elzáró elemre ható erő [N] - a munkahenger mozgó részeinek és a zárólap súlyából származó erő [N] - a maximális működtető erő [N] - a dugattyúra ható rúderő [N] - a dugattyúra ható összes súrlódó erő [N] - az elzáró szegmens súlyából származó erő [N] - az anyag vertikális nyomásából az elzáró elemre ható erő [N] - a gravitációs gyorsulás [m/s2] - az elzáró-szerkezet mozgó elemeinek önsúlya [N] - az ellensúly nagysága [N] - az áramló anyagmagasság [m] - a kifolyógarat kúpmagassága (vagy gúla magassága) [m] - a kifolyógarat magassága [m] - a silószár magassága [m] - a külső csillapítási tényező [Ns/m] - a szemcsék alaktényezője - biztonsági tényező (a szemcsenagyság függvénye) - a silógarat geometriai tényezője - a silógarat kúpszögétől, a belső és a falsúrlódási szögtől függő tényező - a gördülő-ellenállás karja [m] - a maximális lökethossz [m] - a szegmens és a munkahenger csuklópontja közötti távolság [m] - a kifolyónyílás hossza (téglalap alak) [m] - a munkahenger hossza mínusz véghelyzetben [m] - a kiömlőnyílás aktuális hossza [m] - a munkahenger által mozgatott összes tömeg [kg] - a munkahenger mozgó részeinek tömege [kg] - geometriai jellemző - egy szegmens tömege [kg] - a garatgeometriára vonatkozó tényező - a zárólap tömege [kg] - az időegység alatt kiáramló anyagmennyiség [kg/s]
mb m ba m bh
- a kamrákba időegység alatt beáramló levegő mennyisége [kg/s] - a kamrákba hangsebesség alatt beáramló levegő mennyisége [kg/s] - a kamrákba hangsebességgel beáramló levegő mennyisége [kg/s]
-8-
mk
m ka m kh
n p p+ pp0 pA, pB pb pk pt pvkrit qa Qa Qm Qn Qs Qt Qz R Re Res Rz sd sf sz tn tk tkv tsk tt tz T T+ TT0 Tb Tk Tt x xst Uk Us v0
- a kamrákból időegység alatt kiáramló levegő mennyisége [kg/s] - a kamrákból hangsebesség alatt kiáramló levegő mennyisége [kg/s] - a kamrákból hangsebességgel kiáramló levegő mennyisége [kg/s] - a szerkezet és az anyag tulajdonságaitól függő tényező - a nyomás az adott kamrában [Pa] - a nyomás a plusz kamrában [Pa] - a nyomás a mínusz kamrában [Pa] - az atmoszférikus nyomás [Pa] - a nyomás az egyes terekben [Pa] - a beáramló levegő nyomása [Pa] - a kiáramló levegő nyomása [Pa] - a működtető nyomás [Pa] - a kritikus nyomásviszony - a kiömlő anyagáram [m3/s] - az adagolás alatt kiáramló anyagmennyiség [m3] - a zárást indító anyagmennyiség [m3] - a nyitási fázis alatt kiáramló anyagmennyiség [m3] - a technológia által igényelt anyagmennyiség [m3] - a nyitvatartási fázis alatt kiáramló anyagmennyiség [m3] - a zárási fázis alatt kiáramló anyagmennyiség [m3] - az egyetemes gázállandó [J/(kg∙K)] - a működtető erő karja [m] - az ellensúly karja [m] - az elzáró szegmens sugara [m] - a pneumatikus munkahenger dugattyújának elmozdulása [m] - a fékezési úthossz [m] - a dugattyú mozgási útja a munkahenger működésének 2. fázisában maximális kiömlés (zárás), illetve zárt szerkezet (nyitás) esetén [m] - a nyitás ideje [s] - a kiömlés ideje [s] - a közvetlen adagolás vezérlési ideje [s] - a sebesség kiegyenlítődési idő [s] - a nyitva tartás ideje [s] - a zárás ideje [s] - a levegő hőmérséklete az adott kamrában [K] - a levegő hőmérséklete a plusz kamrában [K] - a levegő hőmérséklete a mínusz kamrában [K] - külső levegő hőmérséklete [K] - a beáramló levegő hőmérséklete [K] - a kiáramló levegő hőmérséklete [K] - a működtető levegő hőmérséklete [K] - a szemcsék alakjától függő kitevő - az ekvivalens szemcseátmérő [mm] - a kiáramló térfogatelem kerülete [m] - a hengeres térfogatelem kerülete [m] - a kiömlő anyag sebessége a kifolyónyílásban [m/s]
-9-
va vb vba vbh vB vk vka vkh V V+ V+max VV-max V+0 V-0 Vf+0 Vf-0 W Wk
- az áramló anyag átlagsebessége [m/s] - a beáramló levegő sebessége [m/s] - a hangsebesség alatt beáramló levegő sebessége [m/s] - a hangsebességgel beáramló levegő sebessége [m/s] - a B térbe beáramló levegő sebessége [m/s] - a kiáramló levegő sebessége [m/s] - a hangsebesség alatt kiáramló levegő sebessége [m/s] - a hangsebességgel kiáramló levegő sebessége [m/s] - az adott kamra térfogata [m3] - a plusz kamra térfogata [m3] - a plusz kamra maximális térfogata [m3] - a mínusz kamra térfogata [m3] - a mínusz kamra maximális térfogata [m3] - holttér a plusz kamrában [m3] - holttér a mínusz kamrában [m3] - holttér a plusz kamrában fékezéskor [m3] - holttér a mínusz kamrában fékezéskor [m3] - a silószár szélessége (téglalap alak) [m] - a kifolyónyílás szélessége (téglalap alak) [m]
Görög betűvel jelölt változók: e s 1 c d a es mh s 1 c cf+ cff g
- a tolózár dőlésszöge [o] - a szegmens-zár nyitási szöge [o] - az elzáró-szegmens átfogási szöge [o] - az elzáró-szerkezet által meghatározott szög [o] - a szegmens mozgása és a munkahenger helyzete által meghatározott szög [o] - a munkahenger helyzete által meghatározott szög [o] - a szemcsék alaktényezője - az ömlesztett anyag jellege - az elzáró-szerkezet kialakítási jellege - a pneumatikus munkahenger kialakítási, működtetési jellege - a siló kialakítási jellege - az ömlesztett anyag porozitása - a levegő dinamikus viszkozitása [Pas] - a kifolyógarat dőlésszöge [o] - az adiabatikus kitevő - nyomásviszony az ömlesztett anyagban - a zárólap és támasztófelülete közötti súrlódási tényező - a görgőcsapok súrlódási tényezője - a fojtási tényező fékezéskor a plusz kamra esetén - a fojtási tényező fékezéskor a mínusz kamra esetén - a fojtási tényező - a gördülő-ellenállási tényező
- 10 -
h lev s+ ssf+ sfx 0 + a A, B lev p az m m v v0 hs ht vs vt wt wm ws wt e x Ψb Ψba Ψbh Ψk Ψka Ψkh
- súrlódási tényező a hengerben - a szemcsék közötti levegő viszkozitása [m2/s] - a fojtószelep fojtási tényezője nyitott irányban (plusz kamra) - a fojtószelep fojtási tényezője nyitott irányban (mínusz kamra) - a fojtószelep fojtási tényezője fojtott irányban (plusz kamra) - a fojtószelep fojtási tényezője fojtott irányban (mínusz kamra) - falsúrlódási tényező az anyag és a zárólap között - az ömlesztett anyag tulajdonságaitól függő tényező - a levegő sűrűsége atmoszférikus nyomáson [kg/m3] - a levegő sűrűsége a plusz kamrában [kg/m3] - a levegő sűrűsége a mínusz kamrában [kg/m3] - az ömlesztett anyag sűrűsége [kg/m3] - a levegő sűrűsége az egyes terekben [kg/m3] - a szemcsék közötti levegő sűrűsége [kg/m3] - a szemcsesűrűség [kg/m3] - az ömlesztett anyagban ébredő feszültség [Pa] - az álló anyagrészen fellépő horizontális feszültség [Pa] - vertikális feszültség a kiáramló anyag adott térfogatelemén [Pa] - az anyagból az elzáró-szerkezetre ható függőleges irányú nyomás [Pa] - a vertikális feszültség a kifolyógarat és a silószár találkozásánál [Pa] - horizontális feszültség a silószár adott térfogatelemén [Pa] - horizontális feszültség a kifolyógarat adott térfogatelemén [Pa] - vertikális feszültség a silószár adott térfogatelemén [Pa] - vertikális feszültség a kifolyógarat adott térfogatelemén [Pa] - a kifolyógarat oldalfalára merőleges feszültség [Pa] - az álló és a mozgó anyag között ébredő nyírófeszültség [Pa] - nyírófeszültség a siló hengeres oldalfalán [Pa] - nyírófeszültség a kifolyógarat oldalfalán [Pa] - az anyag belső súrlódási szöge [o] - a falsúrlódási szög [o] - a beáramlásra vonatkozó átömlési függvény - a hangsebesség alatti beáramlásra vonatkozó átömlési függvény - a hangsebességgel történő beáramlásra vonatkozó átömlési függvény - a kiáramlásra vonatkozó átömlési függvény - a hangsebesség alatti kiáramlásra vonatkozó átömlési függvény - a hangsebességgel történő kiáramlásra vonatkozó átömlési függvény
- 11 -
1. ÖMLESZTETT ANYAGOK ADAGOLÁSA Ömlesztett anyagok adagolása alatt azt a folyamatot értjük, amikor egy szállító vagy tárolóeszközből, folytonos anyagáram, illetve nagyobb ömlesztett mennyiség (ömlesztett formában tárolva) formájában rendelkezésre álló anyagból, egy előírt recept, vagy más igény szerint kimért mennyiséget egy adott technológiai folyamatba, vagy egy másik szállítóeszközbe juttatunk. Az adagolási folyamat az alkalmazott szállító és tároló eszközöktől, illetve az adagolt ömlesztett anyagtól függ [1.]. A kutatás témájához igazodva, az értekezésben csak a tároló silókból pneumatikus működtetésű, mechanikus elzáró-szerkezettel végzett adagolással foglakozom (1/1. ábra). Ömlesztett anyag feladás
Siló
Ömlesztett anyag
Elzáró-szerkezet
Mérőadagoló berendezés
Pneumatikus munkahenger Ömlesztett anyag adagolás
Beszállítás a technológiai folyamatba
1/1. ábra Az adagolás általános folyamata silóürítés során Az adagolási folyamat leírásához szükség van az egyes folyamatelemek definiálására és a közöttük fennálló kapcsolatok megadására. A folyamat elemei az 1/1. ábra alapján a következők: az adagolt ömlesztett anyag, a siló, az elzáró-szerkezet, a mérő-adagoló szerkezet, a kapcsolódó technológiai, illetve szállító berendezés. Az ömlesztett anyag jelenti az adagolás tárgyát, ennek viselkedése határozza meg a folyamat tényleges formáját [2.]. Mivel az ömlesztett anyagoknak az adagolás szempontjából fontos mechanikai és folyási tulajdonságai nagymértékben függnek az anyag jellemzőitől (szemcse alak, méret, eloszlás, sűrűség, kohézió, belső súrlódás, stb.), az anyag viselkedésének leírása bonyolult feladat. A tárolósiló az egyetlen olyan tárolási módszer (a szabadtéri garmadás, a nyitott, vagy a zárt hombáros tárolással szemben) az ömlesztett anyagok számára, amely könnyen beilleszthető a - 12 -
korszerű technológiai folyamatokba. A siló egy olyan zárt, lábakon álló tartály (henger, vagy hasáb alakú), amelybe az ömlesztett anyagot egy külső rakodó berendezés segítségével (pneumatikus rakodás, rakodógépes rakodás, stb.) juttatják be (feladás), és az ürítés általában külső eszköz nélkül (az anyag tulajdonságainak függvényében), a kifolyónyílásra épített elzáró szerkezettel valósítható meg. A siló tulajdonságai (alak, méret, anyag, a kifolyónyílás alakja és mérete, stb.) befolyással vannak az anyag kiömlési jellemzőire, de ezek általában egyszerű összefüggésekkel figyelembe vehetők [3.]. A silóra épített elzáró-szerkezet nyitja, illetve zárja a kifolyónyílást, amivel szabaddá teszi, elzárja, vagy szabályozza az ömlesztett-anyag kiáramlását [1.]. Az elzáró-szerkezet paraméterei (geometriai méretek, működési sebesség, kialakítás, stb.) közvetlenül befolyásolják az adagolási folyamat jellemzőit (kiömlő mennyiség, időtartam, stb.). Egy adott silóhoz nagyon sokféle elzáró-szerkezetet alkalmazhatunk, amelyek működtetésére számos lehetőség áll rendelkezésre (pl. elektromos, hidraulikus, pneumatikus, stb.). A mérő-adagoló szerkezet határozza meg az adagolt mennyiséget [4.]. Hatása az anyag kifolyási jellemzőire általában közvetett (bizonyos esetekben befolyásolhatja azokat, pl. rezgések, geometriai kialakítás következtében, de ez általában kiküszöbölhető, elhanyagolható), elsősorban (az alkalmazott megoldás függvényében) az adagolási folyamat elindítását, illetve leállítását vezérli. A technológiai folyamatokhoz szükséges pontos anyagmennyiség biztosítása általában térfogat szerinti, vagy tömeg szerinti adagolással történik. A kapcsolódó technológiai, illetve szállító berendezés [2.] közvetlen hatást gyakorol az elzáró-szerkezet és a mérő-adagoló szerkezet kialakítására. Ennek következtében magát az adagolási folyamatot is nagymértékben befolyásolja. Az értekezés korlátai nem teszik lehetővé az alkalmazható eszközök teljes vertikumának áttekintését. Az adagolási folyamat működése függ az alkalmazott megoldásoktól és a figyelembe vett elemektől. A folyamat szempontjából meg kell különböztetni a közvetlen adagolást (az elzáró-szerkezet végzi az adagolást), a közvetett adagolást (elzáró és adagoló berendezést is tartalmaz) és a kombinált adagolást (az adagoló berendezés látja el a silózár szerepét is).
1.1. Közvetlen adagolás Közvetlen adagolásról akkor beszélünk, ha az adagolási folyamat közvetlenül vezérelhető az elzáró elem mozgatásával (1/2. ábra), és annak működése határozza meg a kiömlő anyagmennyiséget. Tárolt ömlesztett anyag Munkahenger
Szállítótartály
Adagolt ömlesztett anyag
1/2. ábra Példa közvetlen adagolásra A szerkezet működése egyszerű, külön adagoló-szerkezet, vagy mérleg alkalmazása nélkül megvalósítható, de rendkívül pontos méretezést és működtetést igényel. A működés során az
- 13 -
adagoló elem mozgása előre meghatározott időkhöz igazodik, így a vezérléséhez speciális vezérlő mechanizmus, vagy időmérő berendezés szükséges. Közvetlen adagolást általában szakaszos anyagáramlás esetén alkalmaznak [4.]. A gyakorlatban a pneumatikus működtetésű közvetlen adagolás csak ritkán fordul elő. Ennek oka, hogy a pneumatikus munkahengerek gyakorlati méretezése a dugattyúsebesség és ciklusidő tekintetében nem egzakt.
1.2. Közvetett adagolás Közvetett adagolásról akkor beszélhetünk, ha az adagolt ömlesztett anyag mennyiségét nem az elzáró-szerkezet mozgatásával állítjuk elő, hanem egy külső adagoló-szerkezettel, amely lehet kihordó berendezés, mérleg-vezérelt adagoló berendezés. A kihordó berendezés [5.] általában egy folyamatos működésű anyagmozgató eszköz (pl. szállítószalag, szállítócsiga, vibrációs csatorna, stb.), amit közvetlenül a siló kifolyónyílása alá építenek be (1/3. ábra). Az adagolást a kihordó berendezés végzi, a mennyiséget a berendezés paraméterei és működési jellemzői határozzák meg. Az elzáró-szerkezetek szerepe ennél az adagolási módnál csak a siló lezárására korlátozódik. A közvetett adagolás előnye, hogy az ömlesztett anyag kiömlési jellemzői kevéssé befolyásolják az adagolási folyamatot. Kihordó berendezések alkalmazásával folyamatos adagolás is megvalósítható. Tárolt ömlesztett anyag
Adagolt ömlesztett anyag Szállítószalag
1/3. ábra Példa kihordó berendezésre Mérleg-vezérelt adagolásnál (1/4. ábra) a silóból kiömlő anyag egy mérlegre kerül, amely a kívánt mennyiség elérésekor, az elzáró-szerkezet mozgatásával bezárja a kiömlőnyílást [4.]. Tárolt ömlesztett anyag Munkahenger Szállítótartály Adagolt ömlesztett anyag Szakaszos működésű anyagmozgató berendezés
Mérleg
1/4. ábra A mérleg-vezérelt adagolás egy lehetséges megoldása
- 14 -
A szerkezet kialakítása és a lemért mennyiség továbbítása a csatlakozó anyagmozgató berendezésektől és a technológiai követelményektől függ. Legfőbb előnye, hogy elméletileg tetszőleges mennyiségek szakaszos adagolására alkalmas, de ehhez vezérlőkör alkalmazását igényli (1/5. ábra). Adagolt anyagmennyiség
Mérleg
Érzékelő szenzor
Adagolást indító jel
1/5. ábra Mérleg-vezérelt adagolás vezérlésének vázlata
1.3. Kombinált adagolás Kombinált adagolásról akkor beszélünk, ha a kiömlőnyílásra szerelt adagolószerkezet egyidejűleg ellátja az elzáró-szerkezet szerepét is. Ez általában azt jelenti, hogy a szerkezet kikapcsolt állapotban teljesen lezárja a siló kiömlőnyílását. A kombinált adagoló berendezések a silóból egy adott térfogatú anyagmennyiséget vesznek ki. Ezekkel az eszközökkel szakaszos és folyamatos adagolást is megvalósíthatunk, a működtetés jellegétől függően. A folyamat során egy adagoló-szekrény, vagy adagolócella választja le a szükséges mennyiséget, forgó (pl. forgócellás adagoló), vagy alternáló (pl. dugattyús adagoló) mozgással [6.]. A kombinált adagolás egy lehetséges megvalósításának elve látható az 1/6. ábrán. Tárolt ömlesztett anyag
1. munkahenger
Adagolószekrény 2. munkahenger Adagolt ömlesztett anyag
1/6. ábra Példa kombinált adagolásra A kombinált adagolás előnye, hogy egyszerű folyamat, nincs szükség mérlegre, vagy változó értékhez igazodó vezérlésre, de szakaszos adagolásnál csak diszkrét, előre meghatározott mennyiség ürítésére alkalmas.
- 15 -
1.4. Az adagolt mennyiség Az adagolási folyamat legfontosabb paramétere az adagolt mennyiség. Ez a teljes adagolási folyamatot tekintve állandó érték és a kapcsolódó technológiai folyamat határozza meg. Az adagolt anyagmennyiség eltérő módon határozható meg szakaszos illetve folyamatos adagolás esetén. Szakaszos adagolásról akkor beszélünk, ha az anyag diszkrét adagokban kerül ki a silóból. Az adagolási folyamat során az egyes adagok mennyisége és követési ideje is eltérhet [1.]. Folyamatos adagolás esetén az anyag megszakítás nélkül áramlik ki a silóból. Az adagolt mennyiség nem adagokban, hanem folytonos anyagáram (pl. [kg/s]) formájában áll rendelkezésre [1.]. A két esetben eltérő módon lehet meghatározni az adagolt mennyiséget. Az elzárószerkezettel végzett adagolás alkalmazása (közvetlen és mérleg vezérelt adagolás) leginkább a szakaszos folyamatok esetén előnyös, ezért az értekezésben csak ennek bemutatására szorítkozom. Az elzáró-szerkezettel végzett adagolási ciklust három, egymástól határozottan elkülönülő fázisra bonthatjuk (1/7. ábra): 1. fázis: az elzáró-szerkezet nyitása 2. fázis: nyitott állapotban történő anyagkiáramlás 3. fázis: az elzáró-szerkezet zárása qa [kg/s]
tn
tt
tz
nyitás
nyitvatartás
zárás
t [s]
1/7. ábra A kiömlő anyagmennyiség alakulása az egyes fázisok alatt Az 1. és 3. fázisban a kiáramló mennyiség folyamatosan változik, nyitás esetén növekszik, zárás esetén csökken az anyagáram, a tényleges kiömlési keresztmetszet függvényében. A 2. fázisban az anyagkiáramlás állandónak tekinthető (ha az anyagok viselkedésének és a siló kialakításának speciális eseteitől eltekintünk). A három fázis alatt összesen kiáramló mennyiségre igaz, hogy Qa Qn Qt Qz Qs ,
ahol Qn Qt Qz Qs
- a nyitás alatt kiáramló anyagmennyiség [m3], - a nyitva tartás alatt kiáramló anyagmennyiség [m3], - a zárás alatt kiáramló anyagmennyiség [m3], - a technológia által igényelt anyagmennyiség [m3].
- 16 -
(1.1)
Az egyes fázisok alatt kiáramló anyagmennyiségek meghatározása: tn
Qn q a t dt ,
(1.2)
0
Qt
tn tt
q a t dt ,
(1.3)
tn
Qz
tn tt tz
q a t dt .
(1.4)
tn tt
ahol qa tn tt tz
- a kiömlő anyagáram [m3/s], - a nyitás ideje [s], - a nyitvatartás ideje [s], - a zárás ideje [s].
A kiömlő anyagáram meghatározása (1/8. ábra):
q a t A n t va t ,
(1.5)
ahol An va
- a kifolyónyílás nyitott keresztmetszete [m2], - az áramló anyag átlagsebessége [m/s].
Fa
va Az
sd
{
An
FR
qa
1/8. ábra A kiömlő anyagáram értelmezése vízszintes tolózárnál Ha a kiömlési sebességet állandónak tekintjük, akkor a nyitvatartási fázis alatt kiáramló anyagmennyiség egyszerűen számítható: Qt qa t t Amax va t t .
(1.6)
ahol Amax
- a kifolyónyílás maximális keresztmetszete [m2].
A nyitás és a zárás alatt kiáramló anyagmennyiség a kiömlési keresztmetszettől függ, amit a munkahenger, az ömlesztett anyag és az elzáró-szerkezet geometriája együtt határoz meg. A kiömlési keresztmetszet figyelembevételére két módszer adódik: 1. Az adagolt mennyiség közelítése 2. Komplex modell felírása
- 17 -
1.4.1. Az adagolt mennyiség közelítése A gyakorlatban elfogadott módszer a kiömlő mennyiség meghatározására az állandó kiömlési keresztmetszet és kiömlési sebesség figyelembevétele [7.]. Ebben az esetben a teljes adagolási folyamat alatt kiömlő mennyiség: Qa
A max va t n t t t z . b
(1.7)
ahol b
- az elzáró-szerkezet működési sebességétől függő tényező.
A b tényező értéke attól függ, hogy az egyes adagolási fázisok ideje hogyan viszonyul egymáshoz. Ha az elzáró-szerkezet működése gyors, tehát teljesül, hogy
tt
tn , tz ,
(1.8)
akkor a maximális kiömlési keresztmetszetet veszik figyelembe (b=1):
Qa Amax va t n t t t z .
(1.9)
Ha az elzáró-szerkezet működése lassú, tehát az (1.8) feltétel nem teljesül, akkor a teljes adagolási folyamat alatt kiömlő mennyiség meghatározásánál a maximális kiömlési keresztmetszet felét veszik figyelembe: Qa
A max va t n t t t z . 2
(1.10)
A b értéke tapasztalati megfontolások alapján más értékeket is felvehet. Hátránya ennek az eljárásnak, hogy adagolt mennyiség meghatározására csak durva becslést ad. 1.4.2. Az adagolás komplex összefüggés rendszere Az adagolási folyamat komplex leírásához át kell tekinteni az egyes elemek közötti kapcsolatokat és a befolyásoló tényezőket. A kifolyónyílás nyitott keresztmetszetét meghatározó tényezők: An t f1 sd t ; es ,
(1.11)
ahol sd es
- a pneumatikus munkahenger dugattyújának elmozdulása [m], - az elzáró-szerkezet kialakítási jellege.
A dugattyúelmozdulást meghatározó tényezők: sd t f 2 FR t ; mh
ahol FR mh
- a dugattyúra ható rúderő [N]. - a pneumatikus munkahenger kialakítási, működtetési jellege.
- 18 -
(1.12)
A dugattyúra ható rúderőt meghatározó tényezők: FR t f3 Fa t ; es
(1.13)
ahol Fa
- az ömlesztett anyag nyomásából a zárólapra ható erő [N].
Az ömlesztett anyag nyomásából a zárólapra ható erőt meghatározó tényezők: Fa t f 4 az t ; A z t
(1.14)
ahol az Az
- az ömlesztett anyagban ébredő feszültség [Pa], - a kifolyónyílás zárt keresztmetszetének mérete [m2].
A folyadékokkal ellentétben, ahol a nyomás minden irányban azonos, ömlesztett anyagok esetén a belső feszültség értéke horizontális és vertikális irányban eltér. Ez okból az (1.14) egyenletben figyelembe vett nyomásérték (az) a horizontális (azh) és a vertikális (azv) feszültségeknek az elzáró-szerkezet geometriája (es) által meghatározott eredője lesz. Az anyag kiáramlási sebessége, az adagolt ömlesztett anyag tulajdonságainak és a siló geometriai jellemzőinek függvénye:
va t f5 a ; s
(1.15)
ahol a s
- az ömlesztett anyag jellege. - a siló kialakítási jellege.
Az ömlesztett anyagok sajátosságainak következtében a kiáramlási sebesség zérus értéket is felvehet - elsősorban kohézív anyagoknál, az adagolási folyamat kezdetén (beboltozódás) -, amit csak speciális eszközök alkalmazásával (pl. bolygatás) lehet megszüntetni. Az értekezésemben kidolgozott adagolási modell olyan anyagokra vonatkozik, amelyeknél nem áll fenn a boltozatképződés veszélye. Figyelembe véve az (1.11-1.15) függvénykapcsolatokat, megállapítható, hogy a kiáramló anyagmennyiségre a következő összefüggés írható fel: qa t f6 Fa (t); es ; mh ; a ; s ; az t ; A n t ; A z t .
(1.16)
Az (1.16) függvénykapcsolat kibontásához külön-külön meg kell határozni az ömlesztettanyagokra, a pneumatikus munkahengerekre és az egyes elzáró-szerkezet kialakításokra vonatkozó összefüggéseket. Értekezésem célja az ehhez szükséges általános modell felállítása.
- 19 -
2. SILÓK ELZÁRÓ-SZERKEZETEI Ömlesztett anyagok silóban történő tárolásának legfőbb előnyei a zárt tárolótér és az egyszerű, jól automatizálható adagolás. Közvetlen és mérleg-vezérelt adagolásnál, külső anyagmozgató vagy adagoló berendezés alkalmazása nélkül, kizárólag az elzáró-szerkezet mozgatásával, a gravitációs erő kihasználásával történik az anyagok ürítése. Az adagolás fontos jellemzője a kiömlési mód, amely lehet [8.]: - szabad kiömlés, - nyitott csatorna (vályú), - zárt csővezeték (függőleges, vagy ferde). Szabad kiömlés akkor valósul meg, ha az áramló anyagban csak vertikális sebességkomponens ébred és horizontális irányú sebesség nem lép fel. A másik két esetben a szerkezet, az irányváltoztatás és a falsúrlódás révén módosítja az áramlást. Kutatásaim során csak a szabad kiömlést biztosító elzáró-szerkezetekkel foglalkozom. Az elzáró-szerkezetek leggyakoribb megoldásai [4.]: - tolózárak, - szegmens-zárak, - csuklós zárak, - billenő surrantók, - függőláncos elzárók, - gereblyés zárak, stb. Mivel csuklós zárak, billenő surrantók, függőláncos elzárók és gereblyés zárak esetén a nyitási, illetve zárási fázis alatt (surrantóknál és gereblyés záraknál teljesen nyitott állapotban is) az elzáró elem kialakítása határozza meg az anyag áramlásának útját, ezért nem beszélhetünk szabad kiömlésről. Emiatt az adagolási folyamat modellezésénél csak a tolózáras és a szegmens-záras elzáró-szerkezeteket vizsgálom. Az elzáró-szerkezetek működtetési módjuk szerint lehetnek [2.] mechanikus, hidraulikus, vagy pneumatikus működtetésűek. A mechanikus működtetés általában egy elektromos motor és egy kényszerkapcsolat segítségével történik (pl. fogaskerék-fogasléc), hidraulikus és pneumatikus működtetés esetén pedig egy munkahenger mozgatja az elzáró elemet. A legtöbb szakaszos elzáró-szerkezet mindhárom módon működtetethető, mindegyik megoldásnak megvannak az előnyei, illetve hátrányai, és ehhez igazodva a jellemző felhasználási területei. Értekezésemben csak a pneumatikus működtetéssel foglalkozom.
2.1. Tolózárak Tolózáraknál a kiömlés nyitása, illetve zárása egy síklap lineáris, alternáló mozgatásával történik és a tényleges kifolyási keresztmetszet a mozgás minden fázisában merőleges a zárólap síkjára (2/1. ábra). A kiömlési jellemzők (elmozdulás, felület, stb.) általában egyszerű módon számíthatók a működtetési mód meghatározó paraméterének (pneumatikus mozgatásnál a munkahenger elmozdulása) függvényében. Technikailag megvalósítható a tolózáras kialakítás két, egymással szemben, szinkronizáltan mozgó síklappal is, de ennek gyakorlati jelentőssége, a bonyolult kialakítás miatt elhanyagolható. A zárólap mozgási síkjának függvényében megkülönböztethetünk vízszintes, függőleges és ferde tolózárat. - 20 -
Zárólap
Siló
Működtető munkahenger Zárólap
Zárólap
2/1. ábra Vízszintes tolózár vázlata Vízszintes tolózáraknál (2/1. ábra) az anyag kiömlése kizárólag függőleges irányú, az anyag által a zárólapot terhelő erő számításánál csak az ömlesztett anyagban ébredő vertikális feszültséget kell figyelembe venni [4.]. Ferde tolózáraknál (2/2. ábra) a zárólap helyzete e szögben döntött a vízszinteshez képest, ezért a mozgás két koordináta irányában történik, ami a szájnyílás kialakításától függően befolyásolja az anyag kiömlését.
Siló Működtető munkahenger
Zárólap
e
2/2. ábra Ferde tolózár vázlata Ferde tolózáraknál a szájnyílás kialakítása többféle lehet, de a gyakorlatban csak két megoldásnak (toldott és a hosszabbított) van jelentősége (2/3. ábra). Toldás esetén egy henger, vagy hasáb alakú csatornával illesztik a tolózárat a kifolyótölcsérhez, hosszabbításnál a kifolyótölcsért „szabják” a tolózár helyzetéhez. A kétféle kialakítás eltérő befolyást gyakorol a kiömlési tulajdonságokra. Hosszabbítás esetén a silógarat nem szimmetrikus, ami a kiömlés során kismértékű (e függvényében) horizontális irányú áramlási sebességet okoz, toldásnál csak vertikális áramlási sebességgel kell számolni. - 21 -
Mivel az értekezésben csak a szabad kiömléssel foglalkozom, ezért a továbbiakban a toldott változatot vizsgálom.
An = Akh
An = A k
Akh
Akh e
e
Ak
Ak
b) hosszabbított
a) toldott
2/3. ábra Ferde tolózár szájnyílás kialakításai A toldás kismértékben megnöveli az anyagban fellépő feszültségeket, de ez nagy anyagmagasságnál - Janssen elmélete szerint [9.] - nem számottevő. Ferde tolózáraknál, a zárólap mozgatása során figyelembe kell venni a horizontális és a vertikális feszültségből adódó terhelőerőt is. Függőleges tolózáraknál (2/4. ábra) a zárólap helyzete függőleges, így az áramló anyag horizontális és vertikális irányú sebességgel is rendelkezik [5.]. Működtető munkahenger
Siló
Zárólap
2/4. ábra Függőleges tolózár vázlata A szájnyílás kialakítását tekintve itt is két megoldásnak (toldott és a hosszabbított) van jelentősége (2/5. ábra), de a gyakorlatban mindkét megoldásnál szükség lehet egy oldalirányú toldat (a siló és a zárólap geometriájának függvényében) alkalmazására (lásd 2/4. ábra), amely a szájnyílásban jelentősen módosíthatja a kiömlési viszonyokat [4.]. Függőleges, illetve ferde tolózáraknál a zárólap a nyitás folyamán általában felfelé mozog. Lefelé, vagy oldalra mozgatás esetén a tolózár akadályozza az anyag áramlását és a szerkezet tömítése nehezebben oldható meg. Emellett oldalra történő mozgatásnál az áramlási viszonyok nem homogének.
- 22 -
b) hosszabbított
a) toldott
2/5. ábra Függőleges tolózár szájnyílás kialakításai
2.2. Szegmens-zárak Szegmens-záraknál (2/6. ábra) a kiömlés nyitása, illetve zárása egy, vagy két íves kialakítású záró-szegmens, adott csuklópont körüli alternáló mozgatásával történik [4.]. A kiömlést meghatározó paraméterek (elmozdulás, felület, stb.) felírásához szükséges a mozgatómechanizmus pontos ismerete és a működtetési mód jellemző paraméterével való kapcsolatának definiálása. A szegmensek felülete általában hengerpalást és nem gömbfelület. A szegmens-zárakon belül három alapvető kialakítási változat fordul elő: egytagú: o szimmetrikus, o aszimmetrikus, kéttagú.
Siló
Működtető munkahenger Zárószegmens a) szegmens elhelyezés
b) szájnyílás kialakítás
2/6. ábra Egytagú, szimmetrikus szegmens-zár Egytagú, szimmetrikus szegmens-záraknál (2/6. ábra) a tényleges kifolyási keresztmetszet nem vízszintes, ezért az anyag kiömlésénél kismértékű horizontális sebesség is fellép (nyitott állapotban nincs horizontális irányú áramlás) [4.]. A horizontális sebesség komponens értéke nem számottevő, nagy szegmens-sugár esetén elhanyagolható, de e tényező hatását értekezésemben nem elemzem. Aszimmetrikus kialakítás esetén (2/7. ábra) a szegmens elhelyezésére számos megoldás elképzelhető toldat alkalmazásával, vagy anélkül (2/8. ábra), illetve a szegmens csuklópontjának helyzete is széles határok között változtatható [2.]. Az adott feladatra
- 23 -
alkalmas megoldás kiválasztása az alkalmazási sajátosságok - elsősorban a geometriai és méretkorlátok - függvénye.
Siló
Működtető munkahenger
Zárószegmens a) szegmens elhelyezés
b) szájnyílás kialakítás
2/7. ábra Egytagú, aszimmetrikus szegmens-zár Aszimmetrikus szegmens-elhelyezésnél a tényleges kifolyási keresztmetszet a kiömlés minden fázisában ferde irányú (2/8. ábra), alakja pedig a szájnyílás kialakításának függvénye (bizonyos esetekben bonyolult felület). A kiömlés paramétereinek meghatározásánál az áramlási sebesség mindkét komponensét figyelembe kell venni.
a) toldat nélkül
b) vertikális toldat
c) oldalsó toldat
2/8. ábra Szájnyílás kialakítási lehetőségek aszimmetrikus szegmenseknél A leírtak alapján az aszimmetrikus szegmens-zár sem alkalmas az általam definiált adagolási modellben történő alkalmazásra. Kéttagú szegmens-záraknál [4.] a tényleges kifolyási keresztmetszet a kiömlés minden fázisában a szájnyílás nyitott részének vízszintes vetületével számítható és az anyag áramlási sebessége csak vertikális irányú komponenssel rendelkezik (2/9. ábra). Pneumatikus munkahengerrel történő működtetés esetén a tényleges nyitási sebesség a szegmensek szimmetrikus mozgása miatt a dugattyúsebesség kétszerese. A két szegmens mozgásának összehangolására különféle megoldások alkalmazhatók (pl. fogaskerekes, fogas íves, stb.), de mindenképpen biztosítani kell a csúszásmentes kényszerkapcsolatot és a szegmensek szinkronizált mozgását. A dugattyúerő meghatározásánál, mindhárom szegmens-zár esetén, figyelembe kell venni az anyagban ébredő horizontális és vertikális feszültségeket is. Az ismertetett megfontolások alapján, vizsgálataimnál a szegmens-zárak közül csak a kéttagú kialakítással foglalkozom részletesen, amely jól szabályozható, centrikus anyagáramot biztosít [4.].
- 24 -
Siló
Működtető munkahenger
Zárószegmensek a) szegmens elhelyezés
b) szájnyílás kialakítás
2/9. ábra Kéttagú szegmens-zár
2.3. Elzáró-szerkezetek méretezése Tároló silók elzáró-szerkezeteit a szakirodalomban három paraméterre vonatkozóan méretezik: geometriai méretek, szerkezeti terhelés, működtető erőszükséglet. A geometriai méreteket a siló kiömlőnyílásának mérete és alakja, valamint az elzáró-szerkezet kialakítása határozza meg. Az elzáró-szerkezetek szilárdsági méretezése általában a tárolt ömlesztett anyag statikus terhelésének felhasználásával történik. Az adagolási folyamat szempontjából a működtető erőszükséglet a legfontosabb paraméter, amelyre vonatkozóan a magyar szakirodalomban Andrási és Polonszky [4.] a ferde tolózárakra, a szegmenszárakra, a gereblyés zárakra és a függőláncos elzárókra, Prezenszki [5.] a vízszintes tolózárakra és a szegmenszárakra, Knoll [7.] csak a függőleges tolózárakra ad meg összefüggéseket. Mindhárom szerző hasonló megközelítést alkalmaz és csak a maximális működtető erő meghatározására szorítkozik. Az általános, ferde tolózárakra vonatkozó összefüggést a maximális működtető erő meghatározására (a zárólap csúszó megtámasztása esetén) Andrási és Polonszky [4.] írta le: Fmax k t Fe x 1 Ge 1 cos e sin e ,
(2.1)
ahol kt Ge e x 1 Fe
- biztonsági tényező (a szemcsenagyság függvénye kt=1,25…1,5), - az elzáró-szerkezet mozgó elemeinek önsúlya [N], - a tolózár dőlésszöge [o], - falsúrlódási tényező az anyag és a zárólap között, - a zárólap és támasztófelülete közötti súrlódási tényező, - az anyagból a zárólapra ható eredő terhelés [N]:
Fe v Ae cos2 e sin 2 e ,
v Ae
- az anyagból az elzáró-szerkezetre ható függőleges irányú nyomás [Pa], - nyomásviszony az ömlesztett anyagban, - az elzáró-szerkezet anyaggal érintkező felülete [m2]. - 25 -
(2.2)
Vízszintes és függőleges kialakítás esetén a (2.1) összefüggés egyszerűsödik, e=0o, illetve e=90o behelyettesítésével. A zárólap görgős megtámasztása esetén a (2.1) összefüggésben a 1 súrlódási tényező helyett, Knoll [7.] nyomán, a gördülő-ellenállási tényezőt kell alkalmazni:
g
c d g 2 lg Dg
,
(2.3)
ahol dg Dg c lg
- a görgőcsapok átmérője [m], - a görgők átmérője [m], - a görgőcsapok súrlódási tényezője, - a gördülő-ellenállás karja [m].
Kézi működtetésű, ellensúlyos szegmenszárak esetén a maximális működtető erő meghatározását részletesen szintén Andrási és Polonszky [4.] adta meg, egy szegmensre vonatkozóan: Fmax
Fe x R z G es R es Fe G e G es Re
dc c 2 ,
(2.4)
ahol Rz Res Ges dc Re
- az elzáró szegmens sugara [m], - az ellensúly karja [m], - az ellensúly nagysága [N], - a csapátmérő [m], - a működtető erő karja [m].
Két szegmens esetén is hasonlóan számítható működtető erő, de az anyagáram jobban kezelhető [4.] A fent leírt, Andrási és Polonszky [4.] által megadott, illetve a többi szakirodalomban található [5., 7.] összefüggések mechanikus működtetésű (fogasleces, csavarorsós, kézi karos, stb.) szerkezetekre vonatkoznak, ezért nem tudják figyelembe venni a pneumatikus működtetés sajátosságait. Mivel pneumatikus működtetésű elzáró-szerkezetekkel kapcsolatos szakirodalom nem ismert, ezért értekezésem feladata az adagolási modell felállításához szükséges, tolózárakra és kéttagú szegmenses-zárakra vonatkozó összefüggések felírása.
- 26 -
3. ÖMLESZTETT ANYAGOK ADAGOLÁSI JELLEMZŐI Ömlesztett anyagoknak nevezzük azokat a rendszerint különböző szemnagyságú, szilárd részeket tartalmazó, de általában egynemű anyagokat, amelyeket nagyobb tömegben, rendezetlenül és csomagolatlanul mozgatnak, tárolnak, szállítanak (pl. homok, ércek, szén, gabonafélék, csavaranyák, stb.) [2.]. Noha első ránézésre az ömlesztett anyagok viselkedése a folyadékokéra hasonlít, számos jellemzőjük jelentősen eltérő jelleget mutat. Ha röviden akarjuk jellemezni az ömlesztett anyagokat, azt mondhatjuk, hogy egyes tulajdonságaikban folyadékként, másokban szilárd anyagokként jellemezhetők. Ennek oka az, hogy az ömlesztett anyagok általában különböző méretű (tulajdonságú) szilárd szemcsékből állnak, amelyek adott feszültségviszonyok között (terhelés) szilárd tömbbé állhatnak össze, más feszültségviszonyok között egymáson könnyen elmozdulhatnak (3/1. ábra).
3/1. ábra Példa egy komplex ömlesztett anyag összetételére Az ömlesztett anyagok legfontosabb tulajdonságai a következők [8.]: mechanikai jellemzők: o belső súrlódás, o kohézió, o adhézió, o külső súrlódás. o rézsűszög, fizikai jellemzők: o fajsúly, o nedvességtartalom, o a szemcsék tulajdonságai: alak, méret, keménység, összetétel, porozitás, o felület, o tixotrópia, o homogenitás.
- 27 -
A különböző mechanikai és fizikai jellemzők együttesen határozzák meg az adott ömlesztett anyag tárolási, illetve kiömlési tulajdonságait, amelyek az anyagban, illetve az anyag és adott felületek között ébredő feszültségekkel jellemezhetők. Bizonyos technikai megoldások segítségével (pl. pneumatikus bolygatás), az ömlesztett anyagok tárolási és kiömlési tulajdonságai külső eszközökkel is befolyásolhatók. Az ömlesztett anyagok paramétereinek leírására, jelentősen eltérő fizikai tulajdonságaik és folyási viselkedésük miatt, nem létezik olyan általános összefüggés-rendszer, amely minden anyagtípus esetén alkalmazható. Az adagolási folyamat modelljéhez szükséges paraméterek meghatározásánál figyelembe kell venni, hogy az ömlesztett anyagok eltérő viselkedést mutatnak statikus (tárolt, nem mozgó anyaghalmaz) és mozgási (kiömlés, mozgatás) állapotban. Ez azt jelenti, hogy az adagolási folyamat során más összefüggéseket kell alkalmazni az elzáró-szerkezet bezárt (statikus állapot) és nyitott állapotában (mozgási állapot). Fontos feltétel, hogy az adagolás során az anyagáramlást folytonosnak és akadálymentesnek tekintem, tehát nem veszem figyelembe a szemcsék, illetve a silófal közötti esetleges összetapadásokat és a beboltozódásból adódó problémákat. Ez a feltétel száraz, kohéziómentes, szemcsés anyagok adagolása esetén teljesül.
3.1. Statikus állapot-jellemzők A szakirodalomban különböző eljárások terjedtek el az ömlesztett anyagok viselkedésének leírására, amelyek eltérő szemszögből közelítik meg a szemcsék közötti kapcsolatot. Széles körben alkalmazott eljárás a Janssen [9.] által felírt modell, amely a siló hengeres részén adja meg az anyagban ébredő statikus feszültségek meghatározására szolgáló összefüggéseket. A kifolyógaratban érvényes összefüggéseket Janssen elméletéből kiindulva Schulze [10.] írta le. A modell alkalmazhatósága korlátozott és nem ad egzakt megoldást, de legnagyobb előnye, hogy analitikus összefüggéseket ad a feszültségek leírására. A magyar szakirodalomban is ez az eljárás a legelterjedtebb [11.]. Egy másik megközelítés, amikor az ömlesztett anyagok viselkedését az egyes szemcsék kapcsolata és egymásra gyakorolt hatása alapján próbálják leírni. Ez az ún. diszkrét elemek módszere (DEM) [12.]. A módszer elméletileg egzakt megoldást eredményez, de a gyakorlatban nagyon nehéz pontosan definiálni az egyes szemcsék közötti kapcsolatokat, különösen különböző szemcséket tartalmazó anyaghalmaz esetén. Másik hátránya a módszernek, hogy az anyagjellemzők meghatározása numerikus megoldást igényel, ami nagy mennyiségű szemcse esetén rendkívül nagy számítási időt igényelhet. Vizsgálataimnál az ömlesztett anyagok statikus jellemzőinek meghatározására a Janssen-féle megközelítést alkalmaztam, mert a DEM eljárás illesztése a pneumatikus munkahenger működését leíró numerikus módszerhez túllép az értekezés célkitűzésein. Janssen [9.] összefüggő halmazként kezeli az ömlesztett-anyagot és e halmaz egy kicsiny térfogatelemének a többi anyagrésszel való kapcsolatát vizsgálja. A modell alapvetően karcsú, hengeres silókra vonatkozik és a siló hengeres részén ad összefüggést az anyagban ébredő feszültségek leírására. A silóban fellépő feszültségek szempontjából, a silómagasság (Hs) és a silóátmérő (Ds) viszonya alapján a Magyar Szabvány (MSZ EN 1991-4:2006) [13.] a silókat különböző kategóriákba sorolja: lapos silók H s / Ds 0,4 , zömök silók 0,4 H s / Ds 1 , közepesen karcsú silók 1 H s / Ds 2 , karcsú silók H s / Ds 2 .
- 28 -
A magasság/átmérő viszony jelentősége abból adódik, hogy a Janssen által felvett görbe alakja lapos silók esetén közelíthető a hidrosztatikus nyomás-görbével, karcsú silóknál a falsúrlódás miatt a fenéknyomás lényegesen kisebb lesz a hidrosztatikus nyomásnál. Az anyagmagasság növelésével állandósuló feszültség miatt a karcsú silók alkalmazása előnyösebb, de a áramlás megindulásakor, a silószár és a kifolyógarat találkozási vonalán fellépő feszültségugrás értéke az átmérő/magasság viszonnyal arányosan nő [11.]. A feszültségek meghatározásához a silót két részre kell osztani (3/2. ábra): silószár (a siló henger vagy hasáb alakú, felső része), kifolyógarat (a siló alsó, általában kúp vagy gúla alakú része).
Silószár Hs
Ds
Kifolyógarat
ht
Dk
h0
3/2. ábra A siló felépítése A kifolyás szempontjából a garat a mérvadó, de a jellemzők meghatározásához a silószárban fellépő feszültségek is szükségesek (3/2. ábra). Janssen elméletét felhasználva a hengeres silószárban fellépő feszültség-viszonyok a 3/3. ábrán látható modell segítségével írhatók le [10.]. vs As
z
ws Us dzs dzs
hs Us dzs (vs + dvs) As g a As dzs
3/3. ábra Korong alakú térfogatelem a silószárban [10.]
- 29 -
A 3/3. ábra alapján felírható a vertikális erőegyensúly a silószárban:
As vs g a As dzs As vs dvs ws Us dzs
(3.1)
ahol vs ws As dzs a g Us
- vertikális feszültség a silószár adott térfogatelemén [Pa], - nyírófeszültség a siló hengeres oldalfalán [Pa], - a silószár keresztmetszete [m2], - a hengeres térfogatelem magassága [m], - az ömlesztett anyag sűrűsége [kg/m3], - a gravitációs gyorsulás [m/s2], - a hengeres térfogatelem kerülete [m].
A (3.1) egyenletből kiindulva felírható a vertikális irányú feszültség (ha vs0=0, azaz az anyag felszínén nincs terhelés): tan x Us z g a As As vs 1 e (3.2) , tan x Us és a horizontális irányú feszültség: g a As hs tan x Us
tan x Us z As 1 e ,
(3.3)
ahol x
- a falsúrlódási szög.
A kúpos kifolyógaratban fellépő feszültség-viszonyokat, szintén a Janssen-féle megközelítésből kiindulva (Schulze nyomán) [10.] lehet meghatározni a 3/4. ábrán látható modell segítségével. vt At + d(At vt) wt dAM
dzt
wt dAM vt At g a At dzt z
3/4. ábra Kúp alakú térfogatelem a kifolyógaratban [10.] A 3/4. ábra alapján felírható a vertikális erőegyensúly a kifolyógaratban:
At vt d At vt g a At dz t At vt dAM wt cos dAM wt sin ahol vt wt wt
- vertikális feszültség a kifolyógarat adott térfogatelemén [Pa], - a kifolyógarat oldalfalára merőleges feszültség [Pa], - nyírófeszültség a kifolyógarat oldalfalán [Pa], - 30 -
(3.4)
At dzt dAM
- a kifolyógarat dőlésszöge [o], - a kifolyógarat keresztmetszete a kúp alakú térfogatelemen [m2], - a kúp alakú térfogatelem magassága [m], - a kúp alakú térfogatelem palástjának felülete [m2].
A (3.4) egyenletből kiindulva felírható a vertikális irányú feszültség értéke: n
g a z g a h 0 z vt v0 , n 1 n 1 h0
(3.5)
és a horizontális irányban ható feszültség: n
g a z g a h 0 z ht v0 , n 1 n 1 h0
n mg 1 wt vt
tgx 1 1 , tg
(3.7)
- a kifolyógarat magassága [m], - geometriai jellemző (ék alakú garatra mg=0, kúp (gúla) alakú garatra mg=1), - a kifolyógarat kúpmagassága (vagy gúla magassága) [m]:
h0
v0
Feszülts égek az anyagban [Pa]
ahol ht mg h0
(3.6)
Ds , 2 tg
(3.8)
- a vertikális feszültség a kifolyógarat és a silószár találkozásánál [Pa]:
v0 vs Ha .
(3.9) v Vertikális feszültség Horizontális feszültség h
h
Silószár
Kifolyógarat
3/5. ábra Példa a feszültségek alakulására karcsú silók esetén - 31 -
Any
Az n kitevő értéke az ömlesztett anyag összenyomhatóságának és az elzáró-szerkezet merevségének függvényében változhat (tapasztalati érték) [10.]: összenyomhatatlan anyag és tökéletesen merev szerkezet esetén n=0, közel összenyomhatatlan anyag és merev szerkezet esetén n=0,1, mérsékelten összenyomható anyag és rugalmas szerkezet esetén n=0,45, összenyomható anyag és rugalmas szerkezet esetén n=0,9.
3.2. Állapot-jellemzők kiömlés közben Silóból történő kiömlésnél az ömlesztett anyagok háromfelé viselkedést mutathatnak [13.]: magkifolyás, tömegkifolyás, vegyes kifolyás. Magkifolyás esetén az elzáró-szerkezet nyitásakor csak a szájnyílás felett elhelyezkedő anyag indul el, a többi nem mozdul, csak később, az anyagmagasság csökkenésével kezd kiáramlani (3/6.a. ábra). Tömegkifolyás esetén a teljes anyagmennyiség egyszerre áramlik ki a szájnyíláson, a kifolyógarat által meghatározott geometria szerint (3/6.b. ábra). Vegyes kifolyás esetén a siló egyes részein magkifolyási (a szájnyílás felett), más részein tömegkifolyási állapot alakul ki (3/6.c. ábra).
a) magkifolyás
b) tömegkifolyás
c) vegyes kifolyás
3/6. ábra Ömlesztett anyagok lehetséges kifolyási esetei [13.] Egy ömlesztett anyag kifolyási viselkedését elsősorban a kifolyógarat hajlásszöge (θ), illetve az ömlesztett anyag és a kifolyógarat anyaga közötti falsúrlódás (x) határozza meg. Minél laposabb a garat oldalfalának kialakítása és minél nagyobb a falsúrlódási szög, annál nagyobb lesz a magkifolyás előfordulásának valószínűsége [14.]. A kifolyási állapotot meghatározó paraméterek értékei és kapcsolataik a 3/7. ábrán láthatók. A szabvány [13.] csak a kifolyógarat hajlásszögének és a falsúrlódás értékének függvényében - 32 -
1 0.9
0.7
3. zóna
2. zóna
0.6 0.5 0.4 0.3
0.7
0.8
0.6
1. zóna
0.1 0
10
0.5
0.6 0.5
0.4
0.4
0.3
0.3
0.2 0.1
0.1 20 30 40 Kifolyótölcsér s zöge [o]
50
60
A kifolyógarat hajlásszöge (θ) [ ] → o
a
0.7
0.2
0.2
0
a 0.9
Falsúrlódási tényezõ
0.8
0.8
1
Falsúrlódási tényezõ
Falsúrlódási tényező (x) →
tényező (x) → Falsúrlódási Falsúrlódási tényezõ
1 ábrázolja a különböző tartományokat (3/7.a. ábra), Schulze [14.] a belső súrlódás hatását is 0.9 figyelembe veszi (3/7.b. ábra).
0
0
0
0 30 40 50 60
1. zóna → tömegkifolyás 2. zóna → átmeneti zóna, tömeg- vagy magkifolyás 3. zóna → magkifolyás
a) MSZ EN 1991-4:2006 [13.]
10
20
30
40
50
60
10 20 30 40 Kifolyótölcsér s zöge [o] A kifolyógarat hajlásszöge (θ) [o] → Kifolyótölcsér s zöge [o] 30e=30o o 40e=40 o 50e=50o =60 e 60
b) Schulze [14.]
3/7. ábra Kifolyási állapotok kúpos, vagy gúla alakú garat esetén A kifolyási jelleg alapvetően befolyásolja a kiömlési paramétereket, közvetlen hatást gyakorolva az adagolt anyagmennyiségre. Tervezési szempontból a magkifolyás a legkedvezőbb megoldás. A silóban tárolt anyag feszültségi állapota jelentősen megváltozik az elzáró-szerkezet kinyitásakor. A silószár és a kifolyógarat találkozásánál a siló geometria és a kifolyási jelleg függvényében jelentős nyomásnövekedés lép fel. A jelenség oka egyelőre nem tisztázott, valószínűsíthetően a kifolyónyílás felett elhelyezkedő, az áramlás megindulásakor létrejövő ún. fellazulási zónában fellépő nyomásátrendeződés következtében áll elő [15.]. A nyomásváltozás következtében az anyag statikus állapotában meghatározott feszültség értékek nem érvényesek a mozgási állapotra. Az elzáró-szerkezetet működtető erő meghatározásához csak a szájnyílás (a kifolyógarat alsó keresztmetszete) mentén fellépő feszültségek értékére van szükség, melyre - egzakt megközelítés hiányában - a kutatók különböző összefüggéseket alkalmaznak. A Magyar Szabvány [13.] a siló szerkezeti méretezését szem előtt tartva a garat fala mentén fellépő érintőirányú feszültség meghatározására ad meg összefüggést. Schulze [10.] egy, a silógarat alá elhelyezett szalagos kihordó-berendezésre az anyagból származó terhelés meghatározásánál a Roberts [16.] féle összefüggést aktualizálja. Sitkei [11.] és Németh [17.] a Janssen-féle megközelítésen alapuló, de a kifolyási állapotra felírt modellt alkalmazza az anyag mozgási jellemzőinek meghatározásánál. Kutatásaim során a Janssen-féle megközelítésen alapuló, magkifolyású állapotra felírt összefüggést alkalmaztam [17.]. Választásom oka kettős: egyrészt ez a módszer vertikális irányú szabad kiömlésre vonatkozik, másrészt a kiömlő mennyiség meghatározásánál is vizsgálom e modellt (többek között), ahol a mérésekkel összevetve jó eredményt ad. Az áramló anyag jellemzőinek meghatározásánál a 3/8. ábra alapján, a mozgó anyag egy kicsiny térfogatelemére felírható erőegyensúlyból kell kiindulni [17.]:
Ak m g a Ak dh m Ak m dm wm Uk dh m ahol m
- vertikális feszültség a kiáramló anyag adott térfogatelemén [Pa],
- 33 -
(3.10)
50
m wm Ak dhm Uk
- az álló anyagrészen fellépő horizontális feszültség [Pa], - az álló és a mozgó anyag között ébredő nyírófeszültség [Pa], - a kiömlési keresztmetszet [m2], - a kiáramló térfogatelem magassága [m], - a kiáramló térfogatelem kerülete [m]. Álló anyag
Mozgó anyag
Álló anyag
m wm
dhm
m
b
a
b
wm m
m+dm
3/8. ábra A kiáramlás modellje magkifolyás esetén [17.] A (3.10) egyenletből kiindulva és figyelembe véve, hogy a sebesség az ömlesztett anyag tulajdonságai miatt eltérő lesz a kiömlési keresztmetszet különböző pontjain, Mohr elmélete alapján az átlagos vertikális feszültségre a következő összefüggést kapjuk [17.]: 2h U k A C A k ca cb k a g ca cb m 1 e , 2 Uk tge
ca
(3.11)
1 2 tg e , tge 4 2
(3.12)
1 2 tge . tge
(3.13)
cb
ahol h C e
- az áramló anyagmagasság [m], - az anyag kohéziójának mértéke [Pa], - az anyag belső súrlódási szöge [o].
3.3. Kiömlő anyagmennyiség Ömlesztett anyagok esetén – hasonlóan a folyadékok kiömléséhez – egy egyszerű összefüggésből indulhatunk ki a silóból kiömlő anyagmennyiség meghatározásánál (3/9. ábra): ma An va a , [kg/s].
(3.14)
A problémát az jelenti, hogy a folyadékok állandó sűrűségével és jól definiálható áramlási sebességével szemben az ömlesztett anyagok e két jellemzője számos nehezen kezelhető paraméter függvénye, melyek közül a legfontosabbak: a szemcsék közötti kohéziós erők, súrlódás a mozgó szemcsék és a silófal között,
- 34 -
a szemcsék közötti levegő nyomásváltozása, a környező levegő beáramlása által a szemcsékre kifejtett húzóerők, vagy ellenállások, a siló állapota (nyitott, vagy zárt fedél, anyagmagasság, tarolási idő, stb.). a
h va An
. ma
3/9. ábra Ömlesztett anyag gravitációs kiömlése Mivel a fenti paraméterek eltérő értékeket vehetnek fel az egyes anyagok tulajdonságainak, illetve az adott környezeti feltételeknek megfelelően, ezért a silóból kiömlő anyagmennyiség meghatározására egzakt összefüggés nem ismert. Ennek következtében számos módszer található szakirodalomban, amely többé-kevésbé pontos közelítéssel, gyakorlati tapasztalatokon alapuló összefüggések segítségével adja meg a kiáramló anyagmennyiséget. Az ömlesztett anyagok szabad kiáramlását leíró ismertebb közelítő módszerek, alkalmazhatóságuk szerint két csoportba sorolhatók: szemcsés (dp≥0,3mm), kohéziómentes ömlesztett-anyagokra felírt módszerek: - Hagen-módszer, - Beverloo-módszer, - Brit szabvány, - Crewdson-módszer, - Williams-módszer, - Zanker-módszer, - Gjacev-Keller módszer, - Oldal-módszer, stb.
finom szemcsés (dp<0,3mm), kohézív porokra felírt módszerek: - Janssen-módszer, - Johanson-módszer, - Tomas-módszer, - Carleton-módszer, stb.
Kutatásom egyik fontos feladata az adagolási folyamat során alkalmazható, a kiömlő anyagmennyiség meghatározására szolgáló módszer kiválasztása. Ehhez először áttekintem az ismertebb módszerek jellemzőit és összefüggéseit. Hagen [18.] írta le az egyik legelső szemcsés, kohéziómentes anyagokra vonatkozó alapösszefüggést, amelynek segítségével, adott kiömlőnyílás esetén, a kifolyó anyagmennyiség meghatározható. A homokon elvégzett kísérletekkel ellenőrzött formula a következő: ma a g Dk k p d p
ahol Dk dp
- a kifolyónyílás átmérője [m], - az átlagos szemcseátmérő [m],
- 35 -
2,5
[kg/s],
(3.15)
kp
- a szemcsék alaktényezője.
Beverloo [19.] tovább pontosította a Hagen által felírt, szemcsés, kohéziómentes, szabad folyású ömlesztett anyagokra vonatkozó összefüggést, homokra és 7 különféle növényi magra elvégzett vizsgálatok alapján: ma Cb a g Dk k p d p
2,5
[kg/s],
(3.16)
ahol Cb
- kifolyási tényező.
A Brit Szabvány (ICS 65.020.40) [20.] szerinti összefüggések szemcsés (dp>500m), szabad folyású anyagokra, tömegkifolyás esetén egyszerűen használhatók és figyelembe veszik a kifolyó tölcsér dőlésszögét is, de csak durva, nagyságrendi közelítést adnak: • kör alakú kifolyónyílásra:
ma 0,58 a g Dk k p d p k [kg/s], 2,5
•
(3.17)
rés (Lk>3∙Wk) alakú kifolyónyílásra:
ma 1,03 a g Lk k p d p Wk k p d p k [kg/s], 1,5
ahol k Wk Lk
(3.18)
- a silógarat geometriai tényezője (k=tan-0,35, ha <45o, egyébként k=1), - a kifolyónyílás szélessége (téglalap alak) [m], - a kifolyónyílás hossza (téglalap alak) [m].
A Crewdson-módszer [21.] szemcsés, kohéziómentes, szabad folyású ömlesztett anyagokra vonatkozik, tömegkifolyás és homogén anyagáramlás mellett, de figyelembe veszi a szemcsék közé beáramló levegő hatását is:
ma Cb a g
2,5 1 dp Dk k p d p [kg/s], a dr
(3.19)
ahol dp/dr - a kifolyógaratban, a szemcsék között kialakuló nyomásváltozás a kifolyógarat csúcspontjától mért távolság függvényében. A Williams-módszer [22.] figyelembe veszi az anyag és a silófal között fellépő súrlódás hatását is, ami kis nyílásméretnél jelentős lehet (Dk<20mm). A módszer nagy szemcséket tartalmazó, homogén anyagok tömegkifolyású ürítésénél alkalmazható hatásosan, de nem számol a szemcsék közötti levegő hatásával. A kifolyó mennyiség meghatározására a következő összefüggést alkalmazza:
ma K h a g Dk2,5 [kg/s],
(3.20)
ahol Kh
- a silógarat kúpszögétől, a belső és a falsúrlódási szögtől függő tényező (tapasztalati diagramokból kiolvasható) [22.].
A Zanker-módszer [23.] granulátumok, illetve szemcsés ömlesztett anyagok kiömlő mennyiségének gyors és egyszerű meghatározására szolgál, magkifolyás és kör alakú
- 36 -
szájnyílás esetén (Dk>100mm). A módszer egy grafikus nomogram segítségével könnyen alkalmazható, amely a következő összefüggésen alapul: ma
p Dkx
6, 288 t ge 23,16 d p 1,889 44,9
[kg/s],
(3.21)
ahol x
- a szemcsék alakjától függő kitevő (x=2,5…3,3) [23.].
A Gjacev-Keller módszer [24.] szemcsés, kohéziómentes ömlesztett anyagok instacionér áramlására vonatkozik és a következő összefüggéssel számítható: ma a A k
4 t g g Dk g t gh t k [kg/s], 4 t g Dk
(3.22)
ahol tk
- a kiömlés ideje [s].
A Janssen-módszer [17.] a magkifolyásra vonatkozó modell alapján (3/8. ábra), a Mohr-féle feszültségkörök felhasználásával, átlagos áramlási sebességgel számolva, a vertikális feszültség átlagos értékéből (3.11) határozza meg a kiáramló anyagmennyiséget, kohézív anyagok esetén: 2h Uk A C 2 A k ca cb k g ca c b ma A k a 1 e [kg/s]. Uk tge a
(3.23)
Kohéziómentes anyagok esetén (C=0) az összefüggés a következőképpen egyszerűsödik: 2h U k A A k ca cb k g ca cb [kg/s], ma A k a 1 e Uk
ahol
(3.24)
- az ömlesztett anyag tulajdonságaitól függő tényező (=0,5…0,65 jól ömleszthető anyagoknál) [17.].
A Johanson-módszer [25.] finom, kohézív ömlesztett anyagok esetén, egy folyamatosan összeomló anyaghídra felírt dinamikus egyensúlyból kiindulva határozza meg a kiáramló anyag mennyiségét. Az áramlás megindulásához szükséges az anyaghidat összetartó erő legyőzése, amely a folyási faktor függvénye. Tömegkifolyás és kúpos silógarat esetén a következő összefüggéssel számítható a kiáramló mennyiség:
D2 D g ff ma a k k 1 4 4 t g ff akt
1/ 2
[kg/s].
(3.25)
Kohéziómentes anyagok esetén (ff/ffakt=0): 1/ 2
D2 D g m a a k k 4 4 t g
ahol ff
- a folyási faktor, - 37 -
[kg/s],
(3.26)
ffakt
- a folyási faktor aktuális értéke.
A Tomas-módszer [26.] finom, kohézív ömlesztett anyagok esetén szintén egy folyamatosan összeomló anyaghídra felírt dinamikus egyensúlyon alapszik, de figyelembe veszi az instacionér áramlást is. Tömegkifolyás esetén a következő összefüggéssel számítható a kiáramló mennyiség: Dk k p d p 1 9 ma a A k (3.27) [kg/s], 2 m t 1 t g t sk f x st2 2 ahol mt - a garatgeometriára vonatkozó tényező [26.], tsk - a sebesség kiegyenlítődési idő [s], - a levegő dinamikus viszkozitása [Pa·s], xst - ekvivalens szemcseátmérő [mm], - az ömlesztett anyag porozitása. A Carleton-módszer [27.] az anyagszemcsék lineáris mozgási sebességének meghatározásán alapul, finom szemcsés anyagok esetén (dp<200m), a szemcsék, a szemcsék közötti levegő és a silószáj tulajdonságainak függvényében. A módszer előnye, hogy egy - a (3.28) egyenlet megoldásaként előállított - diagram segítségével, grafikus úton, egyszerűen meghatározható a szemcsék áramlási sebessége. Hátránya, hogy nem veszi figyelembe a falsúrlódást és a visszaáramló levegő ellenállását, illetve csak atmoszférikus kiáramlás mellett alkalmazható. A szemcsék sebessége a következő egyenletből számítható: 2/3 4 sin v02 15 1/lev3 lev v04 / 3 g, Dk p d5/p 3
(3.28)
ahol p lev lev v0
- a szemcsesűrűség [kg/m3], - a szemcsék közötti levegő sűrűsége [kg/m3], - a szemcsék közötti levegő viszkozitása [m2/s], - a kiömlő anyag sebessége a kifolyónyílásban [m/s].
Az Oldal-módszer [28.] kohéziómentes, szabad folyású, szemcsés ömlesztett anyagokra vonatkozik, magkifolyás és állandó anyagsebesség mellett: ma
ahol
2,5 2 g a Dk d p [kg/s], 6
(3.29)
- a szemcsék alaktényezője (búzára =0,4; PE granulátumra =0,3) [28.].
Az ismertetett kiömlési formulák mindegyike esetén található egy-egy anyagtípus, amelynél jó eredménnyel alkalmazható, de egyik módszer sem kínál univerzális, az ömlesztett anyagok többségénél alkalmazható számítási eljárást. Ennek következtében kutatásom fontos eleme, a bemutatott módszerek közül az adagolási modellben felhasználható összefüggés kiválasztása.
- 38 -
4. PNEUMATIKUS MUNKAHENGEREK MŰKÖDÉSI JELLEMZŐI Pneumatikus működtetésről akkor beszélünk, ha lineáris, vagy forgó mozgást sűrített levegő segítségével állítunk elő. Ezt úgy érhetjük el, ha egy dugattyú, vagy forgólapát két oldalának felületére adott nyomású levegőt juttatunk. A mozgást a felületek és a nyomások különbsége indukálja, amit a bevezetett tápnyomással (üzemi nyomással) lehet szabályozni (mivel a kisebb nyomású tér általában az atmoszférába ürül) [29.]. A pneumatikus végrehajtó elemek kialakításuk függvényében lehetnek: munkahengerek - egyoldali dugattyúrúddal, - átmenő dugattyúrúddal, - dugattyúrúd nélkül, forgatóhengerek, tömlőhengerek, stb. Működési módot tekintve megkülönböztetünk egyoldali és kétoldali működést. Legegyszerűbb megoldás az egyoldali működés, amelynél az alternáló mozgást egyik irányban sűrített levegő bevezetésével, másik irányban egy rúgó segítségével valósítjuk meg. Lényegesen elterjedtebb a kétoldali működés, ahol mindkét irányban sűrített levegő végzi a mozgatást. A kétoldali működésű végrehajtó elemek jobban szabályozhatóak és megbízhatóbbak. Kutatásaimban csak kétoldali működésű, lineáris, dugattyús pneumatikus munkahengerekkel foglalkoztam (4/1. ábra). Silóelzáró-szerkezetek működtetésénél ettől eltérő változatok is szóba jöhetnek (pl. forgatóhengerek), amelyekre az általános modell (a megfelelő összefüggések aktualizálásával) szintén alkalmazható. Hengerfedél
Hengerköpeny
Dugattyú
Dugattyúrúd
Levegő csatlakozás
4/1. ábra Kétoldali működésű munkahenger főbb részei A pneumatikus munkahengerek legfontosabb alkalmazási jellemzői: alkalmazhatók o erőkifejtésre, o adott pálya menti mozgatásra, egyszerű, gazdaságos és könnyen beépíthető építőelemek, általában egyenes-vonalú mozgatáshoz használják, könnyen szabályozható o a kifejthető erő, o az alkalmazott sebesség, nem lehet túlterhelni,
- 39 -
Tömítések
nehéz körülmények között is biztonságosan működik, a hengerátmérőt a maximális erőkifejtés határozza meg, a henger lökethosszát a mozgás maximális úthossza határozza meg, a legnagyobb megengedhető üzemi nyomás a konstrukció függvénye, stb.
Pneumatikus munkahengerek működési paramétereinek meghatározásánál a levegő tulajdonságaiból, illetve a levegőre érvényes hőtani [30.], illetve áramlástani [31.] összefüggésekből kell kiindulni. A pneumatika elméletével foglalkozó szakirodalmak többsége részletesen foglalkozik a csővezetékben történő áramlásokkal [32.], a be-, illetve kiáramló levegőmennyiség meghatározásával tartályok töltése, illetve ürítése esetén [33.]. A gyakorlati pneumatikával foglalkozó szakirodalom [29.] – igazodva ahhoz az alapelvhez, hogy a munkahenger mozgása csak a két véghelyzet előállítására szolgál – általában csak az alapegyenletekre szorítkozik. A munkahengerek kiválasztásához, méretezéséhez gyakorlati összefüggéseket alkalmaznak [34.]. A nemzetközi szakirodalomban régóta ismert a pneumatikus munkahengerek általános működési jellemzőit leíró differenciálegyenlet-rendszer (lásd. Blackburn [35.]), amely az összefüggések bonyolultsága miatt analitikusan nem oldható meg. A magyar szakirodalomban Elek és Hudáky [36.] írja le részletesen az általános és teljes működési modellt, illetve a működési jellemzők közötti kapcsolatot definiáló differenciálegyenlet-rendszert, melynek egyes elemeit a Magyar Szabvány [37.] is tartalmazza. Az Elek és Hudáky [36.] által használt általános módszer alapján külön fel kell írni a mechanikai, a hőtani és az áramlástani modellt.
4.1. Mechanikai modell Pneumatikus munkahengerek mechanikai modelljében a szerkezeti elemeket két részre kell bontani. Az egyik részt (hengerköpeny, hengerfedelek, henger-rögzítés, stb.) a működés során álló helyzetűnek tekintjük, a másik rész (dugattyú, dugattyúrúd, teher-rögzítés, stb.) végzi a mozgást (4/2. ábra). Mínusz kamra
Plusz kamra
Plusz mozgás Fs cd
m
p+ A+ sd
A-
p-
kd Fd Mínusz mozgás
4/2. ábra Munkahengerek mechanikai modellje [32.] A működés szempontjából csak a mozgó részre ható erőknek van jelentősége, amelyek a következők: a dugattyúra ható erők: - a plusz kamrában, - a mínusz kamrában,
- 40 -
a belső súrlódás (a dugattyútömítés és a hengerköpeny között, illetve a dugattyúrúd tömítésénél), a külső súrlódás, a külső terhelés. Plusz kamra alatt azt a hengerteret értjük, amelyet a dugattyú plusz irányú mozgása (a dugattyúrúd kitolódik a munkahengerből) során töltünk (4/2. ábra). Mínusz kamra alatt azt a hengerteret értjük, amelyet a dugattyú mínusz irányú mozgása (a dugattyúrúd visszahúzódik a munkahengerbe) során töltünk (4/2. ábra). A mínusz kamra hasznos térfogata mindig kisebb, mint a plusz kamráé, mert tartalmazza dugattyúrudat is. Ennek következtében a munkahenger plusz mozgási irányban mindig nagyobb erőt tud kifejteni. A felvett mechanikai modellnél az összes fellépő erőhatás iránya megegyezik (rúdirányú), a keresztirányú terheléseket a rendszer szempontjából elhanyagolhatónak tekintjük (a munkahenger csuklós rögzítésével általában kiküszöbölhetők). A mozgó elemekre felírt impulzustétel alapján:
p t A p t A k d sd (t) cd sd (t) Fs Fd t m sd (t) ,
(4.1)
ahol p+ pA+ Am kd cd Fs Fd
- a nyomás a plusz kamrában [Pa], - a nyomás a mínusz kamrában [Pa], - a dugattyú felülete a plusz kamrában [m2], - a dugattyú felülete a mínusz kamrában [m2], - a munkahenger által mozgatott összes tömeg [kg], - a külső csillapítási tényező [Ns/m], - a külső rugalmassági tényező [N/m], - a dugattyúra ható összes súrlódó erő [N], - a dugattyúra ható összes külső terhelőerő [N].
Az adagolási folyamat szempontjából a szerkezet rugalmasságából, illetve a külső csillapításból adódó erőhatások nem jelentősek, illetve a külső terhelőerőben jelennek meg, ezért a vizsgálatok során elhanyagoljuk (kd=0, cd=0). Az elhanyagolások után a (4.1) egyenlet a következőképpen módosul:
p t A p t A Fs Fd t m sd (t) .
(4.2)
4.2. Hőtani modell Pneumatikus munkahengerek működése során a rendszer állapotváltozását egy nyitott rendszer állapotváltozásaként lehet értelmezni, ahol a folyamat kezdetétől a fél ciklus végéig a rendszerbe be-, illetve onnan kiáramló levegőmennyiség okozza a rendszer állapotváltozásait. Mivel a folyamat során a dt időtartam alatt be-, illetve kiáramló levegő mennyisége nem egyezik meg, ezért az állapotváltozás végig instacionér. A 4/3. ábra alapján a rendszerre felírható energiaegyenlet:
dQ cp Tb dmb
v2b v2 dmb dU dW cp Tk dmk k dmk , 2 2
ahol dU dQ dW
- a rendszer belső energia változása [J], - a rendszerbe bevitt hőmennyiség [J], - a rendszer által végzett munka [J], - 41 -
(4.3)
cp Tb Tk pb pk dmb dmk vb vk
- izobár fajhő [J/(kg∙K)], - a beáramló levegő hőmérséklete [K], - a kiáramló levegő hőmérséklete [K], - a beáramló levegő nyomása [Pa], - a kiáramló levegő nyomása [Pa], - az egységnyi idő alatt beáramló levegő mennyisége [kg/s], - az egységnyi idő alatt kiáramló levegő mennyisége [kg/s], - a beáramló levegő sebessége [m/s], - a kiáramló levegő sebessége [m/s].
dmb
dmk
Tb
Tk
dU vb
vk
dQ
dW
4/3. ábra A munkahenger hőtani modellje Pneumatikus munkahengerek esetén a két kamra a közeg szempontjából teljesen elkülönül egymástól (4/4. ábra), az egyes tereket a működés során csak töltjük, vagy ürítjük. Mínusz kamra
Plusz kamra
dugattyú mozgás (plusz irány)
A+ V+ p+
A-
V-
pdmk
dmb
4/4. ábra A munkahenger hengertereinek hőtani elkülönülése Elek és Hudáky [36.] nyomán, a (4.3) egyenlet, az általános gáztörvény, a gázok jellemzői közötti összefüggések, valamint a belső energia és a munka meghatározására vonatkozó összefüggések felhasználásával (ha dQ=0), a következő alakban írható fel az egyes kamrák esetén: ha a hengertérbe csak befelé áramlik a levegő:
v2b R Tb mb 1 mb p V p V , 2
(4.4)
ha a hengertérből csak kifelé áramlik a levegő:
R Tk mk 1
v2k mk p V p V , 2
ahol p
- nyomás az adott kamrában [Pa], - 42 -
(4.5)
R V
- az egyetemes gázállandó [J/(kg∙K)], - az adiabatikus kitevő, - az adott kamra térfogata [m3].
4.3. Áramlástani modell Elek és Hudáky [36.] nyomán, az áramló levegőmennyiség meghatározásához fel kell írni egy, a terek közötti levegőáramlásra vonatkozó általános modellt is (4/5. ábra). A
B
pA > pB pA A
vA
vB
pB B
4/5. ábra Terek közötti áramlás modellje A 4/5. ábrán látható esetre felírva az összenyomható gázokra érvényes Bernoulli-egyenletet [31.], egydimenziós, súrlódásmentes, örvénymentes és stacioner áramlás esetén, ha a levegő tömegét elhanyagoljuk (figyelembe véve, hogy pA>pB és vA=0): pA 1 2 pB vB , 1 A 2 1 B
(4.6)
ahol pA, pB - a nyomás az egyes terekben [Pa], A, B - a levegő sűrűsége az egyes terekben [kg/m3], vB - a B térbe beáramló levegő sebessége [m/s]. A (4.6) összefüggés felhasználásával az áramlási sebességre, hangsebesség alatti áramlás esetén a következő összefüggés adódik: 1 2 pA pB . vB 1 1 A p A
(4.7)
Amennyiben az áramlás sebessége eléri a helyi hangsebesség értékét, akkor az összefüggés a következőképpen módosul: 1 pA p vkrit , A
(4.8)
2 0,52828 . 1
(4.9)
vBhf ahol
p vkrit
1
Fentiek alapján az áramló levegő mennyisége a következő összefüggés szerint számítható: m f Af B vB ,
ahol f – a fojtási tényező, - 43 -
(4.10)
Af – a beömlési keresztmetszet [m2]. A fenti összefüggéseket felhasználva aktualizálhatjuk a munkahenger kamráiba be-, illetve onnan kiáramló levegő mennyiségét és sebességét. A munkahenger kamráinak töltése során (4/4. ábra) az A tér nyomását (működtető nyomás) állandónak tekintjük (pA=pt), a B tér nyomása pedig változik a kiinduló nyomás és a működtető nyomás értéke között (pB=p). Ehhez azt kell feltételeznünk, hogy a működtető nyomás egy megfelelően nagy méretű légtartályon keresztül jut a munkahengerbe (ebben az esetben a kiáramló levegő mennyisége nem okoz nyomás és sűrűségváltozást az állandó térfogatú tartályban). A (4.7) és (4.10) egyenletek felhasználásával, a beáramló levegő mennyisége a következő összefüggéssel számítható [36.]:
m b f Af p t
b , R Tt
(4.11)
ahol pt Tt Ψb
- a működtető nyomás [Pa], - a működtető levegő hőmérséklete [K], - a beáramlásra vonatkozó átömlési függvény.
A beáramlásra vonatkozó átömlési függvényt leíró összefüggés eltér hangsebesség alatti áramlás során: 2 1 2 p p p , ha ba 1 p vkrit , pt 1 pt pt
hangsebességgel történő áramlás során: 2 1 2 p p vkrit 1 p vkrit , ha p vkrit . 1 pt
bh
(4.13)
1
Átömlési függvény
(4.12)
0.8
a
0.6 0.4 0.2 0 5 1 10
2 10
5
5 5 3 10 4 10 Nyomás [P a]
5 10
5
6 10
4/6. ábra Beáramlásra vonatkozó átömlési függvény (3mm csőátmérő, Pt=6∙105Pa és Tt=293,14K esetén)
- 44 -
5
A munkahenger kamráinak ürítése alatt (4/4. ábra), az A tér nyomása változik az atmoszférikus nyomás és a működtető nyomás értéke között (pA=p), a B tér nyomása pedig az atmoszférikus nyomás értékén állandó lesz (pB=p0). Ehhez az szükséges, hogy a kamra a szabadba ürüljön ki. A (4.7) és (4.10) egyenletek felhasználásával, a kamrából az atmoszférába ürülő levegő mennyisége a következő összefüggéssel számítható [36.]:
m k f A f p
k , R T
(4.14)
ahol p0 T Ψk
- az atmoszférikus nyomás [Pa], - a levegő hőmérséklete az adott kamrában [K], - a kiáramlásra vonatkozó átömlési függvény.
Az átömlési függvényt leíró összefüggés itt is eltér hangsebesség alatti áramlás során: 1 p 2 p0 p0 , ha 0 p vkrit , ka 1 1 p p p 2
(4.15)
hangsebességgel történő áramlás során: 2 1 p 2 p vkrit 1 p vkrit , ha 0 p vkrit . 1 p
kh
(4.16)
Átömlési függvény
1 0.8 0.6
a
0.4 0.2 0 5 6 10
5 10
5
5 5 4 10 3 10 Nyomás [P a]
2 10
5
1 10
5
4/7. ábra Kiáramlásra vonatkozó átömlési függvény (3mm csőátmérő, Pt=6∙105Pa és Tt=293,14K esetén)
4.4. Pneumatikus munkahengerek működését leíró egyenletrendszer Összefoglalva a 4.1.-4.3. fejezetben definiált modellek összefüggéseit, megállapíthatjuk, hogy a pneumatikus munkahengerek működése csak akkor írható le matematikailag egzakt módon, ha a differenciálegyenlet-rendszer teljes egészében megoldásra kerül. A munkahenger teljes működési ciklusa egy plusz és egy mínusz mozgást tartalmaz. Mivel a két mozgás teljesen elkülönül egymástól, ezért célszerű ezeket külön vizsgálni. A - 45 -
differenciálegyenletre vonatkozóan ez csak annyit jelent, hogy a plusz és a mínusz kamra jellemzői helyet cserélnek az összefüggésekben. A pneumatikus munkahengerek működési jellemzőit leíró differenciálegyenlet-rendszer elemei (plusz mozgásra):
p t A p t A Fs Fd t m sd (t) .
(4.17)
Plusz kamra töltése:
v2b t mb t R T 1 p t sd t p t sd t , b 2 A p ha p vkrit : pt 2 1 2 p t p t , mba t f Af p t t 1 R Tt 1 pt pt
v ba t
ha
1 2 pt p t , 1 1 t p t
(4.18)
(4.19)
(4.20)
p p vkrit : pt
mbh f Af p t R Tt
vbh
1
2 2 1 const. , 1
pt 2 const. t 1
(4.21)
(4.22)
Mínusz kamra ürítése:
v2k t mk R T 1 p t sd t p t sd t , k 2 A p0 ha p vkrit : p 2 1 p0 p0 2 mka t f Af p t 1 , R T 1 p t p t
v ka t
1 2 p t p0 , 1 1 t p t
- 46 -
(4.23)
(4.24)
(4.25)
ha
p0 p vkrit : p mkh t f Af p t R T
vkh t
p t
t
1
2 2 1 , 1
2 . 1
(4.26)
(4.27)
Mivel a fenti egyenlet-rendszer analitikusan nem oldható meg, a nemzetközi szakirodalomban számos szerző foglalkozik numerikus megoldásával, illetve adott célokra történő alkalmazásával. Az újabb keletű kutatások közül kiemelkednek a munkahengerek két véghelyzet közötti megállításával, pozícionálásával [38.], illetve a tervező szoftverekben történő felhasználással foglalkozó [39.] publikációk. A mozgási fázisokra vonatkozó megoldási sémát ismertet Sorli et al. [40.] lineáris és forgató munkahengerekre vonatkozóan, de az általa felírt modell nem ad információt a mozgásmentes fázisokra vonatkozóan. Az adagolási folyamat szempontjából alapvető fontosságú a munkahenger teljes ciklusának vizsgálata, amelyhez definiálni kell a működés fázisait és az egyes fázisokban érvényes összefüggéseket. Mivel a fenti megoldások egyike sem ad az adagolási folyamat során alkalmazható modellt a pneumatikus munkahenger működésének leírására, ezért kutatásaim egyik feladata a működést leíró általános modell kidolgozása.
- 47 -
5. PNEUMATIKUS MODELLJE
MUNKAHENGEREK
ÁLTALÁNOS
MŰKÖDÉSI
A 4. fejezetben bemutatott differenciálegyenlet-rendszer, az összefüggések bonyolultsága miatt, csak numerikus módszerek felhasználásával oldható meg. A numerikus megoldáshoz olyan modellt kell felírni, amivel a vizsgált folyamat működése leírható.
5.1. A működés fázisai Kétoldali működésű, lineáris dugattyús, véghelyzet fékezésű pneumatikus munkahengerek működési fázisai: 1. Kezdeti töltési-ürítési állapot 2. A dugattyú normál mozgása: a. gyorsuló mozgás b. állandó sebességű mozgás 3. A dugattyú csillapított mozgása 4. Végső töltési-ürítési állapot Az 1. fázis alatt a dugattyú nem mozog. A munkahenger egyik kamrájába a tápnyomás áramlik, a másik az atmoszférába ürül (plusz mozgásirány esetén a plusz kamra töltése és a mínusz kamra ürítése történik). A 2. fázis elején a dugattyú gyorsulni kezd (2.a. fázis), majd állandósult sebességet vesz fel (2.b. fázis). A dugattyúsebesség a munkahenger-geometriától és a terheléstől függően bizonyos esetekben nem állandósul. A 3. fázisban a munkahenger, a löketvégi fékezésnél alkalmazott, csökkentett kiáramlási keresztmetszet miatt lassulni kezd, majd megáll. A 4. fázis alatt a dugattyú szintén nem mozog. A munkahenger kamráiban a nyomások kiegyenlítődnek (plusz mozgásirány esetén a plusz kamra a működtető tápnyomás értékét, a mínusz kamra az atmoszférikus nyomás értékét veszi fel). Az egyes fázisok az alkalmazható összefüggések szempontjából három csoportba sorolhatók: állandó tér feltöltése, illetve ürítése (1. és 4. fázis), ahol dV 0,sd 0,sd 0 ,
állandó sebességű térfogatváltozás (2.b. fázis), ahol dV const.,sd 0 ,
(5.1)
(5.2)
változó sebességű térfogatváltozás (2.a. és 3. fázis), ahol dV 0,sd 0 .
(5.3)
A (5.1) és (5.2) feltételeknek megfelelően a (4.17) egyenlet az 1., 2.b. és 4. fázisban a következő formára egyszerűsödik:
p t A p t A Fs Fd t 0 .
(5.4)
A (5.1) feltételhez igazodva a (4.18) és (4.23) egyenletek az 1. és 4. fázisban a következő módon egyszerűsödnek: v2b t R Tb 1 mb t p t V0 , 2 - 48 -
(5.5)
v2k t R Tk 1 mk t p t V0 . 2
(5.6)
ahol V+0 V-0
- holttér a plusz kamrában [m3], - holttér a mínusz kamrában [m3].
A pneumatikus munkahenger működése során a 2.a. fázis nem minden esetben különíthető el a 2.b. fázistól, de ez nem jelent problémát, mert az összefüggések mindkét fázisban azonosak (a 2.b. fázisban a gyorsulás értéke zérus, de ez nem befolyásolja a megoldhatóságot).
5.2. Működési modell A munkahenger mozgási fázisainak és működési jellemzőinek ismeretében már definiálható a kétoldali működésű, lineáris dugattyús, véghelyzet fékezésű pneumatikus munkahengerek működési folyamatának elvi modellje (5/1. ábra). Folyamat indítása
1. MŰKÖDÉSI FÁZIS
Kezdeti töltési-ürítési állapot
nem teljesül
Indulási feltétel
Bemenő paraméterek Méretek, kiinduló nyomások (pt)
p+ és p- nyomás változik
Plusz mozgás
Mínusz mozgás
p t A p t A Fs Fd t
p t A p t A Fs Fd t
teljesül
2. MŰKÖDÉSI FÁZIS
Gyorsuló és állandó sebességű dugattyú mozgás
nem teljesül
Csillapítási feltétel
p+, p- nyomás, és az s elmozdulás változik
Plusz mozgás
Mínusz mozgás
sd t sf
sd t L sf
teljesül
3. MŰKÖDÉSI FÁZIS Csillapított dugattyú mozgás
nem teljesül
Megállási feltétel
p+, p- nyomás, és az s elmozdulás változik
Plusz mozgás
Mínusz mozgás
sd t 0
sd t L
teljesül
4. MŰKÖDÉSI FÁZIS Végső töltési-ürítési állapot
nem teljesül
Kiegyenlítődési feltétel
p+ és p- nyomás változik
Plusz mozgás
p t pt
Mínusz mozgás
p t p0
teljesül Folyamat leállítása
Kimenő jellemzők p+(t) és p-(t) nyomás, s(t) elmozdulás
5/1. ábra A működési modell általános sémája - 49 -
p t p0
p t pt
A differenciálegyenlet-rendszer összefüggései az egyes fázisokban eltérő alakot vesznek fel, amelyek megadásával rendelkezésre áll a kétoldali működésű, lineáris dugattyús, véghelyzet fékezésű pneumatikus munkahengerek működési folyamatát leíró modell. A modell alapján a folyamat paraméterei numerikus módszerek felhasználásával meghatározhatók. Az összefüggések azonosak az 1. és 4. fázisban (5/2. ábra), a különbség a mínusz kamra levegő kiáramlási keresztmetszeténél és aktív dugattyúfelületénél adódik: az 1. fázisban Af- és A- (plusz mozgásirány), Af- és A- (mínusz mozgásirány), a 4. fázisban Af-=Acf- és A-=Afh- (plusz mozgásirány), Af+=Acf+ és A+=Afh+ (mínusz mozgásirány). A mozgásegyenletnek az 5/2. ábrán szereplő alakja nem befolyásolja a nyomások alakulását, csak az indulási feltétel meghatározásánál játszik szerepet. A) Mozgásegyenlet p t A p t A Fs Fd t 0
B) Plusz kamra (töltés) v 2b t R Tb 1 m b t p t V0 2
p p vkrit pt
igen
nem
C) Mínusz kamra (ürítés)
v2k t R Tk 1 m k t p t V0 2 m k t m ka t
igen
v k t v ka t
nem m t m t p0 k kh p vkrit p v k t v kh t
m b t m ba t
m b t m bh t
v b t v ba t
v b t v bh t
Kimenő paraméterek: p+(t) és p-(t)
5/2. ábra Az 1. és 4. működési fázis esetén érvényes összefüggések (plusz mozgás) A) Mozgásegyenlet p t A p t A Fs Fd t m sd (t)
B) Plusz kamra (töltés)
v2b t mb t V p t sd t p t 0 sd t R Tb 1 2 A A
p p vkrit pt
igen
nem
C) Mínusz kamra (ürítés)
v2k t mk t V p t sd t p t 0 sd t R Tk 1 2 A A m k t m ka t
v k t v ka t
igen
nem m t m t p0 k kh p vkrit p v k t v kh t
m b t m ba t
m b t m bh t
v b t v ba t
v b t v bh t
Kimenő paraméterek: p+(t), p-(t) és sd(t)
5/3. ábra A 2. és 3. működési fázis esetén érvényes összefüggések (plusz mozgás) - 50 -
Hasonlóan azonos összefüggések írhatók fel a 2. és 3. fázis esetén (5/3. ábra), a különbség egyrészt itt is a mínusz kamra levegő kiáramlási keresztmetszeténél és aktív dugattyúfelületénél adódik: a 2. fázisban Af- és A- (plusz mozgásirány), Af+ és A+ (mínusz mozgásirány), a 3. fázisban Af-=Acf- és A-=Afh- (plusz mozgásirány), Af+=Acf+ és A+=Afh+ (mínusz mozgásirány). További eltérés, hogy a gyorsulás értéke a 2.a. fázis esetén nem zérus, a 2.b. fázisban pedig zérus értéket vesz fel. Mínusz mozgásirány esetén az összefüggések nem változnak a plusz mozgásirányhoz képest, csak a nyomások helye felcserélődik (a plusz kamrát ürítjük, a mínusz kamrát töltjük).
5.3. A differenciálegyenlet-rendszer megoldása Bizonyos esetekben – elsősorban a bonyolultabb összefüggéseknél – a differenciálegyenletek megoldására a feladatot diszkretizálni kell, ami azt jelenti, hogy nem keresünk közelítő képletet a megoldás számára, hanem a megoldás közelítő numerikus értékeit határozzuk meg valamilyen kezdeti értékből kiindulva, olyan közelítéssel, amelynek hibája lehetőleg kicsi. A diszkretizált feladatok megoldására rekurzív, vagy iteráló képleteket használunk. A módszer lényege, hogy egy adott pontban a keresendő függvény értékét korábban már meghatározott értékek alapján számítjuk ki. Az eredmények pontosságának fokozása érdekében az értékeket folyamatosan pontosítva ismételgetjük, feltételezve, hogy a folyamat konvergens. A numerikus értékek közötti távolság (ún. alaptávolság) megválasztása kritikus pontja az iteratív módszereknek, mert a megfelelő pontosság eléréséhez kis távolság, és ebből adódóan magas ismétlési szám szükséges, de az ismétlési szám növekedésével növekszik a közelítő értékek hibája is (hibaöröklés miatt) [41.]. Differenciálegyenlet-rendszerek megoldására számos numerikus módszer létezik (pl. Taylormódszer, Runge-Kutta módszer, stb.). Az hogy melyik módszert alkalmazzuk egy adott feladat megoldására, a feladat jellegétől, illetve a rendelkezésre álló peremfeltételektől függ. Az általam felírt modell vizsgálatánál a differenciálegyenlet-rendszer megoldására a RungeKutta módszert választottam. A választásom oka kettős: a Runge-Kutta módszer az ürítési folyamat minden fázisában alkalmazható az aktualizált egyenlet-rendszer megoldására, az alkalmazott szoftver (MathCAD) ezt a módszert alkalmazza a differenciálegyenletrendszerek megoldására. A pneumatikus munkahengerek leírására használt (4.17)-(4.27) általános differenciálegyenletrendszer (5.3) numerikus megoldáshoz szükséges alakja a következő (ha a kiáramlás során a hőmérséklet változásától eltekintünk): p t f1 p t ,sd t ,sd t ,
(5.7)
p t f 2 p t ,sd t ,sd t ,
(5.8)
sd t f3 p t , p t ,sd t .
(5.9)
A numerikus megoldáshoz szükséges a munkahenger mozgási fázisaihoz történő aktualizálás, ami minden alkalmazott fázis esetén külön peremfeltételek definiálását igényli. - 51 -
Az általános alak a következő módon egyszerűsödik állandó sebességű mozgás (5.2): p t f1 p t ,sd t ,
(5.10)
p t f 2 p t ,sd t ,
(5.11)
illetve állandó tér feltöltés, illetve ürítés (5.1) esetén: p t f1 p t ,
(5.12)
p t f 2 p t .
(5.13)
Az differenciálegyenlet-rendszer numerikus megoldásához a MathCAD 2001i Professional szoftvert használtam. A MathCAD egy olyan Microsoft Windows operációs rendszerre telepíthető matematikai szoftver, amely egyszerű, lineáris menetű feladatmegoldást tesz lehetővé egy egyszerű Windows felületen. A feladatok definiálása és megoldása során egy A/4 (vagy más szabványos méretű) formátumú lapon dolgozva kell egymás után definiálni a feladat változóit, állandóit és a közöttük fennálló összefüggéseket. A szoftver beépített függvényei segítségével számos matematikai feladat megoldható, a legegyszerűbb egyenletektől a differenciálegyenletrendszerekig. A kapott eredmények ábrázolása is egyszerűen megvalósítható két- és háromdimenziós, különböző formájú grafikonokon. A szoftver számos más matematikai és grafikai, illetve animálási lehetőséget is tartalmaz, de értekezésemben csak három képességét aknáztam ki: differenciálegyenlet-rendszerek megoldása, kétdimenziós ábrázolás, mérési eredmények kiértékelése és ábrázolása. Az 5/4. ábra szemlélteti a változók és függvények megadásának tényleges megvalósítását a MathCAD szoftver kezelőfelületén.
5/4. ábra Változók és függvények definiálása a MathCAD szoftverrel - 52 -
A differenciálegyenlet-rendszer megoldása után a program által előállított jellemzők: nyomás: - a plusz kamrában, - a mínusz kamrában, dugattyúmozgási jellemzők: - elmozdulás, - sebesség, - gyorsulás.
Nyomások [Pa]
Az 5/5. ábrán látható diagramok pt=6∙105Pa, Dh=25mm, L=0,1m, Ft=19N főbb munkahenger paraméterek esetén mutatnak be egy példát (a dugattyú plusz mozgási irányában) a program működésére. A szoftver segítségével egy diagramban is lehet ábrázolni a teljes ciklust, de ez a munkahenger működésére vonatkozóan nem jelent többletinformációt (az adagolási folyamat szempontjából a teljes ciklus kerül vizsgálatra). 6 10
5
5 10
5
4 10
5
3 10
5
2 10
5
1 10
5
0
0
0.2
0.4 0.6 Idõ [s]
0.8
1
Nyomás a plusz kamrában Nyomás a mínusz kamrában
a) nyomások a hengerterekben
b) dugattyúelmozdulás
c) dugattyúsebesség
d) dugattyúgyorsulás
5/5. ábra A dugattyú plusz irányú mozgásának paraméterei Az 5/5. ábrán látható diagramokon jól megfigyelhetők az egyes működési fázisok: a mozgás nélküli feltöltési és ürítési, a közel állandó nyomáskülönbség által hajtott mozgó fázis, illetve a mozgásnak a levegő összenyomhatóságából adódó dinamikája. - 53 -
5.4. A modell eredményeinek ellenőrzése Az általános modell és az azt megvalósító program eredményeinek ellenőrzésére konkrét méréseket végeztem. 5.4.1. A vizsgált munkahenger A pneumatikus munkahenger paramétereinek kiemelt fontossága van a működés modellezése szempontjából, hiszen ezekből kiindulva határozhatók meg a számítógépes program bemenő adatai. A vizsgálatok, illetve a mérések során egy Mecman 166-55-0210-1 típusú kétoldali működésű, lineáris dugattyús, véghelyzet fékezésű pneumatikus munkahengert alkalmaztam [34.].
Dh
dhf
dh
sf
dc+
L dcf+
dcdcf-
5/6. ábra A vizsgálatoknál alkalmazott pneumatikus munkahenger fontosabb geometriai méretei A vizsgálatok elvégzéséhez szükséges geometriai méretek a következők (5/6. ábra): Dh - a dugattyú átmérője [m], dh - a dugattyúrúd átmérője [m], dhf - a fékeződugattyú átmérője [m], dc1, - a plusz kamra levegőcsatlakozásának átmérője [m], dc2, - a mínusz kamra levegőcsatlakozásának átmérője [m], dcf1, - kiáramlási átmérő fékezéskor a plusz kamránál [m], dcf2, - kiáramlási átmérő fékezéskor a mínusz kamránál [m], sf - a fékezési úthossz [m], L - a maximális lökethossz [m]. A fenti geometriai adatokból számíthatók ki azok a szükséges jellemzők, amelyek a szoftver alkalmazása közben kerülnek felhasználásra. A munkahenger kamráinak térfogata, a dugattyúfelületek és az áramlási keresztmetszetek a dugattyú mozgása során fázisról-fázisra változnak, ezért az egyes mozgásszakaszoknál meg kell határozni az aktuális értékeket. Geometriai szempontból a dugattyú mozgásának szakaszokra bontása és az egyes szakaszokban figyelembevett értékek az 5/1. táblázatban láthatók, ahol sd a dugattyú elmozdulásának pillanatnyi értéke. A 5/1. táblázatban szereplő jellemzők és leíró összefüggéseik a következők: dugattyúfelületek: D2 (5.14) A h , 4
- 54 -
A
D
2 h
d 2h 4
(5.15)
,
dugattyúfelületek fékezéskor: Afh Afh
D
2 h
d hf2 4
(5.16)
,
ahol Vf+0 Vf-0
- holttér a plusz kamrában fékezéskor [m3], - holttér a mínusz kamrában fékezéskor [m3]. 5/1. táblázat A dugattyú mozgási szakaszaiban érvényes geometriai jellemzők
a plusz mozgás során:
PLUSZ KAMRA aktuális térfogat áramlási (V+) keresztm. (beáramlás)
MÍNUSZ KAMRA aktuális térfogat áramlási (V-) keresztm. (kiáramlás)
+ terek feltöltése, ürítése
V+0
Asf+
V-max=V-0+A-L
As-
+ csillapítatlan mozgás
V+0+A+sd
Asf+
V-max-A-sd
As-
+ fékezés a 2. hengertérben
V+0+A+sd
Asf+
Vf-0+Afh-(L-sd)
Acf-
+ terek feltöltése, ürítése
V+max=V+0+A+L
Asf+
Vf-0
Acf-
a mínusz mozgás során
(kiáramlás)
+ terek feltöltése, ürítése
V+max
As+
V-0
Asf-
+ csillapítatlan mozgás
V+max-A+sd
As+
V-0+A-sd
Asf-
+ fékezés az 1. hengertérben
Vf+0+Afh+(L-sd)
Acf+
V-0+A-sd
Asf-
+ terek feltöltése, ürítése
Vf+0
Acf+
V-max
Asf-
Mozgási szakasz
(beáramlás)
A munkahenger holtterei adott állandó értékek, melyek katalógusból [34.] kiolvashatók, illetve manuálisan lemérhetők. A tényleges áramlási keresztmetszetek a mérési folyamat összeállítása (lásd 3.3.2. fejezet) miatt nem a munkahenger csatlakozó méreteivel számíthatók, hanem az alkalmazott fojtóvisszacsapó szelepek határozzák meg és katalógusból felvehetők [42.]. A löketvégi fékezéskor alkalmazott áramlási keresztmetszet a munkahenger jellemzője, de állítható és a katalógus alapján határozható meg [34.]. Az alkalmazott áramlási keresztmetszetek: As+ Asf+ AsAsfAcf+ Acf-
- a fojtószelep átömlési keresztmetszete nyitott irányban (plusz kamra) [m2], - a fojtószelep átömlési keresztmetszete fojtott irányban (plusz kamra) [m2], - a fojtószelep átömlési keresztmetszete nyitott irányban (mínusz kamra) [m2], - a fojtószelep átömlési keresztmetszete fojtott irányban (mínusz kamra) [m2], - az átömlési keresztmetszet fékezéskor a plusz kamra esetén [m2], - az átömlési keresztmetszet fékezéskor a mínusz kamra esetén [m2].
- 55 -
5.4.2. A mérés elrendezési vázlata és elemei A modell és az alkalmazott program eredményeinek ellenőrzéséhez egy olyan rendszert kellett kialakítani, amelyben lehetőség van a működési folyamat minden szükséges paraméterének mérésére (5/7. ábra). A mért paraméterek: nyomás a plusz kamrában, nyomás a mínusz kamrában, a dugattyú elmozdulása. Munkahenger
Útmérő
2. Nyomásmérő
1. Nyomásmérő
A
B
C
D
E
F
G
H SELECTED ON-LINE
1. Fojtóvisszacsapó szelep
2. Fojtóvisszacsapó szelep
Mérési adatgyűjtő
Főszelep
Mérés feldolgozás
Kompresszor
5/7. ábra A mérési folyamat kialakítása A mérés megvalósításához a következő elemeket kellett - a vizsgált munkahenger mellett - a mérési körbe beépíteni: energiaforrás, vezérlő elem, szabályozó elemek, mérőelemek, mérési adatgyűjtő és feldolgozó elemek. A mérésnél alkalmazott energiaforrás egy ABAC 100 HP3 típusú kompresszor, kimeneti nyomástartomány: 0,1-1,2MPa (túlnyomás). A vezérlő elem egy MECMAN MOD E TYPE 463/300 típusú 5/2-es, bistabil, elektromágneses vezérlésű szelep, kézi rásegítéssel. A szabályozó elemek: 1. FESTO D:S-PPV-GR típusú fojtó-visszacsapó szelep, nyomástartomány: 0,03-1MPa (plusz kamra), 2. FESTO D:S-PPV-GRLA típusú fojtó-visszacsapó szelep, nyomástartomány: 0,03-1MPa (mínusz kamra). - 56 -
Az alkalmazott mérőeszközök: 1. HBM P6A típusú nyomásmérő, maximális nyomás: 10MPa, mérési érzékenység: 2mV/V (plusz kamra), 2. HBM P6A típusú nyomásmérő, maximális nyomás: 10MPa, mérési érzékenység: 2mV/V (mínusz kamra), 3. BANNER Q50BVI típusú optikai távadó, mérési tartomány: 0,1-0,3m, legkisebb válaszidő: 0,004 s. A mérési adatgyűjtő egy HBM SPIDER 8 típusú berendezés. A mérési adatok feldolgozása egy ASUS típusú notebook-on történt, Catman 3.1 típusú szoftverrel, az adatokat txt formátumban nyertem ki. 5.4.3. A mérés eredményei Az 5/7. ábrán bemutatott elrendezésben lefolytatott méréseknél a következő beállításokat alkalmaztam: a munkahenger fő paraméterei: o lökethossz: L = 0,1 m, o dugattyúátmérő: Dh = 0,025 m, o a dugattyúrúd átmérője: dh = 0,012 m, o a fékeződugattyú átmérője: dhf = 0,018 m, o a fékezési úthossz: sf = 0,011 m, o áramlási keresztmetszetek: As+ = As- = 2,1 ∙ 10-6 m2, Asf+ = 4,8 ∙ 10-7 m2, Asf- = 3,2 ∙ 10-7 m2, Acf+ = Acf- = 6,3 ∙ 10-7 m2, a dugattyút terhelő erő (állandó): Fd = 20 N, a mozgatott tömeg: m = 0,38 kg, működtető nyomás (abszolút nyomás): pt = 5∙105 - 7∙105 Pa.
a) plusz mozgás
b) mínusz mozgás
Nyomás a plusz kamrában Nyomás a mínusz kamrában
5/8. ábra A nyomások alakulása (pt=5∙105Pa) A vizsgálatok során a működtető nyomást változtattam, 3 diszkrét beállításban (5∙105Pa; 6∙105Pa; 7∙105Pa). - 57 -
A mérési összeállítás nem tette lehetővé a dugattyúra ható erő változtatását, ezért ennek hatását nem vizsgáltam, a terhelést egy állandó terhelőerővel vettem figyelembe. Minden beállításnál 5 azonos mérést végeztem az esetleges működési anomáliák kizárására és ez alapján választottam ki egy jellemző érték-sorozatot. A mérési adatok megjelenítése MathCAD 2001i Professional szoftverrel történt, Microsoft Excel 2003 xls formátumba konvertálás után. Az xls formátumú adatok közvetlenül importálhatók a MathCAD 2001i Professional szoftverbe. A mérési eredmények az M1. mellékletben találhatók, a 5/8. és 5/9. ábrák szemléltetik a nyomások és a dugattyúelmozdulás alakulását egy mérési beállítás esetén. A diagramokon jól megfigyelhetők a munkahenger mozgási fázisai.
a) plusz mozgás
b) mínusz mozgás
5/9. ábra A dugattyú elmozdulása (pt=5∙105Pa) 5.4.4. A mért és számított eredmények összehasonlítása A programban a mérési rendszer jellemzői mellett a következő adatok kerültek felhasználásra: levegő jellemzők: o az atmoszférikus nyomás: p0 = 105 Pa, o az atmoszférikus levegő sűrűsége: 0 = 1,19 kg/m3, o hőmérsékletek: T0 = Tt = T+ = T- = 293 K, o az egyetemes gázállandó: R = 287,14 J/(kg∙K), o az adiabatikus kitevő: = 1,4, o nehézségi gyorsulás: g = 9,81 m/s2, munkahenger paraméterek: o súrlódási tényező: h = 0,1, o holtterek: V+0 = 8,8 ∙ 10-6 m3, V-0 = 9,1 ∙ 10-6 m3, Vf+0 = Vf-0 = 1,24 ∙ 10-6 m3, o fojtási tényezők: s+ = s- = cf+ = cf- = 0,61, sf+ = sf- = 0,51. A súrlódási ellenállás értéke különböző lehet az egyes tömítésektől függően, de hasonló tömítőanyagokat feltételezve egy, az összes súrlódási veszteséget magában foglaló értéket vettem figyelembe (h). Csak a mozgási súrlódással számoltam, a tapadó súrlódás hatását elhanyagoltam.
- 58 -
A számítások során azért volt szükség a fojtó-visszacsapó szelepek áramlási keresztmetszeteinek felhasználására, mert a mérés sajátosságai miatt nem tudtam a nyomásokat a hengerterekben mérni. A nyomásmérő berendezéseket a kapcsolódó légvezeték egy, a hengerhez közeli szakszán helyeztem el. Ennek következtében a vizsgálat során figyelembevett holtterek lényegesen nagyobbak voltak a henger holttereinél, tartalmazták a henger holttereit, az érintett vezetékszakaszokat és a nyomásmérő eszközök tereit is. A munkahenger ki- és beáramlási veszteségeit egy ún. fojtási tényezővel (f) vettem figyelembe, amelynek értéke 0 és 1 között változhat a csatlakozó geometriai kialakításától függően [36.]. Az összehasonlítás során három paraméter (nyomások a hengerterekben és a dugattyúelmozdulás) értékeinek időbeli alakulását elemeztem, igyekezve feltárni az eltéréseket, illetve azok okait (5/10.-5/11. ábrák).
a) plusz mozgás
b) mínusz mozgás
Számított nyomás Mért nyomás
5/10. ábra Nyomás a plusz kamrában (pt=5∙105Pa)
a) plusz mozgás
b) mínusz mozgás
Számított nyomás Mért nyomás
5/11. ábra Nyomás a mínusz kamrában (pt=5∙105Pa)
- 59 -
A 5/10.-5/11. ábrákon látható, hogy a nyomásoknál eltérés tapasztalható a számított és a mért értékek között, mindkét kamra esetén. A folyamat elején az eltérés nem számottevő, a mozgás megindulása után kezd növekedni a különbség. Ennek elsősorban technikai okai vannak, mert a mérés kialakításából adódóan (5/7. ábra) a nyomásmérési pontok helye nem a munkahengerben található. A rendelkezésre álló eszközökkel csak ilyen méréskialakítást tudtam megvalósítani, de célomnak ez is megfelel, mivel a kutatásaim szempontjából az elmozdulás a legfontosabb paraméter (5/12. ábra).
a) plusz mozgás
b) mínusz mozgás
Számított érték Mért érték
5/12. ábra A dugattyú elmozdulása (pt=5∙105Pa) A nyomásoknál jelentkező eltérések ellenére a dugattyúelmozdulásra használható eredményt kaptam. A magyarázat abban keresendő, hogy a számított és a mért adatoknál a figyelembe vett terek (elsősorban a holtterek) és a veszteségek eltérőek, aminek következtében a nyomások nem fedik egymást. Mivel azonban ez a differencia mindkét tér esetén hasonló eltérést okoz, ezért a folyamat jellege azonos lesz. Ez leginkább a nyomáskülönbségek összehasonlításából látható (5/13. ábra).
a) plusz mozgás
b) mínusz mozgás
Számított nyomásoknál Mért nyomásoknál
5/13. ábra Nyomáskülönbségek (pt=5∙105Pa) - 60 -
Mivel a silóelzáró szerkezetek vizsgálatánál az elmozdulás értékének kulcsszerepe van, ezért megvizsgáltam a mért és a számított dugattyú elmozdulás között az eltérés mértékének alakulásait (5/14. ábra). Az összevetés alapja az az idő, amely a dugattyú adott lökethosszhoz tartozó elmozdulásához tartozik.
a) plusz mozgás
b) mínusz mozgás
5/14. ábra Eltérés a mért és a számított dugattyú elmozdulás között (pt=5∙105Pa) A vizsgált esetekben (5/14. ábra és M2. melléklet) megállapítható (a mért értékre vonatkoztatott %-ban), hogy a mért és a számított dugattyúelmozdulás ideje között, a mozgás első szakaszában kb. 15% (plusz mozgás), illetve kb. 25% (mínusz mozgás), a mozgás második szakaszában mindkét irányban ±5% eltérés tapasztalható. A fenti eredmények alapján, a vizsgált működtető nyomások esetén megállapíthatjuk, hogy a modell és az azt megvalósító program eredményei a dugattyú elmozdulást illetően jól közelítik a gyakorlati mérési eredményeket. Mindhárom nyomás esetén, a mozgás második szakaszán (plusz irányú mozgás estén a teljes mozgás 90%-a, mínusz irányú mozgás esetén a teljes mozgás 75%-a) az eltérés a mért és a modell alapján számított értékek között ±5%-on belül marad. A vizsgálatok alapján kijelenthető, hogy a modell a mozgási fázis folyamán - az indulási fázis kivételével - alkalmas a dugattyú elmozdulás meghatározására. Mivel értekezésemben silóelzáróknál használom a modellt, ahol a folyamat kialakítása (amint a későbbiekben látni fogjuk) lehetővé teszi az indulási fázisban jelentkező eltérések hatásának eliminálását (a zárólap nyitása, illetve zárása esetén a kiömlés csak a kisebb eltérést mutató szakasz elérése után indul meg), ezért az eljárás megfelelő a folyamat modellezésére. Emellett az is megállapítható, hogy a munkahenger kamráiban fellépő nyomások vizsgálatára ez a mérés kialakítás nem alkalmas, arra más megoldást kell keresni.
- 61 -
6. PNEUMATIKUS JÁTÉKA
MŰKÖDTETÉSŰ
ELZÁRÓ-SZERKEZETEK
ERŐ-
Az elzáró-szerkezet kialakítása jelentősen befolyásolja az adagolási folyamatot, alapvetően meghatározza: a kapcsolatot a dugattyú mozgása és az elzáró-elem elmozdulása között, az ömlesztett anyagból a dugattyúra ható terhelés nagyságát és jellegét, a tényleges kiömlési keresztmetszet alakját és méretét. A fenti hatások figyelembevételéhez fel kell írni a dugattyú elmozdulása és az elzáró-elem (síklap, vagy ív) által elzárt, illetve szabadon hagyott anyagfelület közötti függvénykapcsolatot. Ehhez meg kell határozni a kiömlési keresztmetszet alakját, melynél számos variáció előfordulhat (négyzet, téglalap, kör, ellipszis, stb.). A kiömlőnyílás alakját alapvetően a silógarat geometriája határozza meg (kúp esetén kör, gúla esetén téglalap), de a különböző járulékos szerkezeti elemekkel (toldatok, fenéklapok, stb.) ez tetszőlegesen módosítható. Az adagolási folyamat szempontjából fontos változatok: Téglalap keresztmetszet - gúla alakú garat esetén jellemző (a kiömlési sajátosságok miatt nem foglalkoztam a résen történő kiömléssel, ahol L>3∙W), illetve L=W esetén a négyzet alak is figyelembe vehető. Kör keresztmetszet - kúp alakú garat esetén jellemző. Terjedelmi és technikai korlátok miatt értekezésemben csak a téglalap keresztmetszetű kiömlőnyílás jellemzőinek meghatározását mutatom be.
6.1. Tolózárak erőjátéka A tolózárak közül két típus alkalmazható az adagolási modell esetén, de geometriai szempontból a vízszintes tolózár a ferde tolózár egy olyan eseteként is felírható, ahol az e szög értéke zérus. Ebből kifolyólag csak ferde tolózárra vonatkozóan írom fel az összefüggéseket, a vízszintes tolózár jellemzői e=0 behelyettesítésével számíthatók. A működtető elem mozgása és az elzáró elem elmozdulása közötti kapcsolat, pneumatikus működtetés esetén, a következő (6/1. ábra):
sd t ez t ,
(6.1)
ahol ez
- a tolózár elmozdulása [m].
FR(t) v h
x
1
sd(t) e
An(t) ez(t)
Lk
(md + mz) . g
6/1. ábra e szögben mozgó tolózár működésének vázlata
- 62 -
Az aktuális kiömlési felület jellemzői a dugattyú elmozdulásának függvényében (téglalap alakú kifolyónyílás, vertikális irányú anyagáramlás): terület: An t L sd t Wk cos e , (6.2) kerület:
K n t L sd t cos e Wk 2 .
(6.3)
A maximális dugattyúlöket a kifolyónyílás hosszának függvényében számítható ki:
L
Lk . cos e
(6.4)
A tolózár dőlésszögének a kiömlési felületre gyakorolt hatása a 6/2. ábrán látható.
e = 0o e = 10o e = 20o e = 30o
6/2. ábra A dugattyúelmozdulás és a kiömlési felület közötti kapcsolat (záráskor) A munkahenger dugattyúját terhelő külső erők komponensei: a munkahenger mozgó részeinek és a zárólap súlyából származó összes erő:
Fm md mz g 1 cos e sin e - zárás alatt,
(6.5)
Fm md mz g 1 cos e sin e - nyitás alatt,
(6.6)
ahol mz md
- a zárólap tömege [kg], - a munkahenger mozgó részeinek tömege [kg],
az anyag vertikális nyomásából a szerkezetre ható erő: Fv t sd t cos e Wk tv 1 x cos e sin e - zárás alatt,
(6.7)
Fv t sd t cos e Wk tv 1 x cos e sin e - nyitás alatt,
(6.8)
- 63 -
az anyag horizontális nyomásából a szerkezetre ható erő: Fh t sd t sin e Wk th 1 x sin e cos e - zárás alatt,
(6.9)
Fh t sd t sin e Wk th 1 x sin e cos e - nyitás alatt.
(6.10)
A munkahenger dugattyúját terhelő összes rúderő:
FR t Fm t Fv t Fh t .
(6.11)
Mivel a dugattyú belső súrlódási tényezője (h) kisebb, mint a zárólap és támasztófelülete közötti súrlódási tényező (1), ezért a súlyerő meghatározásánál a 1 súrlódási tényezővel közelítem a súrlódási viszonyokat.
6.2. Kéttagú szegmens-zárak erőjátéka Kéttagú szegmens-záraknál, a szimmetrikus kialakítás miatt, elegendő csak az egyik szegmensre felírni az összefüggéseket. A kiömlési jellemzők és a terhelő erők esetén az egy szegmensre vonatkozó adatok kétszerezésével kapjuk meg a szükséges értékeket. A munkahengernek a szegmenshez viszonyított helyzete befolyásolja a mozgási paramétereket, de a gyakorlat szempontjából csak a 6/3. ábrán látható elhelyezésnek van jelentőssége.
Lc
Lc
Lm
sd(t)
sd(t)
Rz
Rz
Lm+L
Lk a) a zárási folyamat kezdetén
a) a zárási folyamat végén
6/3. ábra A szegmensek elhelyezésének geometriája A két szegmensre vonatkozó aktuális kiömlési felület jellemzői a dugattyú elmozdulásának függvényében (6/4. ábra): terület: An t 2 Ln Wk 2 Wk R z sin t , (6.12) kerület: K n t 2 Ln Wk 2 4 R z sin t 2 Wk , (6.13) ahol Lc - a szegmens és a munkahenger csuklópontja közötti távolság [m], Ln - a kiömlőnyílás aktuális hossza [m], Lm - a munkahenger hossza mínusz véghelyzetben [m], - a szegmens-zár nyitási szöge [o]:
- 64 -
R 2 L2 L L z c m t arccos 2 R z Lc
2
R 2 L2 L s t c m d arccos z 2 R z Lc
2
.
(6.14)
Lc
s
Lm+ sd(t)
Rz Ln
6/4. ábra Kiömlési felület számítása egy szegmens esetén A munkahenger dugattyúját terhelő külső erők érintőirányú komponensei (6/5. ábra): egy elzáró szegmens súlyából származó erő: Fszeg t ms g 1 cos s t sin s t - zárás alatt, 2 2
(6.15)
Fszeg t ms g 1 cos s t sin s t - nyitás alatt, 2 2
(6.16)
az anyag vertikális nyomásából egy szegmensre ható erő:
Fv t R z sin s t sin t Wk tv - zárás alatt, 1 x cos s t sin s t 2 2
(6.17)
Fv t R z sin s t sin t Wk tv - nyitás alatt, 1 x cos s t sin s t 2 2
(6.18)
az anyag horizontális nyomásából egy szegmensre ható erő:
Fh t Wk R z cos t cos s t th - zárás alatt, 1 x sin s t cos s t 2 2
(6.19)
Fh t Wk R z cos t cos s t th - nyitás alatt. 1 x sin s t cos s t 2 2 ahol ms – egy szegmens tömege [kg]. - 65 -
(6.20)
Lc
s
v
1
FR(t) c
h
d
md . g
s/2
Ft(t) d
ms . g
6/5. ábra Erőjáték egy szegmens esetén A munkahenger mozgó részeinek súlyából származó erő:
ahol
Fdug t md g h sin d t cos d t - zárás alatt,
(6.21)
Fdug t md g h sin d t cos d t - nyitás alatt,
(6.22)
d t 1 t t s ,
(6.23)
R 2 L s t 2 L2 z m d c . 1 t arccos 2 R z Lm sd t
(6.24)
A munkahenger dugattyúját terhelő összes rúderő:
FR t Fdug t 2
Fszeg t Fv t Fh t
ahol c t
cos c t
1 t . 2
- 66 -
,
(6.25)
(6.26)
7. ÖMLESZTETT ANYAGOK VERTIKÁLIS KIÖMLÉSE A 3. fejezetben leírt megfontolások alapján az adagolási folyamat során, a statikus állapotot tekintve a Janssen-féle összefüggéseket [9.], a kiömlési állapotot tekintve az anyagból az elzáró-szerkezetre ható terhelésre a Janssen-féle magkifolyásra vonatkozó összefüggéseket [17.] alkalmaztam. A kiömlő anyagmennyiség meghatározására szolgáló összefüggések alkalmazhatóságának megállapítására vizsgálatokat végeztem, amelyek során az ömlesztett anyagok és az adagolási környezet lehatárolásával, illetve mérésekkel kiválasztottam az adott adagolási folyamat során felhasználható összefüggést.
7.1. Adagolt ömlesztett anyagok A 3. fejezetben felsorolt tulajdonságok alapján az ömlesztett anyagoknak számtalan változata létezik. Az adagolási folyamatnál alkalmazott anyagok kiválasztásánál a következő szempontokat kellett figyelembe venni: az ömlesztett anyag alkalmas legyen silókból történő adagolásra, olyan anyag legyen, amit a gyakorlatban ténylegesen silókban tárolnak, a kiválasztott anyagok jól reprezentálják a vizsgálni kívánt anyagtípust, stb. Az adagolás során fontos követelmény az anyaggal szemben, hogy a kiömlés folyamatosan, akadálymentesen valósuljon meg. Ennek megfelelően olyan anyagokat kell választani, amelyeknél a következő jelenségek nem fordulhatnak elő: a szemcsék összetapadása, az anyag falhoz tapadása, beboltozódás, stb. A fenti esetekben a kiömlőnyíláson átáramló anyagmennyiséget nem a kiömlési keresztmetszet, hanem más, véletlenszerűen bekövetkező jelenségek határozzák meg, ami megnehezíti az adagolási folyamat leírását. Az ömlesztett anyag szemcséinek egymáshoz, illetve a tárolósiló falához tapadását alapvetően két tényező okozhatja: az anyag nedvességtartalma, a szemcsék elektrosztatikus feltöltődése. Ömlesztett anyagok beboltozódása esetén a kifolyónyílás felett a tárolt anyagban egy boltív jön létre, amely statikailag stabil és megakadályozza az anyag kiömlését. A beboltozódás során az egymással érintkező szemcsék a fellépő, rájuk nehezedő nyomás következtében egymásnak feszülnek és stabil anyaghalmazt képeznek. A beboltozódási jelenség létrejöttét befolyásoló tényezők: az anyagszemcsék o alakja, o mérete, nedvességtartalom, a kifolyónyílás mérete, stb. A beboltózódás elkerülésére különböző bolygató-berendezéseket alkalmaznak (pl. mechanikus-, léglazításos bolygatás, stb.), amivel megakadályozható a boltozat keletkezése. A bolygatás megváltoztatja az anyag kiömlési tulajdonságait, ezért az adagolási folyamat
- 67 -
modellezéséhez olyan anyagokat kellett keresni, amelyeknél boltozat nem alakulhat ki. A beboltozódásnak két formája van [8.]: kohéziómentes anyagoknál a szemcsék mérete, alakja és a kiömlőnyílás geometriája által meghatározott esetekben ún. mechanikai boltozat jöhet létre, kohézív anyagoknál, ha az anyagban ébredő vertikális feszültség kisebb, mint az anyagban, adott helyen érvényes folyáshatárhoz tartozó feszültség értéke, akkor kohézív boltozat képződik. Kohéziómentes anyagoknál (C=0), ha a kifolyónyílás átmérője többszöröse a szemcsenagyságnak, nem áll fenn a beboltozódás veszélye. Az adagolási folyamat modellezéséhez, a fenti szempontok megfontolása után, olyan anyagtípust választottam, amely a következő tulajdonságokkal rendelkezik: - szemcsés (dp>300m), - kohéziómentes (C=0), - homogén szemcseeloszlású, - jól folyó, stb. A kiválasztott anyagtípus tulajdonságainak megfelelő anyagok közül háromféle ömlesztett anyagot használtam a vizsgálataimhoz: búza, PVC granulátum, polietilén granulátum. A vizsgálatoknál alkalmazott ömlesztett anyagok jellemzőit az 7/1. táblázat, szemcséik alakját a 7/1. ábra tartalmazza. 7/1. táblázat A vizsgált ömlesztett anyagok jellemzői [43.] Megnevezés
a [kg/m3] 750
dp [mm] 4,0
e [o] 30
x 0,25
0,6
2. PVC granulátum
660
3,0
30
0,4
0,6
3. Polietilén granulátum
580
2,5
30
0,4
0,6
1. Búza
a) búza
b) PVC granulátum
c) polietilén granulátum
7/1. ábra A vizsgált anyagok szemcséi
- 68 -
7.2. Adagolási környezet A siló, a kifolyógarat és a kiömlőnyílás geometriája jelentős hatást gyakorol az anyagkiömlésre, ezért a 3. fejezetben bemutatott kiömlési összefüggések alkalmazása az adagolási jellemzők függvénye. A silógeometriának a kiáramlás jellemzőire gyakorolt hatásának elemzése, nem tartozik jelen értekezés témájához. A modell vizsgálatához a 7/2. ábrán látható geometriával rendelkező karcsú silót választottam.
Hs
W
B
ht Wk
h0
Lk
Wk
Lk
B
W
7/2. ábra A vizsgálatoknál alkalmazott siló geometriája A tényleges geometriai méretek a következők: a silószár magassága: Hs = 1 m, a silószár szélessége: W = 0,194 m,
- 69 -
a silószár hosszúsága: a kifolyógarat magassága: a kifolyógarat dőlésszöge: a kifolyógúla elméleti magassága: a kifolyónyílás szélessége: a kifolyónyílás hossza: az ekvivalens kifolyónyílás átmérő:
B = 0,194 m, ht = 0,13 m, = 30,6o, h0 = 0,164 m, Wk = 0,04 m, Lk = 0,04 m, Dekv = 0,045 m.
Mivel a kiömlési összefüggések többsége kör alakú kiömlőnyílásra vonatkozik, ezért a számításoknál egy ekvivalens átmérővel számolok, amely az általam használt négyzettel megegyező területű kör átmérője:
Dekv
4 Lk Wk [m].
(7.1)
A siló anyaga műanyag (polimetil-metakrilát), tartószerkezete acél és mindhárom kiválasztott anyag tárolására alkalmazható. A siló elzáró-szerkezetének anyaga azonos a siló anyagával, működtetése pneumatikus munkahengerrel történik. A kiválasztás során csak állandó kiömlési keresztmetszetet vizsgáltam. A megfelelő kiömlési módszer kiválasztásához egy egyszerű mérést végeztem, amelyben a cél a kiömlő mennyiség időbeni alakulásának rögzítése volt. Ehhez felépítettem egy egyszerű berendezést, amelynek vázlata a 7/3. ábrán látható. A silóból kiömlő ömlesztett-anyag mennyiségének mérésére egy erőmérő eszközt alkalmaztam. Siló
Pneumatikus munkahenger
Mérés feldolgozás Elzáró-szerkezet
Mérőberendezés
A
B
C
D
E
F
G
H SELECTED ON-LINE
Mérési adatgyűjtő
7/3. ábra A mérési folyamat kialakítása A mérés megvalósításához a következő elemeket kellett a mérési körbe beépíteni: elzáró elem, mérőelem, mérési adatgyűjtő elem, feldolgozó elem.
- 70 -
A siló elzáró-eleme egy pneumatikus működtetésű vízszintes tolózár, amelyhez a 5. fejezetben jellemzett Mecman 166-55-0210-1 típusú kétoldali működésű hengert alkalmaztam. A kiömlő mennyiség mérésére egy HBM TYP S2 típusú erőmérőt alkalmaztam, amelynél a maximális erő: 100N, mérési érzékenység: 2mV/V. A mérési adatgyűjtés és feldolgozás az 5. fejezetben leírt módon történt.
7.3. A kiömlő mennyiség számítására alkalmas összefüggés kiválasztása A 3. fejezetben bemutatott kiömlési összefüggések főbb jellemzőit a 7/2. táblázat tartalmazza, amely alapján leszűkíthetők az adott adagolási modellben alkalmazható módszerek.
Állapot
1.
Hagen-módszer [18.]
x
x
x
x
x
2.
Beverloo-módszer [19.]
x
x
x
x
x
3.
Brit szabvány [20.]
x
x
x
x
x
4.
Crewdson-módszer [21.]
x
x
x
5.
Williams-módszer [22.]
x
x
x
x
x
6.
Zanker-módszer [23.]
x
x
x
x
7.
Gjacev-Keller módszer [24.]
x
x
x
8.
Janssen-féle módszer [17.]
x
x
x
x
x
9.
Johanson-módszer [25.]
x
x
x
x
x
x
x
x
x
x
x
x
x
10. Tomas-módszer [26.] 11. Carleton-módszer [27.]
x
12. Oldal-módszer [28.]
x
x x x
x
x
x
x
instacionér
stacionér
figyelembe veszi
Szemcsék közötti levegőt elhanyagolja
finom
durva
tömegkifolyás
Kifolyási mód magkifolyás
kohézív
Anyagtípus szabad folyású
Modellek
Szemcsenagyság
7/2. táblázat Az egyes kiömlési módszerek jellemzői
x
x
x x x
A 7/2. táblázat alapján, az általam adagolt ömlesztett anyagok (szemcsés, nem kohézív anyagok) kiömlő mennyiségének meghatározására nem használható a Tomas- és a Carletonmódszer. Hasonló megállapítás tehető a Williams- és a Zanker-módszerekkel kapcsolatban, ahol a vizsgálatoknál alkalmazott kifolyónyílás mérete (Lk=40mm) nagyobb, illetve kisebb a figyelembe vehetőnél. A Crewdson-módszer alkalmazásához szükséges a kiömlés alatt a szemcsék közötti levegő nyomásának változása, amelynek értéke az általam alkalmazott mérési eljárás eszközeivel nem határozható meg. Mivel a Gjacev-Keller módszer a kiömlési folyamat instacionér vizsgálatára szolgál és én elsősorban a stacionér szakaszt elemzem, ezért ez a módszer sem ad megfelelő eredményt.
- 71 -
A 7/2. táblázat alapján első lépésben négy módszer (Hagen-, Beverloo- és Oldal-módszer, Britt szabvány) alkalmas az általam vizsgált anyagok és rendszer esetén a kiömlő mennyiség meghatározására, de mivel a Beverloo-módszer a Hagen-módszer továbbfejlesztése és csak egy konstans szorzóban tér el tőle, ezért a kettő közül a továbbiakban csak a Beverloomódszert vizsgálom. A Brit szabvány által megadott összefüggések közül a kör alakú kiömlési nyílásra vonatkozó egyenletet használom, mivel a rés alakra vonatkozó feltétel (Lk>3∙Wk) nem teljesül a vizsgált kiömlőnyílásra (négyzet kiömlési keresztmetszet, Lk=Wk=40mm). A vizsgálataimat kiterjesztem két olyan módszerre (Janssen- és Johanson-módszer) is, amelyek alapvetően kohézív anyagokra vonatkoznak, de adott egyszerűsítések felhasználásával szemcsés, kohéziómentes anyagoknál is alkalmazhatók. A következő lépésben, mérések segítségével meghatároztam - az 5 kiválasztott összefüggés esetén - a különböző kiömlési összefüggésekben szereplő állandók (kp, Cb, és ) értékeit (7/3. táblázat). 7/3. táblázat Az összefüggések állandóinak mért értékei Megnevezés 1. Búza
kp 2,0
Cb 0,65
0,53
0,45
2. PVC granulátum
1,8
0,55
0,53
0,35
3. Polietilén granulátum
2,2
0,55
0,53
0,35
Az alaktényezőre vonatkozóan Bewerloo [19.] kp=1,4…2,9 közötti értéket kapott (kp=2,9-et homok, és kp=1,3…1,5 közötti értéket növényi magokra), a kifolyási tényezőre Cb=0,55…0,65 közötti értéket használt, de nem vizsgálta az általam használt ömlesztett anyagokat. A 7/3. táblázatból kiolvasható, hogy a Bewerloo-féle összefüggés a különböző anyagoknál milyen állandó értékek mellett szolgáltatta a legpontosabb kiömlési mennyiségeket. A Janssen-féle kiömlési összefüggés anyagtényezője jól ömleszthető anyagok esetén =0,5…0,65 [17.]. A mérések során úgy találtam, hogy mindhárom vizsgált anyag esetén a =0,53 érték adja a legjobb kiömlési eredményt. Az Oldal-féle szemcsealak-tényező értéke Oldal vizsgálatai alapján búzára =0,4; PE granulátumra =0,3 [28.]. Az általam elvégzett mérések alapján búzára =0,45; a műanyag granulátumokra =0,35 érték mellet kaptam a legpontosabb kiömlési értékeket. Karcsú silók esetén a kiömlési sebesség, a gyakorlati tapasztalatok alapján, független a tárolt anyag magasságától, de az összefüggések alkalmazhatóságának vizsgálatánál az anyagmagasság hatását is ellenőriztem. A 7/3. ábrán bemutatott mérési elrendezésben lefolytatott méréseknél, mindhárom anyag esetén 3 diszkrét beállítást alkalmaztam (Hs=0,2m; 0,4m; 0,6m). Minden beállításnál, minden anyag esetén 3 azonos mérést végeztem, az esetleges működési anomáliák kizárására. Ez alapján választottam ki egy-egy jellemző érték-sorozatot, amelyekből a 7/4. ábrán látható diagramok mutatnak be anyagonként egy-egy példát (a teljes mérés-sorozat diagramjai az M3. mellékletben találhatók). A módszerek által szolgáltatott eredmények pontosságának összehasonlítására meghatároztam a mért és a számított kiömlési mennyiségek közötti eltéréseket (7/4. táblázat). A töltési magasság hatása a kiömlő mennyiség meghatározásában, az M3. mellékletben szereplő diagramok összehasonlítása alapján - igazolva a gyakorlati tapasztalatokat - az adott siló esetén elhanyagolható, ezért vizsgálatától a későbbiekben eltekintek. - 72 -
0.5 Kiömlõ anyagmennyiség [kg]
Kiömlõ anyagmennyiség [kg]
0.5
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
0.8
0.4
0.3
0.2
0.1
0
1
0
Brit Szabvány szerint
0.2
0.4 0.6 Idõ [s]
0.8
1
0.8
1
Brit Szabvány szerint
a) búza
b) PVC granulátum
Kiömlõ anyagmennyiség [kg]
0.5
Brit Szabvány szerint Johanson szerint Beverloo szerint Janssen szerint Oldal szerint Mért értékek
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
Brit Szabvány szerint
c) polietilén granulátum
7/4. ábra Mérési eredmények (Hs=0,4m) A 7/4. ábra alapján megállapítható, hogy a Beverloo, Janssen és Oldal módszere szerint meghatározott anyagmennyiség lényegesen pontosabb közelítést ad a silóból kiömlő anyagmennyiség meghatározására, mint a Johanson, vagy a Brit Szabvány által meghatározott összefüggések (7/4. táblázat). 7/4. táblázat A kiömlési módszerek jellemző eltérései PVC Polietilén granulátum granulátum 34% 29%
Megnevezés
Búza
1. Brit szabvány (kör alakú kiömlőnyílás)
11%
2. Johanson-módszer
32%
36%
30%
3. Beverloo-módszer
3,5%
6%
1,5%
4. Janssen-módszer
1%
2,5%
-1%
5. Oldal-módszer
2%
-2%
-3%
- 73 -
Számba véve a mérési eredményeket megállapítottam, hogy az öt vizsgált összefüggés közül csak három ad olyan eredményt, amely a gyakorlati mérésekkel összevetve használható a kiömlési mennyiségek meghatározására. Mindhárom vizsgált ömlesztett-anyag tekintetében hasonló megállapítás tehető. Jóllehet a Brit szabvány lényegesen kisebb hibát vétett búza esetén, mint a Johanson-módszer, de a műanyagok esetén nem ad használható eredményt. A három módszer között alapvető különbség, amelynek a gyakorlati felhasználás szempontjából van jelentősége, hogy a Beverloo-módszer két (kp, Cb), a Janssen- és Oldal-féle módszer pedig csak egy (illetve ) tapasztalati konstans értéket használ, ami az eredmények bizonytalanságát fokozhatja (az eredmény függ a konstansok értékeinek megválasztásától). További eltérés a módszerek között, hogy a Janssen-féle módszer esetén a tapasztalati konstans () értéke csak a folyási viselkedés függvénye (tehát az általam vizsgált anyagoknál a hasonló folyási tulajdonságok miatt azonos), míg a Beverloo- és az Oldal-módszernél a Cb, illetve a értéke eltér búzára és a vizsgált műanyagokra. Ennek hatása abban jelentkezik, hogy a vizsgálatokhoz, azonos folyási viselkedésű anyagoknál is külön-külön meg kell határozni a konstansok értékét. A fenti megállapításokkal összhangban, az adagolási modellben, a kiömlő mennyiség számítására a Janssen-féle összefüggést alkalmaztam.
- 74 -
8. AZ ADAGOLÁSI FOLYAMAT ÁLTALÁNOS MODELLJE Az adagolási folyamatnak, az 1. fejezetben felvázolt kapcsolatrendszerében szereplő elemek részletes kifejtése (lásd 5., 6. és 7. fejezetek) után felírható az általános modell pneumatikus elzáró-szerkezettel ürített ömlesztett-anyag tároló silók esetén.
8.1. A modell elemei Az adagolási folyamat működési modelljét úgy állítottam elő, hogy az 5. fejezetben felírt (5/1. ábra), kétoldali működésű, lineáris dugattyús, véghelyzet fékezésű pneumatikus munkahengerek működési modelljébe beillesztettem az elzáró-szerkezet működési fázisait és összefüggéseit (8/1. ábra). 1. MŰKÖDÉSI FÁZIS
Folyamat indítása
Kezdeti töltési-ürítési állapot
nem teljesül
Bemenő paraméterek munkahenger méretek, kiinduló nyomások (pt) elzárószerkezet geometria, anyagjellemzők, silógeometria
Indulási feltétel teljesül
Elzáró-szerkezet nyitása (mínusz dugattyúmozgás)
Munkahenger 2. MŰKÖDÉSI FÁZIS Gyorsuló és állandó sebességű dugattyú mozgás
nem teljesül
Elzáró-szerkezet zárása (plusz dugattyúmozgás)
2/1. FÁZIS: Elzáró-szerkezet zárva
2/1. FÁZIS: Elzáró-szerkezet nyitva
2/2. FÁZIS: Elzáró-szerkezet nyit
2/2. FÁZIS: Elzáró-szerkezet zár
2/3. FÁZIS: Elzáró-szerkezet nyitva
2/3. FÁZIS: Elzáró-szerkezet zárva
Csillapítási feltétel teljesül
3. MŰKÖDÉSI FÁZIS Csillapított dugattyú mozgás
nem teljesül
Megállási feltétel teljesül
4. MŰKÖDÉSI FÁZIS
Kimenő jellemzők adagolási idő, kiömlött mennyiség
Végső töltési-ürítési állapot
nem teljesül
Kiegyenlítődési feltétel
teljesül
Folyamat leállítása
8/1. ábra A pneumatikus munkahengerrel működtetett adagolási folyamat modellje - 75 -
Az adagolási folyamatot, illetve az elzáró-szerkezet geometriát úgy célszerű kialakítani, hogy a záróelem mozgási folyamata teljes egészében a munkahenger 2. működési fázisa alatt menjen végbe. Ezzel a megoldással kiküszöbölhetjük a munkahenger csillapított mozgása során fellépő lengések hatását az adagolt anyagmennyiségre. Ehhez a következő feltételt kell teljesíteni: L Lk 2 sf .
(8.1)
A feltétel teljesülése esetén az adagolási folyamat a munkahenger 2. működési fázisa alatt megy végbe, ezért az 5/1. ábrán látható működési séma csak a 2. fázisban módosul. Az adagolási folyamat működési fázisai alatt a következő összefüggéseket kell figyelembe venni: az 1. fázisban a 5/2. ábrán szereplő összefüggéseket, a 2. fázisban a 5/3. ábrán szereplő és az adagolószerkezet mozgása alatt érvényes, az anyagkiömléshez kötődő jellemzőket leíró összefüggéseket, a 3. fázisban a 5/3. ábrán szereplő összefüggéseket, a 4. fázisban a 5/2. ábrán szereplő összefüggéseket. A 2. működési fázis esetén külön fel kell írni az elzáró-szerkezet nyitási és zárási folyamatának modelljét, mert a működés lépései és az alkalmazott feltételek különböznek (8/2. és 8/3. ábra). Indulási feltétel teljesül 2/1. FÁZIS: Elzáró-szerkezet zárva
p+, p- nyomás, és az s elmozdulás változik, statikus érvényes, kiömlés nincs
nem teljesül Nyitási feltétel
sd t
teljesül 2/2. FÁZIS: Elzáró-szerkezet nyit
p+, p- nyomás, és az s elmozdulás változik, kiömlési érvényes, kiömlőnyílás változik
nem teljesül Maximum kiömlés feltétel
sd t
teljesül 2/3. FÁZIS: Elzáró-szerkezet nyitva
nem teljesül
L Lk dp 2
L Lk Lk 2
p+, p- nyomás, és az s elmozdulás változik, kiömlési érvényes, kiömlőnyílás állandó
Csillapítási feltétel teljesül
8/2. ábra Az adagolás modellje nyitásnál - 76 -
Indulási feltétel teljesül 2/1. FÁZIS: Elzáró-szerkezet nyitva
p+, p- nyomás, és az s elmozdulás változik, kiömlési érvényes, kiömlőnyílás állandó
nem teljesül Zárási feltétel
sd t
teljesül 2/2. FÁZIS: Elzáró-szerkezet zár
p+, p- nyomás, és az s elmozdulás változik, kiömlési érvényes, kiömlőnyílás változik
nem teljesül Nulla kiömlés feltétel
sd t
teljesül 2/3. FÁZIS: Elzáró-szerkezet zárva
L Lk Lk 2
L Lk dp 2
p+, p- nyomás, és az s elmozdulás változik, statikus érvényes, kiömlés nincs
nem teljesül Csillapítási feltétel teljesül
8/3. ábra Az adagolás modellje zárásnál Az adagolási folyamatnak, az adagolószerkezet mozgása alatt érvényes, az anyagkiömléshez kötődő jellemzőket leíró összefüggései azonosak lesznek a zárási (8/2. ábra) és a nyitási (8/3. ábra) folyamat során, csak a fázisok felcserélődnek és a feltételek módosulnak. A nyitási, zárási, maximum és nulla kiömlési feltételek 8/2. és 8/3. ábrákon szereplő összefüggései szimmetrikus működés esetén érvényesek, amelynek feltétele, hogy L Lk 2 s f 2 s z .
(8.2)
ahol sz
- a dugattyú mozgási útja a munkahenger működésének 2. fázisában maximális kiömlés (zárás), illetve zárt szerkezet (nyitás) esetén [m].
A működési fázisok során, az anyagkiömléshez kötődő jellemzőket leíró összefüggések azonosak a nyitás 2.1. és a zárás 2.3. fázisában (8/4. ábra). Hasonlóan azonos összefüggések írhatók fel az anyagkiömléshez kötődően a nyitás 2.2. és 2.3., illetve a zárás 2.1. és 2.2. fázisában (8/5. ábra), a különbség a kiömlőnyílásnál jelentkezik: a nyitás 2.2. és a zárás 2.2. fázisában a kiömlőnyílás mérete (An) és a zárólap ömlesztett anyag által terhelt felülete (Az) változik, a nyitás 2.3. és a zárás 2.1. fázisában a kiömlőnyílás mérete állandó (An=Amax), a zárólap ömlesztett anyag által terhelt felülete zérus (Az=0). - 77 -
D) Feszültségek az anyagban
Statikus feszültség (Az = Amax) vs
tan x Us z g a As As 1 e , hs vs tan x Us n
vt
g a z g a h 0 z v0 , ht vt n 1 n 1 h0
E) Kiömlő mennyiség An = 0 F) Rúderő
Tolózár
FR t Fm t Fv t Fh t
Szegmenszár
FR t Fdug t 2
Fszeg t Fv t Fh t cos c t
8/4. ábra A nyitás 2.1. és a zárás 2.3. fázisában érvényes összefüggések (anyagjellemzők) D) Feszültségek az anyagban
Mozgási feszültség Az = Az(t)
E) Kiömlő mennyiség An = An(t)
Janssen
F) Rúderő
Tolózár
FR t Fm t Fv t Fh t
Szegmenszár
FR t Fdug t 2
2h U k A A c c k m 1 e k a b g ca c b 2 Uk a
2h U k A A k ca cb k g ca c b ma A k a 1 e Uk
Fszeg t Fv t Fh t cos c t
8/5. ábra A nyitás 2.2. és 2.3., illetve a zárás 2.1. és 2.2. fázisában érvényes összefüggések (anyagjellemzők)
8.2. A modell eredményeinek ellenőrzése Az általános adagolási modell számítógépes megvalósításánál az 5. fejezetben leírt, a pneumatikus munkahengernél alkalmazott módszert és szoftvert (MathCAD 2001i Professional) használtam fel, kiegészítve az adatokat és az összefüggéseket az ömlesztettanyag, a siló és az elzáró-szerkezet, valamint a kiömlési mód jellemzőivel. A számítási eredményeket tekintve is hasonló adatok és diagramok állnak rendelkezésre, kiegészülve a kiömlési mennyiségre vonatkozó eredményekkel. Az általános adagolási modellhez kapcsolódóan is gyakorlati mérésekkel igazoltam az elmélet használhatóságát, egy kiválasztott rendszer esetén. 8.2.1. Mérési környezet A vizsgálatok elvégzéséhez a 7. fejezetben már leírt, a 7/2. ábrán látható geometriával rendelkező karcsú silót alkalmaztam, az elzáró-szerkezet egy horizontális tolózár, melynek
- 78 -
mozgatására az 5. fejezetben ismertetett Mecman 166-55-0210-1 típusú kétoldali működésű, lineáris dugattyús, véghelyzet fékezésű pneumatikus munkahengert alkalmaztam. A mérés megvalósításához a korábban ismertetett (lásd 5. és 7. fejezetek) elemeket kellett a mérési körbe beépíteni. A mérés elrendezésének vázlata a 8/6. ábrán, kialakítása az M5. mellékletben látható.
Siló
Pneumatikus munkahenger
Útmérő
Elzárószerkezet Súlymérőberendezés
2. Nyomásmérő
1. Nyomásmérő
2. Fojtóvisszacsapó szelep
1. Fojtóvisszacsapó szelep
A
B
C
D
E
F
G
H SELECTED ON-LINE
Mérési adatgyűjtő Főszelep
Kompresszor
Mérés-feldolgozás
8/6. ábra A mérési folyamat kialakítása 8.2.2. A mérés eredményei A 8/6. ábrán bemutatott elrendezésben lefolytatott méréseknél, az 5. és 7. fejezetekben már ismertetett paraméterek mellett a következő további beállításokat alkalmaztam: a tárolt anyagmagasság: Ha=0,4m, a szerkezet és az anyag tulajdonságaitól függő tényező: n=0,7. A vizsgálatoknál alkalmazott ömlesztett anyagok jellemzőit a 7/1. táblázat (lásd 7. fejezet) tartalmazza. A méréseknél a működtető nyomást változtattam 3 diszkrét beállítással: - 79 -
Nyomások a hengerterekben [Pa]
pt=5∙105Pa; 6∙105Pa; 7∙105Pa. Minden beállításnál 3 azonos mérést végeztem az esetleges működési anomáliák kizárására. Egy jellemző érték-sorozat a teljes ciklusra (nyitás+zárás) vonatkozóan a 8/7-8/9. diagramokon látható. 8 10
5
7 10
5
6 10
5
5 10
5
4 10
5
3 10
5
2 10
5
1 10
5
0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
2.4
2.2
2.4
1.plusz hengertér kamra mínusz kamra
8/7. ábra A nyomások alakulása (búza, pt=7∙105Pa)
Dugattyú elmozdulás [mm]
110 100 90 80 70 60 50 40 30 20 10 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
8/8. ábra A dugattyú elmozdulása (búza, pt=7∙105Pa)
Kiömlött anyagmennyiség [kg]
1 0.83 0.67 0.5 0.33 0.17 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
8/9. ábra A kiömlő mennyiség (búza, pt=7∙105Pa) - 80 -
2.2
2.4
A 8/7. és 8/9. ábrák összevetéséből kitűnik, hogy a zárólap mozgása és a kiömlés között időbeni eltérés van, aminek oka a kiömlő nyílás és a mérleg közötti távolság (az anyag ennyi idő alatt teszi meg az utat a szájnyílástól a mérlegig). A 8/9. diagramon az látható, hogy a mért anyagmennyiség a folyamat kezdetén hirtelen megnő, majd a végén lecsökken. Ennek oka a kiömlő anyag impulzuserejének a mérlegre gyakorolt hatása, amely a folyamat végén megszűnik. A különböző ömlesztett anyagokra elvégzett mérések eredményei az M4. mellékletben találhatók. 8.2.3. A mért és számított eredmények összehasonlítása A számítások során az adagolási folyamat paramétereinek meghatározásához a korábban már felírt alapadatokat (lásd 5. fejezet) használtam fel. Az összehasonlítás során, céljaimhoz igazodva elsősorban a kiömlő mennyiség értékének időbeli alakulását elemeztem, igyekezve feltárni az eltéréseket, illetve azok okait. A 8/10. és 8/11. ábrák egy adott beállítás diagramjait mutatják, a többi beállításnál a kiömlő mennyiségeket tartalmazó diagramok az M4. mellékletben találhatók.
Dugattyú elmozdulás [mm]
110 96.25 82.5 68.75 55 41.25 27.5 13.75 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
2.4
számított értékek mért értékek
8/10. ábra A dugattyú elmozdulásának összevetése (búza, pt=7∙105Pa)
Kiömlött anyagmennyiség [kg]
1 0.88 0.75 0.63 0.5 0.38 0.25 0.13 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
Számított értékek Mért értékek
8/11. ábra A kiömlő mennyiség összevetése (búza, pt=7∙105Pa)
- 81 -
2.4
Az 8/10. ábra és a különböző beállításokban elvégzett mérések eredményeként megállapítható, hogy az eltérések a mért és számított dugattyúelmozdulás értékek között, a kiömlés szempontjából fontos szakaszokon, nem lépik át az 5%-ot (összhangban a 5. fejezetben leírtakkal). A 8/11. ábra (illetve az M4. melléklet diagramjai) alapján kijelenthető, hogy az összes kiömlött mennyiség számított értéke jól közelíti a mért eredményeket. A 8/1. táblázat tartalmazza az adagolási folyamat alatt kiömlött összes mennyiségnek a modell alapján számított és mért értékei közötti eltéréseket a különböző tápnyomásoknál. 8/1. táblázat A számított és mért kiömlő mennyiségek (Qö) összehasonlítása Működtető nyomás (pt)
Eltérés a mért és a számított érték között [%] búza
PVC granulátum
5∙105 Pa
2,6
1,3
polietilén granulátum -4,2
6∙105 Pa
0,3
2,2
-4,0
7∙105 Pa
-0,3
1,8
-1,4
A 8/1. táblázat alapján megállapítható, hogy az eltérések mértéke minden vizsgált esetben 5% alatt marad. Szintén megállapítható, hogy a tápnyomás növekedése kismértékben növeli az eljárás megbízhatóságát. A számított eredmények minősítésére megvizsgáltam azt a lehetőséget (8.2. táblázat), amelynél a munkahenger teljes mozgási ciklusa alatt a gyakorlatban alkalmazott, egyszerű közelítéssel, maximális kiömlési keresztmetszettel számítjuk ki a kiömlő mennyiséget, az (1.9) összefüggéshez igazodva [7.]:
Qa t n t t t z mJanssen .
(8.3)
ahol
mJanssen - a Janssen módszere szerint felírt anyagáram (3.24.) [kg/s]. 8/2. táblázat A közelítő és a mért kiömlő mennyiségek (Qö) összehasonlítása Eltérés a mért és a becsült érték között [%]
Működtető nyomás (pt)
búza
PVC granulátum
5∙105 Pa
21,7
21,4
polietilén granulátum 28,6
6∙105 Pa
18,0
34,9
32,5
7∙105 Pa
14,0
26,3
22,9
A 8/1. és 8/2. táblázat összehasonlításából jól látható, hogy a közelítésnél jelentkező eltérések jelentősen meghaladják a számított értékek eltéréseit (különösen a műanyagok esetén). Vizsgálataim eredményeként megállapítottam, hogy a kidolgozott modell alkalmas az adagolási folyamat leírására, az adott siló, elzáró-szerkezet és az alkalmazott ömlesztettanyagok esetén.
- 82 -
9. AZ ADAGOLÁSI MODELL ALKALMAZÁSA Fontos megállapítani, hogy az ömlesztett anyagok jelentősen eltérő tulajdonságai és a sokféle befolyásoló paraméter miatt, a modell gyakorlati alkalmazhatósággal kapcsolatban általános érvényű megállapításokat nem lehet tenni. Az egyes anyagtípusok esetén külön-külön ellenőrizni kell a modell eredményeit (természetesen a megfelelő kiömlési összefüggések alkalmazásával). Az értekezésben leírt eredmények szemcsés, kohéziómentes, jól folyó anyagokra érvényesek. Hasonló a helyzet a méréshez alkalmazott eszközöknél is. Az értekezésben megjelenő eredmények egy adott (kisméretű) siló és elzáró-szerkezet esetére (korlátozott anyagmagasság, kiömlő nyílás, stb.) vonatkoznak. A modellre vonatkozó általános következtetések csak különböző anyagtípusokra és silógeometriára elvégzett méréssorozatok eredményeképpen vonható le, de ennek elvégzése meghaladja az értekezés terjedelmét. Fontos alkalmazási korlát, hogy az adagolási folyamat csak akkor kezelhető a modell segítségével, ha az adagolt mennyiség sokkal kisebb a silóban tárolt mennyiségnél. Nagy mennyiségek adagolása esetén a nyitás és zárás alatt kiömlő mennyiségek bizonyos esetekben elhanyagolhatók a nyitvatartási fázis alatt kiömlő mennyiséghez képest (9/1. ábra). Ebben az esetben a durva közelítés is használható eredményt ad.
A zárás alatt kiömlõ mennyis ég [%]
20
a
15
10
5
0 0.1
0.2
0.3 0.4 Öss zes adagolt mennyiség [kg]
9/1. ábra A zárás alatt kiömlő mennyiségek aránya
9.1. Az adagolási folyamat paramétereinek meghatározása A modell elsődleges alkalmazási területe az adagolás paramétereinek meghatározása. Adott siló, elzáró-szerkezet és ömlesztett anyag esetén a modell segítségével előállítható a 9/2. ábrán látható diagram, illetve az ahhoz tartozó jellemző paraméterek: az adagolási fázisok időtartamai: - a nyitás ideje (tn), - a nyitva tartás ideje (tt), - a zárás ideje (tz), a kiömlött mennyiségek: - a nyitás alatt (Qn),
- 83 -
Qn
Qt
Qz
- a nyitva tartás alatt (Qt), - a zárás alatt (Qz),
tn
tt
tz
9/2. ábra A kiömlési diagram és a jellemző mennyiségek Közvetlen adagolás esetén a folyamat a nyitás és a zárás megindítása között eltelt idő beállításával szabályozható, amely a következő összefüggéssel számítható:
t kv t n
Qs Q n Q z tt . Qt
(9.1)
Egy egyszerű idő-vezérelt pneumatikus kapcsolás segítségével az adagoló rendszer megbízhatóan működtethető. A folyamat nem tartalmaz közvetlen visszacsatolást, ezért az adagolás pontosságának biztosításához gondoskodni kell a befolyásoló paraméterek állandó értéken tartásáról. Pneumatikus munkahengerek esetén ez gondos tervezést és üzemeltetést igényel (tápnyomás ingadozás, súrlódási viszonyok megváltozása, stb.). Mérleg-vezérelt adagolás esetén a folyamatot egy mérleg szabályozza, ahol a zárást indító mennyiség a következő összefüggéssel számítható: Qm Qs Qz .
(9.2)
Ennél a megoldásnál a mérleg ellenőrzi a kiömlő mennyiséget, amivel nagyrészt kiküszöbölhetők a paraméterváltozások negatív hatásai.
9.2. Elzáró-szerkezetek kiválasztása A modell további fontos alkalmazási területe az adagolási folyamat elemeinek kiválasztása. Különböző silók, elzáró-szerkezetek, munkahengerek esetén a modell segítségével meghatározhatók az adott feladatra alkalmazható legjobb eszközök és megoldás változatok. A különböző adagoló rendszerek azonos beállítások mellett kapott eredményeinek összehasonlításával kiválaszthatók a kívánt célfüggvényeknek leginkább megfelelő változatok. A kiválasztás során vizsgálható paraméterek: a siló geometriája,
- 84 -
az elzáró-szerkezet geometriája, a tolózár, vagy szegmens mérete, elhelyezése, a munkahenger méretei, a tápnyomás, stb.
A kiválasztás célfüggvényei lehetnek: az adagolás időtartama, az adagolás közben fellépő erők nagysága, az adagolási folyamat dinamikája, az adagolt mennyiség, stb. A következőkben egy példát mutatok be az elzáró-szerkezetek összehasonlító vizsgálatára. A vizsgálatokhoz a következő elemeket alkalmaztam: az 5. fejezetben bemutatott munkahenger, a 7. fejezetben bemutatott siló, búza, az alkalmazott elzáró-szerkezetek: - a 2. és 6. fejezetben részletezett horizontális tolózár, - a 2. és 6. fejezetben bemutatott ferde tolózár három változatban (e=10o, 20o, 30o), - a 2. és 6. fejezetben bemutatott kéttagú szegmens-zár. Ferde tolózárak esetén az egyes változatoknál számításba vett maximális lökethosszak: e = 10o smax = 0,041 m, o e = 20 smax = 0,043 m, e = 30o smax = 0,046 m. A kéttagú szegmens-zár paraméterei: Rz = 0,06 m, Lc = 0,1 m, Lm = 0,06 m, s = 19,5o, ms = 0,42 kg, 2 = 0,1. A vizsgálatokat a fent felsorolt típusok esetében, 5 különböző diszkrét tápnyomás érték (4∙105Pa; 5∙105Pa; 6∙105Pa; 7∙105Pa, 8∙105Pa) alkalmazásával végeztem el, búzára. Az eredmények a 9/1-9/3. táblázatokban láthatók. A példában az adagolási folyamathoz alkalmazott munkahenger kiválasztása szempontjából vizsgálom a különböző változatokat. A pneumatikus munkahengerek kiválasztása sok tényező függvénye (pl. alkalmazott tápnyomás, mozgás jellege, környezeti hatások, erőkifejtés, stb.), de az egyik legfontosabb paraméter a dugattyú által kifejthető erő nagysága, ami matematikai úton (közelítőleg) meghatározható, de a legegyszerűbb módszer a katalógusban feltüntetett munkahenger adatok alapján történő kiválasztás. Ehhez ismerni kell az adagolási folyamat során a munkahenger dugattyújára ható erők alakulását és a maximális rúderő függvényében kiválasztható a szükséges munkahenger típusa és paraméterei [38.]. Minél kisebb a maximális rúderő, annál kisebb munkahenger választható a feladat megoldására (a katalógusban diszkrét értékek közül választhatunk). A 9/1. táblázat tartalmazza a lefolytatott vizsgálatok során a maximális rúderők alakulását a különböző beállítások esetén. - 85 -
9/1. táblázat A maximális rúderők (FR) [N] Elzáró-szerkezet típus
pt [Pa] 4∙105
5∙105
6∙105
7∙105
8∙105
1. Horizontális tolózár
48,3
51,1
53,6
55,8
57,9
2. Ferde tolózár (e=10o)
47,6
50,4
52,9
55,2
57,3
3. Ferde tolózár (e=20o)
46,8
49,6
52,2
54,4
56,4
4. Ferde tolózár (e=30o)
45,9
48,7
51,1
53,4
55,4
5. Kéttagú szegmens-zár
30,4
36,6
40,7
43,7
46,4
A 9/1. táblázatból kiolvasható, hogy kéttagú szegmens-zár esetén lesz a legkisebb a maximális rúderő értéke. A tolózáras megoldások esetén megfigyelhető, hogy az emelkedési szög növelése kismértékben csökkenti a rúderőt. Szintén megállapítható, hogy a tápnyomás növelésével a rúderők is növekednek. 9/2. táblázat A kiömlő anyagmennyiségek (Qö) [10-3 Kg] Elzáró-szerkezet típus
pt [Pa] 4∙105
5∙105
6∙105
7∙105
8∙105
1. Horizontális tolózár
447
439
435
433
431
2. Ferde tolózár (e=10o)
447
440
436
435
433
3. Ferde tolózár (e=20o)
450
442
438
436
435
4. Ferde tolózár (e=30o)
454
445
441
439
438
5. Kéttagú szegmens-zár
419
417
417
416
416
9/3. táblázat Az adagolási idők (ta) [s] Elzáró-szerkezet típus
pt [Pa] 5
6∙105
7∙105
8∙105
1,327
1,289
1,273
1,263
1,258
2. Ferde tolózár (e=10 )
1,33
1,292
1,276
1,27
1,265
3. Ferde tolózár (e=20o)
1,346
1,304
1,287
1,279
1,276
4. Ferde tolózár (e=30o)
1,372
1,326
1,306
1,299
1,296
5. Kéttagú szegmens-zár
1,157
1,152
1,15
1,148
1,147
1. Horizontális tolózár o
4∙10
5
- 86 -
5∙10
10. ÚJ TUDOMÁNYOS EREDMÉNYEK Az értekezésben elvégzett kutatások eredményeit a következő tézisekben foglaltam össze: 1.
Kidolgoztam a kétoldali működésű, lineáris dugattyús, véghelyzet fékezésű pneumatikus munkahengerek működési modelljét. [44., 45., 46., 49., 51., 52., 53.] A pneumatikus munkahengerek általános működési jellemzőit leíró differenciálegyenletrendszer az összefüggések bonyolultsága miatt analitikusan nem oldható meg. A modell felírása során definiáltam a működés fázisait. Az egyes fázisok a munkahenger paraméterek függvényében, adott feltételek teljesülése esetén érvényesülnek, ezért megfogalmaztam a fázisváltási feltételeket. A működést leíró differenciálegyenlet-rendszer elemei az egyes fázisokban eltérő alakot vesznek fel, ezért meghatároztam az adott fázisokban érvényes összefüggéseket és felírtam a numerikus megoldáshoz szükséges kapcsolatokat.
2.
Mérésekkel igazoltam a pneumatikus modelljének eredményeit. [44., 46., 49., 52.]
munkahengerek
elméleti
működési
A modell alapján számítógépes eljárást dolgoztam ki MathCAD szoftver környezetben, amellyel a munkahengerek működési jellemzői numerikusan meghatározhatók. A modell alapján számított és a mért eredmények összehasonlításával megállapítottam, hogy az elméleti modell a dugattyú elmozdulást tekintve, jól közelíti a mérési eredményeket (az eltérés minden vizsgált esetben kisebb, mint 5%) plusz mozgás esetén, a mozgási szakasz 90%-án, mínusz mozgás esetén a mozgási szakasz 75%-án. Az eredmények igazolták, hogy a működési modell alkalmas az adagolási folyamatban történő felhasználásra, abban az esetben, ha az elzáró-szerkezet mozgása a dugattyúmozgás második szakaszára esik (ez megfelelő kialakítás esetén biztosítható). 3. Felírtam a dugattyúerő meghatározásához szükséges összefüggéseket az adagolási folyamatban alkalmazott elzáró-szerkezetek esetén. [47., 55., 56.] A silóelzáró-szerkezetek vizsgálatával kiválasztottam a pneumatikus munkahengerrel működtetett adagolási folyamatnál alkalmazható típusokat (vízszintes és ferde tolózár, kéttagú szegmens-zár). Kimutattam, hogy az elzáró-szerkezetek működtetésénél eltérő összefüggéseket kell alkalmazni az anyag nyomásából adódó terhelés meghatározására a statikus és a mozgási állapot alatt. A kiválasztott elzáró-szerkezetek esetén meghatároztam az ömlesztett anyagban ébredő feszültségekből a munkahenger dugattyújára ható erők figyelembevételére alkalmas összefüggéseket: ferde tolózárra: FR t Fm t Fv t Fh t ,
kéttagú szegmens-zárra: FR t Fdug t 2
ahol Fm
Fszeg t Fv t Fh t cos c t
,
- a munkahenger mozgó részeinek és a zárólap súlyából származó erő [N], - 87 -
Fv Fh Fdug Fszeg c
- az anyag vertikális nyomásából az elzáró elemre ható erő [N], - az anyag horizontális nyomásából az elzáró elemre ható erő [N], - a munkahenger mozgó elemeinek súlyából származó erő [N], - az elzáró szegmens súlyából származó erő [N], - a szegmens mozgása és a munkahenger helyzete által meghatározott szög [ o].
Felírtam a dugattyúerő meghatározásánál szerepet játszó erőkre vonatkozó összefüggéseket, valamint definiáltam a dugattyú elmozdulás és az elzáró-szerkezetek mozgása közötti kapcsolatokat. 4. Kidolgoztam a pneumatikus munkahengerrel működtetett silóürítési folyamatok, szemcsés, kohéziómentes, szabad folyású ömlesztett anyagok adagolására alkalmas általános modelljét. [50., 54., 56., 58., 60., 62., 64.] A modell megalkotásához feltártam az adagolási folyamat elemeit, definiáltam az egyes elemek szerepét, jellemzőit és meghatároztam a közöttük fennálló kapcsolatokat. A modell kidolgozása során elemeztem a folyamat működését és definiáltam az adagolás fázisait. A folyamat során az egyes fázisok az elzáró-szerkezet mozgásának függvényében, adott feltételek teljesülésekor érvényesülnek, ezért megfogalmaztam a fázisváltási feltételeket mind a nyitás, mind a zárás esetére. Az adagolási folyamat fázisait beillesztettem, a pneumatikus munkahenger működési folyamatába és meghatároztam az egyes fázisokban érvényes összefüggéseket. 5.
Mérésekkel igazoltam az adagolási folyamat általános elméleti modelljének eredményeit. [62., 63., 64.] A modell alapján számítógépes eljárást dolgoztam ki MathCAD szoftver környezetben, amellyel az adagolási folyamat jellemzői numerikusan meghatározhatók. A számított és a mért eredmények összehasonlításával megállapítottam, hogy az elméleti modell az adagolási folyamat alatt kiömlő mennyiséget tekintve jól közelíti a mérési eredményeket (az eltérés minden vizsgált esetben kisebb, mint 5%). Az eredmények igazolták, hogy az adagolási modell a három vizsgált anyag (búza, PVC granulátum és polietilén granulátum), a vizsgált siló és vízszintes tolózár esetén a kiömlő anyagmennyiség tekintetében megfelelően írja le az adagolási folyamatot.
- 88 -
ÖSSZEFOGLALÁS Értekezésem célja az adagolási folyamat általános leírására alkalmas modell kidolgozása volt pneumatikus elzáró-szerkezettel ürített ömlesztett-anyag tároló silók esetén. Mivel ez a téma eddig sem a hazai, sem a nemzetközi szakirodalomban nem kapott megfelelő publicitást, ezért kutatásom fontos feladata volt a kapcsolódó három szakterület irodalmában történő elmélyülés és a szűkebb témához szükséges szintézis megteremtése. A kutatás során az alábbi eredmények születtek: 1. Kidolgoztam a kétoldali működésű, lineáris dugattyús, véghelyzet fékezésű pneumatikus munkahengerek működési modelljét. 2. Mérésekkel igazoltam a pneumatikus munkahengerek működési modelljének eredményeit. 3. Felírtam az adagolási folyamatban alkalmazható elzáró-szerkezeteknél a dugattyúerő meghatározásához szükséges összefüggéseket. 4. Kidolgoztam a pneumatikus munkahengerrel működtetett silóürítési folyamatok, szemcsés, kohéziómentes, szabad folyású ömlesztett anyagok adagolására alkalmas általános modelljét. 5. Mérésekkel igazoltam az adagolási folyamat elméleti modelljének eredményeit. A kidolgozott modell alkalmazásával: adott siló, elzáró-szerkezet és ömlesztett anyag esetén meghatározható az adagolási folyamat szabályozásához szükséges - ciklusidő (közvetlen adagolás esetén), vagy - mennyiség (mérleg-vezérelt adagolás esetén), a különböző silók, elzáró-szerkezetek és munkahengerek esetén meghatározhatók az adott feladatra leginkább alkalmas eszközök és megoldás változatok. Teljes és mindenre kiterjedő elemzést nem folytattam, mert értekezésem célja elsősorban a modell kidolgozása és alkalmazhatóságának bemutatása volt. Mivel a vizsgálatokat egy adott anyagcsoport, siló-geometria és elzáró-szerkezet esetén végeztem el, ezért az alkalmazhatóság is e környezetben bizonyított. A kidolgozott modell közvetlen és mérleg-vezérelt adagolás esetén alkalmazható hatékonyan, számottevő előny elsősorban olyan ipari feladatoknál (élelmiszeripar, vegyipar, stb.) jelentkezhet, ahol más, bonyolultabb megoldásokat alkalmaznak az ömlesztett anyagok adagolására. Ahhoz, hogy általános érvényű megállapításokat tehessünk, a vizsgálatokat ki kell terjeszteni más anyagokra, silókra és elzáró-szerkezetekre. A legfontosabb továbbfejlesztési irányok: egyéb, kohéziómentes ömlesztett anyagok vizsgálata, kohézív ömlesztett-anyagok ürítésének vizsgálata, valós folyamatokban működő silók vizsgálata, az elzáró-szerkezet geometria és a kiömlési folyamat jellemzői közötti kapcsolatok vizsgálata, DEM módszer alkalmazási lehetőségének vizsgálata.
- 89 -
SUMMARY The aim is my Ph.D. work was to describe the general model of dosing processes for bulk solids stored in silos using pneumatic actuator operated discharging devices. As not too many papers related to this theme were published so far in the international literature, so one of my important tasks was to get knowledge in all of the three fields (pneumatics, bulk materials, handling of bulks) related to my topic and make a suitable synthesis of them required by the specific subject. New results of my research work: 1. I have described a general operation model for double acting pneumatic cylinders with cushioning. 2. I made a computer program to calculate the parameters of the cylinder based on the model. To evaluate the results I also made a measuring system and I have proved that the results are valid. 3. I have described the calculation forms of the forces acted on the piston rod at the applied discharging structures (shutter and double segment gate). 4. I have described a general model for the dosing process of bulk solids stored in silos using double acting pneumatic cylinder operated discharging devices. The model is valid for granulated, cohesion free, free flowing bulk solids. 5. To evaluate the results I made a computer program based on the dosing model and a measuring system. I have proved that the results are valid. Applying the developed dosing model at given silo, discharging device and bulk solid we can determine the control parameters of the dosing process: - cycle time (at direct dosing), or - quantity (at dosing with mass control), at different silos, discharging devices and pneumatic actuators we can determine the optimal equipment and solution varieties for a given task. Results are not valid for all of the varieties of the model elements, because my aim was to develop the model and prove its applicability. As the analysis was made at certain bulk solids, silo geometry and discharging device, the applicability is proved only at this settings. Application of the model has advantages mainly at direct and mass controlled dosing processes and at industrial processes (food industry, chemical industry, etc.) where other, more complex solutions are used for dosing of bulk solids. For general evaluation, the analysis of the model has to be enlarged for other materials, silos and discharging devices. Most important directions for the development in the future: analysis of other cohesion free bulk solids, analysis of cohesive bulk solids, analysis of real-sized industrial silos, analysis of the effects of the geometry of the discharging devices, analysis of the applicability of a DEM method for the dosing process, etc.
- 90 -
KÖSZÖNETNYILVÁNÍTÁS Doktori értekezésem elkészítésében nyújtott segítségéért szeretnék köszönetet mondani témavezetőmnek, Dr. Szaladnya Sándor professzor úrnak, aki a doktori kutatásaim során szakmailag és emberileg is következetesen támogatott. Az Ő vezetése és iránymutatása felbecsülhetetlen segítséget jelentett a munkám során. Szintén köszönetet kell mondanom Dr. Németh János docens úrnak, aki az ömlesztett anyagokkal kapcsolatos kutatásaimban nyújtott nélkülözhetetlen segítséget, és Bozzay Péter tanszéki mérnöknek, aki a gyakorlati mérések végrehajtásában segédkezett. Köszönet illeti még Dr. Gömze A. László egyetemi tanárt, a Kerámia- és Szilikátmérnöki Tanszék vezetőjét, aki rendelkezésemre bocsátotta a mérések során felhasznált modell-silót, és Dr. Barna Balázs tanszéki mérnököt, aki a szükséges pneumatikus elemeket hozzáférhetővé tette. Végül, de nem utolsó sorban köszönetet kell mondanom Dr. Illés Béla egyetemi tanárnak, az Anyagmozgatási és Logisztikai Tanszék vezetőjének és a tanszék minden kollégájának, amiért a doktori munkám során támogattak és segítettek. Külön meg kell említenem feleségemet és családomat, akiktől az értekezés elkészítésére fordított idő egy jelentős részét elvettem, és akikre mindig támaszkodhattam. Az értekezés a TÁMOP-4.2.2/B-10/1-2010-0008 jelű projekt részeként – az Új Magyarország Fejlesztési Terv keretében – az Európai Unió támogatásával, az Európai Szociális Alap társfinanszírozásával valósult meg.
- 91 -
FELHASZNÁLT IRODALMAK JEGYZÉKE [1.]
Pahl, M. H.: Lagern, Fördern und Dosieren von Schüttgüten. Verlag TÜV Rheinland, Köln, 1989. ISBN 3-88585-672-7
[2.]
Felföldi L. (szerk.): Anyagmozgatási kézikönyv. Műszaki Könyvkiadó, Budapest 1975. ISBN 963 10 0423 6
[3.]
Reimbert, M.; Reimbert, A.: Silos. Berechnung, Betrieb und Ausführung. Bauverlag Gmbh. Wiesbaden und Berlin1975. ISBN 3-7625-0369-9
[4.]
Andrási Gy.; Polonszky Gy.: Tárolás. Termékek fogadása, raktározása, megóvása. Műszaki Könyvkiadó, Budapest 1971.
[5.]
Prezenszki J.: Raktározástechnika. Műszaki Könyvkiadó, Budapest, 1984. ISBN 963 10 57364
[6.]
Reisner, W.; Eisenhart Rothe, M. V.: Bins and Bunkers for Handling Bulk Materials. Practical Design and Techniques. Trans Tech Publications, Cleveland, Clausthal-Zellerfeld, 1971. ISBN 0-87849-001-9
[7.]
Knoll I.: Anyagmozgatás a mezőgazdaságban. Mezőgazdasági Kiadó, Budapest, 1983. ISBN 963 231 479 4
[8.]
Woodcock, C. R.; Mason, J. S.: Bulk solids handling. Leonard Hill, Glasgow and London, 1987.
[9.]
Janssen, H. A.: Getreidedruck in silozellen. Zeitschr. d. Vereines deutscher Ingenieure 39. pp.: 1045-1049. 1895.
[10.] Schulze, D.: Untersuchungen zur gegenseitigen Beeinflussung von Silo Austragorgan. Dissertation, Fakultät für Maschinenbau und Elektrotechnik der Technischen Universität Carolo-Wilhelmina zu Braunschweig, 1991. [11.] Sitkei Gy.: A mezőgazdasági anyagok mechanikája. Akadémiai Kiadó, Budapest, 1981. ISBN 963 05 2673 5 [12.] Cundall, P. A.; Strack, O. D. L.: A discrete numerical model for granular assemblies. Geotechnique, 29. (1979) 1., pp.: 47-65. [13.] Magyar Szabványügyi Testület (MSZT): Eurocode 1: A tartószerkezeteket érő hatások. 4. rész: Silók és tartályok. MSZ EN 1991-4:2006. [14.] Schulze, D.: Powders and Bulk Solids. Springer-Verlag, Berlin Heidelberg 2008. ISBN 978-3540-73767-4 [15.] Oldal I.: Szemcsés anyagok kifolyási és boltozódási tulajdonságai. Ph.D értekezés. Szent István Egyetem, Gödöllő, 2007. [16.] Arnold, P. C.; McLean, A. G.: Improved Analytical Flowfactors for Mass-Flow Hoppers. Powder Technology 15. (1976) pp.: 279-281. [17.] Németh J.: Ömlesztett-anyagok rakodása és tárolása. Előadás fóliák. Miskolci Egyetem, 2000. pp.: 1-132. [18.] Hagen, E.: Druck und Bewegung des trockenen Sandes. In: Bericht über die zur Bekanntmachung geeigneten Verhandlungen der Königlich Preuβischen Akademie der Wissenschaften zu Berlin, pp.: 35-42. 1852. [19.] Bewerloo, W. A.; Leniger, H. A.; van der Felde, J.: The flow of granular solids through orifices. Chem. Engineering Sciences 15. pp.: 260-269. 1961. [20.] British Standards Institution (Ed): Draft design code for silos, bins, bunkers and hoppers. ICS 65.020.40 Published by BSI in association with the British Materials Handling Board 1987.
- 92 -
[21.] Crewdson, B. J.; Ormond, A. L.; Nedderman, R. M.: Air-impeded discharge of fine particles from a hopper. Powder Techn. 16. pp.: 197-207. 1977. [22.] Williams, J. C.: The rate of discharge of coarse granular materials from conical massflow hoppers. Chem. Engineering Sciences 32. pp.: 247-255. 1977. [23.] Zanker, A.: Estimating the flow of solids through openings. Process Engineering (July) pp.: 6667. 1975. [24.] Keller H.: Korrekturformeln zur Theorie des Schuttgutausflusses aus Bunkern von L.V.Gjachev / H. Keller // Wissenschaftlice Zeitschrift der Hochschule fur Architectur und Bauwessen. Weimar, 1987. - Raie B.h. 5/6. - S. 292295. [25.] Johanson, J. R.: Method of calculating rate of discharge from hoppers and bins. Trans. Min. Engineers AIME 232. pp.: 69-80. 1965. [26.] Scheibe, M.: Die Fördercharasteristik einer Zellenradschleuse unter Berücksichtigung der Wechselwirkung von Silo und Austragsorgan. Dissertation. Fakultät für Maschinenbau, Verfahrens- und Energietechnik, Technischen Universität Bergakademie Freiberg, 1997. [27.] Carleton, A. J.: The effect of fluid drag forces on the discharge of free flowing solids from hoppers. Powder Technol. 6. pp.:91-96. 1972. [28.] Oldal I.; Keppler I.; Csizmadia B.; Fenyvesi L.: Outflow properties of silos: The effect of arching. Advanced Powder Technology (2011), doi:10.1016/j.apt.2011.03.013 [29.] Barber, A.: Pneumatic Handbook. Elsevier Advanced Technology, Oxford 2006. ISBN 185617249X [30.] Vida Gy.: Műszaki hőtan. Tankönyvkiadó, Budapest 1991. [31.] Czibere T.: Áramlástan. Tankönyvkiadó, Budapest 1990. [32.] Zoebl, H.; Kruschik, J.: Áramlás csövekben és szelepekben. Műszaki Könyvkiadó, Budapest 1986. ISBN 963 10 6580 4 [33.] Baranyi L.: Ideális gáz izentropikus kiáramlása hőszigetelt tartályból. Oktatási segédlet. Nehézipari Műszaki Egyetem, Miskolc 1984. [34.] MECMAN Pneumatika Katalógus. Eurobridge-Typopress Kft. [35.] Blackburn, J. F.; Reethof, G.; Shearer, J. L.: Fluid Power Control, The Technology Press of MIT/Wiley, New York, 1960. [36.] Elek I., Hudáky J.: Az ipari pneumatika alapjai. Interpress, Budapest 1979. [37.] Magyar Szabványügyi Testület (MSZT): Pneumatikus energiaátvitel. Összenyomható munkaközeggel működő elemek. Áramlási jellemzők meghatározása. MSZ ISO 6358: 1995. [38.] Lee, H. K.; Choi, G. S.; Choi, G. H.: A study on tracking position control of pneumatic actuators. Mechatronics 12 (2002) pp.: 813–831. Elsevier Science Ltd. [39.] Lin-Chen, Y.Y.; Wang, J.; Wu, Q. H.: A software tool development for pneumatic actuator system simulation and design. Computers in Industry 51. (2003) pp.: 73–88. [40.] Sorli, M.; Gastaldi, L.; Codina, E.; de las Heras, S.: Dynamic analysis of pneumatic actuators. Simulation Practice and Theory 7. (1999) pp.: 589-602. [41.] Bálint E.: Közelítő matematikai módszerek. Műszaki Könyvkiadó, Budapest 1966. [42.] FESTO Pneumatika gyorskatalógus. 2003. [43.] Martens, M. (edit.): Silo-Handbuch. Ernst&Sohn, Berlin, 1998.
- 93 -
AZ ÉRTEKEZÉSHEZ KAPCSOLÓDÓ PUBLIKÁCIÓK [44.] Szaladnya S., Németh J., Telek P.: Pneumatikus elemek tereinek töltése és ürítése. PneuHidro'96 Konferencia kiadvány, pp.: 83-86., Miskolc, 1996. [45.] Szaladnya S., Németh J., Telek P.: Pneumatikus munkahengerek működésvizsgálata. MicroCAD'96 Konferencia, Anyagmozgatás, logisztika, robotika szekció kiadványa, pp.: 93102., Miskolc, 1996. [46.] Szaladnya S., Németh J., Telek P.: Silók pneumatikus zárószerkezetének tervezése. Géptervezés-termékfejlesztés’96 szeminárium kiadványa, pp.: 239-245. Miskolc 1996. [47.] Szaladnya S., Németh J., Telek P.: Silók pneumatikus zárószerkezetének méretezése. MicroCAD'97 Konferencia, Szilikátipar és porkohászat szekció kiadványa, pp.: 33-41., Miskolc, 1997. [48.] Szaladnya S., Telek P.: Szilikátiparban alkalmazott pneumatikus munkahengerek levegőfelhasználásának optimalizálása. Gép IL évfolyam, 1998/4-5. pp.: 82-83., ISSN 0016-8572 [49.] Szaladnya S., Telek P.: Pneumatikus munkahengerek alkalmazás-specifikus vizsgálata. FMTÜ Konferencia kiadvány, pp.: 241-244., Kolozsvár, 1998. [50.] Szaladnya S., Telek P.: Dimensioning principals of pneumatic system controlling the quantity of outflow material discharging from bulk solid storage silos. Miskolci Beszélgetések ‘99 konferencia kiadványa, pp.: 99-103., Miskolc, 1999. [51.] Szaladnya S., Németh J., Telek P.: Pneumatikus ütőhengerek méretezésének matematikai alapjai. FMTÜ Konferencia kiadvány, pp.: 109-112., Kolozsvár, 2000. [52.] Szaladnya S., Németh J., Telek P.: Pneumatikus ütőhengerek méretezése. Gépgyártástechnológia XL. évfolyam, pp.: 5-8., 2000. szeptember 9. HU ISSN 0016-8580 [53.] Szaladnya S., Telek P.: Effects of random parameter changing of pneumatic devices to the operation of materials handling equipment. MicroCAD 2003. Konferencia, O: Anyagáramlási rendszerek és logisztikai informatika szekció kiadványa. pp.: 57-62., Miskolc, 2003. [54.] Szaladnya S., Telek P.: Theoretical principals of computer analysis of silo discharging devices operated by pneumatic actuators. MicroCAD 2007. Konferencia, N: Anyagáramlási rendszerek és logisztikai informatika szekció kiadványa, pp.: 115-121. Miskolc, 2007. [55.] Telek P.: Silóelzáró szerkezetek elmozdulása és a működtető pneumatikus munkahenger dugattyújának mozgása közötti kapcsolat meghatározása. GÉP LVIII. évfolyam 2007/5-6. ISSN 0016-8572 pp.: 71-77. [56.] Szaladnya S., Telek P.: Comparison analysis of silo discharging devices operated by pneumatic actuator. Budowa Maszyn i Zarzadanie Produkcja No. 8. 2008. ISSN 1733-1919 pp.: 101-111. [57.] Szaladnya S., Telek P.: Pneumatikus automatizálás eszközei, a tervezés módszerei, exkluzív példatár. Építésügyi Tájékoztatási Központ Kft., ISBN 978-963-513-216-4 pp.: 1-173. Budapest, 2009. [58.] Telek P.: Ömlesztett anyagok kiömlési modelljei gravitációs silóürítés esetén. MicroCAD 2009. Konferencia, O: Anyagáramlási rendszerek, Logisztikai informatika szekció kiadványa, pp.: 203-209. Miskolc, 2009. [59.] Telek P.: Logistical criterions and strategies for storage of bulk materials. ADVANCED LOGISTIC SYSTEMS, Theory and Practice, Vol. 3. Edited by B. Illés, B. SkowronGrabowska, P. Telek, pp.: 51-56., HU ISSN 1789-2198, Miskolci Egyetem 2009. [60.] Illés B., Németh J., Telek P.: Evaluierung der Berechnungsmethoden zum Austrag von Massenfluss-Silos. Innovationen gegen die Krise. 14. Fachtagung Schüttgutfördertechnik in Magdeburg, pp.: 175-184. Otto-von-Guericke-Universität Magdeburg, 2009.
- 94 -
[61.] Telek P.: Ömlesztett-anyag tárolás a korszerű termelési folyamatokban. (Korszerű raktározási technikák V.) TRANSPACK Csomagolási, anyagmozgatási, logisztikai szaklap 2010/3. szám pp.: 40-42. HU ISSN 1587-7183, Horizont Média Kft. [62.] Illés B., Németh J., Telek P.: Computer aided calculation method for the silo discharge by a pneumatic actuator. Anlagen – Bauteile - Computersimulationen. 16. Fachtagung Schüttgutfördertechnik Magdeburg, pp.: 221-230. Otto-von-Guericke-Universität Magdeburg, 2011. [63.] Telek P.: Számítógépes módszer pneumatikus munkahengerrel történő silóürítés méretezésére. XXVI. microCAD International Scientific Conference, J. section: Material Flow Systems. Logistical Information Technology and Technical Language, University of Miskolc, 29-30 March, 2012. CD kiadvány. ISBN: 978-963-661-773-8 [64.] Telek P.: General model of silo discharges operated by a pneumatic actuator. ADVANCED LOGISTIC SYSTEMS, Theory and Practice, Vol. 6. Edited by B. Illés, S. Kot, P. Telek, HU ISSN 1789-2198, Miskolci Egyetem 2012. (Megjelenésre elfogadva)
- 95 -
MELLÉKLETEK
M1. Melléklet MUNKAHENGER PARAMÉTEREK MÉRÉSI EREDMÉNYEI M1/1. Mérési eredmények pt=5∙105Pa esetén
6 10
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
Nyomások [Pa]
Nyomások [Pa]
A nyomások alakulása:
0 4.73
4.93
5.13 5.33 Idõ [s]
5.53
0 7.85
5.73
Nyomás az 1. hengertérben Nyomás a plusz kamrában Nyomás a mínusz kamrában
8.05
8.25
8.45 8.65 Idõ [s]
8.85
9.05
Nyomás az 1. hengertérben
a) plusz mozgás
b) mínusz mozgás
0.12
0.12
0.1
0.1 Elmozdulás [m]
Elmozdulás [m]
A dugattyúelmozdulás alakulása:
0.08 0.06 0.04 0.02 0 4.73
0.08 0.06 0.04 0.02
4.93
5.13 5.33 Idõ [s]
5.53
0 7.85
5.73
a) plusz mozgás
8.05
8.25
8.45 8.65 Idõ [s]
b) mínusz mozgás
- 97 -
8.85
9.05
M1/2. Mérési eredmények pt=6∙105Pa esetén
7 10
5
7 10
5
6 10
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
Nyomások [Pa]
Nyomások [Pa]
A nyomások alakulása:
0 4.99
5.19
5.39 5.59 Idõ [s]
5.79
0
5.99
8.2
8.4
Nyomás az 1. hengertérben Nyomás a plusz kamrában Nyomás a mínusz kamrában
8.6
8.8 Idõ [s]
9
9.2
9.4
Nyomás az 1. hengertérben
a) plusz mozgás
b) mínusz mozgás
0.12
0.12
0.1
0.1 Elmozdulás [m]
Elmozdulás [m]
A dugattyúelmozdulás alakulása:
0.08 0.06 0.04 0.02 0 4.99
0.08 0.06 0.04 0.02
5.19
5.39 5.59 Idõ [s]
5.79
0
5.99
a) plusz mozgás
8.2
8.4
8.6
8.8 Idõ [s]
9
b) mínusz mozgás
- 98 -
9.2
9.4
M1/3. Mérési eredmények pt=7∙105Pa esetén
8 10
5
8 10
5
7 10
5
7 10
5
6 10
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
Nyomások [Pa]
Nyomások [Pa]
A nyomások alakulása:
0 4.29
4.49
4.69 4.89 Idõ [s]
5.09
0 7.39
5.29
Nyomás az 1. hengertérben Nyomás a plusz kamrában Nyomás a mínusz kamrában
7.59
7.79
7.99 8.19 Idõ [s]
8.39
8.59
Nyomás az 1. hengertérben
a) plusz mozgás
b) mínusz mozgás
0.12
0.12
0.1
0.1 Elmozdulás [m]
Elmozdulás [m]
A dugattyúelmozdulás alakulása:
0.08 0.06 0.04 0.02 0 4.29
0.08 0.06 0.04 0.02
4.49
4.69 4.89 Idõ [s]
5.09
0 7.39
5.29
a) plusz mozgás
7.59
7.79
7.99 8.19 Idõ [s]
b) mínusz mozgás
- 99 -
8.39
8.59
M2. Melléklet MÉRT ÉS SZÁMÍTOTT MUNKAHENGER PARAMÉTEREK ÖSSZEHASONLÍTÁSA M2/1. Mért és számított eredmények összehasonlítása pt=5∙105Pa esetén Számított érték Mért érték
6 10
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
Nyomások [Pa]
Nyomások [Pa]
Nyomás a plusz kamrában:
0
0
0.2
0.4 0.6 Idõ [s]
0.8
0
1
0
Számított nyomás az 1. hengertérben
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
Nyomások [Pa]
Nyomások [Pa]
6 10
0
0.2
0.4 0.6 Idõ [s]
0.8
0
1
Számított nyomás a 2. hengertérben 0.12
0.1
0.1
0.08 0.06 0.04 0.02 0
0.2
0.8
1
1.2
0
0.2
0.4
0.6 Idõ [s]
0.8
1
1.2
Számított nyomás a 2. hengertérben
0.12
0
0.6 Idõ [s]
tmb j
Dugattyú elmozdulás [m]
Dugattyú elmozdulás [m]
Dugattyúelmozdulás:
0.4
Számított nyomás az 1. hengertérben
Nyomás a mínusz kamrában:
0
0.2
0.4 0.6 Idõ [s]
0.8
0.08 0.06 0.04 0.02 0
1
Számított a) pluszérték mozgás
0
0.2
0.4
0.6 Idõ [s]
0.8
Számított érték b) mínusz mozgás
- 100 -
1
1.2
5 10 5 4 10
5
5 10 5 4 10
5
5
3 10 5 2 10 5 1 10
3 10 5 2 10 5 1 10
Nyomáskülönbség [Pa]
Nyomáskülönbség [Pa]
Nyomáskülönbségek:
0 5
1 10 5 2 10 5 3 10 5 4 10 5
5
0 5
1 10 5 2 10 5 3 10 5 4 10 5
5 10
0
0.2
0.4 0.6 Idõ [s]
0.8
5 10
1
Számított nyomásoknál Mért nyomásoknál
0
0.24
0.48 0.72 Idõ [s]
0.96
1.2
Számított nyomásoknál
a) plusz mozgás
b) mínusz mozgás
Eltérés a mért és a számított dugattyú elmozdulás között 25
25
20
20
15
15
10
10
5
5
0
0
5
5
10
0
0.02 0.04 0.06 0.08 Dugattyú elmozdulás [m]
0.1
a) plusz mozgás
10
0.1
0.08 0.06 0.04 0.02 Dugattyú elmozdulás [m]
b) mínusz mozgás
- 101 -
0
M2/2. Mért és számított eredmények összehasonlítása pt=6∙105Pa esetén Számított érték Mért érték
7 10
5
7 10
5
6 10
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
0
Nyomások [Pa]
Nyomások [Pa]
Nyomás a plusz kamrában:
0
0.2
0.4 0.6 Idõ [s]
0.8
0
1
0
0.2
0.4
0.6 Idõ [s]
0.8
1
1.2
Számított nyomás az 1. hengertérben
Számított nyomás az 1. hengertérben
7 10
5
7 10
5
6 10
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
0
Nyomások [Pa]
Nyomások [Pa]
Nyomás a mínusz kamrában:
0
0.2
0.4 0.6 Idõ [s]
0.8
0
1
0.1
0.1
Dugattyú elmozdulás [m]
Dugattyú elmozdulás [m]
0.12
0.08 0.06 0.04 0.02 0
0.2
0.4
0.6 Idõ [s]
0.8
1
1.2
tmb j
0.12
0
0.2
Számított nyomás a 2. hengertérben
Számított nyomás a 2. hengertérben
Dugattyúelmozdulás:
0
0.4 0.6 Idõ [s]
0.8
0.08 0.06 0.04 0.02 0
1
0
0.2
0.4
0.6 Idõ [s]
0.8
Számított érték b) mínusz mozgás
Számított a) pluszérték mozgás
- 102 -
1
1.2
5 10 5 4 10
5
5 10 5 4 10
5
5
3 10 5 2 10 5 1 10
3 10 5 2 10 5 1 10
Nyomáskülönbség [Pa]
Nyomáskülönbség [Pa]
Nyomáskülönbségek:
0 5
1 10 5 2 10 5 3 10 5 4 10 5
5
0 5
1 10 5 2 10 5 3 10 5 4 10 5
5 10
0
0.2
0.4 0.6 Idõ [s]
0.8
1
5 10
Számított nyomásoknál Mért nyomásoknál
0
0.24
0.48 0.72 Idõ [s]
0.96
1.2
0.08 0.06 0.04 0.02 Dugattyú elmozdulás [m]
0
Számított nyomásoknál
a) plusz mozgás
b) mínusz mozgás
Eltérés a mért és a számított dugattyú elmozdulás között 25
25
20
20
15
15
10
10
5
5
0
0
5
5
10
0
0.02 0.04 0.06 0.08 Dugattyú elmozdulás [m]
0.1
a) plusz mozgás
10
0.1
b) mínusz mozgás
- 103 -
M2/3. Mért és számított eredmények összehasonlítása pt=7∙105Pa esetén Számított érték Mért érték
8 10
5
8 10
5
7 10
5
7 10
5
6 10
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
0
Nyomások [Pa]
Nyomások [Pa]
Nyomás a plusz kamrában:
0
0.2
0.4 0.6 Idõ [s]
0.8
0
1
0
Számított nyomás az 1. hengertérben
0.2
0.4
0.6 Idõ [s]
0.8
1
1.2
Számított nyomás az 1. hengertérben
8 10
5
8 10
5
7 10
5
7 10
5
6 10
5
6 10
5
5 10
5
5 10
5
4 10
5
4 10
5
3 10
5
3 10
5
2 10
5
2 10
5
1 10
5
1 10
5
0
Nyomások [Pa]
Nyomások [Pa]
Nyomás a mínusz kamrában:
0
0.2
0.4 0.6 Idõ [s]
0.8
0
1
Számított nyomás a 2. hengertérben
0.12
0.1
0.1
0.08 0.06 0.04 0.02 0
0.2
0.4
0.6 Idõ [s]
0.8
1
1.2
tmb j
0.12
0
0.2
Számított nyomás a 2. hengertérben
Dugattyú elmozdulás [m]
Dugattyú elmozdulás [m]
Dugattyúelmozdulás:
0
0.4 0.6 Idõ [s]
0.8
0.08 0.06 0.04 0.02 0
1
Számított a) pluszérték mozgás
0
0.2
0.4
0.6 Idõ [s]
0.8
Számított érték b) mínusz mozgás
- 104 -
1
1.2
1 10 5 8 10
6
1 10 5 8 10
6
5
6 10 5 4 10 5 2 10
6 10 5 4 10 5 2 10
Nyomáskülönbség [Pa]
Nyomáskülönbség [Pa]
Nyomáskülönbségek:
0 5
2 10 5 4 10 5 6 10 5 8 10 6
1 10
5
0 5
2 10 5 4 10 5 6 10 5 8 10 6
0
0.2
0.4 0.6 Idõ [s]
0.8
1
1 10
Számított nyomásoknál Mért nyomásoknál
0
0.24
0.48 0.72 Idõ [s]
0.96
1.2
0.08 0.06 0.04 0.02 Dugattyú elmozdulás [m]
0
Számított nyomásoknál
a) plusz mozgás
b) mínusz mozgás
Eltérés a mért és a számított dugattyú elmozdulás között 25
25
20
20
15
15
10
10
5
5
0
0
5
5
10
0
0.02 0.04 0.06 0.08 Dugattyú elmozdulás [m]
0.1
a) plusz mozgás
10
0.1
b) mínusz mozgás
- 105 -
M3. Melléklet ÖMLESZTETT ANYAGOK KIÖMLÉSÉVEL KAPCSOLATOS MÉRÉSEK M3/1. Mérési eredmények búzára
Kiömlõ anyagmennyiség [kg]
0.5
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
0.8
1
0.8
1
0.8
1
Brita) Szabvány szerint Hs = 0,2 m
Kiömlõ anyagmennyiség [kg]
0.5
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
Britb) Szabvány szerint Hs = 0,4 m 0.5
Kiömlõ anyagmennyiség [kg]
Brit Szabvány szerint Johanson szerint Beverloo szerint Janssen szerint Oldal szerint Mért értékek
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
Britc) Szabvány szerint Hs = 0,6 m
- 106 -
M3/2. Mérési eredmények PVC granulátumra
Kiömlõ anyagmennyiség [kg]
0.5
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
0.8
1
0.8
1
0.8
1
Brita) Szabvány Hs = 0,2szerint m
Kiömlõ anyagmennyiség [kg]
0.5
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
Brit b) Szabvány szerint Hs = 0,4 m 0.5 Kiömlõ anyagmennyiség [kg]
Brit Szabvány szerint Johanson szerint Beverloo szerint Janssen szerint Oldal szerint Mért értékek
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
Britc) Szabvány szerint Hs = 0,6 m
- 107 -
M3/3. Mérési eredmények polietilén granulátumra
Kiömlõ anyagmennyiség [kg]
0.5
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
0.8
1
0.8
1
0.8
1
Brita) Szabvány szerint Hs = 0,2 m
Kiömlõ anyagmennyiség [kg]
0.5
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
Britb) Szabvány szerint Hs = 0,4 m 0.5 Kiömlõ anyagmennyiség [kg]
Brit Szabvány szerint Johanson szerint Beverloo szerint Janssen szerint Oldal szerint Mért értékek
0.4
0.3
0.2
0.1
0
0
0.2
0.4 0.6 Idõ [s]
Britc) Szabvány Hs = 0,6szerint m
- 108 -
M4. Melléklet AZ ADAGOLÁSI FOLYAMAT MÉRÉSI ÉS SZÁMÍTÁSI EREDMÉNYEI M4/1. Kiömlő mennyiségek mért és számított értékei búzára
Kiömlött anyagmennyiség [kg]
0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
2.4
Számított értékek a) pt = 5∙105 Pa
Kiömlött anyagmennyiség [kg]
0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
2.4
1.6
1.8
2
2.2
2.4
Számított értékek b) pt = 6∙105 Pa 0.9 Kiömlött anyagmennyiség [kg]
Számított értékek Mért értékek
0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
Számított értékek c) pt = 7∙105 Pa
- 109 -
M4/2. Kiömlő mennyiségek mért és számított értékei PVC granulátumra
Kiömlött anyagmennyiség [kg]
0.6 0.5 0.4 0.3 0.2 0.1 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
2.4
Számított értékek a) pt = 5∙105 Pa
Kiömlött anyagmennyiség [kg]
0.6 0.5 0.4 0.3 0.2 0.1 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
2.4
1.6
1.8
2
2.2
2.4
Számított értékek b) pt = 6∙105 Pa
0.6 Kiömlött anyagmennyiség [kg]
Számított értékek Mért értékek
0.5 0.4 0.3 0.2 0.1 0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
Számított értékek c) pt = 7∙105 Pa
- 110 -
M4/3. Kiömlő mennyiségek mért és számított értékei polietilén granulátumra
Kiömlött anyagmennyiség [kg]
0.4
0.3
0.2
0.1
0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
2.4
Számított értékek a) pt = 5∙105 Pa
Kiömlött anyagmennyiség [kg]
0.4
0.3
0.2
0.1
0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
1.6
1.8
2
2.2
2.4
1.6
1.8
2
2.2
2.4
Számított értékek b) pt = 6∙105 Pa
0.5 Kiömlött anyagmennyiség [kg]
Számított értékek Mért értékek
0.4 0.3 0.2 0.1
0
0
0.2
0.4
0.6
0.8
1
1.2 1.4 Idõ [s]
Számított értékek c) pt = 7∙105 Pa
- 111 -
M5. Melléklet A MÉRÉSEKHEZ HASZNÁLT BERENDEZÉSEK M5/1. A mérőrendszer
M5/2. Az elzáró-szerkezet
- 112 -